衝撃解析方法
【課題】 樹脂部品の複雑な変形過程の衝撃解析を精度よく実現する衝撃解析方法を提供する。
【解決手段】 この衝撃解析方法は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価するものである。そして、樹脂成形体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行うことを特徴とする。延性破壊条件式に応力多軸度を考慮することができる延性破壊条件式を用いるようにしてもよい。また、樹脂特有の粘弾性に対処するために、延性破壊条件式に歪み速度依存性を考慮するようにしてもよい。
【解決手段】 この衝撃解析方法は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価するものである。そして、樹脂成形体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行うことを特徴とする。延性破壊条件式に応力多軸度を考慮することができる延性破壊条件式を用いるようにしてもよい。また、樹脂特有の粘弾性に対処するために、延性破壊条件式に歪み速度依存性を考慮するようにしてもよい。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、樹脂成形体の衝撃負荷時の変形挙動を予測するために用いる衝撃解析方法に関する。
【背景技術】
【0002】
例えば、自動車用の樹脂内装部品を設計する場合、手間の掛かる実地試験の代わりにコンピュータを用いた衝撃解析方法を活用するのが望ましい。しかしながら、自動車用の樹脂内装部品は、リブ構造のような複雑な形状を有しており、さらに、衝突試験等においては、衝撃エネルギーを吸収するため大きく変形し、部分的に破壊に至る。また、樹脂材料は金属材料と異なり、粘弾性特性を有するため、物性値は歪み速度による影響を受ける。従って、樹脂部品の複雑な変形過程の衝撃解析を精度よく実現する物性値を設定することは難しかった。
【0003】
上記のような衝突試験における衝撃特性を解析するための手法として提案された特許文献1の技術は、速度条件ごとに所定形状の樹脂成形体を用いた実地試験と、破断塑性歪みを変えた衝撃解析とを行い、実地試験結果と衝撃解析結果を比較して速度条件ごとの破断塑性歪みを決定するものである。これにより、歪み速度に対応したリブ部の破壊判定が可能となった。しかしながら、樹脂成形体のリブ部と基板部が一体となっていて、リブ部と基板部双方が破壊するような複雑な破壊挙動は、衝撃解析で再現することができなかった。
【0004】
また、衝突試験における衝撃特性を解析するための手法としては特許文献2のように種々の引張試験より得られる歪み速度と破断歪みの関係より破壊判定を行っているものもある。
【0005】
一方、金属に関しては、延性破壊条件式を用いてその変形挙動を定量的に把握する試みがなされてきた。また、特許文献3では、弾塑性体(金属材料)の相当塑性歪み−応力多軸度の関係を考慮し、境界条件として与えられる外力を上下に更新しながら、各要素に対して亀裂発生条件を満足するか判定を行い、亀裂進展予測を行っている。
【0006】
【特許文献1】特開2003-279456号公報
【特許文献2】特開2005-337784号公報
【特許文献3】特開2004-069638号公報
【非特許文献1】後藤 学 「延性破壊条件式について」 塑性と加工(日本塑性加工学会誌) 第38巻 第434号 1997年3月 P200〜205
【非特許文献2】大矢根 守哉 「延性破壊の条件式について」 日本機械学会誌 第75巻 第639号 昭和47年4月 P596〜601
【発明の開示】
【発明が解決しようとする課題】
【0007】
本発明は、前記事情に鑑みて為されたもので、樹脂成形体の複雑な変形過程の衝撃解析を精度よく実現する衝撃解析方法を提供することを目的とする。
【課題を解決するための手段】
【0008】
すなわち本発明は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、樹脂成形体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行うことを特徴とする衝撃解析方法である。
【0009】
また、本発明は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、有限要素モデルの各要素について、延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値を算出する第1の工程と、第1の工程で求めた延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値とを比較して破壊判定を行う第2の工程を有することを特徴とする衝撃解析方法である。
【発明の効果】
【0010】
本発明の衝撃解析方法によれば、樹脂成形体の複雑な変形過程の衝撃解析を精度よく実現することができる。
【発明を実施するための最良の形態】
【0011】
請求項1にかかる発明は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、樹脂成形体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行うことを特徴とする衝撃解析方法である。
【0012】
延性破壊条件式は、金属の破壊挙動を定量的に理解するために、開発されてきた手段である。延性破壊は所定量の変形の後に破壊が起きる点で、脆性破壊等の他の破壊様式と異なっている。一般に延性破壊では、(1)材料内の高い三軸引張場のところでボイドが発生し、(2)新たな発生を伴いながら、成長、合体してボイドが大きくなり、(3)やがて材料が破壊に至るというのがメカニズムである。
【0013】
延性破壊条件式は、主に、塑性加工時の破壊を防ぐことを目的として開発されてきた。延性破壊条件式は、上記のような延性破壊の過程をモデル化し、数式化しようとするもので、非特許文献1に記載されているように、多くのものが提案されている。例えば、非特許文献2において、大矢根は、ボイド欠陥が生成・合体・成長していく過程を多孔質体で近似し、塑性変形によって体積変化をきたすとして、その体積変化の履歴を追跡することで破壊条件式が得られるとして、以下の式(1)を提案した。
【数1】
【0014】
発明者等は、このような延性破壊条件式が、樹脂成形体の破壊に応用できるのではないかと考え、本発明に至った。樹脂成形体は条件によっては延性的に破壊すると考えられるからである。勿論、金属と樹脂では相違点も大きい。例えば、樹脂は、初期状態において既に多くのボイドを有しており、また、粘弾性を有している点も弾性体である金属とは異なる。
【0015】
しかしながら、前者については、例えば式(1)は本来ボイド欠陥を有している多孔質体をモデルとしているので、基本的に問題が無いと考えられる。延性破壊条件式の中でも、式(1)のように、3次元的な変形を扱うのに必要と考えられる応力多軸度を考慮した式を用いることによって、より精度良く解析を行うことができる。
【数2】
【0016】
樹脂の粘弾性を考慮するためには、歪み速度依存性を考慮した式を用いることが好ましい。例えば、式(2−1)または式(2−2)に示すように、式(1)の右辺に、歪み速度に依存する因子を入れることにより、より精度良く解析を行うことができる。式(2−1)や式(2−2)において、左辺を対数とし、右辺の歪み速度を対数としたのは、一軸引張試験において破断歪みと歪み速度の関係が両対数上で直線近似できるからである。また、ある歪み速度以上では、樹脂はガラス状態になるので、その際は右辺をある一定の値とする。
【0017】
【数3】
【0018】
樹脂の粘弾性特性および、歪み速度が大きい際にガラス状態になることを表現するために、例えば、式(2−1)および式(2−2)を示したが、実験事実を表現できるようにした他の延性破壊条件式を用いるようにしても良い。例えば、大矢根の式(1)を右辺の歪み速度のみを対数とするように変形したものも、場合によっては、有効と思われる。
【数4】
【0019】
本発明の衝撃解析方法を適用する樹脂成形体が、リブ部と基板部とを有する樹脂成形体である場合、該樹脂成形体全体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行っても良いが、通常、樹脂成形体のリブ部と基板部では衝撃の加わる方向が異なるため、リブ部と基板部では異なる破壊判定法を用いて衝撃解析を行うことが好ましい。リブ部と基板部とを有する樹脂成形体では、基板部では板面に直交する方向に衝撃力が作用するのに対し、リブ部では、概ね板面に沿って衝撃力が作用する。従って、破断に至る過程も異なるので、それぞれに異なる破壊判定法を用いることが好ましい。特に、基板部では前述した延性破壊方向式を用いて破壊判定を行い、リブ部については異なる破壊判定法を用いることが好ましい。
【0020】
例えば、樹脂成形品にリブ部が有る場合は、リブ部の破壊判定法として、所定形状の樹脂成形体を用いた実地試験と、破断塑性歪みを変えた衝撃解析を複数行い、実地試験結果と衝撃解析結果を比較して決定した破断塑性歪みを用いることが好ましい。
【0021】
また、樹脂成形品にリブ部が有る場合に、リブ部の破壊判定法として、種々の引張速度での一軸引張試験を行うことで得られる破断歪みと歪み速度との関係を用いることも好ましい。
【0022】
また、本発明の他の態様は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、有限要素モデルの各要素について、延性破壊条件式に関する破壊判定値(式(2−1)または式(2−2)の左辺)と延性破壊条件式に関する破壊判定基準値(式(2−1)または式(2−2)の右辺)を算出する第1の工程と、第1の工程で求めた延性破壊条件式に関する破壊判定値(式(2−1)または式(2−2)の左辺)と延性破壊条件式に関する破壊判定基準値(式(2−1)または式(2−2)の右辺)とを比較して破壊判定を行う第2の工程を有する。
【0023】
以下、図面を参照しつつ発明の好適な実施の形態をさらに詳しく説明する。この実施例では、図1に示すような樹脂成形体の衝撃特性および変形挙動を解析プログラムを用いた衝撃解析手法で予測する。図1には、衝撃解析の対象となる樹脂成形体10を正面から見た形状が、コンピュータグラフィックスを用いて描かれている。この樹脂成形体10は、リブ部12と、平板状の基板部14とからなっている。リブ部12は、一定厚さの複数のリブ板が基板部14からこれに対して直交するように延び、規則的な格子を形成するように構成されている。格子の配列は幅方向×長方向が3×9であり、格子間隔は30mmであり、リブの高さは25mmである。リブ部12には型抜きのためのテーパが形成され、板厚は基端側(基板部側)が1.12mm、先端側が0.69mmである。基板部14は幅150mm、長さ310mm、板厚2mmである。解析モデルを作成する際は、リブ部はシェル要素を、基板部はソリッド要素を用いることが望ましい。
【0024】
この樹脂成形体10について、樹脂材料の機械的特性を示す物性値を設定し、図2に示すような衝突試験における衝撃特性、変形挙動を解析プログラムを用いて解析する。この試験では成形体の先端側をバリア(固体障壁)16に向けて対向配置し、基端側からダート18を5.7m/sと4.5m/sで前進させて成形体に衝突させ、その際の成形体の変形と負荷される荷重を測定する。ダート18は、質量6.8kgで直径165mmの半球の衝突体である。
【0025】
樹脂材料の機械的特性を示す物性値には、弾性率、ポアソン比、密度、降伏応力、歪み硬化度、および延性破壊条件式の材料定数が含まれる。これらの物性値のうち、降伏応力については種々の引張速度で一軸引張試験を実施することで決定できる。樹脂素材では、降伏応力には通常歪み速度依存性が顕著に見られ、一般に歪み速度が大きいほど降伏応力も大きくなる。
【0026】
歪み硬化度については応力多軸度性を考慮することができる。歪み硬化度は降伏後真応力が一定である場合を0とし、真応力が一軸引張試験結果の真応力に従う場合を1とする。
【0027】
延性破壊条件式は、上述したように、大矢根の式(1)を変形した式(2−1)および式(2−2)を用いる。
【数5】
【0028】
上記の物性値のうち、延性破壊条件式の材料定数を決定する方法としては、例えば、以下のような方法を採用することができる。
【数6】
【0029】
【数7】
【0030】
解析モデルの要素は無限に小さくすることはできないため、試験時に亀裂先端で発生する応力や歪みなどが解析では再現できない。そのため、割れの起点は再現できても、割れの伝播は再現できない場合がある。その場合、破壊した要素の周囲は破壊が起こりやすくなるように設定して、計算を行っても良い。例えば、破壊した要素の周囲にある要素の破壊判定を行う際は式(2−1)または式(2−2)の左辺の値に重み付けを行って、破壊しやすくなるようにする。
【数8】
【0031】
また、基板部の1つのソリッド要素が破壊している場合に、その要素の基板部の厚み方向にある全てのソリッド要素を削除するように設定して計算を行っても良い。
【0032】
以下、この実施の形態の衝撃解析方法の工程を、図3のフロー図を用いて説明する。
まず、ステップ1において、有限要素モデルを作成し、静的物性値のデータ、降伏応力の歪み速度依存性のデータ、歪み硬化度の応力多軸度性のデータおよび延性破壊条件式の材料定数を入力し、汎用有限要素プログラムLS-DYNAで計算をスタートする。次に、ステップ2において、ある時刻tが終了時間以下である場合には、各時間ステップで決定される時間増分Δtごとに時間を進める。そして、ステップ3において、汎用有限要素プログラムLS-DYNAで節点の変位を計算し、ステップ4において、時刻tにおける要素を順次選択し、ステップ5においてその要素について、応力と延性破壊条件式の破壊判定値と破壊基準値を計算する。
【0033】
次に、ステップ6において、ステップ5で求めた値が式(2−1)および式(2−2)を充足するかどうかを判定する。破壊判定値が破壊基準値より大きいか等しい場合には、ステップ7においてその要素を削除し、破壊判定値が破壊基準値より小さい場合にはそのままで、ステップ4に戻り、未計算の要素が残っている限り、ステップ5において次の要素に関する応力と延性破壊条件式の破壊判定値と破壊基準値を計算する。このようにして、終了時間までステップ3〜ステップ7を繰り返す。
【0034】
図4は、破壊した要素の周囲を破壊しやすくする場合のフロー図である。この方法では、ステップ5でその要素について応力と延性破壊条件式の破壊判定値と破壊基準値を計算した後、ステップ6でその要素と接する要素が破壊しているかどうかを判定する。そして、破壊している場合は、ステップ7において破壊判定値に重み付けを行い、破壊していない場合はそのままで、ステップ8において破壊条件式を満たすかどうかを判定する。そして、破壊判定値が破壊基準値より大きいか等しい場合には、ステップ9においてその要素を削除する。その他の工程は、先の場合と同様である。
【0035】
基板部の1つのソリッド要素が破壊したらその要素の基板部の厚み方向にある全てのソリッド要素を削除する場合は、ステップ8において破壊条件式を満たさない場合において、図に破線で示すように、ステップ10で基板部の同じ厚み方向にある要素が破壊しているかを判定し、破壊している場合はステップ9に戻ってその要素を削除するようにする。
【実施例】
【0036】
樹脂成形体の素材として、住友化学株式会社製ノーブレンAZ864E4を用いた場合の衝撃特性および破壊挙動を評価する衝撃解析方法を以下に示す。
(1)解析に用いたソフトウェア
解析ソフト:LS-DYNA version9.70(Livermore
Software Technology Corporation製)
(2)解析方法
空間の離散化:有限要素法
時間積分:中心差分に基づく陽解法
材料モデル:LS-DYNA物性タイプ44(ユーザー定義物性)
(3)物性値
AZ864E4の物性値としては以下のような値を用いた。
弾性率:1000MPa
降伏応力:22.8MPa
ポアソン比:0.40
比重:0.90
真応力−真塑性歪み曲線:図5
降伏応力−歪み速度:図6
歪み硬化度−応力多軸度:図7
延性破壊条件式の材料定数:C1=0.970、C2=67.2、
C3=−0.128、C4=0.168、
C5=1.51×1014、C6=−1.65
なお、延性破壊条件式(2−1)、(2−2)の材料定数は、上記の材料定数の決定方法の一例で示したように、決めたものである。
【0037】
解析モデルの要素が基板部の割れの伝播を再現するには大きいと考え、破壊した要素の周囲は破壊が起こりやすくした。
重み係数は、W=3.0とした。
また、基板部のソリッド要素が破壊すると、その要素の基板部の厚み方向にある全てのソリッド要素も削除するようにして計算を行った。
【0038】
図1および図2に示す解析対象モデルの成形体10を実際に作成し、衝突速度5.7m/sと4.5m/sにおいて実地試験を行った。次に図2に示す解析モデルに対して上記の物性値を樹脂成形体10に用いて衝突速度5.7m/sと4.5m/sにおいて、本発明の衝撃解析方法で解析を行った(実施例)。具体的には、歪み速度がC5以下の場合は式(2−1)を、C5を超える場合は式(2−2)を用いた。
【数9】
【0039】
また、比較例として、図2に示す解析モデルに対して上記の物性値を樹脂成形体10に用いて衝突速度5.7m/sと4.5m/sにおいて、LS-DYNA物性タイプ24(多直線近似等方弾塑性体)で、解析を行った。物性タイプ24では降伏応力の歪み速度依存性を考慮することができ、破断塑性歪みは歪み速度に寄らず一定であるモデルである。破断塑性歪みのデータとしては比較例1では0.56とし、比較例2では0.57とした。
【0040】
実地試験、実施例、比較例1および2の結果である荷重−変位曲線を図8に、基板部要素の破壊の有無について表1に示す。図8において(a)は、速度5.7m/sの結果、(b)は速度4.5m/sの結果である。
解析対象モデルの成形体10の実施例、比較例1および2の衝撃解析結果後の形状を図9(a)、(b)および(c)に示す。
速度5.7m/sおよび速度4.5m/sにおける基板部の破壊状況を図10および図11に示す。図10,図11において(a)は実地試験の結果、(b)は実施例の方法による解析結果、(c)は比較例1の結果、(d)は比較例2の結果である。
【0041】
【表1】
【0042】
これらの結果より、比較例の方法では、速度によって基板部が破壊するか破壊しないかの判定を行うことができなかったが、実施例の方法では判定を行うことができた。また、図8の荷重−変位曲線もよく一致していることが分かる。
【図面の簡単な説明】
【0043】
【図1】衝撃解析を行う対象の樹脂成形体である。
【図2】衝撃試験の構成の概略を示す概念図である。
【図3】この実施の形態の衝撃解析方法の工程を示すフロー図である。
【図4】この実施の形態の衝撃解析方法の工程を示すフロー図である。
【図5】樹脂成形体の供試材料の真応力−真塑性歪み曲線である。
【図6】樹脂成形体の供試材料の降伏応力−歪み速度曲線である。
【図7】樹脂成形体の供試材料の歪み硬化度−応力多軸度曲線である。
【図8】樹脂成形体の衝撃特性を評価するための荷重−変位曲線である。
【図9】樹脂成形体の破壊状況を示す図である。
【図10】速度5.7m/sでの衝撃試験における樹脂成形体の基板部の破壊状況を示す図である。
【図11】速度4.5m/sでの衝撃試験における樹脂成形体の基板部の破壊状況を示す図である。
【符号の説明】
【0044】
10 樹脂成形体
12 リブ部
14 基板部
16 バリア(固体障壁)
18 ダート
【技術分野】
【0001】
本発明は、樹脂成形体の衝撃負荷時の変形挙動を予測するために用いる衝撃解析方法に関する。
【背景技術】
【0002】
例えば、自動車用の樹脂内装部品を設計する場合、手間の掛かる実地試験の代わりにコンピュータを用いた衝撃解析方法を活用するのが望ましい。しかしながら、自動車用の樹脂内装部品は、リブ構造のような複雑な形状を有しており、さらに、衝突試験等においては、衝撃エネルギーを吸収するため大きく変形し、部分的に破壊に至る。また、樹脂材料は金属材料と異なり、粘弾性特性を有するため、物性値は歪み速度による影響を受ける。従って、樹脂部品の複雑な変形過程の衝撃解析を精度よく実現する物性値を設定することは難しかった。
【0003】
上記のような衝突試験における衝撃特性を解析するための手法として提案された特許文献1の技術は、速度条件ごとに所定形状の樹脂成形体を用いた実地試験と、破断塑性歪みを変えた衝撃解析とを行い、実地試験結果と衝撃解析結果を比較して速度条件ごとの破断塑性歪みを決定するものである。これにより、歪み速度に対応したリブ部の破壊判定が可能となった。しかしながら、樹脂成形体のリブ部と基板部が一体となっていて、リブ部と基板部双方が破壊するような複雑な破壊挙動は、衝撃解析で再現することができなかった。
【0004】
また、衝突試験における衝撃特性を解析するための手法としては特許文献2のように種々の引張試験より得られる歪み速度と破断歪みの関係より破壊判定を行っているものもある。
【0005】
一方、金属に関しては、延性破壊条件式を用いてその変形挙動を定量的に把握する試みがなされてきた。また、特許文献3では、弾塑性体(金属材料)の相当塑性歪み−応力多軸度の関係を考慮し、境界条件として与えられる外力を上下に更新しながら、各要素に対して亀裂発生条件を満足するか判定を行い、亀裂進展予測を行っている。
【0006】
【特許文献1】特開2003-279456号公報
【特許文献2】特開2005-337784号公報
【特許文献3】特開2004-069638号公報
【非特許文献1】後藤 学 「延性破壊条件式について」 塑性と加工(日本塑性加工学会誌) 第38巻 第434号 1997年3月 P200〜205
【非特許文献2】大矢根 守哉 「延性破壊の条件式について」 日本機械学会誌 第75巻 第639号 昭和47年4月 P596〜601
【発明の開示】
【発明が解決しようとする課題】
【0007】
本発明は、前記事情に鑑みて為されたもので、樹脂成形体の複雑な変形過程の衝撃解析を精度よく実現する衝撃解析方法を提供することを目的とする。
【課題を解決するための手段】
【0008】
すなわち本発明は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、樹脂成形体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行うことを特徴とする衝撃解析方法である。
【0009】
また、本発明は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、有限要素モデルの各要素について、延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値を算出する第1の工程と、第1の工程で求めた延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値とを比較して破壊判定を行う第2の工程を有することを特徴とする衝撃解析方法である。
【発明の効果】
【0010】
本発明の衝撃解析方法によれば、樹脂成形体の複雑な変形過程の衝撃解析を精度よく実現することができる。
【発明を実施するための最良の形態】
【0011】
請求項1にかかる発明は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、樹脂成形体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行うことを特徴とする衝撃解析方法である。
【0012】
延性破壊条件式は、金属の破壊挙動を定量的に理解するために、開発されてきた手段である。延性破壊は所定量の変形の後に破壊が起きる点で、脆性破壊等の他の破壊様式と異なっている。一般に延性破壊では、(1)材料内の高い三軸引張場のところでボイドが発生し、(2)新たな発生を伴いながら、成長、合体してボイドが大きくなり、(3)やがて材料が破壊に至るというのがメカニズムである。
【0013】
延性破壊条件式は、主に、塑性加工時の破壊を防ぐことを目的として開発されてきた。延性破壊条件式は、上記のような延性破壊の過程をモデル化し、数式化しようとするもので、非特許文献1に記載されているように、多くのものが提案されている。例えば、非特許文献2において、大矢根は、ボイド欠陥が生成・合体・成長していく過程を多孔質体で近似し、塑性変形によって体積変化をきたすとして、その体積変化の履歴を追跡することで破壊条件式が得られるとして、以下の式(1)を提案した。
【数1】
【0014】
発明者等は、このような延性破壊条件式が、樹脂成形体の破壊に応用できるのではないかと考え、本発明に至った。樹脂成形体は条件によっては延性的に破壊すると考えられるからである。勿論、金属と樹脂では相違点も大きい。例えば、樹脂は、初期状態において既に多くのボイドを有しており、また、粘弾性を有している点も弾性体である金属とは異なる。
【0015】
しかしながら、前者については、例えば式(1)は本来ボイド欠陥を有している多孔質体をモデルとしているので、基本的に問題が無いと考えられる。延性破壊条件式の中でも、式(1)のように、3次元的な変形を扱うのに必要と考えられる応力多軸度を考慮した式を用いることによって、より精度良く解析を行うことができる。
【数2】
【0016】
樹脂の粘弾性を考慮するためには、歪み速度依存性を考慮した式を用いることが好ましい。例えば、式(2−1)または式(2−2)に示すように、式(1)の右辺に、歪み速度に依存する因子を入れることにより、より精度良く解析を行うことができる。式(2−1)や式(2−2)において、左辺を対数とし、右辺の歪み速度を対数としたのは、一軸引張試験において破断歪みと歪み速度の関係が両対数上で直線近似できるからである。また、ある歪み速度以上では、樹脂はガラス状態になるので、その際は右辺をある一定の値とする。
【0017】
【数3】
【0018】
樹脂の粘弾性特性および、歪み速度が大きい際にガラス状態になることを表現するために、例えば、式(2−1)および式(2−2)を示したが、実験事実を表現できるようにした他の延性破壊条件式を用いるようにしても良い。例えば、大矢根の式(1)を右辺の歪み速度のみを対数とするように変形したものも、場合によっては、有効と思われる。
【数4】
【0019】
本発明の衝撃解析方法を適用する樹脂成形体が、リブ部と基板部とを有する樹脂成形体である場合、該樹脂成形体全体の破壊を判定するための破壊判定法として、延性破壊条件式を用いて破壊判定を行っても良いが、通常、樹脂成形体のリブ部と基板部では衝撃の加わる方向が異なるため、リブ部と基板部では異なる破壊判定法を用いて衝撃解析を行うことが好ましい。リブ部と基板部とを有する樹脂成形体では、基板部では板面に直交する方向に衝撃力が作用するのに対し、リブ部では、概ね板面に沿って衝撃力が作用する。従って、破断に至る過程も異なるので、それぞれに異なる破壊判定法を用いることが好ましい。特に、基板部では前述した延性破壊方向式を用いて破壊判定を行い、リブ部については異なる破壊判定法を用いることが好ましい。
【0020】
例えば、樹脂成形品にリブ部が有る場合は、リブ部の破壊判定法として、所定形状の樹脂成形体を用いた実地試験と、破断塑性歪みを変えた衝撃解析を複数行い、実地試験結果と衝撃解析結果を比較して決定した破断塑性歪みを用いることが好ましい。
【0021】
また、樹脂成形品にリブ部が有る場合に、リブ部の破壊判定法として、種々の引張速度での一軸引張試験を行うことで得られる破断歪みと歪み速度との関係を用いることも好ましい。
【0022】
また、本発明の他の態様は、解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、有限要素モデルの各要素について、延性破壊条件式に関する破壊判定値(式(2−1)または式(2−2)の左辺)と延性破壊条件式に関する破壊判定基準値(式(2−1)または式(2−2)の右辺)を算出する第1の工程と、第1の工程で求めた延性破壊条件式に関する破壊判定値(式(2−1)または式(2−2)の左辺)と延性破壊条件式に関する破壊判定基準値(式(2−1)または式(2−2)の右辺)とを比較して破壊判定を行う第2の工程を有する。
【0023】
以下、図面を参照しつつ発明の好適な実施の形態をさらに詳しく説明する。この実施例では、図1に示すような樹脂成形体の衝撃特性および変形挙動を解析プログラムを用いた衝撃解析手法で予測する。図1には、衝撃解析の対象となる樹脂成形体10を正面から見た形状が、コンピュータグラフィックスを用いて描かれている。この樹脂成形体10は、リブ部12と、平板状の基板部14とからなっている。リブ部12は、一定厚さの複数のリブ板が基板部14からこれに対して直交するように延び、規則的な格子を形成するように構成されている。格子の配列は幅方向×長方向が3×9であり、格子間隔は30mmであり、リブの高さは25mmである。リブ部12には型抜きのためのテーパが形成され、板厚は基端側(基板部側)が1.12mm、先端側が0.69mmである。基板部14は幅150mm、長さ310mm、板厚2mmである。解析モデルを作成する際は、リブ部はシェル要素を、基板部はソリッド要素を用いることが望ましい。
【0024】
この樹脂成形体10について、樹脂材料の機械的特性を示す物性値を設定し、図2に示すような衝突試験における衝撃特性、変形挙動を解析プログラムを用いて解析する。この試験では成形体の先端側をバリア(固体障壁)16に向けて対向配置し、基端側からダート18を5.7m/sと4.5m/sで前進させて成形体に衝突させ、その際の成形体の変形と負荷される荷重を測定する。ダート18は、質量6.8kgで直径165mmの半球の衝突体である。
【0025】
樹脂材料の機械的特性を示す物性値には、弾性率、ポアソン比、密度、降伏応力、歪み硬化度、および延性破壊条件式の材料定数が含まれる。これらの物性値のうち、降伏応力については種々の引張速度で一軸引張試験を実施することで決定できる。樹脂素材では、降伏応力には通常歪み速度依存性が顕著に見られ、一般に歪み速度が大きいほど降伏応力も大きくなる。
【0026】
歪み硬化度については応力多軸度性を考慮することができる。歪み硬化度は降伏後真応力が一定である場合を0とし、真応力が一軸引張試験結果の真応力に従う場合を1とする。
【0027】
延性破壊条件式は、上述したように、大矢根の式(1)を変形した式(2−1)および式(2−2)を用いる。
【数5】
【0028】
上記の物性値のうち、延性破壊条件式の材料定数を決定する方法としては、例えば、以下のような方法を採用することができる。
【数6】
【0029】
【数7】
【0030】
解析モデルの要素は無限に小さくすることはできないため、試験時に亀裂先端で発生する応力や歪みなどが解析では再現できない。そのため、割れの起点は再現できても、割れの伝播は再現できない場合がある。その場合、破壊した要素の周囲は破壊が起こりやすくなるように設定して、計算を行っても良い。例えば、破壊した要素の周囲にある要素の破壊判定を行う際は式(2−1)または式(2−2)の左辺の値に重み付けを行って、破壊しやすくなるようにする。
【数8】
【0031】
また、基板部の1つのソリッド要素が破壊している場合に、その要素の基板部の厚み方向にある全てのソリッド要素を削除するように設定して計算を行っても良い。
【0032】
以下、この実施の形態の衝撃解析方法の工程を、図3のフロー図を用いて説明する。
まず、ステップ1において、有限要素モデルを作成し、静的物性値のデータ、降伏応力の歪み速度依存性のデータ、歪み硬化度の応力多軸度性のデータおよび延性破壊条件式の材料定数を入力し、汎用有限要素プログラムLS-DYNAで計算をスタートする。次に、ステップ2において、ある時刻tが終了時間以下である場合には、各時間ステップで決定される時間増分Δtごとに時間を進める。そして、ステップ3において、汎用有限要素プログラムLS-DYNAで節点の変位を計算し、ステップ4において、時刻tにおける要素を順次選択し、ステップ5においてその要素について、応力と延性破壊条件式の破壊判定値と破壊基準値を計算する。
【0033】
次に、ステップ6において、ステップ5で求めた値が式(2−1)および式(2−2)を充足するかどうかを判定する。破壊判定値が破壊基準値より大きいか等しい場合には、ステップ7においてその要素を削除し、破壊判定値が破壊基準値より小さい場合にはそのままで、ステップ4に戻り、未計算の要素が残っている限り、ステップ5において次の要素に関する応力と延性破壊条件式の破壊判定値と破壊基準値を計算する。このようにして、終了時間までステップ3〜ステップ7を繰り返す。
【0034】
図4は、破壊した要素の周囲を破壊しやすくする場合のフロー図である。この方法では、ステップ5でその要素について応力と延性破壊条件式の破壊判定値と破壊基準値を計算した後、ステップ6でその要素と接する要素が破壊しているかどうかを判定する。そして、破壊している場合は、ステップ7において破壊判定値に重み付けを行い、破壊していない場合はそのままで、ステップ8において破壊条件式を満たすかどうかを判定する。そして、破壊判定値が破壊基準値より大きいか等しい場合には、ステップ9においてその要素を削除する。その他の工程は、先の場合と同様である。
【0035】
基板部の1つのソリッド要素が破壊したらその要素の基板部の厚み方向にある全てのソリッド要素を削除する場合は、ステップ8において破壊条件式を満たさない場合において、図に破線で示すように、ステップ10で基板部の同じ厚み方向にある要素が破壊しているかを判定し、破壊している場合はステップ9に戻ってその要素を削除するようにする。
【実施例】
【0036】
樹脂成形体の素材として、住友化学株式会社製ノーブレンAZ864E4を用いた場合の衝撃特性および破壊挙動を評価する衝撃解析方法を以下に示す。
(1)解析に用いたソフトウェア
解析ソフト:LS-DYNA version9.70(Livermore
Software Technology Corporation製)
(2)解析方法
空間の離散化:有限要素法
時間積分:中心差分に基づく陽解法
材料モデル:LS-DYNA物性タイプ44(ユーザー定義物性)
(3)物性値
AZ864E4の物性値としては以下のような値を用いた。
弾性率:1000MPa
降伏応力:22.8MPa
ポアソン比:0.40
比重:0.90
真応力−真塑性歪み曲線:図5
降伏応力−歪み速度:図6
歪み硬化度−応力多軸度:図7
延性破壊条件式の材料定数:C1=0.970、C2=67.2、
C3=−0.128、C4=0.168、
C5=1.51×1014、C6=−1.65
なお、延性破壊条件式(2−1)、(2−2)の材料定数は、上記の材料定数の決定方法の一例で示したように、決めたものである。
【0037】
解析モデルの要素が基板部の割れの伝播を再現するには大きいと考え、破壊した要素の周囲は破壊が起こりやすくした。
重み係数は、W=3.0とした。
また、基板部のソリッド要素が破壊すると、その要素の基板部の厚み方向にある全てのソリッド要素も削除するようにして計算を行った。
【0038】
図1および図2に示す解析対象モデルの成形体10を実際に作成し、衝突速度5.7m/sと4.5m/sにおいて実地試験を行った。次に図2に示す解析モデルに対して上記の物性値を樹脂成形体10に用いて衝突速度5.7m/sと4.5m/sにおいて、本発明の衝撃解析方法で解析を行った(実施例)。具体的には、歪み速度がC5以下の場合は式(2−1)を、C5を超える場合は式(2−2)を用いた。
【数9】
【0039】
また、比較例として、図2に示す解析モデルに対して上記の物性値を樹脂成形体10に用いて衝突速度5.7m/sと4.5m/sにおいて、LS-DYNA物性タイプ24(多直線近似等方弾塑性体)で、解析を行った。物性タイプ24では降伏応力の歪み速度依存性を考慮することができ、破断塑性歪みは歪み速度に寄らず一定であるモデルである。破断塑性歪みのデータとしては比較例1では0.56とし、比較例2では0.57とした。
【0040】
実地試験、実施例、比較例1および2の結果である荷重−変位曲線を図8に、基板部要素の破壊の有無について表1に示す。図8において(a)は、速度5.7m/sの結果、(b)は速度4.5m/sの結果である。
解析対象モデルの成形体10の実施例、比較例1および2の衝撃解析結果後の形状を図9(a)、(b)および(c)に示す。
速度5.7m/sおよび速度4.5m/sにおける基板部の破壊状況を図10および図11に示す。図10,図11において(a)は実地試験の結果、(b)は実施例の方法による解析結果、(c)は比較例1の結果、(d)は比較例2の結果である。
【0041】
【表1】
【0042】
これらの結果より、比較例の方法では、速度によって基板部が破壊するか破壊しないかの判定を行うことができなかったが、実施例の方法では判定を行うことができた。また、図8の荷重−変位曲線もよく一致していることが分かる。
【図面の簡単な説明】
【0043】
【図1】衝撃解析を行う対象の樹脂成形体である。
【図2】衝撃試験の構成の概略を示す概念図である。
【図3】この実施の形態の衝撃解析方法の工程を示すフロー図である。
【図4】この実施の形態の衝撃解析方法の工程を示すフロー図である。
【図5】樹脂成形体の供試材料の真応力−真塑性歪み曲線である。
【図6】樹脂成形体の供試材料の降伏応力−歪み速度曲線である。
【図7】樹脂成形体の供試材料の歪み硬化度−応力多軸度曲線である。
【図8】樹脂成形体の衝撃特性を評価するための荷重−変位曲線である。
【図9】樹脂成形体の破壊状況を示す図である。
【図10】速度5.7m/sでの衝撃試験における樹脂成形体の基板部の破壊状況を示す図である。
【図11】速度4.5m/sでの衝撃試験における樹脂成形体の基板部の破壊状況を示す図である。
【符号の説明】
【0044】
10 樹脂成形体
12 リブ部
14 基板部
16 バリア(固体障壁)
18 ダート
【特許請求の範囲】
【請求項1】
解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、
樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、
樹脂成形体の破壊を判定するための破壊判定法として延性破壊条件式を用いて破壊判定を行うことを特徴とする衝撃解析方法。
【請求項2】
前記延性破壊条件式に応力多軸度を考慮することができる延性破壊条件式を用いたことを特徴とする請求項1に記載の衝撃解析方法。
【請求項3】
前記延性破壊条件式に歪み速度依存性を考慮することを特徴とする請求項1又は請求項2に記載の衝撃解析方法。
【請求項4】
前記樹脂成形体はリブ部および基板部を有し、該基板部については延性破壊条件式を用いて破壊判定を行い、前記リブ部については異なる破壊判定法を用いることを特徴とする請求項1ないし請求項3のいずれかに記載の衝撃解析方法。
【請求項5】
前記リブ部の破壊判定の際に、所定形状の樹脂成形体を用いた実地試験と、破断塑性歪みを変えた衝撃解析を複数行い、実地試験結果と衝撃解析結果を比較して決定した破断塑性歪みを用いることを特徴とする請求項4に記載の衝撃解析方法。
【請求項6】
前記リブ部の破壊判定の際に、種々の引張速度での一軸引張試験を行うことで得られる破断歪みと歪み速度との関係を用いることを特徴とする請求項4に記載の衝撃解析方法。
【請求項7】
解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、
樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、
有限要素モデルの各要素について、延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値を算出する第1の工程と、
第1の工程で求めた延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値とを比較して破壊判定を行う第2の工程を有することを特徴とする衝撃解析方法。
【請求項1】
解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、
樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、
樹脂成形体の破壊を判定するための破壊判定法として延性破壊条件式を用いて破壊判定を行うことを特徴とする衝撃解析方法。
【請求項2】
前記延性破壊条件式に応力多軸度を考慮することができる延性破壊条件式を用いたことを特徴とする請求項1に記載の衝撃解析方法。
【請求項3】
前記延性破壊条件式に歪み速度依存性を考慮することを特徴とする請求項1又は請求項2に記載の衝撃解析方法。
【請求項4】
前記樹脂成形体はリブ部および基板部を有し、該基板部については延性破壊条件式を用いて破壊判定を行い、前記リブ部については異なる破壊判定法を用いることを特徴とする請求項1ないし請求項3のいずれかに記載の衝撃解析方法。
【請求項5】
前記リブ部の破壊判定の際に、所定形状の樹脂成形体を用いた実地試験と、破断塑性歪みを変えた衝撃解析を複数行い、実地試験結果と衝撃解析結果を比較して決定した破断塑性歪みを用いることを特徴とする請求項4に記載の衝撃解析方法。
【請求項6】
前記リブ部の破壊判定の際に、種々の引張速度での一軸引張試験を行うことで得られる破断歪みと歪み速度との関係を用いることを特徴とする請求項4に記載の衝撃解析方法。
【請求項7】
解析プログラムに樹脂材料の機械的特性を示す物性値を設定し、
樹脂成形体の衝撃特性および破壊挙動を評価する衝撃解析方法において、
有限要素モデルの各要素について、延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値を算出する第1の工程と、
第1の工程で求めた延性破壊条件式に関する破壊判定値と、延性破壊条件式に関する破壊判定基準値とを比較して破壊判定を行う第2の工程を有することを特徴とする衝撃解析方法。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【公開番号】特開2006−308569(P2006−308569A)
【公開日】平成18年11月9日(2006.11.9)
【国際特許分類】
【出願番号】特願2006−85267(P2006−85267)
【出願日】平成18年3月27日(2006.3.27)
【出願人】(000002093)住友化学株式会社 (8,981)
【Fターム(参考)】
【公開日】平成18年11月9日(2006.11.9)
【国際特許分類】
【出願日】平成18年3月27日(2006.3.27)
【出願人】(000002093)住友化学株式会社 (8,981)
【Fターム(参考)】
[ Back to top ]