ガスケット要素 — CAE用語解説
ガスケット要素
先生、フランジのボルト締結解析をしているんですけど、ガスケット部分のモデル化がうまくいかなくて。「ガスケット要素」を使えって言われたんですが、普通のソリッド要素と何が違うんですか?
理論と物理
ガスケット要素の物理的役割
ガスケット要素って、普通のソリッド要素やシェル要素と何が根本的に違うんですか?単に薄い材料をモデル化するだけならシェルでいい気がするのですが。
根本的な違いは、ガスケット要素が「面圧依存性」を持つ特殊な界面要素だという点です。シェルは面内の応力・変形を扱いますが、ガスケット要素は主に面に垂直な方向の圧縮と、それに伴う漏れ防止性能のモデル化が目的です。例えば、フランジ間に挟まれた非金属ガスケットは、締付け力が小さいと弾性変形のみ、大きくなると材料が塑性変形して密封面に食い込み、さらに大きくなると材料が破壊します。この非線形な挙動を定義するのがガスケット要素です。
面圧依存性というのは、具体的にどういう数式で表されるんですか?フックの法則みたいな単純な線形関係ではないですよね?
その通り、高度に非線形です。多くの商用ソフトでは、ガスケットの圧縮応力-ひずみ関係を実験データから定義します。一例として、AnsysのGasket要素で使われるモデルの骨子は、荷重-閉そく曲線と、それに連動した漏れ評価です。圧縮時の応力
「漏れ評価」まで計算できるんですか?流体の流れはCFDの分野だと思っていました。
直接的な漏れ流量の計算はしませんが、漏れに直結する「接触圧力」の分布と最小値を計算します。実務では、ROTT (Room Temperature Operational Tightness) テストに基づく経験式が使われることがあります。例えば、ガスケットメーカーのGarlockやFlexitallicが提供する設計ガイドでは、特定のガスケット材料に対して、座面応力から漏れ率を推定するチャートを用意しています。CAEでは、解析で得られたフランジ全面の最小接触圧力が、その材料の要求される最小座面応力(例:スパイラルワウンドガスケットで
数値解法と実装
要素定式化と材料データ入力
ガスケット要素は、有限要素法のマトリックスの中ではどのように扱われるんですか?厚さ方向のバネみたいなものですか?
良い喩えです。多くの実装では、厚さ方向の非線形バネと、面内方向の(通常は低い)剛性を組み合わせた「界面要素」として定式化されます。Abaqusでは「Gasket」要素種別が用意されており、厚さ方向の挙動は「圧縮-閉そく」テーブルで、せん断方向の挙動は「せん断応力-すべり」テーブルで定義します。剛性マトリックスは、この定義に基づいて接線剛性法などで逐次更新されます。重要なのは、引張剛性がゼロまたは非常に小さいこと。現実のガスケットが面から離れる(剥離)挙動を再現するためです。
材料データの「圧縮-閉そく」テーブルは、どうやって準備すればいいですか?自分で実験する必要があると、かなりハードルが高そうです。
まずはガスケットメーカーから提供されるデータシートを探すのが第一歩です。例えば、非アスベストシートガスケット材料の代表的な規格であるJIS R 3453やASTM F 104には、一定の材料クラスとその特性値が規定されています。より精密な解析には、メーカーがCAE用に提供するライブラリを使います。VictaulicやFrenzelitなどは、自社製品のAnsysやAbaqus用ガスケット材料ライブラリを公開しています。それでもデータがなければ、ASTM F 36(ガスケット材料の圧縮特性試験方法)に準拠した単軸圧縮試験を外注するケースもあります。初期設計段階では、公開論文のデータを参考にすることも多いです。
要素のメッシュサイズは、接するフランジのメッシュに合わせる必要がありますか?ガスケットは薄いので、1層の要素でいいんですよね?
基本的には1層で十分ですが、メッシュサイズには注意が必要です。接するフランジ面のメッシュと節点位置を一致させる「マッチドメッシュ」が理想です。一致させない場合、剛性マトリックスの連成が適切に行われず、応力集中がでたらめな場所に発生したり、収束性が悪化したりします。Abaqus/CAEの「メッシュ制約」機能や、Ansys Workbenchの「Mesh Connection」でボルト穴周りなどの詳細形状を揃えることが推奨されます。また、ボルト穴近傍の応力勾配が大きい領域では、フランジと同程度の細かさでメッシュを切る必要があります。
実践ガイド
フランジ接合部解析のワークフロー
実際にフランジとガスケットの組み立て解析をやろうとすると、ボルトの締め付け順序も再現する必要がありますか?それとも一括で荷重をかけますか?
締め付け順序は、特に大型フランジや不均一な座面応力分布が問題となる場合、再現した方が良いです。実務では「2ステップ法」が一般的です。第1ステップでボルト荷重(または締付けトルクから換算した軸力)を適用し、第2ステップで内圧などの作動荷重を加えます。Ansysでは「Bolt Pretension」要素、Abaqusでは「Bolt Load」機能を使い、初期締め付け力を定義します。より精密を期す場合は、ASME PCC-1 Appendix Oに記載の「交互締め付け順序」(例えば、対角線上にあるボルトから順に締める)を複数の荷重ステップで模倣します。ただし、計算コストと精度のトレードオフを考慮します。
解析が終わった後、結果のどこを見て「合格」か「不合格」を判断するんですか?
主に3点をチェックします。1. **最小座面応力**: ガスケット全体、特に内径側で、要求される最小値(材料シートに記載)を下回っていないか。2. **最大座面応力**: ガスケットの圧潰強度(例:グラファイトシートで
ガスケット要素を使う解析では、収束しにくいと聞きました。安定して計算を終了させるコツはありますか?
収束性を改善する確実な方法はいくつかあります。まず、ガスケット材料データの「圧縮-閉そく」曲線の最初の勾配(初期剛性)を極端に大きく設定しないこと。現実的でなめらかな曲線にします。次に、Abaqusであれば「Stabilization」に微小な値を設定し、接触が成立するまでの剛体運動を制御します。また、荷重適用を急峻に行わず、「Ramp」状にゆっくりと増加させます。最も効果的なのは、ボルト締め付けステップを、最初は小さな荷重で接触を確立させ、その後で最終荷重まで増加させる「2段階アプローチ」です。Ansysでは「Auto Time Stepping」を有効にし、ソルバーが自動的に荷重増分を調整できるようにします。
ソフトウェア比較
各ソフトウェアでの実装と特徴
Ansys、Abaqus、COMSOLで、ガスケット要素の機能や使い勝手に大きな違いはありますか?
はい、アプローチがかなり異なります。
Ansys Mechanical: 「Gasket」は一種の特殊な接触要素として実装されています。2D/3D要素があり、圧縮のみならずせん断挙動も定義可能。材料データはテーブル入力で、Workbench環境では専用の「Gasket Material」データセットを使います。ボルトとの連成解析が直感的です。
Abaqus/Standard: ガスケット解析の老舗で機能が豊富。「Gasket」要素種別が明確に定義され、厚さ方向とせん断方向の挙動を独立して詳細に定義できます。「Gasket Behavior」というプロパティ設定が中心です。
COMSOL Multiphysics: 専用の「ガスケット要素」というより、「固体力学」インターフェース内の「関節」機能や「弾性基礎」モデルを組み合わせて等価的なモデリングを行うことが多いです。マルチフィジックス(熱-構造連成など)との結合が容易という強みがあります。
無料や低価格のCAEソフト(FreeCADのFEMモジュールやCalculiXなど)でもガスケット解析はできますか?
CalculiX(Abaqus入力形式互換のオープンソースソルバー)は、Abaqusと同様の「*GASKET」セクションを入力ファイルに記述することでガスケット解析が可能です。ただし、前処理GUIが貧弱で、材料データの入力や結果の可視化が商用ソフトに比べて難易度が高いです。FreeCADのFEMワークベンチでは、現状、専用のガスケット要素はサポートされていません。代替策としては、薄いソリッド要素に非線形弾塑性材料を定義するか、バネ要素を並べることで近似しますが、密封性能の評価までは困難です。実務設計では、信頼性とサポートの面からAnsysやAbaqusがほぼ標準です。
ソフトウェアごとに、推奨される要素タイプ(例えば、AnsysのSolid186とSolid187)はありますか?
ガスケット要素自体が独立した要素種別なので、周りの構造物に使う要素の選択が問題になります。Ansysでは、フランジに**Solid186**(20節点二次六面体)を使うことが多く、これは曲げ変形を精度良く捉え、ガスケット界面との連成も安定するためです。メッシュ生成が難しい部分には**Solid187**(10節点二次四面体)も使いますが、剛性が高めに出る傾向があるので注意。Abaqusでは**C3D20R**(20節点二次六面体、低積分)や**C3D10**(10節点二次四面体)が一般的です。共通のアドバイスは、ガスケットと接触する面で、要素次数を二次要素に統一すること。一次要素では接触圧力の分布が粗すぎて、実際より高い応力集中を生む可能性があります。
トラブルシューティング
よくあるエラーとその対策
解析中に「ガスケット要素が過度に圧縮された」というエラーが出て計算が止まります。原因は何ですか?
このエラーは、定義した「圧縮-閉そく」曲線の範囲を超えるひずみが生じたことを意味します。主な原因は2つ。1. **ボルト締め付け力が大きすぎる**: トルクから計算したボルト軸力が現実的か確認。例えば、M36ボルトで材質SCM440の場合、締付けトルクは約
内圧をかけるステップで、ガスケットが突然大きく開いて(引張方向に変形して)計算が発散します。これは何が起きているんですか?
これは「内圧による分離(Blow-out)」のシミュレーションそのものですが、数値的に不安定で発散している状態です。物理的には、ボルト締め付け力によるクランプ力が、内圧が生み出すフランジを開こうとする力(内圧 × 内径面積)に負けたことを意味します。対策は以下の順で試します。1. **接触安定化の増加**: 微小なダッシュポットを追加し、瞬間的な剛体運動を抑制。2. **荷重増分の縮小**: 内圧をかけるステップの初期増分を、最終値の1%など非常に小さく設定。3. **物理的な対策の検証**: 解析が正しいなら、設計そのものを見直す必要があります。ボルト本数の増加、ボルト径のアップ、あるいはより高い締付け力が許容される**金属ジャケットガスケット**などの使用を検討します。ASME規格では、必要な最小クランプ力の計算式が与えられています。
接触圧力のコンター図を見ると、ボルト位置とボルト間で極端にムラがあり、現実的ではない気がします。メッシュを細かくしても改善しません。
これは「ピーキング」と呼ばれる数値的な問題で、フランジの曲げ剛性が過大評価されている、またはガスケットの面内剛性が過小評価されていることが原因です。対策:1. **フランジモデルの見直し**: フランジと胴板(ハブ)の接合部にフィレットを付け、鋭利な角をなくす。現実の製造形状に近づけるだけで応力集中は緩和されます。2. **ガスケットの面内剛性**: ガスケット要素の面内剛性(横弾性係数)をゼロではなく、材料の実測値(例えば、グラファイトシートで
関連トピック
なった
詳しく
報告