JP2010069654A - 構造解析方法、構造解析装置、構造解析プログラム、構造解析のための物性値算出方法、構造解析のための物性値算出装置および構造解析のための物性値算出プログラム - Google Patents

Abstract

【課題】繊維状物質を含む材料を用いて射出成形した部材の構造解析について、短い時間で精度良く計算することができる構造解析方法を提供する。
【解決手段】細長い繊維状物質を含む材料を用いて射出成形される部材の構造解析をする方法であって、部材の形状を示す形状情報を取得する形状取得ステップと、部材の形状全体を複数の要素に分割し、分割した各要素を示す要素情報を生成する要素分割ステップ（Ｓ２）と、要素情報が示す各要素を用いて、材料を流し込む金型のゲートに対応する部分に熱源があるとして熱解析を行うことによって部材の熱特性を示す熱特性値を算出して、それを含む熱特性情報を生成する熱解析ステップ（Ｓ３）と、構造解析するための物性値を示す物性情報に熱特性情報を変換する変換ステップ（Ｓ４）と、物性情報と各要素を示す情報とを用いて部材の構造解析をする構造解析ステップ（Ｓ５）とを含む。
【選択図】図１１

Description

本発明は、構造解析方法及び構造解析のための物性値算出方法に関し、特に繊維強化樹脂材料を用いた射出成形品の構造解析方法に関する。

射出成形した樹脂部材は、例えば、テレビ、エアコン、携帯電話の筐体、自動車のバンパ等といった日常の様々な部材に使用されている。樹脂部材の用途が広がるとともに、樹脂に要求される性能や性質も高度になってきている。例えば、携帯電話の軽量化の要請に応えるには、薄くて丈夫な材料が必要になる。このように、耐久性や耐衝撃性が要求される部材には、剛性や強度を向上させるために、樹脂にガラスや炭素等の繊維を添加した繊維強化樹脂がしばしば利用される。

部材や部材の強度を検証するために、有限要素法等を用いた構造解析が行われる。数値シミュレーションは、比較的低コストで種々の部材の性能を検証できるため、部材の低コスト化及び信頼性の向上に貢献する。樹脂部材においても構造解析は一般に行われており、特に繊維強化樹脂では、繊維の配向性を考慮した構造解析が提案されている（例えば、特許文献１及び特許文献２参照）。繊維強化樹脂では、繊維方向の強度と繊維に垂直な方向の強度とでは、前者が後者の３〜５倍程度強い。そのため、繊維の配向性を考慮した構造解析は、射出成形される繊維強化樹脂の強度予測の精度を向上させるために、重要な技術である。

具体的には、特許文献１に記載の構造解析は、大きく３つの工程を含む。まず、流動解析用の要素を用いて、金型を流れる樹脂の流動解析を行う。ここでは、各要素の流れ方向を各要素の繊維配向とすることによって、各要素の繊維配向を示す繊維配向情報を取得する。

次に、流動解析工程で取得した繊維配向情報を、構造解析用の繊維配向情報に変換する。流動解析と構造解析とでは、解析に用いるメッシュが異なる。例えば、流動解析では繊維の挙動を正確にシミュレートするために数百万要素の細かいメッシュが用いられるのに対して、構造解析では数十万要素の粗いメッシュが用いられる。そのため、流動解析用の要素に関連付けられた繊維配向情報を構造解析用の要素に関連付けられた繊維配向情報に変換することによって、繊維配向を考慮した物性値を構造解析に利用できる。

最後に、構造解析用の要素及び構造解析用の繊維配向情報を用いて構造解析を行う。構造解析では、繊維配向情報に基づいて、繊維配向を考慮した物性値が利用される。これにより、高精度の構造解析が可能になる。
特開２００６−２７２９２８号公報 特開２００４−２５７９６号公報

しかしながら、従来の繊維配向を考慮した構造解析では、構造解析全体の計算に多くの時間を要するという問題がある。

すなわち、まず、流動解析に多くの時間が掛かる。解析ではメッシュを細かくする程、計算時間が多くなるところ、流動解析では、上述のように比較的細かいメッシュを用いて計算するからである。また、流動解析用の繊維配向情報を構造解析用の繊維配向情報に変換するために時間が掛かる。流動解析と構造解析とでは解析に用いるメッシュが異なるため、この変換では、流動解析用の要素と構造解析用の要素の対応付け、その後に、対応付けた各要素に適合するように、繊維配向情報を流動解析用のものから構造解析用のものへと変換する必要がある。上述のように流動解析用の要素が数百万要素程度、構造解析用の要素が数十万要素程度であるため、この計算に要する時間は大きい。

本発明は、上述の課題を解決するものであり、繊維状物質を含む材料を用いて射出成形した部材の構造解析について、従来の構造解析の精度を維持又は向上しながら、より短い時間で計算することができる構造解析方法、構造解析装置、構造解析プログラム、構造解析のための物性値算出方法、構造解析のための物性値算出装置および構造解析のための物性値算出プログラムを提供する。

上述の課題を解決するために、本発明に係る構造解析方法は、
繊維状物質を含む材料を用いて射出成形される部材の構造解析方法であって、
上記部材の形状情報を取得する形状取得ステップと、
上記部材を複数の要素に分割し、分割した各要素を示す要素情報を生成する要素分割ステップと、
上記要素情報が示す上記各要素を用いて、上記材料を流し込む金型のゲートに対応する部分に熱源があるとして熱解析を行うことによって上記部材の熱特性を示す熱特性値を算出し、算出した上記熱特性値を含む熱特性情報を生成する熱解析ステップと、
上記熱特性情報を取得し、取得した上記熱特性情報を構造解析に用いる物性値を含む物性情報に変換する変換ステップと、
上記物性情報に含まれる上記物性値と、上記要素情報が示す上記各要素とを用いて上記部材の構造解析を行う構造解析ステップとを含む。

なお、本発明は、構造解析方法として実現されるだけではなく、構造解析方法に含まれるステップを実行する機能を有する処理部を備える構造解析装置として実現することも、構造解析方法に含まれるステップをコンピュータに実行させるプログラムとして実現することもできる。また、本発明は、構造解析のための物性値算出方法として実現されるだけではなく、構造解析のための物性値方法に含まれるステップを実行する機能を有する処理部を備える構造解析のための物性値算出装置として実現することも、構造解析のための物性値算出方法に含まれるステップをコンピュータに実行させるプログラムとして実現することもできる。

本発明では、繊維状物質を含む材料を用いて射出成形した部材の構造解析について、従来の構造解析の精度を維持又は向上しながら、より短い時間で計算することができる。

本発明の実施の形態について図を参照して説明する。

（実施の形態１）
図１は、本発明の実施の形態１に係る構造解析装置の外観を示す。構造解析装置２０ａは、細長い繊維状物質を含む材料を用いて射出成形される部材の構造解析を行う装置であって、例えばソフトウェアプログラムを実行可能なコンピュータである。構造解析装置２０ａには、外部に、液晶ディスプレー等の表示部２２、キーボード、マウス等の入力部２４、ＣＡＤ（Ｃｏｍｐｕｔｅｒ−ａｉｄｅｄ ｄｅｓｉｇｎ，図２参照）装置等が接続されている。本実施の形態では、構造解析を行う部材の一例としてファンを用いて説明する。図１の表示部２２には、ファンを後述するように要素分割した状態が示されている。

なお、部材は、繊維状物質を含む樹脂を用いて作られるものであればよく、製品、半製品、部品等を含み、具体例としては上記ファンの他に、自動車のバンパ、テレビの外枠体、携帯電話の筐体等を挙げることができる。

図２は、実施の形態１に係る構造解析装置２０ａのブロック図である。本図に示すように、構造解析装置２０ａは、ユーザにより入力された情報を入力部２４から取得する入力情報取得部５２と、ハードディスク等の記憶媒体により実現される複数の記憶部とを備える。各記憶部には、詳細は後述するように、各処理部によって生成された情報、ユーザの入力に基づき入力情報取得部５２が取得した情報等が所定の形式で保持されている。

要素記憶部３２は、解析のために分割された各要素の位置関係を特定するための情報である要素情報３２ａを記憶している。図３は、実施の形態１に係る要素情報３２ａの具体例を示す。要素情報３２ａは、各要素を特定するための「要素番号」と、各要素の位置を示す情報とを含む。本実施の形態では、各要素は、「頂点１」から「頂点４」までの４つの頂点を含む四面体である。そのため、各要素の位置を示す情報には、部材に対して固定された座標系での各頂点の座標を示す情報が保持されている。なお、本実施の形態では、各要素が四面体の例を示すが、各要素の形状は四面体に限らず、直方体等が部材の形状に応じて適宜選択されてよい。

熱解析条件記憶部３４は、熱解析のための境界条件を示す情報、熱源の位置を示す情報等の熱解析を行うために必要な条件を示す情報である熱解析条件情報３４ａを記憶している。熱解析条件情報３４ａは、入力情報取得部５２を介してユーザが記憶させる情報である。図４は、実施の形態１に係る熱解析条件情報３４ａの具体例を示す。熱解析条件情報３４ａは、熱解析に用いるための熱伝導率、比熱、ゲート温度、節点温度、ゲート位置等を含む。本発明では、熱解析は部材の熱特性を得るために行われるため、熱伝導率、比熱は部材の材料の正確な値である必要はなく、任意の値、例えば１であってよい。ゲート温度は、射出成形のために樹脂が流し込まれる金型のゲートに対応する部分の温度であり、節点温度は、各要素の境界が交わる節点の温度である。これらの温度も適宜設定される任意の値である。図４では、ゲート温度を１とし、節点温度を０とする例を示している。ゲート位置は、金型のゲートに対応する部分の位置を示す。

熱特性記憶部３６は、構造解析装置２０ａが行う熱解析によって算出される熱特性値を示す熱特性情報３６ａ，３６ｂを記憶している。本実施の形態の熱特性情報３６ａ，３６ｂは、熱流束ベクトル情報３６ａと温度情報３６ｂとを含む。図５は、実施の形態１に係る熱流束ベクトル情報３６ａの具体例を示す。熱流束ベクトル情報３６ａは、熱解析によって算出された各要素の熱流束ベクトルを表す情報である。図６は、実施の形態１に係る温度情報３６ｂの具体例を示す。温度情報３６ｂは、熱解析によって算出された各要素の温度を表す情報である。

変換条件記憶部３８は、熱解析によって算出された熱特性値を、部材の構造解析に用いる物性値に変換するための情報である変換条件情報３８ａを記憶している。変換条件情報３８ａは、入力情報取得部５２を介してユーザが記憶させる情報である。図７は、実施の形態１に係る変換条件情報３８ａの具体例を示す。変換条件情報３８ａは、同図に示されるように、配向性を考慮した物性値に熱特性値を変換するための条件を含む。変換条件情報３８ａは、具体的には、繊維方向の弾性率Ｅ_ＭＤ、繊維方向に垂直な方向の弾性率Ｅ_ＴＤ及び密度ρのそれぞれについて、部材の中の局所的な最大値及び最小値を示す情報を含む。変換条件情報３８ａは、例えば、繊維方向の最大弾性率Ｅ_{ＭＤｍａｘ}を１４００００（ｋｇ／ｃｍ^２）とする。

物性記憶部４０は、構造解析に用いる各要素の物性値を示す情報である物性情報４０ａを記憶している。図８は、実施の形態１に係る物性情報４０ａの具体例を示す。物性情報４０ａは、各要素の熱特性値を変換することによって得られる各要素の物性値を含む。熱特性値及び物性値が算出される要素は、いずれも同一であり、両者は、本実施の形態では要素番号によって関連付けられる。

構造解析条件記憶部４２は、構造解析をするための条件を示す情報である構造解析条件情報４２ａを記憶している。構造解析条件情報４２ａは、入力情報取得部５２を介してユーザが記憶させる。構造解析条件情報４２ａの詳細は、図示しないが、例えば構造解析に用いる境界条件等を含む。

構造解析結果記憶部４４は、構造解析によって算出された結果を示す情報である構造解析結果情報４４ａを記憶している。図９は、実施の形態１に係る構造解析結果情報４４ａの具体例を示す。構造解析結果情報４４ａは、構造解析の結果としての各要素の応力と変位とを示す情報を含む。ここでの要素は、上述の熱特性値及びそれに基づく特性値のそれぞれが算出される要素と同一であり、これらは、要素番号によって関連付けられている。

再び図２を参照すると、構造解析装置２０ａは、形状取得部５６と、要素分割部５８と、熱解析部６０ａと、変換部６４ａと、構造解析部６６と、表示制御部７０とを更に備える。ユーザによる各処理部への指示は、入力情報取得部５２を介して引き渡される。

形状取得部５６は、部材の形状を示す形状情報を取得する。本実施の形態の形状取得部５６は、ＣＡＤ装置２６から接続回線を通じて形状情報を取得する。

要素分割部５８は、形状取得部５６から形状情報を取得すると、部材の全体を複数の要素に分割し、それによって、要素情報３２ａを生成する。

熱解析部６０ａは、熱解析条件記憶部３４から熱解析条件情報３４ａを取得し、要素分割部５８から部材の各要素を示す情報を取得し、これらの取得した情報に基づいて熱解析を行う。また、熱解析部６０ａは、熱解析により得られる解析結果を含む熱特性情報３６ａ，３６ｂを熱特性記憶部３６に記憶させる。

変換部６４ａは、熱解析部６０ａによる熱解析処理が終了すると、変換条件記憶部３８から変換条件情報３８ａを取得し、熱特性記憶部３６から熱特性情報３６ａ，３６ｂを取得し、これらの取得した情報に基づいて、熱特性値を部材の物性値に変換する。

図１０は、実施の形態１に係る変換部６４ａの詳細を示すブロック図である。図示するように、変換部６４ａは、熱特性情報取得部７２ａと、変換条件情報取得部７８と、物性値算出部８０とを有する。

熱特性情報取得部７２ａは、熱特性記憶部３４から熱特性情報３４ａを取得する処理部であり、熱流束ベクトル情報３６ａを取得する熱流束ベクトル取得部７４ａと温度情報３６ｂを取得する温度情報取得部７４ｂとを有する。変換条件情報取得部７８は、変換条件記憶部３８から変換条件情報３８ａを取得する。物性値算出部８０は、熱流束ベクトル情報３６ａと、温度情報３６ｂと、変換条件情報３８ａとに基づいて、各要素の物性値を算出する。物性値算出部８０は、変換して得られる物性値を含む物性情報４０ａを物性記憶部４０に記憶させ、変換処理が終了したことを構造解析部６６に通知する。

構造解析部６６は、変換部６４ａによる変換処理が終了すると、要素記憶部３２から要素情報３２ａを取得し、物性記憶部４０から物性情報４０ａを取得し、構造解析条件記憶部４２から構造解析条件情報４２ａを取得し、これらの取得した情報に基づいて、構造解析を行う。構造解析部６６は、構造解析により得られる解析結果を含む構造解析結果情報４４ａを構造解析結果記憶部４４に記憶させる。また、構造解析部６６は、構造解析結果情報４４ａを表示処理部７０に引き渡す。

表示処理部７０は、構造解析結果情報４４ａに基づいて、構造解析の結果を表示部２２に表示させる。

次に、実施の形態１に係る構造解析装置２０ａが実行する処理の詳細について、フローチャートを参照して説明する。

図１１は、実施の形態１に係る構造解析装置２０ａが実行する処理のフローチャートである。形状取得部５６は、ＣＡＤ装置２６から通信回線を介して形状情報を取得する（Ｓ１）。要素分割部５８は、形状情報に基づいて、部材全体を重複のない複数の要素に分割し、それによって、要素情報３２ａを生成し、生成した要素情報３２ａを要素記憶部３２に記憶させる（Ｓ２）。部材としてのファン９０ａを要素分割した例を図１２に示す。

熱解析部６０ａは熱解析処理を実行する（Ｓ３）。図１３は、熱解析処理（Ｓ３）の詳細を示すフローチャートである。熱解析部６０ａは、要素分割部５８から要素情報３２ａを取得する（Ｓ１１）。熱解析部６０ａは、熱解析条件記憶部３４から熱解析条件情報３４ａを取得する（Ｓ１２）。熱解析部６０ａは、各取得処理（Ｓ１１及びＳ１２）において取得した情報を用いて、例えば有限要素法による熱解析を行い、それによって、各要素の熱流束ベクトル及び温度を含む熱特性値を算出する（Ｓ１３）。ここで、有限要素法は、熱解析及び後述する構造解析等に利用される数値解析手法の１つであって、物体を複数の小さい要素に分割し、各要素に単純化した方程式を組み合わせて適用し、適用した方程式を満足する解を数値的に得る手法である。有限要素法は、よく知られた手法であるため、これに関する詳細な説明は省略する。

熱解析部６０ａは、熱特性値の算出処理（Ｓ１３）において算出した熱特性値を含む熱特性値情報（熱流束ベクトル情報３６ａ，温度情報３６ｂ）を熱特性記憶部３６に格納して記憶させる（Ｓ１４）。

ここで、熱解析部６０ａが熱解析に用いる熱伝導の支配方程式を式（１）に示す。

・・・式（１）

式（１）において、Κは熱伝導率、Ｔは温度，（上に点を３つ付した）ｑは単位体積あたりの内部発熱量、ρは密度、ｃは比熱、ｔは時間、Ｖ_ｘ，Ｖ_ｙ，Ｖ_ｚは伝導媒体中の熱の速度をそれぞれ表す。

これに対して、流動解析に用いられる支配方程式を式（２）に示す。

・・・式（２）

式（２）において、Ｖ_ｘ，Ｖ_ｙ，Ｖ_ｚは樹脂の流速、Ｐは圧力、τは表面力をそれぞれ表す。

式（２）の右辺第１項〜第５項は、それぞれ、熱伝導、圧縮による内部エネルギー、法線応力による仕事、せん断応力による仕事、内部発熱を表す。また、射出成形ＣＡＥ（Ｃｏｍｐｕｔｅｒ Ａｉｄｅｄ Ｅｎｇｉｎｅｅｒｉｎｇ）において、式（２）は、以下に示す（ａ）〜（ｇ）の仮定により、式（３）のように簡略化できる。
（ａ）流動塲は薄肉である。すなわち、肉厚に比べて流動長が十分に長く、そのため、肉厚方向の流動成分Ｖ_ｚは無視できる。
（ｂ）樹脂は高粘性流体であり、そのため、慣性項及び体積力項は表面力に比べて小さく無視できる。
（ｃ）平面方向の速度勾配が肉厚方向の速度勾配に比べて十分に小さく、そのため、面内の表面力は無視できる。
（ｄ）モデル形状が薄肉であるため、面方向の熱伝導率は肉厚方向の熱伝導に比べて十分に小さく、そのため、面方向の熱伝導率は無視できる。
（ｅ）圧縮による内部エネルギーは、速度勾配が小さいので無視できる。
（ｆ）法線応力による仕事は、断面急変部のみで内部発熱として考慮する。
（ｇ）剪断応力による仕事は、肉厚方向の速度勾配成分のみを対象とする。

・・・式（３）

式（３）の右辺は、熱伝導と、肉厚方向の速度勾配成分のせん断応力による仕事と、内部発熱の項とに簡略化される。これを式（１）と比較すると、熱伝導の支配方程式は、簡略化した流動解析の支配方程式において、肉厚方向の速度勾配成分のせん断応力の影響を無視したものである。従って、せん断応力を無視できる場合、すなわち、せん断応力を決定する粘度およびせん断速度を無視できる場合に、流動解析を伝熱解析で簡易的に近似できる。ここで、粘度は、せん断速度、温度及び圧力の関数で表されるため、これらが大きく変化しない単純な形状を部材が有している場合には、材料に含まれる繊維状物質の配向性、すなわち、金型に流し込まれた繊維の方向を解析するために熱伝導の支配方程式を利用することは有用である。

せん断速度、温度及び圧力が大きく変化しない形状の具体例としては、ゲート数が少ない形状、板厚の増減が少ない形状等を挙げることができる。また、流動過渡部（流動先端、金型から遠い成形品内部等）も、せん断速度、温度及び圧力が大きく変化しない部位であり、当該部分についても熱解析を用いても十分な精度で、繊維の配向性を取得することができる。この場合、熱流束ベクトルが繊維の配向性を示す値になる。

このように、熱解析を行うことによって、部材に含まれる樹脂の配向性を求める。式（２）で表される支配方程式の複雑さから分かるように、流動解析の計算は複雑であり、計算時間が多く掛かる。これに対して、熱解析の計算負荷は、流動解析の場合に比べて軽く、そのため、樹脂の配向性を算出するために要する時間を短縮することが可能になる。

図１４は、ファン９０ａの熱解析によって得られた各要素の熱流束ベクトルの一部を示す。図１４に示す各矢印の方向及び長さが、それぞれ、各要素の熱流束ベクトルの方向及び大きさを示す。成形時に金型のゲートが位置する部分であるゲート部分９２ａに樹脂が流し込まれると、樹脂がゲート部分９２ａから羽根の遠方へ向かって放射状に広がり、羽根の輪郭に沿って流れる様子が示されている。樹脂に含まれる繊維の向きは樹脂の流れ方向に沿うため、図１４の矢印は繊維の配向性を示す。

図１５は、ファン９０ａの正面を示す写真であり、図１６は、図１５に示される矩形部分９４ｂのＸ線ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）装置による画像を拡大して示しており、図１４の矩形部分９４ａに対応する部分である。図１６に示す画像において、黒色の線状の物質が樹脂に含まれる繊維である。図１６から分かるように、繊維は概ね右上から左下を向いている。

図１７は、Ｘ線ＣＴ装置による画像に基づいて作成したファン９０ａの繊維の配向性を示す図である。図１４の全体と比較すると、熱解析によって得られた繊維配向性、すなわち、図１４に示す矢印の方向と、Ｘ線ＣＴ装置による画像に基づく実際の繊維の配向性とは、ほぼ一致していることがわかる。従って、熱解析によって、樹脂に含まれる繊維の配向性を精度よくシミュレートすることができる。

再び図１１を参照する。変換部６４は、熱解析で得られた熱特性値を構造解析に用いる物性値に変換する（Ｓ４）。図１８は、実施の形態１の変換処理（Ｓ４）の詳細を示すフローチャートである。

熱特性情報取得部７２ａは、熱特性情報３６ａ，３６ｂを取得する（Ｓ２１）。より詳細には、熱流束ベクトル取得部７４ａが熱流束ベクトル情報３６ａを、温度情報取得部７４ｂが温度情報３６ｂをそれぞれ取得する。変換条件情報取得部７８は、変換条件記憶部３８から変換条件情報３８ａを取得する（Ｓ２２）。物性値算出部８０は、熱流束ベクトル情報３６ａと、温度情報３６ｂと、変換条件情報３８ａとに基づいて、物性値を算出する（Ｓ２３）。

ここで、物性値算出部８０による物性値算出処理（Ｓ２３）について詳しく説明する。

まず、物性値算出部８０は、各要素の温度を各要素の弾性率の大きさに変換する。この変換では、以下に示す式（４）を用いる。式（４）では、変換対象となる要素の温度をＴ、最高温度及び最低温度をそれぞれＴ_ｍａｘ及びＴ_ｍｉｎ、要素の繊維方向の弾性率をＥ_ＭＤ、繊維方向の最大弾性率及び最小弾性率をそれぞれＥ_{ＭＤｍａｘ}及びＥ_{ＭＤｍｉｎ}とする。

・・・式（４）

なお、繊維と垂直な方向の弾性率Ｅ_ＴＤは、式（４）のＥ_ＭＤをＥ_ＴＤに、Ｅ_{ＭＤｍａｘ}をＥ_{ＴＤｍａｘ}に、Ｅ_{ＭＤｍｉｎ}をＥ_{ＴＤｍｉｎ}にそれぞれ置き換えた式によって変換する。

このように、物性値算出部８０は、温度情報３６ｂに含まれる最高温度を変換条件情報３８ａに含まれる弾性率の最大値（例えば、Ｅ_{ＭＤｍａｘ}）に対応付け、温度情報３６ｂに含まれる最低温度を変換条件情報３８ａに含まれる弾性率の最小値（例えば、Ｅ_{ＭＤｍｉｎ}）に対応付け、温度情報３６ｂに含まれる温度が高い要素である程、当該要素の上記弾性率が大きくなるように、温度情報３６ｂに含まれる各要素の温度を各要素の弾性率に一対一で対応付ける。本実施の形態では、上記最高温度と上記最低温度とのそれぞれに重み付けをして各要素の温度が得られる場合に、上記弾性率の最大値と上記弾性率の最小値とに上記各要素の温度を得た場合と同じ重み付けをして得られる値を前記各要素の弾性率として算出する。

一般に、１つの部材においてゲートの近傍は弾性率が大きく、ゲートから離れるに従って弾性率が小さくなる。繊維を含む樹脂の弾性率の違いは、主に樹脂に含まれる繊維の長さの違いによって生じる。具体的には、繊維の長さが長い程、樹脂の弾性率は大きくなる。樹脂を金型に流す場合、長い繊維は、流れ難いためにゲート近傍に多くとどまる。これに対して、短い繊維は、流れ易いためにゲートから離れた部分にまで流れる。そのため、上述のような、ゲート近傍では弾性率が大きく、そこから離れる程小さくなるという弾性率の傾向が生じる。

他方、ゲートを熱源とした場合の部材の温度分布も、繊維の流れ易さと同様に、ゲートの近傍は温度が高く、ゲートから離れるに従って温度が低くなる。そのため、繊維の流動性の違いによって生じる弾性率の分布と、熱伝達の結果である温度分布とは、類似した傾向を示す。そして、各要素の温度と弾性率とは、上記式（４）のような簡単な変換によってほぼ対応付けられる。

次に、物性値算出処理（Ｓ２３）において物性値算出部８０は、各要素の熱流束ベクトルに基づいて弾性率の配向性を算出する。各要素の繊維方向が熱流束ベクトルの方向と一致するように、全体の座標系を各要素の座標系である要素座標系に変換するための係数（マトリクス）を算出する。この係数（マトリクス）が各要素の配向性を示す値（マトリクス）となる。

最後に、物性値算出部８０は、各要素の温度を各要素の密度に変換する。この変換では、以下に示す式（５）を用いる。式（５）では、要素の温度をＴ、温度情報３６ｂに含まれる最高温度及び最低温度をそれぞれＴ_ｍａｘ及びＴ_ｍｉｎ、要素の密度をρ、変換条件情報３８ａの最大密度及び最小密度をそれぞれρ_ｍａｘ及びρ_ｍｉｎとする。

・・・式（５）

このように、物性値算出部８０は、温度情報３６ｂに含まれる最高温度を変換条件情報３８ａに含まれる密度の最大値ρ_ｍａｘに対応付け、温度情報３６ｂに含まれる最低温度を変換条件情報３８ａに含まれる密度の最小値ρ_ｍｉｎに対応付け、温度情報３６ｂに含まれる温度が高い要素である程、当該要素の上記密度が大きくなるように、温度情報３６ｂに含まれる各要素の温度を各要素の弾性率に一対一で対応付ける。本実施の形態では、上記最高温度と上記最低温度とのそれぞれに重み付けをして各要素の温度が得られる場合に、上記密度の最大値と上記密度の最小値とに上記各要素の温度を得た場合と同じ重み付けをして得られる値を前記各要素の弾性率として算出する。

一般に、均一な厚さを有する板のように単純な形状を有する部材の場合、金型のゲートからの距離が遠ざかる程、樹脂は流れにくくなる。金型のゲート付近では樹脂の温度は高いが、ゲートから離れるに従って樹脂の温度は低下し、それに伴って、樹脂の粘度が高くなるためである。このような金型を流れる樹脂の粘度が変化するために、部材の密度は、ゲート部分で大きく、ゲート部分離れるに従って小さくなる。従って、部材の密度分布は、ゲート部分を熱源とした場合の部材の温度分布と同様になる。そして、各要素の温度と密度とは、上記式（５）のような簡単な変換によってほぼ対応付けられる。

図１８の変換処理の説明に戻って、最後に、物性値算出部８０は、算出した物性値を示す物性情報４０ａを物性記憶部４０に格納する（Ｓ２４）。ここでの物性値は、上述のような繊維方向及びそれに垂直な方向の弾性値、配向性、密度等を含む（図８には、要素番号「１０００」である要素に関する繊維方向及びそれに垂直な方向の弾性値と、密度とを例示する。）。繊維方向及びそれに垂直な方向の弾性値、並びに、配向性を示す値をまとめて、以下、配向性弾性値という。

このように、変換処理において、分割要素の違いを考慮しておらず、すなわち、分割要素が異なるために生じる要素のマッピング、各要素の値の統合や分配等を行わない。後述するように、本発明では熱解析と構造解析とで同一の分割要素を用いることができるためである。従って、変換処理を簡単にすることができ、処理時間を短くすることが可能になる。

再び図１１を参照する。構造解析部６６は、変換処理（Ｓ４）において算出された物性値を用いて部材の構造解析を行う（Ｓ５）。図１９は、構造解析処理（Ｓ５）の詳細を示すフローチャートである。

構造解析部６６は、要素記憶部３２から要素情報３２ａを取得する（Ｓ３１）。ここで取得する要素情報３２ａは、要素分割処理（Ｓ２）において要素分割部５８により生成された要素情報である。すなわち、熱解析処理（Ｓ３）に利用された要素情報と同一である。

構造解析部６６は、物性記憶部４０から物性情報４０ａを取得する（Ｓ３２）。ここで取得する物性情報４０ａは、変換処理（Ｓ４）において算出された物性値であり、すなわち、本実施の形態では配向性弾性値及び密度を含む。

構造解析部６６は、構造解析条件記憶部４２から構造解析条件情報４２ａを取得する（Ｓ３３）。構造解析条件情報４２ａは、構造解析に必要な境界条件等の条件を含む情報であり、詳しくは、分割された各要素又は各要素の節点に、例えば変位及び回転の拘束、荷重、応力等の解析条件を想定する状況に応じてユーザが適宜設定する条件を含む情報である。

ファン９０ａの場合、構造解析条件情報４２ａは、例えば、ファン９０ａの中心に挿入される回転軸と嵌り合う部分を固定端とし、その他の表面を自由端とする情報を含む。また、この場合の構造解析条件情報４２ａは、回転速度に応じた遠心力及びファンが受ける空気圧を示す情報を含んでもよい。

構造解析部６６は、熱解析処理（Ｓ３）と同様に、要素情報３２ａに含まれる要素を用いて、例えば有限要素法等の手法によって構造解析を行う。それによって、構造解析部６６は、構造解析結果を算出し、算出した構造解析結果を含む構造解析結果情報４４ａを生成する（Ｓ３４）。本実施の形態では、構造解析結果は部材の応力及び変位を含む。構造解析結果によって、例えば、最大回転速度に耐えられるファン９０の羽根の厚さ等を繊維の配向性を考慮した構造解析に基づいて検討できる。従って、

このように、構造解析と熱解析とで、同一の分割要素を用いる。そのため、異なる要素間でのマッピング、すなわち、熱解析の結果に基づいて算出された各要素の物性値を異なる要素に変換する必要がなくなる。そのため、構造解析処理を簡単にすることができ、当該処理に要する計算時間を短縮することが可能になる。また、構造解析には、繊維の配向性を考慮した物性情報が用いられる。そのため、高精度の構造解析が可能になり、部材及び射出成形のための金型の緻密で高度な設計が可能になる。

図１９の構造解析処理の説明に戻って、最後に、構造解析部６６は、構造解析結果情報４４ａを構造解析結果記憶部４４に格納する（Ｓ３５）。これにより、ユーザは構造解析の結果が保持され、それを随時参照できる。

再び図１１に戻る。表示制御部７０は、構造解析部６６から構造解析結果情報４４ａを取得し、構造解析結果情報４４ａを示す画像を表示部２２に表示する。これによって、構造解析装置２０ａは処理を終了する。

以上のように、繊維特性情報の算出工程及び繊維特性情報を物性情報に変換する工程で計算負荷を軽減できる。そのため、各工程に要する処理時間が短縮される。また、繊維の配向性を考慮した物性情報に基づいて構造解析を行うため、高精度の構造解析の結果を得ることができる。従って、高精度な構造解析を、従来の技術によって計算する場合よりも短時間で行うことが可能になる。

（実施の形態２）
実施の形態２では、熱解析によって算出される熱特性値を変換することによって構造特性値を算出する。ここで、構造特性値とは、部材の構造上の特性を示す値であって、本実施の形態では部材の各要素の強度である。

また、実施の形態２では、ウェルドが発生する部材において、ウェルドが発生する位置（ウェルド位置）を考慮して強度を算出する。

ウェルドとは、金型に流し込まれた材料の流れが合流して融着する部分であり、線状又は帯状に表れることが多いが形状はこれに限られない。ウェルドは、複数のゲートから材料が流し込まれる場合に、各ゲートから流入する材料が衝突する位置で生じる。また例えば、部材の中に窓や孔が設けられる場合に、ウェルドはその周辺で生じる。ウェルドは部材の美観を損ない、またウェルド位置では強度が低下する。そのため、ウェルドが部材の不良の原因となることがある。

ウェルドによる不良を防止するために、成形用の金型や部材には種々の工夫がなされる。例えば、部材の強度に与える影響が少ない位置にウェルドが生じるように、金型のゲートが設けられ、また金型及び部材の形状が設計される。また、ウェルドの発生を抑えられるように、流し込まれる材料の温度は調整される。このように、ウェルドが生じる部材では、ウェルドによる不良を防ぐ前提として、ウェルド位置の予測が重要である。

以下、本発明の実施の形態２について図を参照して説明する。本実施の形態では、全図を通して、実施の形態１と同一の部位及び処理部には同じ参照符号を付している。

図２０は、本発明の実施の形態２に係る構造解析装置の外観を示す。構造解析装置２０ｂは、細長い繊維状物質を含む材料を用いて射出成形される部材の強度を算出する装置である。本実施の形態では、図２０に示す表示部２２に表示されている板９０ｂを部材の例として説明する。板９０ｂは、図示するように、金型のゲートに対応する２つのゲート部分９２ｂ，９２ｃを有する。部材が複数のゲート部分を有するため、板９０ｂにはウェルドが発生する。

図２１は、実施の形態２に係る構造解析装置２０ｂのブロック図である。以下、本実施の形態に特徴的な情報を記憶している記憶部及び本実施の形態に特徴的な処理を行う処理部について説明し、実施の形態１と同様の情報を記憶している記憶部及び処理を行う処理部に関するここでの説明は省略する。

要素記憶部３２は、実施の形態１と同様に要素情報３２ｂを記憶している。要素情報３２ｂは、図２２に示すように、各要素のｘ方向、ｙ方向及びｚ方向の大きさと、各要素を特定するための「要素番号」に対応付けられた各要素の基準位置とを含む。本実施の形態の各要素は、図２１の表示部２２に示すように、立方体であるため、要素情報３２ｂによって各要素の位置関係が特定できる。本実施の形態では、ｘ方向の大きさ（Δｘ）、ｙ方向の大きさ（Δｙ）及びｚ方向の大きさ（Δｚ）を全て「０．５」とし、また、例えば要素番号が「１」である要素の基準位置を「（０，０，０）」としている。なお、この要素は、実施の形態１と同様に、熱解析等の計算に用いられるが、要素の形状は、四面体、立方体等が適宜選択的に使用されてよい。

熱解析条件記憶部３４は、熱解析条件情報３４ｂを記憶している。本実施の形態の熱解析条件情報３４ｂは、図２３に示すように、実施の形態１に係る熱解析条件情報３４ａと異なり、２つのゲート位置を含む。

熱特性記憶部３６は、熱特性情報３６ｃ〜３６ｅを記憶している。図２４〜図２６は、それぞれ熱特性情報のうち、温度勾配情報３６ｃ、熱流束ベクトル情報３６ｄ及び温度情報３６ｅの具体例を示す。温度勾配情報３６ｃは、各要素における温度勾配を示す情報を含む。熱流束ベクトル情報３６ｄ及び温度情報３６ｅはそれぞれ、実施の形態１と同様に、各要素の熱流束ベクトル及び温度を示す情報である。

変換条件記憶部３８は、変換条件情報３８ｂを記憶している。変換条件情報３８ｂは、入力情報取得部５２を介してユーザが記憶させる情報である。図２７は、実施の形態２に係る変換条件情報３８ｂの具体例を示す。変換条件情報３８ｂは、同図に示されるように、熱特性値を強度に変換するための条件を含む。強度には、引張り強さ、降伏点強度、引張強さの７０％か降伏点強度のうち小さい方等が用いられてよい。繊維を含む樹脂の強度は、繊維の方向とそれに垂直な方向とでは大きく異なる。そのため、繊維を含む樹脂の強度は配向性を有する。以下、繊維に沿った方向の強度をσ_ＭＤ、繊維に垂直な方向の強度をσ_ＴＤとする。

変換条件情報３８ｂは、具体的には、各要素の温度を各要素の強度（繊維方向の強度σ_ＭＤ及び繊維方向に垂直な方向の強度σ_ＴＤ）を得るために、部材中の強度の最大値及び最小値（最大強度σ_{ＭＤｍａｘ}、最小強度σ_{ＭＤｍｉｎ}）を示す情報を含む。図２７に示す変換条件情報３８ｂは、繊維方向の最大強度σ_{ＭＤｍａｘ}を１８００（ｋｇ／ｃｍ^２）とする例を示す。また、変換条件情報３８ｂは、ウェルド位置での強度を算出するための係数であるウェルド強度倍率も含む。

ウェルド強度倍率は、材料自体の強度に対するウェルド部分での強度が倍率であり、すなわち、材料でウェルドが発生した部分の強度が材料自体の強度に比べてどの程度であるかを示す。図２７に示す変換条件情報３８ｂは、ウェルド強度倍率を１／３とする例を示す。

再び図２１を参照する。本実施の形態の熱解析部６０ｂは、熱解析条件記憶部３４から熱解析条件情報３４ｂを取得し、要素分割部５８から部材の各要素を示す情報を取得し、これらの取得した情報に基づいて熱解析を行い、これによって、温度勾配情報３６ｃと熱流束ベクトル情報３６ｄと温度情報３６ｅとを含む熱特性情報３６ｃ〜３６ｅを算出する。

変換部６４ｂは、熱解析部６０ｂによる熱解析処理が終了すると、変換条件記憶部３８から変換条件情報３８ｂを取得し、熱特性記憶部３６から熱特性情報３６ｃ〜３６ｅを取得し、これらの取得した情報に基づいて、熱特性値を強度に変換する。

図２８は、実施の形態２に係る変換部６４ｂの詳細を示すブロック図である。図示するように、変換部６４ｂは、熱特性情報取得部７２ｂと、変換条件情報取得部７８と、ウェルド位置特定部８２と、強度算出部８４ａとを有する。

熱特性情報取得部７２ｂは、熱特性記憶部３６から熱流束ベクトル情報３６ｄを取得する熱流束ベクトル情報取得部７４ａと、熱特性記憶部３６から温度情報３６ｅを取得する温度情報取得部７４ｂと、熱特性記憶部３６から温度勾配情報３６ｃを取得する温度勾配情報取得部７４ｃとを有する。変換条件情報取得部７８は、実施の形態１と同様に、変換条件記憶部３８から変換条件情報３８ｂを取得する。

ウェルド位置特定部８２は、温度勾配情報取得部７４ｃを介して温度勾配情報７４ｃを取得し、取得した温度勾配情報７４ｃに基づいて部材のウェルド位置を特定する。図２９は、２つのゲート部分を熱源として熱解析によって算出された温度分布を示す図である。２９を参照すると、部材９０ｂにおけるウェルド位置１００は、ゲート９２ｂ，９２ｃから流入した樹脂等の材料の流れが合流する場所であり、熱の伝わりに置き換えて考えると、ゲート９２ｂ，９２ｃを熱源として隣接する要素間で温度が等しい場所である。そのため、温度勾配が閾値以下になる要素が、ウェルド位置１００であるとして特定される。

ウェルド位置特定部８２は、特定したウェルド位置の要素を示す情報をウェルド位置情報として生成し、強度算出部８４ａに引き渡す。

なお、ウェルド位置を特定するための温度勾配の閾値は、「０」であってもよい。この場合、温度勾配が０である要素が、ウェルド位置として特定される。

なお、温度勾配ではなく、例えば温度情報を用いてもよい。例えば図２９の等しい温度を結んだ等温線９６ａ〜９６ｅから分かるように、ゲートから金型に流入した材料の温度は、ゲート部分９２ｂ及び９２ｃから離れるに従って次第に低くなる。そのため、ウエルド位置は、温度が低い位置に対応する。従って、ウエルド位置は、閾値よりも温度の低い位置として特定することもできる。

強度算出部８４ａは、変換条件情報３８ｂと、温度情報３６ｅと、熱流束ベクトル情報３６ｄと、ウェルド位置情報とを取得し、これらの取得した情報に基づいて強度を算出する。

次に、実施の形態２に係る構造解析装置２０ｂが実行する処理の詳細について、フローチャートを参照して説明する。

図３０は、実施の形態２に係る構造解析装置２０ｂが実行する処理のフローチャートである。形状取得部５６及び要素分割部５８は、それぞれ、実施の形態１の形状情報取得処理（Ｓ１）及び要素分割処理（Ｓ２）と同様に、形状情報取得処理（Ｓ４１）及び要素分割処理（Ｓ４２）を実行する。本実施の形態に係る板９０ｂの要素分割の例を図２０に示す。

熱解析部６０ｂは熱解析処理を実行する（Ｓ４３）。図３１は、実施の形態２に係る熱解析処理（Ｓ４３）の詳細を示すフローチャートである。熱解析部６０ｂは、要素分割部５８から要素情報３２ｂを取得する（Ｓ５１）。熱解析部６０ａは、熱解析条件記憶部３４から熱解析条件情報３４ａを取得する（Ｓ５２）。熱解析部６０ａは、各取得処理（Ｓ５１及びＳ５２）において取得した情報を用いて、例えば有限要素法による熱解析を行い、それによって、各要素の熱流束ベクトル、温度勾配及び温度を含む熱特性値を算出する（Ｓ５３）。熱解析部６０ａは、熱特性値の算出処理（Ｓ５３）において算出した熱特性値を含む熱特性情報３６ｃ〜３６ｅを熱特性記憶部３６に格納する（Ｓ１４）。

再び図３０を参照して、変換部６４ｂは、熱特性値を強度に変換する（Ｓ４４）。図３２は、実施の形態２の変換処理（Ｓ４４）の詳細を示すフローチャートである。

熱特性情報取得部７２ｂは、熱特性情報３６ｃ〜３６ｅを取得する（Ｓ６１）。より詳細には、熱流束ベクトル取得部７４ａが熱流束ベクトル情報３６ｄを、温度情報取得部７４ｂが温度情報３６ｅを、温度勾配情報取得部７４ｃが温度勾配情報３６ｃをそれぞれ取得する。

ウェルド位置特定部８２は、温度勾配情報３６ｃに基づいてウェルド位置を特定する（Ｓ６２）。ここで、ウェルド位置特定部８２は、上述のように、各要素の温度勾配が閾値以下であるか否かを判断し、閾値以下である要素をウェルド位置と特定する。

なお、上述のように、ウェルド位置は、閾値よりも温度が低い位置として特定することもできる。この場合、ウェルド位置特定部８２は、温度勾配情報３６ｃに基づいてウェルド位置を特定する処理（Ｓ６２）に代えて、温度分布（各要素の温度）からウエルド位置を特定する。すなわち、ウェルド位置特定部８２は、各要素の温度を含む温度情報３６ｅを参照して、板９０ｂの中の最低温度を特定するとともに、各要素の温度が最低温度と等しいか又は閾値以下であるか否かを判断し、最低温度と等しい又は閾値以下である要素をウェルド位置として特定する。ここで、板９０ｂの中の最低温度とは、板９０ｂに含まれる各要素の温度の中で最も低い温度をいう。

変換条件情報取得部７８は、変換条件記憶部３８から変換条件情報３８ｂを取得する（Ｓ６３）。強度算出部８４ａは、熱流束ベクトル情報３６ａと、温度情報３６ｂと、変換条件情報３８ａとに基づいて強度を算出する（Ｓ６４）。

ここで、強度算出部８４ａによる強度算出処理（Ｓ６４）について詳しく説明する。

強度算出部８４ａは、各要素の温度を各要素の強度に変換する。ここで、強度には、繊維方向の強度の大きさ、繊維と垂直な方向の強度の大きさ、及び、繊維方向に従って強度の配向性を示す値を含む。

強度算出処理（Ｓ６４）では、強度の各方向の大きさを算出するために、以下に示す式（６）を用いる。式（６）は、式（４）の弾性値Ｅを強度σに置き換えたものであり、要素の温度をＴ、最高温度及び最低温度をそれぞれＴ_ｍａｘ及びＴ_ｍｉｎ、要素の繊維方向の強度をσ_ＭＤ、繊維方向の最大強度及び最小強度のそれぞれの大きさをσ_{ＭＤｍａｘ}及びσ_{ＭＤｍｉｎ}とする。

・・・式（６）

なお、繊維と垂直な方向の弾性率σ_ＴＤは、式（６）のσ_ＭＤをσ_ＴＤに、σ_{ＭＤｍａｘ}をσ_{ＴＤｍａｘ}に、σ_{ＭＤｍｉｎ}をσ_{ＴＤｍｉｎ}にそれぞれ置き換えた式によって変換する。

このように、強度算出部８４ａは、温度情報３６ｂに含まれる最高温度を変換条件情報３８ｂに含まれる最大強度（例えば、繊維方向の最大強度σ_{ＭＤｍａｘ}）に対応付け、温度情報３６ｂに含まれる最低温度を変換条件情報３８ａに含まれる最小強度（例えば、繊維方向の最小強度σ_{ＭＤｍｉｎ}）に対応付け、温度情報３６ｂに含まれる温度が高い要素である程、当該要素の上記強度が大きくなるように、温度情報３６ｂに含まれる各要素の温度を各要素の強度に一対一で対応付ける。本実施の形態では、上記最高温度と上記最低温度とのそれぞれに重み付けをして各要素の強度が得られる場合に、上記最大強度と上記最小強度とに上記各要素の温度を得た場合と同じ重み付けをして得られる値を前記各要素の強度として算出する。

上述のように、１つの部材においてゲートの近傍は長い長さを有する繊維を多く含み、そのため樹脂の強度は強い。そして、ゲートから離れるに従って、含まれる繊維の長さが短くなり、そのため強度は弱くなる。また、ゲートを熱源とした場合の部材の温度分布は、ゲート近傍程温度が高く、ゲートから離れる程温度が低くなる。そのため、繊維の流動性の違いによって生じる樹脂の強度の大きさの分布と、熱伝達の結果である温度分布とは、類似した傾向を有しており、上記式（６）のような簡易な変換によってほぼ対応付けられる。

次に、強度算出処理（Ｓ６４）において強度算出部８４ａは、各要素の熱流束ベクトルに基づいて強度の配向性を算出する。各要素の繊維方向が熱流束ベクトルの方向と一致するように、全体の座標系を各要素の座標系である要素座標系に変換するための係数（マトリクス）を算出する。この係数（マトリクス）が各要素の強度の配向性を示す値（マトリクス）となる。

強度算出部８４ａは、ウェルド位置の強度を算出する（Ｓ６５）。強度算出部８４ａは、ウェルド位置特定処理（Ｓ６２）において特定されたウェルド位置に対応する要素について、強度算出処理（Ｓ６４）において算出された繊維方向及びそれに垂直な方向の強度のそれぞれの大きさにウェルド強度倍率を掛けることによって、ウェルド位置に対応する要素の各方向の強度の大きさを修正する。

強度算出部８４ａは、ウェルド位置に対応する要素については、ウェルド位置の強度算出処理（Ｓ６５）において修正された強度を含み、又その他の要素については、強度算出処理（Ｓ６４）において算出された強度を含む強度情報を構造解析結果情報４４ｂとして構造解析結果記憶部４４に格納する（Ｓ６６）。

以上のように、熱解析によって得られる温度分布を変換することによって、温度分布から直接的に部材の強度を算出する。しかも、変換は上記式（６）のように簡単な式に従って処理され、そのため、変換に要する処理時間は短い。従って、部材の強度を十分な精度でかつ非常に短い時間で行うことが可能になる。

（実施の形態３）
本実施の形態では、実施の形態１と同様の方法によって応力分布を算出し、また実施の形態２と同様の方法によって強度分布を算出し、これらの算出された応力分布と強度分布をと用いて部材の安全率分布を作成する。

以下、本発明の実施の形態３について図を参照して説明する。本実施の形態では、全図を通して、実施の形態１及び実施の形態２と同一の部位及び処理部には同じ参照符号を付している。

実施の形態３に係る構造解析装置２０ｃの外観は、図１に示す実施の形態１に係る構造解析装置２０ａの外観と同様であり、実施の形態３に係る構造解析装置２０ｃは、例えばコンピュータ及びそれに接続されたモニタ等の表示部２２とキーボード、マウス等の入力部２４により実現される。

図３３は、本発明の実施の形態３に係る構造解析装置２０ｃのブロック図である。

変換条件記憶部３８は、図３４に示すように、弾性率Ｅ、密度ρ及び強度σを算出するための条件を含む変換条件情報３８ｃを記憶している。弾性率Ｅを算出するための変換条件として変換条件情報３８ｃは、繊維方向及びそれに垂直な方向の最大弾性率Ｅ_{ＭＤｍａｘ，}Ｅ_{ＴＤｍａｘ}及び最小弾性率Ｅ_{ＭＤｍｉｎ，}Ｅ_{ＴＤｍｉｎ}を含む。密度ρを算出するための変換条件として変換条件情報３８ｃは、最大密度ρ_ｍａｘと最小密度ρ_ｍｉｎを含む。強度σを算出するための変換条件として変換条件情報３８ｃは、繊維方向及びそれに垂直な方向の最大弾性率σ_{ＭＤｍａｘ，}σ_{ＴＤｍａｘ}及び最小弾性率σ_{ＭＤｍｉｎ，}σ_{ＴＤｍｉｎ}を含む。

構造解析結果記憶部４４は、構造解析の結果としての各要素の応力を示す情報と、各要素の強度を示す情報とを含む構造解析結果情報４４ｃを記憶している。本実施の形態では、ウェルド位置を考慮しない強度を例に説明するが、実施の形態２と同様にウェルド位置を考慮した強度を算出し、それによって、ウェルド位置を考慮した強度を用いて、後述する安全率を算出してもよい。また、構造解析結果情報４４ｃには、変位を示す情報が含まれてもよい。

また、構造解析結果記憶部４４は、応力及び強度を用いて算出された各要素の安全率を示す情報を更に含む構造解析結果情報４４ｃを記憶している。

変換部６４ｃは、熱解析部６０ａによる熱解析処理が終了すると、変換条件記憶部３８から変換条件情報３８ａを取得し、熱特性記憶部３６から熱特性情報３６ａ，３６ｂを取得し、これらの取得した情報に基づいて、熱特性値を物性値及び強度に変換する。

図３５は、実施の形態３に係る変換部６４ｃの詳細を示すブロック図である。変換部６４ｃは、図示するように、実施の形態１に係る変換部６４ａが有する各処理部に加えて、強度算出部８４ｂを有する。

強度算出部８４ｂは、熱特性情報取得部７２ａから熱流束ベクトル情報と温度情報と取得し、変換条件情報取得部７８から変換条件情報の強度σに関する情報を取得し、取得したこれらの情報に基づいて部材の各要素の強度を算出する。

再び図３３を参照する。安全率算出部８６は、変換部６４ａで算出された各要素の強度と構造解析部により算出された各要素の応力とを構造解析部６６から取得し、取得した各要素の強度と各要素の応力とを用いて、対応する要素の応力に対する強度の倍率（強度／応力）を各要素について算出する。

次に、実施の形態３に係る構造解析装置２０ｃが実行する処理について説明する。図３６は、実施の形態３に係る構造解析装置２０ｃが実行する処理のフローチャートである。

形状取得部５６及び要素分割部５８は、それぞれ、実施の形態１の形状情報取得処理（Ｓ１）及び要素分割処理（Ｓ２）と同様に、形状情報取得処理（Ｓ７１）及び要素分割処理（Ｓ７２）を実行する。

熱解析部６０ａは、実施の形態１の熱解析処理（図１１のＳ３）と同様に、熱解析処理（Ｓ７３）を実行する。熱解析処理（Ｓ７３）の詳細は、実施の形態１の熱解析処理の詳細（図１３参照）と同じであるため、ここでの詳細な説明は省略する。

変換部６４ｃは、図１８に詳細を示す熱特性値−物性値の変換処理（図１１のＳ４）と同様の処理を行い、これによって、熱特性値を物性値に変換する（Ｓ７４）。

更に、変換部６４ｃは、熱特性値を強度に変換する（Ｓ７５）。実施の形態３に係る熱特性値−強度変換処理（Ｓ７５）の詳細を図３７に示すフローチャートを参照して、説明する。

熱特性情報取得部７２ａは、熱特性情報３６ａ，３６ｂを取得する（Ｓ８１）。より詳細には、熱流束ベクトル取得部７４ａが熱流束ベクトル情報３６ａを、温度情報取得部７４ｂが温度情報３６ｂをそれぞれ取得する。

変換条件情報取得部７８は、変換条件記憶部３８から変換条件情報３８ｂを取得する（Ｓ８２）。強度算出部８４ｂは、熱流束ベクトル情報３６ａと、温度情報３６ｂと、変換条件情報３８ａとに基づいて、各要素の強度を算出する（Ｓ８３）。ここで算出される強度は、実施の形態２の強度算出処理（Ｓ６４）において算出される強度と同じであり、繊維方向の強度の大きさ、繊維と垂直な方向の強度の大きさ、及び、繊維方向に従って強度の配向性を示す値を含む。続けて、強度算出部８４ｂは、算出した各要素の強度を構造解析結果情報４４ｃとして構造解析結果記憶部４４に格納する（Ｓ８４）。

再び図３６を参照すると、構造解析部６６は、熱特性値−物性値変換処理（Ｓ７４）において算出された物性値を用いて部材の構造解析を行い、それによって、ユーザが想定する条件下での部材の応力及び変位を算出する（Ｓ７６）。構造解析処理（Ｓ７６）の詳細は、図１９に示す実施の形態１の構造解析処理と同様であるため、ここでの詳細な説明は省略する。

安全率算出部８６は、各要素の安全率を算出する（Ｓ７７）。実施の形態３に係る安全率算出処理（Ｓ７７）の詳細なフローチャートを図３８に示す。安全率算出部８６は、構造解析部６６から各要素の強度を含む強度情報を取得する（Ｓ９１）。安全率算出部８６は、構造解析結果記憶部４４から各要素の応力情報を取得する（Ｓ９２）。

安全率算出部８６は、各要素の強度を対応する要素の応力で割ることによって、各要素の安全率を算出する（Ｓ９３）。

ここで、強度は、上述のように配向性を有する。また、応力も各方向について算出できる。そのため、各要素の安全率は、例えば、各要素の繊維の方向とそれに垂直な方向とのそれぞれについて算出され、小さい方を当該要素の安全率とする。

なお、安全率の算出には、例えば、各方向の強度を１つの大きさに変換した強度を用いてもよい。この場合、安全率の算出に用いる応力には、各要素に生じる応力の絶対値を用いるとよい。

安全率算出部８６は、算出した各要素の安全率を構造解析結果記憶部４４に格納する（Ｓ９４）。

再び図３６を参照する。表示制御部７０は、例えば各要素の安全率の大きさの違いを色の変化に対応付けて表すコンタ図、各要素における安全率の絶対値のリスト等の形式で表示部２２に表示させ（Ｓ７８）、処理を終了する。

以上のように、本実施の形態では、繊維を含む部材の要素毎の安全率を算出する。そのため、安全率の高精度なシミュレーションが実現できる。また、このような安全率の算出に用いる要素毎の弾性率及び強度は熱解析に基づいて算出されるため、その算出に要する時間は、流動解析等によって算出する場合よりも短い。従って、部材の安全率について、短時間で算出することが可能になる。さらに、各要素の繊維の配向性を考慮した安全率を算出できる。そのため、高精度な安全率の算出が可能になり、製品及び製品を射出成形するための金型の高精度な設計が可能になる。

なお、本発明の各実施の形態において、解析の結果は、要素毎に算出されることとしたが、各要素を分割する境界の交点である節点毎に算出されてもよい。これにより、適切な計算方法を適宜選択することが可能になる。

本発明は、繊維状の物質を含み、射出成形される材料を用いて作られる部材等の構造解析に利用することができる。

本発明の実施の形態１に係る構造解析装置の外観を示す図。 実施の形態１に係る構造解析装置のブロック図。 実施の形態１に係る要素情報の具体例を示す図。 実施の形態１に係る熱解析条件情報の具体例を示す図。 実施の形態１に係る熱流束ベクトル情報の具体例を示す図。 実施の形態１に係る温度情報の具体例を示す図。 実施の形態１に係る変換条件情報の具体例を示す図。 実施の形態１に係る物性情報の具体例を示す図。 実施の形態１に係る構造解析結果情報の具体例を示す図。 実施の形態１に係る変換部の詳細を示すブロック図。 実施の形態１に係る構造解析装置が実行する処理のフローチャート。 部材としてのファン９０ａを要素分割した例を示す図。 熱解析処理Ｓ３の詳細を示すフローチャート。 ファンの熱解析によって得られた各要素の熱流束ベクトルの一部を示す図。 実施の形態１で例示するファンの正面を示す写真。 図１５に示される矩形部分のＸ線ＣＴ装置による画像を拡大して示す図。 Ｘ線ＣＴ装置による画像に基づいて作成したファンの繊維の配向性を示す図。 実施の形態１の変換処理の詳細を示すフローチャート。 構造解析処理Ｓ５の詳細を示すフローチャート。 本発明の実施の形態２に係る構造解析装置の外観を示す図。 実施の形態２に係る構造解析装置のブロック図。 実施の形態２に係る要素情報の具体例を示す図。 実施の形態２に係る熱解析条件情報の具体例を示す図。 実施の形態２に係る温度勾配情報の具体例を示す図。 実施の形態２に係る熱流束ベクトル情報の具体例を示す図。 実施の形態２に係る温度情報の具体例を示す図。 実施の形態２に係る変換条件情報の具体例を示す図。 実施の形態２に係る変換部の詳細を示すブロック図。 ２つのゲートを熱源として熱解析によって算出された温度分布を示す図。 実施の形態２に係る構造解析装置が実行する処理のフローチャート。 実施の形態２に係る熱解析処理Ｓ４３の詳細を示すフローチャート。 実施の形態２の変換処理Ｓ４４の詳細を示すフローチャート。 実施の形態３に係る構造解析装置２０ｃのブロック図。 実施の形態３に係る変換条件情報の具体例を示す図。 実施の形態３に係る変換部の詳細を示すブロック図。 実施の形態３に係る構造解析装置が実行する処理のフローチャート。 実施の形態３に係る熱特性値−強度変換処理の詳細を示すフローチャート。 実施の形態３に係る安全率算出処理Ｓ７７の詳細を示すフローチャート。

符号の説明

２０ａ，２０ｂ，２０ｃ 構造解析装置、２２ 表示部、２４ 入力部、２６ ＣＡＤ装置、３２ 要素記憶部、３４ 熱解析条件記憶部、３６ 熱特性記憶部、３８ 変換条件記憶部、４０ 物性記憶部、４２ 構造解析条件記憶部、４４ 構造解析結果記憶部、５２ 入力情報取得部、５６ 形状取得部、５８ 要素分割部、６０ａ，６０ｂ 熱解析部、６４ａ，６４ｂ，６４ｃ 変換部、６６ 構造解析部、７０ 表示制御部、７２ａ，７２ｂ 熱特性情報取得部、７４ａ 熱流束ベクトル情報取得部、７４ｂ 温度情報取得部、７４ｃ 温度勾配情報取得部、７８ 変換条件情報取得部、８０ 物性値算出部、８２ ウェルド位置特定部、８４ａ，８４ｂ 強度算出部、８６ 安全率算出部。

Claims (13)

1. 繊維状物質を含む材料を用いて射出成形される部材の構造解析方法であって、
   上記部材の形状情報を取得する形状取得ステップと、
   上記部材を複数の要素に分割し、分割した各要素を示す要素情報を生成する要素分割ステップと、
   上記要素情報が示す上記各要素を用いて、上記材料を流し込む金型のゲートに対応する部分に熱源があるとして熱解析を行うことによって上記部材の熱特性を示す熱特性値を算出し、算出した上記熱特性値を含む熱特性情報を生成する熱解析ステップと、
   上記熱特性情報を取得し、取得した上記熱特性情報を構造解析に用いる物性値を含む物性情報に変換する変換ステップと、
   上記物性情報に含まれる上記物性値と、上記要素情報が示す上記各要素とを用いて上記部材の構造解析を行う構造解析ステップとを含むことを特徴とする構造解析方法。
2. 上記熱解析ステップにおいて、上記各要素の熱流束ベクトルを示す熱流束ベクトル情報を含む上記熱特性情報を生成し、
   上記変換ステップにおいて、
   上記熱流束ベクトルを含む上記熱特性情報を取得するステップと、
   上記熱流束ベクトルを上記繊維状物質の配向性とすることによって、上記繊維状物質の配向性を考慮した物性値を含む物性情報に変換するステップとを含むことを特徴とする請求項１に記載の構造解析方法。
3. 上記熱解析ステップにおいて、上記各要素の温度を示す温度情報を更に含む上記熱特性情報を生成し、
   上記変換ステップは、
   上記温度情報を含む上記熱特性情報を取得するステップと、
   上記部材における物性値の最大値及び最小値を含む変換条件情報を取得するステップと、
   上記温度情報に含まれる温度が高い要素である程、当該要素の上記物性値が大きくなるように、上記温度を上記物性値に一対一で対応付けて変換するとを含むことを特徴とする請求項２に記載の構造解析方法。
4. 上記構造解析ステップにおいて、ユーザによって設定された境界条件及び算出された上記各要素の物性値を参照して、上記各要素の応力及び／又は変位を算出することを特徴とする請求項１から請求項３のいずれか１項に記載の構造解析方法。
5. 上記変換ステップは、
   上記部材における弾性率の最大値及び最小値と、密度の最大値及び最小値とを含む変換条件情報を取得するステップと、
   上記温度情報に含まれる最高温度を上記弾性率の最大値に、最低温度を上記弾性率の最小値に対応付け、上記最高温度と上記最低温度とのそれぞれに重み付けをして各要素の温度が得られる場合に、上記弾性率の最大値と上記弾性率の最小値とに上記各要素の温度を得た場合と同じ重み付けをして得られる値を前記各要素の弾性率として算出するステップと、
   上記熱流束ベクトルを上記繊維状物質の配向性とすることによって、上記繊維状物質の配向性を考慮した前記各要素の弾性率を含む物性情報に変換するステップと、
   上記温度情報に含まれる最高温度を上記密度の最大値に、最低温度を上記密度の最小値にそれぞれ対応付け、上記最高温度と上記最低温度とのそれぞれに重み付けをして各要素の温度が得られる場合に、上記密度の最大値と上記密度の最小値とに上記各要素の温度を得た場合と同じ重み付けをして得られる値を前記各要素の密度として算出するステップとを含むことを特徴とする請求項３又は４に記載の構造解析方法。
6. 上記変換ステップは、
   上記部材における物性値の最大値及び最小値と、上記部材における強度の最大値及び最小値とを含む変換条件情報を取得するステップと、
   上記温度情報に含まれる最高温度を上記物性値の最大値に、最低温度を上記物性値の最小値にそれぞれ対応付け、上記最高温度と上記最低温度とのそれぞれに重み付けをして各要素の温度が得られる場合に、上記物性値の最大値と上記物性値の最小値とに上記各要素の温度を得た場合と同じ重み付けをして得られる値を前記各要素の物性値として算出するステップと、
   上記温度情報に含まれる最高温度を上記強度の最大値に、最低温度を上記強度の最小値にそれぞれ対応付け、上記最高温度と上記最低温度とのそれぞれに重み付けをして各要素の温度が得られる場合に、上記構造特性値の最大値と上記強度の最小値とに上記各要素の温度を得た場合と同じ重み付けをして得られる値を前記各要素の強度として算出するステップとを含み、
   上記構造解析ステップにおいて、予め設定された境界条件及び算出された上記各要素の各物性値を参照して、上記各要素の応力を算出し、
   更に、
   上記各要素について上記強度を上記応力で割ることによって、上記各要素の安全率を算出するステップとを含むことを特徴とする請求項５に記載の構造解析方法。
7. 上記変換ステップは、
   上記材料でウェルドが発生した部分の強度が上記材料の強度に比べてどの程度であるがを示す強度倍率を更に含む変換条件情報を取得するステップと、
   上記温度情報に含まれる上記各要素の温度が上記最低温度と等しいか又は閾値以下であるかを判断し、上記最低温度と等しい又は上記閾値以下である要素をウェルド位置として特定するステップと、
   ウェルド位置であると特定された要素について、上記温度情報に基づいて得られた強度に上記強度倍率を掛けることによって上記ウェルド位置の強度を修正し、修正したウェルド位置の強度を含む強度情報を生成するステップとを含むことを特徴とする請求項６に記載の構造解析方法。
8. 上記熱解析ステップにおいて、上記各要素の温度勾配を示す温度勾配情報を含む上記熱特性情報を生成し、
   上記変換ステップは、
   上記温度勾配情報を含む上記熱特性情報を取得するステップと、
   上記材料でウェルドが発生した部分の強度が上記材料の強度に比べてどの程度であるがを示す強度倍率を更に含む変換条件情報を取得するステップと、
   上記温度勾配情報に含まれる上記各要素の温度勾配が閾値以下であるかを判断し、上記閾値以下である要素をウェルド位置として特定するステップと、
   ウェルド位置であると特定された要素について、上記温度情報に基づいて得られた強度に上記強度倍率を掛けることによって上記ウェルド位置の強度を修正し、修正したウェルド位置の強度を含む強度情報を生成するステップとを含むことを特徴とする請求項６に記載の構造解析方法。
9. 射出成形された部材を複数の要素に分割し、分割した当該複数の要素を用いて当該部材について構造解析を行うための物性値を算出する物性値算出方法であって、
   上記構造解析と同一の上記複数の要素を用いた熱解析によって得られた部材の熱特性を示す熱特性情報を、上記構造解析をするための上記各要素の物性値を示す物性情報に変換することを特徴とする構造解析のための物性値算出方法。
10. 繊維状物質を含む材料を用いて射出成形される部材の構造解析装置であって、
    上記部材の形状情報を取得する形状取得手段と、
    上記部材を複数の要素に分割し、分割した各要素を示す要素情報を生成する要素分割手段と、
    上記要素情報が示す上記各要素を用いて、上記材料を流し込む金型のゲートに対応する部分に熱源があるとして熱解析を行うことによって、上記部材の熱特性を示す熱特性情報を生成する熱解析手段と、
    上記熱特性情報を取得し、取得した上記特性情報を構造解析に必要な物性値を含む物性情報に変換する変換手段と、
    上記物性情報に含まれる上記物性値と、上記要素情報が示す上記各要素とを用いて上記部材の構造解析を行う構造解析手段とを含むことを特徴とする構造解析装置。
11. 繊維状物質を含む材料を用いて射出成形される部材の構造解析プログラムであって、
    上記部材の形状情報を取得する形状取得ステップと、
    上記部材を複数の要素に分割し、分割した各要素を示す要素情報を生成する要素分割ステップと、
    上記要素情報が示す上記各要素を用いて、上記材料を流し込む金型のゲートに対応する部分に熱源があるとして熱解析を行うことによって、上記部材の熱特性を示す熱特性情報を生成する熱解析ステップと、
    上記熱特性情報を取得し、取得した上記特性情報を構造解析に必要な物性値を含む物性情報に変換する変換ステップと、
    上記物性情報に含まれる上記物性値と、上記要素情報が示す上記各要素とを用いて上記部材の構造解析を行う構造解析ステップとをコンピュータに実行させることを特徴とする構造解析プログラム。
12. 射出成形された部材を複数の要素に分割し、分割した当該複数の要素を用いて当該部材を構造解析するための物性値を算出する物性値算出装置であって、
    上記構造解析と同一の上記複数の要素を用いた熱解析によって得られた部材の熱特性を示す熱特性情報を、上記構造解析をするための上記各要素の物性値を示す物性情報に変換することを特徴とする構造解析のための物性値算出装置。
13. 射出成形された部材を複数の要素に分割し、分割した当該複数の要素を用いて当該部材を構造解析するための物性値を算出する物性値算出プログラムであって、
    上記構造解析と同一の上記複数の要素を用いた熱解析によって得られた部材の熱特性を示す熱特性情報を、上記構造解析をするための上記各要素の物性値を示す物性情報に変換することをコンピュータに実行させることを特徴とする構造解析のための物性値算出プログラム。

Images