变步长差分法求解薄板小挠度弯曲问题.pdfVIP

一 第 12卷 第 3期 漳州师范学院学报 (自然科学版) Vo1．12No．3 1999年 7月 JournalofZhangzhouTeachersCollege (Nat．Sci．) July 1999 变步长差分法求解薄板小挠度弯曲问题 []；杉 ，f 杨 文元 (漳州农业 槭【 学柱 ，潭州 363000) 摘要；本文舟绍斜弯薄板弹性阶段小挠度变形的也移、应力的一种数值解法一一盎步长差分法，为瞻蛙精确较快得到 问题的 解 ，文 中推导了任意变步长网括划分的差分方程及越界条件 ． 求解薄板小挠度弯 曲问厦。 关链词：小挠度 ，差分方程 ，变步长，变步长差分方程，曲界条件 0引言 j：反， _ 靛藿步苏，j芋性誓， 薄板在自重．和外力作用下产生的小挠度变形，所引起的 力是在弹性阶段范围内， 可用弹性力学进行分析，建立弹性力学基本方程一般均为微分方程，而且边界条件也是 微分方程，对于工程上许多问题，并不能够得出函数式精确解答，基本都是采用各种数 值解法，差分法是其中一种常用解法，本文拟用变步长差分法来分析斜弯薄板小挠度弯 曲的变形。 1 等间距网格的差分方程 差分法是将弹性力学基本微分方程和边界条件微分方程用差分方程这种代数方程来近 似表示，将求解微分方程的问题转变为求解代数方程的问题 差分方程的建立首先需要 导出弹性力学上常用的差分公式，进而推导差分方程 ，[1]给 出如下的推导过程： 如图 l所示，在弹性体上用平行于坐标轴 z、y的两组平行线组成的网格进行划分，各网 格之间的间距都相等为 h。设 f(x，Y)为弹性体内某一连续函数，它可以是某一位移分量 或应力分量，也可 温度函数，f(x，y)在平行于 x轴的任一根网线上的值只随 x坐标改 变而变化，在 y向上也同样 。设在结点 0处，函数 f (x，Y)的值为 fo，则邻近结点 0处 函数 f(x，y)展开的台劳级数 的形式如下； X X 12 8 4 5 8 4 5 儿 3 0 1 l1 3 0 1 9 7 2 6 7 2 6 10 10 ry y 图2不等间距网格划分 本文1999-06·O8收到 第 3期 漳州师范学院学报 (自然科学版) ，=，。+(誓)。。一z。)+1azf)。一z。)+击()。。一 对结点 3及 l有 ：zj= z0一 h，z1= ￡04-^所 以 ￡3一 zo一一 ^，￡ 代入 上面台劳级数展开式得 ： ，一，。一(薯)。+萼()。一等()。+ ()_．-… ，一，。+^(薯)。+萼()。一譬(器)。+24r坐axe、+…… 一 般 ^取值较小，可以略去三次幂及 以上的各项，联立求解得差分格式 ： ( = 及 (a,f 同理可得：(薯)。= 及()。= 二垫 这是基本差分公式，可以利用它们来推出(菩)。及()。