CAEで最大応力が収束しないのはなぜ?|応力集中と応力特異点の見分け方

CAEで最大応力が収束しない理由と、応力集中・応力特異点の見分け方を示す表紙

メッシュを細かくするたびに最大応力が上がる。解析は正常終了し、変形図にも大きな違和感はない。それでも最大応力だけが安定しない――構造解析では珍しくない現象です。

原因として応力特異点が疑われますが、「最大応力が上がったから特異点」とは断定できません。有限の応力集中をまだ十分に表現できていない場合、評価する位置が変わっている場合、要素形状や応力の平均化方法が結果へ影響している場合にも、最大値は動きます。

まず、比較条件を固定します。位置、評価量、要素形式、メッシュの変え方、後処理条件をそろえ、ピーク以外の変位・反力・応力分布と一緒に確認します。

この記事では、均質・等方な延性金属を対象とする線形静解析、小変形を前提に、最大応力を設計判断へ使えるかどうかを切り分けます。疲労、脆性破壊、き裂、溶接部、接触端、異材界面は評価方法が異なるため、同じ判断をそのまま当てはめません。

解析結果全体を点検する順番は、線形静解析の結果を検証する方法|最大応力より先に見る7項目にまとめています。本記事では、その確認を済ませた後にも残る局所ピークへ焦点を絞ります。

目次

最初に結論|応力集中と応力特異点は同じ「赤」でも意味が違う

応力集中と応力特異点は、どちらも小さな高応力域として表示されます。一枚のコンター図だけでは区別できません。モデルの成り立ちと、メッシュを変えたときの傾向を確認します。

応力集中と応力特異点を、形状とメッシュ細分化時の傾向で比較した概念図
応力集中と応力特異点では、メッシュ細分化時の傾向が異なる
状態数学モデル上のピークメッシュを細かくしたとき最大値の扱い
通常の応力集中有限条件が適切なら一定値へ近づく収束を確認して採用する
応力特異点上限を持たない特異点へ近い評価値が上がり、有限値への接近を確認できないその一点の最大値は採用しない
未解像・比較条件の不一致本来は有限の場合もある不規則に変わる、最大位置が移る、表示方法で大きく変わる条件をそろえて再確認する

特異点を含む有限要素モデルでも、一回の計算が返す応力値は有限です。要素寸法を小さくしていくと、理想化した連続体モデルの解が有限の上限へ近づかない、という意味で「特異」と呼びます。

逆に、メッシュ細分化の途中で最大値が上がっただけでは、まだ細分化が足りない有限の応力集中と区別できません。複数段階の傾向と、特異性を生み得る形状・荷重・拘束の有無を合わせて見ます。

「解析が収束した」と「最大応力が収束した」は別の話

CAE結果を判断するときは、次の三つを分けます。

  1. 計算の完了と数値解:計算が正常終了し、現在の離散化モデルの方程式が所定の精度で解かれたか。反復法を使う場合は反復収束も確認する
  2. メッシュ収束:メッシュを変えても、設計判断に使う評価量が許容範囲内で安定したか
  3. モデルの妥当性:形状、材料、荷重、拘束などが、解析目的に対して実物を十分に表しているか

正常終了だけで直接確認できるのは、主に一つ目です。粗いメッシュでも線形静解析は解けますし、現実には有限面積で伝わる荷重を一点へ与えたモデルも計算できます。

メッシュ収束は「選んだ数学モデルの解へ近づいた」ことの確認です。実物にある最小フィレットを省略したまま収束を調べても、実物の局所応力を確認したことにはなりません。NASA-STD-7009Bでも、数値モデルを正しく解いているかという検証と、実世界を適切に表しているかという妥当性確認を分けて扱っています。

変位法の有限要素解析では、主な未知量は変位です。ひずみは変位の勾配から、応力はひずみと材料則から求めます。このため、一般に応力は変位よりメッシュへの感度が高くなります。特異点が局所的なら、全体変位や離れた位置の応力が安定していても、角部の最大応力だけが動くことがあります。

したがって、「変位が安定したから最大応力も正しい」「最大応力が安定しないから解析全体が無効」のどちらも早計です。何を設計判断へ使うのかを先に決め、その評価量ごとに確認します。

応力特異点を生みやすい形状・荷重・拘束

半径ゼロとして表した凹角

切欠きや角穴など、材料側から見た凹角は、線形弾性モデルで特異性を生じることがあります。ただし、すべての鋭い角や90度角が必ず特異点になるわけではありません。角度に加え、自由端・固定端などの境界条件も関係します。M. L. Williamsの古典的な研究も、角部の特異性が境界条件に依存することを示しています。

図面で最小半径が指定された重要部を半径ゼロへ簡略化すると、本来は有限である応力集中を、上限のないモデルへ変えてしまう可能性があります。局所強度を判定する場所なら、製造上保証される最小半径をモデルへ戻す必要があります。

ソリッドモデルの一点・一線へ与えた荷重

実物の力は、接触面、ボルト座面、治具などの有限領域を通じて伝わります。これを連続体の一点へ集中させると、荷重点の近傍に特異応力を導入する場合があります。モデルの次元や要素形式によって挙動は異なるため、「点荷重はどのモデルでも同じ」とは扱いません。

関心領域が荷重点から十分に離れ、合力とモーメントによる全体挙動だけが目的なら、簡略荷重が役立つこともあります。一方、荷重導入部そのものの強度を評価するなら、実際の荷重伝達面や締結・接触のモデル化を検討します。

点拘束や拘束条件の端部

完全固定は、相手側を変形しない環境として扱う理想化です。一点だけの拘束、固定領域と自由領域が急に切り替わる端部、剛体領域と柔らかい領域の接続部には、局所的な高応力が現れることがあります。

対策は、固定端の最大値を機械的に除くことではありません。モデル範囲を延長して境界を関心領域から遠ざける、支持剛性をばね等で表す、荷重や反力を実際の取付面へ分散するなど、解析目的に合う境界へ置き換えます。

このほか、き裂先端、接触端、異材界面の端部・角部、ゼロ半径で理想化した溶接止端・ルート部にも、特異な応力場が生じる場合があります。ただし、これらは最大ミーゼス応力の除外だけでは評価できません。本記事の範囲外です。それぞれに適した評価法へ切り替えます。

メッシュ比較は位置・評価量・変更方法をそろえる

図2は典型的な傾向を説明する概念図で、実解析データではありません。実際の結果は単調に変化するとは限らず、要素品質や最大位置の移動によって上下することもあります。

メッシュ細分化で一定値へ近づく評価量と上がり続ける評価量を対比した概念グラフ
概念図(実解析データではありません)。最大値だけでなく、固定位置の応力や変位も併せて確認する

1.単位・荷重・反力・変形を確認する

単位系、荷重方向、合力と反力を確認します。変位や回転によって形状・荷重方向の変化を無視できない場合は、メッシュ以前に小変形解析の前提を見直します。

2.評価する領域と量を先に決める

モデル全体の最大値だけでは、別の拘束端へ最大位置が移ったときに比較できません。対象の角部、穴周辺、断面など、設計上の関心領域を定めます。延性金属の静的な降伏確認ならミーゼス応力が候補になりますが、荷重経路を理解するために応力成分や主応力も確認します。

3.特異点候補からの距離を実寸で固定する

「角から2要素離れた位置」のような決め方は避けます。要素寸法を変えると、2要素分の物理距離も変わるためです。角から同じ実距離にある点、同じ座標、同じパス上の分布を比較します。除外範囲も要素数ではなく、形状寸法や荷重導入範囲との関係で説明します。

4.同じ要素形式で段階的に細分化する

粗・中・細の少なくとも3水準を用意し、3点で傾向を判断できなければ段階を追加します。各段階の細分化率をおおむねそろえ、要素形式、要素次数、荷重、拘束、接触条件、評価位置、後処理設定は変えません。局所細分化を使う場合は、粗密の遷移部を関心領域へ重ねないようにします。

有限Rを一次要素で表す場合、細分化に伴って曲面の近似形状そのものが変わることがあります。これは、離散化誤差とは別の形状近似の影響です。元のR寸法を保ち、曲辺を表現できる要素または十分な周方向分割を使って、形状近似とメッシュ細分化の影響を切り分けます。

NAFEMSも、特定位置の重要な結果をメッシュ密度に対して追跡し、複数回の解析で収束曲線を確認する方法を示しています。ただし、何%なら十分かは解析目的と設計余裕によって変わります。収束判定の許容差は、設計判断に必要な精度から事前に決めます。

5.ピーク以外の量も同じ表へ記録する

記録する量確認したいこと
関心領域の最大応力ピークが有限値へ近づくか
固定した実距離・座標の応力特異点候補から離れた場が安定するか
パス上の応力分布高応力域の広さと勾配が安定するか
変位全体剛性の評価が安定するか
合力と反力荷重と境界条件の釣り合いが取れているか
平均化あり・なし、または積分点値後処理と要素間ばらつきの影響が大きくないか

変化率を使う場合は、比較式、分母、判定基準を報告書に残します。最大値だけ別の抽出方法に変えたり、最も都合のよいメッシュだけを採用したりすると、比較そのものが崩れます。

結果の組み合わせから「特異性の疑い」を判断する

一つの兆候ではなく、複数の結果の組み合わせで判断します。

観察された傾向考えられる状態次に行うこと
最大応力が一定値へ近づき、分布も安定する有限の応力集中を解像しつつある設定した収束基準と設計余裕を確認する
同じ箇所の最大値は上昇傾向だが、固定位置の応力・変位は安定する局所的な特異性が疑われる形状、荷重、拘束の原因と評価目的を確認する
最大位置が移り、固定位置の応力や変位も変わる未収束、要素品質、モデル条件の問題が残るメッシュとモデルを見直し、断定を保留する
平均化値だけ安定し、非平均化値や分布が不安定平滑化が粗いメッシュを隠している可能性積分点値、要素間差、局所メッシュを確認する

「最大値が上がり続ける」は特異点を疑う重要な所見ですが、それだけで確定とはしません。ピークが同じモデル化上の不連続に留まり、比較条件をそろえても有限値への接近が見えず、その一方で周辺や全体の必要量が安定している――ここまで確認して、局所的な特異性の可能性が高いと判断します。

平均化あり・なしの片方を「真値」にしない

一般的な変位型ソリッド要素の多くでは、応力を積分点で評価し、表示時に節点へ外挿・平均化します。ただし、出力位置と処理方法は要素やソフトの設定で異なります。共有節点には隣接要素から複数の値が届くため、それらを平均して滑らかに表示する機能があります。

平均化した図が滑らかでも、数学モデルの特異性が消えたわけではありません。反対に、平均化していない値が常に真値または安全側とも限りません。粗いメッシュでは、節点への外挿によって積分点値より大きな値が表示されることがあります。Abaqusの公式資料も、外挿した節点応力と積分点応力が大きく異なる場合は、急な応力変化を粗いメッシュで表している可能性を示しています。

ミーゼス応力のような不変量は、応力成分を平均してから計算するか、各要素で計算してから平均するかでも値が変わります。ソフト名だけで設定を推測せず、出力位置、外挿、平均化の範囲と順序を確認します。

実務では、同じ設定の平均化値と非平均化値を各メッシュで記録します。両者の差が細分化で小さくなることは解像度の参考になりますが、差が小さいことだけで収束や妥当性を証明するものではありません。材料境界や意図した不連続をまたいで平均化しないことも必要です。

架空のブラケットで比較実験を設計するなら

ここで示すのは実解析結果ではなく、確認手順を説明するための架空の実験設計です。

L字ブラケットの内側コーナーを対象に、二つの教育用モデルを用意します。一方は図面で保証される有限の内Rを持たせ、もう一方は比較のためR=0へ理想化します。関心部以外の形状、材料、分布荷重、境界条件、要素形式はそろえます。拘束端や点荷重が別の最大値を作らないよう、境界は関心部から離すか、荷重伝達を面で表します。

局所要素寸法だけを一定比率で細かくし、各段階で次を記録します。

  • コーナーを含む指定領域の最大ミーゼス応力
  • コーナーから同じ実距離に置いた複数点の応力
  • 同じパス上の応力分布
  • 同じ位置の変位と全体の反力
  • 平均化あり・なし、可能なら積分点での値

有限Rモデルのピークが一定値へ近づくか、R=0モデルではどの量だけが安定しないかを比較します。ここで得る数値は、実製品の強度値や一般的な許容値には流用できません。要素数や「何要素離す」といった普遍則を作るための実験でもありません。

公開記事へ実際の解析結果を示すグラフを追加する場合は、使用ソルバー、要素形式、メッシュ条件、抽出位置、平均化設定、元データを残し、概念図と区別します。数値やグラフを記事へ載せる前の確認方法は、製造業の技術記事で行う事実確認にもまとめています。

特異点が疑われた後は、解析目的で対応を分ける

特異点があっても、モデル全体が直ちに無効になるとは限りません。ただし、その一点の最大応力を降伏判定や安全率へそのまま使うことはできません。

入力確認からモデル修正または根拠付き評価へ進む応力特異点の判断フロー
特異点の疑いがあるときの確認と処置の流れ

局所強度が設計判断を左右する場合

製造上保証される最小R、実際の荷重伝達面、支持部の剛性をモデルへ反映し、有限の応力集中としてメッシュ収束を確認します。収束した線形弾性のミーゼス応力が降伏応力を超えるなら、線形解析だけで安全・破損を確定せず、材料の応力―ひずみ関係を用いた弾塑性解析や、目的に合う限界状態の評価を検討します。

境界の簡略化が関心領域から離れている場合

モデルを延長する、荷重・拘束を分散する、支持剛性を表す、サブモデルを使うといった方法で、関心領域への影響を確認します。理想化を残す場合は、除外した位置と理由、代わりに採用した評価量、その収束結果を報告します。

疲労・脆性材・き裂・接触端を扱う場合

「局所的だから無視する」と判断しません。疲労では実際の形状や表面、応力振幅が寿命へ影響します。き裂では応力拡大係数やJ積分など、目的に対応した破壊力学量を使います。接触端や異材界面も、対象分野の評価法へ切り替えます。

解析報告書には、採用しない値と採用する値を両方書く

最大応力を採用しない場合、コンターの表示範囲を変えて赤い部分を見えなくするだけでは説明になりません。例えば、次のように記録します。

対象角部はR=0として理想化している。比較条件をそろえて局所メッシュを段階的に細分化したところ、同一角部の最大ミーゼス応力は有限値への接近を確認できなかった。一方、指定した実距離の応力、対象位置の変位および反力は、本解析で事前に定めた変化許容範囲に収まった。このため、角部一点の最大応力は局所降伏判定に使用しない。本解析の判定範囲を[全体剛性/遠方部の応力など、案件固有の目的]に限定し、その目的に対応する[変位/指定位置の応力など]を採用した。角部の局所強度は本モデルでは未評価とし、必要な場合は、製造上保証される最小半径と実際の荷重伝達を反映して再解析する。

角括弧内と「変化許容範囲」は、案件ごとに決めて記入します。特異点という言葉だけで除外せず、原因候補、比較条件、採用した評価量、その評価量で判定できる範囲まで残すと、別の担当者が判断を追跡できます。

まとめ|最大値ではなく、設計判断に必要な量を収束させる

最大応力が動いたときに最初に行うのは、値を採用するか捨てるかの即断ではありません。

  • 同じ位置、同じ評価量、同じ後処理条件でメッシュを比較する
  • ピークに加え、固定した実距離の応力、分布、変位、反力を見る
  • 形状・荷重・拘束の理想化が局所強度の判断へ影響するなら、物理的なモデルへ戻す

応力特異点を適切に切り分けるとは、都合の悪い赤い点を消すことではありません。採用できない最大値と、代わりに設計判断へ使う収束量を明確にすることです。

解析担当者の頭の中にある判断手順を、社内資料や顧客向け記事として残したい場合は、サービス資料をご覧ください。案件固有の条件や公開範囲を含むご相談は、お問い合わせから承ります。関連する制作実績の一部は、ご相談内容を確認したうえで個別にご案内します。

参考にした資料

  • URLをコピーしました!
目次