




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
FiniteElementMethod18.图示一结构网络,共有10个结点,9个单元和四个支杆,在结点10作用竖向荷载10KN,按平面应力问题计算各结点位移和应力.设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')难草NP悔J0=N骡PJ莫培IF悉(NPJ耽.EQ.莲0)NP骆J=1户剃CA瓶LL艇INPU汗T惭(洁NJ,N醒E,NZ圣,ND,俘NPJ,宴NPJ0秆,IPS上,E,P走R,泽别&T,文V,LN禽D,X,呆Y,JZ蚀,PJ)按探N=袋2*NJ扮捏DO射40队I=1,篇N核熊DO咏40介J=1,组ND售40水AK航(I,J勉)=0.糊0披登DO则50乔IE=1灯,NE哄类CA屠LLM拢KE(I华E,NJ享,LND朴,X,Y君,E,P藏R,T,黄AKE)沾柔DO前50献I=1,卵3务谜DO窃50跟II=1晕,2绝队IH蛾=2*(鲁I-1)睁+II作乐ID碰H=2*兰(LND进(IE,像I)-1祝)+II续篮DO随50咐J=1,雹3没召DO订50冰JJ=1科,2激它L=纺2*(J览-1)+滑JJ惜庙IL秤=2*(劲LND(式IE,J脂)-1)功+JJ班歼ID让L=IL酬-IDH躬+1门的IF悔(IDL涂.LE.记0)GO卵TO5嫂0定霜AK盛(IDH制,IDL受)=AK富(IDH参,IDL蜓)+AK会E(IH覆,L)柄50婆CO沿NTIN糠UE帝输CA腿LLM章F(NJ寺,NE,楼NPJ,咳NPJ0结,N,T避,V,L恭ND,X御,Y,P耗J,P)君斥CA色LLR恋KR(N大Z,ND葵,N,J抬Z,AK魄,P)检围CA赚LLS床LOV(奸NJ,N余,ND,答AK,P皮)剥烂CA荡LLM闷ADE(讯NE,N选J,N,戚E,PR照,LND弄,X,Y绸,P)鉴屿CL酱OSE(冬6)拣结CL缝OSE(组8)酿牛ST南OP筛棋EN燃D介奥SU车BROU朵TINE卫炮INPU明T合(励NJ,N高E,NZ眼,ND,彼NPJ,娃NPJ0惭,IPS卖,E,P某R,T,责V,如北&LN躺D,X,哑Y,JZ拨,PJ)塔仍DI采MENS予ION承LND(棋500,惹3),X进(200临),Y(靠200)效,JZ(典50,3乳),PJ吴(50,屿3)流厚RE膛AD(6回,*)E源,PR,派T,V盆唇RE镰AD(6衰,*)(没(LND股(I,J腾),J=溜1,3)蚀,I=1糕,NE)混粘RE尽AD(6购,*)(包X(I)层,Y(I对),I=商1,NJ傍)狭六RE蹈AD(6雁,*)(影(JZ(汉I,J)捡,J=1狐,3),匹I=1,行NZ)够合IF语(NPJ秧0.EQ嚷.0)拘GOTO姓10拘贼RE迎AD(6污,*)(抄(PJ(温I,J)剧,J=1别,3),殃I=1,咱NPJ)袖10淋ND希=0诊坟DO泥20晨IE=1兆,NE精遇DO漫20翅I=1,舱3份响DO框20迟J=1,晒3钥够IW矮=IAB侍S(LN协D(IE盐,I)-扩LND(御IE,J羊))狼保IF迟(ND.约GE.I答W)G英OTO右20肿县ND音=IW惹20摘CO鸽NTIN事UE掏雄ND耽=(ND酸+1)*现2仪云IF毫(IPS吧.EQ.委0)GO消TO3洁0评租E=检E/(1奸.0-P棕R)艘吃PR蜻=PR/赠(1.0弟-PR)症30贫RE剩TURN业阿EN君D收阶SU恶BROU除TINE酒MKE锐(IE,仅NJ,L肚ND,X毒,Y,E此,PR,导T,AK心E)肯胸DI臂MENS惜ION计LND(白500,脑3),X然(200织),Y(铲200)违,B(3产,6),电D(3,甲3),S雀(3,6湾),肿肉&AK访E(6,讯6)兽循CA根LLM础A(IE订,NJ,奉LND,捡X,Y,芹AE)葛四CA银LLM办D(E,爪PR,D狐)吼龄CA格LLM姓B(IE纤,NJ,象LND,稠X,Y,垫AE,B四)疗掘DO叫10头I=1,陕3驾典DO穿10查J=1,怕6服伸S(慰I,J)稀=0.0燕肥DO殊10矩K=1,木3殃10关S(慎I,J)薯=S(I潜,J)+欣D(I,捎K)*B棋(K,J化)国喂DO船20赢I=1,所6满叙DO樱20催J=1,势6辅信AK腐E(I,拨J)=0室.0认就DO锈20饮K=1,闷3第20纠AK彩E(I,鼓J)=A临KE(I羡,J)+秘S(K,仁I)*B牺(K,J励)*AE蓝*T臂芦RE货TURN野秆EN华D必竞SU隔BROU梢TINE擦MA(布IE,N信J,LN旗D,X,称Y,AE久)征尘DI钻MENS亚ION东LND(伪500,菊3),X令(200督),Y(序200)瓦鼓I=迷LND(江IE,1堆)处似J=妄LND(服IE,2却)彻挨K=苗LND(杰IE,3嫩)财棚XI架J=X(腿J)-X访(I)宇停YI薪J=Y(设J)-Y蔑(I)尾险XI弦K=X(狸K)-X缴(I)天臂YI马K=Y(受K)-Y闪(I)掩急AE毙=0.5嫂*(XI鸦J*YI叮K-XI四K*YI腊J)嫩似RE冠TURN妹奖EN恢D看厘SU微BROU番TINE捎MD(并E,PR劈,D)德顾DI勿MENS扇ION销D(3,抄3)担绒DO耀10赖I=1,脆3剩讨DO赵10娇J=1,卧3单10烫D(熔I,J)禾=0.0躁孕D(指1,1)吨=E/(娇1.0-妹PR*P狱R)铁浊D(首1,2)再=E*P绘R/(1洁.0-P矛R*PR框)狸针D(衫2,1)定=D(1辛,2)徒贿D(诱2,2)傅=D(1稠,1)筹哈D(侄3,3)楼=0.5鞠*E/(黎1.0+收PR)雕允RE材TURN话吹EN硬D无筑SU忙BROU籍TINE牵MB(渗IE,N凶J,LN毅D,X,蚁Y,AE啊,B)邮射DI聪MENS罪ION梢LND(度500,木3),X堆(200谊),Y(变200)猾,B(3温,6)仓添I=袋LND(丸IE,1丢)惩娃J=轻LND(煌IE,2能)朗斯K=筛LND(嚼IE,3牙)电DO1扔0II惰=1,3念DO1春0JJ铃=1,6筐10凝B(气II,J岂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)底B(2,晴4)=X义(I)-点X(K)捧B(2,呢6)=X缝(J)-摇X(I)江B(3,卧1)=B景(2,2拿)理B(3,蒙2)=B赤(1,1想)菜B(3,哑3)=B束(2,4绿)猪B(3,截4)=B肚(1,3赚)蛇B(3,巷5)=B查(2,6她)汪B(3,朗6)=B玉(1,5勉)背DO2毛0I1曲=1,3框DO2呆0J1列=1,6氧20黎B(永I1,J浸1)=0想.5/A屡E*B(良I1,J踏1)训RETU秋RNEND净SUBR异OUTI未NEM立F(NJ防,NE,录NPJ,监NPJ0嫂,N,T闭,V,L盯ND,X烦,Y,P用J,P)床DIME决NSIO贯NLN斜D(50渗0,3)腔,X(2片00),斗Y(20存0),P耕J(50爷,3),涉P(50查0)劝DO1疏0I=衣1,N夏10提P(喘I)=0秩.0哗IF(N艇PJ0.蒸EQ.0称)GO卖TO3躺0岔DO2债0I=说1,NP型J肿II=P次J(I,冈1)他P(2*菊II-1讽)=PJ游(I,2垒)佳20咐P(谈2*II捕)=PJ辟(I,3镜)呼30惨IF课(V.E誉Q.0)以GOT批O50罢DO4识0IE尿=1,N桃E借CALL酒MA(兽IE,N纱J,LN毛D,X,孕Y,AE经)职PE=-虾V*AE麦*T/3秒.0您DO4暑0I=乱1,3乖II=L海ND(I奔E,I)懒40证P(夫2*II楚)=P(然2*II音)+PE续50兼RE助TURNEND棋SUBR矮OUTI庸NER授KR(N秃Z,ND晶,N,J揭Z,AK挨,P)锡DIME研NSIO卡NP(怀500)势,JZ(非50,3乳),AK舌(500停,100圾)彻DO3蝴0I=排1,NZ外IR=J拘Z(I,羊1)冒DO3刊0J=翼2,3宴IF(J最Z(I,厌J).E至Q.0)框GOTO驼30隔II=2龙*IR+糕J-3掌AK(I槽I,1)梳=1.0柴DO1广0JJ供=2,N震D泉10肌AK鞠(II,促JJ)=息0.0绘IF(I市I.GT焰.ND)趣JO=N捆D票IF(I挤I.LE京.ND)竿JO=I著I革DO2伯0JJ倍=2,J近O刮20得AK忠(II-逗JJ+1壁,JJ)奸=0.0彼P(II赵)=0.吩0样30够CO生NTIN独UE根RETU帐RNEND纱SUBR准OUTI蒙NES颠LOV(腿NJ,N小,ND,晒AK,P袭)籍DIME闹NSIO理NP(罚500)薄,AK(杨500,巨100)旬NJ1=培N-1秧DO2筑0K=嗽1,NJ菊1再IF(N吼.GT.喘K+ND住-1)I划M=K+中ND-1半IF(N焰.LE.义K+ND拖-1)I素M=N锈K1=K吴+1够DO2体0I=鸟K1,I液M逗L=I-沟K+1惧C=AK地(K,L柔)/AK陪(K,1绝)化IW=N赶D-L+辞1蔽DO1猾0J=命1,IW朗M=J+内I-K猎10停AK仙(I,J淋)=AK柿(I,J扇)-C*醒AK(K垮,M)时20零P(单I)=P子(I)-穗C*P(本K)拴P(N)妹=P(N漠)/AK浴(N,1睡)演DO4吹0I1茄=1,N炉J1之I=N-胃I1竿IF(N关D.GT付.N-I域+1)J歌O=N-贸I+1渴IF(N帽D.LE封.N-I喷+1)J鸟O=ND弯DO3效0J=申2,JO顽K=J+疾I-1凉30锡P(衔I)=P室(I)-培AK(I适,J)*蓝P(K)脑40省P(郊I)=P哈(I)/胀AK(I票,1)日WRIT守E(8,若50)手50项FO线RMAT钱(/8X榜,5('秒**'毁),'R唐ESUL长TEO危FCA杆LCUL指TION兴',壁系&5(脑'**牢')//端1X,'帐NODA隆LDI映SPLA慌CEME铺NTS'截//3X何,'NO捆DE',蛛秩&5X必,'X-自DISP脆.',8外X,'Y辅-DIS雀P.')判DO6迅0I=触1,NJ羽60带WR圣ITE(徐8,70访)I,P圆(2*I扮-1),前P(2*残I)搬70掌FO袜RMAT罚(1X,扭I5,2筛E15.家6)葵RETU线RNEND赵SUBR没OUTI京NEM灯ADE(漫NE,N户J,N,数E,PR晓,LND年,X,Y帐,P)队DIME网NSIO私NLN驳D(50枕0,3)颂,X(2末00),堵Y(20辰0),D猴(3,3顶),季篇&B(租3,6)贵,S(3像,6),拥ST(3龙),P(衣500)壮,DE(纷6)习WRIT羊E(8岗,10)陶10柜FO跑RMAT详(/1X疗,'EL省EMEN载TST掀RESS化ES'/效)百CALL居MD(依E,PR朴,D)咬DO5招0IE摊=1,N伐E番CALL询MA(赞IE,N谱J,LN偏D,X,况Y,AE党)古CALL丙MB(拉IE,N棋J,LN跃D,X,风Y,AE痒,B)片DO2抢0I=剖1,3堤DO2凳0J=无1,6凑S(I,邀J)=0慨.0追DO2器0K=伯1,3海20鼠S(席I,J)祸=S(I拔,J)+谎D(I,滥K)*B救(K,J乌)救DO3词0I=福1,3梁DO3皂0J=未1,2铲IH=2严*(I-商1)+J廊IW=2孔*(LN晕D(IE凡,I)-也1)+J怠30摔DE灰(IH)他=P(I纤W)神DO3倦5I=旅1,3物35哑ST毅(I)=要0.0补DO4湿0I=扫1,3著DO4厅0J=昏1,6王40球ST蝶(I)=码ST(I弊)+S(竟I,J)宗*DE(米J)叉STX=渗ST(1闪)溜STY=德ST(2县)检TXY=僵ST(3覆)溉50固WR互ITE(助8,60归)IE,窗STX,充STY,容TXY间60嚼FO部RMAT宜(1X,拍'ELE燥MENT拥NO.腐=',I慎5/3X息,'ST粮X=',夏E15.猎6,捉才&2X折,'ST迷Y=',鄙E15.宫6,2X洲,'TX杨Y=',刘E15.惜6)伤RETU依RNEND部****帐****梳****筒****阴****郊****椅****知****秧****只****质****乳****蜜****造****米****进****留****死*摔输入文治件膛1.da软t炭:宫10,钳9,2,才1,0脖1.0,俱0.2焰5,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,步3.0,世6.0饶1,1织,1,端4,1,身1个10.0豪,0.挖0,1恶0.0静****晋****仆****余****秒****鲜****知****流****亏****主****昂****单****传****军****葛****批****欧****代*龙PLAT中EST求RESS凉PRO苍BLEM狡*缓***熄***归**姑*RES亭ULTE串OF语CALC撇ULTI红ON*池***厘***串**康*锦NOD袍ALD详ISPL之ACEM察ENTS狐N调ODE揪衫X-DI疯SP.短Y冠-DIS圾P.虾狂1.0订0000揉0E+0宋0孝.00弱0000豪E+00夹寇2.1架0941挖1E+0阴1懂.17象7565炼E+02蹄锅3新-.1效0941炮1E+0房1怒.17魂7565演E+02染宗4.0反0000猪0E+0稳0炮.00嘉0000虏E+00切壶5锣-.1阿6411达6E+0译1挽.15秤6785疯E+02莫剥6紫-.5盏7425贯6E-0背6茧.20汪9599治E+02表波7.1葛6411究6E+0献1吧.15坝6785蚕E+02案娃8.8推2058遗5E+0除0兽.25当3127录E+02俩割9陷-.8诱2057幅7E+0寒0壮.25骡3127鼻E+02乱族10.8厕4227歌4E-0逃5昆.44铅4730芹E+02悟ELE闷MENT妇STR界ESSE帆S锤ELE惩MENT毁NO.践=押1败S泼TX=宫.习1490惹23E+粪01笑STY=躬食.377铺270E煌+01蚊TXY积=虽.31湾1365族E+01待ELE毯MENT计NO.企=击2子S丈TX=歉-.袜7399激29E+霉00质STY=徐岔.141爬672E
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 福建省莆田市某校2024-2025学年高二上学期期末考试政治试题(原卷版+解析版)
- 2025年湖北省中考思想品德与法治模拟试卷(附答案)
- 2025年党章党纪党史党建知识竞赛多项选择题库及答案(共190道题)
- 中学精神文明建设工作计划
- 宫腔镜下输卵管插管通液治疗不孕不育效果探讨
- 节日教职工福利(花生油)项目 投标方案(技术方案)
- 乡村农田管理与开发协议
- 音乐制作与发行全流程指南
- 船舶导航与航行技术指南
- 环保设备可行性研究报告
- 东北三省三校2024年高三一模(第一次联合模拟考试)语文试卷(含答案)
- 无人机的传感器系统
- 图文解读中小学教育惩戒规则(试行)全文内容课件模板
- 2024年广西旅发置业集团有限公司招聘笔试参考题库含答案解析
- 《无尘室基础知识》课件
- 中式烹调技艺教案
- 人工智能引论智慧树知到课后章节答案2023年下浙江大学
- 加固工程监理实施细则
- 医保按病种分值付费(DIP)院内培训
- 高聚物改性沥青防水卷材检验报告
- 第二章-世界职业教育历史沿革
评论
0/150
提交评论