有限元法大作業(yè)

《有限元法大作業(yè)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、有限元法大作業(yè)一平面剛架的程序用VisualC++編制的平面剛架的源程序如下：///////////////////////////////////////////////////////程序開始//////////////////////////////////////////////////////////////////#include"iostream.h"#include"math.h"#include"stdlib.h"#include"conio.h"http://*****************//聲明必要變量//*****************#definePI3.14159

2、2654intNE;//單元數(shù)intNJ;//節(jié)點(diǎn)數(shù)intNZ;//支承數(shù)intNPJ;//有節(jié)點(diǎn)載荷作用的節(jié)點(diǎn)數(shù)intNPF;//非節(jié)點(diǎn)載荷數(shù)intHZ;//載荷碼intE;//單元碼intfangchengshu;doubleF[303];//各節(jié)點(diǎn)等效總載荷數(shù)值intdym_jdm[100][2];//單元碼對應(yīng)的節(jié)點(diǎn)碼:dym_jdm[][0],dym_jdm[][1]分為前后節(jié)點(diǎn)總碼intzhichengweizhi[300];//記錄支持節(jié)點(diǎn)作用點(diǎn)的數(shù)組intfjzhzuoyongdanyuan[100];//非節(jié)點(diǎn)載荷作用單元intfjzhleixing[100];//非節(jié)

3、點(diǎn)載荷類型:1-均布,2-垂直集中,3-平行集中,4-力偶,5-角度集中doublefjzhzhi[100];//非節(jié)點(diǎn)載荷的值doublefjzhzuoyongdian[100];//非節(jié)點(diǎn)載荷在各竿的作用點(diǎn)doublefjzhjiaodu[100];//非節(jié)點(diǎn)載荷作用角度intjdzhzuoyongdian[100];//節(jié)點(diǎn)載荷作用的節(jié)點(diǎn)數(shù)組doublejiedianzaihe[101][3];//節(jié)點(diǎn)載荷值,其jiedianzaihe[][0]--jiedianzaihe[][2]分別為U,V,Mdoublezhengtigangdu[303][303];//整體剛度數(shù)組doub

4、lechangdu[100];//各單元竿長數(shù)組doublejiaodu[100];//各單元角度數(shù)組doubletanxingmoliang[100];//各單元彈性模量數(shù)組doubleJ_moliang[100];//各單元J模量數(shù)組doublemianji[100];//各單元面積數(shù)組doubleweiyi[303];//記錄各個(gè)節(jié)點(diǎn)位移的數(shù)組doubledy_weiyi[100][6];//各個(gè)單元在局部坐標(biāo)系中的位移數(shù)組dy_weiyi[i][0]-dyweiyi[i][6]分別為第i+1單元的u1,v1,@1,u2,v2,@2doubledy_neili[100][6];//

5、各個(gè)單元在局部坐標(biāo)系中的固端內(nèi)力dy_weiyi[i][0]-dyweiyi[i][6]分別為第i+1單元的U1,V1,M1,U2,V2,M2doublegan_neili[100][6];//各個(gè)單元的竿端內(nèi)力數(shù)組,gan_neili[i][6]表示第i+1單元的6內(nèi)力.//*******************//一些處理過程的函數(shù)31//*******************voidinput();//數(shù)據(jù)的輸入voidzonggang();//計(jì)算總剛度,存放于zhengtigangdu[][]數(shù)組中voidzongzaihe();//計(jì)算等效總節(jié)點(diǎn)載荷voidzhichengy

6、inru();//引入支承條件voidjsweiyi();//求各個(gè)節(jié)點(diǎn)位移voidjs_dy_weiyi();//求局部坐標(biāo)系中的位移voidganduanneili();//求竿端內(nèi)力voiddy_gangdu(inti,doubledg[6][6]);//求單元在局部坐標(biāo)系中的單剛voidjs_T_T1(inti,doubleT[6][6],doubleT1[6][6]);//求單元的轉(zhuǎn)換矩陣及其逆陣//************//主函數(shù)//************voidmain(){input();cout<<"輸出結(jié)果"<

7、＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝”<