结构力学程序方法.doc

下载文档 降价啦

39
0
约6.69千字
约 11页
2016-12-31 发布于重庆
举报
版权申诉
保障服务

结构力学程序方法.doc

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

MODULE DATATYPE IMPLICIT NONE INTEGER(KIND(1)),PARAMETER::IKIND=KIND(1),RKIND=KIND(0.0D0) REAL(RKIND),PARAMETER::ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0,FOUR=4.D0,FIVE=5.D0, SIX=6.D0,SEVEN=7.D0,EIGHT=8.D0,NINE=9.D0,TEN=10.D0 TYPE PJ_TYPE INTEGER(IKIND)::LOCNo REAL(RKIND)::VALUE END TYPE TYPE PF_TYPE INTEGER(IKIND)::NO INTEGER(IKIND)::TYP REAL(RKIND)::LOC REAL(RKIND)::VALUE END TYPE END MODULE DATATYPE MODULE SUBROUTINES USE DATATYPE PRIVATE PUBLIC::EJC,SCL,ESM,CTM,EFX CONTAINS SUBROUTINE EJC(M,JE,JN,JC) !子程序：形成单元定位向量 IMPLICIT NONE INTEGER (IKIND),INTENT(IN)::M,JE(:,:),JN(:,:) INTEGER (IKIND),INTENT(OUT)::JC(:) INTEGER (IKIND)::I,J1,J2 J1=JE(1,M) J2=JE(2,M) DO I=1,3 JC(I)=JN(I,J1) JC(I+3)=JN(I,J2) END DO RETURN END SUBROUTINE EJC SUBROUTINE SCL(M,BL,SI,CO,JE,X,Y) !子程序:求单元常数 IMPLICIT NONE INTEGER(IKIND),INTENT(IN)::M,JE(:,:) REAL(RKIND),INTENT(IN OUT)::X(:),Y(:) REAL(RKIND),INTENT(OUT)::BL,SI,CO INTEGER(IKIND)::J1,J2 REAL(RKIND)::DX,DY J1=JE(1,M) J2=JE(2,M) DX=X(J2)-X(J1) DY=Y(J2)-Y(J1) BL=SQRT(DX*DX+DY*DY) SI=DY/BL CO=DX/BL RETURN END SUBROUTINE SCL SUBROUTINE ESM(M,BL,EA,EI,KD) !子程序:单元刚度矩阵 IMPLICIT NONE INTEGER(IKIND),INTENT(IN) ::M REAL(RKIND),INTENT(IN OUT) ::BL,EA(:),EI(:) REAL(RKIND),INTENT(OUT) ::KD(:,:) INTEGER(IKIND)::I,J,I1 REAL(RKIND) ::G,G1,G2,G3 G=EA(M)/BL G1=TWO*EI(M)/BL G2=THREE*G1/BL G3=TWO*G2/BL KD=ZERO KD(1,1)=G KD(1,4)=-G KD(4,4)=G KD(2,2)=G3 KD(5,5)=G3 KD(2,5)=-G3 KD(2,3)=G2 KD(2,6)=G2 KD(3,5)=-G2 KD(5,6)=-G2 KD(3,3)=TWO*G1 KD(6,6)=TWO*G1 KD(3,6)=G1 DO I=1,5 I1=I+1 DO J=I1,6 KD(J,I)=KD(I,J) END DO END DO RETURN END SUBROUTINE ESM SUBROUTINE CTM(SI,CO,T) !子程序:形成单元坐标转换矩阵 IMPLICIT NONE REAL(RKIND),INTENT(IN OUT)::SI,CO REAL(RKIND),INTENT(OUT)::T(:,:) INTEGER(IKIND)::I,J T=ZERO