構造力学　構造工学特論

FORTRANプログラム

平板（４辺単純支持、等分布荷重）

C     plate 4-simple-supported
C     uniform load
C     1998.7.2
C
      REAL E,H,PNU,A,B,X,Y,Q,D,PAI,W,MPY,MPX,ALM
      INTEGER K
C
      READ(5,100) E,H,PNU
  100 FORMAT(F10.0,2F5.0)
      READ(5,110) A,B,Q,X,Y
  110 FORMAT(5F5.0)
C
      K=10
      D=E*H**3/12.0/(1.0-PNU**2)
      PAI=3.14159
      W=0.0
C
      DO 10 M=1,K,2
         MPY=M*PAI*Y/A
         MPX=M*PAI*X/A
         ALM=M*PAI*B/2/A
         W=W+(1.0-(ALM*TANH(ALM)+2.0)/2/COSH(ALM)*COSH(MPY)
     &      +1.0/2.0/COSH(ALM)*MPY*SINH(MPY))*SIN(MPX)/M**5
   10 CONTINUE
C
      W=W*4.0*Q*A**4/PAI**5/D
C
      WRITE(6,120) D
  120 FORMAT(1H ,'D= ',F10.3)
      WRITE(6,130) A,B,X,Y,Q
  130 FORMAT(1H ,'a=',F8.3,'cm  b=',F8.3,'cm  x=',F8.3,
     &             'cm  y=',F8.3,'cm  q=',F8.3,'kg/cm2')
      WRITE(6,140) K
  140 FORMAT(1H ,'項数=',I3,'項まで')
      WRITE(6,*) 'w=',W,'cm'
C
      STOP
      END