有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)

有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)

ID:11934968

大小:40.50 KB

页数:10页

时间:2018-07-15

有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)_第1页
有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)_第2页
有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)_第3页
有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)_第4页
有限元法基础与程序设计 truss.for(basic and program design of finite element method truss.for)_第5页
资源描述:

《有限元法基础与程序设计 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-

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。