物理法則:バネを変形させるのに必要な力は距離に比例する
フックの法則:力は伸びに比例する
ブルドン管はフックの法則に基づいています。上部のコイル状の金属管内のガス 圧 によって生じる力により、 圧力に比例した量だけ管がほどけます。
多くの機械式時計や腕時計の中核をなすテンプは フック の法則に基づいています。コイルばねによって生成されるトルクはテンプの回転角度に比例するため、その振動周期はほぼ一定です。
物理学 において 、 フックの法則は、 バネ をある距離 ( x )だけ 伸ばすか縮めるのに必要な 力 ( F ) はその距離に対して 直線的に 増減するという 経験法則 である 。つまり、 F s = kx であり 、 k は バネに特有の定数 ( つまりバネの 剛性 ) であり、 x は バネの可能な総変形量に比べて小さい。この法則は、17 世紀のイギリスの物理学者 ロバート・ フックにちなんで名付けられた。彼は 1676 年にこの法則をラテン語の アナグラム として初めて述べた。 [1] [2] 彼は 1678 年にこのアナグラムの解法を 次のように発表した 。 [3] ut tensio, sic vis (「伸びた分だけ力も大きくなる」または「伸びた分は力に比例する」)。フックは 1678 年の著作で、1660 年からこの法則に気づいていたと述べている。
フックの方程式は、高層ビルに吹く風や、ギターの 弦を弾く音楽家など、 弾性 体が 変形する 他の多くの状況でも(ある程度)成り立ちます。この方程式が仮定できる弾性体または材料は、 線形弾性 または フック型 と呼ばれます 。
フックの法則は、バネやその他の弾性体が加えられた力に対して実際に示す反応の 一次線形近似にすぎ ません。力が何らかの限界を超えると、この法則は最終的に破綻します。なぜなら、永久変形や状態変化を起こさずに、ある最小サイズを超えて圧縮したり、最大サイズを超えて引き伸ばしたりすることはできないからです。多くの材料は、 弾性限界 に達するずっと前から、フックの法則から著しく逸脱します。
一方、フックの法則は、力と変形が十分に小さい限り、ほとんどの固体に対して正確な近似値となります。このため、フックの法則は科学と工学のすべての分野で広く使用されており、 地震学 、 分子力学 、 音響学など多くの分野の基礎となっています。また、 バネ秤 、 圧力計 、 検流計 、 機械式時計 の テンプ の基本原理でもあります 。
現代の 弾性理論で は、フックの法則を一般化して、 弾性物体または材料の ひずみ(変形) は、それに加えられた 応力 に比例すると述べています。ただし、一般的な応力とひずみには複数の独立した要素があるため、「比例係数」は単一の実数ではなく、 実数の
行列 で表すことができる 線形マップ ( テンソル) になる可能性があります。
この一般的な形式では、フックの法則により、複雑な物体のひずみと応力の関係を、その物体を構成する材料の固有の特性の観点から推測することが可能になります。たとえば、均一な 断面 を持つ 均質な 棒は、伸ばされると単純なバネのように動作し、剛性 k は断面積に正比例し、長さに反比例することが推測できます。
リニアスプリング
バネの伸びと圧縮
単純な らせんばね を考えてみましょう。このばねの片端は固定された物体に取り付けられており、自由端は F s の大きさの力で引っ張られています。ばねは平衡 状態に達し 、長さはそれ以上変化していないものとします。 ばねの自由端が「リラックス」した位置(伸びていないとき)から変位した量を xとします。フックの法則では、
F s = k x {\displaystyle F_{s}=kx}
または、同等の式が成り立ちます。
ここで、 k は ばねの特性を表す正の実数です。コイル間に隙間があるばねは圧縮することができ、圧縮についても同じ式が成り立ちます 。その場合、 F s と x は両方とも負になります。 [4]
x
=
F
s
k
{\displaystyle x={\frac {F_{s}}{k}}}
グラフ導出
この式によれば、変位 x の関数として 作用する力 F s のグラフは、 原点 を通る直線となり 、その 傾きは k に なります 。
バネに関するフックの法則は、バネの自由端を引っ張る物体にバネが及ぼす 復元力 が F s で あるという慣例に基づいても述べられます。その場合、復元力の方向は変位の方向と反対なので、式は 次のようになります。
F
s
=
−
k
x
{\displaystyle F_{s}=-kx}
ねじりバネ
フックの法則の ねじり 類似法則は、 ねじりバネ に適用されます。これは、物体を回転させるために必要なトルク (τ) が、平衡位置からの角度変位 (θ) に正比例することを述べています。これは、物体に加えられるトルクと、ねじりによる結果として生じる角度 変形 との関係を説明しています。数学的には、次のように表現できます。
τ
=
−
k
θ
{\displaystyle \tau =-k\theta }
どこ:
τ はニュートンメートルまたは N·m で測定される トルク です。
k は ねじり定数 (N·m/ラジアンで測定)であり、ねじりバネの剛性または角度変位に対する抵抗を特徴付けます。
θ は平衡位置からの 角度変位 (ラジアンで測定)です。
線形の場合と同様に、この法則はトルクが角変位に比例することを示しており、負の符号はトルクが角変位と反対方向に作用し、システムを平衡状態に戻す復元力を提供することを示しています。
一般的な「スカラー」スプリング
フックのバネの法則は、変形と応力の両方が正と負の両方をとることができる単一の数値で表現できる限り、任意の複雑さを持つあらゆる弾性物体に通常適用されます。
たとえば、2 枚の平行なプレートに取り付けられたゴムのブロックが、伸張や圧縮ではなく せん断 によって変形する場合、せん断力 F s とプレートの横方向の変位 x は フックの法則に従います (変形が十分に小さい場合)。
フックの法則は、両端で支えられた真っ直ぐな鉄筋またはコンクリートの梁(建物で使用されているものなど)が、中間点に置かれた重り F によって曲げられる場合にも適用されます。この場合の変位 x は 、梁の横方向の偏差であり、荷重がかかっていない形状に対して測定されます。
軸 に沿って伸張または圧縮されるらせんばねの場合 、加えられた(または復元する)力と結果として生じる伸張または圧縮は同じ方向(前述の軸の方向)になります。したがって、 F s と x が ベクトル として定義されている場合でも 、フックの 方程式は有効であり、力のベクトルは 伸張ベクトル に 固定スカラーを乗じたもの であることがわかります 。
弾性体の中には、方向の異なる力を受けると、一方向に変形するものがあります。一例として、垂直でも水平でもない横方向の荷重によって曲がる、正方形でない長方形の断面を持つ水平の木製梁が挙げられます。このような場合、 変位 xの 大きさは 、力 F s の方向が同じである限り(そしてその値が大きすぎない限り)、力 F s の大きさに比例します。したがって、フックの法則のスカラーバージョン F s = − kx が成り立ちます。ただし、力のベクトルと変位の ベクトルは 方向が異なるため、互いのスカラー倍数にはなりません。さらに、それらの大きさの比 k は、ベクトル F s の方向によって異なります 。
しかし、このような場合、 力と変形ベクトルが十分に小さい限り、それらの間には固定された 線形関係 が存在することがよくあります。つまり、ベクトルからベクトルへの 関数 κ が存在し、任意の実数 α 、 β および任意の変位ベクトル X 1 、 X 2に対して、 F = κ ( X ) 、および κ ( α X 1 + β X 2 ) = α κ ( X 1 ) + β κ ( X 2 ) が成り立ちます。このような関数は (2 次) テンソル と呼ばれます。
任意の直交座標系 に関して 、力と変位ベクトルは、実数の 3 × 1 行列 で表すことができます。次に、それらを結ぶテンソル κ は、実係数の 3 × 3 行列 κ で表すことができ 、変位ベクトルを 乗じる と、力のベクトルが得られます。
つまり、 i = 1, 2, 3 の場合です 。したがって、フックの法則 F = κ X は、 X と F が 可変方向のベクトルである場合にも成立すると言えます 。ただし、オブジェクトの剛性は、 単一の実数 kではなく、テンソル κ です。
F
=
[
F
1
F
2
F
3
]
=
[
κ
11
κ
12
κ
13
κ
21
κ
22
κ
23
κ
31
κ
32
κ
33
]
[
X
1
X
2
X
3
]
=
κ
X
{\displaystyle \mathbf {F} \,=\,{\begin{bmatrix}F_{1}\\F_{2}\\F_{3}\end{bmatrix}}\,=\,{\begin{bmatrix}\kappa _{11}&\kappa _{12}&\kappa _{13}\\\kappa _{21}&\kappa _{22}&\kappa _{23}\\\kappa _{31}&\kappa _{32}&\kappa _{33}\end{bmatrix}}{\begin{bmatrix}X_{1}\\X_{2}\\X_{3}\end{bmatrix}}\,=\,{\boldsymbol {\kappa }}\mathbf {X} }
F
i
=
κ
i
1
X
1
+
κ
i
2
X
2
+
κ
i
3
X
3
{\displaystyle F_{i}=\kappa _{i1}X_{1}+\kappa _{i2}X_{2}+\kappa _{i3}X_{3}}
(a) ポリマーナノスプリングの概略図。コイル半径 R、ピッチ P、スプリングの長さ L、巻き数 N は、それぞれ 2.5 μm、2.0 μm、13 μm、4 です。荷重前 (be)、伸長 (f)、圧縮 (g)、曲げ (h)、回復 (i) のナノスプリングの電子顕微鏡写真。スケールバーはすべて 2 μm です。スプリングは加えられた力に対して線形応答を示し、ナノスケールでのフックの法則の妥当性を実証しています。 [5]
連続した 弾性材料(ゴムのブロック、 ボイラー の壁、鉄の棒など)内部の材料の応力とひずみは 、数学的にはフックのバネの法則に似た線形関係で結びついており、その名前で呼ばれることが多いです。
しかし、ある点の周りの固体媒体のひずみ状態は、単一のベクトルでは説明できません。同じ物質の塊は、それがどんなに小さくても、異なる方向に沿って同時に圧縮、伸張、せん断される可能性があります。同様に、その塊の応力は、同時に押す、引っ張る、せん断することもあります。
この複雑さを捉えるためには、点の周りの媒体の関連状態を、2 つの 2 次テンソル、つまりひずみ テンソル ε (変位 X の代わり ) と 応力テンソル σ (復元力 F の代わり) で表す必要があります。連続媒体に対するフックのバネの法則の類似物は、次のようになります。ここで、 c は 4 次テンソル (つまり、2 次テンソル間の線形写像) で、通常は剛性テンソルまたは弾性テンソルと呼ばれます 。次のように書くこともできます。ここで 、 テンソル s は コンプライアンス テンソル と 呼ばれ 、前述の線形写像の逆を表します。
σ
=
c
ε
,
{\displaystyle {\boldsymbol {\sigma }}=\mathbf {c} {\boldsymbol {\varepsilon }},}
ε
=
s
σ
,
{\displaystyle {\boldsymbol {\varepsilon }}=\mathbf {s} {\boldsymbol {\sigma }},}
デカルト座標系では、応力テンソルとひずみテンソルは 3 × 3 行列で表すことができます。9
つの数値 σ ij と 9 つの数値 ε kl の 間の線形マッピングであるため、剛性テンソル c は3 × 3 × 3 × 3 = 81 個 の実数 c ijkl の行列で表されます 。フックの法則によれば、
i 、 j = 1,2,3 です。
ε
=
[
ε
11
ε
12
ε
13
ε
21
ε
22
ε
23
ε
31
ε
32
ε
33
]
;
σ
=
[
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
]
{\displaystyle {\boldsymbol {\varepsilon }}\,=\,{\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{31}&\varepsilon _{32}&\varepsilon _{33}\end{bmatrix}}\,;\qquad {\boldsymbol {\sigma }}\,=\,{\begin{bmatrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\end{bmatrix}}}
σ
i
j
=
∑
k
=
1
3
∑
l
=
1
3
c
i
j
k
l
ε
k
l
{\displaystyle \sigma _{ij}=\sum _{k=1}^{3}\sum _{l=1}^{3}c_{ijkl}\varepsilon _{kl}}
これら 3 つのテンソルは、一般に媒体内の点ごとに変化し、時間とともに変化することもあります。ひずみテンソル ε は 、点の近傍における媒体粒子の変位を指定するだけですが、応力テンソル σ は 、媒体の隣接する小片が互いに及ぼす力を指定します。したがって、これらは材料の組成や物理的状態とは無関係です。 一方、剛性テンソル c は、材料の特性であり、温度、圧力 、 微細構造 などの物理的状態変数に依存することがよくあります。
σ 、 ε 、 c の固有の対称性により 、後者の弾性係数のうち独立したものは 21 個のみです。 [6] この数は、材料の対称性によってさらに減少します。 斜方 晶の場合は 9 個、 六方 晶の場合は 5 個、 立方 晶の場合は 3 個です。 [7] 等方性 媒体(どの方向でも同じ物理的特性を持つ) の場合、 c は 、体積弾性率 K と せん断弾性率 G という2 つの独立した数値にのみ減少します 。これらは、それぞれ体積変化とせん断変形に対する材料の抵抗を定量化します。
類似の法律
フックの法則は2つの量の間の単純な比例関係であるため、その公式と結果は、流体 の運動や 電界 による 誘電体 の 分極を 記述する法則など、他の多くの物理法則と数学的に類似しています 。
特に、弾性応力とひずみを関連付けるテンソル方程式 σ = cε は 、粘性 流体 の流れにおける 粘性応力テンソル τ と ひずみ速度テンソル ε̇ を関連付ける 方程式 τ = με̇ と完全に類似しています。ただし、前者は静的 応力 (変形 量 に関連) に関連し、後者は 動的 応力 ( 変形
速度に関連) に関連します。
測定単位
SI 単位 では 、変位はメートル (m) で測定され、力は ニュートン (N または kg·m/s 2 ) で測定されます。したがって、バネ定数 k とテンソル κ の各要素は、ニュートン毎メートル (N/m) またはキログラム毎秒の 2 乗 (kg/s 2 ) で測定されます。
連続媒体の場合、応力テンソル σ の各要素は力を面積で割った値です。したがって、圧力の単位、つまり パスカル (Pa、N/m 2 、kg/(m·s 2 )) で測定されます。ひずみテンソル ε の要素は無次元 です(変位を距離で割った値)。したがって、 c ijkl の要素 も圧力の単位で表されます。
弾性材料への一般的な応用
低炭素鋼の応力-ひずみ曲線。 応力 (単位面積あたりの力)と ひずみ (結果として生じる圧縮/伸張、変形として知られる)の関係を示します 。フックの法則は、曲線の原点と降伏点(2)の間の部分に対してのみ有効です。 見かけの応力( F / A 0 ) 実応力( F / A )
力によって変形された後、その物質の分子または原子が安定した平衡の初期状態に戻り、すぐに元の形状に戻る物体は、多くの場合、フックの法則に従います。
フックの法則は、特定の荷重条件下では一部の材料にのみ適用されます。ほとんどの工学用途では、鋼鉄は線形弾性挙動を示します。フックの法則は、鋼鉄の 弾性範囲 全体にわたって有効です(つまり、 降伏強度 以下の応力に対して)。アルミニウムなどの他の材料の場合、フックの法則は弾性範囲の一部に対してのみ有効です。これらの材料には 比例限界 応力が定義されており、それ以下では線形近似に関連する誤差は無視できます。
ゴムは、その弾性が応力に依存し、温度と荷重率に敏感であるため、一般に「非フック」材料と見なされます。
大変形 の場合のフックの法則の一般化は、 新フック固体 と ムーニー・リブリン固体 のモデルによって提供されます 。
弾性 材料の棒は 線形 バネ として考えることができます。棒の長さは L 、断面積は A です。その 引張応力 σ は、 弾性係数 E によって、その分数伸びまたはひずみ ε に比例します 。
σ
=
E
ε
.
{\displaystyle \sigma =E\varepsilon .}
弾性率は、多くの場合、一定であると考えられます。つまり、
(つまり、長さの変化率)であり、
次のようになります。
ε
=
Δ
L
L
{\displaystyle \varepsilon ={\frac {\Delta L}{L}}}
σ
=
F
A
,
{\displaystyle \sigma ={\frac {F}{A}}\,,}
ε
=
σ
E
=
F
A
E
.
{\displaystyle \varepsilon ={\frac {\sigma }{E}}={\frac {F}{AE}}\,.}
長さの変化は次のように表される。
Δ
L
=
ε
L
=
F
L
A
E
.
{\displaystyle \Delta L=\varepsilon L={\frac {FL}{AE}}\,.}
春のエネルギー
バネに蓄えられた 位置エネルギー U el ( x ) は、バネを徐々に圧縮するのに必要なエネルギーを加算することで得られます 。つまり、力と変位の積分です。外力は変位と同じ方向を向いているため、バネの位置エネルギーは常に非負です。代入すると、
U
e
l
(
x
)
=
1
2
k
x
2
{\displaystyle U_{\mathrm {el} }(x)={\tfrac {1}{2}}kx^{2}}
x
=
F
/
k
{\displaystyle x=F/k}
U
e
l
(
F
)
=
F
2
2
k
.
{\displaystyle U_{\mathrm {el} }(F)={\frac {F^{2}}{2k}}.}
このポテンシャル U el は、 Ux 平面上の 放物線 として視覚化することができ 、 U el ( x ) = 1 / 2 kx 2 。バネが正の x 方向に伸びると、位置エネルギーは放物線状に増加します (バネが圧縮されるときも同じことが起こります)。位置エネルギーの変化は一定の割合で変化するため、
変位と加速度がゼロの場合でも
U
の変化は一定であることに注意してください
d
2
U
e
l
d
x
2
=
k
.
{\displaystyle {\frac {d^{2}U_{\mathrm {el} }}{dx^{2}}}=k\,.}
緩和された力の定数(一般化されたコンプライアンス定数)
緩和された力の定数(一般化されたコンプライアンス定数 の逆 )は、通常の「固定された」力の定数とは対照的に、分子系に対して一意に定義されるため、それを使用することで、 化学反応の 反応物 、 遷移状態 、および生成物に対して計算された力場の間で意味のある相関関係を作成できます。 位置エネルギーが 内部座標で2次形式で記述できるのと同様に 、一般化された力で記述することもできます。 結果として得られる係数は コンプライアンス定数 と呼ばれます。 標準モード解析を行う必要がなく、分子の任意の内部座標のコンプライアンス定数を直接計算する方法があります。 [8]緩和された力の定数(逆コンプライアンス定数)が共有 結合 強度記述子として適していることは、 1980年に早くも実証されました。 最近、非共有結合強度記述子としての適性も実証されました。 [9]
調和振動子
バネで吊り下げられた質量は調和振動子の典型的な例である。
バネの端に取り付けられた 質量 m は、 調和振動子 の典型的な例です。質量をわずかに引っ張ってから放すと、システムは平衡位置を中心に 正弦波振動運動を開始します。バネがフックの法則に従い、 摩擦 とバネの質量を無視できる限り 、振動の振幅は一定のままです。その周波数 f は 振幅とは無関係で、バネの質量と剛性のみによって決まります。
この現象により 、船や人々のポケットに入れて持ち運べる
正確な 機械式時計 や腕時計
の製造が可能になりました。
f
=
1
2
π
k
m
{\displaystyle f={\frac {1}{2\pi }}{\sqrt {\frac {k}{m}}}}
無重力空間での回転
質量 m が力定数 k のバネに取り付けられ 、自由空間で回転する場合、バネの張力 ( F t ) によって必要な向心力 ( F c )が供給されます 。
F t = F c かつ x = r で
あるため 、次のようになります。
ω = 2π f
とすると 、上記と同じ周波数方程式が導かれます。
F
t
=
k
x
;
F
c
=
m
ω
2
r
{\displaystyle F_{\mathrm {t} }=kx\,;\qquad F_{\mathrm {c} }=m\omega ^{2}r}
k
=
m
ω
2
{\displaystyle k=m\omega ^{2}}
f
=
1
2
π
k
m
{\displaystyle f={\frac {1}{2\pi }}{\sqrt {\frac {k}{m}}}}
等方性材料
等方性材料は、空間の方向に依存しない特性を特徴とする。従って、等方性材料を含む物理方程式は、それらを表すために選択された座標系に依存しない必要がある。ひずみテンソルは対称テンソルである。テンソルの トレース はどの座標系からも独立しているため、対称テンソルの最も完全な座標フリー分解は、定数テンソルとトレースレス対称テンソルの和として表すことである。 [10] したがって、 インデックス表記 では次の式で表される。
ここで、 δ ij は クロネッカーのデルタ である 。直接テンソル表記では次の式で表される。
ここで、 I は2次単位テンソルである。
ε
i
j
=
(
1
3
ε
k
k
δ
i
j
)
+
(
ε
i
j
−
1
3
ε
k
k
δ
i
j
)
{\displaystyle \varepsilon _{ij}=\left({\tfrac {1}{3}}\varepsilon _{kk}\delta _{ij}\right)+\left(\varepsilon _{ij}-{\tfrac {1}{3}}\varepsilon _{kk}\delta _{ij}\right)}
ε
=
vol
(
ε
)
+
dev
(
ε
)
;
vol
(
ε
)
=
1
3
tr
(
ε
)
I
;
dev
(
ε
)
=
ε
−
vol
(
ε
)
{\displaystyle {\boldsymbol {\varepsilon }}=\operatorname {vol} ({\boldsymbol {\varepsilon }})+\operatorname {dev} ({\boldsymbol {\varepsilon }})\,;\qquad \operatorname {vol} ({\boldsymbol {\varepsilon }})={\tfrac {1}{3}}\operatorname {tr} ({\boldsymbol {\varepsilon }})~\mathbf {I} \,;\qquad \operatorname {dev} ({\boldsymbol {\varepsilon }})={\boldsymbol {\varepsilon }}-\operatorname {vol} ({\boldsymbol {\varepsilon }})}
右側の最初の項は定数テンソル( 体積ひずみテンソル とも呼ばれる)であり、2 番目の項はトレースレス対称テンソル( 偏差ひずみテンソル またはせん断テンソルとも呼ばれる)です。
等方性材料に対するフックの法則の最も一般的な形式は、次の 2 つのテンソルの線形結合として表すことができます。
ここで、 K は 体積弾性係数 、 G は せん断弾性係数 です 。
σ
i
j
=
3
K
(
1
3
ε
k
k
δ
i
j
)
+
2
G
(
ε
i
j
−
1
3
ε
k
k
δ
i
j
)
;
σ
=
3
K
vol
(
ε
)
+
2
G
dev
(
ε
)
{\displaystyle \sigma _{ij}=3K\left({\tfrac {1}{3}}\varepsilon _{kk}\delta _{ij}\right)+2G\left(\varepsilon _{ij}-{\tfrac {1}{3}}\varepsilon _{kk}\delta _{ij}\right)\,;\qquad {\boldsymbol {\sigma }}=3K\operatorname {vol} ({\boldsymbol {\varepsilon }})+2G\operatorname {dev} ({\boldsymbol {\varepsilon }})}
弾性係数 の関係を利用すると 、これらの式はさまざまな方法で表現することもできます。等方性材料に対するフックの法則の一般的な形式は、直接テンソル表記で表現され、
[11]
で、 λ = K −
σ
=
λ
tr
(
ε
)
I
+
2
μ
ε
=
c
:
ε
;
c
=
λ
I
⊗
I
+
2
μ
I
{\displaystyle {\boldsymbol {\sigma }}=\lambda \operatorname {tr} ({\boldsymbol {\varepsilon }})\mathbf {I} +2\mu {\boldsymbol {\varepsilon }}={\mathsf {c}}:{\boldsymbol {\varepsilon }}\,;\qquad {\mathsf {c}}=\lambda \mathbf {I} \otimes \mathbf {I} +2\mu {\mathsf {I}}}
2 / 3 G = c 1111 − 2 c 1212 および μ = G = c 1212 は ラメ定数 、 I は第 2 階の恒等テンソル、 I は 第 4 階の恒等テンソルの対称部分です。 指数表記では、
逆の関係は [12]です。 したがって、関係 ε = s : σ
におけるコンプライアンステンソルは です。
ヤング率 と ポアソン比
を使用すると 、等方性材料に対するフックの法則は次のように表すことができます。
これは、工学においてひずみが応力テンソルで表される形式です。 展開された形式の式は、
ここで E は ヤング率 、 νは ポアソン比 です 。 ( 3-D 弾性 を参照)。
σ
i
j
=
λ
ε
k
k
δ
i
j
+
2
μ
ε
i
j
=
c
i
j
k
l
ε
k
l
;
c
i
j
k
l
=
λ
δ
i
j
δ
k
l
+
μ
(
δ
i
k
δ
j
l
+
δ
i
l
δ
j
k
)
{\displaystyle \sigma _{ij}=\lambda \varepsilon _{kk}~\delta _{ij}+2\mu \varepsilon _{ij}=c_{ijkl}\varepsilon _{kl}\,;\qquad c_{ijkl}=\lambda \delta _{ij}\delta _{kl}+\mu \left(\delta _{ik}\delta _{jl}+\delta _{il}\delta _{jk}\right)}
ε
=
1
2
μ
σ
−
λ
2
μ
(
3
λ
+
2
μ
)
tr
(
σ
)
I
=
1
2
G
σ
+
(
1
9
K
−
1
6
G
)
tr
(
σ
)
I
{\displaystyle {\boldsymbol {\varepsilon }}={\frac {1}{2\mu }}{\boldsymbol {\sigma }}-{\frac {\lambda }{2\mu (3\lambda +2\mu )}}\operatorname {tr} ({\boldsymbol {\sigma }})\mathbf {I} ={\frac {1}{2G}}{\boldsymbol {\sigma }}+\left({\frac {1}{9K}}-{\frac {1}{6G}}\right)\operatorname {tr} ({\boldsymbol {\sigma }})\mathbf {I} }
s
=
−
λ
2
μ
(
3
λ
+
2
μ
)
I
⊗
I
+
1
2
μ
I
=
(
1
9
K
−
1
6
G
)
I
⊗
I
+
1
2
G
I
{\displaystyle {\mathsf {s}}=-{\frac {\lambda }{2\mu (3\lambda +2\mu )}}\mathbf {I} \otimes \mathbf {I} +{\frac {1}{2\mu }}{\mathsf {I}}=\left({\frac {1}{9K}}-{\frac {1}{6G}}\right)\mathbf {I} \otimes \mathbf {I} +{\frac {1}{2G}}{\mathsf {I}}}
ε
i
j
=
1
E
(
σ
i
j
−
ν
(
σ
k
k
δ
i
j
−
σ
i
j
)
)
;
ε
=
1
E
(
σ
−
ν
(
tr
(
σ
)
I
−
σ
)
)
=
1
+
ν
E
σ
−
ν
E
tr
(
σ
)
I
{\displaystyle \varepsilon _{ij}={\frac {1}{E}}{\big (}\sigma _{ij}-\nu (\sigma _{kk}\delta _{ij}-\sigma _{ij}){\big )}\,;\qquad {\boldsymbol {\varepsilon }}={\frac {1}{E}}{\big (}{\boldsymbol {\sigma }}-\nu (\operatorname {tr} ({\boldsymbol {\sigma }})\mathbf {I} -{\boldsymbol {\sigma }}){\big )}={\frac {1+\nu }{E}}{\boldsymbol {\sigma }}-{\frac {\nu }{E}}\operatorname {tr} ({\boldsymbol {\sigma }})\mathbf {I} }
ε
11
=
1
E
(
σ
11
−
ν
(
σ
22
+
σ
33
)
)
ε
22
=
1
E
(
σ
22
−
ν
(
σ
11
+
σ
33
)
)
ε
33
=
1
E
(
σ
33
−
ν
(
σ
11
+
σ
22
)
)
ε
12
=
1
2
G
σ
12
;
ε
13
=
1
2
G
σ
13
;
ε
23
=
1
2
G
σ
23
{\displaystyle {\begin{aligned}\varepsilon _{11}&={\frac {1}{E}}{\big (}\sigma _{11}-\nu (\sigma _{22}+\sigma _{33}){\big )}\\\varepsilon _{22}&={\frac {1}{E}}{\big (}\sigma _{22}-\nu (\sigma _{11}+\sigma _{33}){\big )}\\\varepsilon _{33}&={\frac {1}{E}}{\big (}\sigma _{33}-\nu (\sigma _{11}+\sigma _{22}){\big )}\\\varepsilon _{12}&={\frac {1}{2G}}\sigma _{12}\,;\qquad \varepsilon _{13}={\frac {1}{2G}}\sigma _{13}\,;\qquad \varepsilon _{23}={\frac {1}{2G}}\sigma _{23}\end{aligned}}}
3次元におけるフックの法則の導出
フックの法則の 3 次元形式は、ポアソン比とフックの法則の 1 次元形式を使用して次のように導き出すことができます。ひずみと応力の関係を、荷重方向の伸張 (1) と垂直方向の収縮 (荷重による) (2 と 3) という 2 つの効果の重ね合わせとして考えます。
ここで、 ν はポアソン比、 E は ヤング率です。
ε
1
′
=
1
E
σ
1
,
ε
2
′
=
−
ν
E
σ
1
,
ε
3
′
=
−
ν
E
σ
1
,
{\displaystyle {\begin{aligned}\varepsilon _{1}'&={\frac {1}{E}}\sigma _{1}\,,\\\varepsilon _{2}'&=-{\frac {\nu }{E}}\sigma _{1}\,,\\\varepsilon _{3}'&=-{\frac {\nu }{E}}\sigma _{1}\,,\end{aligned}}}
方向2と3の荷重についても同様の式が得られ
、
ε
1
″
=
−
ν
E
σ
2
,
ε
2
″
=
1
E
σ
2
,
ε
3
″
=
−
ν
E
σ
2
,
{\displaystyle {\begin{aligned}\varepsilon _{1}''&=-{\frac {\nu }{E}}\sigma _{2}\,,\\\varepsilon _{2}''&={\frac {1}{E}}\sigma _{2}\,,\\\varepsilon _{3}''&=-{\frac {\nu }{E}}\sigma _{2}\,,\end{aligned}}}
ε
1
‴
=
−
ν
E
σ
3
,
ε
2
‴
=
−
ν
E
σ
3
,
ε
3
‴
=
1
E
σ
3
.
{\displaystyle {\begin{aligned}\varepsilon _{1}'''&=-{\frac {\nu }{E}}\sigma _{3}\,,\\\varepsilon _{2}'''&=-{\frac {\nu }{E}}\sigma _{3}\,,\\\varepsilon _{3}'''&={\frac {1}{E}}\sigma _{3}\,.\end{aligned}}}
3つの場合を合計すると( ε i = ε i ′ + ε i ″ + ε i ‴ )
、 νσを
1つ加算して減算すると、 σ 1を
解くと次の式が得られます。
ε
1
=
1
E
(
σ
1
−
ν
(
σ
2
+
σ
3
)
)
,
ε
2
=
1
E
(
σ
2
−
ν
(
σ
1
+
σ
3
)
)
,
ε
3
=
1
E
(
σ
3
−
ν
(
σ
1
+
σ
2
)
)
,
{\displaystyle {\begin{aligned}\varepsilon _{1}&={\frac {1}{E}}{\big (}\sigma _{1}-\nu (\sigma _{2}+\sigma _{3}){\big )}\,,\\\varepsilon _{2}&={\frac {1}{E}}{\big (}\sigma _{2}-\nu (\sigma _{1}+\sigma _{3}){\big )}\,,\\\varepsilon _{3}&={\frac {1}{E}}{\big (}\sigma _{3}-\nu (\sigma _{1}+\sigma _{2}){\big )}\,,\end{aligned}}}
ε
1
=
1
E
(
(
1
+
ν
)
σ
1
−
ν
(
σ
1
+
σ
2
+
σ
3
)
)
,
ε
2
=
1
E
(
(
1
+
ν
)
σ
2
−
ν
(
σ
1
+
σ
2
+
σ
3
)
)
,
ε
3
=
1
E
(
(
1
+
ν
)
σ
3
−
ν
(
σ
1
+
σ
2
+
σ
3
)
)
,
{\displaystyle {\begin{aligned}\varepsilon _{1}&={\frac {1}{E}}{\big (}(1+\nu )\sigma _{1}-\nu (\sigma _{1}+\sigma _{2}+\sigma _{3}){\big )}\,,\\\varepsilon _{2}&={\frac {1}{E}}{\big (}(1+\nu )\sigma _{2}-\nu (\sigma _{1}+\sigma _{2}+\sigma _{3}){\big )}\,,\\\varepsilon _{3}&={\frac {1}{E}}{\big (}(1+\nu )\sigma _{3}-\nu (\sigma _{1}+\sigma _{2}+\sigma _{3}){\big )}\,,\end{aligned}}}
σ
1
=
E
1
+
ν
ε
1
+
ν
1
+
ν
(
σ
1
+
σ
2
+
σ
3
)
.
{\displaystyle \sigma _{1}={\frac {E}{1+\nu }}\varepsilon _{1}+{\frac {\nu }{1+\nu }}(\sigma _{1}+\sigma _{2}+\sigma _{3})\,.}
合計を計算し
、それを σ 1
について解いた方程式に代入すると、次の ようになります
。
ここで、 μ と λ は ラメパラメータ です 。
ε
1
+
ε
2
+
ε
3
=
1
E
(
(
1
+
ν
)
(
σ
1
+
σ
2
+
σ
3
)
−
3
ν
(
σ
1
+
σ
2
+
σ
3
)
)
=
1
−
2
ν
E
(
σ
1
+
σ
2
+
σ
3
)
σ
1
+
σ
2
+
σ
3
=
E
1
−
2
ν
(
ε
1
+
ε
2
+
ε
3
)
{\displaystyle {\begin{aligned}\varepsilon _{1}+\varepsilon _{2}+\varepsilon _{3}&={\frac {1}{E}}{\big (}(1+\nu )(\sigma _{1}+\sigma _{2}+\sigma _{3})-3\nu (\sigma _{1}+\sigma _{2}+\sigma _{3}){\big )}={\frac {1-2\nu }{E}}(\sigma _{1}+\sigma _{2}+\sigma _{3})\\\sigma _{1}+\sigma _{2}+\sigma _{3}&={\frac {E}{1-2\nu }}(\varepsilon _{1}+\varepsilon _{2}+\varepsilon _{3})\end{aligned}}}
σ
1
=
E
1
+
ν
ε
1
+
E
ν
(
1
+
ν
)
(
1
−
2
ν
)
(
ε
1
+
ε
2
+
ε
3
)
=
2
μ
ε
1
+
λ
(
ε
1
+
ε
2
+
ε
3
)
,
{\displaystyle {\begin{aligned}\sigma _{1}&={\frac {E}{1+\nu }}\varepsilon _{1}+{\frac {E\nu }{(1+\nu )(1-2\nu )}}(\varepsilon _{1}+\varepsilon _{2}+\varepsilon _{3})\\&=2\mu \varepsilon _{1}+\lambda (\varepsilon _{1}+\varepsilon _{2}+\varepsilon _{3})\,,\end{aligned}}}
方向 2 と 3 を同様に処理すると、3 次元でのフックの法則が得られます。
行列形式では、等方性材料に対するフックの法則は次のように表すことができます
。ここで、 γ ij = 2 ε ij は工学的せん断ひずみ です 。逆の関係は次のように表すことができます
。これは
、ラメ定数のおかげで簡略化できます。
ベクトル表記では、これは次のように表されます。
ここで、 I は恒等テンソルです。
[
ε
11
ε
22
ε
33
2
ε
23
2
ε
13
2
ε
12
]
=
[
ε
11
ε
22
ε
33
γ
23
γ
13
γ
12
]
=
1
E
[
1
−
ν
−
ν
0
0
0
−
ν
1
−
ν
0
0
0
−
ν
−
ν
1
0
0
0
0
0
0
2
+
2
ν
0
0
0
0
0
0
2
+
2
ν
0
0
0
0
0
0
2
+
2
ν
]
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
{\displaystyle {\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}}\,=\,{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\\gamma _{23}\\\gamma _{13}\\\gamma _{12}\end{bmatrix}}\,=\,{\frac {1}{E}}{\begin{bmatrix}1&-\nu &-\nu &0&0&0\\-\nu &1&-\nu &0&0&0\\-\nu &-\nu &1&0&0&0\\0&0&0&2+2\nu &0&0\\0&0&0&0&2+2\nu &0\\0&0&0&0&0&2+2\nu \end{bmatrix}}{\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}}}
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
=
E
(
1
+
ν
)
(
1
−
2
ν
)
[
1
−
ν
ν
ν
0
0
0
ν
1
−
ν
ν
0
0
0
ν
ν
1
−
ν
0
0
0
0
0
0
1
−
2
ν
2
0
0
0
0
0
0
1
−
2
ν
2
0
0
0
0
0
0
1
−
2
ν
2
]
[
ε
11
ε
22
ε
33
2
ε
23
2
ε
13
2
ε
12
]
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}}\,=\,{\frac {E}{(1+\nu )(1-2\nu )}}{\begin{bmatrix}1-\nu &\nu &\nu &0&0&0\\\nu &1-\nu &\nu &0&0&0\\\nu &\nu &1-\nu &0&0&0\\0&0&0&{\frac {1-2\nu }{2}}&0&0\\0&0&0&0&{\frac {1-2\nu }{2}}&0\\0&0&0&0&0&{\frac {1-2\nu }{2}}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}}}
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
=
[
2
μ
+
λ
λ
λ
0
0
0
λ
2
μ
+
λ
λ
0
0
0
λ
λ
2
μ
+
λ
0
0
0
0
0
0
μ
0
0
0
0
0
0
μ
0
0
0
0
0
0
μ
]
[
ε
11
ε
22
ε
33
2
ε
23
2
ε
13
2
ε
12
]
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}}\,=\,{\begin{bmatrix}2\mu +\lambda &\lambda &\lambda &0&0&0\\\lambda &2\mu +\lambda &\lambda &0&0&0\\\lambda &\lambda &2\mu +\lambda &0&0&0\\0&0&0&\mu &0&0\\0&0&0&0&\mu &0\\0&0&0&0&0&\mu \end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}}}
[
σ
11
σ
12
σ
13
σ
12
σ
22
σ
23
σ
13
σ
23
σ
33
]
=
2
μ
[
ε
11
ε
12
ε
13
ε
12
ε
22
ε
23
ε
13
ε
23
ε
33
]
+
λ
I
(
ε
11
+
ε
22
+
ε
33
)
{\displaystyle {\begin{bmatrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{12}&\sigma _{22}&\sigma _{23}\\\sigma _{13}&\sigma _{23}&\sigma _{33}\end{bmatrix}}\,=\,2\mu {\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{12}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{13}&\varepsilon _{23}&\varepsilon _{33}\end{bmatrix}}+\lambda \mathbf {I} \left(\varepsilon _{11}+\varepsilon _{22}+\varepsilon _{33}\right)}
平面応力
平面応力 条件下では 、 σ 31 = σ 13 = σ 32 = σ 23 = σ 33 = 0 となる。この場合、フックの法則は次のようになる。
[
σ
11
σ
22
σ
12
]
=
E
1
−
ν
2
[
1
ν
0
ν
1
0
0
0
1
−
ν
2
]
[
ε
11
ε
22
2
ε
12
]
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}\,=\,{\frac {E}{1-\nu ^{2}}}{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&{\frac {1-\nu }{2}}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\2\varepsilon _{12}\end{bmatrix}}}
ベクトル表記では、これは次のようになる。
[
σ
11
σ
12
σ
12
σ
22
]
=
E
1
−
ν
2
(
(
1
−
ν
)
[
ε
11
ε
12
ε
12
ε
22
]
+
ν
I
(
ε
11
+
ε
22
)
)
{\displaystyle {\begin{bmatrix}\sigma _{11}&\sigma _{12}\\\sigma _{12}&\sigma _{22}\end{bmatrix}}\,=\,{\frac {E}{1-\nu ^{2}}}\left((1-\nu ){\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}\\\varepsilon _{12}&\varepsilon _{22}\end{bmatrix}}+\nu \mathbf {I} \left(\varepsilon _{11}+\varepsilon _{22}\right)\right)}
逆関係は通常、簡約形式で表されます。
[
ε
11
ε
22
2
ε
12
]
=
1
E
[
1
−
ν
0
−
ν
1
0
0
0
2
+
2
ν
]
[
σ
11
σ
22
σ
12
]
{\displaystyle {\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\2\varepsilon _{12}\end{bmatrix}}\,=\,{\frac {1}{E}}{\begin{bmatrix}1&-\nu &0\\-\nu &1&0\\0&0&2+2\nu \end{bmatrix}}{\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}}
平面ひずみ
平面ひずみ 条件では 、 ε 31 = ε 13 = ε 32 = ε 23 = ε 33 = 0 となる。この場合、フックの法則は次のようになる。
[
σ
11
σ
22
σ
12
]
=
E
(
1
+
ν
)
(
1
−
2
ν
)
[
1
−
ν
ν
0
ν
1
−
ν
0
0
0
1
−
2
ν
2
]
[
ε
11
ε
22
2
ε
12
]
{\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}\,=\,{\frac {E}{(1+\nu )(1-2\nu )}}{\begin{bmatrix}1-\nu &\nu &0\\\nu &1-\nu &0\\0&0&{\frac {1-2\nu }{2}}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\2\varepsilon _{12}\end{bmatrix}}}
異方性材料
コーシー応力テンソル の対称性 ( σ ij = σ ji ) と一般化フックの法則 ( σ ij = c ijkl ε kl )は、 c ijkl = c jikl を意味します。同様に、微小ひずみテンソル の対称性は、 c ijkl = c ijlk を意味します。これらの対称性は、 剛性テンソル cの マイナー対称性 と呼ばれます 。これにより、弾性定数の数が 81 から 36 に減少します。
さらに、変位勾配とコーシー応力は仕事共役であるため、応力-ひずみ関係はひずみエネルギー密度汎関数 ( U ) から導くことができ、
微分次数の任意性は c ijkl = c klij を意味します。これらは剛性テンソルの主要な対称性 と呼ばれます 。これにより、弾性定数の数が 36 から 21 に減ります。主要な対称性とマイナーな対称性は、剛性テンソルに 21 個の独立したコンポーネントしかないことを示しています。
σ
i
j
=
∂
U
∂
ε
i
j
⟹
c
i
j
k
l
=
∂
2
U
∂
ε
i
j
∂
ε
k
l
.
{\displaystyle \sigma _{ij}={\frac {\partial U}{\partial \varepsilon _{ij}}}\quad \implies \quad c_{ijkl}={\frac {\partial ^{2}U}{\partial \varepsilon _{ij}\partial \varepsilon _{kl}}}\,.}
行列表現(剛性テンソル)
フックの法則の異方性形式を行列表記( フォークト表記 とも呼ばれる)で表現すると便利な場合が多い。これを行うには、応力テンソルとひずみテンソルの対称性を利用し、それらを正規直交座標系( e 1 、 e 2 、 e 3 )
の6次元ベクトルとして表現する。
すると、剛性テンソル( c )は次のように表され
、フックの法則は次のように書かれる。
同様に、コンプライアンステンソル( s )は次のように書かれる。
[
σ
]
=
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
≡
[
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
]
;
[
ε
]
=
[
ε
11
ε
22
ε
33
2
ε
23
2
ε
13
2
ε
12
]
≡
[
ε
1
ε
2
ε
3
ε
4
ε
5
ε
6
]
{\displaystyle [{\boldsymbol {\sigma }}]\,=\,{\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}}\,\equiv \,{\begin{bmatrix}\sigma _{1}\\\sigma _{2}\\\sigma _{3}\\\sigma _{4}\\\sigma _{5}\\\sigma _{6}\end{bmatrix}}\,;\qquad [{\boldsymbol {\varepsilon }}]\,=\,{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}}\,\equiv \,{\begin{bmatrix}\varepsilon _{1}\\\varepsilon _{2}\\\varepsilon _{3}\\\varepsilon _{4}\\\varepsilon _{5}\\\varepsilon _{6}\end{bmatrix}}}
[
c
]
=
[
c
1111
c
1122
c
1133
c
1123
c
1131
c
1112
c
2211
c
2222
c
2233
c
2223
c
2231
c
2212
c
3311
c
3322
c
3333
c
3323
c
3331
c
3312
c
2311
c
2322
c
2333
c
2323
c
2331
c
2312
c
3111
c
3122
c
3133
c
3123
c
3131
c
3112
c
1211
c
1222
c
1233
c
1223
c
1231
c
1212
]
≡
[
C
11
C
12
C
13
C
14
C
15
C
16
C
12
C
22
C
23
C
24
C
25
C
26
C
13
C
23
C
33
C
34
C
35
C
36
C
14
C
24
C
34
C
44
C
45
C
46
C
15
C
25
C
35
C
45
C
55
C
56
C
16
C
26
C
36
C
46
C
56
C
66
]
{\displaystyle [{\mathsf {c}}]\,=\,{\begin{bmatrix}c_{1111}&c_{1122}&c_{1133}&c_{1123}&c_{1131}&c_{1112}\\c_{2211}&c_{2222}&c_{2233}&c_{2223}&c_{2231}&c_{2212}\\c_{3311}&c_{3322}&c_{3333}&c_{3323}&c_{3331}&c_{3312}\\c_{2311}&c_{2322}&c_{2333}&c_{2323}&c_{2331}&c_{2312}\\c_{3111}&c_{3122}&c_{3133}&c_{3123}&c_{3131}&c_{3112}\\c_{1211}&c_{1222}&c_{1233}&c_{1223}&c_{1231}&c_{1212}\end{bmatrix}}\,\equiv \,{\begin{bmatrix}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{16}\\C_{12}&C_{22}&C_{23}&C_{24}&C_{25}&C_{26}\\C_{13}&C_{23}&C_{33}&C_{34}&C_{35}&C_{36}\\C_{14}&C_{24}&C_{34}&C_{44}&C_{45}&C_{46}\\C_{15}&C_{25}&C_{35}&C_{45}&C_{55}&C_{56}\\C_{16}&C_{26}&C_{36}&C_{46}&C_{56}&C_{66}\end{bmatrix}}}
[
σ
]
=
[
C
]
[
ε
]
or
σ
i
=
C
i
j
ε
j
.
{\displaystyle [{\boldsymbol {\sigma }}]=[{\mathsf {C}}][{\boldsymbol {\varepsilon }}]\qquad {\text{or}}\qquad \sigma _{i}=C_{ij}\varepsilon _{j}\,.}
[
s
]
=
[
s
1111
s
1122
s
1133
2
s
1123
2
s
1131
2
s
1112
s
2211
s
2222
s
2233
2
s
2223
2
s
2231
2
s
2212
s
3311
s
3322
s
3333
2
s
3323
2
s
3331
2
s
3312
2
s
2311
2
s
2322
2
s
2333
4
s
2323
4
s
2331
4
s
2312
2
s
3111
2
s
3122
2
s
3133
4
s
3123
4
s
3131
4
s
3112
2
s
1211
2
s
1222
2
s
1233
4
s
1223
4
s
1231
4
s
1212
]
≡
[
S
11
S
12
S
13
S
14
S
15
S
16
S
12
S
22
S
23
S
24
S
25
S
26
S
13
S
23
S
33
S
34
S
35
S
36
S
14
S
24
S
34
S
44
S
45
S
46
S
15
S
25
S
35
S
45
S
55
S
56
S
16
S
26
S
36
S
46
S
56
S
66
]
{\displaystyle [{\mathsf {s}}]\,=\,{\begin{bmatrix}s_{1111}&s_{1122}&s_{1133}&2s_{1123}&2s_{1131}&2s_{1112}\\s_{2211}&s_{2222}&s_{2233}&2s_{2223}&2s_{2231}&2s_{2212}\\s_{3311}&s_{3322}&s_{3333}&2s_{3323}&2s_{3331}&2s_{3312}\\2s_{2311}&2s_{2322}&2s_{2333}&4s_{2323}&4s_{2331}&4s_{2312}\\2s_{3111}&2s_{3122}&2s_{3133}&4s_{3123}&4s_{3131}&4s_{3112}\\2s_{1211}&2s_{1222}&2s_{1233}&4s_{1223}&4s_{1231}&4s_{1212}\end{bmatrix}}\,\equiv \,{\begin{bmatrix}S_{11}&S_{12}&S_{13}&S_{14}&S_{15}&S_{16}\\S_{12}&S_{22}&S_{23}&S_{24}&S_{25}&S_{26}\\S_{13}&S_{23}&S_{33}&S_{34}&S_{35}&S_{36}\\S_{14}&S_{24}&S_{34}&S_{44}&S_{45}&S_{46}\\S_{15}&S_{25}&S_{35}&S_{45}&S_{55}&S_{56}\\S_{16}&S_{26}&S_{36}&S_{46}&S_{56}&S_{66}\end{bmatrix}}}
座標系の変更
線形弾性材料を基準構成から別の構成に回転させた場合、回転構成における剛性テンソルの成分が基準構成の成分と関係 [13]
によって関連付けられると、材料は回転に対して対称になります。
ここで、 l abは 直交回転行列 [ L ] の成分です 。同じ関係は反転にも当てはまります。
c
p
q
r
s
=
l
p
i
l
q
j
l
r
k
l
s
l
c
i
j
k
l
{\displaystyle c_{pqrs}=l_{pi}l_{qj}l_{rk}l_{sl}c_{ijkl}}
行列表記では、変換された基底(回転または反転)が参照基底と
次
のように関係付けられている場合、さらに、物質が変換 [ L ]
に関して対称である場合 、
[
e
i
′
]
=
[
L
]
[
e
i
]
{\displaystyle [\mathbf {e} _{i}']=[L][\mathbf {e} _{i}]}
C
i
j
ε
i
ε
j
=
C
i
j
′
ε
i
′
ε
j
′
.
{\displaystyle C_{ij}\varepsilon _{i}\varepsilon _{j}=C_{ij}'\varepsilon '_{i}\varepsilon '_{j}\,.}
C
i
j
=
C
i
j
′
⟹
C
i
j
(
ε
i
ε
j
−
ε
i
′
ε
j
′
)
=
0
.
{\displaystyle C_{ij}=C'_{ij}\quad \implies \quad C_{ij}(\varepsilon _{i}\varepsilon _{j}-\varepsilon '_{i}\varepsilon '_{j})=0\,.}
直交異方性材料
直交異方性材料は3つの 直交対称 面 を持つ 。基底ベクトル( e 1 、 e 2 、 e 3 )が対称面の法線である場合、座標変換関係は次式を意味する。
この関係の逆は一般に次のように表される [14] [ ページ必要 ]
。
[
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
]
=
[
C
11
C
12
C
13
0
0
0
C
12
C
22
C
23
0
0
0
C
13
C
23
C
33
0
0
0
0
0
0
C
44
0
0
0
0
0
0
C
55
0
0
0
0
0
0
C
66
]
[
ε
1
ε
2
ε
3
ε
4
ε
5
ε
6
]
{\displaystyle {\begin{bmatrix}\sigma _{1}\\\sigma _{2}\\\sigma _{3}\\\sigma _{4}\\\sigma _{5}\\\sigma _{6}\end{bmatrix}}\,=\,{\begin{bmatrix}C_{11}&C_{12}&C_{13}&0&0&0\\C_{12}&C_{22}&C_{23}&0&0&0\\C_{13}&C_{23}&C_{33}&0&0&0\\0&0&0&C_{44}&0&0\\0&0&0&0&C_{55}&0\\0&0&0&0&0&C_{66}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{1}\\\varepsilon _{2}\\\varepsilon _{3}\\\varepsilon _{4}\\\varepsilon _{5}\\\varepsilon _{6}\end{bmatrix}}}
[
ε
x
x
ε
y
y
ε
z
z
2
ε
y
z
2
ε
z
x
2
ε
x
y
]
=
[
1
E
x
−
ν
y
x
E
y
−
ν
z
x
E
z
0
0
0
−
ν
x
y
E
x
1
E
y
−
ν
z
y
E
z
0
0
0
−
ν
x
z
E
x
−
ν
y
z
E
y
1
E
z
0
0
0
0
0
0
1
G
y
z
0
0
0
0
0
0
1
G
z
x
0
0
0
0
0
0
1
G
x
y
]
[
σ
x
x
σ
y
y
σ
z
z
σ
y
z
σ
z
x
σ
x
y
]
{\displaystyle {\begin{bmatrix}\varepsilon _{xx}\\\varepsilon _{yy}\\\varepsilon _{zz}\\2\varepsilon _{yz}\\2\varepsilon _{zx}\\2\varepsilon _{xy}\end{bmatrix}}\,=\,{\begin{bmatrix}{\frac {1}{E_{x}}}&-{\frac {\nu _{yx}}{E_{y}}}&-{\frac {\nu _{zx}}{E_{z}}}&0&0&0\\-{\frac {\nu _{xy}}{E_{x}}}&{\frac {1}{E_{y}}}&-{\frac {\nu _{zy}}{E_{z}}}&0&0&0\\-{\frac {\nu _{xz}}{E_{x}}}&-{\frac {\nu _{yz}}{E_{y}}}&{\frac {1}{E_{z}}}&0&0&0\\0&0&0&{\frac {1}{G_{yz}}}&0&0\\0&0&0&0&{\frac {1}{G_{zx}}}&0\\0&0&0&0&0&{\frac {1}{G_{xy}}}\\\end{bmatrix}}{\begin{bmatrix}\sigma _{xx}\\\sigma _{yy}\\\sigma _{zz}\\\sigma _{yz}\\\sigma _{zx}\\\sigma _{xy}\end{bmatrix}}}
E i は 軸 iに沿った ヤング率 である。
G ij は 、法線が方向i である平面上の 方向 jの せん断弾性係数 である。
ν ij は、 方向i に拡張が適用されたときに方向 j の収縮に対応する ポアソン比 です 。
平面応力 条件下 では、 σ zz = σ zx = σ yz = 0 、直交異方性材料に対するフックの法則は次の形式になります
。
逆の関係は、
上記の剛性マトリックスの転置形式もよく使用されます。
[
ε
x
x
ε
y
y
2
ε
x
y
]
=
[
1
E
x
−
ν
y
x
E
y
0
−
ν
x
y
E
x
1
E
y
0
0
0
1
G
x
y
]
[
σ
x
x
σ
y
y
σ
x
y
]
.
{\displaystyle {\begin{bmatrix}\varepsilon _{xx}\\\varepsilon _{yy}\\2\varepsilon _{xy}\end{bmatrix}}\,=\,{\begin{bmatrix}{\frac {1}{E_{x}}}&-{\frac {\nu _{yx}}{E_{y}}}&0\\-{\frac {\nu _{xy}}{E_{x}}}&{\frac {1}{E_{y}}}&0\\0&0&{\frac {1}{G_{xy}}}\end{bmatrix}}{\begin{bmatrix}\sigma _{xx}\\\sigma _{yy}\\\sigma _{xy}\end{bmatrix}}\,.}
[
σ
x
x
σ
y
y
σ
x
y
]
=
1
1
−
ν
x
y
ν
y
x
[
E
x
ν
y
x
E
x
0
ν
x
y
E
y
E
y
0
0
0
G
x
y
(
1
−
ν
x
y
ν
y
x
)
]
[
ε
x
x
ε
y
y
2
ε
x
y
]
.
{\displaystyle {\begin{bmatrix}\sigma _{xx}\\\sigma _{yy}\\\sigma _{xy}\end{bmatrix}}\,=\,{\frac {1}{1-\nu _{xy}\nu _{yx}}}{\begin{bmatrix}E_{x}&\nu _{yx}E_{x}&0\\\nu _{xy}E_{y}&E_{y}&0\\0&0&G_{xy}(1-\nu _{xy}\nu _{yx})\end{bmatrix}}{\begin{bmatrix}\varepsilon _{xx}\\\varepsilon _{yy}\\2\varepsilon _{xy}\end{bmatrix}}\,.}
横方向等方性材料
横方向等方性 材料は、 対称軸の 周りの回転に関して対称である 。このような材料の場合、 e 3 が 対称軸であれば、フックの法則は次のように表される。
[
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
]
=
[
C
11
C
12
C
13
0
0
0
C
12
C
11
C
13
0
0
0
C
13
C
13
C
33
0
0
0
0
0
0
C
44
0
0
0
0
0
0
C
44
0
0
0
0
0
0
C
11
−
C
12
2
]
[
ε
1
ε
2
ε
3
ε
4
ε
5
ε
6
]
{\displaystyle {\begin{bmatrix}\sigma _{1}\\\sigma _{2}\\\sigma _{3}\\\sigma _{4}\\\sigma _{5}\\\sigma _{6}\end{bmatrix}}\,=\,{\begin{bmatrix}C_{11}&C_{12}&C_{13}&0&0&0\\C_{12}&C_{11}&C_{13}&0&0&0\\C_{13}&C_{13}&C_{33}&0&0&0\\0&0&0&C_{44}&0&0\\0&0&0&0&C_{44}&0\\0&0&0&0&0&{\frac {C_{11}-C_{12}}{2}}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{1}\\\varepsilon _{2}\\\varepsilon _{3}\\\varepsilon _{4}\\\varepsilon _{5}\\\varepsilon _{6}\end{bmatrix}}}
より一般的には、 x ≡ e 1 軸が対称軸とされ、逆フックの法則は次のように記述される
[15]。
[
ε
x
x
ε
y
y
ε
z
z
2
ε
y
z
2
ε
z
x
2
ε
x
y
]
=
[
1
E
x
−
ν
y
x
E
y
−
ν
z
x
E
z
0
0
0
−
ν
x
y
E
x
1
E
y
−
ν
z
y
E
z
0
0
0
−
ν
x
z
E
x
−
ν
y
z
E
y
1
E
z
0
0
0
0
0
0
1
G
y
z
0
0
0
0
0
0
1
G
x
z
0
0
0
0
0
0
1
G
x
y
]
[
σ
x
x
σ
y
y
σ
z
z
σ
y
z
σ
z
x
σ
x
y
]
{\displaystyle {\begin{bmatrix}\varepsilon _{xx}\\\varepsilon _{yy}\\\varepsilon _{zz}\\2\varepsilon _{yz}\\2\varepsilon _{zx}\\2\varepsilon _{xy}\end{bmatrix}}\,=\,{\begin{bmatrix}{\frac {1}{E_{x}}}&-{\frac {\nu _{yx}}{E_{y}}}&-{\frac {\nu _{zx}}{E_{z}}}&0&0&0\\-{\frac {\nu _{xy}}{E_{x}}}&{\frac {1}{E_{y}}}&-{\frac {\nu _{zy}}{E_{z}}}&0&0&0\\-{\frac {\nu _{xz}}{E_{x}}}&-{\frac {\nu _{yz}}{E_{y}}}&{\frac {1}{E_{z}}}&0&0&0\\0&0&0&{\frac {1}{G_{yz}}}&0&0\\0&0&0&0&{\frac {1}{G_{xz}}}&0\\0&0&0&0&0&{\frac {1}{G_{xy}}}\\\end{bmatrix}}{\begin{bmatrix}\sigma _{xx}\\\sigma _{yy}\\\sigma _{zz}\\\sigma _{yz}\\\sigma _{zx}\\\sigma _{xy}\end{bmatrix}}}
普遍弾性異方性指数
あらゆるクラスの異方性の程度を把握するために、 普遍弾性異方性指数 (AU) [16] が策定されました。これは、 立方晶系結晶 に適した ツェナー 比に代わるものです。
熱力学的基礎
弾性体の線形変形は、断熱変形 として近似できます 。これらの条件下および準静的プロセスでは、変形した物体に対する 熱力学の第一法則は 次のように表すことができます
。
ここ
で、 δUは 内部エネルギー の増加 、 δW は 外部 力による仕事 です 。仕事は 2 つの項に分割できます。δW
s は表面力 による仕事で 、 δW b は 体積力 による仕事です 。δ u が 物体内の変位場 u の変化 である場合 、2 つの外部仕事項は
次のように 表すことができます。ここで、 t は表面 牽引 ベクトル、 b は体積力ベクトル、 Ω は物体、 ∂ Ω は 物体の表面を表します。 コーシー応力 と表面牽引力 の関係 t = n · σ (ここで n は∂ Ω の外向きの法線の単位 )
を用いると、次の式 が得られます
。 発散定理 を用いて
表面積 分を 体積積分 に変換すると、次の式 が得られます。
コーシー応力の対称性と恒等式を用いると
、
次の式が得られます。
ひずみ
の定義と 平衡 方程式から、次の式が 得られます
。
したがって、次のように書くことができます
。 したがって、 内部エネルギー 密度
の変化は次のように表されます。
弾性 材料は、総内部エネルギーが 内部力の 位置エネルギー( 弾性ひずみエネルギー とも呼ばれる) に等しい材料として定義されます。 したがって、内部エネルギー密度はひずみの関数 U 0 = U 0 ( ε ) であり、内部エネルギーの変化は次のように表すことができます。
ひずみの変化は任意であるため、弾性材料の応力とひずみの関係は次のように表されます。
線形弾性材料の場合、量
δ
W
=
δ
U
{\displaystyle \delta W=\delta U}
δ
W
=
δ
W
s
+
δ
W
b
{\displaystyle \delta W=\delta W_{\mathrm {s} }+\delta W_{\mathrm {b} }}
δ
W
s
=
∫
∂
Ω
t
⋅
δ
u
d
S
;
δ
W
b
=
∫
Ω
b
⋅
δ
u
d
V
{\displaystyle \delta W_{\mathrm {s} }=\int _{\partial \Omega }\mathbf {t} \cdot \delta \mathbf {u} \,dS\,;\qquad \delta W_{\mathrm {b} }=\int _{\Omega }\mathbf {b} \cdot \delta \mathbf {u} \,dV}
δ
W
=
δ
U
=
∫
∂
Ω
(
n
⋅
σ
)
⋅
δ
u
d
S
+
∫
Ω
b
⋅
δ
u
d
V
.
{\displaystyle \delta W=\delta U=\int _{\partial \Omega }(\mathbf {n} \cdot {\boldsymbol {\sigma }})\cdot \delta \mathbf {u} \,dS+\int _{\Omega }\mathbf {b} \cdot \delta \mathbf {u} \,dV\,.}
δ
U
=
∫
Ω
(
∇
⋅
(
σ
⋅
δ
u
)
+
b
⋅
δ
u
)
d
V
.
{\displaystyle \delta U=\int _{\Omega }{\big (}\nabla \cdot ({\boldsymbol {\sigma }}\cdot \delta \mathbf {u} )+\mathbf {b} \cdot \delta \mathbf {u} {\big )}\,dV\,.}
∇
⋅
(
a
⋅
b
)
=
(
∇
⋅
a
)
⋅
b
+
1
2
(
a
T
:
∇
b
+
a
:
(
∇
b
)
T
)
{\displaystyle \nabla \cdot (\mathbf {a} \cdot \mathbf {b} )=(\nabla \cdot \mathbf {a} )\cdot \mathbf {b} +{\tfrac {1}{2}}\left(\mathbf {a} ^{\mathsf {T}}:\nabla \mathbf {b} +\mathbf {a} :(\nabla \mathbf {b} )^{\mathsf {T}}\right)}
δ
U
=
∫
Ω
(
σ
:
1
2
(
∇
δ
u
+
(
∇
δ
u
)
T
)
+
(
∇
⋅
σ
+
b
)
⋅
δ
u
)
d
V
.
{\displaystyle \delta U=\int _{\Omega }\left({\boldsymbol {\sigma }}:{\tfrac {1}{2}}\left(\nabla \delta \mathbf {u} +(\nabla \delta \mathbf {u} )^{\mathsf {T}}\right)+\left(\nabla \cdot {\boldsymbol {\sigma }}+\mathbf {b} \right)\cdot \delta \mathbf {u} \right)\,dV\,.}
δ
ε
=
1
2
(
∇
δ
u
+
(
∇
δ
u
)
T
)
;
∇
⋅
σ
+
b
=
0
.
{\displaystyle \delta {\boldsymbol {\varepsilon }}={\tfrac {1}{2}}\left(\nabla \delta \mathbf {u} +(\nabla \delta \mathbf {u} )^{\mathsf {T}}\right)\,;\qquad \nabla \cdot {\boldsymbol {\sigma }}+\mathbf {b} =\mathbf {0} \,.}
δ
U
=
∫
Ω
σ
:
δ
ε
d
V
{\displaystyle \delta U=\int _{\Omega }{\boldsymbol {\sigma }}:\delta {\boldsymbol {\varepsilon }}\,dV}
δ
U
0
=
σ
:
δ
ε
.
{\displaystyle \delta U_{0}={\boldsymbol {\sigma }}:\delta {\boldsymbol {\varepsilon }}\,.}
δ
U
0
=
∂
U
0
∂
ε
:
δ
ε
.
{\displaystyle \delta U_{0}={\frac {\partial U_{0}}{\partial {\boldsymbol {\varepsilon }}}}:\delta {\boldsymbol {\varepsilon }}\,.}
σ
=
∂
U
0
∂
ε
.
{\displaystyle {\boldsymbol {\sigma }}={\frac {\partial U_{0}}{\partial {\boldsymbol {\varepsilon }}}}\,.}
∂U0 の / ∂ε は ε の線形関数な ので、次のように表すことができます
ここで、 c は 材料定数の 4 階テンソルで、 剛性テンソル とも呼ばれます。 線形弾性材料の場合、 c が 4 階テンソルでなければならない理由は、次の点に注目することで
インデックス表記では
、右側の定数には 4 つのインデックスが必要で、4 階の量です。また、この量は、ひずみテンソルを応力テンソルに変換する線形変換であるため、テンソルでなければならないこともわかります。また、定数が 4 階テンソルのテンソル変換規則に従っていることも示せます。
σ
=
c
:
ε
{\displaystyle {\boldsymbol {\sigma }}={\mathsf {c}}:{\boldsymbol {\varepsilon }}}
∂
∂
ε
σ
(
ε
)
=
constant
=
c
.
{\displaystyle {\frac {\partial }{\partial {\boldsymbol {\varepsilon }}}}{\boldsymbol {\sigma }}({\boldsymbol {\varepsilon }})={\text{constant}}={\mathsf {c}}\,.}
∂
σ
i
j
∂
ε
k
l
=
constant
=
c
i
j
k
l
.
{\displaystyle {\frac {\partial \sigma _{ij}}{\partial \varepsilon _{kl}}}={\text{constant}}=c_{ijkl}\,.}
参照
注記
^ アナグラムはアルファベット順に、 ceiiinosssttuv と 書かれており、 Ut tensio, sic vis (「拡張すれば力も広がる」) を表しています。Petroski , Henry (1996)。Invention by Design: How Engineers Get from Thought to Thing . Cambridge, MA: Harvard University Press. p. 11. ISBN 978-0674463684 。
^ http://civil.lindahall.org/design.shtml を参照してください。ここでは catenary のアナグラムも見つかります 。
^ ロバート・フック 、 『復元力について』、または『バネの力について。バネの力を説明する』 、ロンドン、1678年。
^ ヤング、ヒュー・D.、フリードマン、ロジャー・A.、フォード、A.ルイス(2016年)。 シアーズとゼマンスキーの大学物理学:現代物理学とともに (第14版)。ピアソン。209ページ。
^ 牛場 翔太; 増井 京子; 田口 夏生; 浜野 智樹; 河田 聡; 庄司 悟 (2015). 「コイルスプリング形状ポリマーナノワイヤのサイズ依存ナノメカニクス」. Scientific Reports . 5 :17152. Bibcode :2015NatSR...517152U. doi :10.1038/srep17152. PMC 4661696. PMID 26612544 .
^ ベレンキー;サラエフ(1988)。 「層状結晶の変形効果」。 ウスペキ・フィジチェスキフ・ナウク 。 155 (5): 89. 土井 : 10.3367/UFNr.0155.198805c.0089 。
^ Mouhat, Félix; Coudert, François-Xavier (2014 年 12 月 5 日). 「さまざまな結晶系における必要かつ十分な弾性安定条件」. Physical Review B. 90 ( 22): 224104. arXiv : 1410.0065 . Bibcode :2014PhRvB..90v4104M. doi :10.1103/PhysRevB.90.224104. ISSN 1098-0121. S2CID 54058316.
^ Vijay Madhav, M.; Manogaran, S. (2009). 「冗長内部座標におけるコンプライアンス定数の再検討といくつかの新しい知見」 J. Chem. Phys . 131 (17): 174112–174116. Bibcode :2009JChPh.131q4112V. doi :10.1063/1.3259834. PMID 19895003.
^ Ponomareva, Alla; Yurenko, Yevgen; Zhurakivsky, Roman; Van Mourik, Tanja; Hovorun, Dmytro (2012). 「潜在的な HIV-1 逆転写酵素阻害剤 d4U および d4C の完全なコンフォメーション空間。量子化学研究」 。Phys . Chem. Chem. Phys . 14 (19): 6787–6795。Bibcode : 2012PCCP ...14.6787P。doi :10.1039/C2CP40290D。PMID 22461011 。
^ Symon, Keith R. (1971). 「第10章」. 力学 . マサチューセッツ州レディング: Addison-Wesley. ISBN 9780201073928 。
^ Simo, JC; Hughes, TJR (1998). 計算非弾性 . Springer. ISBN 9780387975207 。
^ ミルトン、グレアム W. (2002)。 複合理論 。応用数学と計算数学に関するケンブリッジモノグラフ。ケンブリッジ大学出版局 。ISBN 9780521781251 。
^ スローター、ウィリアム S. (2001)。 弾性の線形化理論 。ビルクハウザー 。ISBN 978-0817641177 。
^ Boresi, AP; Schmidt, RJ; Sidebottom, OM (1993). Advanced Mechanics of Materials (第 5 版). Wiley. ISBN 9780471600091 。
^ Tan, SC (1994). 積層複合材の応力集中 . ランカスター、ペンシルバニア州: Technomic Publishing Company. ISBN 9781566760775 。
^ Ranganathan, SI; Ostoja-Starzewski, M. (2008). 「普遍弾性異方性指数」. Physical Review Letters . 101 (5): 055504–1–4. Bibcode :2008PhRvL.101e5504R. doi :10.1103/PhysRevLett.101.055504. PMID 18764407. S2CID 6668703.
参考文献
フックの法則 - ファインマン物理学講義
フックの法則 - 古典力学 - 物理学 - MIT OpenCourseWare
外部リンク
スプリングとフックの法則を説明する JavaScript アプレット
スプリングフォースをデモンストレーションする JavaScript アプレット