動的弾性率のソースを表示

'''動的弾性率'''（どうてきだんせいりつ、{{Lang-en-short|Dynamic modulus, Dynamic Elastic Modulus}}）<ref name=Meyers>Meyers and Chawla (1999): "Mechanical Behavior of Materials," 98-103.</ref>は、[[物体]]の[[粘弾性]]を記述する[[物理量]]の一つで、弾性率（[[ヤング率]]）を拡張した概念である。「振動する[[応力]]」と、それによって生じた[[ひずみ|歪み]]の[[フェーザ表示]]による「比」として定義される。[[複素数]]で表わされることが多いため、'''複素弾性率'''（ふくそだんせいりつ、{{Lang|en|Complex modulus}}）と呼ばれることもある。動的弾性率は、動的な粘弾性特性を、粘弾性物質の応力-ひずみ特性の[[位相]]遅れに着目して複素弾性率の[[複素数#極形式|偏角]]として表現したものであり、複素弾性率の実数部にあたる「[[貯蔵弾性率]]」、虚数部にあたる「[[損失弾性率]]」の2つの項に分解できる。同様に、[[剛性率]]についても同様に[[動的剛性率]]が定義されるので、併せて説明する。

数式的な取扱いにおいて、電気工学で用いられる[[インピーダンス]]や、[[制御工学]]の[[周波数伝達関数]]に良く似ており、一種の[[アナロジー]]が成立する。

==定義==
===動的弾性率===
以下の{{EquationNote|式1-1}}のような応力が物体に印加されたとき、以下の{{EquationNote|式1-2}}のような変形が生じた、即ち、応力に対し、その応答たる歪みに位相遅れ {{Mvar|δ}} が生じたとした場合を考える。

{{NumBlk|*|応力:<math>\sigma = \sigma_0 \cos(t\omega) \,</math><ref name="Meyers"/>|{{EquationRef|式1-1}}}}
{{NumBlk|*|歪み:<math>\varepsilon = \varepsilon_0 \cos(t\omega  + \delta)</math>|式1-2}}

貯蔵弾性率 {{Mvar|E&prime;}} ({{Lang|en|storage modulus}})と損失弾性率 {{Mvar|E″}} ({{Lang|en|loss modulus}}) を以下のように導入する。

{{NumBlk|*|貯蔵弾性率: <math> E'  = \frac {\sigma_0} {\varepsilon_0} \cos \delta </math>|式1-3}}
{{NumBlk|*|損失弾性率: <math> E'' =  \frac {\sigma_0} {\varepsilon_0} \sin \delta </math><ref name=Meyers/>|式1-4}}

但し、{{Math|1=''&omega;'' = 2''&pi;f''}} であり、{{Mvar|f}} は周期的応力の[[振動数]]、{{Mvar|t}} は時刻、 {{Mvar|δ}} は応力 ({{Lang|en|stress}})と歪み ({{Lang|en|strain}}) の間の位相遅れを意味する。

さらに、複素弾性率 ({{Lang|en|complex modulus}}) を、以下のように定義する。
{{NumBlk|:|<math>{E}^{*} = E' + i E'' </math>|式1-5}}
但し、{{Mvar|i}} は、[[虚数単位]]を表す。動的弾性率とは、この複素弾性率を表す場合が多いが、複素弾性率、貯蔵弾性率、損失弾性率の少なくともいずれかを指すこともある。

応力、歪の両方に、虚数項を付け加え、
{{NumBlk|*|複素化応力:<math> {\sigma}^{*} = \sigma_0 \exp(it\omega) \,</math> <ref name=Meyers/>|式1-6}}
{{NumBlk|*|複素化歪　:<math> {\varepsilon}^{*} = \varepsilon_0 \exp(i(t\omega  + \delta))</math>|式1-7}}
とすると、

{{NumBlk|:|<math> {\sigma}^{*}={E}^{*} {\varepsilon}^{*}</math>|式1-8}}

が成り立つ。この意味で、{{Math|''E''<sup>*</sup>}}は、応力と歪みの比となっている。

===動的剛性率===
剛性率についても、同様の
:<math>G^* = G' + iG'' \,</math> <ref name=Meyers/>
が定義される。

==粘弾性体による応力-ひずみ間の位相遅れ==
日常的に目にする粘弾性体の例としては[[ゼリー]]がある。ゼリーの"ぷるるん"とした性質に粘弾性の本質がある。この”ぷるるん”とした性質は、粘弾性体による応力-ひずみ間の位相遅れの帰結であり、これを定量的かつシステマティックに表現したものが、動的弾性率である。従って、動的弾性率を理解する前に、粘弾性体による応力-ひずみ間の位相遅れについて理解しておく必要がある。

皿の上にゼリーを置いて、少し揺らすと、ゼリーは”ぷるるん”と震える （ここでは、[[インパルス応答]]を与えたとしても、[[正弦波]]的な振動（後述の{{EquationNote|式2-1}}を参照）を与えてその応答をみたとしても、どちらでもかまわない）。ここで、より「シャキッとした」（[[弾性体]]に近い）[[ゼリー]]は、”ゆらし”に対して、より機敏に応答するであろう。より「フニャっとした」（[[粘性]]体に近い）ゼリーは、”揺らし”に対して少し遅れた（位相の遅れを伴った）応答するであろう。このように、我々は、経験的に食べる前にゼリーの食感を推定する方法を知っている<ref name=zeri>http://www2.kobe-u.ac.jp/~komoda/kyougen/komoda_05.pdf</ref>。

簡単に言えば、上記の日常的な、食感推定法をより精密化したものの一つが、「{{仮リンク|動的粘弾性測定|en|Dynamic mechanical analysis}}」と言われる測定手法である。動的粘弾性測定では、粘弾性物質に対し、以下の{{EquationNote|式2-1}}のような[[正弦波]]状の"振動する応力"を印加をし、それに応答して生じた変形（歪み）の大きさが、リアルタイムで測定される<ref name=pe>[http://las.perkinelmer.com/content/ApplicationNotes/APP_FilmsandCoatings.pdf PerkinElmer "Mechanical Properties of Films and Coatings"]</ref>。

このとき、測定されるひずみの位相遅れについては、以下のような、2種類の極端なケースが考えられる。
*完全な[[弾性体]]においては、印加された応力と、それに対する応答として生じたひずみの間に位相の遅れが生じない。(後述の{{EquationNote|式2-2}}において、{{Mvar|δ}}=0の場合に相当する)
*完全な[[粘性]]体においては、印加された応力と、それに対する応答として生じたひずみの間に位相の遅れが90度(π/2 rad)分生じる。(後述の{{EquationNote|式2-2}}において、{{Mvar|δ}}=90度(π/2 rad)の場合に相当する)日常的な現象においては、「シャキッとした」ゼリーは前者に近く、「フニャっとした」ゼリーは後者に近いということになる。より一般的に、物質が弾性体に近いときは上記の「位相の遅れ」（{{Mvar|δ}}）が0に近く、粘性体に近いときには、上記の「位相の遅れ」が、90度に近いということになる。このように、物質が粘性体または弾性体に近いのかは、その物質に一定のひずみを与えたときの応力緩和（応力の時間変化）の位相遅れや、緩和時間をから判別できる。

粘弾性物質の性質を、物質の弾性(変形を元に戻そうとする性質)による効果と、粘性（変形を抑制する性質）の”重ね合わせ”として考えることにすると（同一周期の正弦波の重ね合わせは位相の変化として表現される）、理想的な（物質からの応答によって振動周期が変化しない）粘弾性物質においては、以下の{{EquationNote|式2-1}}のような応力 ({{Lang|en|stress}}) に対して、物質の変形は、{{EquationNote|式2-2}}のように表現されることになる。

{{NumBlk|*|応力:<math> \sigma = \sigma_0 \cos(t\omega) \,</math> <ref name=Meyers/>|{{EquationRef|式2-1}}}}
{{NumBlk|*|歪:<math> \varepsilon = \varepsilon_0 \cos(t\omega  + \delta)</math>|{{EquationRef|式2-2}}}}

ここで、{{Math|1=''&omega;'' = 2''&pi;f''}} であり、{{Mvar|f}} は周期的応力の振動数、{{Mvar|t}} は時刻、 {{Mvar|δ}} は応力 ({{Lang|en|stress}})と歪み ({{Lang|en|strain}}) の間の位相遅れを意味する。

さて、位相の遅れは、エネルギー損失の帰結である。実際、微小な変形 {{Math|d''ε''}} の間にうける間に加わる力は {{Math|''σ''(''t'')}} であるため、この間のエネルギー変化は、
{{NumBlk|:|<math>\mathrm dU=\sigma (t) d\epsilon = \sigma (t) \frac{\mathrm d\epsilon}{\mathrm dt}\mathrm dt =  -\varepsilon_0  \sigma_0 \omega \cos(t\omega) \sin(t\omega + \delta)\mathrm dt</math>|式2-3}}
となり、加法定理
{{NumBlk|:|{{Math|1=sin(''a''+''b'') = sin(''a'')cos(''b'') + cos(''a'')sin(''b'')}}|式2-4a}}
{{NumBlk|:|{{Math|1=cos(''a''+''b'') = cos(''a'')cos(''b'') &minus; sin(''a'')sin(''b'')}}|式2-4b}}
より、
{{NumBlk|:|{{Math|1=2sin(''a'')cos(''b'') = (sin(''a''+''b'') +sin(''a''&minus;''b''))}}|式2-4c}}
であるため、
{{NumBlk|:|<math>\mathrm dU=-\frac{1}{2}\varepsilon_0  \sigma_0 \omega 
(\sin(2t\omega + \delta )  + \sin(-\delta)) \mathrm dt</math>|式2-5}}

であり、一方で[[三角関数]]を一周期 {{Math|1=''T''=1/''f''}} に渡って積分すると０になるため、一周期毎のエネルギー損失 {{Mvar|U}} は、
{{NumBlk|:|<math> U = {\int}_{0}^{T} \mathrm dU = -\frac{1}{2}\varepsilon_{0} \sigma_0 \omega 
{\int}_{0}^{T}  \left[\sin(2t\omega + \delta) + \sin(-\delta)\right] \mathrm dt
=\varepsilon_{0}\sigma_{0} \pi \sin(\delta)</math>|式2-6}}

となる。

==脚注==
<references />

==関連項目==
*{{仮リンク|動的粘弾性測定|en|Dynamic mechanical analysis}}
*[[粘弾性]]
*[[弾性率]]
*[[ヤング率]]
*[[非ニュートン流体]]
*[[デボラ数]]

{{DEFAULTSORT:とうてきたんせいりつ}}
[[Category:物理量]]
[[Category:固体力学]]
[[Category:粘度]]
[[Category:レオロジー]]
[[Category:連続体力学]]
[[Category:弾性]]
[[Category:変形]]