有限单元法程序报告.doc

目录
程序一:平面刚架静力分析程序() 1
程序二:平面三节点有限元程序 17
程序三:四节点矩形薄板单元程序 24
程序一:平面刚架静力分析程序()
已知各杆截面均为矩形,,, ,,各杆E=×104 MPa。
图2
节点、单元编号如下图3,…..为节点号,①②③……为单元号:
图3
总共有13个节点,13个单元。
计算源程序如下:
! (A program for analysis of plane frame)
! Version 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() 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() 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 & print them