弹性力学第五章差分法.pptVIP

第一页，共二十二页，2022年，8月28日 差分法简介 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 差分法是微分方程的一种近似数值解法。它不是去寻求函数式的解答，而是寻求函数在一些网格结点上的数值。 差分法就是把微分用有限差分代替，把导数用有限差商代替，从而把基本方程和边界条件（一般均为微分方程）近似地改用差分方程（代数方程）来表示，从而把求解微分方程的问题转化成求解代数方程的问题。 第二页，共二十二页，2022年，8月28日 差分公式的推导（1） 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 设函数f为弹性体内的某一个连续函数（可以是应力函数，应力分量函数，位移函数…），将函数f在0点处沿3-0-1这条平行于x轴的直线展开： 第三页，共二十二页，2022年，8月28日 差分公式的推导（2） 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 网格间距h很小： 联立求解： 第四页，共二十二页，2022年，8月28日 差分公式的推导（3） 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 第五页，共二十二页，2022年，8月28日 差分公式的推导（4） 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 第六页，共二十二页，2022年，8月28日 应力函数的差分解（1） 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 对于边界内距边界距离大于h的结点可建立如下方程： 第七页，共二十二页，2022年，8月28日 应力函数的差分解（2） 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 对于距边界内一行的结点建立的方程须用到边界上的结点以及边界外一行上的结点的应力函数值。因此首先我们要用边界条件确定边界上的点的应力函数值。 第八页，共二十二页，2022年，8月28日 应力函数的差分解（3） s -dx dy o (xB, yB) x y ds n A B 第九页，共二十二页，2022年，8月28日 应力函数的差分解（4） s -dx dy o (xB, yB) x y ds n A B 从基点A到边界上任意点B对s积分 第十页，共二十二页，2022年，8月28日 应力函数的差分解（5） s -dx dy o (xB, yB) x y ds n A B 由于： 对应力函数在s上利用分部积分做积分： 联合下式： 得到： 第十一页，共二十二页，2022年，8月28日 应力函数的差分解（6） s -dx dy o (xB, yB) x y ds n A B 代入 得到： （1） （2） （3） 第十二页，共二十二页，2022年，8月28日 应力函数的差分解（7） s -dx dy o (xB, yB) x y ds n A B 由于在应力函数中加上一个线性函数不影响应力的解，因此我们可以假想通过在应力函数中加上一个特殊的线性函数使得应力函数在A点的值以及对x和y的一阶偏导都为零。从而使1-3式有如下简化形式： 第十三页，共二十二页，2022年，8月28日 应力函数的差分解（8） s -dx dy o (xB, yB) x y ds n A B 第一式：表示A与B之间的x方向的面力之和。 第二式：表示A与B之间的y方向的面力之和的负数。 第三式：表示A与B之间的面力对B点的力矩之和，在如右图的坐标系中，顺时针为正。 第十四页，共二十二页，2022年，8月28日 应力函数的差分解（9） s -dx dy o (xB, yB) x y ds n A B 对于多连体问题，在一条边界上可以采用以上三个公式，对另一条边界则需要通过位移单值条件求得边界上一个基点的应力函数及其对应的偏导值，然后再利用未简化的1-3式求解。 第十五页，共二十二页，2022年，8月28日 应力函数的差分解（10） 1 2 3 4 0 5 6 7 8 9 10 11 12 13 14 h h y x B A 对于边界外一行，利用下面的差分公式： 可得到： 第十六页，共二十二页，2022年，8月28日