资源描述：
《有限元法基础与程序设计framefor》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、C框架。(框架问题)-维LOC(3,2)IFIXo(6),(3)区域,弗林特(3),残雪(4),CY(4),1F(12),FR(12)、铁(3、6)常见的神经网络、NE、NDNFIXEC打开文件(5=“F1.TXT”，地位二’老')打开文件(6二0UTFE)阅读(5*)神经网络,NE,NDNFIX,E写(6105)神经网络,NENDNFIXE105格式(4X,nneNFIXE/4I5,Ell)阅读(5*)(LOC(我,1),代码行(我,2),(我)，弗林特(1),1=1,NE)写(6108)108格式(/IX,“元素N0DE1N0DE2
2、区域FINT”)写(6110)(L0C(我,1),我疯狂的(我,2),(我)，弗林特(1),1=1,NE)110格式(1X3162el2.4)阅读(5*)(CX(J),CY(J),J二1,NN)写(6112)112格式(/4x”节点X-C00RDY-COORD”)写(6115)(J,残雪(J),CY(J),J二1,NN)115格式(4x,13,3x2el2.4)阅读(5*)(IFIXo(K),K二1,NFIX)写(6118)(IFIXo(K),K二1,NFIX)118格式(/3x的IFIXo二6预告)C调用春秋国旅(LOC,IFIX用于
3、区域，弗林特,残雪,CY,F,FR,FE)C写(6130)130格式(/4x,“节点”,5x,X-DISP,8x,Y-DISP,9x,THITA)写(6135)(F(3*我2),F(3张)，F(3*我)，我二1,NN)135格式(2x、153el5.6)写(6140)140种格式(/4X,反应节点的力，1/4x“节点”5x,X-L0AD,8x,Y-L0AD,9x,“时刻”)写(6145)(FR(3*我2),FR(3张)，FR(3*我)，我二1,NN)145格式(2X、153el5.6)写(6150)150格式(/4x,“元素”,4x,A
4、XI-F,7x,SHEAR-Q,7x,MOMENT-M)10我二1,不10写(6155)、铁俄，1)，铁俄,2)、铁(我,3)、铁(我,4)、铁俄,5)、铁俄,6)155格式(2x,预告3el5・6/6x3el5.6)停止结束子例程春秋国旅(LOCIFIX用于区域，弗林特,残雪,CY,F,FR,FE)维LOC(NE,2),IFIXo(NFIX)、区域(NE),弗林特(NE),残雪(NN),CY(NN),1F(ND),FR(ND)>铁(NE,6),2Pl(3、6),P2(3、6),门将(12日12),GK1(12日12),(3),TT(3
5、、6、6),食品(3、6、6),3T(6,6),CSS(3),SNN(3),埃克(6,6),EKL(6,6),EKK(3、6、6),XX(6),英航(6,6)常见的神经网络、NE、NDNFIXEC不，5I二1吗11二LOC(我,1)12二LOC(我，2)X12二CX(I2)CX(I1)日元二CY(I2)CY(I1)艾尔(I)=VX12**2+日元**2)CSS⑴二X12/AL俄)SNN⑴二日元/AL俄)5继续C对于中心力，设为一个节点阅读(5*)NP,NVD如果(np.eq.0)到15做102二1,NP阅读(5*),F(3*我2),(3
6、张),F(3*我)10继续15日继续对于一个元素的垂直分布加载如果(nvd.eq.0)到25做202二1,NVD吗阅读(5*)即问II二LOC(l)12二LOC(2)Pl(3*Il-1)=Q*艾尔(IE)/2OPl(3*11)二Q*艾尔(IE)**2/12oPl(3*12-1)二Q*艾尔(IE)/2oPl(3*12)二-q*艾尔(IE)**2/12o20日继续25日继续C做30I二1,ND做30J二1,ND吗30g-(I,J)二0.0C不，100I二1吗6,35J二1吗6,35K二1吗EKL(J,K)二0.0埃克(J,K)二0.035T
7、二0.0(J,K)C宝莲寺二E*区(我)/AL(我)EIL1二E*弗林特(我)/AL(我)EIL2二E*弗林特(我)/AL(I)**2EIL3二E*弗林特(我)/AL(I)**3EKL(1,1)二宝莲寺EKL⑵2)二12o*EIL3EKL(3)=4o*EIL1EKL(4,4)二宝莲寺EKL(5,5)=12o*EIL3EKL(6,6)二4*EIL1EKL(2)=6o*EIL2EKL(4,1)=宝莲寺EKL(5,2)=-12*EIL3EKL(5,3)=6o*EIL2EKL⑵二6o*EIL2EKL(6,3)=2*EIL1EKL(6,5)=6o
8、*EIL25,笫二50二1吗做50个JJ二2+1,650EKL(II,JJ)=EKL(JJ,(二)6,55K二1吗6,55L二1吗EKK(K,L)二EKL(K,L)C写(6140),K,L,EKL(K,L)
