|
外部接触条件を拘束条件と内部接触面の非貫入条件を課した勾配流モデルを用いて、弾性体の変形において変分問題と拘束条件としてモデル化し、FreeFEM と IPOPT パッケージを用いてその数値シミュレーションコードを開発した。また、裁断機による断裂のように移動する外部拘束条件を持つ場合について、さらに内部接触面の非貫入条件を課した2次元のき裂進展現象のシミュレーションを行い、ブレードの間隔がある値より大きい場合にき裂断面形状が大きく変形することがわかった。これらの結果を、国際会議ICCCM2023において発表した。また国際会議ICIAM2023においてミニシンポジウム"FreeFEM software package for finite element modeling of PDEs"をオーガナイズし、FreeFEMによってこのアルゴリズムをどのように再現したかについて発表した。
接触面における摩擦のモデル化についても検討を行った。現象をできるだけ取り込もうとする転がりすべり接触による摩擦熱の発生について考慮したモデルも提案されているが、柏原らによって研究が進められているTresca摩擦条件下での線形動弾性体方程式を元にモデル化できないか調査を始め、検討している。
更に、モデル化にあたり、き裂進展の方程式に立ち戻り、1次元平衡解の解の性質とフェーズフィールドモデルにおける時定数の物理的意味について検討を行った。前者については1次元複数き裂解のエネルギー比較から単一き裂解への遷移が予想されることを日本応用数理学会第20回研究部会連合発表会において発表した。後者については、時定数と粘弾性との関係性について現象・理論・シミュレーションで解析を行い、木村、田中らとプレプリントを発行した（雑誌投稿中）。
|
