微小変形線形弾性静解析

微小変形線形弾性静解析

ここでは微小変形理論に基づく弾性静解析についての定式化を示す。応力・ひずみ関係として線形弾性を仮定している。

基礎方程式

固体力学の平衡方程式、力学的境界条件、幾何学的境界条件(基本境界条件)は次式で与えられる(図2.1.1参照)。

ここで、は応力、は表面力、は物体力であり、は力学的境界、は幾何学的境界を表す。

固体力学における境界値問題(微小変形問題)

図 2.1.1 固体力学における境界値問題(微小変形問題)

微小変形問題におけるひずみ・変位関係式は次式で与えられる。

また、線形弾性体での応力・ひずみ関係式(構成式)は次式で与えられる。

ここでは4階の弾性テンソルである。

仮想仕事の原理

基礎方程式、式 、式と等価である、微小変形線形弾性問題についての仮想仕事の原理は次式のように表される。

さらに構成式を考慮して式は次式のように表される。

において、はひずみテンソル、は4階の弾性テンソルである。 ここで、応力テンソルとひずみテンソルを、 それぞれベクトル形式で と表すと、構成式 は次式のように表される。

ここで、 は弾性マトリクスである。

ベクトル形式で表された応力 および式 を考慮して、式は次式のように表わされる。

および式が、本開発コードにおいて離散化される仮想仕事の原理である。

定式化

仮想仕事の原理式を有限要素ごとに離散化して次式を得る。

要素ごとに、要素を構成する節点の変位を用いて変位場を次式のように内挿する。

このときひずみは、式を用いて次式のように与えられる。

、式を式に代入して、次式を得る。

は次式のようにまとめることができる。

ここで、 式、式で定義されるマトリクスおよびベクトルの成分は、有限要素ごとに計算し、重ねあわせることができる。

が、任意の仮想変位について成立することにより次式を得る。

一方、変位境界条件式は次式のように表される。

を拘束条件式のもとで解くことにより、節点変位を決定することができる。