FEM バー要素理論
局所座標系での剛性マトリクス:$k = \frac{EA}{L}\begin{bmatrix}1 & -1 \\ -1 & 1\end{bmatrix}$
部材力:$F_{bar}= \frac{EA}{L}(u_2 - u_1)$
応力:$\sigma = F_{bar}/A$
プラット・ワーレン・ハウ型トラス橋をFEMで解析。荷重をかけると各部材が赤(引張)・青(圧縮)で表示され、変形図がアニメーションされる。
トラス構造の各部材は、軸方向の力(引張・圧縮)のみを伝える「バー要素」としてモデル化されます。1本の部材の剛性(変形のしにくさ)は、断面積A、材質の硬さE、長さLで決まり、以下の「要素剛性マトリクス」で表されます。
$$k = \frac{EA}{L}\begin{bmatrix}1 & -1 \\ -1 & 1\end{bmatrix}$$ここで、$E$はヤング率(弾性係数)[Pa]、$A$は断面積[m²]、$L$は部材長さ[m]です。このマトリクスは、部材の両端(節点)の変位と力の関係を表しています。
全ての部材の剛性マトリクスを組み合わせて「全体剛性マトリクスK」を作り、支点条件を考慮した上で、以下の連立一次方程式を解くことで、各節点の変位uが求まります。
$$K \cdot u = F$$$K$は全体剛性マトリクス、$u$は全節点の変位ベクトル、$F$は全節点に作用する外力ベクトルです。変位uが分かれば、個々の部材の伸び縮みから部材力と応力を計算できます。
鋼道路橋の設計:実際のプラット・トラス橋の設計では、このFEM解析を用いて、大型トラックなどの活荷重に対して各部材に生じる応力を評価します。部材の断面積を最適化し、軽量化と強度確保の両立を図ります。
既存橋梁の健全性評価:老朽化した橋梁の補修・補強計画策定に利用されます。材料の劣化(Eの低下)や断面欠損(Aの減少)をモデルに反映させ、現在の荷重条件で安全かどうかをシミュレーションで確認します。
鉄道橋の動的解析への発展:トラス橋の基本静解析を発展させ、列車が走行する時の繰返し荷重(疲労)や振動(動的応答)を評価する基礎となります。特に連結部の詳細な応力集中評価に繋がります。
建築構造物への応用:橋梁だけでなく、大型スポーツ施設の屋根トラスや鉄塔、クレーンなどの骨組構造の設計にも全く同じ原理が応用されています。複雑な形状の骨組の力の流れを可視化するのに役立ちます。
このシミュレーターで遊んでいると、いくつか勘違いしやすいポイントがあるんだ。まず「断面積Aを大きくすれば、絶対に変形も応力も減る」と思いがちだけど、そう単純じゃない場合がある。例えば、プラットトラスの斜材のうち、引張力がかかっている部材(赤い部材)の断面積だけを極端に太くしても、橋全体のたわみはあまり減らないことがある。なぜなら、変形は「最も弱いリンク」に引っ張られるからだ。圧縮材(青い部材)や他の部材が柔らかいままなら、そこが変形のボトルネックになる。実務では、部材ごとに応力の大きさを見て、効率的に断面を決める「最適化」の考え方が必要になるよ。
次に、「シミュレーション結果の数値は、そのまま実物の安全判定に使える」という誤解。これは絶対にダメだ。このツールは「線形静解析」で、材料は永遠に弾性変形(元に戻る変形)しかせず、座屈も破壊も起こらない理想化されたモデルだ。実際の設計では、許容応力度や座屈耐力といった安全率を考慮する。例えば、計算で100MPaの応力が出ても、材料の降伏強度が235MPaなら、安全率を考慮して「許容応力度は140MPa」といった別の基準値と比較する。シミュレーション結果は判断材料の一つでしかないことを覚えておこう。
最後に、支点のモデル化の重要性。このツールでは橋脚との接合部は「ピン支承」(回転自由)として固定されている。でも実際の橋梁は、溶接やボルトでがっちり固定されていたり、ローラー支承で温度膨張を吸収したりする。支承の条件を「完全固定」に変えるだけで、部材に曲げモーメントが生じ、この単純なトラスモデルでは捉えられない応力が発生する。FEMで一番怖いのは「入力(境界条件)を間違えること」だから、モデルが現実のどの部分をどう単純化しているか、常に意識することがプロの第一歩だ。
このトラス解析の原理は、橋梁以外の様々なモノづくりの根幹を支えている。まず航空宇宙工学だ。飛行機の主翼や機体の骨組(フレーム)は、まさに複雑な3次元トラス構造と言える。軽量化が命だから、各部材の応力を精密に計算し、必要最小限の材料で設計する「トポロジー最適化」や「軽量トラス構造」の技術に直結する。例えば、人工衛星の展開トラスアンテナも、宇宙空間での熱膨張による応力を、このFEMで事前にシミュレーションしている。
次にロボティクス。産業用ロボットアームは、リンク(アーム)とジョイント(関節)で構成される。このリンク部分の変形を抑えて先端位置精度を高めるためには、剛性設計が必須だ。まさに「断面積A」と「材質E」をどう選ぶかという問題になる。さらに、アームの動きに伴って重心が変わるので、それによる動的荷重をトラスモデルに載せることで、ロボットの構造設計の基礎が学べる。
意外なところでは生体力学にも応用されている。人間の大腿骨などの骨は、内部が海綿骨というトラス状(梁構造)をしていて、効率的に体重を支えている。この構造を工学的に分析する「バイオメカニクス」では、骨の微小構造をバー要素の集合体としてモデル化し、力の伝達経路を解析する。つまり、自然界が長い時間をかけて最適化した「構造」を、我々が工学の言葉で理解するためのツールでもあるんだ。
このシミュレーターに慣れたら、次は「連立方程式 $K u = F$ の向こう側」を覗いてみよう。まず手を動かすなら、簡単な2部材や3部材のトラスを、自分でマトリクスを組み立てて解いてみるのが一番の近道だ。例えば、縦棒1本と斜め棒1本がくっついたL字型の簡単なトラスを考え、支点と荷重を設定する。要素剛性マトリクスを書き、全体剛性マトリクス$K$を組み立て、支点条件を適用して…というプロセスを体験すると、「シミュレーターがブラックボックス」ではなくなる。
数学的な背景としては、線形代数、特に行列とベクトルの演算、連立一次方程式の解法(直接法、反復法)の理解が役立つ。なぜならFEMの核心は、複雑な微分方程式を、最終的には$K u = F$という大きな行列問題に帰着させてコンピュータに解かせる技術だからだ。また、変分法やエネルギー原理(最小ポテンシャルエネルギーの原理)を学ぶと、「なぜ剛性マトリクスで表せるのか?」という根本的な理解が深まる。
次のステップのトピックとしては、「平面ラーメン構造の解析」がおすすめだ。トラスは「部材がピン接合で曲げモーメントが生じない」という特殊なケース。一般的な骨組構造では、接合部が剛接合で、部材には曲げとせん断も加わる。その解析には「梁要素」を使うことになる。トラス解析でFEMの流れを理解した後で梁要素を学ぶと、要素剛性マトリクスがより複雑になる理由や、節点の自由度(変位と回転)が増えることの意味が、対比しながらよくわかるはずだ。