動的解析手法

動的解析手法

本節では直接時間積分法を適用した動的問題解析手法について示す。以下に示すように、本開発コードでは、陰解法及び陽解法による時刻歴応答解析が可能である。

陰解法の定式化について

動的問題を対象として、下式に示す運動方程式の解法に直接時間積分法を適用した。

ここでは、は質量マトリクス,は減衰マトリクス、は内力ベクトル,は外力ベクトルである。 なお、本ソフトは質量の変化を考慮せず、質量マトリクスは非線形において変形によらず一定となる。

時間増分内での変位、速度及び加速度の変化は、Newmark-法を用いて式及び式に示すように近似している。

ここで、はNewmark-法のパラメータである.

よく知られているように、及びを以下の値にした場合、線形加速度法あるいは台形則に一致する。

(線形加速度法)

(台形則)

、式を式に代入すると次式が得られる。

特に、線形問題に対してはは線形剛性マトリクスとし、となり、この式を上式に代入すると次式が得られる。

なお、幾何学的境界条件として加速度が指定されている箇所では、式から次式の変位を得る。

同様に、速度が指定されている箇所では、式から次式の変位を得る。

ここで、 は時刻における節点変位、は時刻の節点速度、 は時刻の節点加速度、は節点自由度番号、は節点番号である. また、質量項及び減衰項の取り扱いは次のとおりとした。

質量項の取り扱い

質量マトリックスについては原則として集中質量マトリックスとして扱っている。

減衰項の取り扱い

減衰項については式で表されるRayleigh減衰として扱っている。

ここで、はRayleigh減衰のパラメータである。

陽解法の定式化について

陽解法では下式に示す時刻tにおける運動方程式を基にする。

ここでは、時刻及び時刻における変位を時刻$t$におけるTaylor展開により表し、に関する2次項までとると、次のようになる。

、式の差及び和から次式が得られる。

、式を式に代入すると次式が得られる。

特に、線形問題に対してはとなり、上式は以下になる

ここで、質量マトリックスを集中質量マトリックス、減衰マトリックスを比例減衰マトリックス のようにおくと、式は連立方程式の求解操作を不要とする。

従って、式からは次式により求めることができる。