有限元方法编程-第5章new

《有限元方法编程-第5章new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第5章平面桁架分析程序§1平面桁架计算模型•1、平面桁架单元刚度矩阵平面桁架单元1020u1u2v1v2平面桁架单元刚度矩阵EAEA0000ll12EI6EI12EI6EIX1032032u1llllYEIEIEIEIv10640621M1l2ll2l1XEAEAu200002Yllv22M2012EI6EI012EI6EI23232llll6EI2EI6EI4EI0022llllEAE

2、A00lle0000kEAEA00ll0000cossin00sincos00T00cossin00sincoseTekTkT§2平面桁架分析程序设计•1.变量说明(1)基本信息NJ--结点总数NE–单元总数NS–约束结点总数NP--结点荷载总数NW–总刚矩阵半带宽（NW=2*(单元结点的最大号差+1））NN–总自由度数（NN=2*NJ）(2)单元信息EF(NE,1)–单元始端结点号EF(NE,2)–单元末端结点号EF(NE,3)--单元杆长EF(NE,4)--单元截面面积EF(NE,5)

3、--单元弹性模量EF(NE,6)--单元局部坐标系夹角(3)约束结点信息JS(NS,1)–约束结点号JS(NS,2)–x方向约束信息JS(NS,3)–y方向约束信息(4)结点荷载信息PJ(NP,1)--荷载作用结点号PJ(NP,2)–荷载作用方向PJ(NP,3)--荷载值(5)其它及结果KS(NN,NW)–结构总刚度矩阵PS(NN)–结构荷载向量，后存储位移向量EG(NE)–单元抗拉线刚度PN(NE)–单元轴力PR(NJ,2)–支座反力KE(4,4)–单元刚度矩阵IJ(4)–单元定位向量开始主程序数组说明输入原始数据子程序INDATA回送原始数据子程序PTDATA计算半带宽子

4、程序BANDNW形成总刚度矩阵子程序TOTLKS形成荷载向量子程序TOTLPS引入支承条件子程序RESTCD解方程组子程序GAUSS2计算轴力及支座反力子程序INTEPN输出位移、轴力及反力子程序PTDPPN结束子程序名INDATA输出变量NJ,NE,NS,NP,EF,JS,PJ,EG输入基本信息NJ,NE,NS,NP输入单元信息EF(NE,6)＝＞EG(NE)输入约束结点信息JS(NS,3)输入结点荷载PJ(NP,3)子程序结束子程序名PIDATA输入变量NJ,NE,NS,NP,EF,JS,PJ输出基本信息NJ,NE,NS,NP输出单元信息EF(NE,6)输出约束结点信息J

5、S(NS,3)输出结点荷载PJ(NP,3)子程序结束子程序名BANDNW输入变量NE,EF输出变量NWNW<=0循环N=1,NEIJ<=ABS(EF(N,1)-EF(N,2))是判断IFIJ＞NW否NW<=IJNW<=2＊(NW+1)子程序结束子程序名TOTLKS输入变量NE,NW,EF,EG输出变量KS循环N=1,NE子程序TOTLKE(N,EF,EG,KE)IE,JE<=EF(N,1~2)IJ(1~4)<=(IE,JE)循环I=1,4循环J=1,4判断IJ(J)＜IJ(I)是IW<=IJ(I)JW<=IJ(J)-IJ(I)+1KS(WI,JW)<=KS(IW,JW)+KE

6、(I,J)子程序结束子程序名TOTLKS输入变量N,EF,EG输出变量KEEH<=EF(N,6)＊3.14159/180CX,CY<=(EH)KE(1~4,1~4)<=(CX,CY)KE(1~4,1~4)<=KE()＊EG(N)子程序结束子程序名TOTLPS输入变量NP,PJ输出变量PS循环N=1,NPJJ,JD<=PJ(N,1~2)判断JD=1是否JJ<=2＊JJ-1JJ<=2＊JJPS(JJ)<=PS(JJ)+PJ(N,3)子程结束子程序名RESTCD输入变量NS,NW,JS,KS,PS输出变量KS,PS循环N=1,NSJJ<=JS(N,1)循环J=1,2判断IFJS(N

7、,J+1)=0是JD<=JJ＊2-2+J接下页接上页循环I=1,NWKS(JD,I)<=0循环I=1,JDJW<=JD-I+1是判断IFJW＞NWKS(I,JW)<=0PS(JD)<=0KS(JD,1)<=1子程结束子程序名GAUSS2输入变量NN,NW,KS,PS输出变量PS循环N=1,NN-1判断IFNN≤N+NW-1是否IM<=NNIM<=N+NW-1循环I=N+1,IMIW<=I-N+1C<=KS(N,IW)/KS(N,1)循环J=1,IN-I+1JW<=J+I-N接下页接上页KS(I,J)<=