NAFEMS LE5: Z断面片持ち梁のねじれ
ベンチマークの仕様と理論的背景
LE5は何を検証する問題か
NAFEMS LE5はシェル要素の検証問題で、Z形断面の片持ち梁の自由端にねじりの偶力を与えます。面白いのは、ねじり荷重なのに合否判定が根元付近の「軸方向応力」で行われることです。Z形のような開断面はサンブナンねじり剛性が極めて小さく、ねじりをそり(warping)の拘束——フランジが互いに逆向きに曲げられる機構——で受け持ちます。その結果、フランジには大きな軸方向の膜応力が立つ。この「開断面の拘束ねじり」という構造力学の要所を、シェル要素の膜挙動が正しく再現できるかを試すのがLE5です。
問題設定
| 項目 | 内容 |
|---|---|
| 形状 | Z断面の片持ち梁:長さ10 m、ウェブ高さ2 m、フランジ幅1 m、板厚0.1 m(薄肉・シェルでモデル化) |
| 材料 | \( E = 210 \) GPa、\( \nu = 0.3 \) |
| 荷重 | 自由端にねじりモーメント1.2 MN·m相当:上下フランジ端辺に一様分布のせん断力0.6 MNずつを逆向きに(偶力として)与える |
| 拘束 | 根元端面を固定 |
| 目標値 | 点A(根元近くの規定点)の中立面(膜)軸応力 \( \sigma_{xx} = -108 \) MPa |
なぜシェルの良問なのか
ねじりの問題なのに軸応力で判定って、ちょっと不思議です。板曲げの検証じゃないんですか?
そこがLE5の設計の妙でね。目標値が中立面の膜応力だから、検証されるのは板の曲げ性能ではなく膜としての面内挙動+断面全体の力の流れなんだ。シェル要素は「膜+曲げ+横せん断」の合成体で、要素定式化の差(低減積分・ドリリング自由度・ロッキング対策)は複雑な力の流れで初めて露呈する。Z断面のそり拘束ねじりは、フランジとウェブの間で膜力が受け渡される絶妙に複雑な流れを作る——単純な引張や純曲げでは差が出ない要素の実力が、ここで数字になるんだよ。
実施手順と読み取りの要点
標準的な実施手順
- モデル化 — 中立面でシェルモデルを作成(板厚0.1 mをプロパティで付与)。フランジ・ウェブの接続線の共有を確認
- 荷重 — 上下フランジの自由端の辺に、辺に沿って一様な分布せん断(エッジロード)として0.6 MNずつを逆向きに付与。合力・合モーメントを必ず検算
- 拘束 — 根元端面の全自由度固定
- メッシュ系列 — 仕様の推奨分割から始め、2倍・4倍で点Aの膜応力の収束を確認
- 読み取り — 点Aの中立面(membrane/mid-surface)応力の軸方向成分を抽出し、−108 MPaへの収束を報告
最頻出ミス——表・裏・中立面の取り違え
シェルの応力出力は上面(top)・下面(bottom)・中立面(mid)の3層があります。LE5の目標値は膜応力(mid)——板曲げ成分を含むtop/bottomを読むと、曲げ分だけ系統的にずれます。ツールの既定出力がどの層か(そして「表」がどちら側か=シェル法線の向き)を確認するのが、この問題の最初の関門です。実務でも「シェルの応力はどの層の値か」は報告書の必須記載事項——LE5はその規律を体に入れる問題でもあります。
荷重の等価置換の作法
「ねじりモーメント1.2 MN·m」を1点のトルクとして与えるのは誤りで、仕様通りフランジ端辺への一様分布せん断の偶力として与えます。集中トルクは局所変形と特異性を持ち込み、根元の応力にも影響します。等価置換の検算として、①付与した荷重の合力がゼロ(純偶力)、②合モーメントが1.2 MN·m、を荷重サマリで確認する——この「荷重を置換したら合力・合モーメントを検算する」習慣は、境界条件検証の荷重版として実務全般に効きます。
実務での活用法
シェル要素の受け入れ試験として
LE5は要素タイプ・定式化の比較に最適です。同じメッシュ密度で、①1次低減積分シェル、②1次完全積分、③2次シェル、を並べて点A応力の収束を比較すると、自社標準要素の実力(と粗メッシュでの癖)が定量化できます。ソルバー更新・要素設定変更(ドリリング剛性・砂時計制御パラメータ)の回帰確認にも、軽量なこの問題が向いています。T4・LE11と並べ、「熱・熱応力・シェル」の3点セットを検証台帳の基本枠にするのが推奨構成です。
構造設計への教訓——開断面はねじりに弱い
LE5の物理は実設計の重要教訓と直結しています。開断面(Z・C・L形)はサンブナンねじり剛性が閉断面の数十〜数百分の1で、ねじりをそり拘束の軸応力で受けるため、「ねじったら曲げ応力で壊れる」挙動をします。薄肉部材のねじり設計で断面を閉じる(箱にする)ことの効果、そり拘束部(根元・ダイアフラム)に応力が集中する理由——LE5を解いた経験は、これらを定量的な直感として与えてくれます。ベンチマークは検証のためだけでなく、構造力学の生きた教材でもあります。
梁理論(Vlasov)との突き合わせ
拘束ねじりはVlasovの薄肉梁理論(そり関数・バイモーメント)で手計算でき、シェル解と突き合わせると「シェルFEM vs 古典理論」の相互検証になります。根元近傍はサンブナンの原理が効かない(そり拘束の影響が支配的な)領域なので、理論との一致は完全にはなりませんが、応力のオーダーと分布傾向の確認には十分——「FEMの答えを別の理論体系で検算する」訓練問題としても価値があります。
ツール別の設定要点
設定の対比
| ツール | エッジ荷重 | 膜応力の出力 |
|---|---|---|
| Ansys Mechanical | エッジへのForce(分布指定) | Shell内の Membrane Stress / ミッドサーフェス指定 |
| Abaqus | *SHELL EDGE LOAD(分布せん断) | 断面点(SPOS/SNEG/中央)指定、または膜力SFから換算 |
| Nastran系 | エッジ節点への等価節点力(一様分布から換算) | 要素力(膜力)出力F1から σ=F/t で換算が確実 |
| シェル対応の汎用FEM全般 | 辺荷重機能がない場合は等価節点力(端部節点は半分) | 「Z1/Z2/mid」等の層指定を確認 |
等価節点力を手作りする場合の注意
辺への一様分布荷重を節点力に落とす場合、1次要素なら内部節点に等分・両端節点に半分(2次要素では中間節点に重みが乗る配分)という整合節点力の規則に従います。全節点に同じ値を与える「均等割り」は端部で分布が崩れ、局所応力を歪めます。ここも「荷重置換の合力・合モーメント検算」で守れる範囲ですが、分布形の崩れは検算をすり抜けるため、配分規則そのものを意識する必要があります。
シェル検証の現在
シェル定式化の系譜とベンチマークの役割
シェル要素は、せん断ロッキング・膜ロッキング対策(低減積分+砂時計制御、MITC系の混合補間、EAS法)の発展史そのものであり、LE5を含む標準ベンチマーク群(ピンチド円筒・Scordelis-Lo屋根・捩り板等の「シェル義務テスト」)がその審判役を務めてきました。新しい要素・ソルバーを評価するとき、この古典セットを一巡させるのが要素開発コミュニティの慣行で、ユーザー側も同じセットで「自分の使う要素の系譜と癖」を把握できます。
ソリッドシェル・中立面レス化の潮流
CADソリッドから中立面を抽出する工数を嫌い、薄肉をソリッドシェル要素や厚み方向1要素のソリッドで直接解く流れが強まっています。LE5のようなベンチマークは、「その簡略化・要素選択で膜挙動の精度が保てるか」を確かめる物差しとして、モデリング方針の意思決定にも使えます。中立面シェルとソリッドシェルでLE5を解き比べる社内スタディは、方針転換時の説得力ある根拠になります。
自動検証への組み込み
T4・LE11と同様、LE5も軽量ゆえCI・回帰テストに向きます。要素ライブラリの更新・既定パラメータ変更は膜挙動に静かに影響し得るため、「シェルの基準器」として±1%許容の自動判定を仕込んでおくと、環境更新の安心感が段違いになります。
トラブル対応
「−108 MPaにならない」ときの診断表
| 症状 | 考えられる原因 | 対策 |
|---|---|---|
| ±数十%ずれる | top/bottom応力を読んでいる(曲げ成分混入) | 中立面(膜)応力を指定して再抽出 |
| 符号が逆 | 点Aの取り違え(上下フランジ)、荷重の向き | 仕様の点A座標と荷重方向を再確認 |
| 倍・半分のずれ | 偶力の片側しか掛けていない/0.6 MNを両側合計と誤解 | 荷重サマリで合モーメント=1.2 MN·mを検算 |
| 粗メッシュで大きく外れ、収束も遅い | 1次低減積分の粗メッシュ癖、フランジ幅方向の分割不足 | 幅方向に最低4分割、2次要素で比較 |
| 局所的に応力が暴れる | 集中トルク・不正な節点力配分、根元の特異性 | 分布エッジ荷重へ。根元角の値でなく点A規定位置で評価 |
| ツール間で数%差が残る | ドリリング剛性・せん断補正・応力回復法の差 | 要素定式化の設定を明記して比較。収束値同士で判定 |
検証台帳「基本3点セット」の完成
T4、LE11、LE5と3つ通しました。この後は何を足していけばいいでしょう?
関連トピック
なった
詳しく
報告