特許 7302112 架橋反応シミュレーション装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B1)

(11)【特許番号】

(24)【登録日】2023-06-23

(45)【発行日】2023-07-03

(54)【発明の名称】架橋反応シミュレーション装置

(51)【国際特許分類】

   G06F 30/20 20200101AFI20230626BHJP

【ＦＩ】

G06F30/20

【請求項の数】 36

(21)【出願番号】P 2023052538

(22)【出願日】2023-03-29

【審査請求日】2023-03-29

【早期審査対象出願】

(73)【特許権者】

【識別番号】000219602

【氏名又は名称】住友理工株式会社

(74)【代理人】

【識別番号】110000648

【氏名又は名称】弁理士法人あいち国際特許事務所

(72)【発明者】

【氏名】住田 弘志

(72)【発明者】

【氏名】酒井 潔

(72)【発明者】

【氏名】杉浦 隆典

(72)【発明者】

【氏名】木村 憲仁

【審査官】合田 幸裕

(56)【参考文献】

【文献】特開２０１３－１１６５８３（ＪＰ，Ａ）

【文献】特開２００８－００１８７８（ＪＰ，Ａ）

【文献】特開２０１７－０７１０５７（ＪＰ，Ａ）

【文献】中国特許出願公開第１１４８３３９７３（ＣＮ，Ａ）

【文献】中川泰忠，部品の設計及び製造工程を適正化するシミュレーション技術，東芝レビュー，日本，株式会社東芝，2008年08月01日，Vol. 63, No. 8，pages 20-24，ISSN: 0372-0462

【文献】雁部竜也 ほか，熱硬化性樹脂成形品の残留応力分布・接着界面強度解析，富士電機技報，日本，富士電機株式会社，2016年03月30日，Vol. 89, No. 1，pages 30-34，ISSN: 2187-1817

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｆ ３０／２０

ＩＥＥＥ Ｘｐｌｏｒｅ

ＪＳＴＰｌｕｓ（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

シミュレーションに用いるデータを記憶する記憶部と、
対象ワークモデルのポリマー部の架橋反応時の伝熱解析を行う伝熱解析部と、
前記伝熱解析の結果を用いて前記ポリマー部の架橋反応の反応率の解析を行う架橋反応解析部と、
を備える架橋反応シミュレーション装置であって、
前記記憶部は、
原料ポリマーを含んで構成される前記ポリマー部を有する前記対象ワークモデルと、
基準反応温度での基準反応時間のときの架橋反応の反応量に対する、対象反応温度での対象反応時間のときの架橋反応の反応量の比を等価反応量として定義し、アレニウスプロットの傾きを表す傾き係数を含んで定義される等価反応量算出モデルと、
前記ポリマー部の架橋反応の進行度に応じて設定された前記傾き係数と、を記憶し、
前記架橋反応解析部は、
前記伝熱解析の結果として、架橋反応における前記対象ワークモデルの前記ポリマー部の各要素について時刻毎の温度を取得する温度取得部と、
取得した前記架橋反応における前記ポリマー部の各要素についての時刻毎の温度、前記等価反応量算出モデル、対象時刻における前記架橋反応の前記進行度に応じた前記傾き係数に基づいて、時刻毎の前記ポリマー部の前記等価反応量を算出し、算出した前記等価反応量に基づいて前記ポリマー部の架橋反応の反応率を算出する反応率算出処理部と、を備える、架橋反応シミュレーション装置。

【請求項2】

前記ポリマー部は、防振性能を発揮するように構成されたゴム部であり、
前記対象ワークモデルは、防振ゴム装置のモデルである、請求項１に記載の架橋反応シミュレーション装置。

【請求項3】

前記傾き係数は、架橋反応の反応率がピークに到達するまでの反応進行期と、反応率がピークを越えた反応戻り期とで、異なる値に設定される、請求項１または２に記載の架橋反応シミュレーション装置。

【請求項4】

前記傾き係数は、前記反応進行期および前記反応戻り期を含む全期間を等間隔に区分して、各区分に応じた値に設定される、請求項３に記載の架橋反応シミュレーション装置。

【請求項5】

前記傾き係数は、前記反応率がピークに到達するまでの架橋反応において、架橋反応が進行し始める反応進行初期、反応進行速度が反応進行初期よりも早くなる反応促進期、反応進行速度が反応促進期よりも遅くなる反応進行後期のそれぞれで、異なる値に設定される、請求項３に記載の架橋反応シミュレーション装置。

【請求項6】

前記傾き係数は、前記反応進行期より前記反応戻り期の方が小さな値に設定される、請求項３に記載の架橋反応シミュレーション装置。

【請求項7】

前記傾き係数は、
前記反応進行期より前記反応戻り期の方が小さな値に設定され、
前記反応進行期は、架橋反応が進行し始める反応進行初期、反応進行速度が反応進行初期よりも早くなる反応促進期、反応進行速度が反応促進期よりも遅くなる反応進行後期に区分され、
前記反応促進期が前記反応進行初期および前記反応進行後期より大きな値に設定される、請求項３に記載の架橋反応シミュレーション装置。

【請求項8】

前回時刻の反応率に基づいて、対象時刻で使用する前記傾き係数を決定し、
対象時刻における温度、等価反応量算出モデル、前記傾き係数に基づいて、対象時刻の微小時間における等価反応量増加量を算出し、
前記等価反応量増加量に基づいて、架橋反応開始から対象時刻までの等価反応量積算値を算出し、
架橋反応開始から対象時刻までの前記等価反応量積算値に基づいて、対象時刻における反応率を算出する、請求項１または２に記載の架橋反応シミュレーション装置。

【請求項9】

前記ポリマー部の架橋反応の反応率と前記傾き係数との対応関係を定義した第一関係データマップを記憶し、
前回時刻の反応率および前記第一関係データマップを用いて、前記傾き係数を決定する、請求項８に記載の架橋反応シミュレーション装置。

【請求項10】

前記ポリマー部の架橋反応の反応率と前記傾き係数との対応関係を定義した第一関数を記憶し、
前回時刻の反応率および前記第一関数を用いて、前記傾き係数を決定する、請求項８に記載の架橋反応シミュレーション装置。

【請求項11】

前記第一関数は、前記ポリマー部の架橋反応の進行度に応じて設定された複数の反応率区分に応じて、異なる関数に設定されている、請求項１０に記載の架橋反応シミュレーション装置。

【請求項12】

前記等価反応量積算値と対象時刻における反応率との対応関係を定義した第二関係データマップを記憶し、
対象時刻における前記等価反応量積算値および前記第二関係データマップを用いて、対象時刻における反応率を決定する、請求項８に記載の架橋反応シミュレーション装置。

【請求項13】

前記等価反応量積算値と対象時刻における反応量との対応関係を定義した第二関数を記憶し、
対象時刻における前記等価反応量積算値および前記第二関数を用いて、対象時刻における反応率を決定する、請求項８に記載の架橋反応シミュレーション装置。

【請求項14】

前記第二関係データマップは、前記対象ワークモデルの前記ポリマー部に相当する試験対象ポリマー材料を用いて、架橋反応特性試験機により測定された基準反応温度において架橋反応時間と前記試験対象ポリマー材料に発生したトルクとの関係に基づく基準反応曲線より設定されている、請求項１２に記載の架橋反応シミュレーション装置。

【請求項15】

前記第二関数は、前記対象ワークモデルの前記ポリマー部に相当する試験対象ポリマー材料を用いて、架橋反応特性試験機により測定された基準反応温度において架橋反応時間と前記試験対象ポリマー材料に発生したトルクとの関係に基づく基準反応曲線より設定されている、請求項１３に記載の架橋反応シミュレーション装置。

【請求項16】

前記第二関数は、前記ポリマー部の架橋反応の進行度に応じて設定された複数の反応率区分に応じて、異なる関数に設定されている、請求項１５に記載の架橋反応シミュレーション装置。

【請求項17】

架橋反応の反応率は、
架橋反応の反応率がピークに到達するまでの反応進行期において、架橋反応の開始時を０％とし、架橋反応の反応率の増加度を０％～１００％までの範囲で定義し、
反応率がピークを越えた反応戻り期において、反応率ピークからの反応率の低下度を１００％に加算した値として定義する、請求項１または２に記載の架橋反応シミュレーション装置。

【請求項18】

さらに、前記架橋反応解析部の解析結果に基づいて、前記対象反応時間に応じた前記反応率から得られる値を表示する表示部を備える、請求項１または２に記載の架橋反応シミュレーション装置。

【請求項19】

さらに、前記架橋反応解析部により解析された前記ポリマー部の架橋反応の反応率を用いて構造解析を行う構造解析部を備え、
前記構造解析部は、
前記架橋反応解析部の前記反応率算出処理部により算出された前記ポリマー部の各要素における前記反応率を取得する反応率取得部と、
前記ポリマー部において、取得した前記反応率に応じた弾性率を割り当てる弾性率割当部と、
前記ポリマー部に前記弾性率を割り当てた状態で構造解析を行うことにより、前記対象ワークモデルの特性を取得する特性取得部と、を備える、請求項１または２に記載の架橋反応シミュレーション装置。

【請求項20】

さらに、前記構造解析の結果に基づいて、前記ポリマー部の架橋反応に用いられる成形型モデルの温度と、架橋反応開始から成形型モデルの脱型までの型内反応時間と、に応じた前記特性を表示する表示部を備える、請求項１９に記載の架橋反応シミュレーション装置。

【請求項21】

前記対象ワークモデルの前記特性は、前記対象ワークモデルの弾性率である、請求項１９に記載の架橋反応シミュレーション装置。

【請求項22】

前記記憶部は、さらに、
成形型モデルと、
さらにカーボンブラックを含んで構成される前記ポリマー部を有する前記対象ワークモデルと、
前記原料ポリマーに対する前記カーボンブラックの質量比と前記ポリマー部の熱拡散率との関係を表す熱拡散率特性と、を記憶し、
前記伝熱解析部は、
前記対象ワークモデルにおける前記原料ポリマーに対する前記カーボンブラックの質量比、および、前記成形型モデルの温度条件を入力する条件入力部と、
前記条件入力部により入力された前記質量比と前記記憶部に記憶された前記熱拡散率特性とに基づいて、前記対象ワークモデルの前記ポリマー部の熱拡散率であるポリマー熱拡散率を決定するポリマー熱拡散率決定部と、
前記成形型モデルに前記対象ワークモデルを配置した状態で、前記ポリマー熱拡散率決定部により決定された前記ポリマー熱拡散率および前記記憶部に記憶された前記温度条件を用いて、伝熱解析を行う解析部と、を備える、請求項１または２に記載の架橋反応シミュレーション装置。

【請求項23】

前記熱拡散率特性は、前記原料ポリマーの種類に依存しない関係である、請求項２２に記載の架橋反応シミュレーション装置。

【請求項24】

前記熱拡散率特性は、複数種の原料ポリマーで、異なるカーボンブラックの質量比で実験した場合の実測値より設定された関係である、請求項２２に記載の架橋反応シミュレーション装置。

【請求項25】

前記熱拡散率特性は、前記カーボンブラックの質量比に対して線形の関係を有する、請求項２２に記載の架橋反応シミュレーション装置。

【請求項26】

前記対象ワークモデルは、
前記ポリマー部と、
前記ポリマー部に接合された接合部材と、を備え、
前記記憶部は、さらに、
前記成形型モデルに前記対象ワークモデルを配置した状態で、前記成形型モデルと前記接合部材との間の熱伝達率である接触熱伝達率を記憶し、
前記解析部は、前記成形型モデルに前記対象ワークモデルを配置した状態で、前記ポリマー熱拡散率決定部により決定された前記ポリマー熱拡散率、前記記憶部に記憶された前記温度条件および前記接触熱伝達率を用いて、前記伝熱解析を行う、請求項２２に記載の架橋反応シミュレーション装置。

【請求項27】

前記条件入力部は、
前記成形型モデルから脱型された前記対象ワークモデルの周囲の気温を含む外気条件が入力されるようになっており、
前記記憶部は、さらに
空気の熱伝達係数である空気熱伝達係数を記憶し、
前記解析部は、前記対象ワークモデルの前記ポリマー部を前記成形型モデル内にて架橋反応させた後に前記成形型モデルから前記対象ワークモデルを脱型した状態で、前記ポリマー熱拡散率決定部により決定された前記ポリマー熱拡散率、前記記憶部に記憶された前記空気熱伝達係数および前記外気条件を用いて、前記伝熱解析を行う、請求項２２に記載の架橋反応シミュレーション装置。

【請求項28】

前記記憶部は、さらに、成形型モデルを記憶し、
さらに、前記対象ワークモデルの前記ポリマー部を前記成形型モデル内にて架橋反応させた後に前記成形型モデルを脱型する架橋反応工程に適用され、前記対象ワークモデルの前記ポリマー部の内部において前記成形型モデルの脱型に伴う気泡の発生を推定する内部気泡推定部を備え、
前記内部気泡推定部は、
前記対象ワークモデルの前記ポリマー部の架橋反応の反応率を取得する反応率取得部と、
脱型時における前記反応率に基づいて、前記対象ワークモデルの前記ポリマー部の内部における気泡の発生を推定する推定部と、
を備える、請求項１または２に記載の架橋反応シミュレーション装置。

【請求項29】

前記記憶部は、さらに
架橋反応開始からの経過時間と、前記対象ワークモデルの前記ポリマー部における架橋反応の進行度に対応する値であって前記ポリマー部に相当する試験対象ポリマー材料を用いて架橋反応特性試験機により測定可能なトルクと、の関係を定義する架橋反応曲線を記憶しており、
前記内部気泡推定部は、さらに、
前記反応率取得部により取得された前記反応率と前記記憶部に記憶された前記架橋反応曲線とに基づいて、取得された前記反応率に対応する前記トルクを算出するトルク算出部を備え、
前記推定部は、
脱型時における前記トルクに基づいて、前記対象ワークモデルの前記ポリマー部の内部における気泡の発生を推定する、請求項２８に記載の架橋反応シミュレーション装置。

【請求項30】

前記推定部は、脱型時における前記トルクが予め設定された閾値を超えていない場合に前記対象ワークモデルの前記ポリマー部の内部において気泡が発生したと推定する、請求項２９に記載の架橋反応シミュレーション装置。

【請求項31】

前記閾値は、複数種の原料ポリマーに対して統一した値に設定されている、請求項３０に記載の架橋反応シミュレーション装置。

【請求項32】

前記閾値は、ポリマー材料のブローポイントにおける架橋反応の反応率よりも大きな反応率に対応する前記トルクである、請求項３０に記載の架橋反応シミュレーション装置。

【請求項33】

前記反応率取得部は、前記ポリマー部の部位ごとの前記反応率を取得し、
前記トルク算出部は、前記ポリマー部の部位ごとの前記トルクを算出し、
前記推定部は、前記ポリマー部の部位ごとに気泡の発生を推定する、請求項２９に記載の架橋反応シミュレーション装置。

【請求項34】

さらに、前記内部気泡推定部による推定結果に基づいて、前記成形型モデルの温度と架橋反応開始から前記成形型モデルの脱型までの型内反応時間と応じて、前記気泡の発生の有無を表示する表示部を備える、請求項２９に記載の架橋反応シミュレーション装置。

【請求項35】

前記表示部は、前記架橋反応解析部による解析結果に基づいて、前記成形型モデルの温度と、架橋反応開始から前記成形型モデルの脱型までの型内反応時間と、に応じた前記反応率から得られる値を、前記気泡の発生の有無に合わせて表示する、請求項３４に記載の架橋反応シミュレーション装置。

【請求項36】

さらに、前記架橋反応解析部により解析された前記反応率を用いて構造解析を行うことにより、前記対象ワークモデルの特性を取得する構造解析部を備え、
前記表示部は、前記構造解析の結果に基づいて、前記成形型モデルの温度と、架橋反応開始から前記成形型モデルの脱型までの型内反応時間と、に応じた前記特性を、前記気泡の発生の有無に合わせて表示する、請求項３４に記載の架橋反応シミュレーション装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、架橋反応シミュレーション装置に関する。

【背景技術】

【0002】

従来より、ポリマーを架橋反応させることにより化学物質の物性を向上させることが行われている。このようなポリマーとしてゴムが例示される。ゴムは、原料ゴムに硫黄やその他の架橋剤，加硫促進剤等を加え、加熱を行なうことによって、ゴム分子鎖間あるいはその分子鎖の中に三次元網目状の架橋構造が形成されている。

【0003】

ゴムの加硫度、すなわち、ゴムの架橋反応の反応率を推定する技術として、非特許文献１に記載のものがある。上記の技術は、アレニウスの式に基づいて、ゴムの温度履歴から加硫度を推定し、適切な加硫条件（金型温度、加硫時間）を推定する。

【先行技術文献】

【非特許文献】

【0004】

【文献】有松利雄、“加硫工程設計の実際”、日本ゴム協会誌、第５９巻，第３号，（１９８６）

【発明の概要】

【発明が解決しようとする課題】

【0005】

しかしながら、上記の技術は精度が十分に高くないため、ゴムの架橋反応の反応率を精度よく推定することができなかった。

【0006】

本発明は、かかる背景に鑑みてなされたものであり、架橋反応の反応率を精度よく推定可能な架橋反応シミュレーション装置を提供しようとするものである。

【課題を解決するための手段】

【0007】

本発明の一態様は、
シミュレーションに用いるデータを記憶する記憶部と、
対象ワークモデルのポリマー部の架橋反応時の伝熱解析を行う伝熱解析部と、
前記伝熱解析の結果を用いて前記ポリマー部の架橋反応の反応量の解析を行う架橋反応解析部と、
を備える架橋反応シミュレーション装置であって、
前記記憶部は、
原料ポリマーを含んで構成される前記ポリマー部を有する前記対象ワークモデルと、
基準反応温度での基準反応時間のときの架橋反応の反応量に対する、対象反応温度での対象反応時間のときの架橋反応の反応量の比を等価反応量として定義し、アレニウスプロットの傾きを表す傾き係数を含んで定義される等価反応量算出モデルと、
前記ポリマー部の架橋反応の進行度に応じて設定された前記傾き係数と、を記憶し、
前記架橋反応解析部は、
前記伝熱解析の結果として、架橋反応における前記対象ワークモデルの前記ポリマー部の各要素について時刻毎の温度を取得する温度取得部と、
取得した前記架橋反応における前記ポリマー部の各要素についての時刻毎の温度、前記等価反応量算出モデル、対象時刻における前記架橋反応の前記進行度に応じた前記傾き係数に基づいて、時刻毎の前記ポリマー部の前記等価反応量を算出し、算出した前記等価反応量に基づいて前記ポリマー部の架橋反応の反応率を算出する反応率算出処理部と、を備える、架橋反応シミュレーション装置にある。

【発明の効果】

【0008】

本発明の一態様によれば、対象時刻における架橋反応の進行度に応じた傾き係数に基づいて、時刻毎のポリマー部の等価反応量を算出し、算出した等価反応量に基づいてポリマー部の架橋反応の反応率を算出する。これにより、架橋反応の進行度を考慮せず算出した等価反応量に基づいてポリマー部の架橋反応の反応率を算出した場合に比べて、架橋反応の反応率を精度よく推定することができる。

【図面の簡単な説明】

【0009】

【図1】実施形態１の架橋反応シミュレーション装置を示すブロック図。

【図2】実施形態１において、成形型モデルに対象ワークモデルを配置した状態を示す図。

【図3】実施形態１の対象ワークモデルを示す図。

【図4】実施形態１の架橋反応シミュレーション装置の動作を示すフローチャート。

【図5】実施形態１において、成形型モデルに対象ワークモデルを配置した状態の各部の温度の経時変化を示す図。

【図6】従来技術において、ポリマー部の温度の実測値と予測値とを示すグラフ。

【図7】実施形態１において、ポリマー部の温度の実測値と予測値とを示すグラフ。

【図8】実施形態１において、カーボンブラックの含有量に対する熱伝導率の変化量を示すグラフ。

【図9】実施形態１において、成形型モデルに対象ワークモデルが配置された状態の解析メッシュの様子を説明するための図。

【図10】実施形態１の加硫試験機の構造を示す模式図。

【図11】実施形態１の加硫試験機によって測定されたトルクの、反応率への変換方法を説明するための図であって、（ａ）は、トルクの、反応時間に対する変化を示すグラフであり、（ｂ）は、反応率の、反応時間に対する変化を示すグラフである。

【図12】実施形態１の架橋反応解析処理における、傾き係数の算出方法を説明するための図。

【図13】実施形態１の架橋反応解析処理における、等価反応量の算出方法を説明するための図。

【図14】実施形態１の架橋反応解析処理における、等価反応量増加量から等価反応量を算出する方法を説明するための図。

【図15】実施形態１の架橋反応解析処理における、等価反応量から反応率への変換方法を説明するための図。

【図16】従来技術において、ポリマー部の反応率の、実測値と予測値とを示す図。

【図17】実施形態１において、架橋反応の反応進行度に応じて区分し、活性化エネルギーを、各区分ごとに算出する方法を説明するための図

【図18】実施形態１において、ポリマー部の反応率の、実測値と予測値とを示す図。

【図19】実施形態１において、反応率の、反応時間に対する変化量の表現方法を説明するための図。

【図20】実施形態１において、反応率の、反応時間に対する変化量を示す図。

【図21】実施形態１の第二関数を示す図。

【図22】実施形態１における、対象ワークモデルの解析メッシュの様子を説明するための図。

【図23】ポリマー部に内部気泡が発生した状態を示す図。

【図24】架橋反応における、ブローポイントにおける反応率、脱型時における反応率および最終反応率を示す図。

【図25】（ａ）複数のサンプルのブローポイントにおける反応率を示すグラフ。（ｂ）複数のサンプルのブローポイントにおける、加硫試験機によって測定されたトルクを示すグラフ。

【図26】内部気泡を抑制するメカニズムを説明するための図。

【図27】実施形態１において、ポリマー部の反応率分布を示すヒートマップ図。

【図28】実施形態１において、金型温度が低温の場合の、反応時間に対する、対象ワークモデルの軸直角方向の弾性率の変化量を示す図。

【図29】実施形態１において、金型温度が中温の場合の、反応時間に対する、対象ワークモデルの軸直角方向の弾性率の変化量を示す図。

【図30】実施形態１において、金型温度が高温の場合の、反応時間に対する、対象ワークモデルの軸直角方向の弾性率の変化量を示す図。

【図31】実施形態１において、ポリマー部の弾性率と、内部気泡の発生状況の、反応時間に対する変化を示す図。

【図32】実施形態１において、内部気泡の発生予想結果を視覚的に表示する図。

【発明を実施するための形態】

【0010】

（実施形態１）
１．架橋反応シミュレーション装置１の構成
１－１．架橋反応シミュレーション装置１の全体構成
実施形態１の架橋反応シミュレーション装置１の全体構成について、図１を参照して説明する。本形態の架橋反応シミュレーション装置１は、原料ポリマーの分子鎖同士を架橋させる架橋反応についてのシミュレーションを行う。

【0011】

原料ポリマーは、分子鎖同士が架橋反応可能であれば特に限定されず、フェノール樹脂、エポキシ樹脂、メラミン樹脂、尿素樹脂、不飽和ポリエステル樹脂、アルキド樹脂、シリコーン樹脂、ポリウレタン等の熱硬化性樹脂、架橋ポリエチレン、架橋ポリプロピレン等の熱可塑性樹脂、天然ゴム、スチレン・ブタジエンゴム、ブタジエンゴム、イソプレンゴム等の合成ゴム、などの任意のポリマーを適宜に選択できる。本形態においては、天然ゴムおよび合成ゴムを含むゴムがポリマーとして用いられる。なお、ゴムの場合、硫黄等の架橋剤が添加された状態で加熱されることにより、ゴムを構成する分子鎖が架橋し、いわゆる加硫反応が起こる。

【0012】

図１に示すように、架橋反応シミュレーション装置１は、記憶部２と、伝熱解析部３と、架橋反応解析部４と、構造解析部５と、内部気泡推定部６と、表示部７と、を備える。

【0013】

１－２．記憶部２の構成
記憶部２は、シミュレーションに用いるデータを記憶する。データは、成形型モデルＭＭと、対象ワークモデルＷＭと、伝熱解析部３が用いる伝熱解析用データＴＤと、架橋反応解析部４が用いる架橋反応解析用データＲＤと、内部気泡推定部６が用いる内部気泡推定用データＢＤと、を含む。

【0014】

成形型モデルＭＭについて図２を参照して説明する。図２に示すように、成形型モデルＭＭは、金型１０と、金型１０に取付けられる熱板１１と、を含んで構成される。金型１０は、下側に位置する下型１０Ａと、下型１０Ａに上方から組付けられる上型１０Ｂと、を備える。熱板１１は、下型１０Ａの下面に取付けられる下熱板１１Ａと、上型１０Ｂの上面に取付けられる上熱板１１Ｂと、を含んで構成される。下型１０Ａには、上方に開口する下キャビティ１００Ａが形成されている。また、上型１０Ｂには、下方に開口する上キャビティ１００Ｂが形成されている。下型１０Ａと上型１０Ｂとが組付けられた状態で、下キャビティ１００Ａと上キャビティ１００Ｂとによって形成される空間内には、対象ワークモデルＷＭが配置される。

【0015】

対象ワークモデルＷＭについて図３を参照して説明する。図３に示すように、対象ワークモデルＷＭは、ポリマー部１２を備える。ポリマー部１２は、原料ポリマーを含む。さらに、ポリマー部１２は、酸化防止剤等の添加剤や、カーボンブラック等を含んでも良い。本形態のポリマー部１２は、原料ポリマーおよびカーボンブラックを含んで構成される。原料ポリマーは特に限定されず、ゴム、熱硬化性樹脂等、任意の材料を適宜に選択できる。本形態では、原料ポリマーはゴムにより構成される。本形態のポリマー部１２は、防振性能を発揮するように構成されたゴム部であり、対象ワークモデルＷＭは、防振ゴム装置のモデルである。ただし、ポリマー部１２は、原料ポリマーとして熱硬化性樹脂を含む構成としてもよいし、また、カーボンブラックを含まない構成としても良い。

【0016】

本形態のポリマー部１２は、上下方向に延びる軸線Ａに沿った筒状に形成されている。ポリマー部１２の外周部には、上下方向に延びる円筒形状に形成された外側接合部材１３（接合部材の一例）が配置されている。外側接合部材１３は、金属製、樹脂製、または金属と樹脂との複合体として構成される。本形態の外側接合部材１３は金属製であって、ポリマー部１２の外周部と接合されている。ポリマー部１２の内周部には、上下方向に延びる円筒形状に形成された内側接合部材１４（接合部材の一例）が配置されている。内側接合部材１４は、金属製、樹脂製、または金属と樹脂との複合体として構成される。本形態の内側接合部材１４は金属製であって、ポリマー部１２の内周部と接合されている。内側接合部材１４の上下方向の長さ寸法は、外側接合部材１３の上下方向の長さ寸法よりも大きく形成されている。ポリマー部１２の上面、および下面は凹形状に形成されている。ただし、ポリマー部１２の形状は上記の形状に限定されない。また、対象ワークモデルＷＭは、外側接合部材１３および内側接合部材１４の、双方または一方を備えない構成としても良い。

【0017】

図２に戻って、成形型モデルＭＭに対象ワークモデルＷＭを配置した状態で、成形型モデルＭＭと外側接合部材１３および内側接合部材１４との間には、それぞれ、隙間１５が形成されている。

【0018】

図１に戻って、伝熱解析用データＴＤは、原料ポリマーに対するカーボンブラックの質量比とポリマー部１２の熱拡散率との関係を示す熱拡散率特性ＴＳと、成形型モデルＭＭに対象ワークモデルＷＭを配置した状態で、成形型モデルＭＭと外側接合部材１３および内側接合部材１４との間の熱伝達率である接触熱伝達率ＣＨと、成形型モデルＭＭから対象ワークモデルＷＭを脱型した状態で、対象ワークモデルＷＭの周囲の空気の熱伝達係数である空気熱伝達係数ＡＣと、後述する条件入力部３０から入力される成形型モデルＭＭの温度条件ＴＣと、脱型した対象ワークモデルＷＭの周囲の気温を含む外気条件ＯＣと、を含む。

【0019】

架橋反応解析用データＲＤは、基準反応温度での基準反応時間のときの架橋反応の反応量に対する、対象反応温度での対象反応時間のときの架橋反応の反応量の比を等価反応量として定義し、アレニウスプロットの傾きを表す傾き係数を含んで定義される等価反応量算出モデルＥＭと、ポリマー部１２の架橋反応の進行度に応じて設定された傾き係数ＳＣと、を含む。架橋反応解析用データＲＤは、さらに、後述する第一関係データマップＤＭ１、第二関係データマップＤＭ２、第一関数Ｆ１、第二関数Ｆ２および基準反応曲線ＲＣを含む。

【0020】

内部気泡推定用データＢＤは、架橋反応開始からの経過時間と、対象ワークモデルＷＭのポリマー部１２における架橋反応の進行度に対応する値であってポリマー部１２に相当する試験対象ポリマー材料を用いて架橋反応特性試験機により測定可能なトルクと、の関係を定義する架橋反応曲線ＣＣを含む。

【0021】

１－３．伝熱解析部３の構成
伝熱解析部３は、対象ワークモデルＷＭのポリマー部１２の架橋反応時における伝熱解析を行う。図１に示すように、伝熱解析部３は、条件入力部３０と、ポリマー熱拡散率決定部３１と、解析部３２と、を備える。

【0022】

条件入力部３０は、対象ワークモデルＷＭにおける原料ポリマーに対するカーボンブラックの質量比、および、成形型モデルＭＭの温度条件ＴＣを入力する。条件入力部３０としては、キーボード、マウス、トラックボール、ジョイスティック等の入力装置でもよいし、半導体メモリ、ハードディスクメモリ等の外部記憶媒体でもよい。

【0023】

ポリマー熱拡散率決定部３１は、条件入力部３０により入力されたカーボンブラックの質量比と、記憶部２に記憶された熱拡散率特性ＴＳと、に基づいて、対象ワークモデルＷＭのポリマー部１２の熱拡散率であるポリマー熱拡散率を決定する。

【0024】

解析部３２は、成形型モデルＭＭに対象ワークモデルＷＭを配置した状態で、ポリマー熱拡散率決定部３１により決定されたポリマー熱拡散率および記憶部２に記憶された温度条件ＴＣを用いて、対象ワークモデルＷＭのポリマー部１２の架橋反応時における伝熱解析を行う。

【0025】

１－４．架橋反応解析部４の構成
架橋反応解析部４は、伝熱解析部３による伝熱解析の結果を用いてポリマー部１２の架橋反応の反応率の解析を行う。図１に示すように、架橋反応解析部４は、温度取得部４０と、反応率算出処理部４１と、を備える。

【0026】

温度取得部４０は、伝熱解析部３による伝熱解析の結果として、架橋反応における対象ワークモデルＷＭのポリマー部１２の各要素について時刻毎の温度を取得する。

【0027】

反応率算出処理部４１は、取得した架橋反応におけるポリマー部１２の各要素についての時刻毎の温度、等価反応量算出モデルＥＭ、対象時刻における架橋反応の進行度に応じた傾き係数ＳＣに基づいて、時刻毎のポリマー部１２の等価反応量を算出し、算出した等価反応量に基づいてポリマー部１２の架橋反応の反応率を算出する。

【0028】

１－５．構造解析部５の構成
構造解析部５は、架橋反応解析部４により解析されたポリマー部１２の架橋反応の反応率を用いて構造解析を行う。図１に示すように、構造解析部５は、温度取得部５０と、反応率取得部５１と、弾性率割当部５２と、特性取得部５３と、を備える。ただし、温度取得部５０は省略しても良い。

【0029】

温度取得部５０は、伝熱解析部３による伝熱解析の結果として、架橋反応における対象ワークモデルＷＭのポリマー部１２の各要素について時刻毎の温度を取得する。

【0030】

反応率取得部５１は、架橋反応解析部４の反応率算出処理部４１により算出されたポリマー部１２の各要素における反応率を取得する。

【0031】

弾性率割当部５２は、ポリマー部１２において、取得した反応率に応じた弾性率を割り当てる。

【0032】

特性取得部５３は、ポリマー部１２に弾性率を割り当てた状態で構造解析を行うことにより、対象ワークモデルＷＭの特性を取得する。

【0033】

１－６．内部気泡推定部６の構成
内部気泡推定部６は、対象ワークモデルＷＭのポリマー部１２を成形型モデルＭＭ内にて架橋反応させた後に成形型モデルＭＭを脱型する架橋反応工程に適用され、対象ワークモデルＷＭのポリマー部１２の内部において成形型モデルＭＭの脱型に伴う内部気泡の発生を推定する。図１に示すように、内部気泡推定部６は、反応率取得部６０と、トルク算出部６１と、推定部６２と、を備える。

【0034】

反応率取得部６０は、対象ワークモデルＷＭのポリマー部１２の架橋反応の反応率を取得する。上記の反応率は、金型１０を型開きした時点、すなわち型締め圧力が解放される時点での反応率である。なお、ポリマー部１２がゴムである場合、上記の反応率は脱型時加硫度とも呼ばれる。

【0035】

トルク算出部６１は、反応率取得部６０により取得された反応率（脱型時反応率）と記憶部２に記憶された架橋反応曲線ＣＣとに基づいて、取得された反応率に対応するトルク（脱型時トルクともいう）を算出する。

【0036】

推定部６２は、トルク算出部６１が算出したトルク（脱型時トルク）と、ポリマー部１２におけるブローポイント加硫度（後に詳述する）から算出したブローポイントトルク（後に詳述する）と、を比較して、ポリマー部１２の内部に気泡が発生するか否かを推定する。

【0037】

１－７．表示部７の構成
表示部７は、架橋反応解析部４の解析結果に基づいて、反応時間に応じた反応率から得られる値を表示する。また、表示部７は、構造解析部５の構造解析の結果に基づいて、ポリマー部１２の架橋反応に用いられる成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じた特性を表示する。

【0038】

また、表示部７は、内部気泡推定部６の推定結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じて、内部気泡の発生の有無を表示する。また、表示部７は、架橋反応解析部４の解析結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じた反応率から得られる値を、内部気泡の発生の有無に合わせて表示する。また、表示部７は、構造解析部５の構造解析の結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じた特性を、内部気泡の発生の有無に合わせて表示する。

【0039】

２．架橋反応シミュレーション装置１の全体の動作
図４を参照しつつ、本形態に係る架橋反応シミュレーション装置１の全体の動作を説明する。ただし、以下の説明は架橋反応シミュレーション装置１の動作の一例であり、架橋反応シミュレーション装置１の動作は以下の記載に限定されない。

【0040】

架橋反応シミュレーション装置１が起動されると、伝熱解析処理Ｓ１が実行される。伝熱解析処理Ｓ１により、架橋反応における対象ワークモデルＷＭのポリマー部１２の各要素について時刻毎の温度が得られる。詳細については後述する。

【0041】

次に、架橋反応解析処理Ｓ２が実行される。架橋反応解析処理Ｓ２においては、伝熱解析の結果として得られた、架橋反応における対象ワークモデルＷＭのポリマー部１２の各要素についての時刻毎の温度、等価反応量算出モデルＥＭ、対象時刻における架橋反応の進行度に応じた傾き係数ＳＣに基づいて、時刻毎のポリマー部１２の等価反応量を算出し、算出した等価反応量に基づいて、ポリマー部１２の架橋反応の反応率が算出される。詳細については後述する。

【0042】

次に、構造解析処理Ｓ３が実行される。構造解析処理Ｓ３においては、架橋反応解析部４により算出されたポリマー部１２の各要素における反応率に応じて弾性率を割り当て、ポリマー部１２に弾性率を割り当てた状態で構造解析を行うことにより、対象ワークモデルＷＭの特性を取得する。詳細については後述する。

【0043】

次に、内部気泡推定処理Ｓ４が実行される。内部気泡推定処理Ｓ４においては、架橋反応解析処理により算出されたポリマー部１２の反応率と記憶部２に記憶された架橋反応曲線ＣＣとに基づいて、ポリマー部１２の反応率に対応するトルクを算出し、脱型時におけるトルクに基づいて、対象ワークモデルＷＭのポリマー部１２の内部における内部気泡の発生を推定する。詳細については後述する。

【0044】

次に、表示処理Ｓ５が実行される。表示処理Ｓ５においては、内部気泡推定部６による推定結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じて、内部気泡の発生の有無を表示する。また、架橋反応解析部４による解析結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じた反応率から得られる値を、内部気泡の発生の有無に合わせて表示する。また、構造解析部５による構造解析の結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じた特性を、内部気泡の発生の有無に合わせて表示する。

【0045】

表示処理Ｓ５が終了すると、架橋反応シミュレーション装置１の動作が終了する。

【0046】

３．伝熱解析処理Ｓ１
３－１．カーボンブラックの質量比
次に、伝熱解析処理Ｓ１の詳細について、図２および図５～図８を参照して説明する。図２に示すように、下型１０Ａと上型１０Ｂとが組合された状態で、下キャビティ１００Ａと上キャビティ１００Ｂとによって形成される空間内に、外側接合部材１３と内側接合部材１４とが配置される。さらに、下キャビティ１００Ａと上キャビティ１００Ｂとによって形成される空間内に原料ポリマーが射出される。その後、金型１０内において、ポリマー部１２を構成する分子鎖の架橋反応が進行する。

【0047】

図５に示すように、熱板１１の温度は一定に設定されている。金型１０の温度は、熱板１１と同じ温度に設定されている。射出時におけるポリマー部１２の温度は、金型１０の設定温度よりも低い。このため、金型１０の温度は原料ポリマーが射出された後、一時的に下がる。その後、金型１０の温度は上昇し、熱板１１の温度と同じになる。ポリマー部１２の温度は、金型１０に射出された後、上昇し、十分に時間が経過した場合には金型１０の温度と同じになる。内側接合部材１４および外側接合部材１３の温度は、金型１０に配置される前の状態においては、外気温と同じであってもよいし、また、予熱されていても良い。外気温と同じである場合には、例えば冬季には０℃程度まで下がる場合もある。また、予熱される場合には、任意の温度に昇温してもよく、例えば１００℃付近まで昇温してもよい。内側接合部材１４および外側接合部材１３が金型１０に配置された後、内側接合部材１４および外側接合部材１３の温度は金型１０からの熱伝導により上昇し、原料ポリマーが射出された後は、さらに上昇し、十分に時間が経過した場合には金型１０の温度と同じになる。

【0048】

図５に示すように、ポリマー部１２の温度が熱板１１および金型１０の温度に達して平衡状態になるまでにある程度の時間が必要となる。平衡状態に達するまでの非定常な過程を計算するために下記の式（１）が用いられる。式（１）は、いわゆる、非定常熱伝導方程式である。

【0049】

【数1】

ここで、
ｔ：時間変数
ｘ：位置変数
Ｔ：温度
α：熱拡散率

【0050】

上記の熱拡散率αは、下記の式（２）で表される。

【0051】

【数2】

ここで、
λ：熱伝導率
ρ：密度
Ｃ：比熱

【0052】

熱拡散率αおよび熱伝導率λが、熱拡散率特性ＴＳに相当する。ただし、熱拡散率特性ＴＳは、熱拡散率αおよび熱伝導率λに限定されない。

【0053】

熱拡散率αが既知である場合には、熱拡散率αを直接用いて、非定常熱伝導方程式（１）によりポリマー部１２の温度を計算することができる。熱拡散率αが既知でない場合には、式（２）から、熱伝導率λ、密度ρおよび比熱Ｃを用いて熱拡散率αを計算し、算出されたαを用いて、式（１）によりポリマー部１２の温度を計算することができる。

【0054】

しかしながら、式（１）を用いてポリマー部１２の温度を計算する際に、一般的なＪＩＳ等においてゴムの熱拡散率αとして記載されている値を用いたり、または式（２）を用いる場合においても、ＪＩＳ等に熱伝導率λとして記載されている値を用いたりした場合には、図６に示すように、金型１０内に配置した対象ワークモデルＷＭのポリマー部１２の温度を実測したときの値と、式（１）に基づいてポリマー部１２の温度を予測したときの値が、十分な精度で一致しないという問題が生じた。発明者らが鋭意検討した結果、ＪＩＳ等に記載されたポリマー部１２の原料ポリマーの熱拡散率αまたは熱伝導率λの値と、実際のポリマー部１２の熱拡散率αまたは熱伝導率λの値とが、異なっていたことが原因であることが分かった。

【0055】

そこで、測定対象となる対象ワークモデルＷＭのポリマー部１２の温度を実測し、ポリマー部１２の熱拡散率αまたは熱伝導率λを同定することにより、測定対象となる対象ワークモデルＷＭのポリマー部１２の熱拡散率αまたは熱伝導率λを得た。

【0056】

図７には、対象ワークモデルＷＭのポリマー部１２の温度の実測値と、上記のようにして得られた熱拡散率αまたは熱伝導率λの実測値を用いて計算した温度と、を示す。すると、ポリマー部１２の温度について、実測値と予測値とが高い精度で一致することが分かった。

【0057】

しかしながら、対象ワークモデルＷＭを構成するポリマー部１２の熱拡散率αまたは熱伝導率λは、ポリマー部１２の原料ポリマーの種類や、原料ポリマーに添加される添加材料の配合量または種類等により、変化する。このため、異なる配合のすべてのポリマー部１２について、熱拡散率αまたは熱伝導率λを実測することは煩雑である。

【0058】

そこで、発明者らは、ポリマー部１２に添加される添加材料のうち、カーボンブラックに着目した。ポリマー部１２を構成するゴムは、熱が伝わりにくい材料である。一方、カーボンブラックは熱が伝わりやすい材料である。このため、ポリマー部１２の熱拡散率αまたは熱伝導率λは、ポリマー部１２を構成する原料ポリマーに対するカーボンブラックの質量比と、関連性が高いと考えられる。カーボンブラックとしては、ファーネスブラック、アセチレンブラック、ケッチェンブラック等が挙げられる。

【0059】

図８に、複数種の原料ポリマーを用い、異なるカーボンブラックの質量比で実験した場合における、原料ポリマーに対するカーボンブラックの質量比と、ポリマー部１２の熱伝導率λとの関係性を示す。図８においては、原料ポリマーに対するカーボンブラックの質量比は、質量部とした。ただし、原料ポリマーに対するカーボンブラックの質量比は、質量％としてもよい。

【0060】

図８に示すように、ポリマー部１２の熱伝導率λは、原料ポリマーに対するカーボンブラックの質量比に対して、線形の関係を有する。また、ポリマー部１２の熱伝導率λは、原料ポリマーの種類に依存しない関係であることが分かった。これにより、原料ポリマーの種類によらず、カーボンブラックの配合量に基づいて、ポリマー部１２の熱伝導率λを精度よく予測することが可能になった。また、式（２）を用いて熱伝導率λから熱拡散率αを算出することができるので、熱拡散率αについても精度よく予測することが可能となった。

【0061】

３－２．接触熱伝達率ＣＨ
図２に示すように、金型１０内に外側接合部材１３が配置された状態で、金型１０の内面と、外側接合部材１３との間には、外側隙間１５Ａが形成される場合がある。この場合、熱は、金型１０の内面から、外側隙間１５Ａ内の空気を介して、外側接合部材１３に伝わる。伝熱解析処理Ｓ１においては、金型１０の内面と外側接合部材１３とが接触した状態を仮定し、金型１０および外側接合部材１３において外側隙間１５Ａを考慮する面を選択し、外側隙間１５Ａの大きさに応じて、接触熱伝達率ＣＨを設定する。接触熱伝達率ＣＨは、テストピースを用いた実験において、外側隙間１５Ａの大きさを変化させたときの温度挙動を測定し、それぞれの温度挙動を再現可能な値を解析で求めたものが採用される。

【0062】

また、金型１０内に内側接合部材１４が配置された状態で、金型１０の内面と内側接合部材１４との間には、内側隙間１５Ｂが形成される場合がある。この場合、熱は、金型１０の内面から、内側隙間１５Ｂ内の空気を介して、内側接合部材１４に伝わる。伝熱解析処理Ｓ１においては、金型１０の内面と内側接合部材１４とが接触した状態を仮定し、金型１０および内側接合部材１４において内側隙間１５Ｂを考慮する面を選択し、内側隙間１５Ｂの大きさに応じて、接触熱伝達率ＣＨを設定する。接触熱伝達率ＣＨは、テストピースを用いた実験において、内側隙間１５Ｂの大きさを変化させたときの温度挙動を測定し、それぞれの温度挙動を再現可能な値を解析で求めたものが採用される。

【0063】

上記したように、記憶部２は、成形型モデルＭＭに対象ワークモデルＷＭを配置した状態で、成形型モデルＭＭと金具との隙間１５Ａ，１５Ｂにおける接触熱伝達率ＣＨを記憶している。

【0064】

伝熱解析処理Ｓ１においては、解析部３２は、成形型モデルＭＭに対象ワークモデルＷＭを配置した状態で、ポリマー熱拡散率決定部３１により決定されたポリマー熱拡散率、記憶部２に記憶された温度条件ＴＣおよび接触熱伝達率ＣＨを用いて、伝熱解析を行う構成となっている。これにより、精度よく、ポリマー部１２の温度を予測することができる。

【0065】

３－３．空気熱伝達係数ＡＣ
成形型モデルＭＭから対象ワークモデルＷＭが脱型された後も、対象ワークモデルＷＭのポリマー部１２の温度はすぐには下がらないので、ポリマー部１２において架橋反応が進行する。そこで、本形態においては、成形型モデルＭＭから対象ワークモデルＷＭが脱型された後においても、ポリマー部１２の温度を予測するために、条件入力部３０は、成形型モデルＭＭから脱型された対象ワークモデルＷＭの周囲の気温を含む外気条件ＯＣが入力されるようになっている。さらに記憶部２は、成形型モデルＭＭから脱型された対象ワークモデルＷＭの周囲の空気の熱伝達係数である空気熱伝達係数ＡＣを記憶する。

【0066】

本形態の解析部３２は、対象ワークモデルＷＭのポリマー部１２を成形型モデルＭＭ内にて架橋反応させた後に成形型モデルＭＭから対象ワークモデルＷＭを脱型した状態で、ポリマー熱拡散率決定部３１により決定されたポリマー熱拡散率、記憶部２に記憶された空気熱伝達係数ＡＣおよび外気条件ＯＣを用いて、伝熱解析を行う。これにより、成形型モデルＭＭから対象ワークモデルＷＭが脱型された後の状態において、ポリマー部１２の温度を精度よく予測することができる。

【0067】

図９に示すように、成形型モデルＭＭに対象ワークモデルＷＭを配置した状態について解析メッシュを作成する。解析メッシュのサイズは、予め定められていても良いし、条件入力部３０から入力されても良い。解析メッシュの作成には、成形型モデルＭＭおよび対象ワークモデルＷＭの形状の特性（例えば対称性）を考慮する。メッシュタイプとしては、８節点ソリッド要素を採用した。ただし、メッシュタイプは任意のモデルを採用可能であり、例えば４節点ソリッド要素を用いても良い。

【0068】

解析メッシュを用いることにより、対象ワークモデルＷＭの各節点における温度を予測することができる。つまり、対象ワークモデルＷＭのポリマー部１２の温度分布を予測することができる。特に、各時刻における、ポリマー部１２の各節点の温度を予測することができる。

【0069】

４．架橋反応解析処理Ｓ２
次に架橋反応解析処理Ｓ２について説明する。架橋反応解析処理Ｓ２については、まず反応率について説明する。次に従来技術について説明し、続いて従来技術の課題について説明する。その後、本形態の架橋反応解析処理について説明する。
４－１．反応率について
架橋反応の反応率は、以下のようにして測定される。図１０に示すように、本形態に係るゴムの反応率については、ＪＩＳ Ｋ６３００－２に準拠し、加硫試験機７０（架橋反応試験機の一例）を使用して、得られたトルクと反応時間との関係性から、反応率（加硫度）を算出する。加硫試験機７０としては、例えばキュラストメーター（登録商標）を用いることができる。また、反応率（加硫度）を考慮した振動特性を出力する場合には、例えばゴム加工性試験機（ＲＰＡ）を用いて、所望の振動条件での反応挙動（加硫挙動）を測定することができる。ＲＰＡは、広範な範囲で振動数、振幅の設定が可能に構成されている。

【0070】

図１０に示すように、加硫試験機７０は、キャビティ７１を有する下ダイ７２と上ダイ７３の間に、ロータ７４の先端に取付けられたディスク７５が配置されている。下ダイ７２と上ダイ７３は所定の温度に設定することができる。ロータ７４は、ねじり角が＋３°～－３°、または＋１°～－１°に選択可能になっている。キャビティ７１とディスク７５との間に測定試料７６が配置され、一定温度下でディスク７５から回転トルクを与えながら架橋反応させ、トルクの変化から架橋反応の進行を検出するようになっている。

【0071】

図１１に測定結果の一例を示す。図１１（ａ）は、所定の反応温度における、反応時間に対するトルクの変化量のグラフ（架橋反応曲線ＣＣ）である。一方、反応率は、最小トルクを反応率０％とし、最大トルクを反応率１００％と定義される。図１１（ｂ）に、所定の反応温度における、反応時間に対する反応率の変化を示すグラフ（基準反応曲線ＲＣ）を示す。本形態に係るゴムの場合、反応率は加硫度と呼称される。加硫度は、加硫ゴムの物性（弾性率、伸び、引張強さ、硬さなど）を指標としてみた場合の加硫の程度と定義される。なお、図１１（ａ）において、反応時間０から最小トルクを示す反応時間までの領域は、ポリマー部１２の架橋反応が未反応であり、可塑的な領域となっている。このため、図１１（ｂ）においては、当該領域においては、反応率は０％とする。

【0072】

４－２．反応率の算出方法
次に、反応率の算出方法について説明する。上記した非特許文献１においては、アレニウスの式（３）を変形することにより得られた等価反応量に係る式（４）を用いて、架橋反応の等価反応量を推定している。なお、等価反応量に係る式（４）は、反応温度Ｔ、反応時間ｔにおける反応量が、ある基準温度Ｔ_０、基準時間ｔ_０における反応量の何倍になっているかを表す。なお、式（４）は、ゴムの加硫反応においては、等価加硫の式と呼称される。この等価加硫の式を用いることにより、後述する反応戻り期における反応率を予測することができる。後に詳述するが、本形態においては、反応戻り期における予測精度を、単に等価加硫の式を用いた場合に比べて、さらに向上させることが可能となる。なお、一般的な解析ソフトにおいて、ゴムやエポキシ等の熱硬化反応の挙動を予測する際に用いられる式（例えば、Ｋａｍａｌの反応速度式）によっては、反応戻り期における反応率を予測することはできない。

【0073】

【数3】

ここで、
ｋ：反応速度
Ａ：頻度因子
Ｅ：活性化エネルギー
Ｒ：ガス定数
Ｔ：反応温度

【0074】

【数4】

ここで、
Ｕ：等価反応量
Ｅ：活性化エネルギー
Ｒ：ガス定数
Ｔ：反応温度
Ｔ_０：基準温度
ｔ：反応時間

【0075】

アレニウスの式（３）から、以下のように式変形して、等価反応量の式（４）が導出される。

【0076】

反応時間ｔにおける架橋反応の反応量Ｚは、式（５）で表される。

【0077】

【数5】

ここで、
Ｚ：反応量
ｔ：反応時間

【0078】

反応量Ｚは、アレニウスの式（３）と、上記の式（５）から、下記の式（６）のようになる。

【0079】

【数6】

【0080】

ここで、ある基準になる反応温度Ｔ_０で、反応時間ｔ_０のときの反応量Ｚ_０を考えると、Ｚ_０は、下記の式（７）で表される。

【0081】

【数7】

【0082】

反応温度Ｔ_０および反応時間ｔ_０における反応量Ｚ_０に対する、反応温度Ｔおよび反応時間ｔにおける反応量Ｚの比Ｕをとると、比Ｕは下記の式（８）で表される。

【0083】

【数8】

【0084】

式（８）に、式（６）および式（７）を代入すると比Ｕについて式（９）を得る。

【0085】

【数9】

【0086】

式（９）において、ｔ_０＝１のとき、式（９）は上記した式（４）のように表される。

【0087】

反応温度Ｔが刻々と変わるときは、対象時刻における温度、上記の式（４）、傾き係数ＳＣに基づいて、対象時刻の微小時間における等価反応量増加量を算出し、等価反応量増加量に基づいて、架橋反応開始から対象時刻までの等価反応量積算値を算出し、架橋反応開始から対象時刻までの前記等価反応量積算値に基づいて、対象時刻における反応率を算出する。

【0088】

式（４）に基づいて、対象時刻の微小時間Δｔにおける等価反応量増加量ΔＵ_ｉは、下記の式（１０）により表される。

【0089】

【数10】

ここで、
ΔＵ_ｉ：対象時刻の微小時間における等価反応量増加量
Δｔ：対象時刻の微小時間

【0090】

反応温度Ｔが刻々と変わるときには、下記の式（１１）のように単位時間毎に累積することにより、等価反応量Ｕを求めることができる。

【0091】

【数11】

ここで、
Ｕ：等価反応量
Ｅ：活性化エネルギー
Ｒ：ガス定数
Ｔ：反応温度
Ｔ_０：基準温度
Δｔ：対象時刻の微小時間

【0092】

上記した式（４）は、上記の式（１１）のうち、Σの内部の式であって、累積される個々の等価反応量に対応する。

【0093】

式（４）または式（１１）を計算するためには、（－Ｅ／Ｒ）の値が必要となる。従来技術においては、（－Ｅ／Ｒ）を下記のようにして求めていた。

【0094】

アレニウスの式（３）の両辺の自然対数をとることにより、下記の式（１２）を得る。

【0095】

【数12】

ここで、
ｋ：反応速度
Ａ：頻度因子
Ｅ：活性化エネルギー
Ｒ：ガス定数
Ｔ：反応温度

【0096】

上記した反応量について、所定の反応率γにおける反応時間ｔ_γと、反応速度kとの間に下記の式（１３）が成り立つと考える。反応率γは、例えば５％、９０％等、０％～１００％の任意の値を適宜に採用することができる。

【0097】

【数13】

ここで、
Ｂ：任意の定数
ｔ_γ：所定の反応率γにおける反応時間

【0098】

式（１３）の両辺の自然対数をとり、式（１２）に代入して式変形することにより、式（１４）を得る。

【0099】

【数14】

【0100】

式（１４）から、ｌｎ（１／ｔ_γ）を縦軸に、（１／Ｔ）を横軸にして、反応率がα％に到達する時間と、各測定温度をプロットすると、直線に近似することが可能であり、この直線の傾きが（－Ｅ／Ｒ）になることがわかる（図１２参照）。

【0101】

実際の反応においては、図１３に示すように、反応時間の経過とともに温度が変化する。このような場合、等価反応量の計算には式（４）および式（１１）を用い、下記のように計算する。

【0102】

まず、図１３に示すように、反応時間を所定の間隔（例えば、１分）で区切る。次に、所定間隔ごとの温度を読取る。

【0103】

次に、所定間隔当たりの等価反応量を、式（４）を用いて計算する。式（４）において、所定間隔Δｔを代入して計算する。所定間隔として、例えば、Δｔ＝１とすることにより、計算を簡略化できる。

【0104】

次に、式（４）によって計算された所定間隔Δｔ当たりのΔＵ_ｉを足し合わせる。具体的な計算方法について、図１４を参照して説明する。

【0105】

図１４に示すように、反応時間ｔ_ｉ－１、反応温度Ｔ_ｉ－１における等価反応量Ｕ_ｉ－１を式（４）で計算する。同様にして、反応時間ｔ_ｉ、反応温度Ｔ_ｉにおける等価反応量Ｕ_ｉを算出する。

【0106】

次に、Ｕ_ｉ－１×Δｔと、Ｕ_ｉ×Δｔとを足す。これにより、反応時間ｔ_ｉ－１～ｔ_ｉにおける、等価反応量の積算値が得られる。なお、本形態ではΔｔ=１／２ｔである。同様にして、式（１１）に従い、ｔ＝１～ｔ＝ｎまで所定時間当たりの等価反応量ΔＵ_ｉを積算することにより、等価反応量Ｕを得ることができる。

【0107】

続いて、等価反応量Ｕから反応率を算出する方法について説明する。まず、反応温度Ｔ_０、反応時間tにおける反応量Zは式（５）で表される。次に、式（５）から、反応温度Ｔ_０で、反応時間ｔ_０のときの反応量Ｚ_０は、下記の式（１５）で表される。ここで、式（５）と式（１５）の反応温度は同一の値Ｔ_０であるため、式（３）からそれぞれのkは同一の値となる。

【0108】

【数15】

【0109】

上記した式（８）に、上記した式（５）および式（１５）を代入すると、反応速度ｋが約分されて下記の式（１６）を得る。

【0110】

【数16】

【0111】

式（１６）において、ｔ_０＝１とすると、式（１６）は下記の式（１７）のように表される。

【0112】

【数17】

式（１７）によれば、ｔ_０＝１としたとき、反応時間ｔを、等価反応量Ｕと読み替えることができる。これにより、図１１（ｂ）に示した、反応温度Ｔ_０での、反応率の、反応時間に対する変化を示すグラフにおいて、反応時間ｔを等価反応量Ｕと読み替えることができる。図１５に、縦軸を反応率とし、横軸を等価反応量としたグラフを示す。このグラフに基づいて、計算によって得られた等価反応量Ｕから反応率を算出することができる。

【0113】

４－３．反応率算出方法の課題
しかしながら、上記した方法に基づいて反応率を推定した場合、実測値との差が大きく、十分な精度が得られないという問題があった。

【0114】

図１６に、反応時間に対する反応率の変化を示す。式（４）を用いるに際して、基準温度を１６０℃としたので、反応温度１６０℃の場合については、実測値と予測値とは精度よく一致している。

【0115】

しかし、反応温度１７０℃および１８０℃については、実測値と予測値とが十分に一致していない。概ね、実測値よりも予測値の方が大きな値となっている。なお、反応率の極大値が一致しているのは、上記したように、反応率の極大値を１００％と定義したためである。

【0116】

図１７を参照して、従来技術において、実測値と予測値とが十分に一致しない理由について説明する。図１７には、反応時間に対する反応率の変化を表すグラフを示す。本形態に係る架橋反応は、以下の過程を経て進行する。

【0117】

反応が開始した反応進行初期においては、徐々に反応率が増加する。ある程度、時間が経過して反応進行促進期に達すると、反応率が急激に増加する。その後、反応進行後期になると、反応率の増加量は緩やかになり反応率は極大値に達する。反応率が極大値に達した後は、戻り反応が進行する反応戻り期となる。この反応戻り期においては、反応率は漸減する。

【0118】

上記した技術においては、反応進行初期から反応戻り期に至るまでを、いわば、１つの架橋反応として捉えていた。このため、１つの活性化エネルギーＥを用いて、反応率の予測を行っていた。この結果、反応率の実測値と、予測値とが十分に一致しなかったと考えられる。

【0119】

４－４．本形態の架橋反応解析処理
本形態においては、反応の進行状態に対応して複数の活性化エネルギーを算出し、得られた活性化エネルギーに基づいて等価反応量を予測する構成とした。この考え方に基づいて、式（４）を下記の式（１８）のように変形した。

【0120】

【数18】

ここで、
Ｅ_Ｎ：Ｎ段階に分けられた、反応の進行段階ごとの活性化エネルギー

【0121】

図１７には、反応の進行状態を、反応進行初期、反応進行促進期、反応進行後期および反応戻り期の４段階に分ける例を示す。ただし、反応の進行状態は２つ～３つ、または５つ以上の任意の段階に分けてもよい。また、反応進行期および反応戻り期を含む全期間を等間隔に区分してもよい。等間隔に区分する指標として、架橋反応の全期間の進行度を基準として、例えば架橋反応の進行度を５％ごとに区分してもよい。ただし、区分の指標は５％に限られず、１％～４％または、６％以上でも良い。

【0122】

本形態においては、反応率が０％～５％までを反応進行初期とし、反応率が５％～５０％までを反応進行促進期とし、反応率が５０％～１００％までを反応進行後期とし、反応進行後期よりも後の段階において、反応率が１００％から９０％にまで減少した範囲を反応戻り期とした。

【0123】

各反応段階において、上記した式（１４）を用いて、活性化エネルギーＥ_Ｎを含む項（Ｅ_Ｎ／Ｒ）の値を求めた。図１７に示すように、反応進行初期（反応率が０％～５％）における（Ｅ_１／Ｒ）は約１２０００であり、反応進行促進期（反応率が５％～５０％）における（Ｅ_２／Ｒ）は約９０００であり、反応進行後期（反応率が５０％～１００％）における（Ｅ_３／Ｒ）は約１３０００であり、反応戻り期（反応率が１００％～９０％）における（Ｅ_４／Ｒ）は約１６０００であった。ただし、（Ｅ_Ｎ／Ｒ）の数値は上記の値に限定されない。なお、反応進行初期から反応戻り期までを１つの反応段階として捉えたときの（Ｅ／Ｒ）は１１０００であった。

【0124】

図１８に、上記のようにして得られた活性化エネルギーＥ_１～Ｅ_４を用いて反応率を予測した結果を、実測値と合わせて示す。基準温度である１６０℃において、実測値と予測値とが一致するだけでなく、反応温度１７０℃および１８０℃においても、実測値と予測値とが精度よく一致した。これにより、精度よく反応率を予測することができる。

【0125】

本形態においては、傾き係数ＳＣは、架橋反応の反応率がピークに到達するまでの反応進行期と、反応率がピークを越えた反応戻り期とで、異なる値に設定される。また、本形態においては、傾き係数ＳＣは、反応進行期より反応戻り期の方が小さな値に設定される。

【0126】

また、本形態においては、傾き係数ＳＣは、反応進行初期、反応進行促進期、反応進行後期および反応戻り期のそれぞれで、異なる値に設定される。また、本形態においては、反応促進期における傾き係数ＳＣが、反応進行初期および反応進行後期より大きな値に設定される。

【0127】

また、傾き係数ＳＣは、架橋反応の全期間が等間隔に区分された場合には、各区分に応じた値に設定される。

【0128】

４－５．架橋反応の反応率の表現方法
図１９を参照して、本形態における架橋反応の反応率の表現方法について説明する。図１９は、反応時間に対する、架橋反応の反応率の変化を示すグラフである。実線で示されたグラフは、反応率の実測値である。

【0129】

実線で示す架橋反応の反応率は、極大値に達した後は、戻り反応が進行することによって反応率は漸減する。このため、例えば、反応率が極大値に達する前に反応率が８０％なる点と、反応率が極大値に達した後に反応率が８０％になる点と、が存在することになる。この場合、単に反応率の数字だけで比較すると、反応率が極大値に達する前の８０％の数値と、反応率が極大値に達した後の８０％の数値を区別することができない。

【0130】

そこで、本形態においては、架橋反応の反応率は、架橋反応の反応率がピークに到達するまでの反応進行期において、架橋反応の開始時を０％とし、架橋反応の反応率の増加度を０％～１００％までの範囲で定義し、反応率がピークを越えた反応戻り期において、反応率のピークからの反応率の低下度を１００％に加算した値として定義する。

【0131】

図１９のグラフで説明すると、反応率が１００％の点から横軸に平行に延びる直線を対称軸として、反応率がピークを越えた反応戻り期のグラフを、矢線Ｂで示すように反転させる。これにより、反応率が極大値に達した後の反応率８０％の点は、反応率１２０％と表現することができる。この結果、反応率が極大値に達する前における反応率と、反応率が極大値に達した後における反応率と、を明確に識別することができる（図２０参照）。

【0132】

４－６．シミュレーションにより算出される値の決定方法
次に、シミュレーションにより算出される値の決定方法について説明する。本形態においては、傾き係数ＳＣおよび反応率は、データベースおよび関数の、一方または双方を用いて決定される。データベースを用いることにより、推定値の精度を向上させることができる。一方、関数を用いることにより、演算速度を向上させることができる。

【0133】

（１）データベースを用いる方法
図１に示すように、記憶部２は、ポリマー部１２の架橋反応の反応率と傾き係数ＳＣとの対応関係を定義した第一関係データマップＤＭ１を記憶する。反応率算出処理部４１は、前回時刻の反応率および第一関係データマップＤＭ１を用いて、傾き係数ＳＣを決定する。

【0134】

記憶部２は、等価反応量積算値と対象時刻における反応率との対応関係を定義した第二関係データマップＤＭ２を記憶する。反応率算出処理部４１は、対象時刻における等価反応量積算値および第二関係データマップＤＭ２を用いて、対象時刻における反応率を決定する。

【0135】

第二関係データマップＤＭ２は、対象ワークモデルＷＭのポリマー部１２に相当する試験対象ポリマー材料を用いて、加硫試験機７０により測定された基準反応温度において架橋反応時間と試験対象ポリマー材料に発生したトルクとの関係に基づく基準反応曲線ＲＣより設定されている。

【0136】

（２）関数を用いる方法
図１に示すように、記憶部２は、ポリマー部１２の架橋反応の反応率と傾き係数ＳＣとの対応関係を定義した第一関数Ｆ１を記憶する。反応率算出処理部４１は、前回時刻の反応率および第一関数Ｆ１を用いて、傾き係数ＳＣを決定する。

【0137】

第一関数Ｆ１は、ポリマー部１２の架橋反応の進行度に応じて設定された複数の反応率区分に応じて、異なる関数に設定されている。架橋反応の進行度の区分数は任意であり、１つまたは２つ以上の複数に区分してもよい。関数については、線形補完を含むｎ次関数、スプライン曲線等、任意の関数を適宜に選択できる。本形態では、例えば、６次関数が好適に用いられる。なお、第一関数Ｆ１の区分数は、Ｎ段階に分けられた反応の進行段階ごとの活性化エネルギーＥ_Ｎの値を求める際の区分とは独立して設定することができる。

【0138】

記憶部２は、等価反応量積算値と対象時刻における反応率との対応関係を定義した第二関数Ｆ２を記憶する。反応率算出処理部４１は、対象時刻における等価反応量積算値および第二関数Ｆ２を用いて、対象時刻における反応率を決定する。

【0139】

第二関数Ｆ２は、対象ワークモデルＷＭのポリマー部１２に相当する試験対象ポリマー材料を用いて、加硫試験機７０により測定された基準反応温度において架橋反応時間と前記試験対象ポリマー材料に発生したトルクとの関係に基づく基準反応曲線ＲＣより設定されている。

【0140】

第二関数Ｆ２は、ポリマー部１２の架橋反応の進行度に応じて設定された複数の反応率区分に応じて、異なる関数に設定されている。架橋反応の進行度の区分数は任意であり、１つまたは２つ以上の複数に区分してもよい。関数については、線形補完を含むｎ次関数、スプライン曲線等、任意の関数を適宜に選択できる。本形態では、例えば、６次関数が好適に用いられる。なお、第二関数Ｆ２の区分数は、Ｎ段階に分けられた反応の進行段階ごとの活性化エネルギーＥ_Ｎの値を求める際の区分とは独立して設定することができる。

【0141】

図２１に、第二関数Ｆ２について、架橋反応の進行度を４つに区分した場合の関数を例示する。本形態では、反応進行初期、反応進行促進期、反応進行後期および反応戻り期に区分されている。各区分において、反応率と、反応時間とが、互いに異なる関数により近似されている。

【0142】

次に、対象ワークモデルＷＭの解析メッシュを作成する。解析メッシュのサイズは、予め定められていても良いし、条件入力部３０から入力されても良い。解析メッシュの作成には、対象ワークモデルＷＭの形状の特性（例えば対称性）を考慮して、図２２に示すような解析メッシュを作成する。メッシュタイプとしては、８節点ソリッド要素を採用した。ただし、メッシュタイプは任意のモデルを採用可能であり、例えば４節点ソリッド要素を用いても良い。

【0143】

解析メッシュを用いて、対象ワークモデルＷＭのポリマー部１２のみの各時刻の各節点の温度から、各節点の反応率を算出することができる。

【0144】

５．構造解析処理Ｓ３
５－１．初期処理
次に、構造解析処理Ｓ３について説明する。構造解析部５は、対象ワークモデルＷＭの解析メッシュを作成する。対象ワークモデルＷＭの解析メッシュは、架橋反応解析処理Ｓ２で作成した解析メッシュと同じでも良いし、異なっていても良い。本形態では構造解析処理Ｓ３で用いられる解析メッシュは、架橋反応解析処理Ｓ２で用いられた解析メッシュと同じものを使用した。

【0145】

５－２．温度および反応率取得
次に、構造解析部５は、伝熱解析部３によって解析された対象ワークモデルＷＭのポリマー部１２の温度分布についての解析結果を取得する。また、構造解析部５は、架橋反応解析部４によって解析された対象ワークモデルＷＭのポリマー部１２の反応率についての解析結果を取得する。

【0146】

５－３．弾性率予測
次に、構造解析部５は、対象ワークモデルＷＭのポリマー部１２の、各時刻における各節点の弾性率を算出する。まず、架橋反応の反応率が１００％のテストピースを作成し、このテストピースの弾性率が実測される。次に、対象とする反応率におけるトルクを、反応率１００％におけるトルクで除した商に基づいて、補正係数が算出される。上記のトルクは、加硫試験機７０で測定可能なトルクである。構造解析部５は、実測された弾性率に、算出された補正係数を乗じることにより、各時刻における各節点の弾性率を算出する。

【0147】

５－４．特性予測
次に、構造解析部５は、対象ワークモデルＷＭの特性を予測する。有限要素解析に用いる入力条件は、各時刻の各節点の温度、各時刻の各節点の弾性率、並びに、外側接合部材１３および内側接合部材１４の境界条件である。

【0148】

入力条件が設定された後、対象ワークモデルＷＭの軸線Ａに対して直角な方向の力または歪を、対象ワークモデルＷＭに付与する。

【0149】

構造解析部５は、対象ワークモデルＷＭの全体の弾性率を出力する。本形態では、対象ワークモデルＷＭの軸線Ａに直角な方向に静的に加えられた力に対する弾性率が、特性として予測される。ただし、予測される特性は、軸線Ａに直角な方向の力または歪に対する弾性率に限られず、軸線Ａに直角な方向に振動する力に対する弾性率、軸線Ａに沿う方向に静的に加えられた力または歪に対する弾性率、軸線Ａに沿う方向に振動する力に対する弾性率等、任意の物性値を特性として予測することができる。

【0150】

６．内部気泡推定処理Ｓ４
次に内部気泡推定処理Ｓ４について説明する。本形態の一例であるゴム製品は、成形加硫して作られる。このため、元々ゴムに溶解していた気体、または加硫反応によって発生した気体が、加硫の高温高圧条件下でゴム中に溶解した状態になっている。金型１０を開放すると、ゴムへかかる圧力が低下し、ゴムへの気体の溶解度が低下するため、加硫が十分に進行していない状態では、ゴムにより構成されるポリマー部１２の中に内部気泡８０が発生する場合がある（図２３参照）。

【0151】

金型１０内にゴムを配置した状態で加硫時間を長くし、加硫を進行させると内部気泡８０が発生しなくなる。実際の製品の加硫においては、熱源である金型１０からの距離が異なる製品内部の点において、温度上昇履歴の違いから、受ける総熱量が異なり、同じ加硫時間であっても、加硫進行状況が異なる。そのため、最遅加硫部分において、内部気泡８０の発生が見られなくなる時間まで金型１０内で加熱して加硫する必要がある。この最遅加硫部分において内部気泡８０の発生が見られなくなるまでの加硫時間で成形した時の、金型１０を開放した時点での、最遅加硫部分のゴムの加硫度をブローポイント加硫度（以下、ブローポイント）という（図２４参照）。

【0152】

内部気泡８０を抑制するためには、金型１０内での加熱時間はできるだけ長くした方が好ましい。しかし、金型１０内での加熱時間を過剰に伸ばすと製品の製造効率が低下する。このため、製品の加硫時間を決めるためには、上記のブローポイントを推定することが重要となる。

【0153】

図２５(ａ)に、複数種のポリマーに対して、複数種のフィラーを、異なる配合量で添加した場合の、ブローポイント反応率を示す。ブローポイント反応率は、ブローポイントにおけるポリマー部１２の架橋反応の反応率である。

【0154】

原料ポリマーは、Ａ～Ｈの８種類のポリマーを用いた。サンプルＡ１～Ａ３は、原料ポリマーＡを用いて、フィラーの配合量を変化させたものである。同様に、サンプルＢ１～Ｂ３は、原料ポリマーＢに対して、フィラーの配合量を変化させたものであり、サンプルＣ１～Ｃ３は、原料ポリマーＣに対して、フィラーの配合量を変化させたものである。サンプルは市販品を用いたため、フィラーの種類および配合量について正確には分からない。

【0155】

図２５（ａ）に示すように、ブローポイント反応率は、原料ポリマーの種類、フィラーの配合量によって大きく異なる。このため、ブローポイント反応率を用いて、ポリマー部１２のブローポイントを予測することは困難であることが分かった。

【0156】

図２６に示すように、内部気泡８０が膨張しようとする圧力Ｐ１と、ポリマー部が内部気泡８０を抑え込もうとする圧力Ｐ２とが、Ｐ１＜Ｐ２となったときに、内部気泡８０がゴムによって抑え込まれて消失すると考えられる。そこで、ブローポイント反応率を、当該ブローポイント反応率に対応する、加硫試験機７０により測定されたトルクに変換して、各サンプルで比較した。

【0157】

図２５（ｂ）に、複数種のポリマーに対して、複数種のフィラーを、異なる配合量で添加した場合の、ブローポイントトルクを示す。ブローポイントトルクは、上記したように、ブローポイント反応率を、当該ブローポイント反応率に対応する、加硫試験機７０により測定されたトルクに変換したものである。

【0158】

図２５（ｂ）に示すように、原料ポリマーの種類や、フィラーの配合量によらず、内部気泡８０の発生が見られなくなるトルク（ブローポイントトルク）は、ほぼ同等の値となった。つまり、容易に測定可能なトルクを用いることにより、ブローポイントを容易に予測可能となる。図２５（ｂ）の縦軸に、閾値を記載した。内部気泡推定部６の推定部６２は、対象ワークモデルＷＭが脱型されるときにおけるトルクが閾値を超えていない場合に、対象ワークモデルＷＭのポリマー部１２において内部気泡８０が発生したと推定する。

【0159】

閾値は、複数種の原料ポリマーに対して統一した値に設定されている。また、閾値は、ポリマー部１２のブローポイントにおける架橋反応の反応率よりも大きな反応率に対応するトルクである。ただし、閾値は、複数のグループを設定し、グループごとに閾値を設定してもよい。複数のグループは任意の基準で区分することができる。

【0160】

反応率取得部５１は、ポリマー部１２の部位ごとの反応率を取得し、トルク算出部６１は、ポリマー部１２の部位ごとのトルクを算出し、推定部６２は、ポリマー部１２の部位ごとの内部気泡８０の発生を推定する。

【0161】

７．表示処理Ｓ５
次に、表示処理Ｓ５について説明する。図１に示すように、表示部７は、架橋反応解析部４の解析結果に基づいて、反応時間に応じた反応率から得られる値を表示する。

【0162】

さらに、表示部７は、構造解析部５の解析結果に基づいて、ポリマー部１２の架橋反応に用いられる成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じた特性を表示する。

【0163】

さらに、表示部７は、内部気泡推定部６の推定結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じて、内部気泡８０の発生の有無を表示する。

【0164】

さらに、表示部７は、架橋反応解析部４の解析結果に基づいて、成形型モデルＭＭの温度と、架橋反応開始から成形型モデルＭＭの脱型までの型内反応時間と、に応じた反応率から得られる値を、内部気泡８０の発生の有無に合わせて表示する。ただし、表示部７が表示する内容は上記に限定されない。

【0165】

７－１．表による表示方法
（１）反応率を表示する形態
表１～表４を参照して、架橋反応の反応率を、表を用いて表示する形態について説明する。

【0166】

【表1】

【0167】

表１は、対象ワークモデルＷＭのうち、架橋反応の反応率が最小である要素における架橋反応の反応率について、金型温度と、反応時間との関係についてまとめた表である。

【0168】

表１に記載された反応時間は、金型１０の内部に対象ワークモデルＷＭが配置された状態における反応時間を表す。また、表１に記載された反応率の数値は、金型１０の内部に対象ワークモデルＷＭが配置された状態における反応時間と、対象ワークモデルＷＭが金型１０から脱型された後の時間と、の双方の温度履歴に基づいて算出されたものである。対象ワークモデルＷＭが金型１０から脱型された後も架橋反応が進行することから、上記のような考慮が必要となる。以下の説明において、表２～表４に記載された反応時間および反応率についても同様である。

【0169】

表１のうち、右上部において太い罫線で囲まれた領域は、反応率が１１５％より大きい値となっている。この領域においては、架橋反応の戻り反応が過剰に進行しているので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、１１５％と異なる値を閾値としても良い。

【0170】

表１の左下部において太い罫線で囲まれた領域は、反応率が９５％よりも小さな値となっている。この領域においては、架橋反応が十分に進行していないので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、９５％と異なる値を閾値としても良い。

【0171】

表１のうち、太い罫線で囲まれた、右上部の領域および左下部の領域を除く領域は、架橋反応の反応率の観点からは好ましい領域である。なお、表１には、後述する内部気泡８０に関する情報は記載されていない。

【0172】

表２は、対象ワークモデルＷＭのすべての要素の架橋反応の平均値について、金型温度と、反応時間との関係についてまとめた表である。

【0173】

【表2】

【0174】

表２のうち、右上部において太い罫線で囲まれた領域は、反応率が１１５％より大きい値となっている。この領域においては、架橋反応の戻り反応が過剰に進行しているので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、１１５％と異なる値を閾値としても良い。

【0175】

表２の左下部において太い罫線で囲まれた領域は、反応率が９５％よりも小さな値となっている。この領域においては、架橋反応が十分に進行していないので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、９５％と異なる値を閾値としても良い。

【0176】

表２のうち、右上部の領域と、左下部の領域と、を除く領域は、架橋反応の反応率の観点からは好ましい領域である。なお、表２には、後述する内部気泡８０に関する情報は記載されていない。

【0177】

表３は、対象ワークモデルＷＭのうち、架橋反応の反応率が最大である要素における架橋反応の反応率について、金型温度と、反応時間との関係についてまとめた表である。

【0178】

【表3】

【0179】

表３のうち、右上部において太い罫線で囲まれた領域は、反応率が１１５％より大きい値となっている。この領域においては、架橋反応の戻り反応が過剰に進行しているので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、１１５％と異なる値を閾値としても良い。

【0180】

表３の左下部において太い罫線で囲まれた領域は、反応率が９５％よりも小さな値となっている。この領域においては、架橋反応が十分に進行していないので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、９５％と異なる値を閾値としても良い。

【0181】

表３のうち、右上部の領域と、左下部の領域と、を除く領域は、架橋反応の反応率の観点からは好ましい領域である。なお、表３には、後述する内部気泡８０に関する情報は記載されていない。

【0182】

表４は、表２に、内部気泡８０に関する情報を追記したものである。

【0183】

【表4】

【0184】

表４のうち、右上部において太い罫線で囲まれた領域は、反応率が１１５％より大きい値となっている。この領域においては、架橋反応の戻り反応が過剰に進行しているので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、１１５％と異なる値を閾値としても良い。

【0185】

表４のうち、二重線の罫線で囲まれた部分においては、対象ワークモデルＷＭのポリマー部１２において内部気泡８０が発生したことを示す。

【0186】

表４の左下部において二重線の罫線で囲まれた部分のうち破線よりも下方の領域は、反応率が９５％よりも小さな値となっている。この領域においては、架橋反応が十分に進行していないので、好ましくない。ただし、反応率の範囲は任意に設定可能であり、９５％と異なる値を閾値としても良い。

【0187】

表４において、二重線の罫線で囲まれた領域のうち、破線よりも上方の領域は、架橋反応の反応率は９５％～１１５％の範囲内であったが、対象ワークモデルＷＭのポリマー部１２に内部気泡８０が発生したものであることを示す。この領域については、ポリマー部１２に内部気泡８０が発生したので、製品として好ましくない。

【0188】

本形態では、表２に記載された予測結果に、内部気泡８０に関する情報を追加することにより表４を作成したが、これに限られず、内部気泡８０に関する情報は、任意の温度における反応率の予測結果に追加することができる。

【0189】

（２）弾性率を表示する形態
表５は、対象ワークモデルＷＭの弾性率について、金型温度と、反応時間との関係についてまとめた表である。

【0190】

表５に記載された反応時間は、金型１０の内部に対象ワークモデルＷＭが配置された状態における反応時間を表す。また、表５に記載された弾性率の値は、金型１０の内部に対象ワークモデルＷＭが配置された状態における反応時間と、対象ワークモデルＷＭが金型１０から脱型された後の時間と、の双方の温度履歴に基づいて算出されたものである。対象ワークモデルＷＭが金型１０から脱型された後も架橋反応が進行することから、上記のような考慮が必要となる。なお、以下の説明において、図２８～図３１に記載された反応時間および弾性率についても同様である。

【0191】

【表5】

【0192】

弾性率は、対象ワークモデルＷＭの軸線Ａに直角な方向について静的な力を加えたときの弾性率である。表中の数値は、架橋反応の反応率が１００％であるときの弾性率に対する比の値である。ただし、対象ワークモデルＷＭに加えられる力の方向は限定されず、対象ワークモデルＷＭの軸線Ａに対して平行な方向について静的な力を加えた場合、対象ワークモデルＷＭの軸線Ａを軸にして回転する方向について静的な力を加えた場合、対象ワークモデルＷＭの軸線Ａに対して傾ける方向について静的な力を加えた場合等、任意の方向の力を加えた場合の弾性率を採用しうる。また、対象ワークモデルＷＭに加えられる力は静的な力に限定されず、例えば、対象ワークモデルＷＭに振動が加えられた場合の弾性率であってもよい。

【0193】

表５は、対象ワークモデルWMの各要素に、加硫度に応じた弾性率を割り当てて構造解析を行った結果得られる弾性率を、金型温度と、反応時間との関係についてまとめた表である。表５に記載された、太い罫線、二重線の罫線および破線の罫線の意味は、表４に記載された内容と同じなので、重複する説明を省略する。なお、表５において、「－」が記入されたセルについては、弾性率のシミュレーションにおいて計算エラーが発生した箇所である。

【0194】

さらに、表５には、例えば内部気泡８０に関する情報等、任意の情報を追加することができる。

【0195】

本形態においては、表に記載された数値を区分するために罫線を用いたが、これに限られず、表に記載された数値の色を変えたり、フォントを変えたり、書体を斜字体や太字に変えたり、数値が記載されたセルの背景色を変えたり、セルの背景のパターンを変えたりする等、任意の手法により、表中の数値を区分することができる。

【0196】

７－２．ヒートマップ図９０による表示方法
図２７を参照して、架橋反応の反応率を、ヒートマップ図９０を用いて表示する形態について説明する。図２７には、対象ワークモデルＷＭについて、架橋反応の反応率の分布を濃淡で表示したヒートマップ図９０を示す。

【0197】

ヒートマップ図９０においては、対象ワークモデルＷＭのうち、外側接合部材１３および内側接合部材１４は白抜きで表示されている。また、対象ワークモデルＷＭのうち、ポリマー部１２は架橋反応の反応率の大小にしたがって領域９１～領域９４に分割して表示され、分割された各領域９１～領域９４が、架橋反応率の大小に従って濃淡表示されている。

【0198】

領域９１は、架橋反応の反応率が最も高く、最も濃い色調で表示されている。領域９２は、架橋反応の反応率が二番目に高く、二番目に濃い色調で表示されている。領域９３は、架橋反応の反応率が三番目に高く、三番目に濃い色調で表示されている。領域９４は、架橋反応の反応率が最も低く、最も薄い色調で表示されている。

【0199】

表示部７はヒートマップ図９０を作業者に表示する。作業者は、対象ワークモデルＷＭのポリマー部１２における架橋反応の反応率について、濃淡表示された領域９１～領域９４を視認することにより、ポリマー部１２の架橋反応の反応率の分布について直感的に理解することができる。

【0200】

ヒートマップ図９０は、対象ワークモデルＷＭに関するすべての反応条件（例えば、温度、時間等）について、作製し、出力することができる。また、ヒートマップ図９０は、所定の反応条件についてのみ作製し、出力する態様としても良い。

【0201】

表示部７によるヒートマップ図９０の表示方法は特に限定されず、例えば、表示部７が上記した表５を表示し、作業者が、表示された表５の各セルを選択することにより、当該セルに対応する反応条件におけるヒートマップ図９０が表示される形態としてもよい。また、作業者が架橋反応の反応条件を入力することにより、入力された反応条件に対応するヒートマップ図９０が表示される態様としても良い。

【0202】

本形態においては、ヒートマップ図９０は、ポリマー部１２が領域９１～領域９４に分割して表示される構成としたが、これに限られず、ポリマー部１２が２個～３個、または５個以上の領域に分割されて表示される構成としても良い。また、ヒートマップ図９０において、反応率の高い領域を薄い色調で表示し、反応率の低い領域を濃い色調で表示しても良い。

【0203】

また、ヒートマップ図９０において、反応率の高い領域を赤色などの暖色系の色で表示し、反応率の低い領域を青色などの寒色系の色で表示してもよいし、反応率の高い領域を青色などの寒色系の色で表示し、反応率の低い領域を赤色などの暖色系の色で表示してもよい。

【0204】

７－３．グラフによる表示方法
（１）弾性率を表示する形態
図２８～図３０を参照して、対象ワークモデルＷＭの弾性率について、グラフを用いて表示する形態について説明する。図２８～図３０は、対象ワークモデルＷＭの弾性率の、架橋反応の反応時間に対するグラフである。図２８は、最も低い金型１０の温度で、対象ワークモデルＷＭのポリマー部１２の架橋反応を行ったときのグラフであり、図２９は、中程度の金型１０の温度で、対象ワークモデルＷＭのポリマー部１２の架橋反応を行ったときのグラフであり、図３０は、最も高い金型１０の温度で、対象ワークモデルＷＭのポリマー部１２の架橋反応を行ったときのグラフである。図２８～図３０において、実線は弾性率の実測値を示し、破線は予測値を示す。

【0205】

図２８～図３０に示すように、弾性率は、架橋反応が開始してから上昇して極大値に達し、極大値に達した後は減少する。

【0206】

金型１０の温度が最も低い図２８においては、弾性率の、実測値と予測値とは精度よく一致した。

【0207】

金型１０の温度が中程度である図２９においては、反応開始後２分までは、弾性率の実測値と予測値とが一致しなかったが、２分以後は、実測値と予測値とは精度よく一致した。反応開始後２分までは、弾性率は、実測値よりも小さく予測されていた。しかしながら、反応開始後２分以内の領域については、製品として利用されることは考えられないので、製品への影響は少ない。

【0208】

金型１０温度が最も高い図３０においては、反応開始後２分までは、弾性率の実測値と予測値とが一致しなかったが、２分以後は、実測値と予測値とは精度よく一致した。反応開始後２分までは、弾性率は、実測値よりも小さく予測されていた。実測値と予測値との差は、金型１０温度が中程度である場合よりも大きかった。しかしながら、上記したように、反応開始後２分以内の領域については、製品として利用されることは考えられないので、製品への影響は少ない。

【0209】

（２）弾性率と内部気泡８０とを併せて表示する形態
図３１を参照して、対象ワークモデルＷＭのポリマー部１２の弾性率と、内部気泡８０の発生についての予測結果を、グラフを用いて表示する形態について説明する。金型１０の温度は、上記した図２９と同様に、中程度の温度に設定されている。

【0210】

図３１には、対象ワークモデルＷＭの軸線Ａに直角な方向に力が加えられたときの弾性率の、反応時間に対するグラフを示す。四角形状のシンボルは、対象ワークモデルＷＭのポリマー部１２の温度が最高値である要素における弾性率を示す。丸形状のシンボルは、対象ワークモデルＷＭのポリマー部１２の温度が中央値である要素における弾性率を示す。三角形状のシンボルは、対象ワークモデルＷＭのポリマー部１２の温度が最低値である要素における弾性率を示す。また、図３１において、白抜きのシンボルは、内部気泡８０が発生すると予測された場合を示し、黒塗りのシンボルは、内部気泡８０が発生しないと予測された場合を示す。

【0211】

なお、ポリマー部１２の温度が最高値である要素と、中央値である要素についてのグラフは、反応進行期を過ぎて、反応戻り期におけるグラフとなっている。

【0212】

ポリマー部１２の温度が最高値である要素においては、架橋反応の進行が最も早い。このため、反応時間が３分を超えて、３分３０秒経過したときに、内部気泡８０が発生しないと予測された。ポリマー部１２の温度が最も高い要素においては、架橋反応の戻り反応の進行も早くなるので、反応時間６分経過時の弾性率は、最も低くなっている。

【0213】

ポリマー部１２の温度が中央値である要素においては、架橋反応の進行も中程度である。このため、反応時間が３分３０秒を超えて４分経過したときに、内部気泡８０が発生しないと予測された。反応時間６分経過時の弾性率は、ポリマー部１２の温度が最高値である要素よりも大きい。

【0214】

ポリマー部１２の温度が最低値である要素においては、架橋反応の進行が最も遅い。このため、反応時間が４分より前の状態においては、反応進行期となっており、弾性率は反応時間の経過とともに増加すると予測された。反応時間が４分経過した後は、架橋反応の戻り反応が進行し、弾性率が漸減すると予測された。ポリマー部１２の温度が最低値である要素においては、反応時間が５分を超えて、５分３０秒経過したときに、内部気泡８０が発生しないと予測された。この時点で、対象ワークモデルＷＭのポリマー部１２の全領域において内部気泡８０が消失したと予測された。反応時間６分経過時の弾性率は、ポリマー部１２の全領域の中で最も大きいと予測された。

【0215】

上記したように、反応時間が５分３０秒経過したときに、対象ワークモデルＷＭのポリマー部１２の全領域において内部気泡８０が消失したと予測されたが、危険率を考慮して、反応時間６分経過時を、この対象ワークモデルＷＭの反応時間として採用した。反応時間６分経過時において、ポリマー部１２の弾性率のばらつきは、１２％であると予測された。

【0216】

このように、本形態によれば、対象ワークモデルＷＭについて、内部気泡８０が発生しなくなる反応時間を予測可能であるとともに、当該反応時間におけるポリマー部１２の弾性率のばらつきも予測することができる。

【0217】

（３）内部気泡８０を視覚的に表示する形態
図３２を参照して、対象ワークモデルＷＭのポリマー部１２における、内部気泡８０の発生の予測結果を視覚的に表示する形態について説明する。

【0218】

図３２は、対象ワークモデルＷＭの断面図を利用して、ポリマー部１２における内部気泡８０の発生可能性の高さを濃淡で表示した図である。対象ワークモデルＷＭの外側接合部材１３および内側接合部材１４は、図３２において最も濃いパターンで示されている。

【0219】

対象ワークモデルＷＭのポリマー部１２は、図３２において最も薄いパターンで示されている。ポリマー部１２の内部には、楕円形状の領域Ｐと、領域Pの内部に位置する楕円形状の領域Ｑが表示されている。領域Ｐは、ポリマー部１２を表す薄いパターンよりもやや濃いパターンで表示されている。領域Ｑは、領域Ｐよりもやや濃いパターンで表示されている。

【0220】

内部気泡推定処理Ｓ４の結果、例えば、領域Ｐが比較的に内部気泡８０の発生しやすい領域と判定され、領域Ｑが領域Ｐよりも内部気泡８０の発生しやすい領域と判定された場合を例にして説明する。領域Ｐおよび領域Ｑは、実際に内部気泡８０の発生しやすい領域に対応している（図２３参照）。

【0221】

内部気泡推定処理Ｓ４が終了すると、表示部７は、図３２に示す対象ワークモデルＷＭを表示する。作業者は、表示部７に表示された図３２を見ることにより、比較的に濃いパターンで示される外側接合部材１３および内側接合部材１４と、比較的に薄いパターンで示されるポリマー部１２と、を視認する。作業者は、濃いパターンで示される外側接合部材１３および内側接合部材１４が、ポリマー部１２と異なる部分であることを直感的に理解することができる。これにより、作業者は、内部気泡８０を探すためには、比較的に薄いパターンで示されるポリマー部１２を確認すればよいことを、直感的に理解することができる。ただし、外側接合部材１３および内側接合部材１４を白抜きで表示しても良い。

【0222】

次に、作業者は、薄いパターンで示されるポリマー部１２の中に、比較的に濃いパターンで示される領域Ｐおよび領域Ｑが表示されていることを視認する。これにより、作業者は、ポリマー部１２のうち、領域Ｐおよび領域Ｑにおいて、内部気泡８０が発生しやすいことを直感的に理解することができる。さらに、領域Ｑは領域Ｐよりも濃いパターンで表示されているので、領域Ｐにおいて最も内部気泡８０が発生しやすいことを直感的に理解することができる。これにより、対象ワークモデルＷＭの製造条件を検討する際に、内部気泡８０が発生する可能性の高い領域において内部気泡８０を探せばよいので、対象ワークモデルＷＭの製造条件を効率的に検討することができる。ただし、内部気泡８０の発生する可能性の高い領域を比較的に薄いパターンで表示しても良い。

【0223】

本発明は上記各実施形態に限定されるものではなく、その要旨を逸脱しない範囲において種々の実施形態に適用することが可能である。

【0224】

本形態においては、内部気泡推定部６は、トルク算出部６１を備える構成としたが、これに限られず、対象ワークモデルＷＭのポリマー部１２の架橋反応の反応率（加硫度）を取得する反応率取得部６０と、脱型時における反応率に基づいて、対象ワークモデルＷＭのポリマー部１２における気泡の発生を推定する推定部６２と、を備える構成としても良い。この場合、加硫反応の反応解析を行う前に、基準反応曲線ＲＣから、閾値となるトルクにおける加硫度を求め、この加硫度をブローポイント加硫度とすることができる。また、ブローポイント加硫度は、例えば、ブローポイントアナライザや、ブローポイントテスタ等の公知の測定機器を用いて直接計測しても良い。その後、加硫反応の反応解析を行い、反応率取得部６０の演算結果（加硫度）と、上記のブローポイント加硫度と、を比較して、内部気泡の発生を推定することができる。例えば、加硫度が、ブローポイント加硫度を超えていない場合に、対象ワークモデルＷＭのポリマー部１２において内部気泡８０が発生したと推定する。

【符号の説明】

【0225】

１：架橋反応シミュレーション装置、２：記憶部、３：伝熱解析部、４：架橋反応解析部、５：構造解析部、６：内部気泡推定部、７：表示部、１０：金型、１１：熱板、１２：ポリマー部、１３：外側接合部材、１４：内側接合部材、１５：隙間、１５Ａ：外側隙間、１５Ｂ：内側隙間、３０：条件入力部、３１：ポリマー熱拡散率決定部、３２：解析部、４０：温度取得部、４１：反応率算出処理部、５０：温度取得部、５１：反応率取得部、５２：弾性率割当部、５３：特性取得部、６０：反応率取得部、６１：トルク算出部、６２：推定部、７０：加硫試験機、８０：内部気泡、Ａ：軸線、ＡＣ：空気熱伝達係数、ＢＤ：内部気泡推定用データ、Ｃ：比熱、ＣＣ：架橋反応曲線、ＤＭ１：第一関係データマップ、ＤＭ２：第二関係データマップ、Ｅ：活性化エネルギー、ＥＭ：等価反応量算出モデル、Ｅ_Ｎ：活性化エネルギー、Ｆ１：第一関数、Ｆ２：第二関数、ＣＨ：接触熱伝達率、ｋ：反応速度、ＭＭ：成形型モデル、ＯＣ：外気条件、ＲＣ：基準反応曲線、ＲＤ：架橋反応解析用データ、Ｓ１：伝熱解析処理、Ｓ２：架橋反応解析処理、Ｓ３：構造解析処理、Ｓ４：内部気泡推定処理、Ｓ５：表示処理、ＳＣ：係数、ｔ：反応時間、Ｔ：反応温度、ｔ_０：基準時間、Ｔ_０：基準温度、ＴＣ：温度条件、ＴＤ：伝熱解析用データ、ＴＳ：熱拡散率特性、Ｕ：等価反応量、ＷＭ：対象ワークモデル、Ｚ：反応量、α：熱拡散率、Δｔ：微小時間、ΔＵｉ：等価反応量増加量、λ：熱伝導率

【要約】

【課題】架橋反応の反応率を精度よく推定可能な架橋反応シミュレーション装置を提供する。
【解決手段】記憶部２と、対象ワークモデルＷＭのポリマー部１２の架橋反応時の伝熱解析を行う伝熱解析部３と、ポリマー部１２の架橋反応の反応率の解析を行う架橋反応解析部４と、を備える架橋反応シミュレーション装置１であって、記憶部２は、ポリマー部１２を有する対象ワークモデルＷＭと、等価反応量算出モデルＥＭと、ポリマー部１２の架橋反応の進行度に応じて設定された傾き係数ＳＣと、を記憶し、架橋反応解析部４は、伝熱解析の結果として、架橋反応における対象ワークモデルＷＭのポリマー部１２の各要素について時刻毎の温度を取得する温度取得部４０と、時刻毎のポリマー部１２の等価反応量を算出し、算出した等価反応量に基づいてポリマー部１２の架橋反応の反応率を算出する反応率算出処理部４１と、を備える。
【選択図】図１

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】

【図22】

【図23】

【図24】

【図25】

【図26】

【図27】

【図28】

【図29】

【図30】

【図31】

【図32】