本文主要是介绍有限元变分形式,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
有如下的微分方程
a d 2 u d x 2 + b = 0 , 0 ≤ x ≤ 2 L u ∣ x = 0 = 0 a d u d x ∣ x = 2 L = R \begin{aligned} & a\frac{d^2u}{dx^2}+b=0, \quad 0 \leq x \leq 2L \\ & \left. u \right|_{x=0}=0 \\ & a\left. \frac{du}{dx} \right|_{x=2L}=R \end{aligned} adx2d2u+b=0,0≤x≤2Lu∣x=0=0adxdu∣∣∣∣x=2L=R
其中 a = E A a=EA a=EA,该微分方程在固体力学中拥有具体的物理意义。它描述了由弹性模量为 E E E的材料制成的截面积为 A A A的一维杆,在右端承受分布载荷 b b b和集中载荷 R R R的张力,如图1.3所示。
此问题可以用最小化势能泛函 Π \Pi Π来表示
Π = ∫ L 1 2 a ( d u d x ) 2 d x − ∫ L b u d x − R u ∣ x = 2 L u ∣ x = 0 = 0 \begin{aligned} & \Pi=\int_L\frac{1}{2}a \left(\frac{du}{dx}\right)^2dx-\int_Lbudx-\left.Ru\right|_{x=2L}\\ & \left. u \right|_{x=0}=0 \end{aligned} Π=∫L21a(dxdu)2dx−∫Lbudx−Ru∣x=2Lu∣x=0=0
我们同样使用上一节中使用的线性拉格朗日插值基函数来表示 { u } \{u\} {u},则可以写出第二个单元上的势能的值
Π e
这篇关于有限元变分形式的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!