超弾性体解析で現れる剛性行列の性質とその解法に関し.pptVIP

本発表の概要 Aが有限要素離散化においてどのようなプロセスを経て構成されるか Aの正定値性に関する議論 全体行列がどのようなプロセスを経て構成されるか 反復解法における前処理の構成法について 収束性について ゴム状物質の弾性ポテンシャル例 ポテンシャルWのLandscape(想像） まとめ 変位部剛性行列の安定性について 　　初期応力項が不安定性を引き起こす可能性あり 　　非圧縮性制約条件が大変形に対する安定性に貢献？ 　　Eを基準としない妥当な弾性ポテンシャルの定義法はあるか？ 線形化方程式の解法に関して 　行列のブロック構造および符号を考慮して前処理行列を構築することが重要。 　収束性の理論的評価の残された問題。 * * * * * * * * * * * * * * * * * * * * * 超弾性体解析で現れる剛性行列の性質とその解法に関して 鷲尾　巧　　久田　俊明 東京大学 日本応用数理学会「行列?固有値の解法とその応用」研究部会 第６回研究会 ほとんど話題は、不思議に思うこと、自分ではわかっていないことです!! ゴムと金属の違い 鉄鋼 2x106 ゴム 1～3 自然状態 伸ばす エントロピーSの減少 結晶構造 ゴム構造 ポテンシャルEの増加 ヘルムホルツの自由エネルギー ヤング率　(MPa) 鋳鉄 0.6x10－11 ゴム 53.7x10－11 圧縮率　(Pa－1) 弾性変形の限界の伸び 金属 1%以下 ゴム 1000% 統計力学的ミクロモデルにおける非圧縮性の関与（参照：ゴムの弾性（裳華房）　久保亮五　著） １ １ １ 伸ばす 高さ方向のエントロピー 水平方向のエントロピー 非圧縮的伸長 応力と歪みの関係式 連続体の変形を表すテンソル Lagrange座標 （変形前の位置を表す） Euler座標 （現配置での位置を表す） 変形勾配テンソル 右極分解　本質的な変形を取り出す R: 回転行列, U:正定値対称行列 歪みテンソル 右Cauchy-Green変形テンソル Green-Lagrangeひずみ 引張り 圧縮 せん断 弾性ポテンシャルWとその変分 歪みEを与えるために要する単位体積あたりのエネルギー 第１変分式 第２変分式 初期変位項 初期応力項 第２Piola-Kirchhoff応力テンソル 応力テンソルの物理的意味 参照配置 現配置 剛性行列の正定値性 初期変位項　: Wが凹関数ならば正定値 初期応力項　: Sが正定値ならば正定値 初期応力項の正定値性は不確定： Sが正定値　　?　主応力方向がすべって引張り 主不変量 等方的連続体の弾性ポテンシャルは主不変量の関数として表わされる。 不変量から定まる剛性行列の符号 第１不変量　　初期応力項：正定値　　初期変位項：ゼロ 第２不変量　　初期応力項：正定値　　初期変位項：ゼロ３負３ 第３不変量　　初期応力項：正定値　　初期変位項：ゼロ３負３ 第２、３不変量の剛性行列（初期応力項＋初期変位項）の符号判定は困難 あるいは など Mooney-Rivlin体 非圧縮性制約条件付き ペナルティ方的なポテンシャル 第２不変量の凹関数性を保証できないため、Newton-Raphson解法の収束性を 一般に保証できない。 境界および制約条件を満たす解の集合 初期状態 心筋の弾性ポテンシャルの例 初期応力項：　歪みに依存　　　初期変位項：　正定値 自然基底 異方性に従った基底 fiber sheet Usyk, Mazhari and McCulloch, Journal of Elasticity 61,2000 異方性連続体モデルの不安定性 　異方性連続体モデルは、初期変位項が正定値となるように設計される。 　しかし、安定性議論においては初期応力項の影響が無視されているようである。 　異方性モデルは、変形過程の初期において、解けなくなるケースが多々生じる。 　安定な第１不変量を少し加えることにより、不安定性を回避できる。 変位 変位勾配テンソル 歪みテンソル 弾性ポテンシャル 線形 非線形 （２次項あり） 凹関数 　歪みテンソルが変位の２次関数となっていることが問題。 　しかし、回転不変性などを考慮しつつ弾性ポテンシャルを変位勾配Zの凹関数として直接的に定義することは困難。 体積弾性の取扱について 未定乗数の導入 第１変分式　（制約条件式も付加） 第２変分式　（制約条件式も付加） 不定値か正定値悪条件か？ 未定乗数導入時 直接的変分 剛性行列 剛性行列 不定値 正定値　しかし悪条件 固有値分布 (1) #positive eigenvalues = dim (A) #negative eigenvalues =dim(C) c d 0 a b 固有値分布 (2)