構造力学　構造工学特論

FORTRANプログラム

平板の曲げの計算データ例（剛性マトリックス法）
FEMPL2.FOR

C     READ A(I),B(I): REAL LENGTH OF PLATE 2*A(I),2*B(I)
C     弾性床上の平板　ＦＥＭ
C     "FEMPL2.FOR"  台形分布荷重に対応　1990.9.29
C                   Ｙ方向変化とＸ方向変化にそれぞれ対応　1991.7.22
C                    いちおうＹ方向変化に対応する
C     1979,11,5 MIYAMOTO, DANSEI SYOJOO NO PLATE
　　−　　以下に主な入力データの説明　　−

C       NODT: 全節点数　　　NEPL: 全要素数
C       NF: 全節点荷重数
C       KH=1: 分布荷重　　　KH<>1: 集中荷重
      READ(ICARD,10,END=90) NODT,NEPL,NF,KH,NCASE
   10 FORMAT(5I5)

C       KAKOM(I,J): 矩形要素の節点番号　A(I): 辺ａの長さ   B(I):辺ｂの長さ
C       　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)

C       L: 各節点
C       IZR: Z方向拘束　　　ISXR: X方向回りの拘束　　　ISYR: Y方向回りの拘束
C       1: 拘束あり　　　　0: 拘束なし
      DO 200 I=1,NODT
      READ(ICARD,30) L,IZR,ISXR,ISYR
   30 FORMAT(I5,5X,3I5)

C       PTH: 板厚
C       E: ヤング係数　　　G: せん断弾性係数
      READ(ICARD,40) PTH,E,G
   40 FORMAT(3E10.0)

C       NUM: 各要素
C       Q1: 台形分布荷重q1　　　Q2: 台形分布荷重q2
      DO 420 I=1,NFT
      READ(ICARD,460) NUM,Q1,Q2
  460 FORMAT(I5,5X,2F10.0)

C       NO: 各節点
C       FZ: Z方向荷重　　　FMX: X方向回り荷重　　　FMY: Y方向回り荷重
  330 DO 80 I=1,NF
   80 READ(ICARD,95) NO(I),FZ(I),FMX(I),FMY(I)
   95 FORMAT(I5,5X,3F10.0)

このプログラムに関心のある方は宮本まで連絡してください。