FEABeRjp / FEABeR

https://feaber.readthedocs.io/ja/latest/
Creative Commons Attribution 4.0 International
3 stars 0 forks source link

📝 Cook’s membrane を追加 (GetFEM) #20

Closed tkoyama010 closed 2 years ago

tkoyama010 commented 2 years ago

GetFEMを使ったCook’s membrane を追加します。main のドキュメント構造が整備されたらそれに合わせて変更をします。

tkoyama010 commented 2 years ago

Draft設定にしてあるので現時点ではマージできません。変更を加えたら再度レビューリクエストを送りますね。

tkoyama010 commented 2 years ago

Cook's membraneを解いた論文です。 https://www.researchgate.net/publication/269705642_Mixed_Galerkin_and_Least-Squares_formulations_for_Isogeometric_analysis

tkoyama010 commented 2 years ago

こちらも参照できます。( @mmer547 様より情報をいただきました。) https://www.academia.edu/7059312/A_class_of_mixed_assumed_strain_methods_and_the_method_of_incompatible_modes

tkoyama010 commented 2 years ago

こちらの論文には参照値23.9の問題があります。 https://journal.augc.asso.fr/index.php/ajce/article/view/ajce.34.1.9

tkoyama010 commented 2 years ago

こちらの論文に解析解がありました。 https://www.researchgate.net/publication/265733095_A_mixed_three-field_FE_formulation_for_stress_accurate_analysis_including_the_incompressible_limit

tkoyama010 commented 2 years ago

解析解はないのでメッシュを細かくした収束値が既存の論文と整合しているか確認しましょう。

tkoyama010 commented 2 years ago

こちらの論文には三角形メッシュも含めてより詳しい値が掲載されていました。 https://zh.booksc.eu/book/3082842/fe9b18

tkoyama010 commented 2 years ago

こちらも詳しそうです。 http://symech.fgg.uni-lj.si/Library_10.0/Tests/Cook.html#

tkoyama010 commented 2 years ago

勉強会の方から頂いた情報でCook's membraneの原本の論文が見つかりました。 https://cybertycoons.files.wordpress.com/2014/04/robert_d-_cook_david_s-_malkus_michael_e-_pleshbookos-org-fem.pdf

tkoyama010 commented 2 years ago

https://www.dealii.org/current/doxygen/deal.II/code_gallery_Quasi_static_Finite_strain_Compressible_Elasticity.html

mmer547 commented 2 years ago

これモデル横方向44㎜でやられてますか?

mmer547 commented 2 years ago

Abaqusでやっているやつはこんな感じです image

tkoyama010 commented 2 years ago

GetFEMの結果です👇 displacement

mmer547 commented 2 years ago

AbaqusSEでも同じくらいの結果になりました。荷重の与え方がうまくいってなかったです。 image

tkoyama010 commented 2 years ago

@mmer547 よかったです。ところで左端の灰色になっているところはどのような状況でしょう?

mmer547 commented 2 years ago

表示バグっぽいですね。他のコンタは連続で出てるのですが、コンタのU2だけさっきの感じになってます。 image

mmer547 commented 2 years ago

拘束条件でY方向は0にしてるんですが、数値拾うと‐38乗のオーダーで値が拾えるので、最小値以下のコンタに割り当てられる色だと思います。 image

tkoyama010 commented 2 years ago

ありがとうございます。納得しました。

tkoyama010 commented 2 years ago

@mmer547 ちなみに最大値はいくらでしょうか?