有限要素法の1次要素、2次要素とは?同じ様に要素分割しても計算精度や計算時間が大きく変わることをフリーソフトで確認しよう!

zairiki
自宅で無料で使える有限要素法(FEM)
  • 自宅パソコンで無料で有限要素法を使ってみたい人向け
  • 有限要素法を使いこなして生活や仕事に活かそう!

有限要素法(FEM)の計算結果は、どの種類の要素を使うかどれだけ細かく分割するかによって大きく変わります。今回は、片持ち梁のたわみを例に、1次要素と2次要素で要素分割を変えた計算をして、結果への影響を確認しました。

この記事はYoutube動画で紹介した内容の概要です。詳細は動画をご覧ください。

1. 1次要素と2次要素の違い(復習)

FEMで形状を分割するときに使う要素には、1次要素2次要素があります。

1次要素と2次要素の特性比較表(形状・変形・計算コスト)
1次要素と2次要素の特性比較。1次要素は計算コストが小さいが曲げ変形に不向き。2次要素は曲げ変形も含めて精度が良いが、計算コストが大きい。
  • 1次要素:要素の角にのみ節点がある。形状・変形を1次式、応力・ひずみを0次(要素内で一定)で表す。直線形状や垂直・せん断変形には向くが、曲げ変形には不向き。計算コストが小さい。
  • 2次要素:辺の中間にも節点がある。形状・変形を2次式、応力・ひずみを1次(要素内で直線近似)で表す。曲線形状や曲げ変形も精度よく扱える。ただし計算コストが大きい

2. 今回扱う例題

1次要素が不向きとされる曲げ変形を例題として取り上げます。片方の端部を固定し、先端に集中荷重を与えた片持ち梁のたわみを計算します。計算には、FEMのフリーソフトCalculixとPrePoMaxを用います。

片持ち梁の例題:形状・材料・荷重・理論解の一覧表
例題の条件。長さ 100 mm・断面 10×10 mm、縦弾性係数 200 GPa・ポアソン比 0.3、先端荷重 10 N。梁モデルの理論解(最大たわみ)は FL³/3EI = 0.02 mm。

この問題を、要素の粗さを5段階に変えながら、1次要素と2次要素のそれぞれで計算し、理論解 0.02 mm と比較します。

3. 1次要素の計算結果

1次要素で要素分割の粗さを変えながら計算した結果を示します。

1次要素の要素数と最大変形量の関係グラフ(収束グラフ)
1次要素の計算結果。要素数を増やすにつれて最大変形量が理論値に近づくが、要素数が少ない条件では変形量が大幅に小さく(剛性が過大に)計算されている。

要素数が最も少ない条件では、変位の計算結果が理論解の わずか 17% とかけ離れた値になり、非常に計算精度が悪くなりました。1次要素は曲げ変形をせん断変形で近似するため、要素が粗いと誤差が非常に大きくなります。

4. 1次要素と2次要素の比較

要素数で比較すると

1次要素と2次要素の最大変形量と要素数の比較グラフ(理論値付き)
1次要素(青)と2次要素(赤)の比較(横軸:要素数)。2次要素は要素数が少なくても理論値(緑線)にほぼ一致する。

2次要素では、最も粗い条件(要素数 33)でも最大変形量は 0.0194 mm(理論値の 97%)と高精度です。要素数を増やすと少ない要素数でも早く理論値に収束します。

節点数で比較しても2次要素が有利

1次要素と2次要素の最大変形量と節点数の比較グラフ(理論値付き)
1次要素(青)と2次要素(赤)の比較(横軸:節点数)。節点数が同程度でも2次要素の方が少ない節点数で精度の良い結果が得られる。

1次要素で細かく分割した条件と2次要素で荒く分割した条件を比較するため、横軸を節点数で整理しました。節点数が同程度でも2次要素の方が精度が高く、曲げ変形への適性が確認できます。

5. 注意点:分布図だけでは計算精度の悪さに気づけない

変形量の分布比較:同じ要素数での1次要素(最大0.011mm)と2次要素(最大0.020mm)の色コンター図
同じ要素数(約500)での変形量分布の比較。1次要素(左)の最大値は約0.011 mmで理論値の半分以下だが、カラーコンターの見た目はどちらも「それらしい」分布になっている。

精度が悪い1次要素の結果であっても、変形量の分布コンターを見ると「それらしい」分布になっており、見た目だけでは精度が悪いことに気づきにくいです。有限要素法の計算結果を使うときは、理論解や近似解と比較して妥当性を検証することが不可欠です。

まとめ

  • 曲げ変形の計算では、1次要素を用いると変形量が実際より小さく(剛性が過大に)計算される。要素分割が粗いほど誤差が大きく、今回の最粗条件の変位は理論値の17%とかけ離れた値になった。
  • 2次要素は少ない要素数でも高精度な結果が得られ、曲げ変形を含む問題に適している。ただし節点数が多くなるため計算コストは大きい。
  • 精度の悪い計算結果であっても、分布図のカラーコンターの見た目はそれらしく見える。有限要素法を使うときは、理論解・近似解などと比較して計算結果を検証することが重要。

このチャンネルでは、材料力学を学んだり生活に役立てるための様々なコンテンツを用意しています。Youtubeではさらに詳しく解説していますので、ぜひあわせてご覧ください。

再生できない方はこちらからどうぞ

関連する内容

自宅で無料で使える有限要素法(FEM)のシリーズ

材料力学の教室

材料力学用語辞典

身近な材料力学

HOME

ABOUT ME
生活に役立つ材料力学
生活に役立つ材料力学
材料力学分野の仕事に20年以上従事する博士(工学)が運営しています。
2022年にYoutubeチャンネルを開始し、動画数が増えてきたので探しやすくするために本サイトを開設しました。
ここで紹介する動画以外にもYoutubeには多くの動画が有るので、ぜひチャンネル登録もお願いします。
記事URLをコピーしました