特許 6061982 異方性樹脂成形体の構造解析方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6061982

(24)【登録日】2016年12月22日

(45)【発行日】2017年1月18日

(54)【発明の名称】異方性樹脂成形体の構造解析方法

(51)【国際特許分類】

   B29C 45/76 20060101AFI20170106BHJP

   G06F 17/50 20060101ALI20170106BHJP

   B29K 105/12 20060101ALN20170106BHJP

【ＦＩ】

   B29C45/76

   G06F17/50 612J

   G06F17/50 680C

   B29K105:12

【請求項の数】11

【全頁数】30

(21)【出願番号】特願2015-91622(P2015-91622)

(22)【出願日】2015年4月28日

(65)【公開番号】特開2016-203584(P2016-203584A)

(43)【公開日】2016年12月8日

【審査請求日】2016年8月30日

【早期審査対象出願】

(73)【特許権者】

【識別番号】390006323

【氏名又は名称】ポリプラスチックス株式会社

(74)【代理人】

【識別番号】100106002

【弁理士】

【氏名又は名称】正林 真之

(74)【代理人】

【識別番号】100120891

【弁理士】

【氏名又は名称】林 一好

(72)【発明者】

【氏名】小林 憲一郎

(72)【発明者】

【氏名】青木 現

【審査官】 田代 吉成

(56)【参考文献】

【文献】 特開２０１４−２２６８７１（ＪＰ，Ａ）

【文献】 特開２０１４−１９９２１９（ＪＰ，Ａ）

【文献】 特開２０１１−７５８（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｂ２９Ｃ ４５／７６

Ｇ０６Ｆ １７／５０

Ｂ２９Ｋ １０５／１２

(57)【特許請求の範囲】

【請求項1】

異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルを作成する構造解析用モデル作成方法であって、
前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、
前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、
前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、
前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、
前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとを含む、異方性樹脂成形体の構造解析用モデル作成方法。

【数1】

（式（１）において、
ν_２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
ν_ｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
ν_ｍは、対象材料を構成する樹脂組成物のポアソン比であり、
Ｖ_ｆは、対象材料に含まれる異方性フィラーの対象材料に対する体積含有率である。）

【数2】

（式（２）において、
α_νは、対象材料の体積膨張率であり、
α_１１は、対象材料の主軸方向の線膨張率であり、
α_２２は、材料主軸方向に直交する方向の線膨張率であり、
α_３３は、材料主軸方向に直交する方向の線膨張率であり、
α_２２≧α_３３という関係を有する。）

【請求項2】

前記α_１１は、下記式（３）によって求められ、前記α_２２は、下記式（４）によって求められる、請求項１に記載の構造解析用モデル作成方法。

【数3】

（式（３）において、
α_ｓは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率であり、
α_ｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、上記式（３）’で規定され、
λ_１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（３）’において、α_νは、対象材料の体積膨張率である。）

【数4】

（式（４）において、
α_ｂは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率であり、
α_ｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、下記式（４）’で規定され、
λ_１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（４）’において、α_νは、対象材料の体積膨張率である。）

【請求項3】

前記物性情報は、下記式（５）で規定される弾性率をさらに含む、請求項１又は２に記載の構造解析用モデル作成方法。

【数5】

（式（５）において、
Ｅ_１１は、対象材料の主軸方向の弾性率であり、
Ｅ_ｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率であり、
ζ_１１は、式（５）’で示される値であり、
ηは、式（５）’’で示される値であり、
Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。）

【数6】

（式（５）’において、
λ_１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）

【数7】

（式（５）’’において、Ｅ_ｆは、対象材料に含まれるフィラーの弾性率であり、Ｅ_ｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。）

【請求項4】

前記物性情報は、式（６）で規定されるせん断弾性率をさらに含む、請求項１から３のいずれかに記載の構造解析用モデル作成方法。

【数8】

（式（６）において、
Ｇ_１２は、対象材料の主軸方向のせん断弾性率であり、
Ｇ_ｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率であり、
ζ_１１は、式（６）’で示される値であり、
η_ｇは、式（６）’’で示される値であり、
Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。）

【数9】

（式（６）’において、
λ_１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）

【数10】

（式（６）’’において、Ｇ_ｆは、対象材料に含まれるフィラーのせん断弾性率であり、Ｇ_ｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。）

【請求項5】

前記α_νは、ＰＶＴ樹脂特性解析によって得られる対象材料の体積膨張率の実測値である、請求項１から４のいずれかに記載の構造解析用モデル作成方法。

【請求項6】

前記異方性樹脂成形体がウエルド部を含む、請求項１から５のいずれかに記載の構造解析用モデル作成方法。

【請求項7】

前記探索ステップは、
前記複数の第１領域の各々について第１重心位置を導出する第１重心位置導出ステップと、
前記複数の第２領域の各々について第２重心位置を導出する第２重心位置導出ステップと、
前記第２領域の各々について、前記第２重心位置から最も近い第１重心位置を探索し、前記第２重心位置から最も近い第１重心位置を有する第１領域を前記第２領域から最も近い第１領域とする最短第１領域設定ステップとを有する、請求項１から６のいずれかに記載の構造解析用モデル作成方法。

【請求項8】

前記配向状態は、対象材料の主軸方向の配向度であり、
前記複数種類の仮想成形体は、前記配向度の範囲に応じて１０種類以上１０００種類以下の範囲で設定される、請求項１から７のいずれかに記載の構造解析用モデル作成方法。

【請求項9】

請求項１から８のいずれかに記載の構造解析用モデル作成方法を使用して異方性樹脂成形体を構造解析する構造解析方法であって、
前記異方性樹脂成形体について樹脂流動解析を行い、要素分割された前記複数の第１領域の各々における配向状態を算出するための配向状態算出基礎情報を含む第１領域配向情報を取得する第１領域配向情報取得ステップと、
前記第１領域配向情報取得ステップとは別に、前記異方性樹脂成形体の構造解析のために前記異方性樹脂成形体を前記複数の第２領域に要素分割する構造解析用要素分割ステップと、
前記第２領域配向状態の情報及び前記第２領域物性情報を含む構造解析モデル情報に基づいて前記異方性樹脂成形体の構造解析を行う構造解析ステップとを含み、
前記第１領域配向状態算出ステップは、前記配向状態算出基礎情報に基づいて前記複数の第１領域の各々における配向状態に相当する第１領域配向状態を算出するステップであり、
前記探索ステップは、前記構造解析用モデルにおける複数の第２領域の各々について、第２領域から最も近い第１領域を探索するステップである、異方性樹脂成形体の構造解析方法。

【請求項10】

請求項９に記載の構造解析方法を使用することで、前記異方性樹脂射出成形体の反り変形を予測する反り変形予測方法。

【請求項11】

異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルの作成をコンピュータに実行させるための構造解析用モデル作成プログラムであって、
前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、
前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、
前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、
前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、
前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとをコンピュータに実行させるためのプログラム。

【数11】

（式（１）において、
ν_２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
ν_ｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
ν_ｍは、対象材料を構成する樹脂組成物のポアソン比である。）

【数12】

（式（２）において、
α_νは、対象材料の体積膨張率であり、
α_１１は、対象材料の主軸方向の線膨張率であり、
α_２２は、材料主軸方向に直交する方向の線膨張率であり、
α_３３は、材料主軸方向に直交する方向の線膨張率であり、
α_２２≧α_３３ という関係を有する。）

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、異方性樹脂成形体の構造解析用モデル作成方法、構造解析方法及び構造解析用モデル作成プログラムに関する。

【背景技術】

【0002】

熱可塑性樹脂と異方性フィラーとを含有する異方性樹脂組成物を射出成形等することによって得られる異方性樹脂成形体は、工業製品として広く使用されており、数値解析による形状検討が多用されている。ところで、数値解析による樹脂成形体の強度設計には、弾性率、ポアソン比、線膨張率等の物性値が必要である。樹脂成形体が異方性フィラーを含有する場合、これらの物性値を実測すること、及び予測することのいずれも容易ではない。

【0003】

特許文献１に記載の弾性率計算方法は、対象物の繊維配向方向の第１弾性率と、繊維配向と直交する方向の第２弾性率と、構造解析を行う際の荷重方向とを予め格納しておき、対象物の任意の要素ついて繊維配向の向きと荷重方向とが平行している場合には第１弾性率を、繊維配向の向きと荷重方向とが直交している場合には第２弾性率をそれぞれ設定し、さらに、繊維配向の向きと荷重方向とが平行も直交もしていない場合には荷重方向と繊維配向の向きとの角度と第１及び第２弾性率のいずれかに基づいて弾性率を算出して弾性率を設定することを特徴とする。

【0004】

特許文献２には、各繊維単体の物性と含有繊維をマトリクスと等価な等価介在物として取り扱う為の仮想歪テンソルを加算演算することにより、複合材全体の物性を表現するテンソルを求めるという手法が紹介されている。

【0005】

特許文献３では、含有する強化繊維及びベースレジン単体の機械物性を使用して評価形状における流動配向解析を実施し、算出した樹脂成形品の配向パラメータより算出した機械物性を、実際に評価形状をした成形品の機械物性を測定した数値と比較し、計算値が実測定値と等しくなるようなベースレジン、および強化繊維単体の機械物性を同定し、同定した機械物性を使用して成形部品の変形量を高精度に評価する方法が記載されている。

【0006】

しかしながら、いずれの方法も、繊維配向による成形品の機械物性予測方法による予測値と実際の成形品の物性値には数値的な不一致が存在しており、さらなる改善が求められていた。

【0007】

この不一致の原因としては、計算にて求める場合の仮定と実際との違いが挙げられる。例えば、射出充填される際に異方性フィラーが破損して、計算に使用している異方性フィラーのアスペクト比の値と実際の値が異なってしまうこと、繊維含有率の不均一、異方性フィラー同士の干渉等による異方性フィラー含有時の流動配向の予測力不足が考えられる。すなわち、計算上の仮定が実現象を十分表現できないことに加え、各構成要素の物性値が不正確であることで、解析モデル構築時の誤差に入力物性値の誤差が加算されて複合された誤差が存在しており、フィラー配向計算結果から機械物性を予測しても、この予測では不十分であった。

【0008】

予測の精度を高めるため、特許文献４では、充填剤を所定の割合で含む樹脂組成物を成形してなる樹脂成形品についてＸ線ＣＴ測定等により、複数のスライス画像を取得し、スライス画像を微小面積の画素に分割し、その情報を元に樹脂成形品内の充填剤の配向状態を反映させた有限要素法モデルを作成し有限要素法解析を介して機械物性を同定し、成形部品の強度、変形量を高精度に評価する方法が記載されている。

【0009】

また、特許文献５では、フィラーを含む射出成形品内のフィラーの３次元配向のデータをＸ線ＣＴによる３次元画像から取得し、せん断応力の積分値および分子配向状態のデータに基づいて均質化法を介して線膨張係数を求め、構造解析により、射出成形品の熱間反りを求めている。

【0010】

しかしながら、特許文献４及び５に記載の方法では、各構成要素の物性値を得るために多大な実験、測定が必要である上、それを求める際の計算時間も多大となる。それにより、物性値を決定する際の計算に要するコストが膨大となっている。

【先行技術文献】

【特許文献】

【0011】

【特許文献1】特開平９−２３７２６７号公報

【特許文献2】特開平７−３０４０５６号公報

【特許文献3】特開２００４−２５７９６号公報

【特許文献4】特開２０１１−７５８号公報

【特許文献5】特開２０１４−１００８７９号公報

【発明の概要】

【発明が解決しようとする課題】

【0012】

本発明は、以上のような課題を解決するためになされたものであり、構造解析が容易ではないと考えられている異方性樹脂成形体を構造解析するにあたり、構造解析の精度を高めるとともに、解析結果を短時間で、かつ、コンピュータにかかる負荷を大幅に軽減可能な構造解析用モデルを提供することを目的とする。

【課題を解決するための手段】

【0013】

本発明者らは、上記課題を解決するために鋭意研究を重ねた。その結果、異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体を想定し、これら複数種類の仮想成形体の各々について、ポアソン比と線膨張率の少なくとも一方を含む物性情報を予め作成し、この物性情報を構造解析のために要素分割された複数の第２領域の各々に当てはめることで、上記の課題を解決できることを見出し、本発明を完成するに至った。より具体的に、本発明は、以下のものを提供する。

【0014】

［１］本発明は、異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルを作成する構造解析用モデル作成方法であって、前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとを含む、異方性樹脂成形体の構造解析用モデル作成方法である。

【数1】

（１）
（式（１）において、
ν_２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
ν_ｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
ν_ｍは、対象材料を構成する樹脂組成物のポアソン比であり、
Ｖ_ｆは、対象材料に含まれる異方性フィラーの対象材料に対する体積含有率である。）

【数2】

（式（２）において、
α_νは、対象材料の体積膨張率であり、
α_１１は、対象材料の主軸方向の線膨張率であり、
α_２２は、材料主軸方向に直交する方向の線膨張率であり、
α_３３は、材料主軸方向に直交する方向の線膨張率であり、
α_２２≧α_３３という関係を有する。）

【0015】

［２］また、本発明は、前記α_１１は、下記式（３）によって求められ、前記α_２２は、下記式（４）によって求められる、［１］に記載の構造解析用モデル作成方法である。

【数3】

（式（３）において、
α_ｓは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率であり、
α_ｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、上記式（３）’で規定され、
λ_１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（３）’において、α_νは、対象材料の体積膨張率である。）

【数4】

（式（４）において、
α_ｂは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率であり、
α_ｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、下記式（４）’で規定され、
λ_１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（４）’において、α_νは、対象材料の体積膨張率である。）

【0016】

［３］また、本発明は、前記物性情報は、下記式（５）で規定される弾性率をさらに含む、［１］又は［２］に記載の構造解析用モデル作成方法である。

【数5】

（式（５）において、
Ｅ_１１は、対象材料の主軸方向の弾性率であり、

Ｅ_ｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率であり、
ζ_１１は、式（５）’で示される値であり、
ηは、式（５）’’で示される値であり、
Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。）

【数6】

（式（５）’において、
λ_１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）

【数7】

（式（５）’’において、Ｅ_ｆは、対象材料に含まれるフィラーの弾性率であり、Ｅ_ｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。）

【0017】

［４］また、本発明は、前記物性情報は、式（６）で規定されるせん断弾性率をさらに含む、［１］から［３］のいずれかに記載の構造解析用モデル作成方法である。

【数8】

（式（６）において、
Ｇ_１２は、対象材料の主軸方向のせん断弾性率であり、
Ｇ_ｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率であり、
ζ_１１は、式（６）’で示される値であり、
η_ｇは、式（６）’’で示される値であり、
Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。）

【数9】

（式（６）’において、
λ_１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）

【数10】

（式（６）’’において、Ｇ_ｆは、対象材料に含まれるフィラーのせん断弾性率であり、Ｇ_ｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。）

【0018】

［５］また、本発明は、前記α_νは、ＰＶＴ樹脂特性解析によって得られる対象材料の体積膨張率の実測値である、［１］から［４］のいずれかに記載の構造解析用モデル作成方法である。

【0019】

［６］また、本発明は、前記異方性樹脂成形体がウエルド部を含む、［１］から［５］のいずれかに記載の構造解析用モデル作成方法である。

【0020】

［７］また、本発明は、前記探索ステップは、前記複数の第１領域の各々について第１重心位置を導出する第１重心位置導出ステップと、前記複数の第２領域の各々について第２重心位置を導出する第２重心位置導出ステップと、前記第２領域の各々について、前記第２重心位置から最も近い第１重心位置を探索し、前記第２重心位置から最も近い第１重心位置を有する第１領域を前記第２領域から最も近い第１領域とする最短第１領域設定ステップとを有する、［１］から［６］のいずれかに記載の構造解析用モデル作成方法である。

【0021】

［８］また、本発明は、前記配向状態は、対象材料の主軸方向の配向度であり、前記複数種類の仮想成形体は、前記配向度の範囲に応じて１０種類以上１０００種類以下の範囲で設定される、［１］から［７］のいずれかに記載の構造解析用モデル作成方法である。

【0022】

［９］また、本発明は、［１］から［８］のいずれかに記載の構造解析用モデル作成方法を使用して異方性樹脂成形体を構造解析する構造解析方法であって、前記異方性樹脂成形体について樹脂流動解析を行い、要素分割された前記複数の第１領域の各々における配向状態を算出するための配向状態算出基礎情報を含む第１領域配向情報を取得する第１領域配向情報取得ステップと、前記第１領域配向情報取得ステップとは別に、前記異方性樹脂成形体の構造解析のために前記異方性樹脂成形体を前記複数の第２領域に要素分割する構造解析用要素分割ステップと、前記第２領域配向状態の情報及び前記第２領域物性情報を含む構造解析モデル情報に基づいて前記異方性樹脂成形体の構造解析を行う構造解析ステップとを含み、前記第１領域配向状態算出ステップは、前記配向状態算出基礎情報に基づいて前記複数の第１領域の各々における配向状態に相当する第１領域配向状態を算出するステップであり、前記探索ステップは、前記構造解析用モデルにおける複数の第２領域の各々について、第２領域から最も近い第１領域を探索するステップである、異方性樹脂成形体の構造解析方法である。

【0023】

［１０］また、本発明は、［９］に記載の構造解析方法を使用することで、前記異方性樹脂射出成形体の反り変形を予測する反り変形予測方法である。

【0024】

［１１］また、本発明は、異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルの作成をコンピュータに実行させるための構造解析用モデル作成プログラムであって、前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとをコンピュータに実行させるためのプログラムである。

【数11】

（式（１）において、
ν_２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
ν_ｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
ν_ｍは、対象材料を構成する樹脂組成物のポアソン比である。）

【数12】

（式（２）において、
α_νは、対象材料の体積膨張率であり、
α_１１は、対象材料の主軸方向の線膨張率であり、
α_２２は、材料主軸方向に直交する方向の線膨張率であり、
α_３３は、材料主軸方向に直交する方向の線膨張率であり、
α_２２≧α_３３ という関係を有する。）

【発明の効果】

【0025】

本発明によると、異方性樹脂成形体を構造解析するにあたり、高精度の結果を短時間で得られるとともに、コンピュータにかかる負荷を大幅に軽減できる。

【図面の簡単な説明】

【0026】

【図1】本発明に係る異方性樹脂成形体の構造解析方法の一例を示すフローチャートである。

【図2】図１における構造解析用モデル作成ステップの一例を示すフローチャートである。

【図3】図２における探索ステップの一例を示すフローチャートである。

【図4】一連の解析プログラムを実現するためのハードウェア資源の一例を示す。

【図5】樹脂流動解析における要素分割後の試験片Ａの形状を示す。

【図6】樹脂流動解析によって得られる試験片Ａにおける異方性フィラーの配向方向分布を示す。

【図7】構造解析を行うために要素分割した後の試験片Ａの形状を示す。

【図8】流動解析によって得られる試験片Ａにおける異方性フィラーの配向方向分布を、構造解析の用に供する有限要素分割モデルに当てはめたときの例を示す。

【図9】試験片Ａを実施例１の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。

【図10】試験片Ｂを比較例１の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。

【図11】樹脂流動解析における要素分割後の試験片Ｂの形状を示す。

【図12】樹脂流動解析によって得られる試験片Ｂにおける異方性フィラーの配向方向分布を示す。

【図13】構造解析を行うために要素分割した後の試験片Ｂの形状を示す。

【図14】ガラス繊維３０重量％含有強化ポリブチレンテレフタレートの温度と比容積との関係を示す。

【図15】流動解析によって得られる試験片Ｂにおける異方性フィラーの配向方向分布を、構造解析の用に供する有限要素分割モデルに当てはめたときの例を示す。

【図16】試験片Ｂにおいて線膨張率を評価した位置を示す。

【図17】試験片Ｂを実施例２の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。

【図18】試験片Ｂを比較例２の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。

【発明を実施するための形態】

【0027】

以下、本発明の具体的な実施形態について、詳細に説明するが、本発明は、以下の実施形態に何ら限定されるものではなく、本発明の目的の範囲内において、適宜変更を加えて実施することができる。

【0028】

＜異方性樹脂成形体の構造解析方法＞
図１は、本発明に係る異方性樹脂成形体の構造解析方法の一例を示すフローチャートである。本発明に係る構造解析方法は、異方性フィラーを含有する異方性樹脂成形体について樹脂流動解析を行い、要素分割された複数の第１領域の各々における配向状態を算出するための配向状態算出基礎情報を含む第１領域配向情報を取得する第１領域配向情報取得ステップ（Ｓ１）と、第１領域配向情報取得ステップ（Ｓ１）とは別に、異方性樹脂成形体の構造解析のために異方性樹脂成形体を複数の第２領域に要素分割する構造解析用要素分割ステップ（Ｓ２）と、異方性樹脂成形体の構造解析を行うための構造解析用モデルを作成する構造解析用モデル作成ステップ（Ｓ３）と、構造解析モデルに基づいて異方性樹脂成形体の構造解析を行う構造解析ステップ（Ｓ４）とを含む。

【0029】

〔第１領域配向情報取得ステップ（Ｓ１）〕
まず、第１領域配向情報取得ステップ（Ｓ１）について説明する。異方性樹脂成形体の配向状態を取得するにあたり、Ｘ線ＣＴ等を用いて、射出成形等の加工によって得られる実際の異方性樹脂成形体を観察することによって取得することも考えられる。しかしながら、試作の前段階である設計段階にある場合、実際の異方性樹脂成形体が試作されていないことから、配向状態を実測するのは難しい。そのため、樹脂流動解析を用い、要素分割された複数の第１領域の各々における配向状態をシミュレートするのが好ましい。

【0030】

樹脂流動解析の手法は特に限定されるものでなく、例えば、有限要素法を用いた樹脂流動解析等が挙げられる。有限要素法を用いた樹脂流動解析により配向状態を取得する場合、異方性樹脂成形体の肉厚方向に対し、５分割以上で要素分割することが好ましい。要素分割の程度が４分割以下である場合、肉厚方向の繊維配向が平均化され、構造解析の精度が低下し得るため、好ましくない。

【0031】

配向状態の種類は特に限定されるものでなく、例えば、配向度、配向方向等が挙げられる。

【0032】

配向度に関する配向度情報は、通常３×３のマトリクスにて表現される。このマトリクスから、３方向の固有ベクトル、固有値を求める。それぞれの固有ベクトルから配向方向が得られる。またそれぞれの固有値から配向度が得られる。

【0033】

必須ではないが、射出成型の際に２点ゲートを用いると、ウエルドが発生し、ウエルド部に応力が集中するため、異方性樹脂成形体がウエルド部を含むものとしてモデル化することが好ましい。

【0034】

また、流動解析の精度、ひいては、構造解析の最終的な精度を高めるため、ゲート、ランナー、インサート部等についてもモデル化することが好ましい。

【0035】

〔構造解析用要素分割ステップ（Ｓ２）〕
続いて、構造解析用要素分割ステップ（Ｓ２）について説明する。構造解析用要素分割ステップ（Ｓ２）は、第１領域配向情報取得ステップ（Ｓ１）とは別に、異方性樹脂成形体の構造解析のために異方性樹脂成形体を複数の第２領域に要素分割するステップである。

【0036】

構造解析用要素分割ステップ（Ｓ２）は、第１領域配向情報取得ステップ（Ｓ１）での要素分割とは別に要素分割を行うものであれば特に限定されるものでなく、第１領域配向情報取得ステップ（Ｓ１）の後に行ってもよいし、第１領域配向情報取得ステップ（Ｓ１）の前に行ってもよいし、第１領域配向情報取得ステップ（Ｓ１）と並行して行ってもよい。

【0037】

要素分割する手法は特に限定されるものでない。一例として、まず、ＣＡＤインターフェース等を利用して、異方性樹脂成形体の形状をパソコン等に取り込むか、あるいは、ＣＡＤシステムにより異方性樹脂成形体の形状を作成し、モデル化範囲を設定する。次いで、要素分割プリプロセッサ等で有限要素法等の要素分割を行うことで異方性樹脂成形体を複数の領域に分割する。

【0038】

要素の形状は特に限定されるものでなく、四面体１次要素、四面体２次要素、六面体１次要素、六面体２次要素等が選択可能であり、有限要素法ソフトウェアの仕様、計算するコンピュータシステムの仕様、計算コスト等に応じて適宜選択すればよい。

【0039】

要素の数も特に限定されるものでなく、計算精度、計算時間等を考慮して適宜選択すればよい。

【0040】

〔構造解析用要素分割ステップ（Ｓ３）〕
構造解析用モデル作成ステップ（Ｓ３）は、異方性樹脂成形体の構造解析を行うための予備段階として構造解析用モデルを作成するステップである。以下、図２を参照しながら、構造解析用モデル作成ステップ（Ｓ３）を詳しく説明する。

【0041】

図２は、図１における構造解析用モデル作成ステップ（Ｓ３）の一例を示すフローチャートであり、本発明に係る構造解析用モデル作成方法に相当する。構造解析用モデル作成ステップ（Ｓ３）は、異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、所定のポアソン比及び所定の線膨張率の少なくとも一方を含む物性情報を作成する物性情報作成ステップ（Ｓ３１）と、異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップ（Ｓ３２）と、異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップ（Ｓ３３）と、探索ステップ（Ｓ３３）で最も近いとされた第１領域における第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップ（Ｓ３４）と、物性情報を参照し、複数の第２領域の各々について、第２領域配向状態に応じた物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップ（Ｓ３５）とを含んで構成される。

【0042】

［物性情報作成ステップ（Ｓ３１）］
物性情報作成ステップ（Ｓ３１）は、配向度等で例示される配向状態に応じた物性情報（機械物性値）を予め作成するステップであり、具体的には、異方性樹脂成形体と配向状態は異なるが、異方性樹脂成形体を構成する樹脂の組成や樹脂成形体の形状等、他の条件は同一である複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成するステップである。物性情報（機械物性値）を予め作成することにより、入力する物性情報の確認、修正が容易となるという利点がある。

【0043】

物性情報の作成情報は特に限定されるものでなく、剛性マトリクスを構成する物性情報を作成する場合、Ｈａｌｐｉｎ−Ｔｓａｉ式を元に、一部を改良したモデルを用いることが挙げられる。また、作成された物性情報は、ハードディスク、メモリ等に例示される、ハードウェアの記録領域に格納される。

【0044】

【数13】

【0045】

式（１）において、ν_２３は、対象材料の主軸方向に直交する２方向に対するポアソン比である。また、ν_ｆは、対象材料に含まれる異方性フィラーのポアソン比である。また、ν_ｍは、対象材料を構成する樹脂組成物のポアソン比である。また、Ｖ_ｆは、対象材料に含まれる異方性フィラーの対象材料に対する体積含有率である。

【0046】

本発明では、ポアソン比を求める方法を改良している。一般的に、各方向のポアソン比を実験的に求める手法として、特開２００９−１２８０３３号公報に記載の手法等が知られているが、配向度依存性を求めることは極めて困難であるため、実際はモデル計算等で求めている。その際、算出方法によっては実際に不可能な値となりやすい。その場合、正定値が負となり、構造解析ソフトウェアによっては計算自体が不可能になる場合がある。式（１）を採用することで、材料物性における正定値が正になりやすくなるため、汎用構造計算ソフトウェアの取扱いが容易になる利点がある。

【0047】

【数14】

式（２）において、α_νは、対象材料の体積膨張率である。また、α_１１は、対象材料の主軸方向の線膨張率である。また、α_２２、α_３３は、材料主軸方向に直交する方向の線膨張率であり、α２２≧α３３という関係性がある。

【0048】

式（２）において、計算に各方向の線膨張率と体積膨張率との関係を用いたが、この関係を満たさないと、体積の膨張・収縮量と各方向の膨張・収縮量が矛盾する。しかしながら、線膨張率を求める方式にて一般的なＳｃｈａｐｅｒｙモデル、Ｈａｓｈｉｎ−Ｓｈｔｒｉｋｍａｎモデル等ではこの関係を必ずしも満たさないため、計算結果の精度が悪い原因になっていた。

【0049】

α_１１は、下記式（３）によって求められることが好ましく、α_２２は、下記式（４）によって求められることが好ましい。

【数15】

【0050】

式（３）において、α_ｓは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率である。また、α_ｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、式（３）’で規定される。λ_１１は、対象材料の主軸方向の配向度である。ｘは、配向度に対する各方向の非線形度を示すパラメーターである。

【0051】

また、式（３）’において、α_νは、対象材料の体積膨張率である。

【0052】

【数16】

【0053】

式（４）において、α_ｂは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率である。α_ｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、式（４）’で規定される。λ_１１は、対象材料の主軸方向の配向度である。ｘは、配向度に対する各方向の非線形度を示すパラメーターである。

【0054】

また、式（４）’において、α_νは、対象材料の体積膨張率である。

【0055】

本発明の実施形態では、仮想異方性樹脂成形体の主軸方向の配向度に応じて連続的に材料物性値を与えるのではなく、区分的に扱う。具体的に、仮想異方性樹脂成形体の主軸方向の配向度λ_１１は、１／３以上１以下の値をとり得るため、例えば、λ_１１が０．９以上１以下の範囲内にある場合を一種類の材料物性を持つものとする。この区分については、区間が細かいほど連続的になり、解析精度がより高まるが、区分数が多くなるため、材料物性の確認／変更は困難になる。構造解析での解析精度、及び材料物性の確認／変更の容易性の双方を考慮すると、複数種類の仮想成形体は、前記配向度の範囲に応じて１０種類以上１０００種類以下の範囲で設定されることが好ましい。

【0056】

物性情報として、上記式（１）で表されるポアソン比、及び式（２）で表される線膨張率のほか、下記式（５）で規定される弾性率、下記式（６）で規定されるせん断弾性率等が挙げられる。

【0057】

【数17】

【0058】

【数18】

【0059】

式（５）において、Ｅ_１１は、対象材料の主軸方向の弾性率であり、Ｅ_ｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率であり、ζ_１１は、式（５）’で示される値であり、ηは、式（５）’’で示される値であり、Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。

【0060】

【数19】

【0061】

式（５）’において、λ_１１は、対象材料の主軸方向の配向度であり、Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。

【0062】

【数20】

式（５）’’において、Ｅ_ｆは、対象材料に含まれるフィラーの弾性率であり、Ｅ_ｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。

【0063】

式（６）において、Ｇ_１２は、対象材料の主軸方向のせん断弾性率であり、Ｇ_ｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率であり、ζ_１１は、式（６）’で示される値であり、η_ｇは、式（６）’’で示される値であり、Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。

【0064】

【数21】

【0065】

式（６）’において、λ_１１は、対象材料の主軸方向の配向度であり、Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。

【0066】

【数22】

式（６）’’において、Ｇ_ｆは、対象材料に含まれるフィラーのせん断弾性率であり、Ｇ_ｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。

【0067】

式（５）’及び式（６）’に関し、配向度等で例示される配向状態を元に弾性率等を求める手法は、これまでも均質化法として知られている。しかしながら、均質化法は、弾性率等を求める際の計算方法が複雑であるため、計算エラー、プログラムコーディングエラー等を生じ得る。計算手順、計算プログラム作成を容易にするため、微小領域における配向度を考慮できるモデルとして、上記の式（５）’及び式（６）’を採用することが好ましい。

【0068】

式（５）’及び式（６）’において、アスペクト比Ｌ／ｄは、実測値を採用することが好ましいが、異方性樹脂成型体に含まれるフィラーのアスペクト比は、フィラーの破損等により分布があるとともに、射出成型前の状態からは変化する。したがって、式（５）’及び式（６）’の運用に関しては、アスペクト比をパラメーターとして扱ってもよい。

【0069】

その他、物性情報として、以下の式（７）から式（９）を例示できる。

【0070】

【数23】

【0071】

式（７）において、ν_１２及びν_１３は、対象材料の主軸方向に直交する方向に対するポアソン比である。また、Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率であり、ν_ｆは、対象材料に含まれる異方性フィラーのポアソン比である。また、Ｖ_ｍは、樹脂組成物の対象材料に対する体積含有率であり、ν_ｍは、対象材料を構成する樹脂組成物のポアソン比である。

【0072】

【数24】

【0073】

式（８−１）及び式（８−２）において、Ｅ_２２及びＥ_３３は、対象材料の主軸方向に直交する方向の弾性率であり、Ｅ_ｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。また、ζ_２２は、２であり、ηは、上記式（５）’’で示される値であり、Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。

【0074】

【数25】

【0075】

式（９−１）及び式（９−２）において、Ｇ_２３は、対象材料の主軸方向に直交する方向のせん断弾性率であり、Ｇ_ｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。また、ζ_２２は、２であり、η_ｇは、式（６）’’で示される値であり、Ｖ_ｆは、異方性フィラーの対象材料に対する体積含有率である。

【0076】

［第１領域配向状態算出ステップ（Ｓ３２）］
第１領域配向状態算出ステップ（Ｓ３２）は、異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出するステップである。本実施形態では、配向状態が配向方向及び配向度であるものとして説明するが、これに限るものではない。

【0077】

配向度に関する配向度情報は、通常３×３のマトリクスにて表現される。このマトリクスから、固有ベクトル、固有値を求める。固有ベクトルから配向方向が得られる。また固有値から配向度が得られる。固有ベクトルの３方向を座標系に変換し、後のステップにて各要素に与える。また、配向度に応じて物性情報作成ステップ（Ｓ３１）にて得られた物性情報を後のステップにて各要素に設定する。これらの設定は、構造解析の前後において確認しやすいため、物性設定における不具合を確認しやすくなるという効果がある。

【0078】

［探索ステップ（Ｓ３３）］
探索ステップ（Ｓ３３）は、異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索するステップである。探索ステップ（Ｓ３３）の態様は特に限定されるものでないが、例えば、図３に示すように、複数の第１領域の各々について第１重心位置を導出する第１重心位置導出ステップ（Ｓ３３１）と、複数の第２領域の各々について第２重心位置を導出する第２重心位置導出ステップ（Ｓ３３２）と、第２領域の各々について、第２重心位置から最も近い第１重心位置を探索し、第２重心位置から最も近い第１重心位置を有する第１領域を第２領域から最も近い第１領域とする最短第１領域設定ステップ（Ｓ３３３）とを有する態様を例示できる。

【0079】

（第１重心位置導出ステップ（Ｓ３３１））
図３を参照しながら、探索ステップ（Ｓ３３）の具体例について説明する。第１重心位置導出ステップ（Ｓ３３１）では、第１領域配向情報取得ステップ（Ｓ１）にて計算に用いた第１領域に係る要素毎に、要素を構成する節点の位置から第１重心位置を求める。第１重心位置は、次に計算する構造解析用モデルの第２重心位置と比較される。

【0080】

（第２重心位置導出ステップ（Ｓ３３２））
第２重心位置導出ステップ（Ｓ３３２）では、第１重心位置導出ステップ（Ｓ３３１）と同様に、構造解析用要素分割ステップ（Ｓ２）で分割された第２領域の要素毎に、要素を構成する節点の位置から重心位置を求める。

【0081】

（最短第１領域設定ステップ（Ｓ３３３））
最短第１領域設定ステップ（Ｓ３３３）では、構造解析用要素分割ステップ（Ｓ２）で分割された第２領域の要素毎に、第２重心位置導出ステップ（Ｓ３３２）で求めた第２重心位置と、第１重心位置導出ステップ（Ｓ３３１）で求めた第１重心位置との間の距離を計算し、最も近い位置の情報を当てはめる。この最短第１領域設定ステップ（Ｓ３３３）を終えると、図２の第２領域配向状態設定ステップ（Ｓ３４）に移る。

【0082】

［第２領域配向状態設定ステップ（Ｓ３４）］
図２に戻る。第２領域配向状態設定ステップ（Ｓ３４）は、探索ステップ（Ｓ３３）で最も近いとされた第１領域における第１領域配向状態を各々の第２領域配向状態とするステップである。この方式の場合、第１領域配向情報取得ステップ（Ｓ１）にて用いた樹脂流動解析での要素分割が粗く、構造解析用要素分割ステップ（Ｓ２）で分割された第２領域の要素の第２重心位置に第１領域配向情報取得ステップ（Ｓ１）にて用いた第１領域に係る要素がない場合でも、第１領域配向状態の情報を第２領域に当てはめることができる。

【0083】

この方式に従えば、要素分割状態によらず第１領域配向状態の情報を第２領域に当てはめることが可能である。この当てはめは、要素の形状関数を用いて配向状態の情報を補間して求めることも可能であるが、要素分割が異なる場合には補間できない状況になる場合があるため、配向状態の情報を求める際の情報密度、構造解析の要素分割に制約が生じる場合があり、その制約により、構造解析自体の精度に影響を及ぼす場合がある。

【0084】

［第２領域物性情報設定ステップ（Ｓ３５）］
第２領域物性情報設定ステップ（Ｓ３５）は、物性情報作成ステップ（Ｓ３１）で予め作成され、ハードウェアの記録領域に格納されている物性情報を参照し、複数の第２領域の各々について、第２領域配向状態に応じた物性情報を各々の第２領域物性情報とするステップである。このステップでは、第２領域配向状態設定ステップ（Ｓ３４）にて求めた第２領域配向状態の情報に基づいた座標系、及び配向度に応じた材料物性を構造解析の要素として設定する。この第２領域物性情報設定ステップ（Ｓ３５）を終えると、構造解析の実行が可能になるため、図１の構造解析ステップ（Ｓ４）に移る。

【0085】

〔構造解析ステップ（Ｓ４）〕
構造解析ステップ（Ｓ４）は、構造解析モデルに基づいて異方性樹脂成形体の構造解析を行うステップである。このステップでは、第２領域配向状態設定ステップ（Ｓ３４）にて設定した第２領域配向状態の情報と、第２領域物性情報設定ステップ（Ｓ３５）にて設定した第２領域物性情報とを含む構造解析モデル情報を用い、構造解析ソフトウェアのプログラムを実行して構造解析を実行する。構造解析における拘束条件、荷重条件、計算条件等の各種条件は、構造解析ソフトェアの仕様にしたがえばよく、特に限定されない。

【0086】

構造解析により、異方性樹脂成形体の変形、応力等をシミュレートできる。このシミュレートの結果から、異方性樹脂成形体における破損の有無等を判定できる。

【0087】

具体的態様として、異方性樹脂射出成形体の反り変形は、構造解析における荷重条件を温度荷重及び残留応力荷重に設定することで解析可能である。そのため、本発明に係る構造解析方法を使用することで、異方性樹脂射出成形体の反り変形を予測できる。

【0088】

＜解析プログラム＞
本発明に係る構造解析用モデル作成プログラムを含む一連の解析プログラムは、ソフトウェアとハードウェア資源とが協働することによって実現される。

【0089】

図４は、一連の解析プログラムを実現するためのハードウェア資源Ｈの一例を示す。ハードウェア資源Ｈは、情報処理装置１と、設計者からの各種要求を受け付ける入力装置２と、情報処理装置１が行った解析結果を出力する出力装置３とを備える。また、情報処理装置１は、ＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）等のネットワークＮＷを介して、ＣＡＤ装置４に接続されている。

【0090】

情報処理装置１は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１０と、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）等により構成される主記憶装置２０と、Ｉ／Ｏインタフェース３０と、ハードディスク等により構成される補助記憶装置４０と、ネットワークＮＷに接続されている装置との間で行うデータ授受の制御を行うネットワークインタフェース５０とを備える。

【0091】

補助記憶装置４０には、上述した一連のステップをハードウェア装置１に実行させるための解析プログラム４１が格納されている。解析プログラム４１は、流動解析による第１領域配向情報取得ステップ（Ｓ１）をハードウェア装置１に実行させるための第１領域配向情報取得プログラム４１Ａと、構造解析用要素分割ステップ（Ｓ２）をハードウェア装置１に実行させるための構造解析用要素分割プログラム４１Ｂと、構造解析用モデル作成ステップ（Ｓ３）をハードウェア装置１に実行させるための構造解析用モデル作成プログラム４１Ｃと、構造解析ステップ（Ｓ４）をハードウェア装置１に実行させるための構造解析プログラム４１Ｄとを含んで構成される。本発明に係る構造解析方法は、ＣＰＵ１０が補助記憶装置４０に格納されている解析プログラム４１を主記憶装置２０にロードして実行することにより実現される。特に、本発明に係る異方性樹脂成形体の構造解析用モデル作成方法は、ＣＰＵ１０が補助記憶装置４０に格納されている構造解析用モデル作成プログラム４１Ｃを主記憶装置２０にロードして実行することにより実現される。

【0092】

上記の実施形態では、第１領域配向情報取得ステップ（Ｓ１）から構造解析ステップ（Ｓ４）に至るまでの各ステップは、複数のプログラムが組み合わせられて実行されることで行われているが、これに限られることなく、最初から一体として構築されたプログラムでもよく、また、実行されるコンピュータの形態や規模、設置場所等も限定されるものではない。

【実施例】

【0093】

以下、実施例及び比較例を示し、本発明を具体的に説明するが、本発明はこれらの実施例に限定されるものではない。

【0094】

＜試験例１＞
試験例１では、平板部分の長片側にリブを設置した形状であって、射出成型にて得られる試験片Ａを用いて、３点曲げ試験時における耐荷重を検討した。

【0095】

〔実施例１〕
まず、この試験片Ａに対し、ＪＩＳ１７８に準拠した曲げ試験を行った。固定冶具間の間隔を６２ｍｍに設定し、試験片中央部に圧子を配置し、圧子をひずみ速度が１％／ｍｉｎにて移動させ、試験機器に搭載のロードセンサーにて荷重を測定した。その結果、試験片は荷重２７２Ｎにて破損した。

【0096】

［第１領域配向情報取得ステップ（Ｓ１）］
続いて、試験片Ａの形状をＣＡＤデータとして作成し、ランナー、ゲート等を加え有限要素分割を実行した。樹脂流動解析には、ＡＵＴＯＤＥＳＫ社ＡＵＴＯＤＥＳＫ ＳＩＭＵＬＡＴＩＯＮ ＭＯＬＤＦＬＯＷ ＩＮＳＩＧＨＴを用いた。樹脂流動解析の条件を以下に示す。
樹脂：ガラス繊維３０重量％含有強化ポリブチレンテレフタレート
樹脂温度：２６０℃
金型温度：８０℃
射出流量：５７．７ｃｍ^３／ｓ
保圧圧力：７０ＭＰａ
保圧時間：１５秒
冷却時間：１０秒
要素：四面体１次要素（分割数 ４５０５３８）

【0097】

樹脂流動解析における要素分割後の試験片Ａの形状を図５に示し、樹脂流動解析によって得られる試験片Ａにおける異方性フィラーの配向方向分布を図６に示す。

【0098】

［構造解析用要素分割ステップ（Ｓ２）］
流動解析用の有限要素モデルとは別に、異方性樹脂成形体の構造解析のために、試験片Ａを複数の第２領域に要素分割した。要素分割後の試験片Ａの形状を図７に示す。なお、図７においては、ランナー、ゲート等を削除し、代わりに試験片支持台、および圧子を加えた。要素については四面体２次要素を用い、分割数は６６５６５であった。試験片支持台を完全拘束し、圧子に、試験片の耐荷重である荷重２７２Ｎを与え、最大主応力および発生位置を検討した。

【0099】

［構造解析用モデル作成ステップ（Ｓ３）］
まず、物性情報作成ステップ（Ｓ３１）において、試験片Ａと配向状態が異なる複数種類の仮想成形体の各々について、上記式（１）〜（９）にしたがって物性情報を作成し、コンピュータの所定の記憶領域に格納した。式（１）〜（９）にしたがって得られる物性情報は、以下のとおりであった。
対象材料に含まれるフィラーの弾性率Ｅ_ｆ ：７２０００ＭＰａ
対象材料を構成する樹脂組成物からなる成形体の弾性率Ｅ_ｍ：２５００ＭＰａ
異方性フィラーの対象材料に対する体積含有率Ｖ_ｆ ：１８ｖｏｌ％
樹脂組成物の対象材料に対する体積含有率Ｖ_ｍ ：８２ｖｏｌ％
対象材料に含まれるフィラーのせん断弾性率Ｇ_ｆ ：２９５００ＭＰａ
対象材料を構成する樹脂組成物からなる成形体のせん断弾性率Ｇ_ｍ：９２６ＭＰａ
対象材料に含まれる異方性フィラーのポアソン比ν_ｆ ：０．２２
対象材料を構成する樹脂組成物のポアソン比ν_ｍ ：０．３５
対象材料に含まれる異方性フィラーのアスペクト比Ｌ／ｄ ：１４

【0100】

そして、試験片Ａの樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップ（Ｓ３２）と、試験片Ａの構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップ（Ｓ３３）と、探索ステップ（Ｓ３３）で最も近いとされた第１領域における第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップ（Ｓ３４）と、物性情報を参照し、複数の第２領域の各々について、第２領域配向状態に応じた物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップ（Ｓ３５）とを経て構造解析モデルを完成させる。完成した構造解析モデルを図８に示す。元となる図６の分布状況と比較すると、方向を示す位置が異なるが、両者の方向が一致している。

【0101】

［構造解析ステップ（Ｓ４）］
得られた構造解析用のモデルを用いて計算する。計算には、構造解析ソフトＮＸ Ｉ−ＤＥＡＳ（シーメンスＰＬＭソフトウェア社製）を用いた。

【0102】

図９は、試験片Ａを実施例１の手法にて構造解析することによって得られる主応力分布を示す。最大応力を示す位置は、試験片の中央部から若干ずれている。試験片の成形の際には２点ゲートを用いたので、ウエルドが発生し、ウエルド部に応力が集中した結果と考えられる。また、実施例１における最大主応力は１３９ＭＰａであった（表１）。

【0103】

また、試験片の引張破壊強度の実測値は１４０ＭＰａであり、ウエルド部が破壊の起点となっていた（表１）。

【0104】

このように、実施例１の結果は、実測に極めて近く、実施例１において本発明の精度は極めて良好といえる。

【0105】

〔比較例１〕
実施例１と同一の試験片について、弾性率として、流動方向の弾性率と流動垂直方向の弾性率とを平均化した値である６３２０ＭＰａを用い、ポアソン比として、既存のモデル計算によって求められる０．３９を用いたこと以外は、実施例１と同じ手法にて、試験片Ａの構造解析を行った。

【0106】

図１０は、試験片Ａを比較例１の手法にて構造解析することによって得られる主応力分布を示す。最大応力を示す位置は試験片中央部であり、ウエルド発生位置とは異なる位置であった。比較例１における最大主応力は１０６ＭＰａであった（表１）。

【0107】

このように、比較例１の結果は、実測と異なり、比較例１の精度は、実施例１に比べて劣るといえる。

【0108】

【表1】

【0109】

＜試験例２＞
試験例２では、平板部分に箱型形状、リブ、円筒等を設置した形状であって、射出成型にて得られる試験片Ｂを用いて、熱膨張挙動を検討した。

【0110】

〔実施例２〕
まず、射出成型にて得られた試験片Ｂを切り出し、１２０℃環境にて恒温槽にてアニールを２時間程実施し、その後放冷した。その後、線膨張率測定装置に設置し、−３０℃から室温まで温度を上昇させ、４か所において寸法変化量から、線膨張率を求めた。その結果は後述の表２にまとめた。

【0111】

［繊維配向情報取得ステップ（Ｓ１）］
続いて、実施例１と同様の手法にて、試験片Ｂの有限要素分割を実行した。樹脂流動解析の条件を以下に示す。
樹脂：ガラス繊維３０重量％含有強化ポリブチレンテレフタレート
樹脂温度：２５０℃
金型温度：６０℃
射出流量：３２ｃｍ^３／ｓ
保圧圧力：７８．４ＭＰａ
保圧時間：１０秒
冷却時間：１５秒
要素：四面体１次要素（分割数 ４８６７４９）

【0112】

樹脂流動解析における要素分割後の試験片Ｂの形状を図１１に示し、樹脂流動解析によって得られる試験片Ｂにおける異方性フィラーの配向方向分布を図１２に示す。

【0113】

［構造解析用要素分割ステップ（Ｓ２）］
流動解析用の有限要素モデルとは別に、異方性樹脂成形体の構造解析のために、試験片Ｂを複数の第２領域に要素分割した。要素分割後の試験片Ｂの形状を図１３に示す。なお、図１３においては、切り出し後の形状を用いた。要素については四面体２次要素を用い、分割数は６６５６５であった。モデルのコーナー部３ヶ所にＸＹＺ方向、ＹＺ方向、Ｚ方向を拘束し、温度荷重１００℃を与えた。

【0114】

［材料物性の設定ステップ（Ｓ３）］
まず、物性情報作成ステップ（Ｓ３１）において、試験片Ａと配向状態が異なる複数種類の仮想成形体の各々について、上記式（１）〜（９）にしたがって物性情報を作成し、コンピュータの所定の記憶領域に格納した。式（１）〜（９）にしたがって得られる物性情報は、以下のとおりであった。
完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率α_ｓ ：２．４×１０^−５／℃
完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率α_ｂ ：６．２×１０^−５／℃
対象材料の体積膨張率α_ν ：１．５３×１０^−４／℃
配向度に対する各方向の非線形度ｘ ：１５．５（無次元）

【0115】

体積膨張率α_νは、ＰＶＴ測定装置を用い、圧力が０ＭＰａであるときのガラス繊維３０重量％含有強化ポリブチレンテレフタレート（強化ＰＢＴ）の体積変化を求め、強化ＰＢＴの温度に対する比容積の変化率を計算することによって求めた値である。

【0116】

そして、実施例１と同様に、第１領域配向状態算出ステップ（Ｓ３２）と、探索ステップ（Ｓ３３）と、第２領域配向状態設定ステップ（Ｓ３４）と、第２領域物性情報設定ステップ（Ｓ３５）とを経て構造解析モデルを完成させる。完成した構造解析モデルを図１５に示す。元となる図１２の分布状況と比較すると、方向を示す位置が異なるが、両者の方向が一致している。

【0117】

［構造解析ステップ（Ｓ４）］
続いて、実施例１と同様の手法によって試験片Ｂを構造解析した。

【0118】

図１６に示す位置Ａ，Ｂ，Ｃ，Ｄにおける寸法変化から、各位置における線膨張率を求めた。結果を表２に示す。また、図１７は、試験片Ｂを実施例２の手法にて構造解析することによって得られる変形量分布を示す。図１７から、試験片Ｂは、平面部がゆがむように変形することを示唆している。

【0119】

また、試験片Ｂにおいて、位置Ａ，Ｂ，Ｃ，Ｄにおける線膨張率を実測した。結果を表２に示す。この実測の際、試験片Ｂの平面部がゆがむように変形することを目視にて確認した。

【0120】

このように、実施例２における位置Ａ，Ｂ，Ｃ，Ｄ間の線膨張率差の変化は、実測値を反映しているといえ、実際の試験片Ｂのゆがみを反映しているといえる。したがって、実施例２の結果は、実測に極めて近く、実施例２において本発明の精度は極めて良好といえる。

【0121】

〔比較例２〕
実施例２と同一の試験片Ｂについて、線膨張率として、流動方向の線膨張率と流動垂直方向の線膨張率とを平均化した値である４．３×１０^−５／℃を用いたこと以外は、実施例２と同じ手法にて、試験片Ｂの構造解析を行った。

【0122】

図１６に示す位置Ａ，Ｂ，Ｃ，Ｄにおける寸法変化から、各位置における線膨張率を求めた。結果を表２に示す。比較例２の場合、位置Ａ，Ｂ，Ｃ，Ｄ間で線膨張率の有意差は認められない。また、図１８は、試験片Ｂを比較例２の手法にて構造解析することによって得られる変形量分布を示す。図１８から、試験片Ｂの平面部が平坦のままであり、実際の試験片Ｂのゆがみを反映していないことが明らかである。

【0123】

このように、比較例２の結果は、実測と異なり、比較例２の構造解析精度は、実施例２に比べて劣るといえる。

【0124】

【表2】

【符号の説明】

【0125】

Ｈ ハードウェア資源
１ 情報処理装置
２ 入力装置
３ 出力装置
４ ＣＡＤ装置
１０ ＣＰＵ
２０ 主記憶装置
３０ Ｉ／Ｏインタフェース
４０ 補助記憶装置
４１ 解析プログラム
４１Ａ 第１領域配向情報取得プログラム
４１Ｂ 構造解析用要素分割プログラム
４１Ｃ 構造解析用モデル作成プログラム
４１Ｄ 構造解析プログラム
５０ ネットワークインタフェース
ＮＷ ネットワーク
Ｓ１ 第１領域配向情報取得ステップ
Ｓ２ 構造解析用要素分割ステップ
Ｓ３ 構造解析用モデル作成ステップ
Ｓ４ 構造解析ステップ
Ｓ３１ 物性情報作成ステップ
Ｓ３２ 第１領域配向状態算出ステップ
Ｓ３３ 探索ステップ
Ｓ３４ 第２領域配向状態設定ステップ
Ｓ３５ 第２領域物性情報設定ステップ
Ｓ３３１ 第１重心位置導出ステップ
Ｓ３３２ 第２重心位置導出ステップ
Ｓ３３３ 最短第１領域設定ステップ

【図1】

【図2】

【図3】

【図4】

【図14】

【図16】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図15】

【図17】

【図18】