变厚度矩形弹性薄板的内力分析研究.pdfVIP

第九届全国建筑工程计算机应用学末会议论文集 黑龙x ·哈尔滨 1998 变厚度矩形弹性薄板的内力分析 江耀铃 徐汝新 (深圳市物资集团公司) (深圳市鹏基工业发展总公司) 摘【要】本文用离散型最小二乘法，通过计算机编程对三边固定、一边自由变厚度矩形弹性薄板在静水 压力作用下的情形进行内力分析，解决了此种板一直无可信性成果的问题。所采用的试函数无需满足边界 条件或定解微分方程式，计算结果完全满足工程设计需要而且误差可知。对于屋顶泳池、地下室贮水池等， 池壁按变厚度板设计更能充分利用建筑材料力学性能，节约材料。 一、概 述 在现代多功能建筑中，对于屋顶泳池、地下室贮水池等，传统设计中的结构计算取值 多按 建《筑结构静力计算手册》w七查得，因此设计中基本上按等厚度平板进行计算和配 筋。然而，我们均知道，池壁是受静水压力 (三角形荷载)作用，设计中池面与池底取相 同的壁厚显然未能充分利用建筑材料力学性能。 对于泳池池壁，在传统设计中的结构计算取值其实是按三边固定一边自由的等厚度矩 形弹性薄板受三角形荷载作用下的内力分析来考虑的，而本文探讨的问题则归结为三边固 定一边自由的变厚度矩形弹性薄板受三角形荷载作用下的内力分析问题。由于变厚度板的 问题，到目前为止尚没有能应用于工程设计的成果，而本文作离散型最小二乘法解法，无 需找到满足自由边边界条件，或者满足定解微分方程的函数，同样能得出满意的结果，只 要在边界上或内部配足够的点就行了，大量复杂的计算通过编程由计算机去完成。 二、离散型最小二乘法方程 设有一边值问题待解，它的微分方程及边界条件分别为; Fu一f=0 (于域V内) (1) Gu一g=0 (于域SE.) (2) 式中u为待求的函数，而F及G为算子，f、g为不含u的函数。上两式可以为线性 的或非线性的微分方程，我们若假设一个近似AF数uc(,x)称为试函数((c为待定函数，x 为独立变量)代入上式，由于函数是假设的，一般不能满足方程而出现了误差，称为残数: R,(c,x)=而一f (V域中) (3) RB(c,x)=Gu一g (S边上) (4) 这里R，和R。分别为内部残数及边界残数，于域V中及S上选择有限个点x, 第九届全国建筑工程计算机应用学术会议论文集 黑龙江 ·哈尔滨 1998 (i=1,2 ,m)，将其坐标代入残数方程式中，得到残数方程组: R,(-.) 几(二，)一f(.1 R,(}4) FU-k)一Axr) (5) .RR(cx,,, w}Gu(cx,,,)一g(xA .R,(-.) wG(W(cx.)一g(xm)I 式中w为边界残数的权数相对于内部残数权数的比例，此权数一般定为1 写成矩阵形式为: R=AC一B (6) 其中A为系数矩阵(pxl阶)，C为待求系数列阵(数目为n)，且方程个数p_1. 令，=RTR，根据最小二乘法极值原理通一。，得: cl了 07一a[RTRI一2ATAC一2:ATB=。 GC 厌 了“