(57)【要約】【課題】単量体の逐次付加反応をシミュレーションする。【解決手段】所定の単量体の逐次付加反応をシミュレーションする方法は、所定の付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定するステップと、特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して所定の反応速度定数を割り当てるステップと、各反応に割り当てられた反応速度定数と所定の反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、特定した物質の範囲で、各物質の生成量を算出するステップと、を備える。【選択図】なし【特許請求の範囲】
【請求項１】
所定の単量体の逐次付加反応をシミュレーションするシミュレーション装置であって、
所定の付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定する物質特定手段と、
上記物質特定手段で特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して所定の反応速度定数を割り当てる反応速度定数割当手段と、
上記反応速度定数割当手段によって各反応に割り当てられた反応速度定数と所定の反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、上記物質特定手段で特定した物質の範囲で、各物質の生成量を算出する生成量算出手段と、
を備えたシミュレーション装置。
【請求項２】
単量体の種類を入力する単量体種入力手段をさらに備えた請求項１に記載されたシミュレーション装置。
【請求項３】
付加単量体数を入力する付加単量体数入力手段をさらに備えた請求項１又は２に記載されたシミュレーション装置。
【請求項４】
反応パターンに対応した反応速度定数を入力する反応速度定数入力手段をさらに備えた請求項１乃至３のいずれかに記載されたシミュレーション装置。
【請求項５】
実験データを入力する実験データ入力手段と、
上記実験データ入力手段から入力された実験データと上記生成量算出手段で算出する生成物の生成量から導出したその実験データに対応する計算データとが最小二乗法により最も近似するように上記反応速度定数入力手段から入力された反応速度定数を補正する反応速度定数補正手段と、
をさらに備えた請求項４に記載されたシミュレーション装置。
【請求項６】
反応初期物質量を入力する反応初期物質量入力手段をさらに備えた請求項１乃至５のいずれかに記載されたシミュレーション装置。
【請求項７】
所定の単量体の逐次付加反応をシミュレーションする方法であって、
所定の付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定するステップと、
上記特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して所定の反応速度定数を割り当てるステップと、
上記各反応に割り当てられた反応速度定数と所定の反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、上記特定した物質の範囲で、各物質の生成量を算出するステップと、
を備えたシミュレーション方法。
【請求項８】
所定の単量体の逐次付加反応をコンピュータを用いてシミュレーションするコンピュータプログラムであって、
コンピュータに、
所定の付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定する手順と、
上記特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して所定の反応速度定数を割り当てる手順と、
上記各反応に割り当てられた反応速度定数と所定の反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、上記特定した物質の範囲で、各物質の生成量を算出する手順と、
を実行させるコンピュータプログラム。
【請求項９】
請求項８に記載されたコンピュータプログラムを記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、所定の逐次付加反応をシミュレーションするシミュレーション装置、シミュレーション方法、それを実行するためのコンピュータプログラム及びそれを記録した記録媒体に関する。
【背景技術】
【０００２】
反応速度論に基づいた反応シミュレーションは一般的に行われており、特許文献１及び２に記載されている。
【０００３】
また、特許文献３には、複数の反応が同時に進行する化学反応系において、所定の反応時間における反応物中の各成分濃度を分析手段によって測定し、そして、コンピュータによって、所定の反応時間における各成分濃度を、ルンゲ・クッタ・ギル法により複数のパラメータを含む微分方程式(反応速度式)を数値積分することによって推算し、さらに、コンピュータによって、各成分濃度の実測値と推算値の偏差の二乗和を最小にするようなパラメータをマルカート法により演算することにより、反応速度式をシミュレーションにより確定することが開示されている。
【特許文献１】特開２００６−２７６９８６号公報
【特許文献２】特開２００３−２１３００３号公報
【特許文献３】特開平０６−０５５０６２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
例えば、グリシドールを逐次付加反応させた場合、グリシドールの付加位置によっては分岐鎖を形成するため、付加単量体数が増加すると、得られる生成物の種類が指数関数的に増大する。従って、このような場合、特定の生成物がどの程度の割合で生成しているのかといった詳細な情報を得ることは非常に困難である。
【０００５】
本出願は、所定の単量体の逐次付加反応をシミュレーションするシミュレーション装置、シミュレーション方法、それを実行するためのコンピュータプログラム及びそれを記録した記録媒体を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明のシミュレーション装置は、所定の単量体の逐次付加反応をシミュレーションするものであって、
所定の付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定する物質特定手段と、
上記物質特定手段で特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して所定の反応速度定数を割り当てる反応速度定数割当手段と、
上記反応速度定数割当手段によって各反応に割り当てられた反応速度定数と所定の反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、上記物質特定手段で特定した物質の範囲で、各物質の生成量を算出する生成量算出手段と、
を備える。
【０００７】
本発明のシミュレーション方法は、所定の単量体の逐次付加反応をシミュレーションするものであって、
所定の付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定するステップと、
上記特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して所定の反応速度定数を割り当てるステップと、
上記各反応に割り当てられた反応速度定数と所定の反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、上記特定した物質の範囲で、各物質の生成量を算出するステップと、
を備える。
【０００８】
本発明のコンピュータプログラムは、所定の単量体の逐次付加反応をコンピュータを用いてシミュレーションするものであって、
コンピュータに、
所定の付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定する手順と、
上記特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して所定の反応速度定数を割り当てる手順と、
上記各反応に割り当てられた反応速度定数と所定の反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、上記特定した物質の範囲で、各物質の生成量を算出する手順と、
を実行させるものである。
【０００９】
本発明の記録媒体は、そのコンピュータプログラムを記録したものである。
【発明の効果】
【００１０】
本発明によれば、所定の単量体の逐次付加反応をシミュレーションすることができる。
【発明を実施するための最良の形態】
【００１１】
以下、実施形態について詳細に説明する。
【００１２】
図１は、本実施形態に係るシミュレーション装置１０を示す。
【００１３】
このシミュレーション装置１０は、ＣＰＵ１１、メモリ１２、ハードディスク１３、キーボードやマウスなどの情報入力部１４、及びディスプレイなどの情報表示部１５を備え、ハードディスク１３に逐次付加反応シミュレーション用のコンピュータプログラムがインストールされたコンピュータで構成されている。
【００１４】
そして、この逐次付加反応シミュレーション用のコンピュータプログラムは、大略として、コンピュータに、
情報入力部１４からの情報入力を情報表示部１５への表示により促す手順と、
付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質を特定する手順と、
上記特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応して反応速度定数を割り当てる手順と、
上記各反応に割り当てられた反応速度定数と反応初期物質量とに基づき、反応速度式を所定期間で時間積分することにより、生成物の生成量を算出する手順と、
上記算出した生成物の生成量に基づいて、情報表示部１５に所定の結果情報を表示させる手順と、
を実行させるものである。なお、このコンピュータプログラムには、複数種の単量体の逐次付加反応のそれぞれについて、単量体の採りうる付加構造等の反応情報がデータベースとして含まれている。
【００１５】
このコンピュータプログラムは、ＣＤ−ＲＯＭ等のコンピュータ読み取り可能な記録媒体での商品化が可能なものである。
【００１６】
次に、このシミュレーション装置１０を用いたシミュレーション方法について、図２及び３のフローチャートに従って、グリシドール逐次付加反応を事例としつつ説明する。
【００１７】
なお、グリシドール逐次付加反応は、グリシドールの付加位置によって分岐鎖を形成するため、付加単量体数が増加すると、生成する生成物の種類が指数関数的に増大する。また、グリシドール逐次付加反応には、グリシドールの１位の炭素原子が付加するMarkownikov付加（以下、「正付加」という。）に加えて、２位の炭素原子が付加する逆Markownikov付加（以下、「逆付加」という。）が５％程度含まれるため、表１に示すように、前者の場合だけでも相当数の種類の生成物が生成するが、後者の場合をも含めると、その数はさらに多くなる。
【００１８】
【表１】
2009301379.tif 000002
【００１９】
スタート後のステップＳ１では、情報表示部１５に、情報入力部１４から出発物質を入力することを促す表示をしてステップＳ２に進み、適切な入力があった場合にはステップＳ３に進み、適切な入力がない場合にはステップＳ１に戻る。
【００２０】
本事例では、情報表示部１５に「出発物質？」と表示し、これに対して情報入力部１４から「ＣＯ」と入力することとする。なお、「ＣＯ」は１価アルコールをＳｍｉｌｅｓ表記したものである。
【００２１】
ステップＳ３では、情報表示部１５に、情報入力部１４から単量体の種類を入力することを促す表示をしてステップＳ４に進み、適切な入力があった場合にはステップＳ５に進んで情報表示部１５に単量体の採りうる付加構造を表示した後にステップＳ６に進み、適切な入力がない場合にはステップＳ３に戻る。従って、情報入力部１４が単量体種入力手段を構成する。
【００２２】
本事例では、情報表示部１５に「単量体の種類１，２？」と表示し、これに対して情報入力部１４からグリシドールを特定する「２」を選択して入力することとする。このとき、コンピュータプログラムに含まれたデータベースに基づいて、情報表示部１５に、グリシドールの採りうる付加構造として、「付加構造：ＣＣ（Ｏ）ＣＯ」及び「付加構造：Ｃ（ＣＯ）ＣＯ」と表示する。前者は１位の炭素原子が付加する正付加をＳｍｉｌｅｓ表記したものであり、後者は２位の炭素原子が付加する逆付加をＳｍｉｌｅｓ表記したものである。
【００２３】
ステップＳ６では、情報表示部１５に、情報入力部１４から付加単量体数を入力することを促す表示をしてステップＳ７に進み、適切な入力があった場合にはステップＳ８に進み、適切な入力がない場合にはステップＳ６に戻る。従って、情報入力部１４が付加単量体数入力手段を構成する。
【００２４】
本事例では、情報表示部１５に「考慮する生成物の最大付加単量体数？」と表示し、これに対して情報入力部１４から「７」と入力することとする。
【００２５】
ステップＳ８では、ステップＳ６で入力された付加単量体数以下で構成可能な相互に分子構造の異なる全ての物質をＳｍｉｌｅｓ表記で特定し、ステップＳ９に進む。従って、ＣＰＵ１１によるこのステップＳ８の実行が物質特定手段を構成する。
【００２６】
本事例では、付加単量体数を「７」とするので、付加単量体数７以下で構成可能な相互に分子構造の異なる全ての物質をＳｍｉｌｅｓ表記で特定する。この場合、表１より生成物が１８９９１種類となるので、ここでは、付加単量体数が３以下で構成可能な相互に分子構造の異なる全ての物質をＳｍｉｌｅｓ表記で特定したものを表２に示す。なお、生成物には、付加単量体数が０、つまり、未付加の１価アルコールもその１つに含まれる。
【００２７】
【表２】
2009301379.tif 000003
【００２８】
ステップＳ９では、情報表示部１５に、情報入力部１４から単量体の当量を入力することを促す表示をしてステップＳ１０に進み、適切な入力があった場合にはステップＳ１１に進んで出発物質添加量及び単量体添加量を表示した後にステップＳ１２に進み、適切な入力がない場合にはステップＳ９に戻る。従って、情報入力部１４が反応初期物質量入力手段を構成する。
【００２９】
本事例では、情報表示部１５に「付加する単量体の当量？」と表示し、これに対して情報入力部１４から「４」と入力することとする。このとき、情報表示部１５に「出発物質添加量：１００．０」及び「単量体  添加量：４００．０」と表示する。これは、反応初期条件として、出発物質である１価アルコールの初期添加量として１００モル、及びグリシドールの初期添加量として４００モルがそれぞれ入力されたことを意味する。
【００３０】
ステップＳ１２では、情報表示部１５に、情報入力部１４から逆付加の確率を入力することを促す表示をしてステップＳ１３に進み、適切な入力があった場合にはステップＳ１４に進んで各付加構造の反応性比を表示した後にステップＳ１５に進み、適切な入力がない場合にはステップＳ１２に戻る。
【００３１】
本事例では、情報表示部１５に「逆付加の確率？」と表示し、これに対して情報入力部１４から「０．０５」と入力することとする。このとき、情報表示部１５に「反応性比：０．９５０  ＣＣ（Ｏ）ＣＯ」及び「反応性比：０．０５０  Ｃ（ＣＯ）ＣＯ」と表示する。これは、前者の正付加の起こる確率として９５％、及び後者の逆付加の起こる確率として５％がそれぞれ入力されたことを意味する。
【００３２】
ステップＳ１５では、情報表示部１５に、情報入力部１４から反応速度定数を入力することを促す表示をしてステップＳ１６に進み、適切な入力があった場合にはステップＳ１７に進んで反応速度定数を表示した後にステップＳ１８に進み、適切な入力がない場合にはステップＳ１５に戻る。従って、情報入力部１４が反応速度定数入力手段を構成する。
【００３３】
本事例では、情報表示部１５に「反応速度定数ｋ１？」、「反応速度定数ｋ２？」、及び「反応速度定数ｋ３？」と表示し、これに対して情報入力部１４から順に「０．１」「０．１」、及び「０．０５」と入力することとする。このとき、情報表示部１５に「反応速度定数ｋ１：０．１００００」、「反応速度定数ｋ２：０．１００００」、及び「反応速度定数ｋ３：０．０５０００」と表示する。これは、反応速度定数ｋ１として０．１モル／（単位時間）、反応速度定数ｋ２として０．１モル／（単位時間）、及び反応速度定数ｋ３として０．０５モル／（単位時間）、がそれぞれ入力されたことを意味する。ここで、グリシドール逐次付加反応では、下記に示すように、出発物質である１価アルコールへの付加（［１］）、末端又は末端以外への１級付加（［２］及び［３］）、並びに末端又は末端以外への２級付加（［４］及び［５］）に反応パターンが分かれる。そして、反応速度定数ｋ１は出発物質である１価アルコールへの付加反応のものであり、反応速度定数ｋ２は末端又は末端以外への１級付加反応のものであり、反応速度定数ｋ３は末端又は末端以外への２級付加反応のものである。
【００３４】
【化１】
2009301379.tif 000004
【００３５】
ステップＳ１８では、情報表示部１５に、情報入力部１４から時間ステップを入力することを促す表示をしてステップＳ１９に進み、適切な入力があった場合にはステップＳ２０に進んで時間ステップを表示した後にステップＳ２１に進み、適切な入力がない場合にはステップＳ１８に戻る。
【００３６】
本事例では、情報表示部１５に「時間ステップ？」と表示し、これに対して情報入力部１４から「０．０１」と入力することとする。このとき、情報表示部１５に「時間ステップ：０．０１０」と表示する。これは、時間ステップとして０．０１（単位時間）が入力されたことを意味する。ここで、時間ステップは、時間積分における積分する時間間隔である。
【００３７】
ステップＳ２１では、情報表示部１５に、情報入力部１４から初期物質残留率、分岐物質含有率、及び１，２付加物質含有率の実験データを入力することを促す表示をしてステップＳ２２に進み、適切な入力があった場合にはステップＳ２３に進んで実験データを表示した後にステップＳ２４に進み、適切な入力がない場合にはステップＳ２１に戻る。従って、情報入力部１４が実験データ入力手段を構成する。
【００３８】
本事例では、情報表示部１５に「アルコール残留率  実験データ？」、「分岐物質含有率  実験データ？」、及び「１，２付加物質含有率  実験データ？」と表示し、これに対して情報入力部１４から順に「７．１４」「１２．９４」、及び「１１．０４」と入力することとする。このとき、情報表示部１５に「アルコール残留率  実験データ：７．１４０」、「分岐物質含有率  実験データ：１２．９４０」、及び「１，２付加物質含有率  実験データ：１１．０４０」と表示する。これは、アルコール残留率の実験データとして７．１４０％、分岐物質含有率の実験データとして１２．９４０％、及び１，２付加物質含有率の実験データとして１１．０４０％がそれぞれ入力されたことを意味する。
【００３９】
ステップＳ２４では、情報表示部１５に、シミュレーションを実行するかどうかの回答を促す表示をしてステップＳ２５に進み、所定の入力があった場合にはステップＳ２６に進み、その他の入力があった場合には終了する。
【００４０】
本事例では、情報表示部１５に「シミュレーションを実行しますか？（０／１）」と表示し、これに対して情報入力部１４から「１」と入力することとする。
【００４１】
ステップＳ２６では、ステップＳ８で特定した各物質について、その物質の全ての付加可能位置のそれぞれに単量体を付加する各反応に対し、その反応パターンに対応してステップＳ１５で入力された反応速度定数のいずれかを割り当て、ステップＳ２７に進む。従って、ＣＰＵ１１によるこのステップＳ２６の実行が反応速度定数割当手段を構成する。
【００４２】
本事例では、１８９９１種類の各物質について、その物質の全ての付加可能位置のそれぞれにグリシドールを付加する各反応に対し、その反応パターンに対応して反応速度定数ｋ１〜ｋ３のいずれかを割り当てる。具体的には、出発物質である１価アルコールへの付加反応には反応速度定数ｋ１、末端又は末端以外への１級付加反応には反応速度定数ｋ２、末端又は末端以外への２級付加反応には反応速度定数ｋ３をそれぞれ割り当てる。
【００４３】
ステップＳ２７では、各反応に割り当てた反応速度定数とステップＳ９で入力された単量体添加量（すなわち、反応初期物質量）とに基づき、反応速度式を反応開始から反応終了までステップＳ１８で入力された時間ステップで時間積分することにより、ステップＳ８で特定した物質の範囲で、各物質の生成量を算出し、ステップＳ２８に進む。従って、ＣＰＵ１１によるこのステップＳ２７の実行が生産量算出手段を構成する。
【００４４】
本事例では、出発物質である１価アルコールへの付加には下記（１）の反応速度式、末端又は末端以外への１級付加には下記（２）の反応速度式、及び末端又は末端以外への２級付加には下記（３）の反応速度式をそれぞれ適用し、そして、出発物質である１価アルコールの添加量を１００モル、単量体であるグリシドールの添加量を４００モルとした条件の下、表３に示すように、グリシドールが消費され尽くされるまで時間ステップ０．０５（単位時間）で時間積分を実行し、１８９９１種類の物質の範囲で、各物質の生成量（モル）を算出する。なお、各物質の生成量の算出において、ある物質にグリシドールが付加したときの全生成量を０．９５倍したものを、その物質にグリシドールが正付加した物質の生成量として加算し、全生成量を０．０５倍したものを、その物質にグリシドールが逆付加した物質の生成量として加算する。
【００４５】
【数１】
2009301379.tif 000005
【００４６】
【表３】
2009301379.tif 000006
【００４７】
ステップＳ２８では、ステップＳ２７での計算結果に基づいて、各付加単量体数毎の生成量、分岐物質の生成量、及び１，２付加物質の生成量を計算すると共に、それらから初期物質残留率、分岐物質含有率、及び１，２付加物質含有率を導出し、ステップＳ２９に進んで情報表示部１５にそれらの計算データを表示した後にステップＳ３０に進む。
【００４８】
本事例では、各付加単量体数毎の生成量、分岐物質の生成量、及び１，２付加物質の生成量を計算すると共に、それらからアルコール残留率、分岐物質含有率、及び１，２付加物質含有率を導出する。このとき、情報表示部１５に表４に示す計算データを表示する。
【００４９】
【表４】
2009301379.tif 000007
【００５０】
ステップＳ３０では、情報表示部１５に、フィッティングを実行するかどうかの回答を促す表示をしてステップＳ３１に進み、所定の入力があった場合にはステップＳ３２に進み、その他の入力があった場合には終了する。
【００５１】
本事例では、情報表示部１５に「フィッティングを実行しますか？（０／１）」と表示し、これに対して情報入力部１４から「１」と入力することとする。
【００５２】
ステップＳ３２では、ステップＳ２１で入力された初期物質残留率、分岐物質含有率、及び１，２付加物質含有率の実験データと計算データとが最小二乗法により最も近似するように、ステップＳ１５で入力された反応速度定数を補正して、ステップＳ２７及びステップＳ２８の計算を繰り返し再実行し、ステップＳ３３に進んで情報表示部１５に最終計算データを表示した後にステップＳ３４に進む。従って、ＣＰＵ１１によるこのステップＳ３２の実行が生産反応速度定数補正手段を構成する。なお、反応速度定数の変量は、最小二乗法による計算結果に基づいてコンピュータプログラムのサブルーチンにおいて実行する。
【００５３】
本事例では、反応速度定数ｋ１〜ｋ３をそれぞれ変量し、各物質の生成量の算出、それに基づく各付加単量体数毎の生成量、分岐物質の生成量、及び１，２付加物質の生成量の計算、並びにそれらからのアルコール残留率、分岐物質含有率、及び１，２付加物質含有率の導出を、表５に示すように、実験データと計算データとが最小二乗法により最も近似するまで繰り返し実行する。このとき、情報表示部１５に表５の最下段の太枠に示す最終の計算データ結果を表示する。この場合、反応速度定数ｋ１を０．１０２１０モル／（単位時間）、反応速度定数ｋ２を０．１０５１２モル／（単位時間）、及び反応速度定数ｋ３を０．０４５１８モル／（単位時間）とすることで、実験データと計算データとが最も近似するということである。
【００５４】
【表５】
2009301379.tif 000008
【００５５】
ステップＳ３４では、情報表示部１５に、生成物分布図を表示するかどうかの回答を促す表示をしてステップＳ３５に進み、所定の入力があった場合にはステップＳ３６に進んで情報表示部１５に生成物分布図を表示した後に終了し、その他の入力があった場合には終了する。
【００５６】
本事例では、情報表示部１５に「生成物分布図を表示しますか？（０／１）」と表示し、これに対して情報入力部１４から「１」と入力することとする。このとき、情報表示部１５に図４に示すフィッティング後の各付加単量体数毎の生成物分布図を表示する。
【００５７】
以上のシミュレーション装置１０を用いれば、逐次付加反応のシミュレーションを実行することにより、反応後の各物質の生成量が算出され、実験データでは得られない特性を知ることができ、例えば、これをプラントに適用する触媒の設計等に活用することができる。
【００５８】
なお、本実施形態において、シミュレーションを一通り実行した後、ステップＳ１〜Ｓ８の物質の特定までをスキップし、ステップＳ９から再シミュレーションを実行することができるように構成してもよい。
【００５９】
また、本実施形態において、ステップＳ２９の情報表示部１５への計算データの表示の後、情報表示部１５にフィッティング前の各付加単量体数毎の生成物分布図を表示することができるように構成してもよい。
【００６０】
また、本実施形態では、グリシドール逐次付加反応を事例としたが、特にこれに限定されるものではなく、例えば、アルキレンオキシド逐次付加反応等も同様にシミュレーションすることができる。
【００６１】
また、本実施形態では、単量体の種類、付加単量体数、単量体の当量、及び反応速度定数を入力する構成としたが、特にこれに限定されるものではなく、これらのうちの幾つかが予めコンピュータプログラムにおいてパラメータとして設定されていてもよい。
【産業上の利用可能性】
【００６２】
本発明は、逐次付加反応のシミュレーションに有用である。
【図面の簡単な説明】
【００６３】
【図１】実施形態に係るシミュレーション装置のブロック図である。
【図２】シミュレーション方法の前半部分のフローチャートである。
【図３】シミュレーション方法の後半部分のフローチャートである。
【図４】生成物分布図である。
【符号の説明】
【００６４】
１０  シミュレーション装置
１１  ＣＰＵ
１２  メモリ
１３  ハードディスク
１４  情報入力部
１５  情報表示部
【図１】
2009301379.tif 000009
【図２】
2009301379.tif 000010
【図３】
2009301379.tif 000011
【図４】
2009301379.tif 000012