杆件结构分析程序

杆件结构分析程序

ID:46127800

大小:159.40 KB

页数:36页

时间:2019-11-21

杆件结构分析程序_第1页
杆件结构分析程序_第2页
杆件结构分析程序_第3页
杆件结构分析程序_第4页
杆件结构分析程序_第5页
资源描述:

《杆件结构分析程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、programframeCPROGRAMMAININTEGERRR,AA,GG,EP,RP,TREALKKCHARACTER*12FNAME,FNAOTCOMMON/A1/NE,NP,NR,MT,MXCOMMON/A2/ME(2),RR(20,3),GG(6),TC0MM0N/A3/EP(20,3),RP(20,4),XZ(30,2),AE(5,3)C0MM0N/A4/KI2,KI3,MI2(20),MI3(20),FI2(20,3),FI3(20,6)C0MM0N/A5/AA(200),KK(6000),NX,NHC0MM0N/A6/FF(60),UV(20,6)WRITE(0

2、,101)101FORMAT(///5X,JPLEASEINPUTDATAFTLENAME/)READ(*,100)FNAME100FORMAT(8A)WRITE(0,102)102FORMAT(///5X,?PLEASEINPUTOUTPUTDATAFILENAME/)READ(*,103)FNAOT103FORMAT(8A)OPEN(5,FILE二FNAME,STATUS二'OLD')OPEN(6,F1LE=FNAOT,STATUS二'new')CALLINPUTCALLFORRRCALLFORKKCALLDIVKKDO200N9=1,MXCALLFORFF(N9)CALLFOR

3、DISCALLFORFMCALLFORRC200CONTINUESTOPENDSUBROUTINEINPUTINTEGEREP,RPC0MM0N/A1/NE,NP,NR,MT,MXC0MM0N/A3/EP(20,3),RP(20,4),XZ(30,2),AR(5,3)READ(5,*)NE,NP,NR,MT,MXWRITE(6,50)NE,NP,NR,MT,MXWRITE(0,50)NE,NP,NR,MT,MX50FORMAT(2X,'NUMBEROFELEMENTNE=,,T2/2X,#'NUMBEROFNODENP二',I2/2X,#'NUMBEROFSUPPOTNODENR二'

4、,I2/2X,#'NUMBEROFMATERIALSETMT二',12/2X,#'NUMBEROFWORKINGCASEMX二',12)READ(5,*)((AE(I,J),J二1,3),I二1,MT)WRITE(6,221)221FORMAT(2X,'NO.OFSET',2X,'MODULUSOFELASTICITY',2X,#'CROSS-AREA,,2X,'MOMEMTOFINERTIA,)DO15[二1,MTWRITE(6,222)I,(AE(I,J),J=l,3)222FORMAT(2X,14,1OX,3F16.5)15CONTINUEREAD(5,*)((EP(I,J),

5、J=l,3),I=l,NE)WRITE(6,229)229FORMAT(1OX,'DATAOFELEMENT,)DO171=1,NEWRITE(6,193)1,(EP(I,J),J=l,3)193FORMAT(2X,14,5X,318)17CONTINUEREAD(5,*)((RP(I,J),J二1,4),I二1,NR)WRITE(6,224)224FORMAT(1OX,'DATAOFSUPPORT')WRITE(6,96)((RP(I,J),J二1,4),1=1,NR)WRITE(0,96)((RP(I,J),J二1,4),I二1,NR)96FORMAT(2X,416)READ(5

6、,*)((XZ(I,J),J=l,2),1=1,NP)WRITE(6,226)226FORMAT(1OX,'NODALCOORDINATE,/5X,'NO.NODE',8X,'X-X'#,8X,'Z-Z')DO2271=1,NPWRITE(6,194)1,(XZ(I,J),J=l,2)194FORMAT(5X,14,5X,2F12.6)227CONTINUERETURNENDSUBROUTINEFORRRINTEGERRR,AA,GG,EP,RP,TREALKKCOMMON/A1/NE,NP,NR,MT,MXCOMMON/A2/ME(2),RR(20,3),GG(6),TCOMMON

7、/A3/EP(20,3),RP(20,4),XZ(30,2),AE(5,3)COMMON/A5/AA(200),KK(6000),NX,NHDO483T=1,NPDO483J二1,3RR(I,J)=l483CONTINUEDO4701=1,NRTZ1=RP(T,1)IW1=RP(I,2)IW2=RP(I,3)1W3=RP(1,4)IF(IWl.EQ.0)RR(IZl,1)=0IF(IW2.EQ.0)RR(IZl,2)=0IF(IW3.EQ.0)RR(IZ1

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

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

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