JP3664806B2 - 高分子多相系材料の衝撃強度計測装置及び方法 - Google Patents
高分子多相系材料の衝撃強度計測装置及び方法 Download PDFInfo
- Publication number
- JP3664806B2 JP3664806B2 JP12615496A JP12615496A JP3664806B2 JP 3664806 B2 JP3664806 B2 JP 3664806B2 JP 12615496 A JP12615496 A JP 12615496A JP 12615496 A JP12615496 A JP 12615496A JP 3664806 B2 JP3664806 B2 JP 3664806B2
- Authority
- JP
- Japan
- Prior art keywords
- strain
- energy
- displacement
- release rate
- stress
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000000463 material Substances 0.000 title claims description 70
- 238000000034 method Methods 0.000 title claims description 65
- 229920000642 polymer Polymers 0.000 title claims description 45
- 238000006073 displacement reaction Methods 0.000 claims description 67
- 238000009863 impact test Methods 0.000 claims description 41
- 238000012360 testing method Methods 0.000 claims description 36
- 238000004364 calculation method Methods 0.000 claims description 31
- 229920001971 elastomer Polymers 0.000 claims description 28
- 239000002245 particle Substances 0.000 claims description 25
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 230000001419 dependent effect Effects 0.000 claims description 5
- 229920006127 amorphous resin Polymers 0.000 claims description 4
- 239000000470 constituent Substances 0.000 claims description 3
- 239000000806 elastomer Substances 0.000 claims description 3
- 229920005669 high impact polystyrene Polymers 0.000 description 25
- 239000004793 Polystyrene Substances 0.000 description 24
- 239000004797 high-impact polystyrene Substances 0.000 description 24
- 229920002223 polystyrene Polymers 0.000 description 24
- 239000005060 rubber Substances 0.000 description 20
- 238000010586 diagram Methods 0.000 description 14
- 238000012545 processing Methods 0.000 description 14
- 239000011342 resin composition Substances 0.000 description 10
- 239000006185 dispersion Substances 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 229920005990 polystyrene resin Polymers 0.000 description 7
- 239000005062 Polybutadiene Substances 0.000 description 6
- PPBRXRYQALVLMV-UHFFFAOYSA-N Styrene Chemical group C=CC1=CC=CC=C1 PPBRXRYQALVLMV-UHFFFAOYSA-N 0.000 description 6
- 239000002131 composite material Substances 0.000 description 6
- 229920002857 polybutadiene Polymers 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 150000001875 compounds Chemical class 0.000 description 4
- 239000002861 polymer material Substances 0.000 description 4
- 230000000379 polymerizing effect Effects 0.000 description 4
- 229920001577 copolymer Polymers 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 229920000098 polyolefin Polymers 0.000 description 3
- 239000012779 reinforcing material Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 229920005992 thermoplastic resin Polymers 0.000 description 3
- KAKZBPTYRLMSJV-UHFFFAOYSA-N Butadiene Chemical compound C=CC=C KAKZBPTYRLMSJV-UHFFFAOYSA-N 0.000 description 2
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 2
- 239000004952 Polyamide Substances 0.000 description 2
- 239000004721 Polyphenylene oxide Substances 0.000 description 2
- 229920005601 base polymer Polymers 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001125 extrusion Methods 0.000 description 2
- 239000010419 fine particle Substances 0.000 description 2
- 238000004898 kneading Methods 0.000 description 2
- 229920002647 polyamide Polymers 0.000 description 2
- 229920006380 polyphenylene oxide Polymers 0.000 description 2
- 229920005989 resin Polymers 0.000 description 2
- 239000011347 resin Substances 0.000 description 2
- 238000003756 stirring Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 229920002943 EPDM rubber Polymers 0.000 description 1
- VVQNEPGJFQJSBK-UHFFFAOYSA-N Methyl methacrylate Chemical compound COC(=O)C(C)=C VVQNEPGJFQJSBK-UHFFFAOYSA-N 0.000 description 1
- 229920001328 Polyvinylidene chloride Polymers 0.000 description 1
- GWEVSGVZZGPLCZ-UHFFFAOYSA-N Titan oxide Chemical compound O=[Ti]=O GWEVSGVZZGPLCZ-UHFFFAOYSA-N 0.000 description 1
- 229920001893 acrylonitrile styrene Polymers 0.000 description 1
- 229910000019 calcium carbonate Inorganic materials 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000006229 carbon black Substances 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000635 electron micrograph Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 239000003063 flame retardant Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000001746 injection moulding Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 229920002521 macromolecule Polymers 0.000 description 1
- 239000002480 mineral oil Substances 0.000 description 1
- 235000010446 mineral oil Nutrition 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000004014 plasticizer Substances 0.000 description 1
- 229920000728 polyester Polymers 0.000 description 1
- 229920006149 polyester-amide block copolymer Polymers 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 239000004800 polyvinyl chloride Substances 0.000 description 1
- 229920000915 polyvinyl chloride Polymers 0.000 description 1
- 239000005033 polyvinylidene chloride Substances 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- SCUZVMOVTVSBLE-UHFFFAOYSA-N prop-2-enenitrile;styrene Chemical compound C=CC#N.C=CC1=CC=CC=C1 SCUZVMOVTVSBLE-UHFFFAOYSA-N 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 229920003048 styrene butadiene rubber Polymers 0.000 description 1
- OGIDPMRJRNCKJF-UHFFFAOYSA-N titanium oxide Inorganic materials [Ti]=O OGIDPMRJRNCKJF-UHFFFAOYSA-N 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、複数の成分で一定の３次元単位構造が繰り返される高分子多相系材料即ち所謂複合材料の衝撃試験値を演算によって計測する高分子多相系材料の衝撃強度測定装置及び方法に関する。
【０００２】
【従来の技術】
複数の成分からなる高分子材料とは、高分子化合物の連続相（母材）に繊維、粒子等の連続相を形成する高分子化合物とは異なる性質の強化材が分散、積層した構造体である。
【０００３】
連続層が熱可塑性樹脂で、分散層が有機・無機微粒子又は繊維の組合わせの例として、フィラー強化ポリオレフィンが、連続層が弾性体で分散層が無機微粒子の組合わせの例としてカーボン強化ゴム組成物が、連続層が熱可塑性樹脂で、分散層が弾性体粒子の組合わせとしてゴム強化樹脂組成物が知られている。
【０００４】
これらの複数の成分からなる高分子材料に共通する特徴としては、材料の機械的特性が単に材料を形成する成分単独の性質の和或いは平均によって決まるものではなく、両成分の形状、分散状態といった各成分の分散形態によってて大きく変化する点がある。
【０００５】
したがって、複数の成分からなる高分子材料を開発する場合には、これら各成分単独の性質を調査し、組合わせるだけでは不十分であり、分散状態を変化させた実験を繰り返しながら最適の組合わせと分散状態を見つけ出すといった複雑な手間の係る作業が必要であった。
【０００６】
従来、斯かる煩雑さをさけて最適点を効率的に探し出す試みが、学術的にも工学的にも広く行われていることは周知の通りである。例えば、ゴム強化熱可塑性樹脂の代表的材料であるハイインパクトポリスチレンの定性的設定指標に関するエヒテの論文（Makromol.Chem.58 175 (1977))（以下、第１従来例と称す）が有名である。
【０００７】
また、分散相を応力の集中点に見立ててサイズの異なる複数の粒子周辺の応力を計算し高分子多相系材料の破壊機構を解析する手法（以下、第２従来例と称す）も提案されている（Macromolecules 24-20 5689 (1991)) 。
【０００８】
【発明が解決しようとする課題】
しかしながら、上記第１従来例にあっては、ハイインパクトポリスチレンの定性的設定指標については最適化を行うことができるが、この論文以降にも数多くの論文等が発表されていることから最適化に関係する問題が完全には解決していないことが明らかである。この背景には毎年数多くの高分子化合物が登場し、組合わせられる化合物の数が膨大になったこと、重合技術・押し出し混練技術などが進歩し、形状・分散状態制御の自由度が格段に向上したことが挙げられる。
【０００９】
また、第２従来例にあっては、高分子多相系材料内部の応力分布を推定するには便利であるが、得られる情報は静的な範囲に留まり、破壊という動的な現象を解析し材料強度を推定するには不十分であるという未解決の課題がある。
【００１０】
そこで、本発明は、高分子多相系材料の衝撃強度を衝撃試験機を用いることなく容易に計測することができる高分子多相系材料の衝撃強度計測装置及び方法を提供することを目的としている。
【００１１】
【課題を解決するための手段】
上記目的を達成するために、請求項１に係る発明は、連続相と分散相とが一定の３次元単位構造が繰り返される高分子多相系材料の衝撃強度を計測する高分子多相系材料の衝撃強度計測装置であって、前記連続相及び分散相の構成材料毎のポアソン比及び限界歪みエネルギ開放率、歪み速度依存データとしての応力歪み特性データ等の材料データと、３次元座標で与えられる単位形状データ、及び所定の境界条件を入力して、予め設定した試験モデルを有限要素に分割した各要素について、有限要素法に従って変位を与えたときの歪みエネルギ開放率を算出し、当該歪みエネルギ開放率が限界歪みエネルギ開放率に達するまで前記変位を増加させることにより、前記高分子多相系材料の歪み速度の異なった応力歪み特性を演算する応力歪み特性演算手段と、該応力歪み特性演算手段で演算した応力歪み特性に基づいて衝撃試験モデルを分割した各要素が破断状態に達するまで前記衝撃試験変位を増加させて、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するエネルギ演算手段と、該エネルギ演算手段で算出した各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギに基づいて衝撃値を演算する衝撃値演算手段とを備えたことを特徴としている。
【００１２】
また、請求項２に係る発明は、上記請求項１の発明において、前記応力歪み特性演算手段は、試験モデルに対する変位の増分を設定する変位設定手段と、該変位設定手段で設定された変位を与えたときの各要素の歪みエネルギと外部エネルギとを演算し、両エネルギと各要素の断面積とに基づいて歪みエネルギ開放率を演算する歪みエネルギ開放率演算手段と、該歪みエネルギ開放率演算手段の歪みエネルギ開放率が限界歪みエネルギ開放率に達したか否かを判定し、限界歪みエネルギ開放率に達していないときに前記変位設定手段の変位を増加させ、限界歪みエネルギ開放率に達したときに該当要素の変位、応力、沈み、歪みエネルギ、外部エネルギ、反力及び歪みエネルギ開放率を決定する特性決定手段と、該特性決定手段で決定された反力及び衝撃試験変位に基づいて公称応力及び公称歪みを演算する公称応力歪み算出手段と、該応力歪み演算手段の演算結果に基づいて高分子多相系材料の歪み速度の異なった応力歪み特性を算出する応力歪み特性算出手段とを備えていることを特徴としている。
【００１３】
さらに、請求項３に係る発明は、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みを算出し、当該歪みが破断歪みに達したときに該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されていることを特徴としている。
【００１４】
さらにまた、請求項４に係る発明は、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みエネルギと外部エネルギを算出し、両エネルギと要素断面積とで算出される歪みエネルギ開放率が限界歪みエネルギ開放率に達したときに該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されていることを特徴としている。
【００１５】
なおさらに、請求項５に係る発明は、請求項１乃至４の発明において、前記高分子多相系材料は、連続相が非晶性樹脂で、且つ分散相がエラストマー粒子であることを特徴としている。
【００１６】
また、請求項６に係る発明は、請求項５の発明において、前記材料データが複数の変形速度で求めた各連続相及び分散相の実測値を使用することを特徴としている。
【００１７】
さらに、請求項７に係る発明は、一定の３次元単位構造が繰り返される高分子多相系材料の衝撃強度を計測する高分子多相系材料の衝撃強度計測方法であって、計測対象となる高分子多相系材料のポアソン比、歪み速度依存データとしての応力歪み特性データ等の材料データと、３次元座標で与えられる単位形状データ、及び所定の境界条件とに基づいて、予め設定した試験モデルを有限要素に分割した各要素について、有限要素法に従って変位を与えたときの歪みエネルギ開放率を算出しながら、当該歪みエネルギ開放率が限界歪みエネルギに達するまで前記変位を増加させることにより、高分子多相系材料の歪み速度の異なった応力歪み特性を演算するステップと、演算した応力歪み特性に基づいて衝撃試験変位を与えたときの各要素が破断状態に達するまで前記衝撃試験変位を増加させながら、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するステップと、積算された各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギを前記切欠き部の最小断面積で除して衝撃値を演算することを特徴としている。
【００１８】
【発明の実施の形態】
以下、本発明の実施の形態を図面に基づいて説明する。
図１は、本発明の一実施形態を示すブロック図である。図中、１はデータ入力手段としてのキーボード２、データ格納手段としてのフロッピィデスク装置３、データ表示手段としてのＣＲＴディスプレイ，液晶ディスプレイ等の表示装置４及びデータ出力手段としてのプリンタ５が接続されたコンピュータであって、このコンピュータ１でキーボード２或いはフロッピィディスク装置３から入力される入力される衝撃強度を計測する高分子多相系材料の入力データに基づいて図２〜図４の演算処理を行うことにより、衝撃強度を算出する。
【００１９】
ここで、高分子多相系材料としては、ベースポリマーを汎用ポリスチレンとし、強化材をゴムとした２成分複合材料でなる高衝撃性ポリスチレン（High Impact PolyStyrene)（以下、ＨＩＰＳと略称する）を適用した。
【００２０】
このＨＩＰＳの微視的形態は、図５で模式的に示すようにベースポリマーの汎用ポリスチレンの海のなかにゴムが島状に浮かんでいる、海島構造の形態をとる。すなわち、繊維をゴムに置き換えれば、複合材料と何ら変わりなく、複合材料で既に行われているような素材設計が可能ではないかと推察される。しかし、現在、実用的な力学的素材設計は、その殆どが弾性域での設計であり、ゴムのような非弾性体への実用的な適用例は殆どない。
【００２１】
このＨＩＰＳの強化材であるゴムは厳密には同一の形状、均一分布ではないが、電子顕微鏡による映像を画像処理することにより平均粒径を算出すると、ほぼ重量分立に近似することから、ゴムの形態を球とし且つ島状のゴムに内方される汎用ポリスチレンやゴムでなるオクルードを無視し、ほぼ均一な分布を示すものと仮定する。
【００２２】
そこで、数値解析用モデルとしては、２次元的に簡略化して表すと、図６（ａ）に示す電子顕微鏡による映像を、図６（ｂ）に示すように、汎用ポリスチレンの海にゴムが均一分散しているものと仮定し、各分散領域について図６（ｃ）に示すように幾つかのユニットが連続して存在する系に置換し、次いで図６（ｄ）に示すように連続するユニットとの相互作用は、連続的に存在することを境界条件として与えることで考慮する。
【００２３】
実際の数値解析モデルは、図７に示すように、汎用ポリスチレンの海内に球状の島が存在する３次元モデルとして適用する。
また、衝撃強度としてアイゾット衝撃値を計測するものとし、このアイゾット衝撃値を求めるモデルとして、図８に示すような中央部に切欠き部１０を有する規格化されている樹脂のアイゾット試験片１１の形状をモデル化し、これにハンマ１２を衝接させるようにしている。このハンマは初速３ｍ／秒の速度で試験片に衝突し、試験片から離れるまで（試験片１１が破断するか、極端に変形するまで）を解析する。このモデルの破壊計算に大して適当な分割要素数は１００〜１０００００に選定され、好ましくは１０００〜１００００に選定する。この理由は、分割要素数が１００要素未満では計算精度が充分でなく、１０万を越えると計算に長時間を要し無意味となるからである。
【００２４】
ここで、ＨＩＰＳの衝撃値はクレーズと呼ばれる微視的亀裂によって支配される。このクレーズは、ＨＩＰＳに外荷重を与えた際、ゴムの周囲の汎用ポリスチレンに生じる微少な損傷（亀裂）であり、クレーズが多く発生すれば衝撃値は高くなり、逆に少であれば衝撃値は低下する。本実施形態においては、クレーズ現象を機械的な亀裂として取扱い、微少な欠陥が試験片内部に存在すると仮定することにより、工学的には破壊力学による取扱が可能となる。
【００２５】
そして、ＨＩＰＳの衝撃強度の解析を行う前に、図９に示すように、高分子多相系材料を構成する汎用ポリスチレン及びゴムの異なる歪み速度例えば０．０１２５ｍ／sec ，２ｍ／sec ，６ｍ／sec の応力歪み線図を実験により求めておき、この実験データを例えばフロッピィディスクに格納しておく。この図９のゴムの破断伸びは７以上（７００％以上）であり、ここでは見易くするために歪み０．１までを示す。
【００２６】
また、汎用ポリスチレン及びゴムのポアソン比も求めておく、このポアソン比は、文献的に公知の値を使用しても実験的に求めても、或いは衝撃強度の実験値と計算値のずれを最少にする値として推定して用いても構わない。
【００２７】
そして、コンピュータで図２〜図４に示す演算処理を実行する。
この演算処理は、先ず、ステップＳ１で材料データ、単位形状データ、境界条件を入力する。
【００２８】
ここで、材料データとしては応力歪み線図を格納したフロッピィディスクをフロッピィディスク装置３にセットして、これに格納されている歪み速度の異なる応力歪み線図を入力すると共に、キーボード２から汎用ポリスチレン及びゴムのポアソン比を入力し、さらに実験により求めた材料物性値である限界エネルギ開放率Ｇc を入力する。
【００２９】
また、単位形状データとしては、図７に示す体積分率から算定される形態であって画像処理によって得られる形態を３次元座標で入力する。
さらに、境界条件としては、周期対称条件、引張方向への強制変位を入力する。
【００３０】
次いで、ステップＳ２に移行して、入力された各材料データ、単位形状データ及び境界条件をもとに非線形有限要素法に従って演算を実行して、ＨＩＰＳとしての衝撃試験変位、応力、歪みエネルギ、反力、歪みエネルギ開放率を算出する。
【００３１】
この算出処理の具体例は、図３に示すように、先ず、ステップＳ２ａで初期値として衝撃試験変位即ちハンマによって与えられる変位ｘを設定し、次いでステップＳ２ｂに移行して、対称とする分割要素の１つについて変位ｘを与えたときの歪みエネルギＥstと外部エネルギＥout を算出し、次いでステップＳ２ｃに移行して、該当要素の断面積Ａを算出してからステップＳ２ｄに移行する。
【００３２】
このステップＳ２ｄでは、算出した歪みエネルギＥst及び外部エネルギＥout と断面積Ａとを基に下記（１）式の演算を行って歪みエネルギ開放率Ｇを算出する。
【００３３】
Ｇ＝Ｅout ／Ａ−Ｅst／Ａ …………（１）
次いで、ステップＳ２ｅに移行して、算出した歪みエネルギ開放率Ｇが入力された限界歪み開放率Ｇc に達したか否かを判定し、Ｇ＜Ｇc であるときには限界歪みエネルギ開放率Ｇc に達していないものと判断してステップＳ２ｆに移行して、現在の変位ｘに所定の増加分Δｘを加算した値（ｘ＋Δｘ）を新たな変位ｘとして設定してから前記ステップＳ２ｂに戻り、Ｇ＝Ｇc であるときには限界歪みエネルギ開放率Ｇc に達したものと判断してステップＳ２ｇに移行する。
【００３４】
ここで、歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達した時には、試験片モデルは変形によるエネルギ消費（歪みエネルギ）が不可能となり、破断する状態となったことを判断することができる。
【００３５】
ステップＳ２ｇでは、対象とする有限要素の全ての歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達したか否かを判定し、限界歪みエネルギ開放率Ｇc に達していない有限要素が存在するときには、前記ステップＳ２ａに戻り、限界歪みエネルギ開放率Ｇc に達していない有限要素が存在しないときには処理を終了して、処理結果であるＨＩＰＳの変位、応力、歪み、歪みエネルギ、外部エネルギ、反力及び歪みエネルギ開放率を例えば内蔵する記憶装置の所定記憶領域に格納してから図２のステップＳ３に移行する。
【００３６】
このステップＳ３では、前記ステップＳ２で格納されたＨＩＰＳの反力を変位方向に対して直角の断面積で除して単位形状としての公称応力を算出すると共に、強制変位を変位方向の単位形状長さで除して単位形状としての公称歪みを算出する。
【００３７】
次いで、ステップＳ４に移行して、ステップＳ３で算出したＨＩＰＳの公称応力及び公称歪みをもとに歪み速度の異なったＨＩＰＳの応力歪み特性を算出し、これを表す応力歪み特性線図を作成し、これを記憶装置の所定記憶領域に格納してからステップＳ５に移行する。なお、このステップＳ４の処理は市販されている非線型有限要素法解法プログラム「ＡＢＡＱＵＳ」（商品名）（米国ＨＫＳ社製）を適用することができる。
【００３８】
このステップＳ５では、ステップＳ５で格納された応力歪み特性線図をもとに衝撃試験方法に応じた試験形状データ及び試験方法に応じた境界条件を入力することにより、各有限要素について破断が生じるまでの外部エネルギ、内部エネルギ及び運動エネルギを算出し、これらを試験片全体で積算し、積算結果を記憶装置の所定記憶領域に記憶するエネルギ演算処理を実行する。
【００３９】
このエネルギ演算処理の具体例は、図４に示すように、先ず、ステップＳ５ａで１つの有限要素について初期変位を設定し、次いでステップＳ５ｂに移行して要素各々の歪みεを算出してからステップＳ５ｃに移行する。
【００４０】
このステップＳ５ｃでは、ステップＳ５ｂで算出した歪みεが前述したステップＳ４の応力歪み特性線図から求まる歪み速度に応じた破断歪みεb に達したか否かを判定し、ε＜εb であるときには、破断歪みεb に達していないものと判断して、ステップＳ５ｄに移行して、変位ｘに変位増加分Δｘを加算した値を新たな変位ｘとして設定してから前記ステップＳ５ｂに戻り、ε＝εb であるときには、破断歪みεb に達して、切欠き部から生じる亀裂（クレーズではない巨視的亀裂＝破断）が生じたものと判断してステップＳ５ｅに移行する。
【００４１】
このステップＳ５ｅでは、破断歪みεb に達した有限要素を削除し、次いでステップＳ５ｆに移行して、該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出して、試験片全体での積算を行い、その積算結果を記憶装置の所定記憶領域に更新記憶する。
【００４２】
次いで、ステップＳ５ｇに移行して、対象とする有限要素の全てが破断歪みεb に達したか否かを判定し、破断歪みεb に達していない有限要素即ち削除されていない有限要素について前記ステップＳ５ａ〜Ｓ５ｆの処理を繰り返し、全ての有限要素が破断歪みεb に達したときには破断状態に達したと判断して対象となる要素全てを削除して処理を終了し、図２のステップＳ６に移行する。
【００４３】
このステップＳ６では、記憶装置の所定領域に更新記憶されている試験片全体での外部エネルギ、内部エネルギ及び運動エネルギから試験片に蓄積された吸収エネルギを算出する。
【００４４】
次いで、ステップＳ７に移行して、吸収エネルギを切欠き部の最少断面積で除することにより、アイゾット衝撃値を算出し、これを記憶装置の所定記憶領域に記憶すると共に、表示装置に表示し、必要に応じてプリンタで印字してから処理を終了する。
【００４５】
ここで、図２の処理において、ステップＳ２〜Ｓ４の処理が応力歪み特性演算手段に対応すると共に、ステップＳ２の処理が特性決定手段に対応し、ステップＳ３の処理が公称応力歪み算出手段に対応し、ステップＳ４の処理が応力歪み特性算出手段に対応し、ステップＳ５の処理がエネルギ演算手段に対応し、ステップＳ６及びＳ７の処理が衝撃値演算手段に対応しており、またステップＳ６及びＳ７の処理は市販されている衝撃強度解析プログラム「ＲＡＤＩＯＳＳ」（商品名）（仏ＭＥＣＡＬＯＧ社製）を適用することができる。
【００４６】
このように、図２の処理によれば、高分子多相系材料を構成する連続相及び分散相成分毎のポアソン比、歪み速度を異ならせた応力歪み特性及び限界歪みエネルギ開放率Ｇc とを含む材料データと、ＨＩＰＳの形態を表す単位形状データと、境界条件を入力することにより、複合材料としてのＨＩＰＳの応力歪み特性を演算し、この演算結果から衝撃試験方法に応じた試験形状データと試験方法に応じた境界条件を入力することにより、試験片全体の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これらをエネルギに基づいて試験片に蓄積された吸収エネルギを算出し、これを試験片の切欠き部の最少面積で除することにより、仮想のＨＩＰＳであってもアイゾット衝撃値を算出することができる。
〔実施例〕
下記のように作成したゴム強化ポリスチレン樹脂組成物及びポリスチレンを適宜組み合わせて１１種類の試験片ＴＰ1 〜ＴＰ11を作成し、これをアイゾット衝撃試験機で実測した実測衝撃値と図２の演算処理によって求めた演算衝撃値との比較結果を下記表１に表す。
【００４７】
▲１▼ポリブタジエン（宇部興産製：１３ＨＢ）を溶解したスチレンモノマーを攪拌下で重合する際に、攪拌回転数を変化させ、０．６８ミクロン，１．２ミクロン及び２．７ミクロンの平均粒子径を持ち、ゴム質重合体を１１重量％含有するゴム強化ポリスチレン組成物(I),(II)及び(III) を作成した。
【００４８】
▲２▼ポリブタジエン（旭化成工業製：ＮＦ３５）を溶解したスチレンモノマーを攪拌下で重合し、１．７ミクロンの平均粒子径を持ち、ゴム質重合体を１１重量％含有するゴム強化ポリスチレン樹脂組成物(IV)を作成した。
【００４９】
▲３▼スチレン−ブタジエンブロック共重合体をスチレンモノマーに溶解し攪拌下で重合することにより、ゴム質重合体粒子内部に単一のポリスチレンオクルージョンを有する平均粒子径が０．７ミクロンのゴム質重合体粒子１１重量％含有するゴム強化ポリスチレン樹脂組成物(V) を作成した。
【００５０】
▲４▼スチレンモノマーを重合することにより、樹脂組成物(I) 〜(V) の連続相ポリスチレンと同等の分子量を有するポリスチレン(VI)を作成すると共に、ポリスチレン(VI)に対してミネラルオイル３％を混合し、ポリスチレン(VII) を作成した。
【００５１】
これらの樹脂組成物を単独或いは押し出し混練により適宜組み合わせて試験片ＴＰ1 〜ＴＰ11とし、ＡＳＴＭ Ｄ２５６規格に基づき射出成形法１／４インチ試験片ＩＺＯＤ衝撃強度を測定した。
【００５２】
また、図２の演算処理においては、事前にゴム強化ポリスチレン樹脂組成物の原料に用いたポリブタジエン、ブタジエンスラチレン共重合体、及び別途作成した連続相を形成するポリスチレン相当品を圧縮成形によって１ミリメートル厚さのシートとし、これにより試験片を切り出し、これらを０．０１２５ｍ／sec 、２．０ｍ／sec 及び６．０ｍ／sec の引張り速度で図９に示す応力歪み曲線を測定した。一方、ゴム強化ポリスチレン樹脂組成物の電子顕微鏡写真を画像解析し、繰り返し最少単位の形状を求めた。ゴム強化ポリスチレン樹脂組成物（Ｖ）に関しては、ゴム質重合体粒子内部のポリスチレン粒子を考慮し、同心円（球）内部のポリスチレンの性質は粒子外部を取り巻く連続相ポリスチレンと同一とみなした。これらの実測値を入力し、各成分のポアソン比は、ポリスチレン、ポリブタジエンの文献値を、ポリスチレンの臨界エネルギ開放率も文献値を用いた。
【００５３】
また、図７に示すノッチ（切欠き部）付きアイゾット試験片形状を、ノッチ周辺及び区ラックが進展するノッチの延長線部分は細かく、ノッチから遠ざかるに従って粗く２０００個の有限要素に分割した。
【００５４】
さらに、アイゾット試験片の規定の位置に初速毎秒３メートルでハンマが衝突する条件のもとで、試験片が破断し分離するまでにようする試験片に蓄えられるエネルギをノッチ部分の最小断面積で除してアイゾット衝撃値を求めた。
【００５５】
なおさらに、ノッチ部分から生じるクラックは、解析により求まる歪みが繰り返し単位の数値解析によって得られる速度依存の応力歪み曲線の破断歪みεb を越える時点において生じるとして取り扱った。
【００５６】
【表１】
【００５７】
この表１から明らかなように、試験片ＴＰ1 〜ＴＰ3 及びＴＰ5 〜ＴＰ7 の実測値と演算値の弾性体粒子径及び弾性体量依存性の傾向は、アイゾット衝撃値を最大にする最適粒子径の存在、弾性体量減少によりアイゾット衝撃値が低下する減少を正しく表現し、試験片ＴＰ2 及びＴＰ4 の比較からゴムの種類によって最大値が異なることを正しく表している。
【００５８】
また、試験片ＴＰ7 及びＴＰ8 の比較においても可塑材効果を表現していることが明らかである。さらに試験片ＴＰ1 〜ＴＰ2 とＴＰ9 〜ＴＰ11との比較から、粒子径の小さい単一オクルージョンを有する弾性体粒子からなるゴム強化ポリスチレンは、少量の比較的大きなサイズの一般的弾性体粒子を添加することで飛躍的にそのアイゾット衝撃強度を高める現象が表現されている。
【００５９】
一方、重量分率１２．３％でゴム粒径が０．６８μｍ、１．２５μｍ及び２．７０μｍであるときのＨＩＰＳのステップＳ４で算出される応力歪み特性値（ＦＥＭ）と実験値との関係を図１０〜図１２に示す。
【００６０】
これらの図からも明らかなように、本実施形態で求めたＨＩＰＳの応力歪み特性は、実験値とほぼ一致しており、歪みエネルギ開放率Ｇが限界歪み開放率に達したときに要素の破断状態と判断することが最適であることが立証された。
【００６１】
次に、本発明の第２の実施形態を図１３について説明する。
この第２の実施形態においては、前述した第１の実施形態における図２のステップＳ５のエネルギ演算処理が図１３に示すように変更されていることを除いては前述した図２と同様の処理を実行し、対応するステップには同一ステップ番号を付して、その詳細説明はこれを省略する。
【００６２】
すなわち、第２の実施形態では、エネルギ演算処理が、図１３に示すように、先ず、ステップＳ５ｈで１つの有限要素について初期変位ｘを設定し、次いでステップＳ５ｉに移行して、変位ｘでの歪みエネルギＥstJ 外部エネルギＥout を算出し、次いでステップＳ５ｊに移行して、該当要素の断面積Ａを算出してからステップＳ５ｋに移行する。
【００６３】
このステップＳ５ｋでは、前述した（１）式の演算を行って歪みエネルギ開放率Ｇを算出し、次いでステップＳ５ｌに移行して歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達したか否かを判定し、Ｇ＜Ｇc であるときには限界歪みエネルギ開放率Ｇc に達していないものと判断してステップＳ５ｍに移行して、現在の変位ｘに変位増加分Δｘを加算した値を新たな変位ｘとして設定してから前記ステップＳ５ｉに戻り、Ｇ＝Ｇc であるときには有限要素が破断状態に達したと判断してステップＳ５ｎに移行する。
【００６４】
このステップＳ５ｎでは、破断状態に達した要素を削除し、次いでステップＳ５ｏに移行して、削除した要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これらを試験片全体で積算し、その積算結果を記憶装置の所定記憶領域に更新記憶する。
【００６５】
次いで、ステップＳ５ｐに移行して、対象とする要素全てが限界歪みエネルギ開放率Ｇc に達したか否かを判定し、限界歪みエネルギ開放率Ｇc に達していない要素があるとき即ち削除されずに残っている要素があるときには、前記ステップＳ５ｈに戻って、残りの要素についてステップＳ５ｈ〜Ｓ５ｎの処理を行い、全ての要素が限界歪みエネルギ開放率Ｇc に達したときには、試験片が破断状態に達したものと判断して、ステップＳ５ｑに移行して、対象要素を全て削除してから図２のステップＳ６に移行する。
【００６６】
この図１１のエネルギ演算処理においても、各有限要素について変位を増加させながら各要素の歪みエネルギ開放率Ｇを演算し、この歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達したときに該当有限要素が破断状態となったものと判断して、該当有限要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを試験片全体で積算し、その積算値を記憶装置の所定記憶領域に更新記憶するので、全ての要素の解析が終了した時点で、試験片に蓄積された吸収エネルギを算出し、これを試験片のノッチ部分の最小断面積で除することにより目的とするアイゾット衝撃値を算出することができる。
【００６７】
なお、上記実施例においては、高分子多相系材料がゴム強化ポリスチレン樹脂組成物である場合について説明したが、これに限定されるものではなく、連続相を形成する物質としては、前述したポリスチレンの他、アクリロニトル−スチレン共重合体、ポリフェニレンオキサイド、ポリスチレン−ポリフェニレンオキサイド混合物といったスチレン系樹脂及びこれらに可塑材・難燃剤を添加した組成物があり、他にもポリ塩化ビニル、ポリ塩化ビニリデン、ポリメチルメタクリレート、ポリブタジエンといった非晶性樹脂及びこれらの組成物にも好適に適用することができる。
【００６８】
また、分散相を形成する物質としては、ポリブタジエン、ブタジエン−スチレン共重合体、アクリロニトリル−ブタジエン共重合体、ＥＶＡ、ＥＰＤＭといった連続相より弾性率の低い物質、ボリオレフィン、ボリアミド、ボリエステルといった結晶性高分子、炭酸カルシウム、酸化チタン、カーボンブラックといった弾性率の高い物質があげられる。これらポリオレフィン、ポリエステル、ポリアミドといった結晶性高分子材料の場合には、結晶相を分散相として、非晶相を連続相として取り扱うことで図２の演算処理の適用が可能である。
【００６９】
さらに、分散相の分布形態は、球に限らず、紡錘体、柱状、面等任意に与えることができるし、分布形状は電子顕微鏡写真を画像処理して入力しても、架空の分散形状を定義して入力してもよく、分散相と連続相の比率も任意に設定することができる。
【００７０】
なおさらに、上記実施形態においては、アイゾット衝撃強度を計測する場合について説明したが、これに限定されるものではなく、シャルピー衝撃強度やその他の衝撃強度を計測することができる。この場合にはエネルギ演算手段で衝撃試験方法に応じた衝撃試験モデル及び衝撃試験変位を与えて解析すればよい。
【００７１】
また、上記実施形態においては、衝撃試験片として、切欠き部１１を有するモデルを適用した場合について説明したが、これに限定されるものではなく、切欠き部を有さない衝撃試験モデルについても、破断位置を特定することにより、吸収エネルギに基づいて衝撃値を容易に算出することができる。
【００７２】
【発明の効果】
以上説明したように、請求項１及び請求項７に係る発明によれば、入力手段によって、高分子多相系材料を構成する複合材料毎の応力歪み特性とポアソン比でなる材料データと、高分子多相系材料の形態を表す形状データと所定の境界条件を入力することにより、非線形有限要素法に従って試料モデルを分割した有限要素毎に歪みエネルギ開放率が限界歪みエネルギ開放率に達するまで変位を増加させることにより、高分子多相系材料の破断状態に至る応力歪み特性を演算し、この演算結果から、衝撃試験変位を与えたときの各要素が破断状態に達するまで前記衝撃試験変位を増加させながら、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するステップと、積算された各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギを前記切欠き部の最小断面積で除して衝撃値を演算するように構成したので、高分子多相系材料の各構成材料の歪み速度の異なる応力歪み特性、ポアソン比及び限界歪みエネルギ開放率を与えるだけで、衝撃強度を計測することが可能となり、実際に該当する高分子多相系材料が存在しない場合でも確実に衝撃強度を計測することができ、高分子多相系材料を設計開発する場合に大いに貢献することができるという効果が得られる。
【００７３】
また、請求項２に係る発明によれば、請求項１の発明において、前記応力歪み特性演算手段は、試験モデルに対する変位の増分を設定する変位設定手段と、該変位設定手段で設定された衝撃試験変位を与えたときの各要素の歪みエネルギと外部エネルギとを演算し、両エネルギと各要素の断面積とに基づいて歪みエネルギ開放率を演算する歪みエネルギ開放率演算手段と、該歪みエネルギ開放率演算手段の歪みエネルギ開放率が限界歪みエネルギ開放率に達したか否かを判定し、限界歪みエネルギ開放率に達していないときに前記変位を増加させ、限界歪みエネルギ開放率に達したときに該当要素の変位、応力、沈み、歪みエネルギ、外部エネルギ、反力及び歪みエネルギ開放率を決定する特性決定手段と、該特性決定手段で決定された反力及び衝撃試験変位に基づいて公称応力及び公称歪みを演算するする応力歪み演算手段と、該応力歪み演算手段の演算結果に基づいて歪み速度の異なった応力歪み特性を算出する応力歪み特性算出手段とを備えているので、試験片の歪みエネルギ開放率と限界歪みエネルギ開放率とを比較し、限界歪みエネルギ開放率に達したときに試験片が歪むことによってエネルギ消費をし尽くすことができず破断する状態となったことを確実に判断することができ、構造の定義された高分子多相系材料の応力歪み特性を正確に計測することができるという効果が得られる。
【００７４】
さらに、請求項３に係る発明によれば、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みを算出し、当該歪みが破断歪みに達したときに破断状態と判断して該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されているので、各要素の破断状態を確実に判断することができるという効果が得られる。
【００７５】
さらにまた、請求項４に係る発明によれば、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みエネルギと外部エネルギを算出し、両エネルギと要素断面積とで算出される歪みエネルギ開放率が限界歪みエネルギ開放率に達したときに破綻状態と判断して該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されているので、請求項３の発明と同様に、各要素の破断状態を確実に判断することができるという効果が得られる。
【００７６】
なおさらに、請求項５に係る発明によれば、請求項１乃至４の発明において、前記高分子多相系材料は、連続相が非晶性樹脂で、且つ分散相がエラストマー粒子であるので、弾性率の大きく異なる材料の組み合わされた系の各成分の衝撃強度に対する寄与を確実に判断できるという効果が得られる。
【００７７】
また、請求項６に係る発明によれば、請求項５の発明において、前記材料データは複数の変形速度で求めた各連続相及び分散相の実測値を使用するので、材料データの入力が容易となるという効果が得られる。
【図面の簡単な説明】
【図１】本発明の第１の実施形態を示す概略構成図である。
【図２】第１の実施形態における衝撃強度演算処理手順の一例を示すフローチャートである。
【図３】図２の応力歪み特性演算処理の具体例を示すフローチャートである。
【図４】図２のエネルギ演算処理の具体例を示すフローチャートである。
【図５】ハイインパクトポリスチレンの形態を示す模式図である。
【図６】ハイインパクトポリスチレンの形態をモデル化する場合の説明図である。
【図７】ハイインパクトポリスチレンの３次元モデルを示す説明図である。
【図８】アイゾット衝撃試験片のモデルを示す説明図である。
【図９】ハイインパクトポリスチレンを構成する汎用ポリスチレン及びゴムの応力歪み特性を示す特性線図である。
【図１０】ゴム粒径が０．６８μｍであるときのＨＩＰＳの演算値と実測値との関係を示す特性線図である。
【図１１】ゴム粒径が１．２５μｍであるときのＨＩＰＳの演算値と実測値との関係を示す特性線図である。
【図１２】ゴム粒径が２．７０μｍであるときのＨＩＰＳの演算値と実測値との関係を示す特性線図である。
【図１３】本発明の第２の実施形態を示す図３に対応するエネルギ演算処理の具体例を示すフローチャートである。
【符号の説明】
１ コンピュータ
２ キーボード
３ フロッピィデスク装置
４ 表示装置
５ プリンタ
【発明の属する技術分野】
本発明は、複数の成分で一定の３次元単位構造が繰り返される高分子多相系材料即ち所謂複合材料の衝撃試験値を演算によって計測する高分子多相系材料の衝撃強度測定装置及び方法に関する。
【０００２】
【従来の技術】
複数の成分からなる高分子材料とは、高分子化合物の連続相（母材）に繊維、粒子等の連続相を形成する高分子化合物とは異なる性質の強化材が分散、積層した構造体である。
【０００３】
連続層が熱可塑性樹脂で、分散層が有機・無機微粒子又は繊維の組合わせの例として、フィラー強化ポリオレフィンが、連続層が弾性体で分散層が無機微粒子の組合わせの例としてカーボン強化ゴム組成物が、連続層が熱可塑性樹脂で、分散層が弾性体粒子の組合わせとしてゴム強化樹脂組成物が知られている。
【０００４】
これらの複数の成分からなる高分子材料に共通する特徴としては、材料の機械的特性が単に材料を形成する成分単独の性質の和或いは平均によって決まるものではなく、両成分の形状、分散状態といった各成分の分散形態によってて大きく変化する点がある。
【０００５】
したがって、複数の成分からなる高分子材料を開発する場合には、これら各成分単独の性質を調査し、組合わせるだけでは不十分であり、分散状態を変化させた実験を繰り返しながら最適の組合わせと分散状態を見つけ出すといった複雑な手間の係る作業が必要であった。
【０００６】
従来、斯かる煩雑さをさけて最適点を効率的に探し出す試みが、学術的にも工学的にも広く行われていることは周知の通りである。例えば、ゴム強化熱可塑性樹脂の代表的材料であるハイインパクトポリスチレンの定性的設定指標に関するエヒテの論文（Makromol.Chem.58 175 (1977))（以下、第１従来例と称す）が有名である。
【０００７】
また、分散相を応力の集中点に見立ててサイズの異なる複数の粒子周辺の応力を計算し高分子多相系材料の破壊機構を解析する手法（以下、第２従来例と称す）も提案されている（Macromolecules 24-20 5689 (1991)) 。
【０００８】
【発明が解決しようとする課題】
しかしながら、上記第１従来例にあっては、ハイインパクトポリスチレンの定性的設定指標については最適化を行うことができるが、この論文以降にも数多くの論文等が発表されていることから最適化に関係する問題が完全には解決していないことが明らかである。この背景には毎年数多くの高分子化合物が登場し、組合わせられる化合物の数が膨大になったこと、重合技術・押し出し混練技術などが進歩し、形状・分散状態制御の自由度が格段に向上したことが挙げられる。
【０００９】
また、第２従来例にあっては、高分子多相系材料内部の応力分布を推定するには便利であるが、得られる情報は静的な範囲に留まり、破壊という動的な現象を解析し材料強度を推定するには不十分であるという未解決の課題がある。
【００１０】
そこで、本発明は、高分子多相系材料の衝撃強度を衝撃試験機を用いることなく容易に計測することができる高分子多相系材料の衝撃強度計測装置及び方法を提供することを目的としている。
【００１１】
【課題を解決するための手段】
上記目的を達成するために、請求項１に係る発明は、連続相と分散相とが一定の３次元単位構造が繰り返される高分子多相系材料の衝撃強度を計測する高分子多相系材料の衝撃強度計測装置であって、前記連続相及び分散相の構成材料毎のポアソン比及び限界歪みエネルギ開放率、歪み速度依存データとしての応力歪み特性データ等の材料データと、３次元座標で与えられる単位形状データ、及び所定の境界条件を入力して、予め設定した試験モデルを有限要素に分割した各要素について、有限要素法に従って変位を与えたときの歪みエネルギ開放率を算出し、当該歪みエネルギ開放率が限界歪みエネルギ開放率に達するまで前記変位を増加させることにより、前記高分子多相系材料の歪み速度の異なった応力歪み特性を演算する応力歪み特性演算手段と、該応力歪み特性演算手段で演算した応力歪み特性に基づいて衝撃試験モデルを分割した各要素が破断状態に達するまで前記衝撃試験変位を増加させて、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するエネルギ演算手段と、該エネルギ演算手段で算出した各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギに基づいて衝撃値を演算する衝撃値演算手段とを備えたことを特徴としている。
【００１２】
また、請求項２に係る発明は、上記請求項１の発明において、前記応力歪み特性演算手段は、試験モデルに対する変位の増分を設定する変位設定手段と、該変位設定手段で設定された変位を与えたときの各要素の歪みエネルギと外部エネルギとを演算し、両エネルギと各要素の断面積とに基づいて歪みエネルギ開放率を演算する歪みエネルギ開放率演算手段と、該歪みエネルギ開放率演算手段の歪みエネルギ開放率が限界歪みエネルギ開放率に達したか否かを判定し、限界歪みエネルギ開放率に達していないときに前記変位設定手段の変位を増加させ、限界歪みエネルギ開放率に達したときに該当要素の変位、応力、沈み、歪みエネルギ、外部エネルギ、反力及び歪みエネルギ開放率を決定する特性決定手段と、該特性決定手段で決定された反力及び衝撃試験変位に基づいて公称応力及び公称歪みを演算する公称応力歪み算出手段と、該応力歪み演算手段の演算結果に基づいて高分子多相系材料の歪み速度の異なった応力歪み特性を算出する応力歪み特性算出手段とを備えていることを特徴としている。
【００１３】
さらに、請求項３に係る発明は、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みを算出し、当該歪みが破断歪みに達したときに該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されていることを特徴としている。
【００１４】
さらにまた、請求項４に係る発明は、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みエネルギと外部エネルギを算出し、両エネルギと要素断面積とで算出される歪みエネルギ開放率が限界歪みエネルギ開放率に達したときに該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されていることを特徴としている。
【００１５】
なおさらに、請求項５に係る発明は、請求項１乃至４の発明において、前記高分子多相系材料は、連続相が非晶性樹脂で、且つ分散相がエラストマー粒子であることを特徴としている。
【００１６】
また、請求項６に係る発明は、請求項５の発明において、前記材料データが複数の変形速度で求めた各連続相及び分散相の実測値を使用することを特徴としている。
【００１７】
さらに、請求項７に係る発明は、一定の３次元単位構造が繰り返される高分子多相系材料の衝撃強度を計測する高分子多相系材料の衝撃強度計測方法であって、計測対象となる高分子多相系材料のポアソン比、歪み速度依存データとしての応力歪み特性データ等の材料データと、３次元座標で与えられる単位形状データ、及び所定の境界条件とに基づいて、予め設定した試験モデルを有限要素に分割した各要素について、有限要素法に従って変位を与えたときの歪みエネルギ開放率を算出しながら、当該歪みエネルギ開放率が限界歪みエネルギに達するまで前記変位を増加させることにより、高分子多相系材料の歪み速度の異なった応力歪み特性を演算するステップと、演算した応力歪み特性に基づいて衝撃試験変位を与えたときの各要素が破断状態に達するまで前記衝撃試験変位を増加させながら、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するステップと、積算された各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギを前記切欠き部の最小断面積で除して衝撃値を演算することを特徴としている。
【００１８】
【発明の実施の形態】
以下、本発明の実施の形態を図面に基づいて説明する。
図１は、本発明の一実施形態を示すブロック図である。図中、１はデータ入力手段としてのキーボード２、データ格納手段としてのフロッピィデスク装置３、データ表示手段としてのＣＲＴディスプレイ，液晶ディスプレイ等の表示装置４及びデータ出力手段としてのプリンタ５が接続されたコンピュータであって、このコンピュータ１でキーボード２或いはフロッピィディスク装置３から入力される入力される衝撃強度を計測する高分子多相系材料の入力データに基づいて図２〜図４の演算処理を行うことにより、衝撃強度を算出する。
【００１９】
ここで、高分子多相系材料としては、ベースポリマーを汎用ポリスチレンとし、強化材をゴムとした２成分複合材料でなる高衝撃性ポリスチレン（High Impact PolyStyrene)（以下、ＨＩＰＳと略称する）を適用した。
【００２０】
このＨＩＰＳの微視的形態は、図５で模式的に示すようにベースポリマーの汎用ポリスチレンの海のなかにゴムが島状に浮かんでいる、海島構造の形態をとる。すなわち、繊維をゴムに置き換えれば、複合材料と何ら変わりなく、複合材料で既に行われているような素材設計が可能ではないかと推察される。しかし、現在、実用的な力学的素材設計は、その殆どが弾性域での設計であり、ゴムのような非弾性体への実用的な適用例は殆どない。
【００２１】
このＨＩＰＳの強化材であるゴムは厳密には同一の形状、均一分布ではないが、電子顕微鏡による映像を画像処理することにより平均粒径を算出すると、ほぼ重量分立に近似することから、ゴムの形態を球とし且つ島状のゴムに内方される汎用ポリスチレンやゴムでなるオクルードを無視し、ほぼ均一な分布を示すものと仮定する。
【００２２】
そこで、数値解析用モデルとしては、２次元的に簡略化して表すと、図６（ａ）に示す電子顕微鏡による映像を、図６（ｂ）に示すように、汎用ポリスチレンの海にゴムが均一分散しているものと仮定し、各分散領域について図６（ｃ）に示すように幾つかのユニットが連続して存在する系に置換し、次いで図６（ｄ）に示すように連続するユニットとの相互作用は、連続的に存在することを境界条件として与えることで考慮する。
【００２３】
実際の数値解析モデルは、図７に示すように、汎用ポリスチレンの海内に球状の島が存在する３次元モデルとして適用する。
また、衝撃強度としてアイゾット衝撃値を計測するものとし、このアイゾット衝撃値を求めるモデルとして、図８に示すような中央部に切欠き部１０を有する規格化されている樹脂のアイゾット試験片１１の形状をモデル化し、これにハンマ１２を衝接させるようにしている。このハンマは初速３ｍ／秒の速度で試験片に衝突し、試験片から離れるまで（試験片１１が破断するか、極端に変形するまで）を解析する。このモデルの破壊計算に大して適当な分割要素数は１００〜１０００００に選定され、好ましくは１０００〜１００００に選定する。この理由は、分割要素数が１００要素未満では計算精度が充分でなく、１０万を越えると計算に長時間を要し無意味となるからである。
【００２４】
ここで、ＨＩＰＳの衝撃値はクレーズと呼ばれる微視的亀裂によって支配される。このクレーズは、ＨＩＰＳに外荷重を与えた際、ゴムの周囲の汎用ポリスチレンに生じる微少な損傷（亀裂）であり、クレーズが多く発生すれば衝撃値は高くなり、逆に少であれば衝撃値は低下する。本実施形態においては、クレーズ現象を機械的な亀裂として取扱い、微少な欠陥が試験片内部に存在すると仮定することにより、工学的には破壊力学による取扱が可能となる。
【００２５】
そして、ＨＩＰＳの衝撃強度の解析を行う前に、図９に示すように、高分子多相系材料を構成する汎用ポリスチレン及びゴムの異なる歪み速度例えば０．０１２５ｍ／sec ，２ｍ／sec ，６ｍ／sec の応力歪み線図を実験により求めておき、この実験データを例えばフロッピィディスクに格納しておく。この図９のゴムの破断伸びは７以上（７００％以上）であり、ここでは見易くするために歪み０．１までを示す。
【００２６】
また、汎用ポリスチレン及びゴムのポアソン比も求めておく、このポアソン比は、文献的に公知の値を使用しても実験的に求めても、或いは衝撃強度の実験値と計算値のずれを最少にする値として推定して用いても構わない。
【００２７】
そして、コンピュータで図２〜図４に示す演算処理を実行する。
この演算処理は、先ず、ステップＳ１で材料データ、単位形状データ、境界条件を入力する。
【００２８】
ここで、材料データとしては応力歪み線図を格納したフロッピィディスクをフロッピィディスク装置３にセットして、これに格納されている歪み速度の異なる応力歪み線図を入力すると共に、キーボード２から汎用ポリスチレン及びゴムのポアソン比を入力し、さらに実験により求めた材料物性値である限界エネルギ開放率Ｇc を入力する。
【００２９】
また、単位形状データとしては、図７に示す体積分率から算定される形態であって画像処理によって得られる形態を３次元座標で入力する。
さらに、境界条件としては、周期対称条件、引張方向への強制変位を入力する。
【００３０】
次いで、ステップＳ２に移行して、入力された各材料データ、単位形状データ及び境界条件をもとに非線形有限要素法に従って演算を実行して、ＨＩＰＳとしての衝撃試験変位、応力、歪みエネルギ、反力、歪みエネルギ開放率を算出する。
【００３１】
この算出処理の具体例は、図３に示すように、先ず、ステップＳ２ａで初期値として衝撃試験変位即ちハンマによって与えられる変位ｘを設定し、次いでステップＳ２ｂに移行して、対称とする分割要素の１つについて変位ｘを与えたときの歪みエネルギＥstと外部エネルギＥout を算出し、次いでステップＳ２ｃに移行して、該当要素の断面積Ａを算出してからステップＳ２ｄに移行する。
【００３２】
このステップＳ２ｄでは、算出した歪みエネルギＥst及び外部エネルギＥout と断面積Ａとを基に下記（１）式の演算を行って歪みエネルギ開放率Ｇを算出する。
【００３３】
Ｇ＝Ｅout ／Ａ−Ｅst／Ａ …………（１）
次いで、ステップＳ２ｅに移行して、算出した歪みエネルギ開放率Ｇが入力された限界歪み開放率Ｇc に達したか否かを判定し、Ｇ＜Ｇc であるときには限界歪みエネルギ開放率Ｇc に達していないものと判断してステップＳ２ｆに移行して、現在の変位ｘに所定の増加分Δｘを加算した値（ｘ＋Δｘ）を新たな変位ｘとして設定してから前記ステップＳ２ｂに戻り、Ｇ＝Ｇc であるときには限界歪みエネルギ開放率Ｇc に達したものと判断してステップＳ２ｇに移行する。
【００３４】
ここで、歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達した時には、試験片モデルは変形によるエネルギ消費（歪みエネルギ）が不可能となり、破断する状態となったことを判断することができる。
【００３５】
ステップＳ２ｇでは、対象とする有限要素の全ての歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達したか否かを判定し、限界歪みエネルギ開放率Ｇc に達していない有限要素が存在するときには、前記ステップＳ２ａに戻り、限界歪みエネルギ開放率Ｇc に達していない有限要素が存在しないときには処理を終了して、処理結果であるＨＩＰＳの変位、応力、歪み、歪みエネルギ、外部エネルギ、反力及び歪みエネルギ開放率を例えば内蔵する記憶装置の所定記憶領域に格納してから図２のステップＳ３に移行する。
【００３６】
このステップＳ３では、前記ステップＳ２で格納されたＨＩＰＳの反力を変位方向に対して直角の断面積で除して単位形状としての公称応力を算出すると共に、強制変位を変位方向の単位形状長さで除して単位形状としての公称歪みを算出する。
【００３７】
次いで、ステップＳ４に移行して、ステップＳ３で算出したＨＩＰＳの公称応力及び公称歪みをもとに歪み速度の異なったＨＩＰＳの応力歪み特性を算出し、これを表す応力歪み特性線図を作成し、これを記憶装置の所定記憶領域に格納してからステップＳ５に移行する。なお、このステップＳ４の処理は市販されている非線型有限要素法解法プログラム「ＡＢＡＱＵＳ」（商品名）（米国ＨＫＳ社製）を適用することができる。
【００３８】
このステップＳ５では、ステップＳ５で格納された応力歪み特性線図をもとに衝撃試験方法に応じた試験形状データ及び試験方法に応じた境界条件を入力することにより、各有限要素について破断が生じるまでの外部エネルギ、内部エネルギ及び運動エネルギを算出し、これらを試験片全体で積算し、積算結果を記憶装置の所定記憶領域に記憶するエネルギ演算処理を実行する。
【００３９】
このエネルギ演算処理の具体例は、図４に示すように、先ず、ステップＳ５ａで１つの有限要素について初期変位を設定し、次いでステップＳ５ｂに移行して要素各々の歪みεを算出してからステップＳ５ｃに移行する。
【００４０】
このステップＳ５ｃでは、ステップＳ５ｂで算出した歪みεが前述したステップＳ４の応力歪み特性線図から求まる歪み速度に応じた破断歪みεb に達したか否かを判定し、ε＜εb であるときには、破断歪みεb に達していないものと判断して、ステップＳ５ｄに移行して、変位ｘに変位増加分Δｘを加算した値を新たな変位ｘとして設定してから前記ステップＳ５ｂに戻り、ε＝εb であるときには、破断歪みεb に達して、切欠き部から生じる亀裂（クレーズではない巨視的亀裂＝破断）が生じたものと判断してステップＳ５ｅに移行する。
【００４１】
このステップＳ５ｅでは、破断歪みεb に達した有限要素を削除し、次いでステップＳ５ｆに移行して、該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出して、試験片全体での積算を行い、その積算結果を記憶装置の所定記憶領域に更新記憶する。
【００４２】
次いで、ステップＳ５ｇに移行して、対象とする有限要素の全てが破断歪みεb に達したか否かを判定し、破断歪みεb に達していない有限要素即ち削除されていない有限要素について前記ステップＳ５ａ〜Ｓ５ｆの処理を繰り返し、全ての有限要素が破断歪みεb に達したときには破断状態に達したと判断して対象となる要素全てを削除して処理を終了し、図２のステップＳ６に移行する。
【００４３】
このステップＳ６では、記憶装置の所定領域に更新記憶されている試験片全体での外部エネルギ、内部エネルギ及び運動エネルギから試験片に蓄積された吸収エネルギを算出する。
【００４４】
次いで、ステップＳ７に移行して、吸収エネルギを切欠き部の最少断面積で除することにより、アイゾット衝撃値を算出し、これを記憶装置の所定記憶領域に記憶すると共に、表示装置に表示し、必要に応じてプリンタで印字してから処理を終了する。
【００４５】
ここで、図２の処理において、ステップＳ２〜Ｓ４の処理が応力歪み特性演算手段に対応すると共に、ステップＳ２の処理が特性決定手段に対応し、ステップＳ３の処理が公称応力歪み算出手段に対応し、ステップＳ４の処理が応力歪み特性算出手段に対応し、ステップＳ５の処理がエネルギ演算手段に対応し、ステップＳ６及びＳ７の処理が衝撃値演算手段に対応しており、またステップＳ６及びＳ７の処理は市販されている衝撃強度解析プログラム「ＲＡＤＩＯＳＳ」（商品名）（仏ＭＥＣＡＬＯＧ社製）を適用することができる。
【００４６】
このように、図２の処理によれば、高分子多相系材料を構成する連続相及び分散相成分毎のポアソン比、歪み速度を異ならせた応力歪み特性及び限界歪みエネルギ開放率Ｇc とを含む材料データと、ＨＩＰＳの形態を表す単位形状データと、境界条件を入力することにより、複合材料としてのＨＩＰＳの応力歪み特性を演算し、この演算結果から衝撃試験方法に応じた試験形状データと試験方法に応じた境界条件を入力することにより、試験片全体の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これらをエネルギに基づいて試験片に蓄積された吸収エネルギを算出し、これを試験片の切欠き部の最少面積で除することにより、仮想のＨＩＰＳであってもアイゾット衝撃値を算出することができる。
〔実施例〕
下記のように作成したゴム強化ポリスチレン樹脂組成物及びポリスチレンを適宜組み合わせて１１種類の試験片ＴＰ1 〜ＴＰ11を作成し、これをアイゾット衝撃試験機で実測した実測衝撃値と図２の演算処理によって求めた演算衝撃値との比較結果を下記表１に表す。
【００４７】
▲１▼ポリブタジエン（宇部興産製：１３ＨＢ）を溶解したスチレンモノマーを攪拌下で重合する際に、攪拌回転数を変化させ、０．６８ミクロン，１．２ミクロン及び２．７ミクロンの平均粒子径を持ち、ゴム質重合体を１１重量％含有するゴム強化ポリスチレン組成物(I),(II)及び(III) を作成した。
【００４８】
▲２▼ポリブタジエン（旭化成工業製：ＮＦ３５）を溶解したスチレンモノマーを攪拌下で重合し、１．７ミクロンの平均粒子径を持ち、ゴム質重合体を１１重量％含有するゴム強化ポリスチレン樹脂組成物(IV)を作成した。
【００４９】
▲３▼スチレン−ブタジエンブロック共重合体をスチレンモノマーに溶解し攪拌下で重合することにより、ゴム質重合体粒子内部に単一のポリスチレンオクルージョンを有する平均粒子径が０．７ミクロンのゴム質重合体粒子１１重量％含有するゴム強化ポリスチレン樹脂組成物(V) を作成した。
【００５０】
▲４▼スチレンモノマーを重合することにより、樹脂組成物(I) 〜(V) の連続相ポリスチレンと同等の分子量を有するポリスチレン(VI)を作成すると共に、ポリスチレン(VI)に対してミネラルオイル３％を混合し、ポリスチレン(VII) を作成した。
【００５１】
これらの樹脂組成物を単独或いは押し出し混練により適宜組み合わせて試験片ＴＰ1 〜ＴＰ11とし、ＡＳＴＭ Ｄ２５６規格に基づき射出成形法１／４インチ試験片ＩＺＯＤ衝撃強度を測定した。
【００５２】
また、図２の演算処理においては、事前にゴム強化ポリスチレン樹脂組成物の原料に用いたポリブタジエン、ブタジエンスラチレン共重合体、及び別途作成した連続相を形成するポリスチレン相当品を圧縮成形によって１ミリメートル厚さのシートとし、これにより試験片を切り出し、これらを０．０１２５ｍ／sec 、２．０ｍ／sec 及び６．０ｍ／sec の引張り速度で図９に示す応力歪み曲線を測定した。一方、ゴム強化ポリスチレン樹脂組成物の電子顕微鏡写真を画像解析し、繰り返し最少単位の形状を求めた。ゴム強化ポリスチレン樹脂組成物（Ｖ）に関しては、ゴム質重合体粒子内部のポリスチレン粒子を考慮し、同心円（球）内部のポリスチレンの性質は粒子外部を取り巻く連続相ポリスチレンと同一とみなした。これらの実測値を入力し、各成分のポアソン比は、ポリスチレン、ポリブタジエンの文献値を、ポリスチレンの臨界エネルギ開放率も文献値を用いた。
【００５３】
また、図７に示すノッチ（切欠き部）付きアイゾット試験片形状を、ノッチ周辺及び区ラックが進展するノッチの延長線部分は細かく、ノッチから遠ざかるに従って粗く２０００個の有限要素に分割した。
【００５４】
さらに、アイゾット試験片の規定の位置に初速毎秒３メートルでハンマが衝突する条件のもとで、試験片が破断し分離するまでにようする試験片に蓄えられるエネルギをノッチ部分の最小断面積で除してアイゾット衝撃値を求めた。
【００５５】
なおさらに、ノッチ部分から生じるクラックは、解析により求まる歪みが繰り返し単位の数値解析によって得られる速度依存の応力歪み曲線の破断歪みεb を越える時点において生じるとして取り扱った。
【００５６】
【表１】
この表１から明らかなように、試験片ＴＰ1 〜ＴＰ3 及びＴＰ5 〜ＴＰ7 の実測値と演算値の弾性体粒子径及び弾性体量依存性の傾向は、アイゾット衝撃値を最大にする最適粒子径の存在、弾性体量減少によりアイゾット衝撃値が低下する減少を正しく表現し、試験片ＴＰ2 及びＴＰ4 の比較からゴムの種類によって最大値が異なることを正しく表している。
【００５８】
また、試験片ＴＰ7 及びＴＰ8 の比較においても可塑材効果を表現していることが明らかである。さらに試験片ＴＰ1 〜ＴＰ2 とＴＰ9 〜ＴＰ11との比較から、粒子径の小さい単一オクルージョンを有する弾性体粒子からなるゴム強化ポリスチレンは、少量の比較的大きなサイズの一般的弾性体粒子を添加することで飛躍的にそのアイゾット衝撃強度を高める現象が表現されている。
【００５９】
一方、重量分率１２．３％でゴム粒径が０．６８μｍ、１．２５μｍ及び２．７０μｍであるときのＨＩＰＳのステップＳ４で算出される応力歪み特性値（ＦＥＭ）と実験値との関係を図１０〜図１２に示す。
【００６０】
これらの図からも明らかなように、本実施形態で求めたＨＩＰＳの応力歪み特性は、実験値とほぼ一致しており、歪みエネルギ開放率Ｇが限界歪み開放率に達したときに要素の破断状態と判断することが最適であることが立証された。
【００６１】
次に、本発明の第２の実施形態を図１３について説明する。
この第２の実施形態においては、前述した第１の実施形態における図２のステップＳ５のエネルギ演算処理が図１３に示すように変更されていることを除いては前述した図２と同様の処理を実行し、対応するステップには同一ステップ番号を付して、その詳細説明はこれを省略する。
【００６２】
すなわち、第２の実施形態では、エネルギ演算処理が、図１３に示すように、先ず、ステップＳ５ｈで１つの有限要素について初期変位ｘを設定し、次いでステップＳ５ｉに移行して、変位ｘでの歪みエネルギＥstJ 外部エネルギＥout を算出し、次いでステップＳ５ｊに移行して、該当要素の断面積Ａを算出してからステップＳ５ｋに移行する。
【００６３】
このステップＳ５ｋでは、前述した（１）式の演算を行って歪みエネルギ開放率Ｇを算出し、次いでステップＳ５ｌに移行して歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達したか否かを判定し、Ｇ＜Ｇc であるときには限界歪みエネルギ開放率Ｇc に達していないものと判断してステップＳ５ｍに移行して、現在の変位ｘに変位増加分Δｘを加算した値を新たな変位ｘとして設定してから前記ステップＳ５ｉに戻り、Ｇ＝Ｇc であるときには有限要素が破断状態に達したと判断してステップＳ５ｎに移行する。
【００６４】
このステップＳ５ｎでは、破断状態に達した要素を削除し、次いでステップＳ５ｏに移行して、削除した要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これらを試験片全体で積算し、その積算結果を記憶装置の所定記憶領域に更新記憶する。
【００６５】
次いで、ステップＳ５ｐに移行して、対象とする要素全てが限界歪みエネルギ開放率Ｇc に達したか否かを判定し、限界歪みエネルギ開放率Ｇc に達していない要素があるとき即ち削除されずに残っている要素があるときには、前記ステップＳ５ｈに戻って、残りの要素についてステップＳ５ｈ〜Ｓ５ｎの処理を行い、全ての要素が限界歪みエネルギ開放率Ｇc に達したときには、試験片が破断状態に達したものと判断して、ステップＳ５ｑに移行して、対象要素を全て削除してから図２のステップＳ６に移行する。
【００６６】
この図１１のエネルギ演算処理においても、各有限要素について変位を増加させながら各要素の歪みエネルギ開放率Ｇを演算し、この歪みエネルギ開放率Ｇが限界歪みエネルギ開放率Ｇc に達したときに該当有限要素が破断状態となったものと判断して、該当有限要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを試験片全体で積算し、その積算値を記憶装置の所定記憶領域に更新記憶するので、全ての要素の解析が終了した時点で、試験片に蓄積された吸収エネルギを算出し、これを試験片のノッチ部分の最小断面積で除することにより目的とするアイゾット衝撃値を算出することができる。
【００６７】
なお、上記実施例においては、高分子多相系材料がゴム強化ポリスチレン樹脂組成物である場合について説明したが、これに限定されるものではなく、連続相を形成する物質としては、前述したポリスチレンの他、アクリロニトル−スチレン共重合体、ポリフェニレンオキサイド、ポリスチレン−ポリフェニレンオキサイド混合物といったスチレン系樹脂及びこれらに可塑材・難燃剤を添加した組成物があり、他にもポリ塩化ビニル、ポリ塩化ビニリデン、ポリメチルメタクリレート、ポリブタジエンといった非晶性樹脂及びこれらの組成物にも好適に適用することができる。
【００６８】
また、分散相を形成する物質としては、ポリブタジエン、ブタジエン−スチレン共重合体、アクリロニトリル−ブタジエン共重合体、ＥＶＡ、ＥＰＤＭといった連続相より弾性率の低い物質、ボリオレフィン、ボリアミド、ボリエステルといった結晶性高分子、炭酸カルシウム、酸化チタン、カーボンブラックといった弾性率の高い物質があげられる。これらポリオレフィン、ポリエステル、ポリアミドといった結晶性高分子材料の場合には、結晶相を分散相として、非晶相を連続相として取り扱うことで図２の演算処理の適用が可能である。
【００６９】
さらに、分散相の分布形態は、球に限らず、紡錘体、柱状、面等任意に与えることができるし、分布形状は電子顕微鏡写真を画像処理して入力しても、架空の分散形状を定義して入力してもよく、分散相と連続相の比率も任意に設定することができる。
【００７０】
なおさらに、上記実施形態においては、アイゾット衝撃強度を計測する場合について説明したが、これに限定されるものではなく、シャルピー衝撃強度やその他の衝撃強度を計測することができる。この場合にはエネルギ演算手段で衝撃試験方法に応じた衝撃試験モデル及び衝撃試験変位を与えて解析すればよい。
【００７１】
また、上記実施形態においては、衝撃試験片として、切欠き部１１を有するモデルを適用した場合について説明したが、これに限定されるものではなく、切欠き部を有さない衝撃試験モデルについても、破断位置を特定することにより、吸収エネルギに基づいて衝撃値を容易に算出することができる。
【００７２】
【発明の効果】
以上説明したように、請求項１及び請求項７に係る発明によれば、入力手段によって、高分子多相系材料を構成する複合材料毎の応力歪み特性とポアソン比でなる材料データと、高分子多相系材料の形態を表す形状データと所定の境界条件を入力することにより、非線形有限要素法に従って試料モデルを分割した有限要素毎に歪みエネルギ開放率が限界歪みエネルギ開放率に達するまで変位を増加させることにより、高分子多相系材料の破断状態に至る応力歪み特性を演算し、この演算結果から、衝撃試験変位を与えたときの各要素が破断状態に達するまで前記衝撃試験変位を増加させながら、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するステップと、積算された各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギを前記切欠き部の最小断面積で除して衝撃値を演算するように構成したので、高分子多相系材料の各構成材料の歪み速度の異なる応力歪み特性、ポアソン比及び限界歪みエネルギ開放率を与えるだけで、衝撃強度を計測することが可能となり、実際に該当する高分子多相系材料が存在しない場合でも確実に衝撃強度を計測することができ、高分子多相系材料を設計開発する場合に大いに貢献することができるという効果が得られる。
【００７３】
また、請求項２に係る発明によれば、請求項１の発明において、前記応力歪み特性演算手段は、試験モデルに対する変位の増分を設定する変位設定手段と、該変位設定手段で設定された衝撃試験変位を与えたときの各要素の歪みエネルギと外部エネルギとを演算し、両エネルギと各要素の断面積とに基づいて歪みエネルギ開放率を演算する歪みエネルギ開放率演算手段と、該歪みエネルギ開放率演算手段の歪みエネルギ開放率が限界歪みエネルギ開放率に達したか否かを判定し、限界歪みエネルギ開放率に達していないときに前記変位を増加させ、限界歪みエネルギ開放率に達したときに該当要素の変位、応力、沈み、歪みエネルギ、外部エネルギ、反力及び歪みエネルギ開放率を決定する特性決定手段と、該特性決定手段で決定された反力及び衝撃試験変位に基づいて公称応力及び公称歪みを演算するする応力歪み演算手段と、該応力歪み演算手段の演算結果に基づいて歪み速度の異なった応力歪み特性を算出する応力歪み特性算出手段とを備えているので、試験片の歪みエネルギ開放率と限界歪みエネルギ開放率とを比較し、限界歪みエネルギ開放率に達したときに試験片が歪むことによってエネルギ消費をし尽くすことができず破断する状態となったことを確実に判断することができ、構造の定義された高分子多相系材料の応力歪み特性を正確に計測することができるという効果が得られる。
【００７４】
さらに、請求項３に係る発明によれば、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みを算出し、当該歪みが破断歪みに達したときに破断状態と判断して該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されているので、各要素の破断状態を確実に判断することができるという効果が得られる。
【００７５】
さらにまた、請求項４に係る発明によれば、請求項１又は２の発明において、前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みエネルギと外部エネルギを算出し、両エネルギと要素断面積とで算出される歪みエネルギ開放率が限界歪みエネルギ開放率に達したときに破綻状態と判断して該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されているので、請求項３の発明と同様に、各要素の破断状態を確実に判断することができるという効果が得られる。
【００７６】
なおさらに、請求項５に係る発明によれば、請求項１乃至４の発明において、前記高分子多相系材料は、連続相が非晶性樹脂で、且つ分散相がエラストマー粒子であるので、弾性率の大きく異なる材料の組み合わされた系の各成分の衝撃強度に対する寄与を確実に判断できるという効果が得られる。
【００７７】
また、請求項６に係る発明によれば、請求項５の発明において、前記材料データは複数の変形速度で求めた各連続相及び分散相の実測値を使用するので、材料データの入力が容易となるという効果が得られる。
【図面の簡単な説明】
【図１】本発明の第１の実施形態を示す概略構成図である。
【図２】第１の実施形態における衝撃強度演算処理手順の一例を示すフローチャートである。
【図３】図２の応力歪み特性演算処理の具体例を示すフローチャートである。
【図４】図２のエネルギ演算処理の具体例を示すフローチャートである。
【図５】ハイインパクトポリスチレンの形態を示す模式図である。
【図６】ハイインパクトポリスチレンの形態をモデル化する場合の説明図である。
【図７】ハイインパクトポリスチレンの３次元モデルを示す説明図である。
【図８】アイゾット衝撃試験片のモデルを示す説明図である。
【図９】ハイインパクトポリスチレンを構成する汎用ポリスチレン及びゴムの応力歪み特性を示す特性線図である。
【図１０】ゴム粒径が０．６８μｍであるときのＨＩＰＳの演算値と実測値との関係を示す特性線図である。
【図１１】ゴム粒径が１．２５μｍであるときのＨＩＰＳの演算値と実測値との関係を示す特性線図である。
【図１２】ゴム粒径が２．７０μｍであるときのＨＩＰＳの演算値と実測値との関係を示す特性線図である。
【図１３】本発明の第２の実施形態を示す図３に対応するエネルギ演算処理の具体例を示すフローチャートである。
【符号の説明】
１ コンピュータ
２ キーボード
３ フロッピィデスク装置
４ 表示装置
５ プリンタ
Claims (7)
- 連続相と分散相とが一定の３次元単位構造が繰り返される高分子多相系材料の衝撃強度を計測する高分子多相系材料の衝撃強度計測装置であって、前記連続相及び分散相の構成材料毎のポアソン比及び限界歪みエネルギ開放率、歪み速度依存データとしての応力歪み特性データ等の材料データと、３次元座標で与えられる単位形状データ、及び所定の境界条件を入力して、予め設定した試験モデルを有限要素に分割した各要素について、有限要素法に従って変位を与えたときの歪みエネルギ開放率を算出し、当該歪みエネルギ開放率が限界歪みエネルギ開放率に達するまで前記変位を増加させることにより、前記高分子多相系材料の歪み速度の異なった応力歪み特性を演算する応力歪み特性演算手段と、該応力歪み特性演算手段で演算した応力歪み特性に基づいて衝撃試験モデルを分割した各要素が破断状態に達するまで前記衝撃試験変位を増加させて、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するエネルギ演算手段と、該エネルギ演算手段で算出した各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギに基づいて衝撃値を演算する衝撃値演算手段とを備えたことを特徴とする高分子多相系材料の衝撃強度計測装置。
- 前記応力歪み特性演算手段は、試験モデルに対する変位の増分を設定する変位設定手段と、該変位設定手段で設定された変位を与えたときの各要素の歪みエネルギと外部エネルギとを演算し、両エネルギと各要素の断面積とに基づいて歪みエネルギ開放率を演算する歪みエネルギ開放率演算手段と、該歪みエネルギ開放率演算手段の歪みエネルギ開放率が限界歪みエネルギ開放率に達したか否かを判定し、限界歪みエネルギ開放率に達していないときに前記変位設定手段の変位を増加させ、限界歪みエネルギ開放率に達したときに該当要素の変位、応力、沈み、歪みエネルギ、外部エネルギ、反力及び歪みエネルギ開放率を決定する特性決定手段と、該特性決定手段で決定された反力及び衝撃試験変位に基づいて公称応力及び公称歪みを演算する公称応力歪み算出手段と、該応力歪み演算手段の演算結果に基づいて高分子多相系材料の歪み速度の異なった応力歪み特性を算出する応力歪み特性算出手段とを備えていることを特徴とする請求項１記載の高分子多相系材料の衝撃強度計測装置。
- 前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みを算出し、当該歪みが破断歪みに達したときに破断状態と判断して該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されていることを特徴とする請求項１又は２に記載の高分子多相系材料の衝撃強度計測装置。
- 前記エネルギ演算手段は、各要素について衝撃試験変位を増加させたときの歪みエネルギと外部エネルギを算出し、両エネルギと要素断面積とで算出される歪みエネルギ開放率が限界歪みエネルギ開放率に達したときに破綻状態と判断して該当要素の外部エネルギ、内部エネルギ及び運動エネルギを算出し、これを衝撃試験モデル全体で積算するように構成されていることを特徴とする請求項１又は２に記載の高分子多相系材料の衝撃強度計測装置。
- 前記高分子多相系材料は、連続相が非晶性樹脂で、且つ分散相がエラストマー粒子であることを特徴とする請求項１乃至４の何れかに記載の高分子多相系材料の衝撃強度計測装置。
- 前記材料データは複数の変形速度で求めた各連続相及び分散相の実測値を使用することを特徴とする請求項５記載の高分子多相系材料の衝撃強度計測装置。
- 一定の３次元単位構造が繰り返される高分子多相系材料の衝撃強度を計測する高分子多相系材料の衝撃強度計測方法であって、計測対象となる高分子多相系材料のポアソン比、歪み速度依存データとしての応力歪み特性データ等の材料データと、３次元座標で与えられる単位形状データ、及び所定の境界条件とに基づいて、予め設定した試験モデルを有限要素に分割した各要素について、有限要素法に従って変位を与えたときの歪みエネルギ開放率を算出しながら、当該歪みエネルギ開放率が限界歪みエネルギに達するまで前記変位を増加させることにより、高分子多相系材料の歪み速度の異なった応力歪み特性を演算するステップと、演算した応力歪み特性に基づいて衝撃試験変位を与えたときの各要素が破断状態に達するまで前記衝撃試験変位を増加させながら、各要素毎の外部エネルギ、内部エネルギ及び運動エネルギを演算し、これらを衝撃試験モデル全体で積算するステップと、積算された各エネルギに基づいて吸収エネルギを算出し、該吸収エネルギを前記切欠き部の最小断面積で除して衝撃値を演算することを特徴とする高分子多相系材料の衝撃強度計測方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12615496A JP3664806B2 (ja) | 1996-05-21 | 1996-05-21 | 高分子多相系材料の衝撃強度計測装置及び方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12615496A JP3664806B2 (ja) | 1996-05-21 | 1996-05-21 | 高分子多相系材料の衝撃強度計測装置及び方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09304257A JPH09304257A (ja) | 1997-11-28 |
| JP3664806B2 true JP3664806B2 (ja) | 2005-06-29 |
Family
ID=14928020
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP12615496A Expired - Fee Related JP3664806B2 (ja) | 1996-05-21 | 1996-05-21 | 高分子多相系材料の衝撃強度計測装置及び方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3664806B2 (ja) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4630024B2 (ja) * | 2004-08-31 | 2011-02-09 | 独立行政法人科学技術振興機構 | 実験装置と連動させる仮想実験インタフェース |
| JP4602776B2 (ja) * | 2005-01-18 | 2010-12-22 | 株式会社ブリヂストン | ゴム材料の変形挙動予測方法及びゴム材料の変形挙動予測装置 |
| JP4581758B2 (ja) * | 2005-03-11 | 2010-11-17 | 横浜ゴム株式会社 | ゴムと樹脂との複合材料の評価方法 |
| JP4602929B2 (ja) * | 2006-03-30 | 2010-12-22 | 株式会社ブリヂストン | ゴム材料の変形挙動予測装置及びゴム材料の変形挙動予測方法 |
| JP2008027026A (ja) * | 2006-07-19 | 2008-02-07 | Sumitomo Chemical Co Ltd | 破壊挙動解析方法 |
| JP5169225B2 (ja) * | 2007-05-25 | 2013-03-27 | 横浜ゴム株式会社 | 複合材料モデルの作成装置、シミュレーション装置、複合材料モデルの作成プログラム及び複合材料モデルの作成装置の動作方法 |
| JP5771935B2 (ja) * | 2010-10-04 | 2015-09-02 | 横浜ゴム株式会社 | シミュレーションモデルの作成方法及びそのコンピュータプログラム、シミュレーション方法及びそのコンピュータプログラム、並びにシミュレーションモデルの作成装置 |
| JP7157978B2 (ja) * | 2018-11-07 | 2022-10-21 | 日本メナード化粧品株式会社 | 原料粉末の力学特性から固形粉末化粧料の耐衝撃性を予測する方法 |
| JP7280117B2 (ja) * | 2019-05-31 | 2023-05-23 | 日鉄神鋼建材株式会社 | アイゾット衝撃試験方法 |
| CN117907150B (zh) * | 2024-03-20 | 2024-05-17 | 中国科学院力学研究所 | 一种活性材料冲击释能特性获取方法 |
-
1996
- 1996-05-21 JP JP12615496A patent/JP3664806B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH09304257A (ja) | 1997-11-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3664806B2 (ja) | 高分子多相系材料の衝撃強度計測装置及び方法 | |
| CN113420376B (zh) | 基于多尺度的碳纤维复合材料抗冲击力学性能仿真方法 | |
| Hassanifard et al. | On the strain-life fatigue parameters of additive manufactured plastic materials through fused filament fabrication process | |
| Jendli et al. | Multi-scales modelling of dynamic behaviour for discontinuous fibre SMC composites | |
| JP5548180B2 (ja) | 高分子材料のシミュレーション方法 | |
| Grellmann et al. | Technical material diagnostics-fracture mechanics of filled elastomer blends | |
| Hamdi et al. | Experimental and FEM analysis of scratch behavior on polypropylene thin films: Effect of film orientation and ethylene monomer content | |
| Dundar et al. | Investigation for impact behavior of acrylonitrile-butadiene-styrene amorphous thermoplastic | |
| Jinnai | Electron microscopy for polymer structures | |
| Rohde et al. | Nanostructured thermoset/thermoset blends compatibilized with an amphiphilic block copolymer | |
| Yao et al. | Atomic insights into shock-induced spalling of polyurea by molecular dynamics simulation | |
| Kim et al. | Polyamide-nylon 6 particulate polycarbonate composites with outstanding energy-absorbing properties | |
| Corté et al. | Analysis of polymer blend morphologies from transmission electron micrographs | |
| Muratoǧlu et al. | Microstructural fracture processes accompanying growing cracks in tough rubber-modified polyamides | |
| CN112417724A (zh) | 一种基于反向建模的碳酸钙填充复合材料力学性能预测方法 | |
| Kulmi et al. | A molecular dynamics study of the failure modes of a glassy polymer confined between rigid walls | |
| Cavaillé et al. | Dynamic mechanical spectrometry for the study of multiphase polymer materials | |
| Géhant et al. | Criteria for cavitation of rubber particles: Influence of plastic yielding in the matrix | |
| Yamashita et al. | A study of the microscopic plastic deformation process in poly (methylmethacrylate)/acrylic impact modifier compounds by means of small angle X-ray scattering | |
| Maurya et al. | Free vibration response of agglomerated carbon nanotube-reinforced nanocomposite plates | |
| Inukai et al. | Properties and Interfacial Structure Analysis of MWCNT/ESBS Composites | |
| Chiang | Predictions of elastic properties of reinforced polymers accounting for the agglomeration of fillers | |
| Fard et al. | On the strain-life fatigue parameters of additive manufactured plastic materials through used filament fabrication process | |
| Rizzieri | Rubber particle cavitation in toughened poly (methyl methacrylate) | |
| Pletz et al. | Combined Crack Initiation and Crack Growth Model for Multi-Layer Polymer Materials. Materials 2022, 15, 3273 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20040127 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20041026 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050315 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050330 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080408 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090408 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090408 Year of fee payment: 4 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090408 Year of fee payment: 4 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090408 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100408 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |