応力−ひずみ曲線式を出力するためのプログラム及びその装置
【課題】 粘弾性だけでなく、ゴム弾性を併せて弾性材料の物性を評価できるモデルを提供する。
【解決手段】 弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力する。ひずみと弾性係数との間に相関性があることを見出した点が一つの特徴であり、これにより、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現することが可能となる。
【解決手段】 弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力する。ひずみと弾性係数との間に相関性があることを見出した点が一つの特徴であり、これにより、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現することが可能となる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現するために必要となる応力−ひずみ曲線式を出力するためのプログラム及びその装置に関する。
【背景技術】
【０００２】
従来より、弾性材料の物性を評価する指標として粘弾性があり、その粘弾性を表現するモデルとしては、弾性要素及び粘性要素を直列に配置した一般化Ｍａｘｗｅｌｌモデルが用いられる。一般化Ｍａｘｗｅｌｌモデルでは、実際の応力−時間曲線（応力緩和曲線）（図１（ａ）参照）、即ち、弾性材料の粘弾性を、以下の式（１）に示す応力及び時間の相関式（図１（ｂ）参照）で表現することができるため、以下の式（１）に示す応力及び時間の相関式は、応力−時間曲線式（応力緩和曲線式）として用いられる。
【数１】
ここで、σは応力、ｔは時間、Ｅは一般化Ｍａｘｗｅｌｌモデルにおける弾性要素の引張り弾性係数、τは緩和時間、を表す。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、近年は、弾性材料（例えば、粘着剤）の高品質化乃至高付加価値化に伴い、粘弾性だけでは弾性材料の物性を十分に評価できなくなってきている。具体的に言えば、非線形で且つ大変形を起こすような弾性材料であれば、その物性を評価するためには、粘弾性だけでなく、ゴム弾性を指標に加える必要がある。
【０００４】
ところが、一般化Ｍａｘｗｅｌｌモデルでは、ゴム弾性に関する相関式として、式（１）から以下の式（２）に示す応力及びひずみの相関式を導き出すことができるものの、以下の式（２）に示す応力及びひずみの相関式（図１（ｄ）参照）では、実際の応力−ひずみ曲線（図１（ｃ）参照）、即ち、弾性材料のゴム弾性を表現することができない。
【数２】
ここで、ηは一般化Ｍａｘｗｅｌｌモデルにおける粘性要素の粘性係数、εはひずみ、を表す。
【０００５】
そこで、本発明は、上記問題に鑑みてなされたもので、粘弾性モデルである一般化Ｍａｘｗｅｌｌモデルを発展させ、粘弾性−ゴム弾性両方の物性を評価できるモデル（以下、「一般化Ｍａｘｗｅｌｌ発展モデル」という）を作成し、その定式化を行った上で、それを活用できるツールを提供することを課題とする。
【課題を解決するための手段】
【０００６】
上記課題を解決すべく構成されたもので、ツールとしての本発明に係るプログラムは、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力するためのプログラムであって、コンピュータを、前記パラメータを入力する入力手段、該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段、及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段、として機能させることを特徴とする。
【０００７】
かかる構成によれば、物性に応じて異なるパラメータ（必要に応じてその他のパラメータも含む）を入力することで、弾性材料の物性に応じた応力−ひずみ曲線式が出力される。例えば、適宜、パラメータを変更して入力することで、所望の応力−ひずみ曲線式を追求し、また、それに合うような材料設計を行うといったアプリケーションが期待できる。
【０００８】
また、本発明に係るプログラムにおいては、一例として、前記ひずみ及び弾性係数の相関式は、弾性係数をひずみの指数関数で表したものであり、前記パラメータは、指数関数に乗じられる第一パラメータと、指数においてひずみに乗じられる第二パラメータとからなる。
【０００９】
その場合、さらなる一例として、応力及びひずみの相関式は、Cauchy応力テンソルを太字σ、Jaumann速度であることを右下付き(J)、弾性係数をＧ、ひずみ速度テンソルを太字Ｄ、前記第一パラメータをＡ、前記第二パラメータをＢ、ひずみのスカラー量をハット付きε、弾性要素の偏差応力テンソルを太字σ’、緩和時間をτ、対象を圧縮性線形弾性体とした場合の圧力をｐ、体積弾性係数をＫｖ、演算子traceをtr、単位テンソルを太字Ｉで表した場合、以下の式で規定される。
【数３】
【００１０】
また、上記課題を解決すべく構成されたもので、ツールとしての本発明に係る装置は、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力する装置であって、前記パラメータを入力する入力手段と、該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段と、及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段とを含むことを特徴とする。
【発明の効果】
【００１１】
以上の如く、本発明によれば、粘弾性モデルである一般化Ｍａｘｗｅｌｌモデルを発展させ、粘弾性−ゴム弾性両方の物性を評価できる一般化Ｍａｘｗｅｌｌ発展モデルを作成し、その定式化を行った上で、それを活用できるツールを提供することができる。また、そのツールを用いて、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現することができるようになる。
【図面の簡単な説明】
【００１２】
【図１】（ａ）は、実際の応力−時間曲線図、（ｂ）は、一般化Ｍａｘｗｅｌｌモデルに基づく計算解によって表現される応力−時間曲線図、（ｃ）は、実際の応力−ひずみ曲線図、（ｄ）は、一般化Ｍａｘｗｅｌｌモデルに基づく計算解によって表現される応力−ひずみ曲線図、を示す。
【図２】一般化Ｍａｘｗｅｌｌ発展モデルの概念図を示す。
【図３】１要素モデルの概念図を示す。
【図４】三次元モードでの一般化Ｍａｘｗｅｌｌ発展モデルの概念図を示す。
【図５】一般化Ｍａｘｗｅｌｌ発展モデルを検証するための１要素解析モデルの概念図を示す。
【図６】引張り速度５０ｍｍ／ｍｉｎにて実測した実際の応力−ひずみ曲線図と、一般化Ｍａｘｗｅｌｌ発展モデルに基づく計算解によって表現される応力−ひずみ曲線図とを示す。
【発明を実施するための形態】
【００１３】
以下、本発明の一実施形態について、図面を参酌しつつ説明する。
【００１４】
本実施形態は、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力するためのプログラムに言及する。そこで、以下においては、その応力及びひずみの相関式を導出するまでの流れを説明する。
【００１５】
上記式（１）で表される一般化Ｍａｘｗｅｌｌモデルではゴム弾性を表現できないことは上述の如くである。そこで、ゴム弾性も表現できるよう、一般化Ｍａｘｗｅｌｌモデルにおける弾性要素の引張り弾性係数Ｅ（並びに、一般化Ｍａｘｗｅｌｌモデルにおける粘性要素の粘性係数η）をひずみの関数とするモデル（一般化Ｍａｘｗｅｌｌ発展モデル）を想定した（図２参照）。尚、緩和時間τ（＝η(ε)／Ｅ(ε)）は、便宜上、一定とする。
【００１６】
その上で、まず始めに、粘弾性を表現する上記式（１）に対し、実測の応力緩和試験の値から回帰分析を行い、各ひずみにおける弾性係数Ｅと緩和時間τの組を求めた。より詳しくは、アクリル系ポリマーをシート状にして作製した粘着剤シートから断面積２ｍｍ2の円柱状の試験片を作成し、チャック間１０ｍｍ、速度１０００ｍｍ／ｍｉｎにて、所定のひずみになるまで伸張し、次に、ひずみを一定の状態に保持し、応力緩和試験を行い、上記式（１）を用いて回帰分析を行い、各ひずみにおける計五組の弾性係数Ｅ、緩和時間τの組を求めた。その結果を以下の表１に表す。
【表１】
尚、表１において、ひずみ（strain）は、対数ひずみである。公称ひずみで表せば、左から、１００％、２００％、３００％、４００％、５００％、６００％、となっている。
【００１７】
次に、弾性係数Ｅのひずみ依存性データをプロットし、カーブフィッティングを行うと、以下の表２に示す如く、弾性係数はひずみに対し、以下の式（３）に示す指数関数で近似することで、高い相関係数（表２中のＲ2値が１．０００に近いほど近似していることを示す）を持つことがわかった。即ち、ひずみ及び弾性係数の相関式は、弾性係数Ｅをひずみεの指数関数で表した以下の式（３）が妥当であることが確認できた。
【表２】
【数４】
ここで、Ａは指数関数に乗じられる第一パラメータ、Ｂは指数においてひずみεに乗じられる第二パラメータ、を表す。
従って、一般化Ｍａｘｗｅｌｌモデルでは、取得すべき物性値は、τｉ、Ｅｉである（上記式（１）参照）のに対し、一般化Ｍａｘｗｅｌｌ発展モデルでは、取得すべき物性値は、τｉ、Ａｉ、Ｂｉである。
【００１８】
次に、Ｍａｘｗｅｌｌ要素が一つしかない場合を考える。図３に示すように、応力σ、ひずみε、弾性要素の弾性係数Ｅ、粘性要素の粘性係数ηとすると、
【数５】
【数６】
【数７】
【数８】
となる。そして、上記式（６）を時間で微分すると、以下の式（８）が得られる。
【数９】
また、上記式（５）を時間で微分した式と上記式（７）、上記式（８）を用いて、εa、εbを消去して整理すると、以下の式（９）が得られる。
【数１０】
【００１９】
ここで、上記式（３）において、Ａ，Ｂを定数とし、
【数１１】
とおいて整理すると、以下の式（１０）が得られる。
【数１２】
これが図３に示す１要素モデルの構成式である。
【００２０】
次に、図２に示すような一般型の構成式を誘導する。一般型になっても、各要素では、上記式（１０）が成立するため、以下の式（１１）が成立する。
【数１３】
本モデルでは、各要素のひずみは一定であるため、ひずみに対して添え字ｉは付かない。そこで、一般型モデルの全体の応力σは各要素のσiの和であるから、以下の式（１２）が成立し、
【数１４】
一般化Ｍａｘｗｅｌｌ発展モデルの応力及びひずみの相関式が得られる。
【００２１】
ところで、上記式（１２）は、一次元モードの式であり、応力及びひずみの相関式を概念として表すためのものであるが、シミュレーションを行うためには、三次元モードでの定式化を行わなければならない。そこで、以下においては、三次元モードでの式を導出するまでの流れを説明する。
【００２２】
まず、一次元モードの場合と同様、三次元の観点で一般化Ｍａｘｗｅｌｌ発展モデルを図４のように設定すると、（弾性材料の物性に応じて異なる係数を含む、）ひずみ及び弾性係数の相関式は、以下の式（１３）のようになる。
【数１５】
ここで、Ｇはせん断弾性係数、ハット付きεはひずみのスカラー量、を表す。
尚、一次元モードの場合は、添え字ｉを右下付きで記載していたが、以下の三次元モードにおいては、添え字Ｋを右上付きで記載している。
【００２３】
ひずみのスカラー量（ハット付きε）を微小ひずみテンソル（太字ε）を用いて以下の式（１４）のように置くこととする。
【数１６】
今、Ｋ番目ユニットの弾性要素の変位応力テンソル（太字σ’）は、以下の式（１５）となる。
【数１７】
ここで、添え字がｓｐＫのεはＫ番目ユニットの弾性要素の微小ひずみテンソル、を表す。
【００２４】
上記式（１５）の両辺を物質時間微分すれば、以下の式（１６）となる。
【数１８】
ここで、太字Ｄはひずみ速度テンソル、添え字がｓｐＫの太字ＤはＫ番目ユニットの弾性要素のひずみ速度テンソル、を表す。
上記式（１６）の左辺は客観性のないテンソルであるため、客観応力速度を導入すると、以下の式（１７）を弾性要素の構成方程式として用いることができる。
【数１９】
ここで、右下付きの（＊）は任意の客観応力速度を表す。
【００２５】
一方、Ｋ番目ユニットの粘性要素における偏差成分の構成方程式は、以下の式（１８）で与えられる。
【数２０】
Ｋ番目ユニットの系全体の微小ひずみテンソルを太字ε、Ｋ番目ユニットの微小ひずみテンソルを添え字がＫの太字ε、で表すと、以下の式（１９）となり、
【数２１】
これを物質時間微分すると、以下の式（２０）となる。
【数２２】
【００２６】
上記式（１７）、（１８）、（２０）より、Ｋ番目ユニットの偏差成分の構成方程式は、以下の式（２１）のように導出される。
【数２３】
ここで、上記式（１３）は、以下の式（２２）となることから、
【数２４】
上記式（２１）は、以下の式（２３）となる。
【数２５】
モデルの性質から、モデル全体の応力はユニットＫ個の応力の総和であることから、偏差成分の応力は、以下の式（２４）となる。
【数２６】
【００２７】
次に、体積成分を考える。ここで、対象を圧縮性線形弾性体と仮定すると、圧力ｐは以下の式（２５）となる。
【数２７】
ここで、Ｋｖは体積弾性係数、trは演算子trace、を表し、体積弾性係数Ｋｖは、以下の式（２６）となる。
【数２８】
ここで、νはポアソン比を表す。
【００２８】
上記式（２５）を物質時間微分すると、以下の式（２７）となる。
【数２９】
また、以下の式（２８）であり、
【数３０】
さらに上記式（２２）を代入することにより、以下の式（２９）となる。
【数３１】
従って、体積成分の構成方程式は以下の式（３０）となる。
【数３２】
【００２９】
客観応力速度としてJaumann速度を用いることを考えると、Cauchy応力（太字σ）の物質時間導関数とそのJaumann速度の間には、以下の式（３１）の関係が成り立つ。
【数３３】
ここで、太字Ｗはスピンテンソル、右下付き(J)はJaumann速度であることを表す。
【００３０】
Cauchy応力（太字σ）を偏差成分と体積成分に分離すると、以下の式（３２）となる。
【数３４】
即ち、右辺の第一項が偏差成分であり、第二項が体積成分である。尚、Ｉは単位テンソルを表す。
ここで、上記式（３２）を上記式（３１）の右辺に代入すると、以下の式（３３）が得られる。
【数３５】
従って、客観応力速度としてJaumann速度を採用した場合の全成分の構成方程式は、偏差成分の構成式（上記式（２４））と、体積成分の構成式（上記式（３０））とから、以下の式（３４）が成立し、
【数３６】
ここに、一般化Ｍａｘｗｅｌｌ発展モデルの三次元モードでの応力及びひずみの相関式が得られる。
【００３１】
以下に、上記式（３４）の検証を行った。まず、アクリル系ポリマーをシート状にして作製した粘着剤シートから断面積２ｍｍ2の円柱状の試験片を作成し、チャック間１０ｍｍ、引張り速度５，５０，５００ｍｍ／ｍｉｎの三条件で、応力−ひずみ測定を行った。上記式（３４）の一次元モード式は、以下の式（３５）であり（上記式（１２）参照）、
【数３７】
この式（３５）を用いて非線形最小二乗法により、ＡK，ＢK，τKを決定した。
【００３２】
そして、モデルの妥当性を検証するため、ラグランジェ型解法によって検討を行った。解析モデルは、一辺が１０ｍｍの立方体とした。実際の粘着剤は非圧縮に近いポアソン比を示すが、ポアソン比を０．４９などとすると、ロッキングと思われる現象が生じ、解析が進まず、また、要素分割数を増やしても同様に解析が進まないので、まずはモデルの妥当性を検証することを目的として、１要素、ポアソン比：０．３である図５のような解析モデルを設定して検証実験を行った。
【００３３】
その結果、以下の表３及び図６から明らかなように、実測とシミュレーションの値はよい相関性を示すことがわかった。
【表３】
このように、一般化Ｍａｘｗｅｌｌ発展モデルは、ゴム弾性−粘弾性の物性を示す粘着剤（弾性材料））の変形を表現するモデルとして、妥当性があるモデルであることが確認できたわけである。そして、以上により、粘弾性については従来と同じようにして評価をし、ゴム弾性については、今回提案する方法で評価をすることで、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現することができるようになる。
【００３４】
ここで、本実施形態に係るソフトウェアは、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式（一次元モードにおける上記式（９）、三次元モードにおける上記式（２１））と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式（一次元モードにおける上記式（３）、三次元モードにおける上記式（１３））と、から求められる、上記物性に応じて異なる係数をパラメータ（一次元モードにおけるＡi，Ｂi，τi、三次元モードにおけるＡK，ＢK，τK）して含む、応力及びひずみの相関式（一次元モードにおける上記式（１２）及び式（３５）、三次元モードにおける式（３４））を、応力−ひずみ曲線式として出力するためのプログラムであって、コンピュータを、上記パラメータ（一次元モードにおけるＡi，Ｂi，τi、三次元モードにおけるＡK，ＢK，τK）を入力する入力手段、該入力されたパラメータによって上記応力及びひずみの相関式を特定する特定手段、及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段、として機能させるものである。尚、三次元モードにおいては、上記式（３４）に挙がっているパラメータＡK，ＢK，τK以外のパラメータも必要に応じて入力可能とするか、予め、コンピュータの記憶部（後述するＲＯＭあるいはメモリ）に記憶されるようになっている。
【００３５】
本実施形態に係る装置は、コンピュータによって構成され、それぞれバス７に接続されたＣＰＵ、ＲＯＭ、ワーキングメモリ、フレームメモリ、データ入出力装置、ハードディスク及びディスプレイを備える。ＲＯＭは、上記プログラムや各種パラメータを記憶し、ワーキングメモリは、ＣＰＵが制御を行うために必要なメモリであり、バッファやレジスタ等を含む。ＣＰＵは、ＲＯＭに記憶されたコンピュータプログラムに従って各種演算や処理を行う。データ入出力装置が上記入力手段、ＣＰＵ及びＲＯＭやメモリからなる制御部が上記特定手段を構成する。出力手段については、応力−ひずみ曲線式によって表現される応力−ひずみ曲線をディスプレイに表示させるならば、ディスプレイが出力手段を構成し、また、応力−ひずみ曲線式を外部に送信するならば、入出力装置が出力手段を構成し、また、応力−ひずみ曲線式を制御部内で保有しておくならば（それを用いるアプリケーションが制御部内にインストールされていて、当該パソコン内でアプリケーションを起動させるケースなど）、制御部自体が出力手段を構成する。
【００３６】
また、粘弾性も一緒に評価する観点から、それに関するものもプログラミングされ、装置に保有させる（コンピュータにインストールする）ようにするのが好ましい。
【００３７】
尚、本発明は、上記実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々の変更が可能である。
【００３８】
例えば、応力−ひずみ曲線式並びにそれの前提となる各種の相関式は、以上に挙げられたものに限定されるものではなく、その他の適切な近似式に基づいて各種の相関式及びそれに基づいて導出される応力−ひずみ曲線式であってもよい。
【００３９】
また、上記実施形態においては、一次元モードでは、弾性係数として、引張り弾性係数Ｅを用い、三次元モードでは、弾性係数として、せん断弾性係数を用いていたが、一次元モードでも、引張り弾性係数の代わりにせん断弾性係数を用いるようにしてもよい。この場合、各式におけるＥはＧに置き換えられる。
【産業上の利用可能性】
【００４０】
本発明は、ゴム弾性と粘弾性とを併せ持つ性質を示す高品質乃至高付加価値な弾性材料を設計するために利用することができる。
【技術分野】
【０００１】
本発明は、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現するために必要となる応力−ひずみ曲線式を出力するためのプログラム及びその装置に関する。
【背景技術】
【０００２】
従来より、弾性材料の物性を評価する指標として粘弾性があり、その粘弾性を表現するモデルとしては、弾性要素及び粘性要素を直列に配置した一般化Ｍａｘｗｅｌｌモデルが用いられる。一般化Ｍａｘｗｅｌｌモデルでは、実際の応力−時間曲線（応力緩和曲線）（図１（ａ）参照）、即ち、弾性材料の粘弾性を、以下の式（１）に示す応力及び時間の相関式（図１（ｂ）参照）で表現することができるため、以下の式（１）に示す応力及び時間の相関式は、応力−時間曲線式（応力緩和曲線式）として用いられる。
【数１】
ここで、σは応力、ｔは時間、Ｅは一般化Ｍａｘｗｅｌｌモデルにおける弾性要素の引張り弾性係数、τは緩和時間、を表す。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、近年は、弾性材料（例えば、粘着剤）の高品質化乃至高付加価値化に伴い、粘弾性だけでは弾性材料の物性を十分に評価できなくなってきている。具体的に言えば、非線形で且つ大変形を起こすような弾性材料であれば、その物性を評価するためには、粘弾性だけでなく、ゴム弾性を指標に加える必要がある。
【０００４】
ところが、一般化Ｍａｘｗｅｌｌモデルでは、ゴム弾性に関する相関式として、式（１）から以下の式（２）に示す応力及びひずみの相関式を導き出すことができるものの、以下の式（２）に示す応力及びひずみの相関式（図１（ｄ）参照）では、実際の応力−ひずみ曲線（図１（ｃ）参照）、即ち、弾性材料のゴム弾性を表現することができない。
【数２】
ここで、ηは一般化Ｍａｘｗｅｌｌモデルにおける粘性要素の粘性係数、εはひずみ、を表す。
【０００５】
そこで、本発明は、上記問題に鑑みてなされたもので、粘弾性モデルである一般化Ｍａｘｗｅｌｌモデルを発展させ、粘弾性−ゴム弾性両方の物性を評価できるモデル（以下、「一般化Ｍａｘｗｅｌｌ発展モデル」という）を作成し、その定式化を行った上で、それを活用できるツールを提供することを課題とする。
【課題を解決するための手段】
【０００６】
上記課題を解決すべく構成されたもので、ツールとしての本発明に係るプログラムは、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力するためのプログラムであって、コンピュータを、前記パラメータを入力する入力手段、該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段、及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段、として機能させることを特徴とする。
【０００７】
かかる構成によれば、物性に応じて異なるパラメータ（必要に応じてその他のパラメータも含む）を入力することで、弾性材料の物性に応じた応力−ひずみ曲線式が出力される。例えば、適宜、パラメータを変更して入力することで、所望の応力−ひずみ曲線式を追求し、また、それに合うような材料設計を行うといったアプリケーションが期待できる。
【０００８】
また、本発明に係るプログラムにおいては、一例として、前記ひずみ及び弾性係数の相関式は、弾性係数をひずみの指数関数で表したものであり、前記パラメータは、指数関数に乗じられる第一パラメータと、指数においてひずみに乗じられる第二パラメータとからなる。
【０００９】
その場合、さらなる一例として、応力及びひずみの相関式は、Cauchy応力テンソルを太字σ、Jaumann速度であることを右下付き(J)、弾性係数をＧ、ひずみ速度テンソルを太字Ｄ、前記第一パラメータをＡ、前記第二パラメータをＢ、ひずみのスカラー量をハット付きε、弾性要素の偏差応力テンソルを太字σ’、緩和時間をτ、対象を圧縮性線形弾性体とした場合の圧力をｐ、体積弾性係数をＫｖ、演算子traceをtr、単位テンソルを太字Ｉで表した場合、以下の式で規定される。
【数３】
【００１０】
また、上記課題を解決すべく構成されたもので、ツールとしての本発明に係る装置は、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力する装置であって、前記パラメータを入力する入力手段と、該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段と、及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段とを含むことを特徴とする。
【発明の効果】
【００１１】
以上の如く、本発明によれば、粘弾性モデルである一般化Ｍａｘｗｅｌｌモデルを発展させ、粘弾性−ゴム弾性両方の物性を評価できる一般化Ｍａｘｗｅｌｌ発展モデルを作成し、その定式化を行った上で、それを活用できるツールを提供することができる。また、そのツールを用いて、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現することができるようになる。
【図面の簡単な説明】
【００１２】
【図１】（ａ）は、実際の応力−時間曲線図、（ｂ）は、一般化Ｍａｘｗｅｌｌモデルに基づく計算解によって表現される応力−時間曲線図、（ｃ）は、実際の応力−ひずみ曲線図、（ｄ）は、一般化Ｍａｘｗｅｌｌモデルに基づく計算解によって表現される応力−ひずみ曲線図、を示す。
【図２】一般化Ｍａｘｗｅｌｌ発展モデルの概念図を示す。
【図３】１要素モデルの概念図を示す。
【図４】三次元モードでの一般化Ｍａｘｗｅｌｌ発展モデルの概念図を示す。
【図５】一般化Ｍａｘｗｅｌｌ発展モデルを検証するための１要素解析モデルの概念図を示す。
【図６】引張り速度５０ｍｍ／ｍｉｎにて実測した実際の応力−ひずみ曲線図と、一般化Ｍａｘｗｅｌｌ発展モデルに基づく計算解によって表現される応力−ひずみ曲線図とを示す。
【発明を実施するための形態】
【００１３】
以下、本発明の一実施形態について、図面を参酌しつつ説明する。
【００１４】
本実施形態は、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力するためのプログラムに言及する。そこで、以下においては、その応力及びひずみの相関式を導出するまでの流れを説明する。
【００１５】
上記式（１）で表される一般化Ｍａｘｗｅｌｌモデルではゴム弾性を表現できないことは上述の如くである。そこで、ゴム弾性も表現できるよう、一般化Ｍａｘｗｅｌｌモデルにおける弾性要素の引張り弾性係数Ｅ（並びに、一般化Ｍａｘｗｅｌｌモデルにおける粘性要素の粘性係数η）をひずみの関数とするモデル（一般化Ｍａｘｗｅｌｌ発展モデル）を想定した（図２参照）。尚、緩和時間τ（＝η(ε)／Ｅ(ε)）は、便宜上、一定とする。
【００１６】
その上で、まず始めに、粘弾性を表現する上記式（１）に対し、実測の応力緩和試験の値から回帰分析を行い、各ひずみにおける弾性係数Ｅと緩和時間τの組を求めた。より詳しくは、アクリル系ポリマーをシート状にして作製した粘着剤シートから断面積２ｍｍ2の円柱状の試験片を作成し、チャック間１０ｍｍ、速度１０００ｍｍ／ｍｉｎにて、所定のひずみになるまで伸張し、次に、ひずみを一定の状態に保持し、応力緩和試験を行い、上記式（１）を用いて回帰分析を行い、各ひずみにおける計五組の弾性係数Ｅ、緩和時間τの組を求めた。その結果を以下の表１に表す。
【表１】
尚、表１において、ひずみ（strain）は、対数ひずみである。公称ひずみで表せば、左から、１００％、２００％、３００％、４００％、５００％、６００％、となっている。
【００１７】
次に、弾性係数Ｅのひずみ依存性データをプロットし、カーブフィッティングを行うと、以下の表２に示す如く、弾性係数はひずみに対し、以下の式（３）に示す指数関数で近似することで、高い相関係数（表２中のＲ2値が１．０００に近いほど近似していることを示す）を持つことがわかった。即ち、ひずみ及び弾性係数の相関式は、弾性係数Ｅをひずみεの指数関数で表した以下の式（３）が妥当であることが確認できた。
【表２】
【数４】
ここで、Ａは指数関数に乗じられる第一パラメータ、Ｂは指数においてひずみεに乗じられる第二パラメータ、を表す。
従って、一般化Ｍａｘｗｅｌｌモデルでは、取得すべき物性値は、τｉ、Ｅｉである（上記式（１）参照）のに対し、一般化Ｍａｘｗｅｌｌ発展モデルでは、取得すべき物性値は、τｉ、Ａｉ、Ｂｉである。
【００１８】
次に、Ｍａｘｗｅｌｌ要素が一つしかない場合を考える。図３に示すように、応力σ、ひずみε、弾性要素の弾性係数Ｅ、粘性要素の粘性係数ηとすると、
【数５】
【数６】
【数７】
【数８】
となる。そして、上記式（６）を時間で微分すると、以下の式（８）が得られる。
【数９】
また、上記式（５）を時間で微分した式と上記式（７）、上記式（８）を用いて、εa、εbを消去して整理すると、以下の式（９）が得られる。
【数１０】
【００１９】
ここで、上記式（３）において、Ａ，Ｂを定数とし、
【数１１】
とおいて整理すると、以下の式（１０）が得られる。
【数１２】
これが図３に示す１要素モデルの構成式である。
【００２０】
次に、図２に示すような一般型の構成式を誘導する。一般型になっても、各要素では、上記式（１０）が成立するため、以下の式（１１）が成立する。
【数１３】
本モデルでは、各要素のひずみは一定であるため、ひずみに対して添え字ｉは付かない。そこで、一般型モデルの全体の応力σは各要素のσiの和であるから、以下の式（１２）が成立し、
【数１４】
一般化Ｍａｘｗｅｌｌ発展モデルの応力及びひずみの相関式が得られる。
【００２１】
ところで、上記式（１２）は、一次元モードの式であり、応力及びひずみの相関式を概念として表すためのものであるが、シミュレーションを行うためには、三次元モードでの定式化を行わなければならない。そこで、以下においては、三次元モードでの式を導出するまでの流れを説明する。
【００２２】
まず、一次元モードの場合と同様、三次元の観点で一般化Ｍａｘｗｅｌｌ発展モデルを図４のように設定すると、（弾性材料の物性に応じて異なる係数を含む、）ひずみ及び弾性係数の相関式は、以下の式（１３）のようになる。
【数１５】
ここで、Ｇはせん断弾性係数、ハット付きεはひずみのスカラー量、を表す。
尚、一次元モードの場合は、添え字ｉを右下付きで記載していたが、以下の三次元モードにおいては、添え字Ｋを右上付きで記載している。
【００２３】
ひずみのスカラー量（ハット付きε）を微小ひずみテンソル（太字ε）を用いて以下の式（１４）のように置くこととする。
【数１６】
今、Ｋ番目ユニットの弾性要素の変位応力テンソル（太字σ’）は、以下の式（１５）となる。
【数１７】
ここで、添え字がｓｐＫのεはＫ番目ユニットの弾性要素の微小ひずみテンソル、を表す。
【００２４】
上記式（１５）の両辺を物質時間微分すれば、以下の式（１６）となる。
【数１８】
ここで、太字Ｄはひずみ速度テンソル、添え字がｓｐＫの太字ＤはＫ番目ユニットの弾性要素のひずみ速度テンソル、を表す。
上記式（１６）の左辺は客観性のないテンソルであるため、客観応力速度を導入すると、以下の式（１７）を弾性要素の構成方程式として用いることができる。
【数１９】
ここで、右下付きの（＊）は任意の客観応力速度を表す。
【００２５】
一方、Ｋ番目ユニットの粘性要素における偏差成分の構成方程式は、以下の式（１８）で与えられる。
【数２０】
Ｋ番目ユニットの系全体の微小ひずみテンソルを太字ε、Ｋ番目ユニットの微小ひずみテンソルを添え字がＫの太字ε、で表すと、以下の式（１９）となり、
【数２１】
これを物質時間微分すると、以下の式（２０）となる。
【数２２】
【００２６】
上記式（１７）、（１８）、（２０）より、Ｋ番目ユニットの偏差成分の構成方程式は、以下の式（２１）のように導出される。
【数２３】
ここで、上記式（１３）は、以下の式（２２）となることから、
【数２４】
上記式（２１）は、以下の式（２３）となる。
【数２５】
モデルの性質から、モデル全体の応力はユニットＫ個の応力の総和であることから、偏差成分の応力は、以下の式（２４）となる。
【数２６】
【００２７】
次に、体積成分を考える。ここで、対象を圧縮性線形弾性体と仮定すると、圧力ｐは以下の式（２５）となる。
【数２７】
ここで、Ｋｖは体積弾性係数、trは演算子trace、を表し、体積弾性係数Ｋｖは、以下の式（２６）となる。
【数２８】
ここで、νはポアソン比を表す。
【００２８】
上記式（２５）を物質時間微分すると、以下の式（２７）となる。
【数２９】
また、以下の式（２８）であり、
【数３０】
さらに上記式（２２）を代入することにより、以下の式（２９）となる。
【数３１】
従って、体積成分の構成方程式は以下の式（３０）となる。
【数３２】
【００２９】
客観応力速度としてJaumann速度を用いることを考えると、Cauchy応力（太字σ）の物質時間導関数とそのJaumann速度の間には、以下の式（３１）の関係が成り立つ。
【数３３】
ここで、太字Ｗはスピンテンソル、右下付き(J)はJaumann速度であることを表す。
【００３０】
Cauchy応力（太字σ）を偏差成分と体積成分に分離すると、以下の式（３２）となる。
【数３４】
即ち、右辺の第一項が偏差成分であり、第二項が体積成分である。尚、Ｉは単位テンソルを表す。
ここで、上記式（３２）を上記式（３１）の右辺に代入すると、以下の式（３３）が得られる。
【数３５】
従って、客観応力速度としてJaumann速度を採用した場合の全成分の構成方程式は、偏差成分の構成式（上記式（２４））と、体積成分の構成式（上記式（３０））とから、以下の式（３４）が成立し、
【数３６】
ここに、一般化Ｍａｘｗｅｌｌ発展モデルの三次元モードでの応力及びひずみの相関式が得られる。
【００３１】
以下に、上記式（３４）の検証を行った。まず、アクリル系ポリマーをシート状にして作製した粘着剤シートから断面積２ｍｍ2の円柱状の試験片を作成し、チャック間１０ｍｍ、引張り速度５，５０，５００ｍｍ／ｍｉｎの三条件で、応力−ひずみ測定を行った。上記式（３４）の一次元モード式は、以下の式（３５）であり（上記式（１２）参照）、
【数３７】
この式（３５）を用いて非線形最小二乗法により、ＡK，ＢK，τKを決定した。
【００３２】
そして、モデルの妥当性を検証するため、ラグランジェ型解法によって検討を行った。解析モデルは、一辺が１０ｍｍの立方体とした。実際の粘着剤は非圧縮に近いポアソン比を示すが、ポアソン比を０．４９などとすると、ロッキングと思われる現象が生じ、解析が進まず、また、要素分割数を増やしても同様に解析が進まないので、まずはモデルの妥当性を検証することを目的として、１要素、ポアソン比：０．３である図５のような解析モデルを設定して検証実験を行った。
【００３３】
その結果、以下の表３及び図６から明らかなように、実測とシミュレーションの値はよい相関性を示すことがわかった。
【表３】
このように、一般化Ｍａｘｗｅｌｌ発展モデルは、ゴム弾性−粘弾性の物性を示す粘着剤（弾性材料））の変形を表現するモデルとして、妥当性があるモデルであることが確認できたわけである。そして、以上により、粘弾性については従来と同じようにして評価をし、ゴム弾性については、今回提案する方法で評価をすることで、ゴム弾性と粘弾性とを併せ持つ物性を示す弾性材料の大変形挙動をシミュレーション上で定量性よく表現することができるようになる。
【００３４】
ここで、本実施形態に係るソフトウェアは、弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式（一次元モードにおける上記式（９）、三次元モードにおける上記式（２１））と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式（一次元モードにおける上記式（３）、三次元モードにおける上記式（１３））と、から求められる、上記物性に応じて異なる係数をパラメータ（一次元モードにおけるＡi，Ｂi，τi、三次元モードにおけるＡK，ＢK，τK）して含む、応力及びひずみの相関式（一次元モードにおける上記式（１２）及び式（３５）、三次元モードにおける式（３４））を、応力−ひずみ曲線式として出力するためのプログラムであって、コンピュータを、上記パラメータ（一次元モードにおけるＡi，Ｂi，τi、三次元モードにおけるＡK，ＢK，τK）を入力する入力手段、該入力されたパラメータによって上記応力及びひずみの相関式を特定する特定手段、及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段、として機能させるものである。尚、三次元モードにおいては、上記式（３４）に挙がっているパラメータＡK，ＢK，τK以外のパラメータも必要に応じて入力可能とするか、予め、コンピュータの記憶部（後述するＲＯＭあるいはメモリ）に記憶されるようになっている。
【００３５】
本実施形態に係る装置は、コンピュータによって構成され、それぞれバス７に接続されたＣＰＵ、ＲＯＭ、ワーキングメモリ、フレームメモリ、データ入出力装置、ハードディスク及びディスプレイを備える。ＲＯＭは、上記プログラムや各種パラメータを記憶し、ワーキングメモリは、ＣＰＵが制御を行うために必要なメモリであり、バッファやレジスタ等を含む。ＣＰＵは、ＲＯＭに記憶されたコンピュータプログラムに従って各種演算や処理を行う。データ入出力装置が上記入力手段、ＣＰＵ及びＲＯＭやメモリからなる制御部が上記特定手段を構成する。出力手段については、応力−ひずみ曲線式によって表現される応力−ひずみ曲線をディスプレイに表示させるならば、ディスプレイが出力手段を構成し、また、応力−ひずみ曲線式を外部に送信するならば、入出力装置が出力手段を構成し、また、応力−ひずみ曲線式を制御部内で保有しておくならば（それを用いるアプリケーションが制御部内にインストールされていて、当該パソコン内でアプリケーションを起動させるケースなど）、制御部自体が出力手段を構成する。
【００３６】
また、粘弾性も一緒に評価する観点から、それに関するものもプログラミングされ、装置に保有させる（コンピュータにインストールする）ようにするのが好ましい。
【００３７】
尚、本発明は、上記実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々の変更が可能である。
【００３８】
例えば、応力−ひずみ曲線式並びにそれの前提となる各種の相関式は、以上に挙げられたものに限定されるものではなく、その他の適切な近似式に基づいて各種の相関式及びそれに基づいて導出される応力−ひずみ曲線式であってもよい。
【００３９】
また、上記実施形態においては、一次元モードでは、弾性係数として、引張り弾性係数Ｅを用い、三次元モードでは、弾性係数として、せん断弾性係数を用いていたが、一次元モードでも、引張り弾性係数の代わりにせん断弾性係数を用いるようにしてもよい。この場合、各式におけるＥはＧに置き換えられる。
【産業上の利用可能性】
【００４０】
本発明は、ゴム弾性と粘弾性とを併せ持つ性質を示す高品質乃至高付加価値な弾性材料を設計するために利用することができる。
【特許請求の範囲】
【請求項１】
弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力するためのプログラムであって、
コンピュータを、
前記パラメータを入力する入力手段、
該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段、
及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段、
として機能させることを特徴とするプログラム。
【請求項２】
前記ひずみ及び弾性係数の相関式は、弾性係数をひずみの指数関数で表したものであり、前記パラメータは、指数関数に乗じられる第一パラメータと、指数においてひずみに乗じられる第二パラメータとからなる請求項１に記載のプログラム。
【請求項３】
応力及びひずみの相関式は、Cauchy応力テンソルを太字σ、Jaumann速度であることを右下付き(J)、弾性係数をＧ、ひずみ速度テンソルを太字Ｄ、前記第一パラメータをＡ、前記第二パラメータをＢ、ひずみのスカラー量をハット付きε、弾性要素の偏差応力テンソルを太字σ’、緩和時間をτ、対象を圧縮性線形弾性体とした場合の圧力をｐ、体積弾性係数をＫｖ、演算子traceをtr、単位テンソルを太字Ｉで表した場合、以下の式
【数３８】
で規定される請求項２に記載のプログラム。
【請求項４】
弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力する装置であって、
前記パラメータを入力する入力手段と、
該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段と、
及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段と
を含むことを特徴とする装置。
【請求項１】
弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力するためのプログラムであって、
コンピュータを、
前記パラメータを入力する入力手段、
該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段、
及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段、
として機能させることを特徴とするプログラム。
【請求項２】
前記ひずみ及び弾性係数の相関式は、弾性係数をひずみの指数関数で表したものであり、前記パラメータは、指数関数に乗じられる第一パラメータと、指数においてひずみに乗じられる第二パラメータとからなる請求項１に記載のプログラム。
【請求項３】
応力及びひずみの相関式は、Cauchy応力テンソルを太字σ、Jaumann速度であることを右下付き(J)、弾性係数をＧ、ひずみ速度テンソルを太字Ｄ、前記第一パラメータをＡ、前記第二パラメータをＢ、ひずみのスカラー量をハット付きε、弾性要素の偏差応力テンソルを太字σ’、緩和時間をτ、対象を圧縮性線形弾性体とした場合の圧力をｐ、体積弾性係数をＫｖ、演算子traceをtr、単位テンソルを太字Ｉで表した場合、以下の式
【数３８】
で規定される請求項２に記載のプログラム。
【請求項４】
弾性要素及び粘性要素を直列に配置したＭａｘｗｅｌｌモデルに基づいて求められる、応力、ひずみ、弾性係数及び緩和時間の相関式と、物性に応じて異なる係数を含む、ひずみ及び弾性係数の相関式と、から求められる、前記物性に応じて異なる係数をパラメータとして含む、応力及びひずみの相関式を、応力−ひずみ曲線式として出力する装置であって、
前記パラメータを入力する入力手段と、
該入力されたパラメータによって前記応力及びひずみの相関式を特定する特定手段と、
及び、該特定された応力及びひずみの相関式を、応力−ひずみ曲線式として出力する出力手段と
を含むことを特徴とする装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図２】
【図３】
【図４】
【図５】
【図６】
【公開番号】特開２０１０−２５６２９３（Ｐ２０１０−２５６２９３Ａ）
【公開日】平成２２年１１月１１日（２０１０．１１．１１）
【国際特許分類】
【出願番号】特願２００９−１０９６０５（Ｐ２００９−１０９６０５）
【出願日】平成２１年４月２８日（２００９．４．２８）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 研究集会名：日本機械学会第２１回計算力学講演会 主催者名：社団法人日本機械学会 計算力学部門 開催日：２００８年１１月１日（土）〜３日（月）
【出願人】（０００００３９６４）日東電工株式会社 (5,557)
【出願人】（５０４１３６５６８）国立大学法人広島大学 (924)
【Ｆターム（参考）】
【公開日】平成２２年１１月１１日（２０１０．１１．１１）
【国際特許分類】
【出願日】平成２１年４月２８日（２００９．４．２８）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 研究集会名：日本機械学会第２１回計算力学講演会 主催者名：社団法人日本機械学会 計算力学部門 開催日：２００８年１１月１日（土）〜３日（月）
【出願人】（０００００３９６４）日東電工株式会社 (5,557)
【出願人】（５０４１３６５６８）国立大学法人広島大学 (924)
【Ｆターム（参考）】
[ Back to top ]