用位移变分法求问题.docVIP

用位移变分法求解问题 一 求解问题： 铅直平面内的正方形薄板，边长为2a，四边固定，图1，只受重力作用。设μ=0，试取位移分量的表达式为 图1 图1 O y a a a a x x 用里兹法或伽辽金法求解（在u的表达式中， 布置了因子x和y，因为按照问题的对称条件， u应为x和y的奇函数）。 二 求解步骤： 取坐标如图1所示，位移分量为 （1） （1） 上述位移分量满足位移边界条件 现在，式（1）取各式前两个系数为待定系数，也就是取 （2） （2） 将式（2）代入式（3），即 （3） （3） 其中，μ=0，得 （4） （4） 因为在这里可得 （5） （5） 由式（5）得 （6） （6） 解式（6）得， （7） （7） 将式（7）代入式（2）得位移函数 （8） （8） 将位移函数式（8）代入位移-应力函数中，可得 （9） （9） 三 MATLAB求解 应力空间图 这里取， 由应力函数可到以下图2，图3，图4。分别代表x方向正应力，y方向正应力，剪应力 图2 方向正应力 图2 方向正应力 图3 方向正应力 图3 方向正应力 图4 剪应力 图4 剪应力 由应力图很容易观察到图像关于x轴y轴的对称性，对比应力的函数的奇偶对称性正好相对应。图2知