構造力学 構造工学特論
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 いちおうＹ方向変化に対応する
注） 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): 辺ａの長さ B(I):辺ｂの長さ
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: 垂直変位拘束あれば１ ISXR:ｘ軸回り回転拘束あれば１
ISYR:ｙ軸回り回転拘束あれば１
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:台形分布荷重ｑ１ 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