(57)【要約】【課題】高精度、処理時間短縮、負荷軽減のいずれにも資する異方性樹脂成形体の構造解析用モデルを提供する。【解決手段】本発明は、特定のポアソン比及び／又は特定の線膨張率を含む物性情報を作成する物性情報作成ステップ（Ｓ３１）と、異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の第１領域配向状態を算出する第１領域配向状態算出ステップ（Ｓ３２）と、異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップ（Ｓ３３）と、このステップで最も近いとされた第１領域における第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップ（Ｓ３４）と、上記の物性情報を参照し、複数の第２領域の各々について、第２領域配向状態に応じた物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップ（Ｓ３５）とを含む。【選択図】図２
2016203584.tif 000002
【特許請求の範囲】
【請求項１】
異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルを作成する構造解析用モデル作成方法であって、
前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、
前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、
前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、
前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、
前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとを含む、異方性樹脂成形体の構造解析用モデル作成方法。
【数１】
2016203584.tif 000030
（式（１）において、
ν２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
νｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
νｍは、対象材料を構成する樹脂組成物のポアソン比であり、
Ｖｆは、対象材料に含まれる異方性フィラーの対象材料に対する体積含有率である。）
【数２】
2016203584.tif 000031
（式（２）において、
ανは、対象材料の体積膨張率であり、
α１１は、対象材料の主軸方向の線膨張率であり、
α２２は、材料主軸方向に直交する方向の線膨張率であり、
α３３は、材料主軸方向に直交する方向の線膨張率であり、
α２２≧α３３という関係を有する。）
【請求項２】
前記α１１は、下記式（３）によって求められ、前記α２２は、下記式（４）によって求められる、請求項１に記載の構造解析用モデル作成方法。
【数３】
2016203584.tif 000032
（式（３）において、
αｓは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率であり、
αｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、上記式（３）’で規定され、
λ１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（３）’において、ανは、対象材料の体積膨張率である。）
【数４】
2016203584.tif 000033
（式（４）において、
αｂは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率であり、
αｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、下記式（４）’で規定され、
λ１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（４）’において、ανは、対象材料の体積膨張率である。）
【請求項３】
前記物性情報は、下記式（５）で規定される弾性率をさらに含む、請求項１又は２に記載の構造解析用モデル作成方法。
【数５】
2016203584.tif 000034
（式（５）において、
Ｅ１１は、対象材料の主軸方向の弾性率であり、
Ｅｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率であり、
ζ１１は、式（５）’で示される値であり、
ηは、式（５）’’で示される値であり、
Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。）
【数６】
2016203584.tif 000035
（式（５）’において、
λ１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）
【数７】
2016203584.tif 000036
（式（５）’’において、Ｅｆは、対象材料に含まれるフィラーの弾性率であり、Ｅｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。）
【請求項４】
前記物性情報は、式（６）で規定されるせん断弾性率をさらに含む、請求項１から３のいずれかに記載の構造解析用モデル作成方法。
【数８】
2016203584.tif 000037
（式（６）において、
Ｇ１２は、対象材料の主軸方向のせん断弾性率であり、
Ｇｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率であり、
ζ１１は、式（６）’で示される値であり、
ηｇは、式（６）’’で示される値であり、
Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。）
【数９】
2016203584.tif 000038
（式（６）’において、
λ１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）
【数１０】
2016203584.tif 000039
（式（６）’’において、Ｇｆは、対象材料に含まれるフィラーのせん断弾性率であり、Ｇｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。）
【請求項５】
前記ανは、ＰＶＴ樹脂特性解析によって得られる対象材料の体積膨張率の実測値である、請求項１から４のいずれかに記載の構造解析用モデル作成方法。
【請求項６】
前記異方性樹脂成形体がウエルド部を含む、請求項１から５のいずれかに記載の構造解析用モデル作成方法。
【請求項７】
前記探索ステップは、
前記複数の第１領域の各々について第１重心位置を導出する第１重心位置導出ステップと、
前記複数の第２領域の各々について第２重心位置を導出する第２重心位置導出ステップと、
前記第２領域の各々について、前記第２重心位置から最も近い第１重心位置を探索し、前記第２重心位置から最も近い第１重心位置を有する第１領域を前記第２領域から最も近い第１領域とする最短第１領域設定ステップとを有する、請求項１から６のいずれかに記載の構造解析用モデル作成方法。
【請求項８】
前記配向状態は、対象材料の主軸方向の配向度であり、
前記複数種類の仮想成形体は、前記配向度の範囲に応じて１０種類以上１０００種類以下の範囲で設定される、請求項１から７のいずれかに記載の構造解析用モデル作成方法。
【請求項９】
請求項１から８のいずれかに記載の構造解析用モデル作成方法を使用して異方性樹脂成形体を構造解析する構造解析方法であって、
前記異方性樹脂成形体について樹脂流動解析を行い、要素分割された前記複数の第１領域の各々における配向状態を算出するための配向状態算出基礎情報を含む第１領域配向情報を取得する第１領域配向情報取得ステップと、
前記第１領域配向情報取得ステップとは別に、前記異方性樹脂成形体の構造解析のために前記異方性樹脂成形体を前記複数の第２領域に要素分割する構造解析用要素分割ステップと、
前記第２領域配向状態の情報及び前記第２領域物性情報を含む構造解析モデル情報に基づいて前記異方性樹脂成形体の構造解析を行う構造解析ステップとを含み、
前記第１領域配向状態算出ステップは、前記配向状態算出基礎情報に基づいて前記複数の第１領域の各々における配向状態に相当する第１領域配向状態を算出するステップであり、
前記探索ステップは、前記構造解析用モデルにおける複数の第２領域の各々について、第２領域から最も近い第１領域を探索するステップである、異方性樹脂成形体の構造解析方法。
【請求項１０】
請求項９に記載の構造解析方法を使用することで、前記異方性樹脂射出成形体の反り変形を予測する反り変形予測方法。
【請求項１１】
異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルの作成をコンピュータに実行させるための構造解析用モデル作成プログラムであって、
前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、
前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、
前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、
前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、
前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとをコンピュータに実行させるためのプログラム。
【数１１】
2016203584.tif 000040
（式（１）において、
ν２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
νｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
νｍは、対象材料を構成する樹脂組成物のポアソン比である。）
【数１２】
2016203584.tif 000041
（式（２）において、
ανは、対象材料の体積膨張率であり、
α１１は、対象材料の主軸方向の線膨張率であり、
α２２は、材料主軸方向に直交する方向の線膨張率であり、
α３３は、材料主軸方向に直交する方向の線膨張率であり、
α２２≧α３３  という関係を有する。）
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、異方性樹脂成形体の構造解析用モデル作成方法、構造解析方法及び構造解析用モデル作成プログラムに関する。
【背景技術】
【０００２】
熱可塑性樹脂と異方性フィラーとを含有する異方性樹脂組成物を射出成形等することによって得られる異方性樹脂成形体は、工業製品として広く使用されており、数値解析による形状検討が多用されている。ところで、数値解析による樹脂成形体の強度設計には、弾性率、ポアソン比、線膨張率等の物性値が必要である。樹脂成形体が異方性フィラーを含有する場合、これらの物性値を実測すること、及び予測することのいずれも容易ではない。
【０００３】
特許文献１に記載の弾性率計算方法は、対象物の繊維配向方向の第１弾性率と、繊維配向と直交する方向の第２弾性率と、構造解析を行う際の荷重方向とを予め格納しておき、対象物の任意の要素ついて繊維配向の向きと荷重方向とが平行している場合には第１弾性率を、繊維配向の向きと荷重方向とが直交している場合には第２弾性率をそれぞれ設定し、さらに、繊維配向の向きと荷重方向とが平行も直交もしていない場合には荷重方向と繊維配向の向きとの角度と第１及び第２弾性率のいずれかに基づいて弾性率を算出して弾性率を設定することを特徴とする。
【０００４】
特許文献２には、各繊維単体の物性と含有繊維をマトリクスと等価な等価介在物として取り扱う為の仮想歪テンソルを加算演算することにより、複合材全体の物性を表現するテンソルを求めるという手法が紹介されている。
【０００５】
特許文献３では、含有する強化繊維及びベースレジン単体の機械物性を使用して評価形状における流動配向解析を実施し、算出した樹脂成形品の配向パラメータより算出した機械物性を、実際に評価形状をした成形品の機械物性を測定した数値と比較し、計算値が実測定値と等しくなるようなベースレジン、および強化繊維単体の機械物性を同定し、同定した機械物性を使用して成形部品の変形量を高精度に評価する方法が記載されている。
【０００６】
しかしながら、いずれの方法も、繊維配向による成形品の機械物性予測方法による予測値と実際の成形品の物性値には数値的な不一致が存在しており、さらなる改善が求められていた。
【０００７】
この不一致の原因としては、計算にて求める場合の仮定と実際との違いが挙げられる。例えば、射出充填される際に異方性フィラーが破損して、計算に使用している異方性フィラーのアスペクト比の値と実際の値が異なってしまうこと、繊維含有率の不均一、異方性フィラー同士の干渉等による異方性フィラー含有時の流動配向の予測力不足が考えられる。すなわち、計算上の仮定が実現象を十分表現できないことに加え、各構成要素の物性値が不正確であることで、解析モデル構築時の誤差に入力物性値の誤差が加算されて複合された誤差が存在しており、フィラー配向計算結果から機械物性を予測しても、この予測では不十分であった。
【０００８】
予測の精度を高めるため、特許文献４では、充填剤を所定の割合で含む樹脂組成物を成形してなる樹脂成形品についてＸ線ＣＴ測定等により、複数のスライス画像を取得し、スライス画像を微小面積の画素に分割し、その情報を元に樹脂成形品内の充填剤の配向状態を反映させた有限要素法モデルを作成し有限要素法解析を介して機械物性を同定し、成形部品の強度、変形量を高精度に評価する方法が記載されている。
【０００９】
また、特許文献５では、フィラーを含む射出成形品内のフィラーの３次元配向のデータをＸ線ＣＴによる３次元画像から取得し、せん断応力の積分値および分子配向状態のデータに基づいて均質化法を介して線膨張係数を求め、構造解析により、射出成形品の熱間反りを求めている。
【００１０】
しかしながら、特許文献４及び５に記載の方法では、各構成要素の物性値を得るために多大な実験、測定が必要である上、それを求める際の計算時間も多大となる。それにより、物性値を決定する際の計算に要するコストが膨大となっている。
【先行技術文献】
【特許文献】
【００１１】
【特許文献１】特開平９−２３７２６７号公報
【特許文献２】特開平７−３０４０５６号公報
【特許文献３】特開２００４−２５７９６号公報
【特許文献４】特開２０１１−７５８号公報
【特許文献５】特開２０１４−１００８７９号公報
【発明の概要】
【発明が解決しようとする課題】
【００１２】
本発明は、以上のような課題を解決するためになされたものであり、構造解析が容易ではないと考えられている異方性樹脂成形体を構造解析するにあたり、構造解析の精度を高めるとともに、解析結果を短時間で、かつ、コンピュータにかかる負荷を大幅に軽減可能な構造解析用モデルを提供することを目的とする。
【課題を解決するための手段】
【００１３】
本発明者らは、上記課題を解決するために鋭意研究を重ねた。その結果、異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体を想定し、これら複数種類の仮想成形体の各々について、ポアソン比と線膨張率の少なくとも一方を含む物性情報を予め作成し、この物性情報を構造解析のために要素分割された複数の第２領域の各々に当てはめることで、上記の課題を解決できることを見出し、本発明を完成するに至った。より具体的に、本発明は、以下のものを提供する。
【００１４】
［１］本発明は、異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルを作成する構造解析用モデル作成方法であって、前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとを含む、異方性樹脂成形体の構造解析用モデル作成方法である。
【数１】
2016203584.tif 000003
（１）
（式（１）において、
ν２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
νｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
νｍは、対象材料を構成する樹脂組成物のポアソン比であり、
Ｖｆは、対象材料に含まれる異方性フィラーの対象材料に対する体積含有率である。）
【数２】
2016203584.tif 000004
（式（２）において、
ανは、対象材料の体積膨張率であり、
α１１は、対象材料の主軸方向の線膨張率であり、
α２２は、材料主軸方向に直交する方向の線膨張率であり、
α３３は、材料主軸方向に直交する方向の線膨張率であり、
α２２≧α３３という関係を有する。）
【００１５】
［２］また、本発明は、前記α１１は、下記式（３）によって求められ、前記α２２は、下記式（４）によって求められる、［１］に記載の構造解析用モデル作成方法である。
【数３】
2016203584.tif 000005
（式（３）において、
αｓは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率であり、
αｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、上記式（３）’で規定され、
λ１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（３）’において、ανは、対象材料の体積膨張率である。）
【数４】
2016203584.tif 000006
（式（４）において、
αｂは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率であり、
αｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、下記式（４）’で規定され、
λ１１は、対象材料の主軸方向の配向度であり、
ｘは、配向度に対する各方向の非線形度を示すパラメーターである。）
（式（４）’において、ανは、対象材料の体積膨張率である。）
【００１６】
［３］また、本発明は、前記物性情報は、下記式（５）で規定される弾性率をさらに含む、［１］又は［２］に記載の構造解析用モデル作成方法である。
【数５】
2016203584.tif 000007
（式（５）において、
Ｅ１１は、対象材料の主軸方向の弾性率であり、
Ｅｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率であり、
ζ１１は、式（５）’で示される値であり、
ηは、式（５）’’で示される値であり、
Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。）
【数６】
2016203584.tif 000008
（式（５）’において、
λ１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）
【数７】
2016203584.tif 000009
（式（５）’’において、Ｅｆは、対象材料に含まれるフィラーの弾性率であり、Ｅｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。）
【００１７】
［４］また、本発明は、前記物性情報は、式（６）で規定されるせん断弾性率をさらに含む、［１］から［３］のいずれかに記載の構造解析用モデル作成方法である。
【数８】
2016203584.tif 000010
（式（６）において、
Ｇ１２は、対象材料の主軸方向のせん断弾性率であり、
Ｇｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率であり、
ζ１１は、式（６）’で示される値であり、
ηｇは、式（６）’’で示される値であり、
Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。）
【数９】
2016203584.tif 000011
（式（６）’において、
λ１１は、対象材料の主軸方向の配向度であり、
Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。）
【数１０】
2016203584.tif 000012
（式（６）’’において、Ｇｆは、対象材料に含まれるフィラーのせん断弾性率であり、Ｇｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。）
【００１８】
［５］また、本発明は、前記ανは、ＰＶＴ樹脂特性解析によって得られる対象材料の体積膨張率の実測値である、［１］から［４］のいずれかに記載の構造解析用モデル作成方法である。
【００１９】
［６］また、本発明は、前記異方性樹脂成形体がウエルド部を含む、［１］から［５］のいずれかに記載の構造解析用モデル作成方法である。
【００２０】
［７］また、本発明は、前記探索ステップは、前記複数の第１領域の各々について第１重心位置を導出する第１重心位置導出ステップと、前記複数の第２領域の各々について第２重心位置を導出する第２重心位置導出ステップと、前記第２領域の各々について、前記第２重心位置から最も近い第１重心位置を探索し、前記第２重心位置から最も近い第１重心位置を有する第１領域を前記第２領域から最も近い第１領域とする最短第１領域設定ステップとを有する、［１］から［６］のいずれかに記載の構造解析用モデル作成方法である。
【００２１】
［８］また、本発明は、前記配向状態は、対象材料の主軸方向の配向度であり、前記複数種類の仮想成形体は、前記配向度の範囲に応じて１０種類以上１０００種類以下の範囲で設定される、［１］から［７］のいずれかに記載の構造解析用モデル作成方法である。
【００２２】
［９］また、本発明は、［１］から［８］のいずれかに記載の構造解析用モデル作成方法を使用して異方性樹脂成形体を構造解析する構造解析方法であって、前記異方性樹脂成形体について樹脂流動解析を行い、要素分割された前記複数の第１領域の各々における配向状態を算出するための配向状態算出基礎情報を含む第１領域配向情報を取得する第１領域配向情報取得ステップと、前記第１領域配向情報取得ステップとは別に、前記異方性樹脂成形体の構造解析のために前記異方性樹脂成形体を前記複数の第２領域に要素分割する構造解析用要素分割ステップと、前記第２領域配向状態の情報及び前記第２領域物性情報を含む構造解析モデル情報に基づいて前記異方性樹脂成形体の構造解析を行う構造解析ステップとを含み、前記第１領域配向状態算出ステップは、前記配向状態算出基礎情報に基づいて前記複数の第１領域の各々における配向状態に相当する第１領域配向状態を算出するステップであり、前記探索ステップは、前記構造解析用モデルにおける複数の第２領域の各々について、第２領域から最も近い第１領域を探索するステップである、異方性樹脂成形体の構造解析方法である。
【００２３】
［１０］また、本発明は、［９］に記載の構造解析方法を使用することで、前記異方性樹脂射出成形体の反り変形を予測する反り変形予測方法である。
【００２４】
［１１］また、本発明は、異方性フィラーを含有する異方性樹脂成形体の構造解析を行うための構造解析用モデルの作成をコンピュータに実行させるための構造解析用モデル作成プログラムであって、前記異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成する物性情報作成ステップと、前記異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップと、前記異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップと、前記探索ステップで最も近いとされた第１領域における前記第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップと、前記物性情報を参照し、前記複数の第２領域の各々について、前記第２領域配向状態に応じた前記物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップとをコンピュータに実行させるためのプログラムである。
【数１１】
2016203584.tif 000013
（式（１）において、
ν２３は、対象材料の主軸方向に直交する２方向に対するポアソン比であり、
νｆは、対象材料に含まれる異方性フィラーのポアソン比であり、
νｍは、対象材料を構成する樹脂組成物のポアソン比である。）
【数１２】
2016203584.tif 000014
（式（２）において、
ανは、対象材料の体積膨張率であり、
α１１は、対象材料の主軸方向の線膨張率であり、
α２２は、材料主軸方向に直交する方向の線膨張率であり、
α３３は、材料主軸方向に直交する方向の線膨張率であり、
α２２≧α３３  という関係を有する。）
【発明の効果】
【００２５】
本発明によると、異方性樹脂成形体を構造解析するにあたり、高精度の結果を短時間で得られるとともに、コンピュータにかかる負荷を大幅に軽減できる。
【図面の簡単な説明】
【００２６】
【図１】本発明に係る異方性樹脂成形体の構造解析方法の一例を示すフローチャートである。
【図２】図１における構造解析用モデル作成ステップの一例を示すフローチャートである。
【図３】図２における探索ステップの一例を示すフローチャートである。
【図４】一連の解析プログラムを実現するためのハードウェア資源の一例を示す。
【図５】樹脂流動解析における要素分割後の試験片Ａの形状を示す。
【図６】樹脂流動解析によって得られる試験片Ａにおける異方性フィラーの配向方向分布を示す。
【図７】構造解析を行うために要素分割した後の試験片Ａの形状を示す。
【図８】流動解析によって得られる試験片Ａにおける異方性フィラーの配向方向分布を、構造解析の用に供する有限要素分割モデルに当てはめたときの例を示す。
【図９】試験片Ａを実施例１の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。
【図１０】試験片Ｂを比較例１の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。
【図１１】樹脂流動解析における要素分割後の試験片Ｂの形状を示す。
【図１２】樹脂流動解析によって得られる試験片Ｂにおける異方性フィラーの配向方向分布を示す。
【図１３】構造解析を行うために要素分割した後の試験片Ｂの形状を示す。
【図１４】ガラス繊維３０重量％含有強化ポリブチレンテレフタレートの温度と比容積との関係を示す。
【図１５】流動解析によって得られる試験片Ｂにおける異方性フィラーの配向方向分布を、構造解析の用に供する有限要素分割モデルに当てはめたときの例を示す。
【図１６】試験片Ｂにおいて線膨張率を評価した位置を示す。
【図１７】試験片Ｂを実施例２の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。
【図１８】試験片Ｂを比較例２の手法にて構造解析したときに得られる試験片Ａの主応力分布を示す。
【発明を実施するための形態】
【００２７】
以下、本発明の具体的な実施形態について、詳細に説明するが、本発明は、以下の実施形態に何ら限定されるものではなく、本発明の目的の範囲内において、適宜変更を加えて実施することができる。
【００２８】
＜異方性樹脂成形体の構造解析方法＞
図１は、本発明に係る異方性樹脂成形体の構造解析方法の一例を示すフローチャートである。本発明に係る構造解析方法は、異方性フィラーを含有する異方性樹脂成形体について樹脂流動解析を行い、要素分割された複数の第１領域の各々における配向状態を算出するための配向状態算出基礎情報を含む第１領域配向情報を取得する第１領域配向情報取得ステップ（Ｓ１）と、第１領域配向情報取得ステップ（Ｓ１）とは別に、異方性樹脂成形体の構造解析のために異方性樹脂成形体を複数の第２領域に要素分割する構造解析用要素分割ステップ（Ｓ２）と、異方性樹脂成形体の構造解析を行うための構造解析用モデルを作成する構造解析用モデル作成ステップ（Ｓ３）と、構造解析モデルに基づいて異方性樹脂成形体の構造解析を行う構造解析ステップ（Ｓ４）とを含む。
【００２９】
〔第１領域配向情報取得ステップ（Ｓ１）〕
まず、第１領域配向情報取得ステップ（Ｓ１）について説明する。異方性樹脂成形体の配向状態を取得するにあたり、Ｘ線ＣＴ等を用いて、射出成形等の加工によって得られる実際の異方性樹脂成形体を観察することによって取得することも考えられる。しかしながら、試作の前段階である設計段階にある場合、実際の異方性樹脂成形体が試作されていないことから、配向状態を実測するのは難しい。そのため、樹脂流動解析を用い、要素分割された複数の第１領域の各々における配向状態をシミュレートするのが好ましい。
【００３０】
樹脂流動解析の手法は特に限定されるものでなく、例えば、有限要素法を用いた樹脂流動解析等が挙げられる。有限要素法を用いた樹脂流動解析により配向状態を取得する場合、異方性樹脂成形体の肉厚方向に対し、５分割以上で要素分割することが好ましい。要素分割の程度が４分割以下である場合、肉厚方向の繊維配向が平均化され、構造解析の精度が低下し得るため、好ましくない。
【００３１】
配向状態の種類は特に限定されるものでなく、例えば、配向度、配向方向等が挙げられる。
【００３２】
配向度に関する配向度情報は、通常３×３のマトリクスにて表現される。このマトリクスから、３方向の固有ベクトル、固有値を求める。それぞれの固有ベクトルから配向方向が得られる。またそれぞれの固有値から配向度が得られる。
【００３３】
必須ではないが、射出成型の際に２点ゲートを用いると、ウエルドが発生し、ウエルド部に応力が集中するため、異方性樹脂成形体がウエルド部を含むものとしてモデル化することが好ましい。
【００３４】
また、流動解析の精度、ひいては、構造解析の最終的な精度を高めるため、ゲート、ランナー、インサート部等についてもモデル化することが好ましい。
【００３５】
〔構造解析用要素分割ステップ（Ｓ２）〕
続いて、構造解析用要素分割ステップ（Ｓ２）について説明する。構造解析用要素分割ステップ（Ｓ２）は、第１領域配向情報取得ステップ（Ｓ１）とは別に、異方性樹脂成形体の構造解析のために異方性樹脂成形体を複数の第２領域に要素分割するステップである。
【００３６】
構造解析用要素分割ステップ（Ｓ２）は、第１領域配向情報取得ステップ（Ｓ１）での要素分割とは別に要素分割を行うものであれば特に限定されるものでなく、第１領域配向情報取得ステップ（Ｓ１）の後に行ってもよいし、第１領域配向情報取得ステップ（Ｓ１）の前に行ってもよいし、第１領域配向情報取得ステップ（Ｓ１）と並行して行ってもよい。
【００３７】
要素分割する手法は特に限定されるものでない。一例として、まず、ＣＡＤインターフェース等を利用して、異方性樹脂成形体の形状をパソコン等に取り込むか、あるいは、ＣＡＤシステムにより異方性樹脂成形体の形状を作成し、モデル化範囲を設定する。次いで、要素分割プリプロセッサ等で有限要素法等の要素分割を行うことで異方性樹脂成形体を複数の領域に分割する。
【００３８】
要素の形状は特に限定されるものでなく、四面体１次要素、四面体２次要素、六面体１次要素、六面体２次要素等が選択可能であり、有限要素法ソフトウェアの仕様、計算するコンピュータシステムの仕様、計算コスト等に応じて適宜選択すればよい。
【００３９】
要素の数も特に限定されるものでなく、計算精度、計算時間等を考慮して適宜選択すればよい。
【００４０】
〔構造解析用要素分割ステップ（Ｓ３）〕
構造解析用モデル作成ステップ（Ｓ３）は、異方性樹脂成形体の構造解析を行うための予備段階として構造解析用モデルを作成するステップである。以下、図２を参照しながら、構造解析用モデル作成ステップ（Ｓ３）を詳しく説明する。
【００４１】
図２は、図１における構造解析用モデル作成ステップ（Ｓ３）の一例を示すフローチャートであり、本発明に係る構造解析用モデル作成方法に相当する。構造解析用モデル作成ステップ（Ｓ３）は、異方性樹脂成形体と配向状態が異なる複数種類の仮想成形体の各々について、所定のポアソン比及び所定の線膨張率の少なくとも一方を含む物性情報を作成する物性情報作成ステップ（Ｓ３１）と、異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップ（Ｓ３２）と、異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップ（Ｓ３３）と、探索ステップ（Ｓ３３）で最も近いとされた第１領域における第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップ（Ｓ３４）と、物性情報を参照し、複数の第２領域の各々について、第２領域配向状態に応じた物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップ（Ｓ３５）とを含んで構成される。
【００４２】
［物性情報作成ステップ（Ｓ３１）］
物性情報作成ステップ（Ｓ３１）は、配向度等で例示される配向状態に応じた物性情報（機械物性値）を予め作成するステップであり、具体的には、異方性樹脂成形体と配向状態は異なるが、異方性樹脂成形体を構成する樹脂の組成や樹脂成形体の形状等、他の条件は同一である複数種類の仮想成形体の各々について、下記式（１）で規定されるポアソン比と下記式（２）で規定される線膨張率との少なくとも一方を含む物性情報を作成するステップである。物性情報（機械物性値）を予め作成することにより、入力する物性情報の確認、修正が容易となるという利点がある。
【００４３】
物性情報の作成情報は特に限定されるものでなく、剛性マトリクスを構成する物性情報を作成する場合、Ｈａｌｐｉｎ−Ｔｓａｉ式を元に、一部を改良したモデルを用いることが挙げられる。また、作成された物性情報は、ハードディスク、メモリ等に例示される、ハードウェアの記録領域に格納される。
【００４４】
【数１３】
2016203584.tif 000015
【００４５】
式（１）において、ν２３は、対象材料の主軸方向に直交する２方向に対するポアソン比である。また、νｆは、対象材料に含まれる異方性フィラーのポアソン比である。また、νｍは、対象材料を構成する樹脂組成物のポアソン比である。また、Ｖｆは、対象材料に含まれる異方性フィラーの対象材料に対する体積含有率である。
【００４６】
本発明では、ポアソン比を求める方法を改良している。一般的に、各方向のポアソン比を実験的に求める手法として、特開２００９−１２８０３３号公報に記載の手法等が知られているが、配向度依存性を求めることは極めて困難であるため、実際はモデル計算等で求めている。その際、算出方法によっては実際に不可能な値となりやすい。その場合、正定値が負となり、構造解析ソフトウェアによっては計算自体が不可能になる場合がある。式（１）を採用することで、材料物性における正定値が正になりやすくなるため、汎用構造計算ソフトウェアの取扱いが容易になる利点がある。
【００４７】
【数１４】
2016203584.tif 000016
式（２）において、ανは、対象材料の体積膨張率である。また、α１１は、対象材料の主軸方向の線膨張率である。また、α２２、α３３は、材料主軸方向に直交する方向の線膨張率であり、α２２≧α３３という関係性がある。
【００４８】
式（２）において、計算に各方向の線膨張率と体積膨張率との関係を用いたが、この関係を満たさないと、体積の膨張・収縮量と各方向の膨張・収縮量が矛盾する。しかしながら、線膨張率を求める方式にて一般的なＳｃｈａｐｅｒｙモデル、Ｈａｓｈｉｎ−Ｓｈｔｒｉｋｍａｎモデル等ではこの関係を必ずしも満たさないため、計算結果の精度が悪い原因になっていた。
【００４９】
α１１は、下記式（３）によって求められることが好ましく、α２２は、下記式（４）によって求められることが好ましい。
【数１５】
2016203584.tif 000017
【００５０】
式（３）において、αｓは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率である。また、αｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、式（３）’で規定される。λ１１は、対象材料の主軸方向の配向度である。ｘは、配向度に対する各方向の非線形度を示すパラメーターである。
【００５１】
また、式（３）’において、ανは、対象材料の体積膨張率である。
【００５２】
【数１６】
2016203584.tif 000018
【００５３】
式（４）において、αｂは、完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率である。αｒは、配向がランダム状態の仮想異方性樹脂成形体の主軸方向の線膨張率であり、式（４）’で規定される。λ１１は、対象材料の主軸方向の配向度である。ｘは、配向度に対する各方向の非線形度を示すパラメーターである。
【００５４】
また、式（４）’において、ανは、対象材料の体積膨張率である。
【００５５】
本発明の実施形態では、仮想異方性樹脂成形体の主軸方向の配向度に応じて連続的に材料物性値を与えるのではなく、区分的に扱う。具体的に、仮想異方性樹脂成形体の主軸方向の配向度λ１１は、１／３以上１以下の値をとり得るため、例えば、λ１１が０．９以上１以下の範囲内にある場合を一種類の材料物性を持つものとする。この区分については、区間が細かいほど連続的になり、解析精度がより高まるが、区分数が多くなるため、材料物性の確認／変更は困難になる。構造解析での解析精度、及び材料物性の確認／変更の容易性の双方を考慮すると、複数種類の仮想成形体は、前記配向度の範囲に応じて１０種類以上１０００種類以下の範囲で設定されることが好ましい。
【００５６】
物性情報として、上記式（１）で表されるポアソン比、及び式（２）で表される線膨張率のほか、下記式（５）で規定される弾性率、下記式（６）で規定されるせん断弾性率等が挙げられる。
【００５７】
【数１７】
2016203584.tif 000019
【００５８】
【数１８】
2016203584.tif 000020
【００５９】
式（５）において、Ｅ１１は、対象材料の主軸方向の弾性率であり、Ｅｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率であり、ζ１１は、式（５）’で示される値であり、ηは、式（５）’’で示される値であり、Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。
【００６０】
【数１９】
2016203584.tif 000021
【００６１】
式（５）’において、λ１１は、対象材料の主軸方向の配向度であり、Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。
【００６２】
【数２０】
2016203584.tif 000022
式（５）’’において、Ｅｆは、対象材料に含まれるフィラーの弾性率であり、Ｅｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。
【００６３】
式（６）において、Ｇ１２は、対象材料の主軸方向のせん断弾性率であり、Ｇｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率であり、ζ１１は、式（６）’で示される値であり、ηｇは、式（６）’’で示される値であり、Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。
【００６４】
【数２１】
2016203584.tif 000023
【００６５】
式（６）’において、λ１１は、対象材料の主軸方向の配向度であり、Ｌ／ｄは、対象材料に含まれる異方性フィラーのアスペクト比である。
【００６６】
【数２２】
2016203584.tif 000024
式（６）’’において、Ｇｆは、対象材料に含まれるフィラーのせん断弾性率であり、Ｇｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。
【００６７】
式（５）’及び式（６）’に関し、配向度等で例示される配向状態を元に弾性率等を求める手法は、これまでも均質化法として知られている。しかしながら、均質化法は、弾性率等を求める際の計算方法が複雑であるため、計算エラー、プログラムコーディングエラー等を生じ得る。計算手順、計算プログラム作成を容易にするため、微小領域における配向度を考慮できるモデルとして、上記の式（５）’及び式（６）’を採用することが好ましい。
【００６８】
式（５）’及び式（６）’において、アスペクト比Ｌ／ｄは、実測値を採用することが好ましいが、異方性樹脂成型体に含まれるフィラーのアスペクト比は、フィラーの破損等により分布があるとともに、射出成型前の状態からは変化する。したがって、式（５）’及び式（６）’の運用に関しては、アスペクト比をパラメーターとして扱ってもよい。
【００６９】
その他、物性情報として、以下の式（７）から式（９）を例示できる。
【００７０】
【数２３】
2016203584.tif 000025
【００７１】
式（７）において、ν１２及びν１３は、対象材料の主軸方向に直交する方向に対するポアソン比である。また、Ｖｆは、異方性フィラーの対象材料に対する体積含有率であり、νｆは、対象材料に含まれる異方性フィラーのポアソン比である。また、Ｖｍは、樹脂組成物の対象材料に対する体積含有率であり、νｍは、対象材料を構成する樹脂組成物のポアソン比である。
【００７２】
【数２４】
2016203584.tif 000026
【００７３】
式（８−１）及び式（８−２）において、Ｅ２２及びＥ３３は、対象材料の主軸方向に直交する方向の弾性率であり、Ｅｍは、対象材料を構成する樹脂組成物からなる成形体の弾性率である。また、ζ２２は、２であり、ηは、上記式（５）’’で示される値であり、Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。
【００７４】
【数２５】
2016203584.tif 000027
【００７５】
式（９−１）及び式（９−２）において、Ｇ２３は、対象材料の主軸方向に直交する方向のせん断弾性率であり、Ｇｍは、対象材料を構成する樹脂組成物からなる成形体のせん断弾性率である。また、ζ２２は、２であり、ηｇは、式（６）’’で示される値であり、Ｖｆは、異方性フィラーの対象材料に対する体積含有率である。
【００７６】
［第１領域配向状態算出ステップ（Ｓ３２）］
第１領域配向状態算出ステップ（Ｓ３２）は、異方性樹脂成形体の樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出するステップである。本実施形態では、配向状態が配向方向及び配向度であるものとして説明するが、これに限るものではない。
【００７７】
配向度に関する配向度情報は、通常３×３のマトリクスにて表現される。このマトリクスから、固有ベクトル、固有値を求める。固有ベクトルから配向方向が得られる。また固有値から配向度が得られる。固有ベクトルの３方向を座標系に変換し、後のステップにて各要素に与える。また、配向度に応じて物性情報作成ステップ（Ｓ３１）にて得られた物性情報を後のステップにて各要素に設定する。これらの設定は、構造解析の前後において確認しやすいため、物性設定における不具合を確認しやすくなるという効果がある。
【００７８】
［探索ステップ（Ｓ３３）］
探索ステップ（Ｓ３３）は、異方性樹脂成形体の構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索するステップである。探索ステップ（Ｓ３３）の態様は特に限定されるものでないが、例えば、図３に示すように、複数の第１領域の各々について第１重心位置を導出する第１重心位置導出ステップ（Ｓ３３１）と、複数の第２領域の各々について第２重心位置を導出する第２重心位置導出ステップ（Ｓ３３２）と、第２領域の各々について、第２重心位置から最も近い第１重心位置を探索し、第２重心位置から最も近い第１重心位置を有する第１領域を第２領域から最も近い第１領域とする最短第１領域設定ステップ（Ｓ３３３）とを有する態様を例示できる。
【００７９】
（第１重心位置導出ステップ（Ｓ３３１））
図３を参照しながら、探索ステップ（Ｓ３３）の具体例について説明する。第１重心位置導出ステップ（Ｓ３３１）では、第１領域配向情報取得ステップ（Ｓ１）にて計算に用いた第１領域に係る要素毎に、要素を構成する節点の位置から第１重心位置を求める。第１重心位置は、次に計算する構造解析用モデルの第２重心位置と比較される。
【００８０】
（第２重心位置導出ステップ（Ｓ３３２））
第２重心位置導出ステップ（Ｓ３３２）では、第１重心位置導出ステップ（Ｓ３３１）と同様に、構造解析用要素分割ステップ（Ｓ２）で分割された第２領域の要素毎に、要素を構成する節点の位置から重心位置を求める。
【００８１】
（最短第１領域設定ステップ（Ｓ３３３））
最短第１領域設定ステップ（Ｓ３３３）では、構造解析用要素分割ステップ（Ｓ２）で分割された第２領域の要素毎に、第２重心位置導出ステップ（Ｓ３３２）で求めた第２重心位置と、第１重心位置導出ステップ（Ｓ３３１）で求めた第１重心位置との間の距離を計算し、最も近い位置の情報を当てはめる。この最短第１領域設定ステップ（Ｓ３３３）を終えると、図２の第２領域配向状態設定ステップ（Ｓ３４）に移る。
【００８２】
［第２領域配向状態設定ステップ（Ｓ３４）］
図２に戻る。第２領域配向状態設定ステップ（Ｓ３４）は、探索ステップ（Ｓ３３）で最も近いとされた第１領域における第１領域配向状態を各々の第２領域配向状態とするステップである。この方式の場合、第１領域配向情報取得ステップ（Ｓ１）にて用いた樹脂流動解析での要素分割が粗く、構造解析用要素分割ステップ（Ｓ２）で分割された第２領域の要素の第２重心位置に第１領域配向情報取得ステップ（Ｓ１）にて用いた第１領域に係る要素がない場合でも、第１領域配向状態の情報を第２領域に当てはめることができる。
【００８３】
この方式に従えば、要素分割状態によらず第１領域配向状態の情報を第２領域に当てはめることが可能である。この当てはめは、要素の形状関数を用いて配向状態の情報を補間して求めることも可能であるが、要素分割が異なる場合には補間できない状況になる場合があるため、配向状態の情報を求める際の情報密度、構造解析の要素分割に制約が生じる場合があり、その制約により、構造解析自体の精度に影響を及ぼす場合がある。
【００８４】
［第２領域物性情報設定ステップ（Ｓ３５）］
第２領域物性情報設定ステップ（Ｓ３５）は、物性情報作成ステップ（Ｓ３１）で予め作成され、ハードウェアの記録領域に格納されている物性情報を参照し、複数の第２領域の各々について、第２領域配向状態に応じた物性情報を各々の第２領域物性情報とするステップである。このステップでは、第２領域配向状態設定ステップ（Ｓ３４）にて求めた第２領域配向状態の情報に基づいた座標系、及び配向度に応じた材料物性を構造解析の要素として設定する。この第２領域物性情報設定ステップ（Ｓ３５）を終えると、構造解析の実行が可能になるため、図１の構造解析ステップ（Ｓ４）に移る。
【００８５】
〔構造解析ステップ（Ｓ４）〕
構造解析ステップ（Ｓ４）は、構造解析モデルに基づいて異方性樹脂成形体の構造解析を行うステップである。このステップでは、第２領域配向状態設定ステップ（Ｓ３４）にて設定した第２領域配向状態の情報と、第２領域物性情報設定ステップ（Ｓ３５）にて設定した第２領域物性情報とを含む構造解析モデル情報を用い、構造解析ソフトウェアのプログラムを実行して構造解析を実行する。構造解析における拘束条件、荷重条件、計算条件等の各種条件は、構造解析ソフトェアの仕様にしたがえばよく、特に限定されない。
【００８６】
構造解析により、異方性樹脂成形体の変形、応力等をシミュレートできる。このシミュレートの結果から、異方性樹脂成形体における破損の有無等を判定できる。
【００８７】
具体的態様として、異方性樹脂射出成形体の反り変形は、構造解析における荷重条件を温度荷重及び残留応力荷重に設定することで解析可能である。そのため、本発明に係る構造解析方法を使用することで、異方性樹脂射出成形体の反り変形を予測できる。
【００８８】
＜解析プログラム＞
本発明に係る構造解析用モデル作成プログラムを含む一連の解析プログラムは、ソフトウェアとハードウェア資源とが協働することによって実現される。
【００８９】
図４は、一連の解析プログラムを実現するためのハードウェア資源Ｈの一例を示す。ハードウェア資源Ｈは、情報処理装置１と、設計者からの各種要求を受け付ける入力装置２と、情報処理装置１が行った解析結果を出力する出力装置３とを備える。また、情報処理装置１は、ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）等のネットワークＮＷを介して、ＣＡＤ装置４に接続されている。
【００９０】
情報処理装置１は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１０と、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）等により構成される主記憶装置２０と、Ｉ／Ｏインタフェース３０と、ハードディスク等により構成される補助記憶装置４０と、ネットワークＮＷに接続されている装置との間で行うデータ授受の制御を行うネットワークインタフェース５０とを備える。
【００９１】
補助記憶装置４０には、上述した一連のステップをハードウェア装置１に実行させるための解析プログラム４１が格納されている。解析プログラム４１は、流動解析による第１領域配向情報取得ステップ（Ｓ１）をハードウェア装置１に実行させるための第１領域配向情報取得プログラム４１Ａと、構造解析用要素分割ステップ（Ｓ２）をハードウェア装置１に実行させるための構造解析用要素分割プログラム４１Ｂと、構造解析用モデル作成ステップ（Ｓ３）をハードウェア装置１に実行させるための構造解析用モデル作成プログラム４１Ｃと、構造解析ステップ（Ｓ４）をハードウェア装置１に実行させるための構造解析プログラム４１Ｄとを含んで構成される。本発明に係る構造解析方法は、ＣＰＵ１０が補助記憶装置４０に格納されている解析プログラム４１を主記憶装置２０にロードして実行することにより実現される。特に、本発明に係る異方性樹脂成形体の構造解析用モデル作成方法は、ＣＰＵ１０が補助記憶装置４０に格納されている構造解析用モデル作成プログラム４１Ｃを主記憶装置２０にロードして実行することにより実現される。
【００９２】
上記の実施形態では、第１領域配向情報取得ステップ（Ｓ１）から構造解析ステップ（Ｓ４）に至るまでの各ステップは、複数のプログラムが組み合わせられて実行されることで行われているが、これに限られることなく、最初から一体として構築されたプログラムでもよく、また、実行されるコンピュータの形態や規模、設置場所等も限定されるものではない。
【実施例】
【００９３】
以下、実施例及び比較例を示し、本発明を具体的に説明するが、本発明はこれらの実施例に限定されるものではない。
【００９４】
＜試験例１＞
試験例１では、平板部分の長片側にリブを設置した形状であって、射出成型にて得られる試験片Ａを用いて、３点曲げ試験時における耐荷重を検討した。
【００９５】
〔実施例１〕
まず、この試験片Ａに対し、ＪＩＳ１７８に準拠した曲げ試験を行った。固定冶具間の間隔を６２ｍｍに設定し、試験片中央部に圧子を配置し、圧子をひずみ速度が１％／ｍｉｎにて移動させ、試験機器に搭載のロードセンサーにて荷重を測定した。その結果、試験片は荷重２７２Ｎにて破損した。
【００９６】
［第１領域配向情報取得ステップ（Ｓ１）］
続いて、試験片Ａの形状をＣＡＤデータとして作成し、ランナー、ゲート等を加え有限要素分割を実行した。樹脂流動解析には、ＡＵＴＯＤＥＳＫ社ＡＵＴＯＤＥＳＫ  ＳＩＭＵＬＡＴＩＯＮ  ＭＯＬＤＦＬＯＷ  ＩＮＳＩＧＨＴを用いた。樹脂流動解析の条件を以下に示す。
樹脂：ガラス繊維３０重量％含有強化ポリブチレンテレフタレート
樹脂温度：２６０℃
金型温度：８０℃
射出流量：５７．７ｃｍ３／ｓ
保圧圧力：７０ＭＰａ        
保圧時間：１５秒
冷却時間：１０秒
要素：四面体１次要素（分割数  ４５０５３８）
【００９７】
樹脂流動解析における要素分割後の試験片Ａの形状を図５に示し、樹脂流動解析によって得られる試験片Ａにおける異方性フィラーの配向方向分布を図６に示す。
【００９８】
［構造解析用要素分割ステップ（Ｓ２）］
流動解析用の有限要素モデルとは別に、異方性樹脂成形体の構造解析のために、試験片Ａを複数の第２領域に要素分割した。要素分割後の試験片Ａの形状を図７に示す。なお、図７においては、ランナー、ゲート等を削除し、代わりに試験片支持台、および圧子を加えた。要素については四面体２次要素を用い、分割数は６６５６５であった。試験片支持台を完全拘束し、圧子に、試験片の耐荷重である荷重２７２Ｎを与え、最大主応力および発生位置を検討した。
【００９９】
［構造解析用モデル作成ステップ（Ｓ３）］
まず、物性情報作成ステップ（Ｓ３１）において、試験片Ａと配向状態が異なる複数種類の仮想成形体の各々について、上記式（１）〜（９）にしたがって物性情報を作成し、コンピュータの所定の記憶領域に格納した。式（１）〜（９）にしたがって得られる物性情報は、以下のとおりであった。
対象材料に含まれるフィラーの弾性率Ｅｆ                ：７２０００ＭＰａ
対象材料を構成する樹脂組成物からなる成形体の弾性率Ｅｍ：２５００ＭＰａ
異方性フィラーの対象材料に対する体積含有率Ｖｆ        ：１８ｖｏｌ％
樹脂組成物の対象材料に対する体積含有率Ｖｍ            ：８２ｖｏｌ％
対象材料に含まれるフィラーのせん断弾性率Ｇｆ          ：２９５００ＭＰａ
対象材料を構成する樹脂組成物からなる成形体のせん断弾性率Ｇｍ：９２６ＭＰａ
対象材料に含まれる異方性フィラーのポアソン比νｆ      ：０．２２
対象材料を構成する樹脂組成物のポアソン比νｍ          ：０．３５
対象材料に含まれる異方性フィラーのアスペクト比Ｌ／ｄ  ：１４
【０１００】
そして、試験片Ａの樹脂流動解析において要素分割された複数の第１領域の各々における第１領域配向状態を算出する第１領域配向状態算出ステップ（Ｓ３２）と、試験片Ａの構造解析のために要素分割された複数の第２領域の各々について、第２領域から最も近い第１領域を探索する探索ステップ（Ｓ３３）と、探索ステップ（Ｓ３３）で最も近いとされた第１領域における第１領域配向状態を各々の第２領域配向状態とする第２領域配向状態設定ステップ（Ｓ３４）と、物性情報を参照し、複数の第２領域の各々について、第２領域配向状態に応じた物性情報を各々の第２領域物性情報とする第２領域物性情報設定ステップ（Ｓ３５）とを経て構造解析モデルを完成させる。完成した構造解析モデルを図８に示す。元となる図６の分布状況と比較すると、方向を示す位置が異なるが、両者の方向が一致している。
【０１０１】
［構造解析ステップ（Ｓ４）］
得られた構造解析用のモデルを用いて計算する。計算には、構造解析ソフトＮＸ  Ｉ−ＤＥＡＳ（シーメンスＰＬＭソフトウェア社製）を用いた。
【０１０２】
図９は、試験片Ａを実施例１の手法にて構造解析することによって得られる主応力分布を示す。最大応力を示す位置は、試験片の中央部から若干ずれている。試験片の成形の際には２点ゲートを用いたので、ウエルドが発生し、ウエルド部に応力が集中した結果と考えられる。また、実施例１における最大主応力は１３９ＭＰａであった（表１）。
【０１０３】
また、試験片の引張破壊強度の実測値は１４０ＭＰａであり、ウエルド部が破壊の起点となっていた（表１）。
【０１０４】
このように、実施例１の結果は、実測に極めて近く、実施例１において本発明の精度は極めて良好といえる。
【０１０５】
〔比較例１〕
実施例１と同一の試験片について、弾性率として、流動方向の弾性率と流動垂直方向の弾性率とを平均化した値である６３２０ＭＰａを用い、ポアソン比として、既存のモデル計算によって求められる０．３９を用いたこと以外は、実施例１と同じ手法にて、試験片Ａの構造解析を行った。
【０１０６】
図１０は、試験片Ａを比較例１の手法にて構造解析することによって得られる主応力分布を示す。最大応力を示す位置は試験片中央部であり、ウエルド発生位置とは異なる位置であった。比較例１における最大主応力は１０６ＭＰａであった（表１）。
【０１０７】
このように、比較例１の結果は、実測と異なり、比較例１の精度は、実施例１に比べて劣るといえる。
【０１０８】
【表１】
2016203584.tif 000028
【０１０９】
＜試験例２＞
試験例２では、平板部分に箱型形状、リブ、円筒等を設置した形状であって、射出成型にて得られる試験片Ｂを用いて、熱膨張挙動を検討した。
【０１１０】
〔実施例２〕
まず、射出成型にて得られた試験片Ｂを切り出し、１２０℃環境にて恒温槽にてアニールを２時間程実施し、その後放冷した。その後、線膨張率測定装置に設置し、−３０℃から室温まで温度を上昇させ、４か所において寸法変化量から、線膨張率を求めた。その結果は後述の表２にまとめた。
【０１１１】
［繊維配向情報取得ステップ（Ｓ１）］
続いて、実施例１と同様の手法にて、試験片Ｂの有限要素分割を実行した。樹脂流動解析の条件を以下に示す。
樹脂：ガラス繊維３０重量％含有強化ポリブチレンテレフタレート
樹脂温度：２５０℃
金型温度：６０℃
射出流量：３２ｃｍ３／ｓ
保圧圧力：７８．４ＭＰａ
保圧時間：１０秒
冷却時間：１５秒
要素：四面体１次要素（分割数  ４８６７４９）
【０１１２】
樹脂流動解析における要素分割後の試験片Ｂの形状を図１１に示し、樹脂流動解析によって得られる試験片Ｂにおける異方性フィラーの配向方向分布を図１２に示す。
【０１１３】
［構造解析用要素分割ステップ（Ｓ２）］
流動解析用の有限要素モデルとは別に、異方性樹脂成形体の構造解析のために、試験片Ｂを複数の第２領域に要素分割した。要素分割後の試験片Ｂの形状を図１３に示す。なお、図１３においては、切り出し後の形状を用いた。要素については四面体２次要素を用い、分割数は６６５６５であった。モデルのコーナー部３ヶ所にＸＹＺ方向、ＹＺ方向、Ｚ方向を拘束し、温度荷重１００℃を与えた。
【０１１４】
［材料物性の設定ステップ（Ｓ３）］
まず、物性情報作成ステップ（Ｓ３１）において、試験片Ａと配向状態が異なる複数種類の仮想成形体の各々について、上記式（１）〜（９）にしたがって物性情報を作成し、コンピュータの所定の記憶領域に格納した。式（１）〜（９）にしたがって得られる物性情報は、以下のとおりであった。
完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向の線膨張率αｓ    ：２．４×１０−５／℃
完全配向である仮想異方性樹脂成形体、すなわち配向度が１である仮想異方性樹脂成形体の主軸方向に直交する方向の線膨張率αｂ    ：６．２×１０−５／℃
対象材料の体積膨張率αν    ：１．５３×１０−４／℃
配向度に対する各方向の非線形度ｘ    ：１５．５（無次元）
【０１１５】
体積膨張率ανは、ＰＶＴ測定装置を用い、圧力が０ＭＰａであるときのガラス繊維３０重量％含有強化ポリブチレンテレフタレート（強化ＰＢＴ）の体積変化を求め、強化ＰＢＴの温度に対する比容積の変化率を計算することによって求めた値である。
【０１１６】
そして、実施例１と同様に、第１領域配向状態算出ステップ（Ｓ３２）と、探索ステップ（Ｓ３３）と、第２領域配向状態設定ステップ（Ｓ３４）と、第２領域物性情報設定ステップ（Ｓ３５）とを経て構造解析モデルを完成させる。完成した構造解析モデルを図１５に示す。元となる図１２の分布状況と比較すると、方向を示す位置が異なるが、両者の方向が一致している。
【０１１７】
［構造解析ステップ（Ｓ４）］
続いて、実施例１と同様の手法によって試験片Ｂを構造解析した。
【０１１８】
図１６に示す位置Ａ，Ｂ，Ｃ，Ｄにおける寸法変化から、各位置における線膨張率を求めた。結果を表２に示す。また、図１７は、試験片Ｂを実施例２の手法にて構造解析することによって得られる変形量分布を示す。図１７から、試験片Ｂは、平面部がゆがむように変形することを示唆している。
【０１１９】
また、試験片Ｂにおいて、位置Ａ，Ｂ，Ｃ，Ｄにおける線膨張率を実測した。結果を表２に示す。この実測の際、試験片Ｂの平面部がゆがむように変形することを目視にて確認した。
【０１２０】
このように、実施例２における位置Ａ，Ｂ，Ｃ，Ｄ間の線膨張率差の変化は、実測値を反映しているといえ、実際の試験片Ｂのゆがみを反映しているといえる。したがって、実施例２の結果は、実測に極めて近く、実施例２において本発明の精度は極めて良好といえる。
【０１２１】
〔比較例２〕
実施例２と同一の試験片Ｂについて、線膨張率として、流動方向の線膨張率と流動垂直方向の線膨張率とを平均化した値である４．３×１０−５／℃を用いたこと以外は、実施例２と同じ手法にて、試験片Ｂの構造解析を行った。
【０１２２】
図１６に示す位置Ａ，Ｂ，Ｃ，Ｄにおける寸法変化から、各位置における線膨張率を求めた。結果を表２に示す。比較例２の場合、位置Ａ，Ｂ，Ｃ，Ｄ間で線膨張率の有意差は認められない。また、図１８は、試験片Ｂを比較例２の手法にて構造解析することによって得られる変形量分布を示す。図１８から、試験片Ｂの平面部が平坦のままであり、実際の試験片Ｂのゆがみを反映していないことが明らかである。
【０１２３】
このように、比較例２の結果は、実測と異なり、比較例２の構造解析精度は、実施例２に比べて劣るといえる。
【０１２４】
【表２】
2016203584.tif 000029
【符号の説明】
【０１２５】
Ｈ    ハードウェア資源
１    情報処理装置
２    入力装置
３    出力装置
４    ＣＡＤ装置
１０  ＣＰＵ
２０  主記憶装置
３０  Ｉ／Ｏインタフェース
４０  補助記憶装置
４１  解析プログラム
４１Ａ  第１領域配向情報取得プログラム
４１Ｂ  構造解析用要素分割プログラム
４１Ｃ  構造解析用モデル作成プログラム
４１Ｄ  構造解析プログラム
５０  ネットワークインタフェース
ＮＷ  ネットワーク
Ｓ１  第１領域配向情報取得ステップ
Ｓ２  構造解析用要素分割ステップ
Ｓ３  構造解析用モデル作成ステップ
Ｓ４  構造解析ステップ
Ｓ３１  物性情報作成ステップ
Ｓ３２  第１領域配向状態算出ステップ
Ｓ３３  探索ステップ
Ｓ３４  第２領域配向状態設定ステップ
Ｓ３５  第２領域物性情報設定ステップ
Ｓ３３１  第１重心位置導出ステップ
Ｓ３３２  第２重心位置導出ステップ
Ｓ３３３  最短第１領域設定ステップ
【図１】
2016203584.tif 000042
【図２】
2016203584.tif 000043
【図３】
2016203584.tif 000044
【図４】
2016203584.tif 000045
【図１４】
2016203584.tif 000046
【図１６】
2016203584.tif 000047
【図５】
2016203584.tif 000048
【図６】
2016203584.tif 000049
【図７】
2016203584.tif 000050
【図８】
2016203584.tif 000051
【図９】
2016203584.tif 000052
【図１０】
2016203584.tif 000053
【図１１】
2016203584.tif 000054
【図１２】
2016203584.tif 000055
【図１３】
2016203584.tif 000056
【図１５】
2016203584.tif 000057
【図１７】
2016203584.tif 000058
【図１８】
2016203584.tif 000059