|TOP Page|
有限要素法において複数の節点が自由な値を取ることができずに,まとめて何らかの制約を受けるような場合を考えよう.これは一般的に多点拘束,MPC(Multiple Point Constraint)と呼ばれる.ここではラグランジュ未定乗数法を用いたMPCについて説明する.
このような拘束条件が付加される場合には以下の3つのやり方が考えられる.
自由度の削減は,例えば固体の解析なら変位のxy方向拘束をする場合にx節点の自由度を削減して連立一次方程式を立てることに相当する.ペナルティ法は制約条件から離れてしまうと大きな力が発生するような硬いバネを入れて解く方法である.ラグランジュ未定乗数法は拘束力を未知数として一緒に解く方法である.自由度を削減するやり方は自由度が減りロバストなものの,使える機会が限られる.ペナルティ法は自由度の数はそのままでコーディングも楽だが反復法で連立一次方程式が解きづらくなる.また人工的なパラメータ(硬いバネのバネ定数)が必要であるし,厳密には拘束条件を満たさない.ラグランジュ未定乗数法は自由度が拘束力の分だけ増えるものの,大抵の拘束条件に応用できるので汎用性がある.行列が不定値になり,反復解法でやや解きづらくなるが,ペナルティイ法ほど重大ではない.
変位に対する制約条件が
という風に書けるとすると,このような問題は
というようにかける.
上の制約条件を満たす範囲内での極値問題をそのまま変分的に書くと,解は次を満たす.
但しは
における
の接空間である.
つまりは
を
において満たす
上の定式化では,の接部分空間
を考えなければいけなく,有限要素法では一般的に取り扱いが困難である.そこでこのような空間を考えなくてもよい方法がラグランジュ未定乗数を使った方法である.
上式から,解の周辺においてと
の接空間が一致することがわかる.よって,
と
が平行になる.
つまり,あるスカラー量を使って,
と書くことができる.このがラグランジュ未定乗数である.以上からこの問題の解は
の極値となることがわかる.つまり,
を解は満たすことがわかる.
関数の選び方は幾つも考えられるが
を満たすようにしておくと,
であるから,は制約条件の満たされない方向へのエネルギーの勾配の大きさに等しくなる.エネルギーの変位での微分は力であるから,この場合
は拘束力の大きさを意味している.このようにLagrange未定乗数に拘束力の大きさとしての意味づけをすることができるので便利である.
拘束条件が複数ある場合を考える
この場合,極値が満たす条件は以下のとおり.
但し,は
の接空間の積集合である.
つまり,に対して,
となる.
において,
の接空間と
の接空間の積集合は一致している.よって,それぞれの補空間や直交補空間は一致しているはずである(これは接空間が閉空間である場合にのみ言えることで,非圧縮性拘束条件などの場合は補間方法によっては成り立たない場合がある).直交補空間は補空間に含まれる.
は直交補空間で,
の和集合は補空間である.
よって,ある実数が存在して,
となる.これは次のポテンシャルエネルギーの極値となる.
有限要素法のラグランジュ未定乗数の例として,各節点ごとにベクトル値を持っていて,全体のポテンシャルエネルギーが節点のベクトル値の関数で表されていて,制約条件も各節点のベクトル値について個別に付加される状況を考える.
このときの制約条件化でのポテンシャルエネルギーの仮想変分を考えてみよう.制約条件が課されている節点について変分をとる.制約条件が課されていなければ取り扱いは同じである.
ここで,はそれぞれ,変位と拘束力についての等価節点内力で次の通りである.
これを微分することで接線剛性行列を考えてみよう
係数行列は対称となるこことがわかる.最終的に次の連立一次方程式を解くことで解が得られる
固体の解析にこのラグランジュ未定乗数による拘束を適応してみよう.
節点が直線上にあるという拘束を考える.例えば点を通り,単位法線ベクトルが
である直線より上にあるという条件は
と書ける.但し節点の変形後の位置がで変位が
で,変形前の位置が
である.
であるから,
のようになる.は拘束なしの剛性行列と同一である.
下の動画はラグランジュ未定乗数法を用いて,梁の左端の節点が赤い直線の上に来るように制約条件を付け加えた物.梁は線形弾性体を使い,右端には正弦波状の変位を与えている.
変形前の状態において,剛体の一点(普通は重心)がとして,そこから
だけ離れた点
に節点が固定されているとする.この場合に剛体上に節点が貼り付けられているとしよう.剛体の一点の変位が
として,剛体のこの点周りの回転が
と表されるとする.このとき剛体に貼り付けられた節点の位置
は次のようになる.
ここで,は回転を表し,
である.また,による1回微分,2回微分は
のようになる.
ここで,節点変位は次のようになる.
これは成分で書くと
最終的に拘束条件は,
となる.節点の座標を所定の位置に完全に固定してしまうので,ラグランジュ未定乗数はxy成分に必要である.さて,歪エネルギーポテンシャル関数をとして,この拘束条件を考慮した新たなポテンシャルは次のとおりになる.
但し,は
番目の制約条件の節点番号を表す.
このときの制約条件化でのポテンシャルエネルギーの仮想変分を考えてみよう.制約条件が課されている節点について変分をとる.制約条件が課されていない節点は制約条件が無い場合と取り扱いは同じである.
よって節点についての等価節点内力は次のようになる.
剛体の並進や回転についての最終的な内力は全ての拘束された節点からの内力の重ね合わせであることに注意したい.
接線剛性行列は内力を微分することで得られる.上三角行列は以下のとおりである.対象行列になるから下三角行列は上三角の転置をすればよい.
剛体の回転についての最終的な接線行列は全ての拘束された節点からの接線剛性の重ね合わせであることに注意したい.
非線形有限要素法のためのテンソル解析の基礎 |
久田俊明 著 |