オイラー・ベルヌーイ梁理論
理論と物理
最も基本的な梁理論
先生、「オイラー・ベルヌーイ梁」って何ですか？ オイラーは座屈でも出てきましたよね。
レオンハルト・オイラーとダニエル・ベルヌーイが18世紀に確立した梁の曲げ理論だ。構造力学で最も基本的な理論であり、FEMの梁要素の出発点でもある。
基本仮定
どんな仮定をしているんですか？
3つの基本仮定：
1. 平面保持の仮定 — 変形前に平面だった断面は、変形後も平面のまま
2. 直交の仮定 — 変形後も断面は梁軸（中立軸）に直交する
3. 微小変形 — 変形が十分小さい
仮定2が重要そうですね。断面が傾かないということは…
そう。せん断変形を無視するということだ。曲げだけで梁が変形し、せん断による断面の傾きはゼロ。これがオイラー・ベルヌーイ梁の最大の特徴であり、同時に最大の制約だ。
支配方程式
曲げの微分方程式を教えてください。
たわみ $w(x)$ に関する4階の常微分方程式：
ここで $EI$ は曲げ剛性、$q(x)$ は分布荷重。
4階の微分方程式！ 4回積分すれば解が出る？
そう。4つの積分定数は4つの境界条件で決まる。各端で2つの条件（変位と回転、またはせん断力とモーメント）が必要だ。
各微分の物理的意味：
| 微分 | 物理量 | 式 |
|---|---|---|
| $w$ | たわみ | |
| $w' = dw/dx$ | 回転角 $\theta$ | |
| $w'' = d^2w/dx^2$ | 曲率 $\kappa = M/(EI)$ | $M = EI w''$ |
| $w''' = d^3w/dx^3$ | せん断力 | $V = -EI w'''$ |
| $w'''' = d^4w/dx^4$ | 分布荷重 | $q = EI w''''$ |
たわみを4回微分するだけで荷重に戻る。エレガントですね。
この関係を暗記しておくと、FEMの結果の検証に非常に役立つ。たわみ→回転角→曲率→モーメント→せん断力→荷重、と連鎖する。
せん断変形を無視する影響
せん断変形を無視するとどの程度の誤差が出ますか？
梁のスパン/せい比（$L/h$）が大きいほど誤差が小さい：
| $L/h$ | せん断変形の寄与 | オイラー・ベルヌーイの精度 |
|---|---|---|
| > 20 | < 1% | 十分正確 |
| 10 〜 20 | 1 〜 5% | 実用上問題なし |
| 5 〜 10 | 5 〜 20% | 注意が必要 |
| < 5 | > 20% | 不正確。ティモシェンコ梁を使うべき |
$L/h < 10$ くらいから差が出始めるんですね。
一般的な鉄骨梁（H形鋼）は $L/h = 15 \sim 25$ 程度だからオイラー・ベルヌーイで十分。しかしサンドイッチパネルや短い連結梁（$L/h < 5$）ではティモシェンコ梁理論が必要だ。
FEMでの梁要素
FEMのオイラー・ベルヌーイ梁要素はどんな形ですか？
2節点の梁要素で、各節点に3自由度（2Dの場合）：たわみ $w$、回転角 $\theta$、軸変位 $u$。
重要な特徴は形状関数がエルミート多項式（3次多項式）であること。通常のFEM要素（ラグランジュ多項式）とは異なり、変位だけでなく回転角も節点変数にする。これにより、わずか2節点で曲げの4次の微分方程式を正確に解ける。
2節点で4次の方程式が解ける！ 要素1つで正確なんですか？
一様断面・一定分布荷重の場合、1要素で厳密解が得られる。これはオイラー・ベルヌーイ梁要素の大きな利点だ。集中荷重が作用する位置には節点を置く必要があるが、それ以外は粗いメッシュで十分だ。
まとめ
オイラー・ベルヌーイ梁理論を整理します。
要点：
- せん断変形を無視した古典的梁理論 — 断面は常に中立軸に直交
- $EI w'''' = q$ — 4階の常微分方程式
- $L/h > 10$ なら十分正確 — 細長い梁に適用
- FEMではエルミート補間 — 2節点で一様梁の厳密解
- $L/h < 5$ ではティモシェンコ梁を使う — せん断変形が無視できない
材料力学の授業で解いた片持ち梁のたわみ $\delta = PL^3/(3EI)$ は、まさにこの理論の解だ。FEMの梁要素はこの理論を離散化したもので、原理は同じなんだ。
オイラー梁理論の誕生
オイラー・ベルヌーイ梁理論は1744年にレオンハルト・オイラーが著した「De Curvis Elasticis」に起源を持つ。「断面は変形後も平面を保ち梁軸に直交する」という仮定は当初物議を醸したが、細長い梁（スレンダー比L/h>10）では現在も誤差1%以内の精度を発揮する。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
梁要素の剛性マトリクス
オイラー・ベルヌーイ梁要素の剛性マトリクスを教えてください。
2節点、各節点に曲げ自由度2つ（$w_i, \theta_i$）の場合、4×4の剛性マトリクス：
この $12, 6L, 4L^2$ の数字はどこから来るんですか？
エルミート形状関数を微分して $B$ マトリクスを作り、$\int_0^L B^T EI B \, dx$ を積分すると得られる。エルミート形状関数が3次多項式なので、曲率（2次微分）は1次式になり、積分は厳密に実行できる。
軸力の自由度（$u_i$）を追加すると、曲げと独立な2×2の剛性：
これを曲げの $[K_b]$ と組み合わせて6×6（2D）または12×12（3D）の梁要素剛性マトリクスを構成する。
等価節点荷重
分布荷重を梁要素にどう与えるんですか？
分布荷重 $q$ を等価節点荷重に変換する。一様分布荷重の場合：
各節点に $qL/2$ の力と $\pm qL^2/12$ のモーメント…これは両端固定梁の固定端モーメントですよね！
完璧な理解だ。等価節点荷重は固定端反力の符号を変えたものに等しい。この対応を知っていると、等価節点荷重が正しく計算されているか直感的に検証できる。
ソルバー別の要素名
| 要素 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 2節点梁（EB） | CBAR | B23（2D）, B33（3D） | BEAM3（2D）, BEAM4 |
| 2節点梁（Timoshenko） | CBEAM | B21（2D）, B31（3D） | BEAM188/189 |
NastranのCBARとCBEAMは違うんですか？
CBARはオイラー・ベルヌーイ梁（せん断変形なし）、CBEAMはティモシェンコ梁（せん断変形あり＋ワーピング対応）。CBEAMのほうが高機能だが、細長い梁ではCBARで十分。
AbaqusではB33がオイラー・ベルヌーイ、B31がティモシェンコ。3次元の梁要素ではB31（ティモシェンコ）がデフォルト推奨。B33を使うメリットは稀だ。
断面定義
梁要素の断面はどう定義しますか？
主要なパラメータ：
- $A$ — 断面積（軸力用）
- $I_y, I_z$ — 断面二次モーメント（主軸まわり）
- $J$ — ねじり定数（サン・ブナンねじり）
- $C_w$ — ワーピング定数（CBEAMのみ）
定義方法は2種類：
1. 直接入力 — $A, I_y, I_z, J$ の値を手動で入力
2. 断面形状から自動計算 — H形鋼、箱形、円形等の形状パラメータを与えて断面諸元を自動計算。Abaqusの *BEAM SECTION, Ansysの SECTYPE/SECDATA
手動入力のリスクは？
単位の間違いが最も多い。$I$ は長さの4乗だから、mm系なら mm⁴。mとmmを間違えると $10^{12}$ 倍ずれる。断面形状から自動計算するほうが安全だ。
まとめ
梁要素の数値手法、整理します。
要点：
- エルミート補間で厳密な曲げ解 — 一様梁は1要素で十分
- 等価節点荷重 = 固定端反力の符号反転 — 直感的な検証に使える
- CBARはEB梁、CBEAMはティモシェンコ — Nastranでの使い分け
- 断面諸元は形状から自動計算が安全 — 手動入力は単位ミスのリスク
- 3次元ではティモシェンコ梁が推奨 — EB梁を使うメリットは稀
エルミート補間の威力
オイラー・ベルヌーイ梁要素はたわみwと回転角θの両方に連続性が必要なため、3次エルミート多項式を形状関数に採用する。この4自由度定式化は1965年にTurnerらがNASAの宇宙構造物に適用し、わずか6要素で解析梁との誤差0.3%を達成したと技術報告書に記録されている。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
梁要素の実務適用
梁要素は実務でどう使われていますか？
梁要素は構造解析で最も多く使われる要素タイプだ。ビルのフレーム、橋梁、機械のシャフト、配管…どこにでも使われている。
フレーム解析
ビルの骨組み解析ではどうですか？
鉄骨造や鉄筋コンクリート造のフレーム解析は梁要素の最大の適用分野だ。柱と梁を梁要素でモデル化し、接合部は節点として処理する。
ポイント：
- 柱梁接合部 — 剛接合、ピン接合、半剛接合の区別
- 荷重の流れ — 荷重→床スラブ→梁→柱→基礎の荷重伝達
- ブレースの扱い — トラス要素（軸力のみ）またはバネ要素
床スラブの面剛性はどうモデル化しますか？
3つのアプローチ：
1. 剛床仮定 — 各階のノードを剛体リンクで結合。最も単純
2. シェル要素でスラブをモデル化 — 面剛性を正確に考慮
3. 斜材でブレース効果を入力 — スラブの面内せん断をブレースで代替
梁要素の結果の解釈
梁要素の結果はソリッド要素の結果と同じように読めますか？
全く違う。梁要素の結果は断面力（軸力 $N$、せん断力 $V$、曲げモーメント $M$、ねじりモーメント $T$）として出力される。応力は断面力から計算する：
応力は後処理で計算するんですね。
そう。多くのソルバーは梁要素の応力も出力するが、それは上記の式に基づく公称応力であり、局所的な応力集中（フィレット、ボルト穴、溶接部）は含まない。応力集中の評価にはシェルやソリッドのサブモデルが必要だ。
梁要素の限界
梁要素で解析できないケースは？
梁要素は「全体挙動の把握」に使い、詳細はシェルやソリッドで補う、という使い分けですね。
まさにそう。構造解析のワークフローは「梁要素のフレームモデル → シェルのサブモデル → ソリッドの局所モデル」という階層構造になることが多い。
実務チェックリスト
梁要素を使うときのチェックリストをお願いします。
断面の主軸方向の間違いは深刻そうですね。
H形鋼の強軸と弱軸を間違えると断面二次モーメントが10倍以上変わる。たわみも10倍変わるから、すぐに気づけるはずだが、意外と見落とされる。必ず変位のオーダーを手計算と比較すること。
橋梁設計での実務活用
本州四国連絡橋の明石海峡大橋（1998年完成）の設計解析では、主塔の初期剛性評価にオイラー・ベルヌーイ梁モデルが使われた。主塔高さ298mを約500要素に分割し、風荷重下での水平変位を数時間で試算。3次元シェルモデルとの変位誤差は約2%にとどまった。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
梁要素のツール
梁要素を使ったフレーム解析にはどんなツールがありますか？
汎用FEMの梁要素に加えて、構造設計に特化したフレーム解析ソフトがある。
構造設計専用ソフト
| ソフト | 地域 | 特徴 |
|---|---|---|
| SS3/SuperBuild | 日本 | 建築一貫構造計算。保有水平耐力計算まで |
| MIDAS Gen | 韓国/グローバル | 建築・土木兼用。各国コード対応 |
| ETABS | 米国 | ビルディングフレーム専用。CSIの旗艦製品 |
| SAP2000 | 米国 | 汎用フレーム解析。橋梁にも対応 |
| ROBOT Structural Analysis | Autodesk | Revitとの統合。BIMフロー |
| STAAD.Pro | Bentley | グローバルな一般構造解析 |
SS3は日本の建築設計でよく見ます。
SS3（超高層対応版がSuperBuild/SS7）は日本の建築基準法に完全対応した一貫構造計算プログラムだ。確認申請に必須の保有水平耐力計算を内蔵している。国内の建築設計事務所では事実上の標準だ。
汎用FEMと設計専用ソフトの使い分けは？
設計専用ソフトは設計コードチェックが自動化されている。梁の許容応力度、座屈耐力、たわみ制限…全てソフトが自動で判定する。汎用FEMは応力を出すだけで、コードチェックはエンジニアが手動で行う。
したがって：
- 建築の設計実務 → SS3, ETABS, MIDAS Gen（コードチェック込み）
- 非標準的な構造の詳細解析 → 汎用FEM（Nastran, Abaqus, Ansys）
- 研究・教育 → SAP2000, OpenSees（柔軟性と自由度）
オープンソース
オープンソースのフレーム解析ソフトはありますか？
選定ガイド
まとめると？
梁要素は最もシンプルな要素だけど、ツールの選択肢が最も多い分野ですね。
梁要素ベースのフレーム解析は建築・土木の日常業務だから、効率化の需要が大きい。ツール選びは「どの設計コードを使うか」で決まることが多い。
各ソルバーの梁要素実装
NASTRANのCBEAM要素は1963年のNASAオリジナル版から引き継がれた実装で、断面の剪断センターとニュートラル軸のオフセットを明示的に指定できる。一方、Abaqus B23要素は3次元Euler-Bernoulli梁として3節点（二次補間）を採用し、曲率の表現精度がCBEAMより高いとされる。
選定で最も重要な3つの問い
- 「何を解くか」：オイラー・ベルヌーイ梁理論に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
梁理論の先端研究
オイラー・ベルヌーイ梁は18世紀の理論ですが、まだ研究テーマがあるんですか？
梁理論自体は完成しているが、その拡張と新材料への適用が活発だ。
機能的傾斜材料（FGM）梁
断面方向に材料特性が連続的に変化するFGM梁の曲げ理論。中立軸が幾何学的中心からずれるため、通常の曲げ式 $M = EI\kappa$ が直接使えない。FGM梁の有限要素定式化が研究されている。
ナノ梁（非局所弾性理論）
ナノスケールでも梁理論が使えるんですか？
カーボンナノチューブやグラフェンシートの力学は非局所弾性理論（Eringen理論）に基づくオイラー・ベルヌーイ梁として記述できる。通常の弾性論に「非局所パラメータ」を追加して、分子間力の長距離効果を考慮する。
ここで $e_0 a$ が非局所パラメータ。6階の微分方程式になる。
18世紀の理論がナノテクノロジーに使われているとは！
確率論的梁解析
材料特性や荷重のばらつきを確率変数として扱う確率論的有限要素法（SFEM）の梁要素が研究されている。構造の信頼性評価に使われ、$E$ や $I$ のばらつきがたわみ・応力の分布にどう伝播するかを定量化する。
共回転定式化
梁の大変形はどう扱いますか？
共回転定式化（co-rotational formulation）が大変形梁の標準的手法だ。要素ごとにローカル座標系を回転させて追従し、要素内では小変形の梁理論を使う。大域的には大回転・大変位に対応できる。
OpenSeesの非線形梁要素やAbaqusのB31要素は共回転定式化に基づいている。地震時の鉄骨フレームの大変形崩壊解析に不可欠だ。
まとめ
梁理論の先端研究、まとめます。
18世紀の理論がナノからマクロまで、確定論から確率論まで、今でも構造力学の中心にある。
剪断ロッキング問題の回避
オイラー・ベルヌーイ要素をティモシェンコ梁として使うと、スレンダー比が大きくなるにつれ剛性行列が過剛になる「剪断ロッキング」が発生する。1981年にPrathapとBhashyamはReduced Integration法により2節点Timoshenko要素でもロッキングを排除でき、L/h=1000の細長い梁でも精度を維持できることを示した。
トラブルシューティング
梁要素のトラブル
梁要素の解析でよくあるトラブルを教えてください。
梁要素は設定がシンプルな分、間違いに気づきにくいことがある。
たわみが理論値と合わない
単純梁のたわみが $\delta = 5qL^4/(384EI)$ と一致しません。
確認項目：
1. $I$ の値は正しいか — 強軸/弱軸の間違い、単位の間違い（mm⁴ vs m⁴）
2. 支持条件 — 両端ピンのつもりが片端固定になっていないか。回転自由度の拘束を確認
3. 荷重 — 線荷重（N/m）と集中力（N）の混同。等価節点荷重が正しいか
4. メッシュ — 集中荷重の作用点に節点があるか（節点がないと荷重が正しく伝達されない）
5. せん断変形 — ティモシェンコ梁要素を使っていると、せん断変形分だけ追加のたわみが出る
5番は盲点ですね。ティモシェンコ梁でオイラー・ベルヌーイの理論値と比較したら合わないのは当然…。
その通り。AbaqusのB31（ティモシェンコ）でオイラー・ベルヌーイの理論解と比較すると、数%の差が出る。比較するならB33（オイラー・ベルヌーイ）を使うか、ティモシェンコの理論解と比較すること。
断面力図がおかしい
せん断力図が不連続になります。
集中荷重の作用点でせん断力が不連続になるのは正常だ。分布荷重でも、要素の境界で値が少し変わることがある。梁要素の断面力は要素単位で出力されるため、同じ節点でも左側要素と右側要素で値が異なる。
全体のBMD（曲げモーメント図）やSFD（せん断力図）はどう描けばいいですか？
各要素の端部の値を接続して描く。FEMのポスト処理ツールは通常これを自動で行うが、梁要素のBMD/SFD表示機能があるかどうかはソルバーのポストプロセッサによる。Abaqus/CAEは梁の断面力図を描ける。Nastranのf06ファイルからは手動でプロットする必要がある。
剛体リンクの扱い
梁要素と他の要素を接続するときの注意点は？
梁要素はたわみ $w$ と回転角 $\theta$ の両方の自由度を持つが、ソリッド要素は変位のみ。直接接続すると回転自由度が拘束されない。
対策：
RBE2とRBE3の違いは何ですか？
別のページで詳しく扱うが、簡単に言うと：
- RBE2 — 剛体結合。独立節点の動きが従属節点に剛体的に伝達
- RBE3 — 荷重分配。独立節点の力を従属節点に重み付きで分配。剛性を追加しない
座屈解析で座屈荷重がおかしい
梁要素で座屈解析したら、横座屈が出ません。
オイラー・ベルヌーイ梁要素（CBAR, B33等）はワーピング自由度がないため、横座屈（曲げねじり座屈）を正しく表現できない。横座屈を評価するには：
- ティモシェンコ梁要素（CBEAM, B31OS等） — ワーピング自由度あり
- シェル要素でモデル化 — 横座屈が自動的に出る
まとめ
梁要素のトラブル対処、整理します。
梁要素は設定項目が少ないからこそ、1つの間違いが結果を大きく変えるんですね。
その通り。断面諸元と支持条件の2つが正しければ、梁要素の結果は非常に信頼性が高い。この2つを確実に押さえることが梁要素解析の全てだ。
大変形時の線形梁の限界
線形オイラー梁要素は変位が梁長さの約5%を超えると幾何学的非線形の影響が無視できなくなる。実例として1980年代のNASAシャトル主翼桁の初期解析でこの仮定が適用され、熱変形評価で約8%の過小評価が生じたと1987年のNASA-TM-89865で指摘されている。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——オイラー・ベルヌーイ梁理論の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告