8-1 変位と応力度を計算する位置
要素の20個の節点が図4-3に示した位置関係を保ったまま、図8-1のように、稜の長さが2の立方体の要素に置き換える。置き換えた立方体の座標軸をξ(クサイ)、η(イ-タ)、ζ(ツエ-タ)で表す。このとき要素の20個の節点のξ、η、ζ座標は次のようになる。
節点番号 N31 N32 N33 N34 N35 N36 N37 N38 N39 N40 N41 N42 N43 N44 N45 N46 N47 N48 N49 N50
ξ -1 1 1 -1 -1 1 1 -1 0 1 0 -1 -1 1 1 -1 0 1 0 -1
η -1 -1 1 1 -1 -1 1 1 -1 0 1 0 -1 -1 1 1 -1 0 1 0
ζ -1 -1 -1 -1 1 1 1 1 -1 -1 -1 -1 0 0 0 0 1 1 1 1
図8-1に示したξ、η、ζ座標により、要素内で変位と応力度を計算する位置を図8-2のイ)~ニ)で指定する。イ)~ニ)で指定した位置が、本来の形状のX、Y、Z座標で表すとどうなるかは、図8-1の節点の位置関係から推測できるが、正確な座標は計算結果に出るX、Y、Z座標を見て分かることになる。
図8-1 要素のξηζ座標表示
イ)3×3×3の格子点
ξ、η、ζがそれぞれ-1、0、1になる組み合わせの27点。
ロ)1つの面上の3×3の格子点
ξ、η、ζのいずれかの軸に垂直な面上で、ξ、η、ζのうちの2つが-1、0、1になる組み合わせの9点。
ハ)1つの線上の5点
ξ、η、ζのいずれかの軸に平行な線上で、ξ(またはηかζ)が-1、-0.5、 0、 0.5、 1になる5点。
ニ)任意の一点
ξ、η、ζが任意の値の1点。ただし、-1≦ξ、η、ζ≦1 でないと計算結果は意味を持たないことに注意する。
イ)3×3×3の格子点 ロ)1つの面上の9点 ハ)1つの線上の5点 ニ)任意の1点
図8-2 変位と応力度を計算する位置指定の種類
8-2 位置指定デ-タの作成方法
要素内の変位と応力度を計算する位置を指定するデ-タは、表8-1のような数値の並びで作成する。表8-2は、図4-2に示したモデル401における位置指定デ-タの例である。
N1 整数5字枠、枠の右につめて入れる。要素の番号。
N2 整数5字枠、枠の右につめて入れる。N2の数値で次のように位置を指定する。
1~2000のとき、N1からN2までの要素について図8-2のイ)の27点。
7100のとき、N1の要素についてξ=C1の面上で図8-2のロ)の9点。
7010のとき、N1の要素についてη=C2の面上で図8-2のロ)の9点。
7001のとき、N1の要素についてξ=C3の面上で図8-2のロ)の9点。
8011のとき、N1の要素についてη=C2,ζ=C3の線上で図8-2のハ)の5点。
8101のとき、N1の要素についてξ=C1,ζ=C3の線上で図8-2のハ)の5点。
8110のとき、N1の要素についてξ=C1,η=C2の線上で図8-2のハ)の5点。
9000のとき、N1の要素についてξ=C1、η=C2、ζ=C3の1点。
C1 小数10字枠、小数点付きで入れる。ξの値。
C2 小数10字枠、小数点付きで入れる。ηの値。
C3 小数10字枠、小数点付きで入れる。ζの値。
99999 デ-タの終わりを表すために、最後の行の先頭から99999を入れる。
表8-1 位置指定デ-タの並び
(注:最上行の数字は、数値配列の欄の数で、入力データではない)
5 10 20 30 40
N1 N2 C1 C2 C3
〃 〃 〃 〃 〃
〃 〃 〃 〃 〃
〃 〃 〃 〃 〃
99999
表8-2 モデル401における位置指定デ-タの例
(ファイルD1M401に収納)
2 2 0.0 0.0 0.0
4 7001 0.0 0.0 0.5
5 8011 0.0 -0.5 0.5
1 9000 0.25 -0.25 0.5
2 9000 0.75 0.0 -0.5
99999
8-3 位置指定デ-タの数値の並びの条件
表8-1の位置指定デ-タの数値の並びは必ず次のようにする。
表8-1のN1,N2 5字枠の右端につめて入れる。
表8-1のC1,C2、C3 小数点を付けて10字枠内に入れる。
表8-1の99999 最後の行に先頭から99999を入れる。
8-4 要素内の変位と応力度の計算方法
以下の1)~8)の手順で行う。なおパソコンの操作は、2章を参照。
1)パソコンをWindowsで起動し、ハ-ドディスク上に自分の作業用に作成したフォルダ(ここでは hirai )の窓を開く。
2)フォルダ hirai に節点の変位と節点力の計算結果を入れたファイルをコピ-し、そのファイルの名前をR1とする。
3)フォルダ hirai に8-2節に説明した変位と応力度を計算する位置指定のデ-タを入れたファイルをコピ-し、そのファイルの名前をD1とする。
4)フォルダ hirai の窓の中の F3DSD.EXE のアイコンをダブルクリックする。画面に新しい窓が開き、計算の実行の進行程度の表示がスクリーン上に出て、短時間で計算が終了し窓が消える。
5)計算結果はディスクのファイルSDに収納されている。
表8-3 応力度の符号
8-5 要素内の変位と応力度の計算結果の見方
表6-2に示したモデル401の節点の変位と節点力の計算結果について、表8-2に示したように位置を指定して変位と応力度を計算した結果は、表8-4のような数値の並びでファイルSDに入っている。
表8-4のように、最初にモデルのタイトルと要素数、節点数、材質の数が出た後、位置を指定した要素ごとに結果が出る。要素ごとには、要素番号、図8-2の位置指定の種類とξ(GZ)η(EA)ζ(ZE)の値が出た後、要素内に指定した各点ごとに次の値が出る。
X座標、Y座標、Z座標、X方向変位、Y方向変位、Z方向変位、X方向垂直応力度、Y方向垂直応力度、Z方向垂直応力度、XYせん断応力度、YZせん断応力度、ZXせん断応力度
変位の数値は、例えば 0.778787e-05 は 0.778787×10-5 (=0.00000778787)である。また、応力度は表8-3に示すような符号で表されている。
表8-4 要素内の変位と応力度の計算結果
MODEL TITLE = 401
ELEMENTS = 6 NODES = 80 MATERIALS = 2
ELEMENT = 2
COODINATE X-Y-Z DISPLACEMENT X-Y-Z STRESS X-Y-Z-XY-YZ-ZX
0.00000 0.00000 150.00000 0.11574e-04 -0.18483e-02 -0.53827e-03 0.053 -0.002 0.903 -0.036 -0.629 0.015
0.00000 0.00000 210.00000 0.40483e-04 -0.23431e-02 -0.16330e-03 -0.006 0.285 1.913 -0.031 0.036 0.018
0.00000 0.00000 270.00000 0.39649e-04 -0.15697e-02 0.47478e-03 0.105 -0.112 2.821 0.054 -0.055 0.022
0.00000 30.00000 150.00000 -0.41635e-05 -0.18634e-02 -0.17101e-03 -0.028 -0.089 -0.282 -0.028 -0.482 0.007
0.00000 30.00000 210.00000 -0.49541e-05 -0.23442e-02 -0.23355e-03 0.027 0.048 -0.185 -0.048 -0.091 -0.011
0.00000 30.00000 270.00000 -0.92522e-05 -0.16520e-02 -0.28021e-03 -0.068 -0.578 -0.271 0.013 -0.455 -0.028
0.00000 60.00000 150.00000 -0.21841e-04 -0.18557e-02 0.21343e-03 -0.129 -0.213 -1.631 -0.019 -0.334 0.002
0.00000 60.00000 210.00000 -0.50545e-04 -0.23212e-02 -0.32852e-03 0.019 -0.230 -2.453 -0.064 -0.216 -0.037
0.00000 60.00000 270.00000 -0.56519e-04 -0.17088e-02 -0.11018e-02 -0.304 -1.091 -3.539 -0.029 -0.853 -0.076
20.00000 0.00000 150.00000 0.40353e-14 -0.18473e-02 -0.54395e-03 0.061 0.016 0.928 0.000 -0.593 0.000
20.00000 0.00000 210.00000 0.85902e-14 -0.23313e-02 -0.16366e-03 -0.004 0.276 1.931 0.000 0.051 0.000
20.00000 0.00000 270.00000 0.13145e-13 -0.15472e-02 0.47975e-03 0.101 -0.148 2.832 0.000 -0.060 0.000
20.00000 30.00000 150.00000 0.38193e-14 -0.18608e-02 -0.17041e-03 -0.027 -0.078 -0.284 0.000 -0.452 0.000
20.00000 30.00000 210.00000 0.82859e-14 -0.23342e-02 -0.23430e-03 0.022 0.033 -0.194 0.000 -0.081 0.000
20.00000 30.00000 270.00000 0.12752e-13 -0.16346e-02 -0.28230e-03 -0.080 -0.621 -0.287 0.000 -0.465 0.000
20.00000 60.00000 150.00000 0.36034e-14 -0.18516e-02 0.22030e-03 -0.134 -0.209 -1.661 0.000 -0.308 0.000
20.00000 60.00000 210.00000 0.79816e-14 -0.23130e-02 -0.32966e-03 0.008 -0.253 -2.490 0.000 -0.211 0.000
20.00000 60.00000 270.00000 0.12360e-13 -0.16965e-02 -0.11110e-02 -0.323 -1.141 -3.582 0.000 -0.868 0.000
40.00000 0.00000 150.00000 -0.11574e-04 -0.18483e-02 -0.53827e-03 0.053 -0.002 0.903 0.036 -0.629 -0.015
40.00000 0.00000 210.00000 -0.40483e-04 -0.23431e-02 -0.16330e-03 -0.006 0.285 1.913 0.031 0.036 -0.018
40.00000 0.00000 270.00000 -0.39649e-04 -0.15697e-02 0.47478e-03 0.105 -0.112 2.821 -0.054 -0.055 -0.022
40.00000 30.00000 150.00000 0.41635e-05 -0.18634e-02 -0.17101e-03 -0.028 -0.089 -0.282 0.028 -0.482 -0.007
40.00000 30.00000 210.00000 0.49541e-05 -0.23442e-02 -0.23355e-03 0.027 0.048 -0.185 0.048 -0.091 0.011
40.00000 30.00000 270.00000 0.92522e-05 -0.16520e-02 -0.28021e-03 -0.068 -0.578 -0.271 -0.013 -0.455 0.028
40.00000 60.00000 150.00000 0.21841e-04 -0.18557e-02 0.21343e-03 -0.129 -0.213 -1.631 0.019 -0.334 -0.002
40.00000 60.00000 210.00000 0.50545e-04 -0.23212e-02 -0.32852e-03 0.019 -0.230 -2.453 0.064 -0.216 0.037
40.00000 60.00000 270.00000 0.56519e-04 -0.17088e-02 -0.11018e-02 -0.304 -1.091 -3.539 0.029 -0.853 0.076
ELEMENT = 4 APPOINTMENT = 7001 ZE = 0.50000
COODINATE X-Y-Z DISPLACEMENT X-Y-Z STRESS X-Y-Z-XY-YZ-ZX
0.00000 60.00000 330.00000 0.78040e-05 0.55327e-03 -0.12998e-02 0.008 0.846 -0.357 0.010 0.155 -0.016
0.00000 110.00000 332.50000 0.14275e-04 0.84130e-03 -0.35487e-02 -0.001 0.710 0.108 -0.008 -0.052 0.015
0.00000 160.00000 335.00000 0.24944e-05 0.11044e-02 -0.63933e-02 0.020 0.434 -0.190 -0.023 -0.310 -0.018
20.00000 60.00000 330.00000 0.18813e-13 0.55153e-03 -0.13092e-02 0.011 0.859 -0.356 0.000 0.135 0.000
20.00000 110.00000 332.50000 0.18623e-13 0.84146e-03 -0.35568e-02 0.003 0.721 0.115 0.000 -0.055 0.000
20.00000 160.00000 335.00000 0.18432e-13 0.11065e-02 -0.63999e-02 0.025 0.445 -0.177 0.000 -0.291 0.000
40.00000 60.00000 330.00000 -0.78040e-05 0.55327e-03 -0.12998e-02 0.008 0.846 -0.357 -0.010 0.155 0.016
40.00000 110.00000 332.50000 -0.14275e-04 0.84130e-03 -0.35487e-02 -0.001 0.710 0.108 0.008 -0.052 -0.015
40.00000 160.00000 335.00000 -0.24944e-05 0.11044e-02 -0.63933e-02 0.020 0.434 -0.190 0.023 -0.310 0.018
ELEMENT = 5 APPOINTMENT = 8011 EA = -0.50000 ZE = 0.50000
COODINATE X-Y-Z DISPLACEMENT X-Y-Z STRESS X-Y-Z-XY-YZ-ZX
0.00000 190.00000 335.00000 -0.84801e-06 0.11208e-02 -0.83476e-02 -0.020 -0.047 -0.101 0.007 -0.220 -0.006
10.00000 190.00000 335.00000 -0.42401e-06 0.11223e-02 -0.83494e-02 -0.018 -0.045 -0.096 0.004 -0.198 -0.003
20.00000 190.00000 335.00000 0.18137e-13 0.11229e-02 -0.83501e-02 -0.018 -0.044 -0.094 0.000 -0.191 0.000
30.00000 190.00000 335.00000 0.42401e-06 0.11223e-02 -0.83494e-02 -0.018 -0.045 -0.096 -0.004 -0.198 0.003
40.00000 190.00000 335.00000 0.84801e-06 0.11208e-02 -0.83476e-02 -0.020 -0.047 -0.101 -0.007 -0.220 0.006
ELEMENT COODINATE X-Y-Z DIPLACEMENT X-Y-Z STRESS X-Y-Z-XY-YZ-ZX ( GZ-EA-ZE = 0.25000 -0.25000 0.50000 )
1 25.00000 22.50000 112.50000 0.16899e-05 -0.12185e-02 -0.23029e-03 0.012 -0.040 -0.272 -0.009 -0.171 -0.001
ELEMENT COODINATE X-Y-Z DIPLACEMENT X-Y-Z STRESS X-Y-Z-XY-YZ-ZX ( GZ-EA-ZE = 0.75000 0.00000 -0.50000 )
2 35.00000 30.00000 180.00000 0.30902e-05 -0.22477e-02 -0.20430e-03 0.018 0.074 -0.213 0.036 -0.183 0.001