Salome-Mecaを使用した メッシュ生成（非構造格子）

Transcription

1 Salome-Mecaを使用した熱伝導解析入門 & 解析手法の違いによる熱伝導解析比較 秋山善克 1

2 Salome-Meca とは EDF( フランス電力公社 ) が提供している Linux ベースのオープンソース Code_Aster : 解析ソルバー Salome-Meca : プリポストを中心とした統合プラットフォーム :SALOME Platform に Code_Aster をモジュールとして組み込んだもの Code_Aster は 構造力学 熱力学を中心に非常に高度で多彩な機能と 400 を超える要素 (1 次元 2 次元 3 次元ほか ) を有しています また 2000 以上のテストケースと ページ以上のドキュメント ( 使用方法 テクニック 理論的背景 ) 公式フォーラムなどがあり 他のオープンソース CAE ソフトと較べてサポート体制が充実しているのが特長です より インストール方法 使い方等上記ページを参照してください 3

3 解析の流れ 商用ソフト フリーソフト 設計者 CAD 図面から 3 次元 CAD で形状を作成 CATIA,Pro/E SolidWorks FreeCAD Blender,SALOME メッシュ CAD からメッシュを作成 ICEM,ANSA HyperMesh Femap SALOME Netgen Gmsh 解析技術者 条件設定 解析 解析する条件を設定 解析実行 ABAQUS-CAE MENTAT ABAQUS,MARC FLUENT,STAR-CCM Eficas HelyxOS OpenFOAM Code-Aster Elmer 結果処理 解析結果を評価する HyperView CFD-Post SALOME ParaView 4

4 Salome-Meca の構成 CAD 図面から 3 次元 CAD で形状を作成 Salome メッシュ CAD からメッシュを作成 Salome Salome-Meca 条件設定 解析する条件を設定 Eficas 解析 解析実行 Code-Aster 結果処理 解析結果を評価する Salome

5 本日の演習内容 演習 1: ウィザードによる熱伝導解析手順の確認 演習 2: 有限要素法による1 次元解析 演習 3: 差分法による数値解の比較 ( 未実施のため削除 ) 演習 4: 有限体積法による数値解の比較 6

6 Salome-Meca の起動 デスクトップ上のアイコンをクリック Geometry を選択 7

7 Geometry 起動画面 メニューバー ツールバー オブジェクトブラウザ グラフィックウィンドウ パイソンコンソール 8

8 演習 1 Primitives によるモデル作成 9

9 演習 1 Primitives によるモデル作成 円柱の作成 New Entity>Primitives>Cylinder 座標原点に作成 名前は任意 連続して作成する場合は Apply オブジェクトブラウザに追加される 表示 / 非表示切り替え 10

10 演習 1 Primitives によるモデル作成 点の作成 New Entity>Basic>Point 11

11 演習 1 Primitives によるモデル作成 球の作成 New Entity>Primitives>Sphere 点と半径を指定して作成 矢印を選択するとグラフィックウインドウまたはオブジェクトブラウザから選択可能 12

12 演習 1 Primitives によるモデル作成 13

13 演習 1 Primitives によるモデル作成 球の作成 Operations>Boolean>Cut 14

14 演習 1 Primitives によるモデル作成 15

15 演習 1 グループの作成 グループの作成 New Entity>Group>Create フェースをグループ化 名前は任意 指定のフェースを選択 16

16 演習 1 グループの作成 作成中 選択するとハイライトされる 17

17 演習 1 グループの作成 グループの作成 low up side hole 18

18 Mesh 起動画面 メニューバー ツールバー オブジェクトブラウザ グラフィックウィンドウ パイソンコンソール 19

19 Mesh>Create Mesh 演習 1 メッシュの作成 20

20 演習 1 メッシュサイズの設定 21

21 演習 1 メッシュの作成 メッシュの作成 Mesh>Compute Mesh_1 を選択 22

22 演習 1 メッシュの作成 23

23 グループの作成 Mesh>Create Group 演習 1 メッシュのグループ化 色は任意 24

24 演習 1 メッシュのグループ化 25

25 演習 1 熱伝導解析設定条件 up 温度 :100 熱伝導率 :0.54W/mK low 温度 :0 26

26 演習 1 Aster モジュールの起動 ウィザード Aster>Wizards>Linear thermic 27

27 演習 1 wizard の設定 熱伝導率の指定 温度の指定 熱流束の指定 発熱量の指定 28

28 演習 1 wizard の設定 右クリック 29

29 演習 1 解析の実行 解析実行中 30

30 メニューバー ParaViS 起動画面 ツールバー グラフィックウィンドウ オブジェクトブラウザ プロパティブラウザ パイソンコンソール 31

31 演習 1 解析結果の表示 32

32 演習 1 解析結果の表示 ( 温度 ) 33

33 演習 1 解析結果の表示 ( レンジの変更 ) 34

34 演習 解析解と数値解の比較 解析解を算出する 2-2メッシュサイズの違いによる解析結果への影響 メッシュサイズ 0.01mと0.005mで比較 2-3メッシュタイプの違いによる解析結果への影響 テトラメッシュとヘキサメッシュ ( 各辺 15 分割 ) で比較 2-4 境界条件の違いによる解析結果への影響 端部拘束と法線方向拘束で比較 35

35 演習 2 有限要素法による熱伝導解析 36 右側温度左側温度全長奥行幅 : : : : : 2 1 T T L W B B W L T2 s J W kgk J C m kg mk W s m K p 比熱密度熱伝導率熱拡散率温度 1: 1: 1: 1: : 3 2 K T K T m L m W m B 0 : :100 :10 :1 :1 2 1 T1

36 演習 2 グループの作成 グループの作成 T1 T2 wall 37

37 演習 2 メッシュグループの作成 グループの作成 T1 T2 wall 38

38 温度 演習 2 解析結果 選択状態にする Filters>Data Analysis>Plot Over Line 座標 ラインの始点と終点の座標を入力 39

39 演習 4 有限体積法による熱伝導解析 Tutorials の laplacianfoam の flange をコピー /home/user/openfoam/user-2.2.1/run/tutorials/basic/laplacianfoam/flange フォルダー名を flange に変更 Allclean Allrun Flange.ans Constant>polyMesh 内のファイルは削除する 47

40 メッシュ及びグループを表示状態メッシュを選択状態にする 演習 4 メッシュの出力 File>Export>UNV file ex4 フォルダー直下に保存 48

41 演習 4 メッシュの変換 条件設定 $ ideasunvtofoam Mesh_1.unv メッシュ変換 その他ファイルの編集 /system/controldict の編集 /0/T の編集 開始時間 終了時間 計算時間間隔 結果出力間隔 /constant/transportproperties の編集 熱伝導率 $ parafoam 49

42 温度 演習 4 計算実行 結果処理 $ laplacianfoam 解析実行 $ parafoam 座標 50

43 解析結果の比較 座標 FEM FDM FVM 離散化手法 演習 2 有限要素法 FEM 演習 3 差分法 FDM 演習 4 有限体積法 FVM 51