16.1 シェルの解析-シェルの組み合わせ

トップページは<a href="https://qiita.com/Sagittarius_Chiron/items/fd806361ceca4864bc79">こちら</a>です。
1.0、2.0、4.0の記事の読了を前提として、操作についての記載を省略している場合があります。
<a href="https://qiita.com/Sagittarius_Chiron/items/35e5919069c0ab9db806">16.0</a> の読了も前提にしています。

###解析内容
<dl>
  <dt>使用ファイル</dt>
  <dd>shell_3part.stp</dd>
  <dt>単位</dt>
  <dd>mm-ton</dd>
  <dt>材料物性値:Top</dt>
  <dd>Steel　　(注1)</dd>
  <dd>縦弾性係数：E=210000 [MPa]</dd>
  <dd>ポアソン比：ν=0.3</dd>
</dl>
荷重状態
　屋根(面)： p = 1e-6 [MPa]　青色
　壁の下(エッジ)； Dx = Dy = Dz = 0  [mm]　緑色
形状
　屋根：2000×2000×1 [mm]
　壁　：1000×2×500 [mm]
![Fig16_10_01.PNG](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/888272/ff280dd8-57b7-786f-e7b8-afb67378036b.png)

<ul>
<li>(注1)</li>この値は元の例題での値です。ここでは PrePoMax の Material Library に登録済の近い値の材料を使います。
</ul>

###形状の読込
PrePoMaxの新規起動または新規作業開始をしてください。
unit system:<b>mm,ton,s,℃</b> に設定してください。
<dd>元の例題の記載では、単位系が m-kg-Pa になっています。しかし、Stepファイルの単位はmmでした。</dd>
shell_3part.stpをインポートしてください。
オブジェクト・ブラウザのGeometryタブのツリーにShell_part-1、Shell_part-2、Shell_part-3 が作成され、グラフィック・パネルには部品が表示されます。

Shell_part-1 が屋根で２枚の面から構成されています。Shell_part-2 と Shell_part-3 は壁で１枚の面から構成されています。

Shell_part-1、Shell_part-2、Shell_part-3 を選択状態にして、右クリックメニューの Create Compound Part を左クリックしてください。

![Fig16_10_02.PNG](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/888272/a048d342-78bf-5623-2f9d-cadfe408e1a8.png)

GeometryタブのツリーにCompound-1 が作成されます。

###メッシュ作成
Compound-1だけにメッシュ生成をしてください。メッシュ生成用の設定はデフォルトのままで良いです。

###材料物性値の設定

Material Library から材料 Steel：<b>S185</b> を呼び出してください。

セクションの設定をしてください。

　FE ModelタブのツリーのSectionsの右クリックメニューのCreate。
　メニューバーの場合はProperty→Section→Create。

　<b>Type</b> = Shell section
　<b>Properties</b>
　Name = Shell_selection-1 (デフォルトのまま)
　Material =  S185
　Thickness = <font color="Red">1</font>  [mm]
　Offset 0
　Region type = Selection
　<font color="Red">屋根の２つの面を選択してください</font> 

　「OK-New」を左クリック

　<b>Type</b> = Shell section
　<b>Properties</b>
　Name = Shell_selection-2 (デフォルトのまま)
　Material =  S185
　Thickness = <font color="Red">2</font>  [mm]
　Offset 0
　Region type = Selection
　<font color="Red">壁の２つの面を選択してください</font> 

　「OK」を左クリック

ツリーのSectionsの下に Shell_selection-1 と Shell_selection-2 が表示されます。

###解析設定
解析ステップを新規作成してください。

　FE ModelタブのツリーのStepsの右クリックメニューのCreate。
　メニューバーの場合はStep→Step→Createです。

　<b>Type</b> = Static step
　<b>Data</b>
　Name = Step-1 (デフォルトのまま)
　Solver = Default (デフォルトのまま)
　Nlgeom = <font color="Red">Off</font>　　(注2)
　Incrementation = Default (デフォルトのまま)

<ul>
<li>(注2)</li>板厚に比べて変形量が大きいときは On にする必要があります。今回は、Off で良いです。
</ul>

固定するエッジ(2本)を設定してください。
　FE ModelタブのツリーのStep-1/BCs
　<b>Type</b> = Fixed
　<b>Properties</b>
　Name = Fixed-1 (デフォルトのまま)
　Region type = Selection (デフォルトのまま)

圧力の面を設定してください。
　FE ModelタブのツリーのStep-1/Loads
　<b>Type</b> = Pressure
　Name = Pressure-1 (デフォルトのまま)
　Region type = Selection (デフォルトのまま)
　Magunitude = 1e-6 [MPa]
　<font color="Red">屋根の２つの面を選択してください</font> 

≪ファイル保存≫
いったんファイルを保存してください。名前は任意ですが、ここでは16_10フォルダーの16_10.pmxファイルとしました。

###ソルバー計算
ソルバー計算を実施してください。

表示例を示します。
![Fig16_10_03.PNG](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/888272/6a79e7ea-ad4a-c447-ca90-d497ebd87830.png)

###補足
 複数のシェル部品から Compound(複合部品)を作成するとシェルとして一体になり、もともとのの屋根と壁が接している線分のところで剛体結合されます。今回の屋根の部分は、16.0でのモデルも同じですが、２枚の面のCompoundです。だから「16.0 シェルの解析-基本」の時も一体となっていました。

勘違いしやすい操作を以下に記載します。<br>

【１】
<dd>Gemometryタブで、対象の全てのShell_part を Create Mesh し、FE Modelタブで Merge Mesh する。</dd>
<dd>この操作をした場合はシェル要素が連結されません。また、このあとで「4.0 部品の連結」のようなTieという操作も効きません。TieとContct/Type=Tieは、表面どうしのみで有効です。<br><b><font color="blue">このことは「16.2 シェルとSolidの組み合わせ」で重要です。</font> </b></dd>

![Fig16_10_04.PNG](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/888272/ca4f8804-1341-4568-e1df-5d485e783e77.png)

【２】
接触要素で行った場合(Shell to Edge)
壁が屋根を突き抜けています。また撓んでいる屋根の曲率の方向が逆です。
![Fig16_10_05.PNG](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/888272/de45dae9-0ed4-f7e5-2374-5ea34a6bb5d6.png)

≪接触のマスター面とスレーブ面≫
原則は、マスター側が面(Shell)でスレーブ側が線(Edge)です。
上の図で、左側の壁は原則通りですが右の壁は逆になっています。原則通りの設定となっている左側の壁でも突き抜けています。また線形要素です。デフォルトの２次要素のままで計算すると、変形後の屋根の形状は激しい凹凸になります。
左右の壁とも原則通りにすると、壁と屋根の突き抜けは発生しないものの全体が剛体移動を起こしました。

つまり、接触要素では適切な計算ができません。