资源描述:
《有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、有限元法基础与程序设计truss.for(Basicandprogramdesignoffiniteelementmethodtruss.for)C-------桁架。---(桁架问题)--------------------维禄(3,2),国泰航空(3)、环磷酰胺(3)、(3)、铝等等(3)、地区(3),1GK(6,6),F(6),FF(3)常见的神经网络,NE,Nd,nfix,EC-----------------------------------------------------------打开(5、文件='t1.txt”,地位=老')打开(6、文
2、件='outt1。”)阅读(5,*)NN,NE,Nd,nfix,E写(6105)NN,NE,Nd,nfix,E105格式(2x,'nnNE和E/4i5nfix,e11.4)读(5,*)(LOC(i,1),LOC(i,2),区域(i),i=1,NE)写作(6108)108格式(/1,'elementNode1NODE2地区)写入(6110)(i,LOC(i,1),LOC(i,2),区域(i),i=1,NE)110格式(1X,3i6,e12.4)读(5,*)(CX(j),赛扬(j),j=1,NN)写作(6112)112格式(/3X,'nodex-coordy-c
3、oord”)写入(6115)(J,CX(j),赛扬(j),j=1,NN)115格式(3X,i3,3x,2e12.4)阅读(5,*)(iFIX(K),K=1,nfix)写(6118)(iFIX(K),K=1,nfix)118格式(/1,'ifix=',5i4)C-----------------------------------------------------------做10I=1,ND10f(i)=0f(2)=10000C-----------------------------------------------------------调用CST(
4、LOC,CX,CY,iFix,AL地区,GK,F,FF)写作(6120)120格式(/3X,'node,6X,'x-disp,8X的'y-disp”)写(6125)(I、F(2×I-1),F(2*),i=1,NN)125格式(1X,i5,2e15.4)写作(6130)130格式(/1,'element,5X,力量”)写入(6135)(i,FF(i),i=1,NE)135格式(1X,i5,e15.4)停止结束C-----------------------------------------------------------子程序CST(LOC,CX,CY,
5、iFix,AL地区,GK,F,FF)维禄(NE,2)、CX(NN),CY(NN),iFIX(nfix)、Al(NE),1地区(NE),GK(ND,ND),F(ND)、FF(NE),EK(4,4),XX(4)常见的神经网络,NE,Nd,nfix,EC-----------------------------------------------------------做10I=1,ND做10j=1,ND10GK(i,j)=0C-----------------------------------------------------------做100I=1,NE
6、I1=LOC(I,1)I2=LOC(i,2)X12=cx(I2)-CX(I1)Y12=CY(I2)-CY(I1)Al(i)=SQRT(X12××2+Y12×2)CS=X12/铝(我)Sn=Y12铝(我)A0=区域(i)AL(i)A1=aa*csA2=******A3EK(1,1)=A1EK(2,2)=A2EK(3,3)=A1EK(4,4)=A2EK(1,2)=A3EK(1,3)=1EK(1,4)=A3EK(2,3)=A3EK(2,4)=A2EK(3,4)=A350二=2做50个JJ=1,150EK(II,JJ)=EK(JJ,II)做55k=155=1C---
7、--------------------------------------------------------用C写的(6140)I,K,L,EK(K,L)140格式(1X,我,K,L,EK,3i4,E12.5)55继续C-----------------------------------------------------------做85点=1,2nodei=LOC(我、inode)做85idofn=1,2nrows=(nodei-1)*2+idofnnrowe=(inode-1)*2+idofn做85jnode=1,2nodej=LOC(我jnod
8、e)做85jdofn=1,2军官学校=(nodej-