有限單元法編程之網(wǎng)絡(luò)結(jié)構(gòu)

FiniteElementMethod18．圖示一結(jié)構(gòu)網(wǎng)絡(luò),共有10個(gè)結(jié)點(diǎn),9個(gè)單元和四個(gè)支桿,在結(jié)點(diǎn)10作用豎向荷載10KN,按平面應(yīng)力問(wèn)題計(jì)算各結(jié)點(diǎn)位移和應(yīng)力.設(shè)E=1,μ=0,厚度t=1*******************************************************************************程序如下：cPLANSTRESS/STRAINPROBLAMANALYSISPROGRAMDIMENSIONLND(500,3),X(200),Y(200),&JZ(50,3),PJ(50,3),P(500),&TL(20),AK(500,100),AKE(6,6)OPEN(6,FILE='1.txt',STATUS='OLD')OPEN(8,FILE='2.txt',STATUS='OLD')READ(6,10)TL10FORMAT(20A2)READ(6,*)NJ,NE,NZ,NPJ,IPSWRITE(8,10)TLIF(IPS.EQ.0)WRITE(8,20)IF(IPS.EQ.1)WRITE(8,30)20FORMAT(/1X,'PLATESTRESSPROBLEM')30FORMAT(/1X,'PLATESTRAINPROBLEM')公嫁NPJ0=N電PJ脈育IF(NPJ斧.EQ.0)鑼NPJ=1將慰CALL腸INPUT禿(生NJ,NE,葉NZ文,矛ND堪,NPJ,N融PJ0,IP向S,E,PR事,予&夕T,V,LN剩D,X,Y,慚JZ,PJ)腐燥N=2*NJ扶痛DO40秋I=1,N喂演DO40謝J=1,ND每40囑AK(I,J觀)=0.0剩茂DO50瘡IE=1,N高E色役CALLM媽KE(IE,肅NJ,LND謊,X,Y,E敗,PR,T,服AKE)適役DO50岡I=1,3比艱DO50璃II=1,2軌我IH=2*(途I-1)+I由I微聰IDH=2*丑(LND(I疏E,I)-1本)+II促靈DO50捆J=1,3絞朵DO50花JJ=1,2墊晶L=2*(J伐-1)+JJ剛惡IL=2*(歪LND(IE籮,J)-1)載+JJ程脅IDL=IL瓣-IDH+1柜蒼IF(IDL煌.LE.0)筒GOTO5晶0式晨AK(IDH壺,IDL)=袋AK(IDH嫁,IDL)+性AKE(IH局,L)猶50陶CONTIN珠UE秩令CALLM默F(xiàn)(NJ,N賺E,NPJ,標(biāo)NPJ0,N敞,T,V,L信ND,X,Y楚,PJ,P)作驚CALLR伸KR(星NZ追,鹽ND很,N,畫JZ分,省AK路,P)醫(yī)宵CALLS勻LOV(NJ曠,N,ND,丹AK,P)絹柄CALLM瞧ADE(NE南,NJ,N,拐E,PR,L叼ND,X,Y嚷,P)涂慈CLOSE(錫6)忠申CLOSE(鴨8)牽覺(jué)STOP輛溉END乖拴SUBROU閃TINE抄INPUT析(碧N(yùn)J,NE,箭NZ,ND,沉NPJ,NP隆J0,IPS王,E,PR,她T,V,稀&太LND,X,躬Y,JZ,P妙J)舟暑DIMENS該IONLN頌D(500,資3),X(2流00),Y(愉200),J柏Z(50,3怒),PJ(5竄0,3)錄薄READ(6攝,*)E,P拴R,T,V匹爪READ(6非,*)((L弦ND(I,J弟),J=1,您3),I=1圣,NE)來(lái)藏READ(6冰,*)(X(水I),Y(I鞋),I=1,譽(yù)NJ)裂拿READ(6島,*)((J換Z(I,J)純,J=1,3廣),I=1,愁NZ)妄翁IF(NPJ田鋪)GOTO盯10尺跡READ(6老,*)((P膝J(I,J)蚊,J=1,3損),I=1,允NPJ)暑10購(gòu)ND=0錦史DO20悉IE=1,N螺E棋戴DO20角I=1,3肅鉗DO20漂J=1,3賤油IW=IAB渠S(LND(券IE,I)-懇LND(IE窄,J))產(chǎn)閉IF(ND.報(bào)GE.IW)眾GOTO原20脊紋ND=IW懷20吳CONTIN悄UE漏駕ND=(ND削+1)*2擠倡I(xiàn)F(IPS刑.EQ.0)瘡GOTO3蠻0報(bào)沈E=E/(1滿.0-PR)山旁PR=PR/虧(1.0-P梢R)收30阻RETURN毅容END渣斬SUBROU洪TINEM思KE(IE,尋NJ,LND尊,X,Y,E家,PR,T,獄AKE)棒灌DIMENS腎IONLN怨D(500,擦3),X(2疲00),Y(降200),B生(3,6),能D(3,3)秩,S(3,6闖),勻&穿AKE(6,慎6)合滿CALLM伯A(IE,N刑J,LND,隔X,Y,AE踢)皮跟CALLM多D(E,PR檔,D)更站CALLM臟B(IE,N湊J,LND,萄X,Y,AE鏡,B)傘初DO10要I=1,3濕玻DO10沃J=1,6蛛媽S(I,J)護(hù)=0.0臣飛DO10訂K=1,3要10弦S(I,J)頸=S(I,J誠(chéng))+D(I,妖K)*B(K浪,J)滅鞠DO20眼I=1,6衣戀DO20孝J=1,6族乞AKE(I,吩J)=0.0庸奶DO20駱K=1,3毫20貼AKE(I,惱J)=AKE母(I,J)+睬S(K,I)具*B(K,J縫)*AE*T筆利RETURN艦執(zhí)END捎疫SUBROU茫TINEM歲A(IE,N褲J,LND,威X,Y,AE釀)件涌DIMENS算IONLN免D(500,坐3),X(2攻00),Y(救200)翼容I=LND(盯IE,1)樂(lè)始J=LND(午IE,2)未丘K=LND(賓IE,3)恒何XIJ=X(諸J)-X(I因)兆不YIJ=Y(膏J)-Y(I趨)乒呼XIK=X(真K)-X(I疏)踐車YIK=Y(輝K)-Y(I隆)鴉慈AE=0.5感*(XIJ*后YIK-XI缸K*YIJ)駕座RETURN囑走END江孝SUBROU清TINEM魯D(E,PR詞,D)底糞DIMENS僻IOND(不3,3)尋鈔DO10夢(mèng)I=1,3搭鑄DO10壁J=1,3陜10冬D(I,J)守=0.0鑒哪D(1,1)才=E/(1.診0-PR*P廁R)喊饅D(1,2)寨=E*PR/連(1.0-P稿R*PR)束貿(mào)D(2,1)緊=D(1,2農(nóng))訓(xùn)豎D(2,2)荷=D(1,1爪)郊嬌D(3,3)害=0.5*E槽/(1.0+弟PR)岸披RETURN語(yǔ)犯END撫乖SUBROU躍TINEM鳥B(IE,N帖J,LND,盆X,Y,AE唉,B)郵辨DIMENS森IONLN題D(500,筒3),X(2瞎00),Y(酷200),B叼(3,6)病舟I=LND(匆IE,1)吧閉J=LND(遺IE,2)林單K=LND(歲IE,3)嗎福DO10悄II=1,3溪裂DO10趁JJ=1,6攪10老B(II,J預(yù)J)=0.0栽駐B(1,1)座=Y(J)-宜Y(K)改勞B(1,3)罰=Y(K)-副Y(I)漲途B(1,5)共=Y(I)-歲Y(J)治顛B(2,2)待=X(K)-舌X(J)統(tǒng)蘿B(2,4)蘇=X(I)-柳X(K)寧故B(2,6)災(zāi)=X(J)-榴X(I)拋黑B(3,1)趙=B(2,2犬)庫(kù)也B(3,2)識(shí)=B(1,1拾)礦圖B(3,3)決=B(2,4菜)網(wǎng)豎B(3,4)炮=B(1,3筋)元市B(3,5)鑼=B(2,6載)鬧盾B(3,6)繼=B(1,5葬)蘋果DO20慶I1=1,3驕瞧DO20剝J1=1,6覺(jué)20擇B(I1,J問(wèn)1)=0.5圾/AE*B(鵲I1,J1)排利RETURN鵲久END膀黨SUBROU脂TINEM夾F(NJ,N器E,NPJ,摩NPJ0,N轟,T,V,L訊ND,X,Y擺,PJ,P)炕漏DIMENS叔IONLN擇D(500,練3),X(2許00),Y(中200),P警J(50,3恢),P(50續(xù)0)摩仗DO10化I=1,N協(xié)10災(zāi)P(I)=0轉(zhuǎn).0燥鈔IF(NPJ件時(shí))GOTO拘30元結(jié)DO20嘗I=1,NP走J艱蔑II=PJ(困I,1)屋住P(2*II爛-1)=PJ甜(I,2)挪20撇P(2*II炒)=PJ(I隨,3)宿30遮IF(V.E臘Q.0)G盞OTO50床抹DO40安IE=1,N露E灣為CALLM維A(IE,N愈J,LND,堤X,Y,AE慰)出戰(zhàn)PE=-V*慈AE*T/3浙.0扒祝DO40賽I=1,3慧夕II=LND貫(IE,I)注40癥P(2*II源)=P(2*獎(jiǎng)II)+PE秒50此RETURN耕茄END亂夏SUBROU諷TINER滿KR(蛛NZ鈔,沿ND駁,N,泉JZ泰,粥AK封,P)拳申DIMENS制IONP(罩500),J罪Z(50,3貨),AK(5朱00,100咐)舍嶺DO30厘I=1,NZ贏稅IR=JZ(在I,1)娃市DO30麗J=2,3攔春IF(JZ(貴I,J).E死Q.0)GO靜TO30肚遠(yuǎn)II=2*I勢(shì)R+J-3以根AK(II,帆1)=1.0跨燭DO10畢JJ=2,N般D羅10格AK(II,彈JJ)=0.屢0期剃IF(II.婚GT.ND)閥JO=ND平饑IF(II.職LE.ND)潑JO=II頂萌DO20織JJ=2,J耐O臂20足AK(II-膚JJ+1,J朵J)=0.0循召P(II)=憶0.0全30糞CONTIN少UE貓缺RETURN王辜END斗沒(méi)SUBROU逃TINES籠LOV(NJ拆,N,ND,鄉(xiāng)AK,P)羅往DIMENS量IONP(恭500),A暗K(500,炎100)報(bào)滲NJ1=N-打1樹愉DO20混K=1,NJ員1害韻IF(N.G連T.K+ND愈-1)IM=佛K+ND-1洗眉IF(N.L漆E.K+ND館-1)IM=御N帆國(guó)K1=K+1誰(shuí)狐DO20瓣I=K1,I喝M靠咐L=I-K+敬1詞酸C=AK(K席,L)/AK當(dāng)(K,1)秤荒IW=ND-銹L+1右倉(cāng)DO10萍J=1,IW屆獲M=J+I-毯K始10葵AK(I,J厭)=AK(I將,J)-C*肚AK(K,M煮)絲20柔P(I)=P控(I)-C*計(jì)P(K)瓶志P(N)=P刺(N)/AK搏(N,1)目通DO40物I1=1,N耽J1腥挖I=N-I1原妻IF(ND.猴GT.N-I梢+1)JO=惑N-I+1堆亞IF(ND.公LE.N-I擋+1)JO=駝ND霉添DO30賴J=2,JO潔爬K=J+I-辣1阿30叔P(I)=P像(I)-AK體(I,J)*雖P(K)打40度P(I)=P根(I)/AK角(I,1)銳株WRITE(蹤8,50)掙50占FORMAT病(/8X,5菌('**'巖),'RES鷹ULTEO侄FCALC垃ULTION瓦',尼&我5('**跳')//1X結(jié),'NODA千LDISP蹈LACEME轟NTS'//室3X,'NO企DE',宜&微5X,'X-倉(cāng)DISP.'返,8X,'Y憂-DISP.服')根匆DO60諷I=1,NJ瑞60高WRITE(序8,70)I雞,P(2*I訓(xùn)-1),P(急2*I)睛70納FORMAT賠(1X,I5尸,2E15.高6)住所RETURN忌社END圈犁SUBROU初TINEM凱ADE(NE宣,NJ,N,悠E,PR,L縱ND,X,Y兄,P)蝶珍DIMENS霞IONLN宮D(500,證3),X(2捕00),Y(踢200),D住(3,3),愈&夾B(3,6)慰,S(3,6地),ST(3頭),P(50徹0),DE(喂6)些順WRITE障(8,10)蹲10耳FORMAT餡(/1X,'貪ELEMEN蠶TSTRE屠SSES'/警)聾耽CALLM德D(E,PR艘,D)蛙宋DO50癢IE=1,N消E乎盈CALLM正A(IE,N攔J,LND,介X,Y,AE底)伍孕CALLM本B(IE,N虜J,LND,知X,Y,AE友,B)播茄DO20肌I=1,3僅筑DO20錄J=1,6囑忌S(I,J)廣=0.0報(bào)置DO20蟲K=1,3體20墓S(I,J)夫=S(I,J恥)+D(I,尿K)*B(K擋,J)篩洞DO30吸I=1,3乎捐DO30拌J=1,2譜寒IH=2*(耗I-1)+J玩扯IW=2*(周LND(IE膝,I)-1)估+J脾30懂DE(IH)筍=P(IW)傳衛(wèi)DO35腦I=1,3停35氧ST(I)=鼠0.0吧遍DO40糠I=1,3艦饑DO40行J=1,6乎40脈ST(I)=山ST(I)+哀S(I,J)秀*DE(J)奮檔STX=ST蠅(1)塌續(xù)STY=ST疏(2)蘭嘴TXY=ST墓(3)特50班WRITE(膏8,60)I振E,STX,述STY,TX謎Y猴60等FORMAT舉(1X,'E使LEMENT雀NO.='與,I5/3X撤,'STX=界',E15.映6,添&掠2X,'ST壩Y=',E1厚5.6,2X魄,'TXY=拆',E15.峽6)錫蹤RETURN狗扛END押******角******澇******摩******食******蔽******君******偽******齡******雕******徒******切***振輸入文朝件帆1.da牲t蚊：銳10,9,工2,1,0親1.0,0搜.25,1采.0,0.得0衫1,2,5,隱2,3,6,廟3,4,7,管3,7,6,持2,6,5,呢5,6,8著6,7,9,斧6,9,8,廣8,9,10瓦0.0,0.柿0,2.0,炕0.0,4.遍0,0.0,占6.0,0.俗0,1.0,攝2.0,3.壺0,2.0,各5.0,2.等0幸2.0,4才.0,4.壞0,4.0負(fù),3.0,絮6.0辟1,1,打1,4,1,融1右10.0,暴0.0,1旨0.0岔******霉******丘******喬******譽(yù)******否******駛******借******半******你******川******盲***暮PLATE姑STRESS膽PROBL桃EM嶼轟***鍛****非***R位ESULTE爐OFCA傷LCULTI適ON***冷****靈***菠NODAL屈DISPL鴨ACEMEN幫TS向NOD怕E筍X-DISP姐.富Y-D怨ISP.皇1啊.0粉00000E稠+00惜.0000青00E+00便2翁.1群09411E韻+01撕.1775陸65E+02晴3捷-.1途09411E趨+01引.1775濾65E+02鴉4屠.0述00000E害+00奪.0000談00E+00加5恐-.1砌64116E雖+01疲.1567壤85E+02費(fèi)6鴨-.5住74256E叢-06額.2095裕99E+02法7稿.1塞64116E壘+01隙.1567劑85E+02圖8輛.8察20585E施+00挺.2531匹27E+02況9撞-.8衰20577E伏+00傻.2531籮27E+02膛10遞.8宜42274E翻-05邀.4447逃30E+02物ELEME竄NTSTR井ESSES寬ELEME幻NTNO.棋=1純STX釘=.壯149023塔E+01束STY=游.377嚴(yán)270E+0赤1TXY柳=.脊311365卷E+01訂ELEME趟NTNO.玩=2壘STX徒=-.彈739929堡E+00夫STY=麥.141影672E+0萌1TXY恒=-.皮142490滴E-07殊ELEM