最小ポテンシャルエネルギの原理

内部仮想仕事\(W_i\)は, 仮想変位\(\delta u_i\)によるひずみエネルギの増加 \(\delta U\) と考えられる.

$\displaystyle \delta U (=\delta W_i)$ $\textstyle =$ $\displaystyle \int_V \sigma_{ij} \delta \epsilon_{ij} dV$  
  $\textstyle =$ $\displaystyle \int_V\frac{\partial U_0}{\partial \epsilon_{ij}} \delta \epsilon_{ij} dV
=\int_V\delta U_0 dV = \delta \int_V U_0 dV$ (331)

物体力や表面力が保存力 (conservative force)であり, ポテンシャル \(\Phi_{F}\)\(\Phi_{T}\)であるとすると,
\begin{displaymath}
F_i=-\frac{\partial \Phi_{F}}{\partial u_i}, \qquad
\stackrel{n}{T}_i=-\frac{\partial \Phi_{T}}{\partial u_i}
\end{displaymath} (332)

のように表されるから,仮想仕事の原理は,
\begin{displaymath}
\delta\int_V U_0 dV + \delta \int_V \Phi_{F} dV + \delta \int_S \Phi_{T} dS = 0
\end{displaymath} (333)

すなわち,
\begin{displaymath}
\delta \left(\int_V U_0 dV + \int_V \Phi_{F} dV + \int_S \Phi_{T} dS\right) = 0
\end{displaymath} (334)

ここで,
\begin{displaymath}
\Pi= \int_V U_0 dV + \int_V \Phi_{F} dV + \int_S \Phi_{T} dS
\end{displaymath} (335)

全ポテンシャルエネルギ (total potential energy)という. \(\Pi\)は,弾性体のひずみエネルギと物体力および表面力による ポテンシャルエネルギを加えたものである. 微小な仮想変位\(\delta u_i\)を与えても物体力と表面力が変化しないとすると,
\begin{displaymath}
\Pi= \int_V U_0 dV - \int_V F_iu_i dV - \int_S \stackrel{n}{T}_i u_i dS
\end{displaymath} (336)

と書ける.このように, ポテンシャルエネルギは,位置エネルギともよばれ,外力のなす仕事は, 外力のポテンシャルエネルギの減少を表すので外力のなす仕事に負 の符号をつけたものである.

これを用いると,

\begin{displaymath}
\delta \Pi=0
\end{displaymath} (337)

となる.つまり,物体がつりあい状態にある時に,\(\Pi\)が停留値をとることを 意味している. このことを,停留ポテンシャルエネルギの原理 (stationary principle of potential energy)という.

さて,この停留値が極小値であることを示そう. \(S_u\)\(S_t\) 上の境界条件を満足する実変位\(u_i\)に対する ポテンシャルエネルギ\(\Pi\)と 任意の仮想変位が加えられた場合の変位 \(u_i+\delta u_i\)に対する ポテンシャルエネルギ\(\Pi'\)を比較してみる.

$\displaystyle \Pi'-\Pi$ $\textstyle =$ $\displaystyle \int_V\left\{U_0(\epsilon_{ij}+\delta \epsilon_{ij}) - U_0(\epsilon_{ij})\right\}dV$  
    $\displaystyle -\int_V F_i \delta u_i dV - \int_{S_t} \stackrel{n}{T}_i \delta u_i dS$ (338)

ここで, \(U_0(\epsilon_{ij}+\delta \epsilon_{ij})\)\(\epsilon_{ij}\)のまわりに Taylor 展開すると
\begin{displaymath}
U_0(\epsilon_{ij}+\delta \epsilon_{ij})
= U_0(\epsilon_{ij})...
...psilon_{kl}}
\delta \epsilon_{ij} \delta \epsilon_{kl}
+\cdots
\end{displaymath} (339)

高次の項を無視して,
$\displaystyle \Pi'-\Pi$ $\textstyle =$ $\displaystyle \int_V
\frac{\partial U_0}{\partial \epsilon_{ij}}\delta \epsilon_{ij}dV$  
    $\displaystyle -\int_V F_i \delta u_i dV - \int_{S_t} \stackrel{n}{T}_i \delta u_i dS$  
    $\displaystyle + \int_V \frac{1}{2}
\frac{\partial^2 U_0}{\partial \epsilon_{ij}\partial \epsilon_{kl}}
\delta \epsilon_{ij} \delta \epsilon_{kl}dV$ (340)

となる.第1, 2, 3項の和は仮想仕事の原理式から0となる.

一方, \(\epsilon_{ij}=0\)のとき,\(U_0(0)=0\) であり, \(\partial U_0/\partial \epsilon_{ij}\) は,\(\sigma_{ij}\) を意味するが, 応力ひずみ関係から応力も 0となり, \(\partial U_0/\partial \epsilon_{ij}=0\) となる.したがって,

\begin{displaymath}
U_0(\delta\epsilon_{ij})
=\frac{1}{2}
\frac{\partial^2 U_0}{...
...rtial \epsilon_{kl}}
\delta \epsilon_{ij} \delta \epsilon_{kl}
\end{displaymath} (341)

となるので,
\begin{displaymath}
\Pi'-\Pi = \int_V U_0(\delta \epsilon_{ij})dV
\end{displaymath} (342)

となる. \(U_0\)はひずみエネルギを表わす正値二次形式なので,\(U_0 \ge 0\)である.よって,
\begin{displaymath}
\Pi'\ge\Pi
\end{displaymath} (343)

となる. このことから,無ひずみ状態(これを自然状態 (natural state)という) の近傍では,ポテンシャルエネルギは極小値をもつことがわかる. 後で述べるように,弾性問題では解は唯一であるという原理から, ポテンシャルエネルギの極小値は最小値であることになる. したがって,ポテンシャルエネルギの停留を示す式は 最小値を意味している. これを最小ポテンシャルエネルギの原理 (principle of minimum potential energy)という.

Akihiro Nakatani 2001-06-25