Axisymmetric Elastic FEA of a Pressure Vessel with Salome-Meca 2024 Windows

Salome-Mecaによる 軸対称体の弾性解析 （圧⼒容器の耐圧設計） PE/計算固体⼒学研究室 1

学習⽬標 • 軸対称体の弾性解析を⾏い，簡単な耐圧設計式への 適⽤により整合性を確認する • 2Dスケッチを⽤いて，軸対称体のモデル作成を体 験する • AsterStudy を⽤いて，コマンドファイルを編集し， 簡単な弾性解析におけるコマンド構成を把握する • 解析結果のテキスト印刷出⼒を⽤いて，ワークシー トによる後処理を体験する 2

解析の流れ 問題の設定 Geometry Module • 解析対象の形状，⼨法，材料特性の準備 • 解析結果処理に必要な情報の確認 • 形状の作成（2D Sketchを使⽤） • 境界条件を与える形状グループに名前をつける Mesh Module • メッシュの作成 • 要素に形状グループ名を継承 Aster Module • AsterStudyを⽤いたコマンドファイルの作成 • Studyケースの設定（コマンド＋メッシュ） • 解析の実⾏ ParaVis Module • 変形図によるおおまかな妥当性の確認 • 応⼒分布の確認 結果の評価 • ワークシートによる断⾯平均応⼒計算 3

解析対象（形状） 上下対称な容器 断⾯図 モデル化対象 （1/8容器の断⾯） 4

y (z) y (z) 素材︓炭素鋼（steel） を想定し， 内径 Di = 1000 O (0, 0, 0) 胴部半⾼さ H = 1000 ⾁厚 w = 50 解析対象（境界条件，材料特性） E = 200 GPa  = 0.3 Sm = 100 MPa 対称⾯ uy = 0 x (r) ⼨法の諸元（単位︓mm） O (0, 0, 0) とすれば、 許容圧⼒はどうなるか x (r) 境界条件 5

解析⽅法 • 軸対称2次要素を使⽤（modelisation=AXIS） • Salome-Meca2024 for windows（Code_aster ver. 16.7.9）を使⽤し，AsterStudy の Assistant 機能と編集機能を⽤いてコマンドファイルを作成 • 対称⾯（GROUP_MA=Symm）上の節点に対して、 応⼒成分をテキスト出⼒する • 応⼒成分のテキスト出⼒を Microsoft Excel TM に 読込み，解析結果のグラフ化，分析（断⾯平均算出 など）を⾏う 6

形状データの作成 • Geometryモジュールにある2D Sketchを使⽤して， 容器形状の外形線を描く • 荷重（内圧）をかける⾯，後処理⽤にデータを出⼒ する⾯（対称⾯）に名前をつけておく • 形状作成までの作業を Script に保存する 7

SalomeMecaの起動 (1) run_salome.bat をダブルクリック 8

Geometryの起動 (1) プルダウンメニュー から Geometry を選択 9

2D Sketchの起動 (1) New Entity -> Basic -> 2D Sketch (0, 1550) (0, 1500) (0, 1000) 円弧の (500, 1000) 中⼼ (500, 0) (550, 1000) (550, 0) スケッチする図形 10

2Dスケッチの作成（直線部） (1) グローバル座標系 を確認 (0, 1550) (0, 1500) (0, 1000) 円弧の (500, 1000) 中⼼ (500, 0) (550, 1000) (2) 直線を確認 (550, 0) スケッチする図形 (4) 適⽤する (3) 直線の始点座標を⼊⼒ X=550, Y=0 ※ 反時計回りにスケッチを作成 11

2Dスケッチの作成（直線部） (0, 1550) (0, 1500) (0, 1000) 円弧の (500, 1000) 中⼼ (1) 直線の終点座標を⼊⼒ X=550, Y=1000 (500, 0) (550, 1000) (550, 0) スケッチする図形 (2) 適⽤する 12

2Dスケッチの作成（円弧部） (1) 円弧を選択 (5) Fit All で全体の 形状を確認 (0, 1550) (0, 1500) (2) 中⼼をパラ メータに追加 (4) 適⽤する (0, 1000) 円弧の (500, 1000) 中⼼ (500, 0) (550, 1000) (550, 0) スケッチする図形 (3) 中⼼と終点座標を⼊⼒ X Center = 0, Y Center = 1000 X = 0, Y = 1550 13

2Dスケッチの作成（直線部） (1) 直線にもどす (0, 1550) (3) 適⽤する 3 (0, 1500) (0, 1000) 円弧の (500, 1000) 中⼼ (500, 0) (2) 頭部内側の座標を⼊⼒ X=0, Y=1500 (550, 1000) (550, 0) スケッチする図形 14

2Dスケッチの作成（円弧部） (0, 1550) (1) 円弧を選択 (0, 1500) (0, 1000) 円弧の (500, 1000) 中⼼ (2) 中⼼をパラ メータに追加 (500, 0) (4) 適⽤する (550, 1000) (550, 0) スケッチする図形 (3) 中⼼と終点座標を⼊⼒ X Center = 0, Y Center = 1000 X = 500, Y = 1000 15

2Dスケッチの作成（直線部） (1) 直線にもどす (0, 1550) (3) 適⽤する (0, 1500) (0, 1000) 円弧の (500, 1000) 中⼼ (500, 0) (2) 中⼼内側の座標を⼊⼒ X=500, Y=0 (550, 1000) (550, 0) スケッチする図形 16

2Dスケッチの終了 (1) Sketch を 閉じる 17

2Dスケッチの作成（確認） (1) Geometry の中に Sketch_1 ができていることを確認（⽬の アイコンを点滅させる） (2) Sketch_1 の点滅 を確認 18

フェースの作成 (2) New Entity -> Build -> Face (1) Sketch_1を選択 19

フェースの作成 (2) Face の名前（Face_1） とオブジェクト （Sketch_1）を確認 (1) Face を確認 （⾊が変わる） (3) 適⽤して閉じる 20

形状グループの作成 (1) Face_1を右ク リック (2) グループを作成 21

形状グループの作成（圧⼒⾯） (1) 「線」を選択 (2) 名前は「Press」 (3) 内⾯の円弧 をピック (7) 2つの線の 表⽰を確認 (4) 追加する 5 (5) 内⾯の直線 をピック (6) 追加する (8) 適⽤する 22

形状グループの作成（対称⾯） (1) 名前は 「Symm」 (3) 追加する 4 (4) 線番号を 確認 (5) 適⽤して閉じる (2)下端の線（対 称⾯）をピック 5 23

形状グループの作成（確認） (1) Face_1の中のPress とSymmを確認 24

Script の保存（Geometryまで） (1) File -> Dump Study... 25

作業の保存（Geometryまで） (1) 作業フォルダ を開く (2) ファイル名は 「Geometry1.py」 (3) 保存する 26

Geometry1.py の確認 geomObj_1 = geompy.MakeMarker(0, 0, 0, 1, 0, 0, 0, 1, 0) sk = geompy.Sketcher2D() sk.addPoint(550.000000, 0.000000) ### sk.addSegmentAbsolute(550.000000, 1000.000000) ### This file is generated automatically by SALOME v9.12.0 with dump sk.addArcCenterAbsolute(0.000000, 1000.000000, 0.000000, 1550.000000, 0, python functionality 0) ### sk.addSegmentAbsolute(0.000000, 1500.000000) sk.addArcCenterAbsolute(0.000000, 1000.000000, 500.000000, 1000.000000, import sys 0, 0) import salome sk.addSegmentAbsolute(500.000000, 0.000000) sk.close() salome.salome_init() Sketch_1 = sk.wire(geomObj_1) import salome_notebook Face_1 = geompy.MakeFaceWires([Sketch_1], 1) notebook = salome_notebook.NoteBook() sys.path.insert(0, r'C:/OpenCAE/Salome2024Win/v2024/MyCases/PresVes') Press = geompy.CreateGroup(Face_1, geompy.ShapeType["EDGE"]) geompy.UnionIDs(Press, [10, 12]) Symm = geompy.CreateGroup(Face_1, geompy.ShapeType["EDGE"]) ### geompy.UnionIDs(Symm, [14]) ### GEOM component geompy.addToStudy( O, 'O' ) ### geompy.addToStudy( OX, 'OX' ) geompy.addToStudy( OY, 'OY' ) import GEOM geompy.addToStudy( OZ, 'OZ' ) from salome.geom import geomBuilder geompy.addToStudy( Sketch_1, 'Sketch_1' ) import math geompy.addToStudy( Face_1, 'Face_1' ) import SALOMEDS geompy.addToStudyInFather( Face_1, Press, 'Press' ) geompy.addToStudyInFather( Face_1, Symm, 'Symm' ) geompy = geomBuilder.New() #!/usr/bin/env python O = geompy.MakeVertex(0, 0, 0) OX = geompy.MakeVectorDXDYDZ(1, 0, 0) OY = geompy.MakeVectorDXDYDZ(0, 1, 0) OZ = geompy.MakeVectorDXDYDZ(0, 0, 1) if salome.sg.hasDesktop(): salome.sg.updateObjBrowser() ※ テキストエディタ で開くとこのように なっている 27

メッシュデータの作成 • Meshモジュールを⽤いてメッシュ（有限要素によ るモデル）を作成する • 板厚⽅向にサブメッシュを定義し，分割数を制御す る • メッシュデータをファイル出⼒（エクスポート）す る • メッシュ作成までの作業を Script に保存する 28

Meshの起動 (1) Meshモジュー ルを起動 29

Meshの作成 (3) Fit All で全体表⽰ (2) -OZ の⽅向 に表⽰ (1) Face_1を選択して， 表⽰させる（「⽬」 のマークをオン） 30

Mesh の作成 (2) Mesh -> Create Mesh 1 (1) Face_1 を選択 した状態で 31

Mesh の作成 (1) 分割⽅法を指定 した四辺形分割 32

分割数の指定（全体） (1) セグメントの 分割数は 50 (2) OKする (3) 適⽤して閉じる 33

サブメッシュの作成 (1) Mesh_1 を 右クリック ※ Mesh_1 を選択した状態で Mesh -> Create Sub-mesh で も作成できる 2 (2) サブメッシュ の作成 34

サブメッシュの作成 (1) Mesh_1 の中に SubMesh_1を作成し， Geometry は Symm とする 2 (2) 1Dタブ中のアルゴリ ズムから「線の離散化」 を選択 35

サブメッシュの作成 (1) 設定（⻭⾞マー ク）をクリック 2 (2) セグメントの 分割数を選ぶ 36

サブメッシュの作成 ① 分割数を選ぶ 37

サブメッシュの作成 (1) 5分割にする (2) OKする 38

サブメッシュの作成 (1) 追加の仮定の設定（⻭⾞ マーク）から「反対側の エッジまで貫通」を選択 (2) 適⽤して閉じる 39

メッシュの作成 (1) Mesh_1 を右ク リック 2 (2) メッシュ⽣成 （compute） 40

メッシュの作成（結果） (2) メッシュの 情報を確認 (1) メッシュの確認 (3) 閉じる 41

メッシュの確認 (1) ボックス拡⼤ （Fit Area）する 2 (2) マウスで拡⼤表⽰したい 部分を囲む（囲む四⾓形の 対⾓線を左ボタンでドラッ グする）と拡⼤される 42

⼆次要素化 (2) ⼆次要素への変換 (1) Mesh_1 を右 クリック 43

⼆次要素化 (1) 「中間節点を形状 （例えば円弧）の上に 置く」をチェック (2) ⼆次要素への変換 2 (3) 適⽤して閉じる 44

メッシュ情報の確認 2 (2) メッシュ情報 (1) Mesh_1 を右 クリック 45

メッシュ情報の確認 (1) 線形要素（Linear）がゼロ 個，⼆次要素（Quadratic）が 500個になっている (2) OK する 46

メッシュファイルの保存 (2) Export -> MED file （バイナリファイ ル）を選択 (1) Mesh_1 を 右クリック 47

Script の保存 1 (1) File -> Dump Study 48

Script の保存 (1) 作業フォルダ を開く (2) ファイル名は 「Mesh_1.py」 (3) 保存する 49

Mesh1.py の確認 smesh.SetName(Quadrangle_2D.GetAlgorithm(), 'Quadrangle_2D') smesh.SetName(Number_of_Segments_2, 'Number of Segments_2') smesh.SetName(Propagation_of_1D_Hyp, 'Propagation of 1D Hyp. on Opposite Edges_1') smesh.SetName(Number_of_Segments_1, 'Number of Segments_1') import SMESH, SALOMEDS smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1') from salome.smesh import smeshBuilder smesh.SetName(Press_1, 'Press') smesh.SetName(Symm_1, 'Symm') smesh = smeshBuilder.New() #smesh.SetEnablePublish( False ) # Set to False to avoid publish in study if smesh.SetName(Sub_mesh_1, 'Sub-mesh_1') not needed or in some particular situations: # multiples meshes built in parallel, complex and if salome.sg.hasDesktop(): numerous mesh edition (performance) salome.sg.updateObjBrowser() Mesh_1 = smesh.Mesh(Face_1,'Mesh_1') Regular_1D = Mesh_1.Segment() Number_of_Segments_1 = Regular_1D.NumberOfSegments(50) Quadrangle_2D = Mesh_1.Quadrangle(algo=smeshBuilder.QUADRANGLE) Press_1 = Mesh_1.GroupOnGeom(Press,'Press',SMESH.EDGE) Symm_1 = Mesh_1.GroupOnGeom(Symm,'Symm',SMESH.EDGE) Regular_1D_1 = Mesh_1.Segment(geom=Symm) Number_of_Segments_2 = Regular_1D_1.NumberOfSegments(5) Propagation_of_1D_Hyp = Regular_1D_1.Propagation() isDone = Mesh_1.Compute() [ Press_1, Symm_1 ] = Mesh_1.GetGroups() Mesh_1.ConvertToQuadratic(0) [ Press_1, Symm_1 ] = Mesh_1.GetGroups() Sub_mesh_1 = Regular_1D_1.GetSubMesh() ### ### SMESH component ### ※ テキストエディタで開くとこの ようになっている（Geometry1.py と共通部分は除く） ## Set names of Mesh objects smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D') 50

解析条件データの作成 • AsterStudy の 等⽅線形弾性⽤ Assistant 機能 を⽤いて新しい Stage を作成する • AsterStudy を⽤いてファイルのメッシュファ イルとバイナリ出⼒ファイルの割当を⾏う • AsterStudy を⽤いてコマンドファイルを修正 する • ワークシートで後処理を⾏うデータをテキス トファイルに出⼒する • 解析を実⾏する 51

AsterStudyの起動 1 (1) AsterStudyを起動 52

Code_Aster のバージョンの確認 (1) 安定版を確認 (2) OK する 53

Assistant の起動 (1) Operation -> Add Stage with Assistant -> Isotropic linear elasticity（等⽅線形弾性） 54

Assistant の起動 (1) 注意書きを確認 (2) 次に進む 55

メッシュの選択 (1) メッシュの確認 (2) 次に進む 56

モデルの選択 (1) モデルには AXIS （軸対称）を選ぶ (2) 次に進む 2 57

材料特性の⼊⼒ (1) ヤング率は 200000 MPa， ポアソン⽐は 0.3 2 (2) 次に進む 58

拘束条件の指定 (1) 緑⾊のプラスをクリック すると Group 選択できる (2) 編集する 59

拘束条件の指定 (1) 拘束するグループ は Symm (2) OK する 60

拘束条件の指定 (1) Y ⽅向（軸⽅向） 変位（DY）をゼロに 拘束 (2) 次に進む 61

荷重（内圧）の⼊⼒ ③ Press を選ぶ ① 緑⾊のプラスをクリック して Group 選択を表⽰ ④ OK する ② 編集する ⑤ 圧⼒は 10 MPa ⑥ 次に進む 62

荷重（内圧）の⼊⼒ (1) 緑⾊のプラスをクリック して Group 選択を表⽰ (2) 編集する 63

荷重（内圧）の⼊⼒ (1) 要素グループ 「Press」を選ぶ (2) OK する 64

荷重（内圧）の⼊⼒ (1) 内圧は 10 MPa 2 (2) 次に進む 65

出⼒ファイルの指定 (1) 出⼒ファイル を探す 66

出⼒ファイルの指定 (1) 作業フォルダを 開く (2) ファイル名は 「PresVes.rmed」 (3) 保存する 67

Assistant の終了 (1) 終了 68

メッシュの読込み（LIRE_MAILLAGE） (1) 作成されているコマンド ファイルからLIRE_MAILLAGE （メッシュの読込み）をダブ ルクリック (3) 「...」をクリックして メッシュファイルを探す (2) フォーマットは バイナリ（MED） 69

メッシュファイルの読込み 1 (1) 作業フォルダ を開く (1) 保存してある Mesh_1.med を選ぶ (3) 開く 70

変更の適⽤とコマンドファイルの確認 (2) コマンドファイ ルが更新される (1) 適⽤する (3) OK する 3 71

モデルの割当（AFFE_MODELE） (1) AFFE_MODELE（モデル の割当）をダブルクリック (2) アイテムを開く （▼をクリック） (3) [0] 番⽬のモデル を開く 72

モデルの割当（AFFE_MODELE） (1) モデルの対象は TOUT （=all）を OUI（=yes） (2) PHENOMENE（現象） はMECANIQUE（機械）， MODELISATION（モデル 化）は AXIS（軸対称） 2 (3) OK する 73

モデルの割当（AFFE_MODELE） (1) 適⽤する (3) OK する (2) コマンドファイ ルの更新を確認 74

解析結果の後処理（CALC_CHAMP） (1) Commands -> Post Processing -> CAL_CHAMP （場の計算）を開く 75

解析結果の後処理（CALC_CHAMP） (1) 後処理結果は解析と 同じコンセプト名を継承 (2) CONTRAINTE（応⼒） のアイテムを増やし， SIGM_NOEU（節点上の応 ⼒成分）を選ぶ 76

解析結果の後処理（CALC_CHAMP） (1) 同様にして，CRITERES（基 準）を SIEQ_NOEU（節点上の相 当応⼒）を追加 (3) コマンドファイル が更新される (2) 適⽤する (4) OK する 4 77

結果の印刷出⼒（IMPR_RESU） (3) 「...」をクリックして 結果ファイルを探す (1) IMPR_RESU（結果 の印刷出⼒）をダブル クリック (2) フォーマットは バイナリ（MED） 2 78

結果の印刷出⼒（IMPR_RESU） (1) 作業フォルダ を開く (2) ファイル名を確 認（拡張⼦は rmed） (3) 保存 79

結果の印刷出⼒（IMPR_RESU） (1) 結果を編集 80

結果の印刷出⼒（IMPR_RESU） (1) 出⼒するRESULTAT（結果） は result(CALC_CHAMP) (2) 結果は TOUT_CHAM（=all fields）を OUI（=yes）にする (3) OK する 81

結果の印刷出⼒（IMPR_RESU） (2) コマンドファイル が更新される (1) 適⽤する (3) OK する 82

結果の印刷出⼒（IMPR_RESU）テキスト出⼒ (1) Commands -> Output -> IMPR_RESU（結果の印刷出 ⼒）を開いて，2つ⽬の出⼒ 設定を作る 83

結果の印刷出⼒（IMPR_RESU） (2) 「...」をクリックして 結果ファイルを探す ① フォーマットはテ キスト（RESULTAT） 84

結果の印刷出⼒（IMPR_RESU） 1 (1) 作業フォルダ を開く ファイル名は PresVes.txt (3) 保存 85

結果の印刷出⼒（IMPR_RESU） (1) RESUにアイテム を追加 (2) 編集する 86

結果の印刷出⼒（IMPR_RESU） (1) 出⼒するRESULTAT（結果） は result(CALC_CHAMP) (2) NOM_CHAM（場の名 前）を SIGM_NOEU と SIEQ_NOEU にする (3) 下にスク ロールする 87

結果の印刷出⼒（IMPR_RESU） (2) 出⼒する形状グループ を編集する (1) IMPR_COOR（座標の印 刷）を OUI（=yes）に 88

結果の印刷出⼒（IMPR_RESU） (1) Symmを選ぶ (2) OK する 89

結果の印刷出⼒（IMPR_RESU） (1) OK する 1 90

結果の印刷出⼒（IMPR_RESU） (2) コマンドファイル が更新される (1) 適⽤する (3) OK する 91

コマンドファイルの保存 (1) ステージ名を 右クリック (2) コマンドファイル をエキスポート 92

コマンドファイルの保存 (1) 作業フォルダ を開く (2) ファイル名は Stage_1.comm (3) 保存 93

コマンドファイルの確認 DEBUT(LANG='FR') mesh = LIRE_MAILLAGE(FORMAT='MED', UNITE=2) model = AFFE_MODELE(AFFE=_F(MODELISATION='AXIS', PHENOMENE='MECANIQUE', TOUT='OUI'), MAILLAGE=mesh) mater = DEFI_MATERIAU(ELAS=_F(E=200000.0, NU=0.3)) materfl = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ), TOUT='OUI'), MODELE=model) mecabc = AFFE_CHAR_MECA(DDL_IMPO=_F(DY=0.0, GROUP_MA=('Symm', )), MODELE=model) mecach = AFFE_CHAR_MECA(MODELE=model, PRES_REP=_F(GROUP_MA=('Press', ), PRES=10.0)) result = MECA_STATIQUE(CHAM_MATER=materfl, EXCIT=(_F(CHARGE=mecabc), _F(CHARGE=mecach)), MODELE=model) result = CALC_CHAMP(reuse=result, CONTRAINTE=('SIGM_NOEU', ), CRITERES=('SIEQ_NOEU', ), RESULTAT=result) IMPR_RESU(FORMAT='MED', RESU=_F(RESULTAT=result, TOUT_CHAM='OUI'), UNITE=80) IMPR_RESU(FORMAT='RESULTAT', RESU=_F(GROUP_MA=('Symm', ), IMPR_COOR='OUI', NOM_CHAM=('SIGM_NOEU', 'SIEQ_NOEU'), RESULTAT=result), UNITE=8) FIN() 94

解析ケース（Study）の保存 (1) File -> Save As... を開く 95

作業ケース（Study）の保存 (1) ファイル名は Study1 (2) 保存 96

解析ケース（Study）の実⾏ (4) 解析ケースの追加 （緑⾊のプラス） (1) History タブを開く (2) 解析に使うメモリ， 上限時間を確認 (5) Run (3) ⾃動更新5秒毎に 97

解析ケース（Study）の終了 (1) 緑マークでエラー なく解析が終了したこ とを確認 98

解析結果の確認と利⽤ • ParaVisを⽤いて，解析結果を可視化し，定性的 な妥当性を確認する • 解析結果のテキストファイル（*.txt）を開き， 対称⾯の応⼒分布を確認する • 対称⾯の応⼒を Excel にコピーし，設計評価に 必要な主応⼒，相当応⼒（⼀次⼀般膜応⼒）を 計算する 99

ParaVisの起動 (1) Case View タブ を開く (2) rmed ファイル を右クリック (3) ParaVis で開く 100

結果の読込み (3) 適⽤する (1) 読込むデータに チェックを⼊れる (2) ベクトルを ⽣成する 101

変形図の表⽰ (1) Filters -> Common -> Warp By Vector を開く 102

変形図の表⽰ (3) 変形を確認 (2) 適⽤する (1) Scale Factor は⻘いバーを真 ん中あたりに置く 103

変形図上への相当応⼒のカラーコンターの表⽰ (1) result_SIEQ_NOEU（節点 上の応⼒），VMIS（フォン ミーゼス相当応⼒）を選ぶ 2 (2) 適⽤する 104

結果テキストファイルの処理 -------------------------------------------------------------------------------CONCEPT 00000007 DE TYPE EVOL_ELAS ENTITES TOPOLOGIQUES SELECTIONNEES GROUP_MA : Symm ======> (1) 出⼒されたテキストファイル をエディターで開き，応⼒成分 （SIGM_NOEU）の部分をコピー する ------> CHAMP AUX NOEUDS DE NOM SYMBOLIQUE SIGM_NOEU NUMERO D'ORDRE: 1 INST: 0.00000000000000E+00 NOEUD X Y SIXX SIYY SIZZ SIXY N1 5.50000000000000E+02 2.27373675443232E-15 8.64937663409215E-03 4.76489979183474E+01 9.52397082043099E+01 -2.95585387162626E-06 N6 5.00000000000000E+02 2.27373675443232E-15 -9.98792139387641E+00 4.75893242156517E+01 1.05219224126309E+02 5.65005519440946E-06 N207 5.10000000000000E+02 2.27373675443232E-15 -7.75129026639781E+00 4.76013844800141E+01 1.02985965603772E+02 9.87073265220388E-07 N208 5.20000000000000E+02 2.27373675443232E-15 -5.64262713599956E+00 4.76133821906997E+01 1.00880737677645E+02 1.07138529169108E-06 N209 5.30000000000000E+02 2.27373675443232E-15 -3.65209061925905E+00 4.76253109650605E+01 9.88936369303453E+01 1.03598268374566E-06 N210 5.40000000000000E+02 2.27373675443232E-15 -1.77102201812760E+00 4.76371831992362E+01 9.70159974260350E+01 4.66514293023739E-07 N812 5.05000000000000E+02 2.27373675443232E-15 -8.86946745187538E+00 4.75953467870597E+01 1.04102430602500E+02 3.17983418943223E-06 N813 5.15000000000000E+02 2.27373675443232E-15 -6.69697170466574E+00 4.76073840986936E+01 1.01933366799663E+02 1.11669882345327E-06 N814 5.25000000000000E+02 2.27373675443232E-15 -4.64737028670956E+00 4.76193472059167E+01 9.98872008466797E+01 1.16184310434922E-06 N815 5.35000000000000E+02 2.27373675443232E-15 -2.71156649628747E+00 4.76312475105768E+01 9.79548292349277E+01 8.34687314936836E-07 N816 5.45000000000000E+02 2.27373675443232E-15 -8.81290108867037E-01 4.76430962997634E+01 9.61279763193355E+01 -1.38500723709909E-06 105

ワークシートへのデータのコピー (1) ワークシート（ここで は Excel）を起動し，デー タを貼り付ける ※ 1⾏全体が⽂字列としてA 列にコピーされる 106

ワークシートへのデータのコピー (1) データタブ を開く (3) 固定⻑フィー ルドを選ぶ (2) 区切り位置の 設定 (4) 次へ 107

データのセル区切り設定 (1) 区切り位置を調整 （必要に応じて区切 り線を修正） (3) 完了 (2) 次へ 2 108

ワークシートでの後処理 (1) セルに適切に配分 されていることを確認 109

ワークシートでのグラフ作成 (2) X と SIZZ（周⽅ 向応⼒成分）をグラ フにする (1) X（半径⽅向距 離）が⼩さい順に並 び替える (3) グラフ作成（板厚⽅向平均応⼒ は 100 MPa をわずかに超える） 110

まとめ • 圧⼒容器の弾性解析を例題として，2Dスケッチに よる形状作成，メッシュ作成を体験した • AsterStudy の Assistant を⽤いたコマンドファイ ルを体験した • AsterStudy の GUI でコマンドファイルを修正し， 解析を実⾏した • 解析結果を可視化し，妥当性を確認した • 解析結果のテキストファイルからデータを抽出し， ワークシートを⽤いて結果をグラフ化した 111