第二章-有限差分法.pptVIP

第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法由Taylor级数展开，可得令：第二章有限差分法讨论1.只要初值问题的解对时间的二阶导数和对空间的四阶导数均有界，该差分格式就是收敛的2.?0.5？第二章有限差分法2.差分格式的稳定性定义：对于定解问题的某个差分格式，设初始层上引入误差，是第n层上的误差，如果存在常数K使得是一种范数，比如模等第二章有限差分法该格式在小于1时是稳定的第二章有限差分法§2.5差分解法在力学中的应用举例1.差分法求解梁的弯曲问题梁的弯曲方程差分格式第二章有限差分法例1.等截面超静定梁如图所示，左端绞支右端固支。在绞支端作用一力偶，梁的长度L和截面抗弯刚度EI已知，求梁的挠度和绞支端的转角.第二章有限差分法解：将梁四等分，则步长h＝L/4节点1，2，3的差分方程边界条件的差分格式虚拟节点第二章有限差分法精确解差分解误差最大挠度?13h2/(22EI)5%转角L/(4EI)0.1477L/(EI)40%转角误差大，怎么办？第二章有限差分法例2.文克尔弹性地基上的梁如图所示，已知梁的长度L＝6.096m，地基的等效弹性系数k＝5×104kNm2，抗弯刚度EI＝2.6336×106kNm2，梁上作用均布载荷q＝14.88kN/m。求梁的挠度、弯矩及其地基反力.解：由文克尔假设，地基反力为第二章有限差分法因此弹性地基上的梁的挠曲方程为将地基梁四等分，则h＝L/4，对应的差分方程第二章有限差分法节点1，2，3的差分方程边界条件的差分格式虚拟节点第二章有限差分法节点1、2、3的弯矩为地基反力A端处的支座反力第二章有限差分法2.差分法求解薄板的弯曲问题弹性薄板理论从几何特征上讲，薄板是指tb的板第二章有限差分法弹性薄板小挠度问题常用基本假设1）中法线始终垂直于中面；2）垂直于板中面的应力很小，一般忽略不计；3）弯曲后板的中面无应变.第二章有限差分法根据板的几何变形特征应变与位移间的关系薄板本构关系薄板面内力与位移间的关系第二章有限差分法薄板控制方程边界条件？第二章有限差分法边界条件简支固支第二章有限差分法边界条件自由边滑动边第二章有限差分法2.薄板弯曲问题的有限差分法第二章有限差分法第二章有限差分法薄板弯曲控制方程的差分格式每一个节点的挠度需要涉及13个节点的挠度每一个节点的挠度需要涉及5个节点的挠度和相应的“弯矩和”第二章有限差分法第二章有限差分法例1.如图所示的简支正方形板，承受垂直于板面的横向均布载荷q作用，求板的最大挠度。解：将薄板在x方向和y方向分别四等份，则h=k=a/4，应用弯矩控制方程第二章有限差分法同理，可得挠度的控制方程最大挠度为节点3处的挠度，与精确解0.00406qa4/D,其相对误差为0.79%。如果想减小误差，增加网格数即可，比如将网格数增大到64个.第二章有限差分法例2.边长为a的固定的方形薄板如图所示，承受横向均布载荷q.试用差分法确定各节点的挠度和弯矩（h=a/4）解：第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法第二章有限差分法?定解条件差分格式表示?整理微分方程与定解条件，将微分方程转化为代数方程组其中：第二章有限差分法例1用差分法求解边值问题：解：取步长h＝0.1，节点xn=0.1n(n=0,1,…,10)第二章有限差分法数值计算结果：解析结果：结果比较：第二章有限差分法3.偏微分方程的差分格式构造与求解二阶偏微分方程的分类利用这个判别式，可以判定波动方程、弦振动方程隶属于双曲型方程，静电场、静磁场方程隶属于椭圆型方程，热传导方程隶属于抛物型方程。第二章有限差分法一