#author("2023-07-28T00:29:32+09:00","default:kouzouken","kouzouken")
#author("2023-07-28T00:29:52+09:00","default:kouzouken","kouzouken")
#contents

(注意)応力ひずみ曲線のグラフにひずみ(mm)と書いてありますが、ひずみは無次元なので単位はありません。~
表は訂正しましたが、グラフは修正できない(グラフそのものを編集するファイルを捨ててしまったため)ので、mmは無視してください。~
[[ひずみの詳しい説明:https://www.str.ce.akita-u.ac.jp/~gotou/kouzou/kouzou2_02.html#hizumido]]
[[ひずみの詳しい説明はこちらから:https://www.str.ce.akita-u.ac.jp/~gotou/kouzou/kouzou2_02.html#hizumido]]

*課題で取り組んだこと [#zb051450]
発表日:2023年4月14日(金) 再発表日:2023年4月19日(水)再々発表日:2023年4月26日(水)

片持梁と単純梁を引張と曲げで解析、応力ひずみ曲線を作成して降伏応力が235N/mm^2になることを確かめる。(解析対象はSM400)

*解析モデル [#kd9feada]
引張の場合
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/dansohipari.png

曲げの場合
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/dansomage.png

曲げの際にかけた荷重
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/magekajuu.png

*曲げについての考察(4/26発表用) (Marcで解析した結果も載せました) [#wb1028d5]

・片持ち梁の固定面の数値化

ECRO_LINEの降伏後のヤング率の数値を1000,1,0.001でそれぞれ設定し比較したところ、降伏前までの値(応力,ひずみ)の変化は同じだったので、降伏後のヤング率を0.001で解析したデータを使う。再解析で見た点②を基準として、②が降伏した際の同じpastimeの応力とひずみのデータを取る。降伏点は応力とひずみの変化が一定でなくなったところ(比例関係でなくなったところ)をとる。

8MPaで曲げた際のひずみの様子(手前が固定面)

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/zukataEPSI.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/zusuuji.png

,,ひずみ,応力(MPa)
,①,0.001042,249.267
,②,0.000936,234.025
,③,0.000984,236.751
,④,-0.000013,-1.57226
,⑤,-0.001037,-252.657
,⑥,-0.000935,-233.743
,⑦,-0.001067,-252.314
,⑧,-0.000010,-0.717033
,⑨,0.000101,26.7453

7.5mm(ひずみが最大)でsliceしたところと比較

片持ち梁曲げ(固定面から7.5mm)
,,ひずみ,応力(MPa)
,①,0.001072,221.063
,②,0.001084,223.697
,③,0.001061,218.876
,④,0.000123,25.3584
,⑤,-0.001076,-220.903
,⑥,-0.001111,-227.414
,⑦,-0.001065,-218.376
,⑧,-0.000047,-9.75992
,⑨,-0.000024,-4.57733

分かったこと

前の降伏値の求め方(ヤング率が約1000MPa下がったところで降伏と判断)から今回のやり方に変えたことで固定面でも235MPaに近い値で降伏する結果となった。固定面からz軸正方向に向かうほど、応力は徐々に小さくなっていくので、片持ち梁のモーメント図と同じ変化となったと言えるだろう。(曲げモーメント最大=部材の曲げる力最大)

また、断面の四隅は面の真ん中の点(②,⑥)より大きな値が出ていることから、四隅は外力の影響が受けやすい場所だと考えられる。つまり、最初の発表までに用いていた解析法、面(または線)に存在する全ての点を平均したデータを使って応力ひずみ曲線を描くのは外れ値も含んでしまうためあまりよろしくない手法であったと考えた。

・曲げが引張のときみたいにグラフが曲がらない理由

salomeの設定は問題なかった。また、載荷値を下げてみてもグラフの形状は引張のように曲がることはなく、緩やかな曲線であった。(ECRO_LINEの降伏後のヤング率の数値を小さくしても、解決しなかった。)

曲げの場合、引張のように一気に降伏するのではなく、固定面側からじわじわと部材が降伏する(材質がゆっくり柔らかくなっていく)ため、緩やかな曲線を描くグラフができるのではないかと思われる。

もしかすると、asterstudyで曲げ降伏の解析設定を行う際、引張解析の設定に更に入力する必要のあるコマンドがあるかもしれない。今回、探した範囲では分からなかったので、これからの研究活動で何かしら分かったことがあれば追記しようと思う。

5/2(火)追記:Marc Mentatというツールを使用して片持ち梁の曲げ解析を行い、応力-ひずみ曲線を作成した。グラフを見れば分かると思うが、設定した降伏応力(235MPa)を超えることはなく直線になっている。今回、塑性後のひずみと応力を入力して弾塑性解析を行なった。この方法はSalome-MecaでいうTRACTIONに近いものだと思われるので、ECRO_LINEではなくTRACTIONを使うべきなのだろうか...?

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/Mdanso.png

*発表後の再解析についてのまとめ(4/19発表用) [#kf5c3e65]

**解析手法 [#lfea90fd]

幅10mm✕厚さ10mmのモデルで解析。モデルの設定等は上記の解析モデルと同じだが、単純梁については片方の支承をz方向に自由に動けるように再設定した。(x,y方向は固定)

他のモデルは時間の関係上省略させていただきます。(解析1つに相当の時間がかかるため。)

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/zuTH.png

引張:梁の長さの中間地点でsliceして上側と下側を1点ずつ見る。(片持ち梁:50mm 単純梁:60mm)

これは梁の変形を考慮して、最も変形しにくいと思われる場所が梁の中間だと考えたため。

曲げ:片持ち梁は梁の固定面から少しずれた位置(7mm)でsliceして、引張側の1点で見る。(→ひずみが最大のところで見た。)単純梁は載荷した場所をsliceして、引張側の1点で見る。

・赤丸の点が解析の際に見たところ 

引張の場合、上と下側の赤丸。曲げの場合、片持ち梁は上側の赤丸,単純梁は下側の赤丸。

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/haru1.png

**解析結果 [#iab2c9bf]

片持ち梁引張

降伏時の応力とひずみ 上と下の点、どちらも同じデータが得られました。

,ひずみ,応力(MPa)
,0.00114,233.64

以前の降伏時の応力とひずみデータ

,ひずみ,応力(MPa)
,0.00114,233.64

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/saiKH.png

片持ち梁曲げ

降伏時の応力とひずみ

,ひずみ,応力(MPa)
,0.00114,233.497

以前の降伏時の応力とひずみデータ

,ひずみ,応力(MPa)
,0.00110,229.697

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/saiKM.png

単純梁引張

降伏時の応力とひずみ 上と下の点、どちらも同じデータが得られました。

,ひずみ,応力(MPa)
,0.00114,233.64

以前の降伏時の応力とひずみデータ

,ひずみ,応力(MPa)
,0.00114,234.002


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/saiTH.png

単純梁曲げ

降伏時の応力とひずみ

,ひずみ,応力(MPa)
,0.001139,234.487

以前の降伏時の応力とひずみデータ

,ひずみ,応力(MPa)
,0.00107,218.189

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/saiTM.png


**コメント [#p664fd21]
平均で見る方法からある1点のみを見る方法に変更した結果、SM400の降伏応力である235MPaに近い値で降伏したことが確認できた。引張については以前の結果とほぼ同じ値が得られたが、曲げについては値が大きく変化したので、以前の方法では外れ値の影響を受けていたことが考えられる。
しかし、1点どこを見てもいいわけではなく、影響の受けにくい場所を選んで見る必要があるだろう。

気になった点として、単純梁の引張については片持ち張りの引張と同じ結果が得られたことである。これは片方の支承をz方向に自由に移動できるようになったため、そちら側に大きく引っ張られたことで片持ち張り引張と似たような変形が起きたために同じ結果が得られたと考えられる。

*解析結果の一例(引張のみ) [#y52507b7]
・片持ち梁引張

強制変位(5mm)(DDL_IMPO)を使った解析結果[モデルは10mm✕5mm✕100mm]

,ひずみ,応力(MPa)
,0,0
,0.00025,50.1
,0.00050,100.2
,0.00075,150.3
,0.00100,200.4
,0.00124,234.9
,0.00152,235.2
,0.00179,235.5
,0.00206,235.7
,0.00233,236.0
,0.00259,236.3
,0.00286,236.5
,0.0031,236.8

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/2hipari5mm.png

荷重(FORCE_FACEを236MPaに設定)を使った解析結果(解析値に関しては最初の5点と最後の5点を載せました。)[モデルは10mm✕5mm✕100mm]

,ひずみ,応力(MPa)
,0,0
,1.18001E-05,2.36002
,2.36006E-05,4.72007
,3.54013E-05,7.08015
,4.72024E-05,9.44027
,0.00113418,226.714
,0.001146,229.077
,0.00115783,231.441
,0.00116966,233.804
,0.00284734,236.56

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/hipari236.png

・単純梁の引張[モデルは10mm✕15mm✕120mm]

降伏時の応力とひずみは233.8MPa,0.0011

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/tan1015100.png

まとめてみたもの

Kは片持ち梁、Tは単純梁、Mは曲げ、Hは引張

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/dannsoseikoufuku.png




*メモ [#wcaf6252]

**TRACTIONとECRO_LINEの注意点 [#e973d282]
TRACTIONで設定する場合、AnalysisのCOMPORTEMENT→RELATIONをVMIS_ISOT_TRACに設定すること。
ECRO_LINEの場合、AnalysisのCOMPORTEMENT→RELATIONをVMIS_ISOT_LINEに設定すること。
なお、ECRO_LINEの設定のやり方については[[Salome-Meca_降伏点検証 :https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca_%E9%99%8D%E4%BC%8F%E7%82%B9%E6%A4%9C%E8%A8%BC]]を見てください。

**部材ごとに応力を見る方法 [#p10b3527]
-パイプのブラウザーでファイルの選択(medファイル)>Apply>フィルター >Alphabetical>Extract Group>見たい部材にチェックを入れる(2つあるので一応どっちにも入れる)>Apply>Spread sheet viewのShowingでExtract Groupを選択>select points throughで見たい部分を選択する>選んだポイントがピンク色で表示される>Spread sheet viewで選択した部分のデータが緑色に表示される 
Spread sheet viewで全部選択し、Plot section overtimeを行うことにより各パスタイムごとの平均値を出力することが可能となる。

補足:1つ選択した場合、その点のpastimeごとの変化を見ることができます。どこの点を取ってPlot section overtimeを行うかはお任せします。(1点だけでもいいし、必要な点だけ取ってもいい。全部取るのもあり?)ただ、何故そのように取ったか説明できるようにすること。


*参考資料 [#sd0cd727]
-[[Salome-Meca演習_弾塑性解析(2021) :https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca%E6%BC%94%E7%BF%92_%E5%BC%BE%E5%A1%91%E6%80%A7%E8%A7%A3%E6%9E%90%282021%29]](Salome-Mecaでの設定は基本このページと同じように行った。)
-[[Salome-Meca_降伏点検証 :https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca_%E9%99%8D%E4%BC%8F%E7%82%B9%E6%A4%9C%E8%A8%BC]]
-[[Salome-Meca演習_弾塑性解析:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca%E6%BC%94%E7%BF%92_%E5%BC%BE%E5%A1%91%E6%80%A7%E8%A7%A3%E6%9E%90#s18e3271]]
-[[Salome-Mecaの使用法解説:http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2]] (先頭が07のpdfが参考になるかも)
-[[Salome-Meca非線形解析:https://sites.google.com/site/codeastersalomemeca/home/code_aster-1/kouzou-hisenkei]] (古いver.みたいだけど、参考になるかも)
-[[code-asterの説明:https://code-aster.org/V2/doc/default/en/index.php?man=commande]] (asterstudyで書かれているコードの説明が書かれている。ただ英語で書かれているので、翻訳する必要あり)

トップ   編集 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS