スカラー、ベクトル、および2次テンソルの2次テンソルに対する微分は、連続体力学において非常に有用である。これらの微分は、非線形弾性および塑性理論、特に数値シミュレーションのアルゴリズムの設計において使用される。[1]
方向微分はこれらの微分を求める体系的な方法を提供する。[2]
ベクトルと2階テンソルに関する微分
さまざまな状況での方向微分の定義を以下に示します。関数は十分に滑らかで、微分を取ることができるものと仮定します。
ベクトルのスカラー値関数の微分
f ( v ) をベクトルvの実数値関数とする。すると、 f ( v ) のvに関する(またはvにおける)微分は、任意のベクトルuとのドット積によって定義されるベクトルとなる。
すべてのベクトルuについて。上記のドット積はスカラーを生成し、u が単位ベクトルであれば、u方向のvにおけるfの方向微分を生成します。
プロパティ:
- もしそうなら


- もしそうなら


- もしそうなら


ベクトルのベクトル値関数の微分
f ( v ) をベクトルvのベクトル値関数とする。すると、 f ( v ) のvに関する(またはvにおける)微分は、任意のベクトルuとのドット積で定義される2次テンソルとなる。
すべてのベクトルuに対して、上記のドット積はベクトルを生成し、u が単位ベクトルであれば、方向uにおけるvでのfの方向微分が得られます。
プロパティ:
- もしそうなら


- もしそうなら


- もしそうなら


2階テンソルのスカラー値関数の微分
を2次テンソル の実数値関数とします。 の(または における)方向の微分は、
すべての2次テンソル に対してとして定義される
2次テンソルです。






![{\displaystyle {\frac {\partial f}{\partial {\boldsymbol {S}}}}:{\boldsymbol {T}}=Df({\boldsymbol {S}})[{\boldsymbol {T}}]=\left[{\frac {d}{d\alpha }}~f({\boldsymbol {S}}+\alpha ~{\boldsymbol {T}})\right]_{\alpha =0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3916392ebeebca9ebddb31e76edb1c9f8c08b586)

プロパティ:
- もしそうなら


- もしそうなら


- もしそうなら


2階テンソルのテンソル値関数の微分
を2次テンソル の2次テンソル値関数とします。 の(または における)方向の微分は、
すべての2次テンソル に対してとして定義される
4次テンソルです。






![{\displaystyle {\frac {\partial {\boldsymbol {F}}}{\partial {\boldsymbol {S}}}}:{\boldsymbol {T}}=D{\boldsymbol {F}}({\boldsymbol {S}})[{\boldsymbol {T}}]=\left[{\frac {d}{d\alpha }}~{\boldsymbol {F}}({\boldsymbol {S}}+\alpha ~{\boldsymbol {T}})\right]_{\alpha =0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e2072ee23ca460f6f32851f397b5879a47033c3b)

プロパティ:
- もしそうなら


- もしそうなら


- もしそうなら


- もしそうなら


テンソル場の勾配
任意の定数ベクトルcの方向のテンソル場の勾配 、は次のように定義されます。
次数n
のテンソル場の勾配は、次数n +1のテンソル場です。



デカルト座標
が直交座標系の基底ベクトルで、点の座標が()で示される場合、テンソル場の勾配は次のように表される。



基底ベクトルは直交座標系では変化しないため、スカラー場、ベクトル場v、および2次テンソル場の勾配には次の関係が成り立ちます。


曲線座標
が曲線座標系の反変基底ベクトルであり、点の座標が()で示される場合、テンソル場の勾配は次のように表される(証明については
[3]を参照)。


この定義から、スカラー場、ベクトル場v、および2次テンソル場の勾配について次の関係が得られます。


ここでクリストッフェル記号 は次のように定義される。

円筒極座標
円筒座標では、勾配は次のように表される。
テンソル場の発散
テンソル場の発散は再帰関係を用いて定義される。

ここで、cは任意の定数ベクトル、v はベクトル場です。がn > 1 の位数のテンソル場である場合、場の発散はn − 1 位数のテンソルになります。

デカルト座標
デカルト座標系では、ベクトル場vと 2 次テンソル場の間に次の関係が成り立ちます。

ここで、右端の式では偏微分に対する
テンソル指数表記が使用されている。
対称2次テンソルの場合、発散は次のように表記されることが多い[4]。
上記の式は、直交座標成分形式( と表記されることも多い)における
の定義として使用されることがあります
。このような定義はこの記事の残りの部分と一致しないことに注意してください(曲線座標のセクションを参照)。


違いは、微分が の行に関して実行されるか、列に関して実行されるかによって生じ、 は慣例的です。これは例で示されます。直交座標系では、2 次テンソル (行列)はベクトル関数 の勾配です。



最後の式は、代替定義/解釈[4]と同等である。
曲線座標
曲線座標では、ベクトル場vと2次テンソル場の発散は

より一般的には、
円筒極座標
円筒極座標では
テンソル場の回転
n > 1次テンソル場の回転も、cが任意の定数ベクトル、v がベクトル場
である再帰関係を使用して定義されます
。

1階テンソル(ベクトル)場の回転
ベクトル場vと任意の定数ベクトルc を考えます。インデックス表記では、外積は次のように表されます。
ここでは置換記号、別名レヴィ・チヴィタ記号です。したがって
、



2次テンソル場の回転
2階テンソルの場合
、1階テンソル場の回転の定義を用いると
、


テンソル場の回転に関する恒等式
テンソル場の回転に関する最もよく使われる恒等式は、この
恒等式はあらゆる次数のテンソル場に対して成り立つ。2次テンソルの重要なケースである、この恒等式は、



2階テンソルの行列式の微分
2階テンソルの行列式の導関数は次のように表される。

正規直交基底では、 の成分は行列Aとして表すことができます。その場合、右辺は行列の補因子に対応します。

2階テンソルの不変量の微分
2階テンソルの主な不変量は
これら3つの不変量の微分は

証拠
行列式の導関数から次のことが分かる。
他の2つの不変量の導関数については、特性方程式に戻ってみましょう。
テンソルの行列式と同じアプローチを用いると、次のことが示せる。
左辺は次のように展開できる。
したがって
、
![{\displaystyle {\frac {\partial I_{1}}{\partial {\boldsymbol {A}}}}~\lambda ^{2}+{\frac {\partial I_{2}}{\partial {\boldsymbol {A}}}}~\lambda +{\frac {\partial I_{3}}{\partial {\boldsymbol {A}}}}=\det(\lambda ~{\boldsymbol {\mathit {1}}}+{\boldsymbol {A}})~\left[(\lambda ~{\boldsymbol {\mathit {1}}}+{\boldsymbol {A}})^{-1}\right]^{\textsf {T}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3993a79ec1f95da1b9e188243300861ee7ee2e45)
右辺を展開し、左辺の項を分離すると、
または、
とを定義すると、上記は次のように書ける。


λの様々なべき乗を含む項を集めると、
そしてλの任意性を考慮すると、
これは、
2階単位テンソルの微分
を2階恒等テンソルとします。このテンソルの2階テンソルに対する微分は次のように与えられます。
これは、が に依存しない
からです。





2階テンソルのそれ自身に関する微分
を2階テンソルと
します。すると
したがって、
これは4次の恒等テンソルである。正規直交基底に関する指数表記では

この結果は
、

したがって、テンソルが対称であれば、導関数も対称となり、
対称4次単位テンソルは次の式で表される。


2階テンソルの逆関数の微分
とを2つの2階テンソルと
すると、
正規直交基底に関する指数表記では
次
の式も成り立ちます。
指数表記で
は、テンソルが対称
であれば、






証拠
思い出してください
なので、次のように書くことができる。

2次テンソルの積の法則を使う
あるいは
、

したがって、
部品ごとの統合
ドメイン、その境界、および外向きの単位法線

連続体力学におけるテンソル微分に関連するもう一つの重要な演算は部分積分である。部分積分の式は次のように書ける。
ここで、およびは任意の位数の微分可能なテンソル場であり、はテンソル場が定義されている領域に対する単位外向きの法線であり、一般化されたテンソル積演算子を表し、は一般化された勾配演算子である。が恒等テンソルに等しいとき、発散定理が得られる。




部分積分の式は、直交座標の指数表記で次のように表すことができます。
テンソル積演算が1つのインデックスの縮約であり、勾配演算が発散であり、とが両方とも2階テンソルである特別なケースでは、


インデックス表記では、
参照
参考文献
- ^ JC Simo と TJR Hughes、1998 年、「Computational Inelasticity」、Springer
- ^ JE Marsden と TJR Hughes、2000 年、「弾性の数学的基礎」、ドーバー。
- ^ RW Ogden、2000、「非線形弾性変形」、ドーバー。
- ^ ab Hjelmstad, Keith (2004).構造力学の基礎. Springer Science & Business Media. p. 45. ISBN 9780387233307。