第八讲平面钢架程序设计知识讲解.pptVIP

下载本文档

2
0
约8.94千字
约 27页
2019-11-04 发布于天津
举报
版权申诉

第八讲平面钢架程序设计知识讲解.ppt

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

第八讲平面钢架程序设计; 第一节概述;计算模型与计算方法;在局部坐标系下，单元刚度矩阵如下式所示，即: ;3.支承条件的引入及整体刚度矩阵的组集;;;;第二节;EAI;RT;; 平面刚架静力分析原程序（PFSAP.FOR) 1: C ANALYSIS PROGRAM FOR PLANE FRAME 2: REAL K(200,200), KE(6,6),AKE(6,6),X(100),Y(100),AL(100) 3: EAI(3,100),PJ(100),PF(2,100) ,R(6,6),P(100),FF(6), 4: FE(6),D(100),ADE(6),DE(6),RT(6,6),AFE(6),F(3) 5: INTEGER JE(2,100),JN(3,100),JPJ(100),JPF(2,00),M(6), 6: JEAI(100),NO 7: OPEN(6,FILE=’PFSAP.IN’) 8: OPEN(8,FILE=’PFSAP.OUT’) 9: READ(6,·)NO 10: IF(NO.EQ.0)STOP 11: WRITE(8,(/A5,I3,A1)’)’(NO.=’,NO,’ 12: CALL READ(NJ,N,NE,NM,NPJ,NPF,JN,X,Y,JE, JEAI,EAI,JPJ,PJ,JPF,PF) 13: DO 5 I=1,N 14: P(I)=0 15: DO 5 J=1,N 16: K(I,J)=0 17: DO 10 IE=1,NE 18: CALL MKE(KE,IE,JE,JEAI,EAI,X,Y,AL) 19: CALL MR(R,IE,JE,X,Y) 20: CALL MAKE(KE,R,AKE);21: CALL CALM(M,IE,JN,JE 22: CALL MK(K,AKE,M) 23: 10 CONTINUE 24: DO 20 IP=1,NP 25: CALL MR(R,JPF(1,IP),JE,X,Y） 26: CALL TRAN(R,RT） 27: CALL PE(FE,IP,JPF,PF,AL） 28: CALL MULV6(RT,FE,AFE) 29: CALL CALM(M,JPF(1,IP),JN,JE) 30: CALL MF(P,AFE,M） 31:20 CONTINUE 32: DO 30 I=1,NP 33:30 P(JPJ(I))=P(JPJ(I))+PJ(I) 34: CALL SLOV(K,P,D,N） 35: WRITE(8,’(/2(26(1H·),A))’)’RESULTS OF CALCULATION’ 36: WRITE(8,40) 37:40 FORMAT(/5X,’NO.N’,4X,’X-DISPLACEMENT’,2X, 38: ‘Y-DISPLACEMENT’,3X,’ANG.ROT.(RAD)’) 39: DO 60 KK=1,NJ 40: DO 50 II=1,3 41: F(II)=0 42: I1=JN(II,KK);43: 50 IF(I1.GT.0)F(II)=D(I1) 44: 60 WRITE(8,70)KK,F(1),F(2),F(3) 45: 70 FORMAT(I8,2X,3G16.5) 46: WRITE(8,80) 47: 80 FORMAT(/’NO.E’,5X,’N(1),8X,’Q(1),8X,M(1)’, 48: 8X,’N(2)’,8X,’Q(2),8X,’M(2)’) 49; DO 130 IE=1,NE 50: CALL MADE(IE,JN,JE,D,ADE) 51: CALL MKE(KE,IE,JE,JEAI,EAI,X,Y,AL) 52: CALL MR(R,IE,JE,X,Y) 53: CALL MULV6(R,ADE,DE) 54: CALL MULV6(KE,DE,EF) 55: DO 100 IP=1,NPF 56: IF(JPF(1,IP).EQ.IE) THEN 57: CALL PE(EF,IP,JPF,PF,AL) 58: DO