有限单元法初步 有限单元法是在矩阵位移法基础上发展起来的一种结构 分析方法,用于板壳、实体等结构的分析。 有限元分析的步骤与矩阵位移法基本相同,过程也相似。 离散化: 水坝 单元分析 整体分析: 求应力:
有限单元法初步 有限单元法是在矩阵位移法基础上发展起来的一种结构 分析方法，用于板壳、实体等结构的分析。 有限元分析的步骤与矩阵位移法基本相同，过程也相似。 1 2 4 6 3 5 离散化： 水坝 单元分析： 整体分析： 求应力：
§1杆系结构的有限单元法 §11泛函与变分 “最速落径问题”一-质量为m的小环从A处自由滑下, 试选择一条曲线使所需时间最短。(不计摩擦) 设路径为y=y(x) a A ds=vdx+dy y hgh Y B 2 称T为y(x)的泛函, d√29 y(x)为自变函数。 即以函数作自变量以积 所需时间T[y(x)]= 2dx分形式定义的函数为泛函
§1 杆系结构的有限单元法 §1.1 泛函与变分 “最速落径问题”---质量为m的小环从A处自由滑下， 试选择一条曲线使所需时间最短。（不计摩擦） A B X Y 设路径为 y=y（x） 2 2 ds = dx + dy dx dt y dt ds v − = = 2 1 v = 2gh dx gh y dt − = 2 1 2 − = a dx gh y T y x 0 2 2 1 所需时间 [ ( )] a y 称T为y（x）的泛函， y（x）为自变函数。 即以函数作自变量以积 分形式定义的函数为泛函
§11泛函与变分 y(x)=y(x)+(x) 称S(x)为y(x)的变分,它是一个无穷小的任意函数。 变分运算在形式上与微分运算相同。 y2(x)=2y(x)(x) x x+dx 微分与变分运算次序可以交换。 A (oy)= y=y(x 积分与变分运算次序也可以交换。 Y Sy(x) ∫八x,y(x)x=」/(x,y(x) 5(x)
§1.1 泛函与变分 A X Y ( ) ( ) ( ) * y x = y x +y x ( ) 2 ( ) ( ) 2 y x = y x y x 变分运算在形式上与微分运算相同。 y=y（x） x+dx dy x ( ) * * y = y x 称 y(x) 为y（x）的变分，它是一个无穷小的任意函数。 y(x) 微分与变分运算次序可以交换。 ( ) ( ) dx dy y dx d = 积分与变分运算次序也可以交换。 = 2 1 2 1 [ , ( )] [ ( , ( ))] x x x x f x y x dx f x y x dx
§12变形体虚位移原理 外力虚功 (xd (x)50(x)d 平衡位置 内力虚功 SW=[M(x)ok(x)+O()8y+ N(x)SE]dx 虚功方程 q(x)ay(x)dx=L[M(x)Sk(x)+O(x)8y+N(x)8Eldx §13势能原理 1.应变能 弯曲应变能V=P△/2=Mkt 工△ 拉压应变能V=P△/2 Ndx 2 △ 剪切应变能V=P△/2=2Qr P 工△
§1.2 变形体虚位移原理 = l We q x y x dx 0 ( ) ( ) 外力虚功 = + + l Wi M x k x Q x N x dx 0 [ ( ) ( ) ( ) ( ) ] 内力虚功 虚功方程 Wi = We = + + l l q x y x dx M x k x Q x N x dx 0 0 ( ) ( ) [ ( ) ( ) ( ) ( ) ] §1.3 势能原理 1.应变能 弯曲应变能 P Ve = P/ 2 = l M dx 2 0 1 拉压应变能 Ve = P/ 2 = l N dx 2 0 1 P P 剪切应变能 Ve = P/ 2 = l Q dx 2 0 1 y（x） 平衡位置 q(x) y
§12变形体虚位移原理 虚功方程 9)(k=(x以()+x8y+2(y③) 平衡位置 §13势能原理 1.应变能 弯曲应变能V=P△2=2M △ 拉压应变能v=PA/2=2 Nadx 剪切应变能V=P△2=2grh API 2.外力势能 外力从变形状态退回到无位移的 原始状态中所作的功 P21P3 V:=∑2A=x △ 3.结构势能 Ep=Ve+Vp TIMK+Or+ NEldx-n.ydx yx 0
2.外力势能 §1.3 势能原理 1.应变能 弯曲应变能 P P Ve = P/ 2 = l M dx 2 0 1 拉压应变能 Ve = P/ 2 = l N dx 2 0 1 P 剪切应变能 Ve = P/ 2 = l Q dx 2 0 1 1 2 3 P1 P2 P3 外力从变形状态退回到无位移的 原始状态中所作的功. Ve = −Pi i * y(x) q(x) = − l Ve q x y x dx 0 * ( ) ( ) 3.结构势能 * EP =Ve +VP §1.2 变形体虚位移原理 虚功方程 = + + l l q x y x dx M x k x Q x N x dx 0 0 ( ) ( ) [ ( ) ( ) ( ) ( ) ] y（x） 平衡位置 q(x) y = + + − l l M Q N dx qydx 0 0 [ ] 2 1