共回転定式化
共回転定式化の理論基礎
共回転定式化とは
先生、「共回転定式化」って何ですか?
共回転(co-rotational)定式化は梁やシェルの大回転問題に最適な手法。要素ごとにローカル座標系を回転に追従させ、要素内では小変形の理論を使う。
要素が回転しても、要素内は「小さな変形」として扱える?
そう。全体的には大回転でも、各要素が十分小さければ要素内の変形は微小。「大回転=多数の微小回転の累積」と捉える。
大回転の問題
3次元の有限回転は非可換(回転の順序で結果が変わる)。これが大回転の定式化を複雑にする。
回転の順序が結果に影響する…小さな回転なら無視できるけど、大回転では重要。
共回転定式化は回転マトリクス $[R]$ を各要素で管理し、回転の非可換性を正しく扱う。OpenSeesの非線形梁要素やAbaqusのB31要素は共回転ベース。
まとめ
要点:
- 要素ごとにローカル座標系が回転に追従 — 要素内は小変形
- 3次元の大回転は非可換 — 回転の順序が結果に影響
- 梁・シェルの大変形に最適 — フレーム構造の崩壊解析
- OpenSees, Abaqus B31 — 共回転ベースの梁要素
コロテーショナル座標の発明とArgyris
コロテーショナル定式化は1960〜70年代にArgyris(Stuttgart大)とWempner(Georgia Tech)が独立に発展させた。各要素が自身とともに「回転・並進する局所座標系」を持つことで、大変位・大回転を扱いながら局所変形を線形域に保てる。Argyrisは「わずかに変形した剛体運動」と表現し、大変形解析の幾何学的非線形定式化を大幅に簡素化した。
共回転定式化の数値計算手法
共回転の実装
共回転定式化のアルゴリズム:
1. 各要素の初期ローカル座標系を設定
2. 変形後の要素の節点座標から剛体回転を抽出(極分解)
3. ローカル座標系を剛体回転分だけ更新
4. 更新されたローカル系で微小変形の剛性マトリクスを計算
5. グローバル系に変換
「剛体回転の抽出」が核心ですね。
極分解($[F] = [R][U]$、$[R]$が回転、$[U]$がストレッチ)で剛体回転を分離する。
まとめ
コロテーショナル法の更新手順
コロテーショナル法では各荷重ステップで①現在配置の要素局所座標系更新②局所変位の計算③内力ベクトルの全体座標変換④接線剛性行列の組み立ての順に計算が進む。局所変位をベクトリアル回転ではなくスピノル代数で更新することで、大回転での数値誤差が1/10以下に抑えられる。Abaqusの梁・シェル要素はこの手法を使っている。
共回転定式化の実務適用
共回転の実務
鉄骨フレームの地震崩壊解析、海洋ライザーの大変形、柔軟ロボットの変形追跡で使用。
実務チェックリスト
薄板金属プレスのスプリングバック解析
プレス成形後の「スプリングバック」(弾性回復による形状変化)はコロテーショナル大変形解析の典型適用例だ。高張力鋼(980MPa級)のプレスではパンチ離型後に最大5〜10mmのスプリングバックが生じ、寸法精度に直接影響する。Toyota・Honda・Subaruはいずれもコロテーショナル定式化ベースの板成形解析(Autoform・PAM-STAMP等)でスプリングバック予測を設計の標準としている。
共回転定式化のソフトウェア比較
共回転のツール
選定ガイド
ANSYS Beam188大回転解析の精度
ANSYSのBEAM188はTimoshenko梁理論+コロテーショナル定式化を採用し、細長比L/D>10の梁で大変形・大回転解析が高精度に実行できる。KEYOPT(2)=2の全ニュートン法設定で、90°以上の回転を含む問題でも収束する。ロボットアームの到達点精度解析や展開宇宙構造(折り畳まれた太陽電池パドル展開)の設計検証にも使われている。
共回転定式化の先端研究
共回転の先端研究
コロテーショナル法と一貫性シェル定式化
従来のコロテーショナルシェルは局所座標での「薄板要素」を使うため厚板効果(Timoshenko)が不正確だった。2000年代に開発された「一貫性コロテーショナル(consistent corotational)」定式化では局所変形に完全な3次元シェル理論を適用し、厚板でのせん断変形も正確に再現できる。ANSYSのSHELL181要素はこの定式化を採用している。
共回転定式化のトラブル対応
共回転のトラブル
コロテーショナル解析での過大なドリフト変形
コロテーショナル法で荷重ステップを大きくすると、節点の回転更新でドリフト(回転の累積誤差)が生じることがある。特に局所座標の更新でクォータニオンの正規化を省略するとドリフトが10ステップで1〜2°に達することも。各ステップで回転更新量の正規化チェックを行い、許容値(1e-6以下)を超えたらステップを分割して再計算することを推奨する。
関連トピック
なった
詳しく
報告