變分方法
對連續介質問題，位置函式的\(u\)的泛函為
\[\Pi=\int_{\Omega}\boldsymbol{F}\Bigg(u,\frac{\partial u}{\partial x},\cdots\Bigg)\mathrm{d}\boldsymbol{\Omega}+\int_{\Gamma}\boldsymbol{E}\Bigg(u,\frac{\partial u}{\partial x},\cdots\Bigg)\mathrm{d}\Gamma
\]
其中，F和E是特定的運算元，\(\Omega\) 是求解域,\(\Gamma\) 是 \(\Omega\) 的邊界.
在變分方法中，連續介質問題的解u是使得泛函\(\Pi\) 對於微小變化的\(\delta u\)取駐值，即泛函的變分等於0。
\[\delta\Pi=0
\]
對於可以執行變分原理的問題，
對於可以運用變分原理的問題，可以建立其得到近似解的如下方法. 未知函式的近似解可表示成帶有待定引數的試探函式。
\[u\approx\tilde{u}=\sum_{i=1}^{n}N_{i}\alpha_{i}=N\alpha
\]
式中，\(\alpha_i\)為待定引數，\(N_i\)是已知的函式序列.將式(1-18)代入式(1-16),得到用試探函式\(\tilde{u}\) 和待定引數\(\alpha\) 表示的泛函\(\Pi\).泛函的變分為零相當於將泛函對關於待定引數進行全微分，並令其等於0。即
\[\delta\prod=\frac{\partial\prod}{\partial\alpha_{1}}\delta\alpha_{1}+\frac{\partial\prod}{\partial\alpha_{2}}\delta\alpha_{2}+\cdots+\frac{\partial\prod}{\partial\alpha_{n}}\delta\alpha_{n}=0
\]
由於\(\delta\alpha_1,\delta\alpha_2,\cdotp\cdotp\cdotp,\delta\alpha_n\)是任意的,式(1-19)成立時必有\(\frac{\partial\Pi}{\partial\alpha_1},\frac{\partial\Pi}{\partial\alpha_2},\cdotp\cdotp\cdotp,\frac{\partial\Pi}{\partial\alpha_n}\)都等於零,因而有
\[\frac{\partial\Pi}{\partial\alpha}=\begin{bmatrix}\frac{\partial\Pi}{\partial\alpha_1}\\\frac{\partial\Pi}{\partial\alpha_2}\\\vdots\\\frac{\partial\Pi}{\partial\alpha_n}\end{bmatrix}=\mathbf{0}
\]
由上述與待定引數的數目相等的方程組可求出\(\alpha\).
如果在泛函\(\Pi\)\中u 及其導數的最高階次為二階,則稱泛函 \(\Pi\)為二次泛函,工程中的許多問題都屬於二次泛函.對於二次泛函問題,式(1-20)退化為一線性方程組
\[\frac{\partial\Pi}{\partial\alpha}=K\alpha-P=0
\]
對上式進行變分，得
\[
\delta\biggl(\frac{\partial\Pi}{\partial\alpha}\biggr)=\begin{bmatrix}\frac{\partial}{\partial\alpha_1}\biggl(\frac{\partial\Pi}{\partial\alpha_1}\biggr)\delta\alpha_1+\frac{\partial}{\partial\alpha_2}\biggl(\frac{\partial\Pi}{\partial\alpha_1}\biggr)\delta\alpha_2+\cdots\\\vdots\\\frac{\partial}{\partial\alpha_1}\biggl(\frac{\partial\Pi}{\partial\alpha_n}\biggr)\delta\alpha_1+\frac{\partial}{\partial\alpha_2}\biggl(\frac{\partial\Pi}{\partial\alpha_n}\biggr)\delta\alpha_2+\cdots\end{bmatrix}=K\delta\alpha
\]
由矩陣K的子矩陣
\[K_{ij}=\frac{\partial^{2}\prod}{\partial\alpha_{i}\partial\alpha_{j}} ,\quad K_{ji}=\frac{\partial^{2}\prod}{\partial\alpha_{j}\partial\alpha_{i}}
\]
可知
\[K_{ij}=K_{ji}^{\mathrm{T}}
\]
式\(\frac{\partial\Pi}{\partial\alpha}=K\alpha-P=0\)可以近似泛函表示成
\[\Pi=\frac{1}{2}\alpha^{\mathrm{T}}K\alpha-\alpha^{\mathrm{T}}P
\]
對\(\Pi=\frac{1}{2}\alpha^{\mathrm{T}}K\alpha-\alpha^{\mathrm{T}}P\)進行變分，得
\[\delta \Pi=\frac{1}{2}\delta\boldsymbol{\alpha}^{\mathrm{T}}\boldsymbol{K}\boldsymbol{\alpha}+\frac{1}{2}\boldsymbol{\alpha}^{\mathrm{T}}\boldsymbol{K}\delta\boldsymbol{\alpha}-\delta\boldsymbol{\alpha}^{\mathrm{T}}\boldsymbol{P}
\]
由於K的對稱性，有
\[\delta\alpha^{\intercal}K\alpha=\alpha^{\intercal}K\delta\alpha
\]
因而
\[\delta\Pi=\delta\alpha^{\mathrm{T}}(K\alpha-P)=0
\]
因為\(\delta\boldsymbol{\alpha}\) 是任意的,可得 \(\boldsymbol{K\alpha}-\boldsymbol{P}=0\),即為式(1-25)
例子 對於問題
\[\begin{cases}\frac{\partial^2u}{\partial x^2}+u=-x&(0\leqslant x\leqslant1)\\u(0)=0,&u(1)=0\end{cases}
\]
建立變分原理.
該問題變分可表示為
\[\delta\Pi=\int_0^1\biggl(-\frac{\mathrm{d}^2u}{\mathrm{d}x^2}-u-x\biggr)\delta u\mathrm{d}x+\biggl[\frac{\mathrm{d}u}{\mathrm{d}x}\delta u\biggr]_0^1
\]
對方程中積分的第一項進行分佈積分，得
\[\delta\Pi=\int_0^1\biggl(\frac{\mathrm{d}u}{\mathrm{d}x}\frac{\mathrm{d}(\delta u)}{\mathrm{d}x}-u\delta u-x\delta u\biggr)\mathrm{d}x
\]
利用變分運算元的交換性
\[\delta \Pi=\delta\int_{0}^{1}\Bigg(\frac{1}{2}\Bigg(\frac{\mathrm{d}u}{\mathrm{d}x}\Bigg)^{2}-\frac{1}{2}u^{2}-xu\Bigg)\mathrm{d}x
\]
於是，可得泛函
\[\Pi=\int_{0}^{1}\Bigg(\frac{1}{2}\Bigg(\frac{\mathrm{d}u}{\mathrm{d}x}\Bigg)^{2}-\frac{1}{2}u^{2}-xu\Bigg)\mathrm{d}x
\]
參考文獻
Matlab有限元結構動力學分析與工程應用-徐斌