構造力学 構造工学特論

FORTRANプログラム

平板の曲げの計算データ例(剛性マトリックス法)
FEMDA5.DAT
3辺固定1辺単純支持 24要素 台形分布荷重

以下のプログラム
FEMPL2.FOR
に対応する。
C     READ A(I),B(I): REAL LENGTH OF PLATE 2*A(I),2*B(I)
C     弾性床上の平板 FEM
C     "FEMPL2.FOR"  台形分布荷重に対応 1990.9.29
C                   Y方向変化とX方向変化にそれぞれ対応 1991.7.22
C                    いちおうY方向変化に対応する

注)      NODT: 総節点数   NEPL: 総要素数     NF:集中荷重数
       KH = 1: 分布荷重   KH <> 1: 集中荷重     NCASE:荷重ケース
      READ(ICARD,10,END=90) NODT,NEPL,NF,KH,NCASE
   10 FORMAT(5I5)
   35   24    0    1    1
注)       KAKOM(I,J): 矩形要素の節点番号 A(I): 辺aの長さ   B(I):辺bの長さ
        BK(I):バネ強さ
      DO 100 I=1,NEPL
  100 READ(ICARD,20) (KAKOM(I,J),J=1,4),A(I),B(I),BK(I)
   20 FORMAT(4I5,3E10.3)
    1    2    8    9    0.6125     0.975   0.0E+00
    2    3    9   10    0.6125     0.975   0.0E+00
    3    4   10   11    0.6125     0.975   0.0E+00
    4    5   11   12    0.6125     0.975   0.0E+00
    5    6   12   13    0.6125     0.975   0.0E+00
    6    7   13   14    0.6125     0.975   0.0E+00
    8    9   15   16    0.6125     0.975   0.0E+00
    9   10   16   17    0.6125     0.975   0.0E+00
   10   11   17   18    0.6125     0.975   0.0E+00
   11   12   18   19    0.6125     0.975   0.0E+00
   12   13   19   20    0.6125     0.975   0.0E+00
   13   14   20   21    0.6125     0.975   0.0E+00
   15   16   22   23    0.6125     0.975   0.0E+00
   16   17   23   24    0.6125     0.975   0.0E+00
   17   18   24   25    0.6125     0.975   0.0E+00
   18   19   25   26    0.6125     0.975   0.0E+00
   19   20   26   27    0.6125     0.975   0.0E+00
   20   21   27   28    0.6125     0.975   0.0E+00
   22   23   29   30    0.6125     0.975   0.0E+00
   23   24   30   31    0.6125     0.975   0.0E+00
   24   25   31   32    0.6125     0.975   0.0E+00
   25   26   32   33    0.6125     0.975   0.0E+00
   26   27   33   34    0.6125     0.975   0.0E+00
   27   28   34   35    0.6125     0.975   0.0E+00
注)       L: 節点荷重 IZR: 垂直変位拘束あれば1   ISXR:x軸回り回転拘束あれば1
          ISYR:y軸回り回転拘束あれば1
      DO 200 I=1,NODT
      READ(ICARD,30) L,IZR,ISXR,ISYR
   30 FORMAT(I5,5X,3I5)
    1         1    1    1
    2         1    1    1
    3         1    1    1
    4         1    1    1
    5         1    1    1
    6         1    1    1
    7         1    1    1
    8         1    1    1
    9         0    0    0
   10         0    0    0
   11         0    0    0
   12         0    0    0
   13         0    0    0
   14         0    0    1
   15         1    1    1
   16         0    0    0
   17         0    0    0
   18         0    0    0
   19         0    0    0
   20         0    0    0
   21         0    0    1
   22         1    1    1
   23         0    0    0
   24         0    0    0
   25         0    0    0
   26         0    0    0
   27         0    0    0
   28         0    0    1
   29         1    1    1
   30         1    1    1
   31         1    1    1
   32         1    1    1
   33         1    1    1
   34         1    1    1
   35         1    1    1
注)      PTH:板厚  E:ヤング係数  G:せん断弾性係数
      READ(ICARD,40) PTH,E,G
   40 FORMAT(3E10.0)
   0.4E+00   2.1E+06   807692.
注)     NFT:矩形要素に作用する分布荷重の総数
      IF(KH-1) 330,320,330
  320 READ(ICARD,450) NFT
  450 FORMAT(I5)
   24
注)     NUM:要素番号   Q1:台形分布荷重q1 Q2:台形分布荷重q2
      DO 420 I=1,NFT
      READ(ICARD,460) NUM,Q1,Q2
  460 FORMAT(I5,5X,2F10.0)
    1      0.57      1.35
    2      1.35      2.13
    3      2.13      2.92
    4      2.92      3.70
    5      3.70      4.48
    6      4.48      5.26
    7      0.57      1.35
    8      1.35      2.13
    9      2.13      2.92
   10      2.92      3.70
   11      3.70      4.48
   12      4.48      5.26
   13      0.57      1.35
   14      1.35      2.13
   15      2.13      2.92
   16      2.92      3.70
   17      3.70      4.48
   18      4.48      5.26
   19      0.57      1.35
   20      1.35      2.13
   21      2.13      2.92
   22      2.92      3.70
   23      3.70      4.48
   24      4.48      5.26