有限元分析程序与算例
10 17 18 8.0000 .5000 11 0 19 .0000 .0000 12 20 21 2.0000 .0000 13 22 23 4.0000 .0000 14 24 25 6.0000 .0000 15 26 0 8.0000 .0000
MATERIAL PROPERTIES
N.M. YOUNGS MODULUS POISION RATIO UNIT WEIGHT WIDTH 1 .2000E+08 .2000E+00 .0000E+00 .1000E+01
PARAMETERES OF WATER ANDSILT PRESSURE
N.P. ZERO-PRESSURESURFACE UNIT WEIGHT 1 11.00000 1.00000
NEE= NME= NET= NK= NSF= NST= NOD=
1 1 0 1 0 0 6 7 2 1 SURFACE NEWS NNN= 1 0010
NEE= NME= NET= NK= NSF= NST= NOD=
2 1 0 1 0 0 7 8 3 2 SURFACE NEWS NNN= 1 0010
NEE= NME= NET= NK= NSF= NST= NOD=
3 1 0 1 0 0 8 9 4 3 SURFACE NEWS NNN= 1 0010
NEE= NME= NET= NK= NSF= NST= NOD=
4 1 0 1 0 0 9 10 5 4 SURFACE NEWS NNN= 1 0010
NEE= NME= NET= NK= NSF= NST= NOD=
5 1 0 0 0 0 11 12 7 6 NEE= NME= NET= NK= NSF= NST= NOD=
6 1 0 0 0 0 12 13 8 7 NEE= NME= NET= NK= NSF= NST= NOD=
7 1 0 0 0 0 13 14 9 8 NEE= NME= NET= NK= NSF= NST= NOD=
8 1 0 0 0 0 14 15 10 9 TOTAL DEGREES OF FREEDOM.…………………N= 26 MAX-SEMI-BANDWIDTH ………………………… MX= 13 TOTAL-STORAGE …………………………………… NH= 230
LOAD OF NODES
NODE NO. X-DIRECTION Y-DIRECTION
1 .0000E+00 -.1000E+02 2 .0000E+00 -.2000E+02
- 21 -
有限元分析程序与算例
3 .0000E+00 -.2000E+02 4 .0000E+00 -.2000E+02 5 .0000E+00 -.1000E+02 6 .0000E+00 .0000E+00 7 .0000E+00 .0000E+00 8 .0000E+00 .0000E+00 9 .0000E+00 .0000E+00 10 .0000E+00 .0000E+00 11 .0000E+00 .0000E+00 12 .0000E+00 .0000E+00 13 .0000E+00 .0000E+00 14 .0000E+00 .0000E+00 15 .0000E+00 .0000E+00
DISPLACEMENT OF NODES
NODE NO. X-DIRECTION Y-DIRECTION
1 .0000E+00 -.1905E-02 2 -.6940E-04 -.1764E-02 3 -.1299E-03 -.1358E-02 4 -.1723E-03 -.7407E-03 5 -.1875E-03 -.3574E-05 6 .0000E+00 -.1907E-02 7 .8459E-07 -.1766E-02 8 .2524E-06 -.1359E-02 9 .1587E-06 -.7414E-03 10 .6309E-06 -.2797E-05 11 .0000E+00 -.1905E-02 12 .6961E-04 -.1764E-02 13 .1302E-03 -.1358E-02 14 .1728E-03 -.7412E-03 15 .1894E-03 .0000E+00
ELEMENT MID-POINT STRESS
ELE NO.SIG--XX SIG--YY SIG--XY SIG-- 1 SIG--2 ANGLE X--CO. Y--CO.
1 -348.24 -8.26 10.23 -7.95 -348.55 88.28 1.00 .75 2 -303.17 -8.66 29.51 -5.73 -306.10 84.33 3.00 .75 3 -214.19 -6.33 51.28 5.63 -226.15 76.87 5.00 .75 4 -76.58 -16.76 69.22 28.73 -122.08 56.68 7.00 .75 5 348.24 -1.23 9.78 348.52 -1.50 1.60 1.00 .25 6 303.17 -4.26 30.59 306.18 -7.27 5.63 3.00 .25 7 214.19 8.88 48.82 225.21 -2.14 12.72 5.00 .25
- 22 -
有限元分析程序与算例
8 76.58 -43.42 70.89 109.45 -76.29 24.88 7.00 .25
NODES & PRINCIPAL STRESS
NODE NO. SIG-XX SIG-YY SIG-XY SIG-1 SIG-2 ANGLE
1 -348.242 -8.255 10.233 -7.948 -348.550 88.278 2 -325.706 -8.459 19.869 -7.219 -326.946 86.430 3 -258.680 -7.497 40.392 -1.161 -265.016 81.086 4 -145.387 5 -76.584 6 .000 7 .000 8 -.001 9 -.001 10 .000 11 348.243 12 325.706 13 258.678 14 145.386 15 76.585
-11.548 -16.764 -4.742 -5.601 -2.593 -14.408 -30.090 -1.228 -2.744 2.311 -17.268 -43.417 60.249 69.221 10.004 20.026 40.049 60.051 70.053 9.775 20.183 39.707 59.854 70.886 - 23 -
11.578 28.732 7.910 17.421 38.773 53.278 56.606 348.516 326.942 264.687 165.037 109.454 -168.513 -122.080 -12.652 -23.022 -41.367 -67.686 -86.696 -1.501 -3.979 -3.698 -36.919 -76.286 69.001 56.685 38.334 41.020 44.073 41.580 38.939 1.601 3.503 8.606 18.176 24.877
有限元分析程序与算例
平面刚架静力分析程序
一、 源程序：
! PF.F90(A program form analysis of plane frmae) ! Version 4.3 2007 ! Main program reads the control date & organizes the ! whole calculation by calling subroutines
IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION W(20000) CHARACTER IDFN*12,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
- 1 -
有限元分析程序与算例
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
! Read data of members,joints,supports & print them SUBROUTINE IOMJS (TITLE,E,NM,NJ,NS,NLC,IST,IEN,AR,RI,X,Y,IS,VS) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION IST(NM),IEN(NM),AR(NM),RI(NM),X(NJ),Y(NJ),IS(NS),VS(NS) CHARACTER TITLE(5)*72 WRITE (*,'(/)') WRITE (*,1)(TITLE(M),M=1,5) 1 FORMAT(1X,A72) WRITE (*,2)E,NM,NJ,NS,NLC
2 FORMAT(/13X,'The Input Data'//10X,'The General Information'& //6X,'E',9X,'NM',5X,'NJ',5X,'NS',5X,'NLC' /1X,1PE10.3,4I7) READ (3,*)(IST(M),IEN(M),AR(M),RI(M),M=1,NM) WRITE (*,3)
3 FORMAT(/10X,'The Information of Members'//1X,'member',2X,'start',& 2X,'end',9X,'A',15X,'I') WRITE (*,4)(M,IST(M),IEN(M),AR(M),RI(M),M=1,NM) 4 FORMAT(1X,I4,I8,I6,1P2E16.6) READ (3,*)(X(M),Y(M),M=1,NJ)
- 2 -