hzf第5章用差分法和变分法解平面问题.pptVIP

例题 例题2 用差分法计算图中A和B点的应力分量。 解：为反映对称性，取A为基点。令边界点的应力函数值：边界点的导数值：由上式及，求出边界外一行虚结点的 值：对1点列差分方程： 代入各 值，解出。再求出应力分量：例题3正方形的板块，厚度，受一对集中力F的作用，如图。试 取，应用差分法求解该问题的应力分量。 解：⑴本题具有的两个对称轴，为了反映对称性，在 y 向外荷载作用下，取网格结点编号如图所示。求出边界上各结点的值，如下图所示。 结点 A B CDEGHI者可检验，上述的值反映了边界结点 和边界外一行虚结点上 值的对称性。由得到由得到 ⑷ 对内结点1、2、3、4分别列出下列类型的方程： 0点： ⑸按照应力公式 及，求得AJ及EI截面上的应力分量：例题4试证明，在同样的应变分量 ， 和下，平面应变情况下单位厚度的形变 势能大于平面应力情况下的形变势能。 代入，得 显然，方括号内 将式⑴中的 ， 都作为式(b)的变换，整理后得平面应变情况下的形变势能公式，从式⑶可见，在平面应变情况下，形变势能中的第一、二、三项均大于平面应力情况下的值，而第四项不变。因此，平面应变的形变势能 大于平面应力的形变势能U 。 例题5 图中表示一板块，受到铅直方向均布拉力作用下发生拉伸变形，并使之两端固定下来，若在其中切开一小口AB时，试说明板的形变势能将发生什么变化？ 解： ⑴当AB线切开时，AB线上的应力趋于 0。而形变势能是正定的，，当这部 应力时，相应的形变势能也失去因 此，板的总的形变势能减少。⑵ 当AB线切开后，边界CD和EF仍是固 定的，我们可以比较两种状态： (b) AB线张开，出现裂纹。这是稳定的平衡状态。由于系统的稳定平衡状态与邻近的状态相比，总势能处于极小值，而(a)、(b)两种状态的外力势能不变，因此，(b)的形变势能小于(a)，即形变势能将减少。 例题6 单位厚度的深梁，两侧边固定，上下边受均布荷载q作用，如图所示。试用位移变分法求解其位移。 （取， 并设）。上式已满足两端的约束边界条件， 以及对称和反对称性条件。以下按瑞利- 里茨法进行计算。 假设只取u，v中一项，即 将u和v代入形变势能公式（平面应力问题），得：在本题中体力，在边界上只有的均布荷载，。由此，瑞利-里茨方程成为当取，且时，上两式方程简化为 由此解出, 位移分量的解答是 上式已反映了位移对称于y轴的要求：v为x的 偶函数，u为x的奇函数。仅取各一项进行运算， 由于体力，面力只存在于AB边 （　　　），因此求解的位移变分方程 为： 当，且取泊松系数时，形变 势能简化为 将u、v 代入, 从式(c)、 (d)解出 F/2 F/2 F/2 -Fh/2 -Fh/2 -Fh ⑶ 计算边界外一行结点的 值。 对结点1， 对结点2， 对结点3， 对结点4， 解出 例题 对于平面应变情况，只需将上式中 ， 变换为 解：平面应力情况下，单位厚度的形变势能是： 例题 (a) 例题 (c) 例题 例题 C D E F A B 例题 例题 (a) AB切开后，仍然处于闭合状态，不发生张开。这是不稳定的平衡状态； 例题 q y x b u v b a a o q 解：在图示荷载作用下，深梁的位移应对称于x轴，而反对称于y轴。因此，位移分量u应为 、 的奇函数，而v为 x 、y 的偶函数， x y 如图所示。可以设定位移势函数如下： 例题 例题 例题 再积分求U，边界是，且，从到积分。再将U代入上式，得到两个求的方程： 例题7 图中所示的薄板，厚度，三边固定，一边受到均布压力q的作用。试用瑞利-里茨的位移变分法求解，其中取，。 例题 a a bx y q 解：在瑞利-里茨法中， 设定位移试函数应满 足位移边界条件，并 应反映图示问题的对称性。取 例题 例题 （a） （b） 形变势能U为 将U及代入式(a)，(b)，得 (c) (d) 例题 于是得到位移分量， 再求应力分量，取，得： 例题 在对称轴上，x=0，, 在边界，,本题中，由于u，v中各只取一项，且取，因此，求出的位移解的精度较低；而由近似解的位移求应力时，其应力精度要降低一阶，其精度更差些。对于实际问题，应取更多的项数进行计算。 将几何方程代入，形变势能还可用位移分量来表示 一 变分及其性质高等数学我们学过微分的概念，微分是变量的增量。那么什么是变分呢？变分是函数的增量，通常用δ表示。变分具有以下的性质： §5-5 位移变分方程 二 位移变分方程设弹性体在一定外力作用下，处于平衡状态，发生的