複合材の構成材料の物性値を逆解析により同定する手法 1

本解析では複合材料および複合材料に使用されている樹脂の材料物性値が既知であり、繊維の材料物性値が未知の場合に、リバースエンジニアリングにより繊維の材料物性値を同定することを目的とします。
具体的には、均質化手法を用いたマルチスケール解析ソフトであるMultiscale.Simの線形均質化機能を用いて、繊維のヤング率及びポアソン比を変数として、複合材料の異方性弾性物性値を求めます。この時、既知の複合材料の異方性弾性物性値に合うように繊維のヤング率及びポアソン比をAnsys DesignExplorerを用いた最適化計算により同定します。

本解析で用いる複合材料のミクロモデルは、図1に示すような繊維の体積含有率50%の一方向強化材料を対象とします。また、既知である複合材料および樹脂の材料物性値をそれぞれ表1および表2に示します。なお、本稿では複合材料の実機の試験は実施できなかったので、表1に示す複合材料の物性値については、表3で示す繊維の物性値を用いて予めMultiscale.Simの線形均質化解析によって取得したものを使用します。つまり、繊維の物性値を未知数として最適化解析を実施し、表3に示す物性値が最適値として計算されれば精度よく最適化ができたと判断できます。

解析設定としては、図2に示すようにAとBのシステムを構築します。システムAはMultiscale.Simを使用した線形均質化解析、システムBはDesignExplorerによる最適化解析です。まず、システムAでは繊維のヤング率およびポアソン比を入力変数としてMultiscale.Sim均質化解析を実施し、複合材料の異方性弾性物性値を計算します。ここでは均質化結果のX方向ヤング率(Ex)、XY方向せん断弾性係数(Gxy)およびXY方向ポアソン比(PRxy)を出力変数として設定します。次にシステムBで最適化解析の設定を実施します。ここでは直接最適化を用いて、システムAより計算された複合材料の物性値が、表1と一致するように最適化を実施します。今回は繊維のヤング率の範囲を5e10[Pa]から1e11[Pa]、ポアソン比の範囲を0.15から0.25と指定して最適化を実施しました。なお、物性値の範囲の見当がつかない場合は、幅広い範囲から最適値を1度に探すのは難しいので、まずは応答局面を作成してある程度の範囲の見当をつけ、それから最適化解析を実施するという方法もあります。

最適値として得られた繊維の物性値を表4に、その時の複合材料の物性値を表5に示します。表4より繊維の物性値は表3と比較して非常に近い値が得られています。また、表5に示した最適化結果における複合材料の物性値も、表1と比較して非常に近い値になっており、精度よく繊維の材料物性値が同定できたと言えます。

Ansys Workbench Mechanical 2020R1
Multiscale.Sim 2020R1