材料科学と連続体力学において、粘弾性とは、変形時に粘性と弾性の両方の特性を示す材料の特性のことです。水などの粘性材料は、応力が加えられると、せん断流と歪みの両方に時間とともに直線的に抵抗します。弾性材料は、伸ばされると歪み、応力が除去されるとすぐに元の状態に戻ります。
粘弾性材料は、これら両方の特性の要素を持ち、時間依存の歪みを示します。弾性は通常、秩序だった固体の結晶面に沿った結合の伸張の結果であるのに対し、粘性は非晶質材料内の原子または分子の拡散の結果です。[1]
背景
19 世紀には、ジェームズ・クラーク・マクスウェル、ルートヴィヒ・ボルツマン、ケルビン卿などの物理学者がガラス、金属、ゴムのクリープと回復について研究し、実験しました。粘弾性は、合成ポリマーが設計され、さまざまな用途に使用されるようになった 20 世紀後半にさらに研究されました。 [2]粘弾性の計算は、粘度変数ηに大きく依存します。η の逆数は流動性φとしても知られています。どちらの値も、温度の関数として、または与えられた値 (ダッシュポットの場合など) として導出できます。[1]

材料内部のひずみ速度と応力の変化に応じて、粘度は線形、非線形、または塑性応答に分類できます。材料が線形応答を示す場合、ニュートン流体として分類されます。この場合、応力はひずみ速度に線形比例します。材料がひずみ速度に対して非線形応答を示す場合、非ニュートン流体として分類されます。せん断/ひずみ速度が一定のままで粘度が減少する興味深いケースもあります。このタイプの動作を示す材料はチキソトロピーとして知られています。さらに、応力がこのひずみ速度に依存しない場合、材料は塑性変形を示します。[1]多くの粘弾性材料は、ポリマー弾性の熱力学理論によって説明されるゴムのような動作を示します。
粘弾性材料の例としては、非晶質ポリマー、半結晶性ポリマー、バイオポリマー、非常に高温の金属、ビチューメン材料などがあります。ひび割れは、ひずみが急速に加えられ、弾性限界を超えた場合に発生します。靭帯と腱は粘弾性であるため、靭帯と腱に生じる可能性のある損傷の程度は、靭帯と腱の長さの変化率と加えられた力の両方に依存します。[引用が必要]
粘弾性材料には次の特性があります。
- 応力-ひずみ曲線にヒステリシスが見られる
- 応力緩和が起こる:一定のひずみにより応力が減少する
- クリープが発生する:一定の応力が段階的に加わり、ひずみが増加する
- その剛性はひずみ速度または応力速度に依存する
弾性挙動と粘弾性挙動
純粋に弾性的な物質とは異なり、粘弾性物質には弾性成分と粘性成分があります。粘弾性物質の粘性により、物質は時間に対するひずみ速度依存性を持ちます。純粋に弾性的な物質は、荷重が加えられてから取り除かれたときにエネルギー (熱) を消散しません。しかし、粘弾性物質は、荷重が加えられてから取り除かれたときにエネルギーを消散します。応力-ひずみ曲線にはヒステリシスが見られ、ループの面積は荷重サイクル中に失われたエネルギーに等しくなります。粘性は熱活性化塑性変形に対する抵抗であるため、粘性物質は荷重サイクルを通じてエネルギーを失います。塑性変形はエネルギーの損失をもたらしますが、これは純粋に弾性的な物質の荷重サイクルに対する反応の特徴ではありません。[1]
具体的には、粘弾性は分子の再配置です。ポリマーなどの粘弾性材料に応力がかかると、長いポリマー鎖の一部の位置が変わります。この動きまたは再配置はクリープと呼ばれます。ポリマーは、これらの鎖の一部が応力に対応するために再配置されている場合でも固体材料のままであり、これが発生すると、材料に逆応力が発生します。逆応力が適用された応力と同じ大きさになると、材料はクリープしなくなります。元の応力が取り除かれると、蓄積された逆応力によってポリマーは元の形状に戻ります。材料がクリープすると接頭辞に粘性-が付き、材料が完全に回復すると接尾辞に弾性が付きます。[2]
線形粘弾性と非線形粘弾性
線形粘弾性は、クリープ応答と荷重の両方で関数が分離可能な場合です。すべての線形粘弾性モデルは、応力とひずみを結び付けるボルテラ方程式で表すことができます。 または 、
線形粘弾性は通常、小さな変形にのみ適用できます。
非線形粘弾性は、関数が分離できない場合に発生します。これは通常、変形が大きい場合、または材料が変形によって特性を変える場合に発生します。非線形粘弾性は、粘弾性流体における法線応力、せん断減粘、伸長肥厚などの観察される現象も解明します。[3]
非弾性材料は粘弾性材料の特殊なケースです。非弾性材料は、荷重が除去されると完全に元の状態に戻ります。
弾性、粘性、粘弾性挙動の形態を区別する場合、観察対象物質の緩和時間に対する測定の時間スケール、いわゆるデボラ数(De)を参照すると便利です。[3] ここで 、
- 材料の緩和時間である
- 時間です
動的弾性率
粘弾性は、小さな振動応力を加え、その結果生じる歪みを測定する 動的機械分析を使用して研究されます。
- 純粋に弾性のある材料では、応力と歪みが同位相にあるため、一方が他方に引き起こす反応は即時に起こります。
- 純粋に粘性のある材料では、ひずみは応力より 90 度遅れて発生します。
- 粘弾性材料は、これら 2 種類の材料の中間の挙動を示し、ひずみに若干の遅れが生じます。
複素動的弾性率G は、振動する応力とひずみの関係を表すために使用できます。 ここで、は貯蔵弾性率、は損失弾性率です。 ここで、と はそれぞれ応力とひずみの振幅、はそれらの間の位相シフトです。
線形粘弾性の構成モデル

非晶質ポリマー、半結晶性ポリマー、生体ポリマー、さらには生体組織や細胞などの粘弾性材料[4]は、応力とひずみ、または力と変位の相互作用や時間依存性を決定するためにモデル化することができます。マクスウェルモデル、ケルビン・フォークトモデル、標準線形固体モデル、バーガースモデルなどのこれらのモデルは、さまざまな荷重条件下での材料の応答を予測するために使用されます。
粘弾性挙動には、弾性コンポーネントと粘性コンポーネントがあり、それぞれバネとダッシュポットの線形結合としてモデル化されます。各モデルはこれらの要素の配置が異なり、これらの粘弾性モデルはすべて電気回路として同等にモデル化できます。
同等の電気回路では、応力は電流で表され、ひずみ率は電圧で表されます。バネの弾性率は回路のインダクタンスの逆数(エネルギーを蓄える)に類似しており、ダッシュポットの粘性は回路の抵抗(エネルギーを消散させる)に類似しています。
前述のように、弾性コンポーネントは、次の式で弾性定数 E のスプリングとしてモデル化できます。 ここで、σ は応力、E は材料の弾性率、ε はフックの法則 と同様に、与えられた応力下で発生するひずみです。
粘性成分はダッシュポットとしてモデル化することができ、応力とひずみ速度の関係は次のように表すことができます。 ここで、σ は応力、η は材料の粘度、dε/dt はひずみの時間微分です。
応力とひずみの関係は、特定の応力またはひずみ速度に対して単純化できます。高応力またはひずみ速度/短時間の場合、応力-ひずみ関係の時間微分成分が支配的になります。これらの条件では、ダッシュポットは変形せずに高荷重に耐えることができる剛性ロッドとして近似できます。したがって、ダッシュポットは「短絡」と見なすことができます。[5] [6]
逆に、低応力状態/より長い期間では、時間微分成分は無視でき、ダッシュポットはシステムから効果的に除去できます(「オープン」回路)。[6]その結果、ダッシュポットに並列に接続されたスプリングのみがシステム全体の歪みに寄与します。[5]
マクスウェルモデル

マクスウェル モデルは、図に示すように、純粋に粘性のダンパーと純粋に弾性のスプリングを直列に接続して表すことができます。このモデルは次の式で表すことができます。
このモデルでは、材料が一定のひずみを受けると、応力は徐々に緩和します。材料が一定の応力を受けると、ひずみには 2 つの成分があります。まず、弾性成分が瞬時に発生し、バネに対応し、応力が解放されるとすぐに緩和します。2 つ目は粘性成分で、応力が加えられている限り時間とともに増大します。マクスウェル モデルは、応力が時間とともに指数関数的に減少すると予測しており、これはほとんどのポリマーで正確です。このモデルの限界の 1 つは、クリープを正確に予測できないことです。クリープまたは一定応力条件のマクスウェル モデルは、ひずみが時間とともに直線的に増加すると仮定しています。ただし、ほとんどのポリマーでは、ひずみ率が時間とともに減少することが示されています。[2]
このモデルは、融点付近の熱可塑性ポリマー、生コンクリート(経年変化は無視)、融点に近い温度の多数の金属などの柔らかい固体に適用できます。
しかし、ここで紹介した方程式は、よりミクロなモデルからの一貫した導出を欠いており、観察者に依存しません。上部対流マクスウェルモデルは、コーシー応力テンソルの観点からその妥当な定式化であり、粘弾性の最も単純なテンソル構成モデルを構成します(例えば、[7]または[6]を参照 )。
ケルビン・フォークトモデル

ケルビン-フォークト モデルはフォークト モデルとも呼ばれ、図に示すように、並列に接続されたニュートン ダンパーとフック弾性スプリングで構成されています。このモデルは、ポリマーのクリープ挙動を説明するために使用されます。
構成関係は線形一次微分方程式として表現されます。
このモデルは、可逆的な粘弾性ひずみを受ける固体を表します。一定の応力を加えると、材料は減少率で変形し、定常ひずみに漸近します。応力が解放されると、材料は徐々に変形していない状態に戻ります。一定の応力(クリープ)では、時間が無限に続くにつれてひずみが σ/E に近づくと予測するため、このモデルは非常に現実的です。マクスウェルモデルと同様に、ケルビン-フォークトモデルにも限界があります。このモデルは材料のクリープのモデル化に非常に優れていますが、緩和に関してはモデルの精度ははるかに劣ります。[8]
このモデルは、負荷がそれほど高くない場合、有機ポリマー、ゴム、木材に適用できます。
標準線形ソリッドモデル
標準線形ソリッド モデル (別名 Zener モデル) は、2 つのスプリングと 1 つのダッシュポットで構成されます。これは、粘弾性材料のクリープと応力緩和の両方の挙動を適切に記述する最も単純なモデルです。このモデルでは、支配的な構成関係は次のとおりです。
一定の応力を受けると、モデル化された材料は瞬間的に変形して、ひずみの瞬間弾性部分であるひずみになります。その後、材料は変形を続け、ひずみの遅延弾性部分である定常ひずみに漸近的に近づきます。標準線形ソリッド モデルは、材料の応答を予測する点でマクスウェル モデルやケルビン-フォークト モデルよりも正確ですが、数学的には、特定の荷重条件下でのひずみに対して不正確な結果を返します。
ジェフリーズモデル
ジェフリーズモデルはツェナーモデルと同様に3要素モデルであり、2つのダッシュポットと1つのスプリングで構成されています。[9]

1929年にハロルド・ジェフリーズによって地球のマントルを研究するために提案されました。[10]
バーガーモデル
バーガース モデルは、並列の 2 つのマクスウェル コンポーネント、または直列のケルビン-フォークト コンポーネント、スプリング、ダッシュポットで構成されます。このモデルでは、支配的な構成関係は次のとおりです。
このモデルは、標準的な線形ソリッド モデルに粘性流動を組み込み、固定荷重条件下でのひずみに対して線形に増加する漸近線を提供します。
一般化マクスウェルモデル

一般化マクスウェルモデルは、ヴィーヒャートモデルとも呼ばれ、粘弾性の線形モデルの最も一般的な形式です。このモデルでは、緩和は単一の時点で起こるのではなく、時間の分布で起こることを考慮しています。分子セグメントの長さが異なり、短い分子セグメントの寄与が長い分子セグメントの寄与よりも小さいため、時間の分布は変化します。ヴィーヒャートモデルでは、分布を正確に表すために必要な数のスプリングダッシュポットマクスウェル要素を使用してこれを示します。右の図は、一般化ヴィーヒャートモデルを示しています。[11] 用途: 絶対融点の 4 分の 1 未満の温度での金属および合金 (K で表示)。
非線形粘弾性の構成モデル
非線形粘弾性構成方程式は、法線応力の差、せん断減粘、伸長増粘などの流体の現象を定量的に説明するために必要です。[3] 必然的に、材料が経験した履歴は時間依存の挙動を説明するために必要であり、通常は履歴カーネルKとしてモデルに組み込まれます。[12]
二次流体
2次流体は、通常、最も単純な非線形粘弾性モデルと考えられており、ニュートン流体と他のより複雑な非線形粘弾性流体の間で、高ひずみ振幅とデボラ数で発生する材料挙動の狭い領域で発生します。[3] 2次流体の構成方程式は次のように表されます。
どこ:
- は恒等テンソルである
- は変形テンソルである
- はそれぞれ粘度、第1および第2の法線応力係数を表す。
- は変形テンソルの上部対流微分を表し、は変形テンソルの物質時間微分である。[3]
上対流マクスウェルモデル
上対流型マクスウェルモデルは粘弾性マクスウェルモデルに非線形時間挙動を取り入れたもので、次のように表される。[3]
ここで は応力テンソルを表します。
オールドロイドBモデル
オールドロイドBモデルは上部対流マクスウェルモデルの拡張であり、弾性ビーズとスプリングダンベルで満たされた溶媒として解釈されます。このモデルは、作成者であるジェームズ・G・オールドロイドにちなんで名付けられました。[13] [14] [15]
モデルは 次のように記述できます:
- は応力 テンソルです。
- リラックスタイムです。
- 遅延時間 = ;
- 応力テンソルの上対流時間微分です。
- 流体の速度です。
- は溶媒とポリマー成分からなる総粘度です。
- は変形速度テンソルまたはひずみ速度テンソルです。
このモデルはせん断流における粘弾性流体の近似値として優れているが、伸張流においてはダンベルが無限に引き伸ばされる非物理的な特異点がある。しかし、これは理想化された流れに特有のものであり、クロススロット形状の場合、伸張流は理想的ではないため、応力は特異ではあるが積分可能であり、それに応じて無限に小さい領域では応力は無限大となる。[15]
溶媒粘度がゼロの場合、Oldroyd-B は上部対流マクスウェルモデルになります。
ワーグナーモデル
Wagner モデルは、Bernstein–Kearsley–Zapas モデルの簡略化された実用的な形式であると考えられます。このモデルは、ドイツのレオロジストManfred Wagnerによって開発されました。
等温条件の場合、モデルは次のように記述できます。
どこ:
- は時間tの関数としてのコーシー応力テンソルであり、
- pは圧力である
- は統一テンソルである
- Mはメモリ関数であり、通常は各緩和モードの指数項の合計として表されます。ここで、各緩和モードについて、は緩和係数、は緩和時間です。
- は、フィンガーテンソルの第一および第二不変量に依存するひずみ減衰関数です。
ひずみ減衰関数は通常次のように表される: ひずみ硬化関数の値が1に等しい場合、変形は小さい。0に近づくと、変形は大きくなる。[16] [17]
プロニーシリーズ
1 次元の緩和試験では、材料は試験期間中一定に保たれる突然のひずみを受け、時間の経過とともに応力が測定されます。初期応力は材料の弾性応答によるものです。その後、材料の粘性効果により、時間の経過とともに応力が緩和されます。通常、引張、圧縮、バルク圧縮、またはせん断ひずみのいずれかが適用されます。結果として得られる応力対時間データは、モデルと呼ばれるいくつかの方程式に適合できます。適用されるひずみの種類に応じて表記法のみが変わります。引張圧縮緩和は、せん断は、バルクは と表示されます。せん断緩和の Prony 級数は、
ここで、は材料が完全に緩和された後の長期弾性係数、は緩和時間(図中の と混同しないでください)です。値が高いほど、応力が緩和するのに時間がかかります。データは、予測値とデータ値の間の誤差を最小化するためにパラメータ( )を調整する最小化アルゴリズムを使用して式に適合されます。 [18]
弾性係数が長期係数と次のように関係していることに注目すると、別の形式が得られる。
したがって、
この形式は、弾性せん断弾性率が緩和データとは独立したデータから得られる場合や、Simulia (2010)のように弾性特性を粘性特性とは別に指定するコンピュータ実装の場合に便利です。[19]
クリープ実験は通常、緩和実験よりも簡単に実行できるため、ほとんどのデータは(クリープ)コンプライアンス対時間として入手できます。[20] 残念ながら、プロニー級数の係数に関して、(クリープ)コンプライアンスの閉じた形式は知られていません。そのため、クリープデータがある場合、たとえばで必要な(緩和)プロニー級数の係数を取得するのは簡単ではありません。[19]これらの係数を取得する便利な方法は次のとおりです。まず、コンプライアンスと緩和の両方で閉じた形式のソリューションを持つモデル、たとえばBarbero(2007)[21]のMaxwell-Kelvinモデル(式7.18-7.19)またはBarbero(2007) [21]の標準固体モデル(式7.20-7.21)(セクション7.1.3)でクリープデータをフィッティングします。クリープ モデルのパラメータがわかったら、共役緩和モデルを使用して、元のデータと同じ時間の緩和疑似データを作成します。最後に、疑似データを Prony 級数に適合させます。
温度の影響
ポリマーの二次結合は、熱運動により絶えず切断され、再形成されます。応力を加えると、他の構造よりもいくつかの構造が有利になるため、ポリマーの分子は時間の経過とともに徐々に有利な構造に「流れ」ます。[22] 熱運動はポリマーの変形に寄与する要因の 1 つであるため、粘弾性特性は温度の上昇または低下とともに変化します。ほとんどの場合、適用された応力と時間依存のひずみの比率として定義されるクリープ係数は、温度の上昇とともに減少します。一般的に、温度の上昇は、一定の応力下で等しいひずみを与えるのに必要な時間の対数的な減少と相関しています。言い換えると、粘弾性材料を同じ距離だけ伸ばすのに必要な作業は、温度が高い場合の方が低い場合よりも少なくなります。
ポリマーの粘弾性挙動に対する温度の影響をより詳細にプロットするには、図を参照してください。
典型的なポリマーには主に5つの領域(IVとVを組み合わせた4つと表記されるものもある)が含まれる。[23]
- 領域I:この領域ではポリマーのガラス状態が表れる。この領域ではポリマーの温度が低すぎて分子運動が起こらない。そのため、この領域では分子運動が凍結している。この領域では機械的性質は硬く脆い。[24]
- 領域 II: この領域ではポリマーはガラス転移温度を超えます。Tg を超えると、環境から供給される熱エネルギーは分子の動きを解凍するのに十分です。この領域では分子が局所的に動くことが許されるため、領域 I と比較して剛性が急激に低下します。
- 領域 III: ゴム状プラトー領域。この領域にある材料は、エントロピーによって駆動される長距離弾性を有します。たとえば、この領域の初期状態では、輪ゴムは無秩序です。輪ゴムを伸ばすと、構造も整列してより秩序立ったものになります。したがって、輪ゴムを放すと、輪ゴムは自発的により高いエントロピー状態を求めて、初期状態に戻ります。これが、エントロピー駆動型弾性形状回復と呼ばれるものです。
- 領域IV:ゴム状流動領域における挙動は、時間依存性が高い。この領域のポリマーは、時間と温度の重ね合わせを使用してより詳細な情報を取得し、材料の使用方法を慎重に決定する必要がある。たとえば、材料が短い相互作用時間の目的で使用される場合、それは「硬い」材料として現れる可能性がある。長い相互作用時間の目的で使用する場合、それは「柔らかい」材料として機能する。[25]
- 領域 V: この領域では粘性ポリマーが容易に流れます。剛性がさらに大幅に低下します。

極端に低い温度では、粘弾性材料がガラス相に変化して脆くなることがあります。たとえば、感圧接着剤を極度に低い温度(ドライアイス、凍結スプレーなど)にさらすと、粘着力が失われ、剥離が発生します。
粘弾性クリープ

粘弾性材料は、一定のステップ応力を受けると、時間に依存してひずみが増加します。この現象は粘弾性クリープとして知られています。
時刻 において、粘弾性材料に一定の応力が加わり、この応力は十分に長い時間維持されます。粘弾性液体の場合、材料は応力に応答してひずみが増加し、最終的に破損します。一方、粘弾性固体の場合は、加えられた応力と材料の最終的な抵抗力に応じて破損する場合と破損しない場合があります。応力が短時間維持されると、材料は時刻 まで初期ひずみを受け、その後、ひずみは直ちに減少し (不連続)、その後徐々に減少して残留ひずみになります。
粘弾性クリープデータは、クリープ係数(一定の印加応力を特定の時間における全ひずみで割った値)を時間の関数としてプロットすることで表すことができます。[26] 臨界応力以下では、粘弾性クリープ係数は印加応力とは無関係です。印加応力が材料の臨界応力値を下回る場合、さまざまな印加応力に対するひずみ対時間応答を表す一連の曲線は、単一の粘弾性クリープ係数対時間曲線で表すことができます。
粘弾性クリープは、長期的な構造設計を考える上で重要です。荷重と温度条件を考慮すると、設計者はコンポーネントの寿命に最適な材料を選択できます。
測定
せん断レオメトリー
せん断レオメーターは、測定対象の材料を2枚のプレートの間に置き、片方または両方をせん断方向に動かして材料に応力とひずみを生じさせるというアイデアに基づいています。試験は、一定のひずみ速度、応力、または振動方式(動的機械分析の一形態)で行うことができます。[27] せん断レオメーターは通常、材料が2枚のプレートの間から漏れ出し、材料/プレートの界面で滑るエッジ効果によって制限されます。
伸長レオメトリー
伸張レオメーター(エクステンシオメーターとも呼ばれる)は、粘弾性流体を通常一軸方向に引っ張ることで粘弾性特性を測定します。[28] これは通常毛細管力を利用し、流体を狭い形状に制限するため、この技術は希薄ポリマー溶液や一部の溶融ポリマーなど、比較的粘度の低い流体に限定されることが多いです。[28] 伸張レオメーターは、伸張計の端でのエッジ効果と毛細管の内側と外側の圧力差によっても制限されます。[3]
前述の明らかな制限にもかかわらず、高粘度の流体に対しても伸長レオメトリーを実施できます。これにはさまざまな機器を使用する必要がありますが、これらの技術と装置により、ポリマー溶融物などの材料の伸長粘弾性特性を研究できます。過去 50 年間に開発された最も一般的な伸長レオメトリー機器には、マイスナー型レオメーター、フィラメント伸長レオメーター (FiSER)、およびセントマナット伸長レオメーター (SER) の 3 つがあります。
マイスナーとホステットラーが1996年に開発したマイスナー型レオメーターは、2組の逆回転ローラーを使用してサンプルを一軸方向に歪ませます。[29]この方法では、実験全体を通じて一定のサンプル長を使用し、サンプルをエアクッションでローラー間で支えてサンプルのたるみの影響を排除します。ただし、いくつかの問題があります。1つは、ベルトで流体が滑る可能性があり、予想よりも歪率が低くなることです。さらに、この機器は操作が難しく、購入と維持にコストがかかります。
FiSERレオメーターは、2枚のプレートの間に流体を入れるだけです。実験中、上側のプレートは固定され、下側のプレートに力が加えられ、上側のプレートから離れる方向に動きます。[30]ひずみ速度は、サンプルの中央での半径の変化率によって測定されます。次の式を使用して計算されます。 ここで、は中央半径の値、はひずみ速度です。次に、サンプルの粘度は次の式を使用して計算されます。 ここで、はサンプルの粘度、はサンプルを引き離すために適用される力です。
マイスナー型レオメーターと同様に、SER レオメーターは 2 本のローラーを使用して、サンプルを所定の速度で歪ませます。[31]次に、よく知られている式を使用してサンプルの粘度を計算します。 ここで、は応力、は粘度、は歪速度です。この場合の応力は、機器に搭載されているトルク トランスデューサーによって決定されます。この機器は小型であるため使いやすく、ローラー間のサンプルのたるみがありません。SER 伸長レオメーターの動作を詳しく説明した図は、右側にあります。

その他の方法
材料の機械的および粘弾性応答を試験する機器は多数ありますが、広帯域粘弾性分光法(BVS) と共鳴超音波分光法(RUS) は、周囲温度以上でも以下でも使用でき、粘弾性の試験に特化しているため、粘弾性挙動の試験によく使用されます。これら 2 つの機器は、時間と温度の重ね合わせに頼ることなく、さまざまな周波数と時間範囲で減衰メカニズムを採用しています。BVS と RUS を使用して材料の機械的特性を調べることは、粘弾性を示す材料がどのように機能するかを理解するために重要です。[32]
参照
参考文献
- ^ abcde MeyersとChawla(1999):「材料の機械的挙動」、98-103。
- ^ abc McCrum、Buckley、およびBucknell(2003):「ポリマーエンジニアリングの原理」、117-176。
- ^ abcdefg Macosko, Christopher W. (1994). レオロジー:原理、測定、および応用。ニューヨーク:VCH。ISBN 978-1-60119-575-3. OCLC 232602530.
- ^ ビスワス、アビジット;マニバンナン、M.スリニバサン、マンディヤム A. (2015)。 「パチニ小体のマルチスケール層状生体力学的モデル」。ハプティクスに関するIEEEトランザクション。8 (1): 31-42。土井:10.1109/TOH.2014.2369416。PMID 25398182。S2CID 24658742 。
- ^ ab Van Vliet、Krystyn J. (2006)。「3.032 材料の機械的挙動」
- ^ abc Cacopardo, Ludovica (2019年1月). 「ハイドロゲルの粘弾性工学」. Journal of the Mechanical Behavior of Biomedical Materials . 89 : 162–167. doi :10.1016/j.jmbbm.2018.09.031. hdl : 11568/930491 . PMID 30286375. S2CID 52918639 – Elsevier経由. 引用エラー: 名前付き参照 ":0" が異なる内容で複数回定義されています (ヘルプ ページを参照してください)。
- ^ ラーソン、ロナルド G. (1999 年 1 月 28 日)。複雑な流体の構造とレオロジー (化学工学のトピック): ラーソン、ロナルド G.: 9780195121971: Amazon.com: 書籍。Oup USA。ISBN 019512197X。
- ^ Tanner, Roger I. (1988).エンジニアリングレオロジー. Oxford University Press. p. 27. ISBN 0-19-856197-0。
- ^ バーンズ、ハワード A.; ハットン、ジョン フレッチャー; ウォルターズ、K. (1989)。レオロジー入門。エルゼビア。ISBN 978-0-444-87140-4。
- ^バード、R. バイロン ( 1987-05-27 )。高分子液体のダイナミクス、第1巻:流体力学。Wiley。ISBN 978-0-471-80245-7。
- ^ ロイランス、デイビッド(2001)「粘弾性工学」14-15
- ^ Drapaca, CS; Sivaloganathan, S.; Tenti, G. (2007-10-01). 「粘弾性における非線形構成法則」.固体の数学と力学. 12 (5): 475–501. doi :10.1177/1081286506062450. ISSN 1081-2865. S2CID 121260529.
- ^ Oldroyd, James (1950年2 月)。「レオロジー状態方程式の定式化について」。Proceedings of the Royal Society of London。シリーズ A、数学および物理科学。200 (1063): 523–541。Bibcode : 1950RSPSA.200..523O。doi :10.1098/rspa.1950.0035。S2CID 123239889 。
- ^ Owens, RG; Phillips, TN (2002).計算レオロジー. インペリアル・カレッジ・プレス. ISBN 978-1-86094-186-3。
- ^ ab Poole, Rob (2007年10月). 「純粋に弾性的な流れの非対称性」. Physical Review Letters . 99 (16): 164503. Bibcode :2007PhRvL..99p4503P. doi :10.1103/PhysRevLett.99.164503. hdl : 10400.6/634 . PMID 17995258.
- ^ Wagner, Manfred (1976). 「低密度分岐ポリエチレン溶融物のせん断および伸長流動に関する時間依存非線形応力成長データの分析」Rheologica Acta . 15 (2): 136–142. Bibcode :1976AcRhe..15..136W. doi :10.1007/BF01517505. S2CID 96165087.
- ^ Wagner, Manfred (1977). 「単一積分構成方程式を用いたせん断粘度データからの一次法線応力差の予測」Rheologica Acta . 16 (1977): 43–50. Bibcode :1977AcRhe..16...43W. doi :10.1007/BF01516928. S2CID 98599256.
- ^ EJ Barbero. 「線形粘弾性材料の長期応答を予測するための時間・温度・経年変化の重ね合わせ原理」、ポリマーマトリックス複合材料のクリープと疲労の第 2 章。Woodhead、2011 年。
- ^ ab Simulia. Abaqus 解析ユーザーズマニュアル、19.7.1「時間領域粘弾性」、6.10 版、2010 年
- ^ 統一基準によるコンピュータ支援による素材の事前選択
- ^ ab EJ Barbero. 複合材料の有限要素解析。CRC Press、フロリダ州ボカラトン、2007年。
- ^ SA Baeurle、A. Hotta、AA Gusev、ポリマー47、6243-6253 (2006)。
- ^ Aklonis., JJ (1981). 「ポリマーの機械的性質」J Chem Educ . 58 (11): 892. Bibcode :1981JChEd..58..892A. doi :10.1021/ed058p892.
- ^ IM、Kalogeras (2012)。「ガラス状態の性質:構造とガラス転移」。Journal of Materials Education。34 ( 3): 69 。
- ^ I, Emri (2010).固体ポリマーの時間依存挙動。
- ^ Rosato他(2001):「プラスチック設計ハンドブック」、63-64ページ。
- ^ Magnin, A.; Piau, JM (1987-01-01). 「降伏応力を伴う流体のせん断レオメトリー」.非ニュートン流体力学ジャーナル. 23 : 91–106. Bibcode :1987JNNFM..23...91M. doi :10.1016/0377-0257(87)80012-5. ISSN 0377-0257.
- ^ ab Dealy, JM (1978-01-01). 「溶融ポリマー用伸長レオメーター; レビュー」.非ニュートン流体力学ジャーナル. 4 (1–2): 9–21. Bibcode :1978JNNFM...4....9D. doi :10.1016/0377-0257(78)85003-4. ISSN 0377-0257.
- ^ Meissner, J.; Hostettler, J. (1994-01-01). 「ポリマー溶融物および他の高粘弾性液体用の新しい伸長レオメーター」Rheologica Acta . 33 (1): 1–21. Bibcode :1994AcRhe..33....1M. doi :10.1007/BF00453459. ISSN 1435-1528. S2CID 93395453.
- ^ Bach, Anders ; Rasmussen, Henrik Koblitz; Hassager, Ole (2003 年 3 月)。「フィラメント伸張レオメーターで測定されたポリマー溶融物の伸張粘度」。Journal of Rheology。47 ( 2): 429–441。Bibcode : 2003JRheo..47..429B。doi : 10.1122 /1.1545072。ISSN 0148-6055。S2CID 44889615 。
- ^ Sentmanat, Martin L. (2004-12-01). 「小型ユニバーサル試験プラットフォーム: 伸長溶融レオロジーから固体変形挙動まで」Rheologica Acta . 43 (6): 657–669. Bibcode :2004AcRhe..43..657S. doi :10.1007/s00397-004-0405-4. ISSN 1435-1528. S2CID 73671672.
- ^ ロッド・レイクス(1998年)。粘弾性固体。CRCプレス。ISBN 0-8493-9658-1。
- Silbey および Alberty (2001):物理化学、857 ページ。John Wiley & Sons、Inc.
- Alan S. Wineman および KR Rajagopal (2000):ポリマーの機械的応答: 入門
- アレンとトーマス(1999):『材料の構造』、51。
- クランダル他(1999):固体力学入門348
- J. Lemaitre と JL Chaboche (1994)固体材料の力学
- Yu. Dimitrienko (2011)非線形連続体力学と大規模非弾性変形、Springer、772p
