"繊維強化樹脂射出成形品の固有振動数の推定方法であって、前記推定方法が、前記繊維強化樹脂射出成形品の樹脂流動解析を実行して、前記繊維強化樹脂射出成形品の弾性パラメータ及び繊維配向パラメータを算出するステップI、前記ステップIで算出した弾性パラメータ及び繊維配向パラメータを、前記繊維強化樹脂射出成形品の理論モード解析に導入して、前記理論モード解析を実行して、固有振動数の算出値を得るステップIIを含む繊維強化樹脂射出成形品の固有振動数の推定方法。前記ステップIが、樹脂流動解析モデルを入力するステップI-1、及び前記繊維強化樹脂射出成形品について有限要素法に基づく要素分割、及び、樹脂流動解析モデルに対する境界条件として、前記繊維強化樹脂射出成形品の製造時の射出条件並びに樹脂流動ゲートの種類を入力して前記樹脂流動解析を実行するステップI-2を含む請求項1記載の繊維強化樹脂射出成形品の固有振動数の推定方法。前記樹脂流動解析モデルとして、弾性パラメータを算出するためのTandon-Wengマイクロメカニクスモデル、及び、繊維配向パラメータを算出するためのFloger-Tuckerモデルを入力する請求項2記載の繊維強化樹脂射出成形品の固有振動数の推定方法。前記射出条件として、樹脂温度°C、金型温度°C、射出速度mm/s、保持圧力MPa、保圧時間秒及び冷却時間秒を入力する請求項2又は3記載の繊維強化樹脂射出成形品の固有振動数の推定方法。前記ステップIにおいて、樹脂流動解析を実行して、前記弾性パラメータとして、流れ方向の弾性率E1、垂直方向の弾性率E2、ポアソン比ν、及び、せん断弾性率τ12、τ23並びにτ31を算出し、前記繊維配向パラメータとして、流れ方向角度θを算出する請求項1~4のいずれか1項記載の固有振動数の推定方法。前記樹脂流動解析を、樹脂流動解析ソフトウェアAutodeskMoldflowInsights2010によって実行する請求項1~5のいずれか1項記載の繊維強化樹脂射出成形品の固有振動数の推定方法。前記ステップIIが、前記ステップIで算出された弾性パラメータ及び繊維配向パラメータを入力するステップステップII-1、及び前記繊維強化樹脂射出成形品について有限要素法に基づく要素分割を入力して前記理論モード解析を実行するステップII-2を含む請求項1~6のいずれか1項記載の繊維強化樹脂射出成形品の固有振動数の推定方法。前記理論モード解析を、構造解析ソフトウェアNXNastranによって実行する請求項1~7のいずれか1項記載の繊維強化樹脂射出成形品の固有振動数の推定方法。前記繊維強化樹脂射出成形品が、ポリアミド樹脂及びガラス繊維を含むポリアミド樹脂組成物を射出成形して得られる請求項1~8のいずれか1項記載の繊維強化樹脂射出成形品の固有振動数の推定方法。繊維強化樹脂射出成形品の製造方法であって、請求項1~9のいずれか1項記載の繊維強化樹脂射出成形品の固有振動数の推定方法で繊維強化樹脂射出成形品の固有振動数を算出する工程を含む繊維強化樹脂射出成形品の製造方法。しかし、繊維強化樹脂射出成形品を製造して実験モード解析を行うことは、正確な固有振動数を算出きる一方で、時間と労力を要するので、理論モード解析によって効率的、かつ、正確に固有振動数を算出できることが望ましい。本発明は、射出成形によって得られる繊維強化樹脂射出成形品の固有振動数を理論モード解析によって推定できる解析方法を提供することを課題とする。本発明は、1繊維強化樹脂射出成形品の固有振動数の推定方法であって、前記推定方法が、前記繊維強化樹脂射出成形品の樹脂流動解析を実行して、前記繊維強化樹脂射出成形品の弾性パラメータ及び繊維配向パラメータを算出するステップI、前記ステップIで算出した弾性パラメータ及び繊維配向パラメータを、前記繊維強化樹脂射出成形品の理論モード解析に導入して、前記理論モード解析を実行して、固有振動数の算出値を得るステップIIを含む繊維強化樹脂射出成形品の固有振動数の推定方法、及び2繊維強化樹脂射出成形品の製造方法であって、前記1記載の繊維強化樹脂射出成形品の固有振動数の推定方法で繊維強化樹脂射出成形品の固有振動数を算出する工程を含む繊維強化樹脂射出成形品の製造方法である。本発明によれば、射出成形によって得られる繊維強化樹脂射出成形品の固有振動数を理論モード解析によって推定できる解析方法を提供することができる。本発明の推定方法は、繊維強化樹脂射出成形品の樹脂流動解析を実行して、前記繊維強化樹脂射出成形品の弾性パラメータ及び繊維配向パラメータを算出するステップIを含む。本発明の推定方法は、ステップIで、繊維強化樹脂射出成形品の弾性パラメータと繊維の配向状態を反映する繊維配向パラメータを算出して、これらのパラメータを、後述するステップIIにおける理論モード解析に導入して、繊維強化樹脂射出成形品の固有振動数を算出する点に特徴を有する。本発明の推定方法における繊維強化樹脂射出成形品は、実際に製造した繊維強化樹脂射出成形品でもよいし、射出成形条件と材料の物性を指定することによってモデル化された仮想的な繊維強化樹脂射出成形品でもよい。これらの場合、本発明の推定方法で算出された固有振動数を有する繊維強化樹脂射出成形品をモデルにして、より効果的に、実際に開発可能な繊維強化樹脂射出成形品を設計して製造することができる。具体的には、繊維強化樹脂射出成形品の固有振動数が予測できれば、共振での破壊回避等、効果的に振動、音対策が可能となる。例えば、繊維強化樹脂射出成形品のある振動モードで最も振幅が大きい部位に対し、板厚増し等の剛性を上げるとか金具で固定する等の具体的対策を行うことで、振幅を抑制し振動問題を解決した繊維強化樹脂射出成形品を設計して製造することができる。即ち、本発明は、繊維強化樹脂射出成形品の製造方法であって、本発明の繊維強化樹脂射出成形品の固有振動数の推定方法で繊維強化樹脂射出成形品の固有振動数を算出する工程を含む繊維強化樹脂射出成形品の製造方法を開示する。ステップIは、好ましくは、樹脂流動解析モデルを入力するステップI-1、及び前記繊維強化樹脂射出成形品について有限要素法に基づく要素分割、及び、樹脂流動解析モデルに対する境界条件として、前記繊維強化樹脂射出成形品の製造時の射出条件並びに樹脂流動ゲートの種類を入力して前記樹脂流動解析を実行するステップI-2を含む。ステップI-1における樹脂流動解析モデルは、繊維強化樹脂射出成形品の射出成形条件、射出成形する樹脂組成物の樹脂及び繊維の種類及び物性等に応じて、適したモデルを選択できるが、Tandon-Wengマイクロメカニクスモデルを、弾性パラメータを算出するためのモデルとして、希薄溶液中の楕円体配向挙動の予測モデルであるJefferyモデル非特許文献2を高濃度状態に拡張したFloger-Tuckerモデルを、繊維配向パラメータを算出するためのモデルとして採用することが好ましい。ステップI-2における有限要素法に基づく要素分割は、繊維強化樹脂射出成形品の形状及び大きさに応じて適した要素分割が選択される。ステップI-3で入力する射出条件としては、樹脂温度°C、金型温度°C、射出速度mm/s、保持圧力MPa、保圧時間秒及び冷却時間秒が好ましく挙げられる。ステップI-2における樹脂流動ゲートは、繊維強化樹脂射出成形品の形状及び大きさに応じて、射出成形時の樹脂の流動状態が適切にモデル化されるように選択される。例えば、単純な直方体で近似できる形状で、小型の繊維強化樹脂射出成形品であれば、1点ゲートを選択できる場合があり、複雑な形状で、大型の繊維強化樹脂射出成形品であれば、複数のゲートを選択することが好ましい。樹脂流動解析は、現在では、様々な市販の解析ソフトウェア例えば、樹脂流動解析ソフトウェア3DTIMON東レエンジニアリング社製、樹脂射出成形用流動解析ソフトウェアMoldex3Dセイロジャパン社製、樹脂流動解析ソフトウェアMoldflowInsights2010Autodesk社製等を利用できるが、デファクトスタンダードの観点から、樹脂流動解析ソフトウェアMoldflowInsights2010Autodesk社製によって実行することが好ましい。また、後述するステップIIにおける理論モード解析には、振動特性への繊維配向影響を反映させる観点から、ステップIの樹脂流動解析を実行して、弾性パラメータとして、MD方向の弾性率E1、TD方向の弾性率E2、ポアソン比ν、及び、せん断弾性率τ12、τ23並びにτ31を算出し、繊維配向パラメータとして、MD方向角度θを算出することが好ましい。本発明の推定方法は、ステップIで算出した弾性パラメータ及び繊維配向パラメータを、前記繊維強化樹脂射出成形品の理論モード解析に導入して、前記理論モード解析を実行して、固有振動数の算出値を得るステップIIを含む。ステップIIは、好ましくは、前記ステップIで算出された弾性パラメータ及び繊維配向パラメータを入力するステップステップII-1、及び前記繊維強化樹脂射出成形品について有限要素法に基づく要素分割を入力して前記理論モード解析を実行するステップII-2を含む。ステップIで算出された弾性パラメータ及び繊維配向パラメータは、例えば、後述する構造解析ソフトウェアNXNastranを使用する場合、弾性パラメータとして、MD方向の弾性率E1、TD方向の弾性率E2、ポアソン比ν、及び、せん断弾性率τ12、τ23並びにτ31を、繊維配向パラメータとして、MD方向角度θを入力することが好ましい。ステップII-2における有限要素法に基づく要素分割は、繊維強化樹脂射出成形品の形状及び大きさに応じて適した要素分割が選択される。理論モード解析は、現在では、様々な市販の構造解析ソフトウェア例えば、汎用非線形構造解析ソフトウェアMarcMSCSoftware社製、構造解析ソフトウェアANSYSANSYS社製、汎用非線形有限要素解析ソフトウェアABAQUSSIMULIA社製、構造解析ソフトウェアNXNastranシーメンス社製等を利用できるが、デファクトスタンダードの観点から、構造解析ソフトウェアNXNastranシーメンス社製によって実行することが好ましい。本発明の推定方法の対象となる繊維強化樹脂射出成形品は、通常使用されるエンジニアリングプラスチック用の熱可塑性樹脂と繊維強化樹脂に使用される繊維を含む樹脂組成物を射出成形して得ることができるが、比較的安価で入手しやすく自動車部品でも多用されている観点から、ポリアミド樹脂及びガラス繊維を含むポリアミド樹脂組成物を射出成形して得られるものが好ましい。繊維強化樹脂射出成形品繊維強化樹脂射出成形品図1を、下記繊維強化樹脂組成物を、下記条件で射出成形して試験片として製造した。1繊維強化樹脂組成物:ガラス繊維強化ナイロンガラス30含有宇部興産株式会社製2射出成形機:品名MD100株ニイガタマシンテク社製3繊維強化樹脂成形品の寸法図2:幅31mm、長さ250mm、厚み2.5mm4射出条件:樹脂温度275°C、金型温度55°C、射出速度80mm/s、保持圧力35MPa、保圧時間10sec、冷却時間20sec5樹脂流動ゲートの種類:1点ゲートと3点ゲートの2種類を製造した。1点ゲートは、図3に示される射出成形型枠の長さ方向の端部で幅方向の中央部を、3点ゲートは、図4に示される射出成形型枠の幅方向の端部で長さ方向の側部の3箇所を射出点とする。ステップI繊維強化樹脂射出成形品について、射出成形解析ソフトウェアAutodeskMoldflowInsights2010を使用して、以下の条件を選択又は入力して下記弾性パラメータ及び繊維配向パラメータを算出した。樹脂流動解析モデルとして、弾性パラメータを算出するためのTandon-Wengマイクロメカニクスモデル、及び、繊維配向パラメータを算出するためのFloger-Tuckerモデルを採用した。繊維強化樹脂射出成形品についての有限要素法に基づく要素分割として、3節点三角形要素を使用し、要素数9342、節点数4793の条件で行った。各要素を厚み方向に20層の分割定義を行い、計算モデルを適用して層ごとの物性値を算出した。各要素に対して、厚み方向に対して均一としたものと力学的に等価となるよう、各層の物性値を重み付け平均して、等価物性値を算出した図6。解析モデルに対して、境界条件となる射出条件として、樹脂温度275°C、金型温度55°C、射出速度80mm/s、保持圧力35MPa、保圧時間10sec、冷却時間20sec樹脂流動ゲートの種類1点ゲート又は3点ゲートを、樹脂流動解析ソフトウェアに入力し、9342の要素毎に、弾性及び繊維配向パラメータとして、流れ方向弾性率:E1、TD方向弾性率:E2、ポアソン比:νせん断弾性率:τ12、τ23、τ31、流れ方向角度:θを算出した。なお、流れ方向角度とは、要素座標のX軸からの流れ方向の角度を表す。図6においては、9342の要素のそれぞれは節点例えば、節点121、234、456で構成されており、節点121→234方向がX軸、それに直交するようにY軸が設定されて、要素座標となる。ステップII上記ステップIで算出された弾性パラメータと繊維配向パラメータとを導入して、構造解析ソフトウェアNXNastranを使用して、繊維強化樹脂射出成形品1点ゲート品及び3点ゲート品の固有振動数を算出した。繊維強化樹脂射出成形品についての有限要素法に基づく要素分割として、6節点三角形要素を使用し、要素数9342、節点数18,927の条件で行った。上記要素分割の設定の下で、表1の各パラメータの値を入力して、上記構造解析ソフトウェアで算出した1点ゲート試験片及び3点ゲート試験片の各モードに対応する固有振動数として表1の値を得た。実験モード解析上記ステップI及びIIを実行して算出された固有振動数を、実験モード解析によって算出される固有振動数と比較した。繊維強化樹脂射出成形品について、以下の条件でハンマリング試験を行った。周波数1つの計測点に対して、加振点HITPOINTA~Cを打撃し周波数伝達関数以下、FRFともいうを計測後、次の計測点の加速度計を移動させる加速度計移動方式を採用した図5。加速度計の信号はFFTアナライザLMS社製Pimentを用いて周波数分析した。樹脂の吸水による弾性率変化の影響を排除するため、試験片は成形品の水分率が0重量となるようにハンマリング試験前に繊維強化樹脂成形品を乾燥した。実験モード解析により得られた各モードに対応する固有振動数として表1の値を得た。表2に示される結果から、本発明の推定方法によれば、繊維強化樹脂射出成形品について、その射出成形条件と、樹脂流動解析から算出される弾性パラメータ及び繊維配向パラメータとを導入した理論モード解析を行うことによって、実験モード解析と同様の精度で繊維強化樹脂射出成形品の固有振動数を算出できることがわかる。"