扁壳弯曲问题的优化计算.pdfVIP

强度与环境2000增刊1 扁壳弯曲问题的优化计算 郄禄文 焦永树范慕辉 (河北大学数学系，保定，051002) (河北工业大学机槭工程学院，天津，300130) 文摘本文利用h山ce算子的单调性，论证了扁壳弯曲问题平衡方程的单调性。以一个受均 布载荷作用的正方形底面简支球面扁壳为例，采用数学规划加权残值法得到了壳体中点法向位 移的最小上界和最大下界，精确解位于其问。计算表明，取最小上界和最大下界的平均值作为 近似解，可以用较少的计算工作量得到误差可知的高精度近似解。 一，—、 关t词数学规划加权残值法扁壳位移啦下界j 一…， 1 引言 作为求解微分方程定解问题的一种半解析近似解法，加权残值法具有原理简单，方法简 便，计算工作量小和适用范围广等突出特点。但毋庸讳言，加权残值法的求解精度往往受多 种因素的制约，如待解问题的性质、试函数的完备性、权函数的选取以及求解者的经验等等。 1992年，朱宝安教授【11综合利用加权残值法应用范围广和数学规划可以优化的特点，提出 了“数学规划加权残值法(MP—MWR)”，将具有单调性的一类微分方程问题转化为数学规划 问题，通过对残差方程中的待定系数进行优化，使偏大的近似解尽量小，偏小的近似解尽量 大，从而得到近似解的最小上界和最大下界，真解总是位于这两者之间。这样，便可以用较 少的计算工作量得到误差可知的高精度近似解。 2残差关于解的单调性 数学规划加权残值法是建立在残差关于解的单调性理论基础上的一种近似解法。所谓 残差关于解的单调性是指[2]若微分方程定解问题 舰(』)=0 x=(。I。2，…，Ⅳ。) in D (1a) PM=g On3D (1b) 存在满足边界条件的解u(z)。现有两个近似解u，(z)，u2(z)和关于解的凸集K=(*；“)， 有“1，“2，u∈K 鄞棵文，男．1qq54年8月生，工学硕士，讲师，(∞1∞2)保定市河北大学数学系。 收稿日期：1999年lO月30日 ——64—— “1(≈)≥u(*)≥U2(z) (2) 成立，当且仅当 Rul(z)≥月u(x)=0≥Ru2(z) (3) 成立时，我们称残差算子R关于解u是单调的。 3 扁壳弯曲问题的单调性 考虑一受均布载荷作用的正方形底面简支球面扁壳，其无量纲形式的平衡方程和边界 条件可表示为 V2V2w+船=一口 (4) 。=0， (V2)。=0 (5) 数，q=q004／肌为无量纲的分布载荷，这里E为壳体的弹性模量，D为壳体的弯曲刚度，t 为壳体的厚度，n为壳体底面的边长，qo为壳体所受的分布载荷集度。 考虑在域中连续且至少具有四阶连续偏导数的试函数”，将其代人控制方程(4)，可得 残值方程 Rw=V2V2圳+kw+口≠0 (6) 对于满足边界条件的非常数解”，若不等式 未(鼢)≥o (7) 在域中恒成立时，则残差Rw关于解”是单调的。由文献[3]知，双Laplace算子可2v2”关于 ”是单凋的，由此不难推断边值问题(1)是单调的。 4 数学规划问题的建立 设扁壳弯曲变形的试函数为 (A，*，y)；∑∑^舢sinm丁cxsinnⅡy (8) m=1．3-=I．3 不难验证所选试函数(8)可以满足边界条件(5)。将其代人控制方程(4)，可以得到内部残值 方程