毕业论文：有限单元法程序的报告.doc

下载文档 降价啦

5
0
约3.48万字
约 36页
2017-12-21 发布于浙江
举报
版权申诉
保障服务

毕业论文：有限单元法程序的报告.doc

1、本文档共36页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

毕业论文：有限单元法程序的报告

目录程序一：平面刚架静力分析程序（PF.FOR） 1 程序二：平面三节点有限元程序 17 程序三：四节点矩形薄板单元程序 24 程序一：平面刚架静力分析程序（PF.FOR） …..为节点号，①②③……为单元号：图3 总共有13个节点，13个单元。计算源程序如下： ! PF.FOR (A program for analysis of plane frame) ! Version 6.3 2004 ! Main program reads the control data organizes the whole ! calculation by calling subroutines. DIMENSION W(80000) CHARACTER IDFN*20,TITLE(5)*72 READ(*,(A12))IDFN OPEN(3,FILE=IDFN,STATUS=OLD) READ(3,(A72))(TITLE(M),M=1,5) READ(3,*)E,NM,NJ,NS,NLC L1=1 L2=L1+NM L3=L2+NM L4=L3+NM L11=L4+NM L12=L11+NJ L21=L12+NJ L22=L21+NS L31=L22+NS L41=L31+6*NM CALL IOMJS(TITLE,E,NM,NJ,NS,NLC,W(L1),W(L2),W(L3), W(L4),W(L11),W(L12),W(L21),W(L22)) CALL LCVCT (NM,W(L1),W(L2),W(L31),NJ,N) CALL LCDIA (NM,N,W(L31),W(L41),W(L41),W(L41),MAXBDW,NA) L51=L41+N L52=L51+36 L53=L52+NA*2 L54=L53 L61=L54+N*2 NW=L61+6*NM-1 WRITE (*,1)NA,NW 1 FORMAT(/40X,( NA=,I6, ),/40X,( NW=,I6, )) CALL FORMA (E,NM,NJ,N,NA,W(L1),W(L2),W(L3),W(L4), W(L11),W(L12),W(L31),W(L51),W(L41),W(L52)) CALL AS(NS,N,NA,W(L21),W(L41),W(L52)) CALL LDLT(N,NA,W(L41),W(L52),W(L53)) DO 100 LC=1,NLC READ (3,*)NLJ L62=L61+NLJ L63=L62+NLJ L64=L63+NLJ L71=L61 L81=L71+6*NM CALL B0(LC,N,NLJ,W(L54)) IF(NLJ.NE.0) CALL IOLJB(N,NLJ,W(L61),W(L62), W(L63),W(L64),W(L54)) READ(3,*)NLM L82=L81+NLM L83=L82+NLM L84=L83+NLM CALL F0(NLM,NM,W(L71)) IF(NLM.NE.0) CALL IOLMFB(NM,NJ,N,NLM,W(L81), W(L82),W(L83),W(L84),W(L1),W(L2),W(L11), W(L12),W(L31),W(L71),W(L54)) CALL BS(NS,N,W(L21),W(L22),W(L54)) CALL SLVEQ(N,NA,MAXBDW,W(L41),W(L52),W(L54)) CALL OJD(NJ,N,W(L54)) CALL COTF(E,NM,NJ,N,W(L1),W(L2),W(L3),W(L4), W(L11),W(L12),W(L31),W(L54),W(L71)) NW=L84+NLM-1 WRITE(*,1)NA,NW 100 CONTINUE WRITE(*,(/)) END SUBROUTINE IOMJS(TITLE,E,NM,NJ,NS,NLC,IST,IEN, AR,RI,X,Y,IS,VS) ! Read data of members, joints, supports