Hirai Takayuki

8章 要素内の変位と応力度の計算

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座標を見て分かることになる。

 

3s_8_1

図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点

3s_8_2

図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 応力度の符号

 3s_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

PAGETOP
Copyright © 平居 孝之 All Rights Reserved.
Powered by WordPress & BizVektor Theme by Vektor,Inc. technology.