第五章 有限差分发和变分法.pptx

2016-11-18第五章用差分法和变分法解平面问题5－1差分公式的推导5－1差分公式的推导差分法：是微分方程的近似解法，具体的讲，差分法就是把微分用差分来代替，把导数用差分商来代替，从而把基本方程和边界条件（微分方程）近似用差分方程来表示，把求解微分方程的问题变成求解代数方程问题。差分法的数学基础：泰勒公式5－1差分公式的推导5－1差分公式的推导5－1差分公式的推导注（5－1）－－（5－4）是最基本的差分公式5－1差分公式的推导混合二阶导数的差分公式四阶导数的差分公式5－1差分公式的推导讨论：（1）差分公式是微分方程在数学上的近似；（2）在推导（5－1）－－（5－4）时，略去了三次项及更高阶项；（4）要想求差分解，前提是要有微分方程。5－2应力函数的差分解5－2应力函数的差分解当不计体力时，我们已把弹性力学平面问题归结为在给定边界条件下求解双调和方程的问题。用差分法解平面问题，就应先将双调和方程变换为差分方程，而后求解之。5－2应力函数的差分解1、应力分量（不计体力）5－2应力函数的差分解双调和方程对于弹性体边界以内的每一结点，都可以建立这样一个差分方程。应力函数在域内应该满足上式。整理即得2、差分方程（相容方程）相容方程的差分公式问题：边界上的点（边界附近的点）怎么办？？？？？？5－2应力函数的差分解代入上式，即得：5－2应力函数的差分解由图（5－2）可见因此，式（b）可以改写成5－2应力函数的差分解约去dy、dx得：5－2应力函数的差分解由高等数学可知，将此式亦从A点到B点沿s进行积分，就得到边界上任一点B处的φ值。为此利用分部积分法，得：5－2应力函数的差分解将式(c),(d)代入，整理得：由前知，把应力函数加上一个线性函数，并不影响应力。因此，可设想把应力函数加上a+bx+cy，然后调整a，b，c三个数值，使得5－2应力函数的差分解于是式(d)，式(e)简化为：讨论：（4）以上结果不能用于多连体的情况。5－2应力函数的差分解5－2应力函数的差分解用差分法解弹性平面问题时，可按下列步骤进行：5－2应力函数的差分解（5）按照公式（5－9）计算应力的分量。说明：如果一部分边界是曲线的，或是不与坐标轴正交，则边界附近将出现不规则的内结点。对于这样的结点，差分方程（5－10）必须加以修正。5－4弹性体的变形势能和外力势能5－4弹性体的变形势能和外力势能变分法：主要是研究泛函及其极值的求解方法。泛函：函数是函数的函数；能量法：弹性力学中的变分法；形变势能与弹性体的受力次序无关，也与受力的历史无关完全由应力和变形的最终大小确定－－保守场。设弹性体在一定外力作用下，处于平衡状态，发生的真实位移为u,v,w，它们满足位移分量表示的平衡方程，并满足位移边界条件和用位移表示的应力边界条件。弹性体受力后，发生变形，外力作功，外力功转化为变形能，储存在弹性体内，单元体内的变形能为5－4弹性体的变形势能和外力势能整个弹性体内的变形能5－4弹性体的变形势能和外力势能对应于平面问题，微元的应变能（应变比能）整个弹性体内的变形能把物理方程代入微元的应变能，分别得到用应力应变表示方程5－4弹性体的变形势能和外力势能把几何方程（2－8）代入，得到用位移分量表示的微元变形势能位移分量表示的弹性体变形势能5－4弹性体的变形势能和外力势能讨论（1）变形势能是变形分量或位移分量的二次泛函，叠加原理不再适用；（2）变形或位移发生时，变形势能总是正的；5－4弹性体的变形势能和外力势能在静态或准静态时，外力的势能转化成外力的功，因此弹性体的外力势能5－5位移变分方程5－5位移变分方程对于三维时：一、位移变分方程（拉格朗日变分方程）注：变分和微分都是微量，运算方法相同。5－5位移变分方程给出弹性体的限制条件：（1）没有温度改变（热能没变）；（2）没有速度改变（动能没变）。根据能量守恒，变形势能的增加等于外力势能的减少（外力的虚功）三维：上式：位移变分方程（拉格朗日变分方程）5－5位移变分方程二、虚功方程按照变分原理，变分运算与定积分的运算可以交换次序。利用（5－15）代入位移变分方程5－5位移变分方程对应于二维情况（5－24）就是虚功方程，表示：如果在虚位移发生前，弹性体是处于平衡状态，那么，在虚位移过程中，外力在虚位移上所做的虚功，等于应力在虚应变上所做的虚功。5－5位移变分方程三、极小势能原理令在虚位移过程中，外力的大小和方向保持不变，只是作用点发生了改变将变分与定积分交换次序，移项5－5位移变分方程在给定外力作用下，在满足位移边界条件的所有各组位移中间，实际存在的一组位移应使总势能成为极值，对于稳定平衡状态，这个值是极小值。位移变分方程（极小势能原理或虚功方程）等价于平衡微分方程和应力边界条件。5－6位移变分法5－6位移变分法位移变分法：（1）设定一组包含若干待