13、材料力学中的有限元应用：两端固定梁的欧拉 - 伯努利理论分析

材料力学中的有限元应用：两端固定梁的欧拉 - 伯努利理论分析

1. 问题描述

考虑一根长度为 $L$、横截面积为 $S$、密度为 $\rho$ 的均质梁 $\Omega$，其力学行为为小扰动下的各向同性弹性。梁的材料具有杨氏模量 $E$ 和泊松比 $\nu$。梁的轴线为 $(O;X_1)$，当前横坐标为 $x$，梁的横截面由几何变量 $(y,z)$ 参数化。梁在 $x = 0$ 和 $x = L$ 处两端固定，且沿梁施加横向力的体积密度 $f = -f_2X_2$。

2. 位移变分公式与虚功原理 - 理论部分

2.1 虚位移场结构

选择欧拉 - 伯努利运动学框架来描述实际位移场，虚位移场 $U^ $ 表示为：
$U^ = -y \dot{u}^ (x)X_1 + u^ (x)X_2$
满足边界条件：
$u^ (0) = u^ (L) = 0$
$\dot{u}^ (0) = \dot{u}^ (L) = 0$
其中 $\dot{u}^ \equiv \frac{du^ }{dx}$。$U^*$ 表示满足上述三个条件的所有位移场的集合。

2.2 虚功定义

内力虚功 $T_{Int}^ $ 和外力虚功 $T_{Ext}^ $ 定义如下：
$T_{Int}^ = -\int_{\Omega} \sigma_{ij}\epsilon_{ij}(U^ )d\Omega$
$T_{Ext}^