局部载荷作用下四边固支矩形板弹性解答及其应用

2019-08-16陆锁芳陈伟杰

舰船科学技术 2019年7期

陆锁芳，刘俊，陈伟杰

(1. 上海交通大学海洋工程国家重点试验室，上海 200240；

2. 高新船舶与深海开发装备协同创新中心，上海 200240)

0 引言

对于四边固支矩形板的弹性解答，很多学者已进行了较为深入的研究，主要针对承受均布载荷的矩形板[1]，钟阳[2]提出利用辛几何法推导得到均布载荷下四边固支矩形弹性薄板弯曲问题的解答；李元媛[3]提出采用多项式函数结合Navier 方法来得到近似的矩形薄板弯曲挠度，IMRAK[4]结合三角函数和双曲函数，在忽略高阶小量的前提下求解了四边固支矩形板受均布载荷的最大位移。均布载荷下四周刚性固定矩形板解答可用于对船体中板格进行分析，如承受分布货物荷重的甲板板。但一些特殊情况下，比如海洋平台或军用船舶上直升机的起降，滚装船上轮印载荷的作用[5]，以及极地船舶舷侧加强部分承受冰载荷的作用[6]等，此时分布载荷只作用于板格局部位置，不能再将板格视为承受均布载荷四周刚性固定矩形板求解，将其简化为局部分布载荷下矩形板的问题近似考虑更为合理。对于承受局部分布载荷的刚性固定矩形板，张文福[7]曾尝试给出了板中弯矩和挠度的计算用表，但仅适用于泊松比为0.2，矩形板的长宽比a/b 为1，1.2，1.4，1.6，1.8 和2，局部载荷的边长与对应矩形板边长比c/a 为0.2，0.4，0.6，0.8 和1 的部分情况，应用具有局限性。

本文基于能量法，推导了局部载荷作用下四边固支矩形板弹性阶段的解答方程组，利用Matlab 编程可方便得到挠度计算值。对于冰级船舶，由于与冰的遭遇工况多，如果采用非线性直接计算，在工况繁多的情况下计算工作量会很大。结合极限载荷准则，将本文提出的解答应用于极地船舶极限冰载荷的快速校核中[8]，可提高初步设计阶段的效率。

1 局部载荷作用下四边固支矩形板的理论解答

图 1 局部载荷作用示意图Fig. 1 Diagram of local distributed load

满足位移边界条件的挠度 w表达式为：

其中： M ， N 分别为 x ， y 方向展开的项数；为基函数的系数。

板的应变能为：

力函数为：

将挠度表达式代入总位能公式中得到：

由此得到关于系数amn的线性方程组。

2 局部载荷作用下四边固支矩形板的近似解

随着项数 M ，N 的增加，挠度计算精度也不断增加直至收敛，同时随着项数的增加，求解式（5）的计算量也会陡增，故此借助Matlab 程序进行求解。

将式（5）改写成如式（6）的形式求解，

并有

求出系数 A后，代入式（1）即可求得矩形板的挠度 w，该方法可以求任意载荷面积下的矩形板挠度。

2.1 均布载荷作用下解答验证

为验证本文方法的正确性，与四边刚性固定矩形板受均布载荷作用的经典解答进行对比。弹性解答矩形板中点的挠度为取值见表1[9]。

表 1 四边刚性固定的矩形板在均布载荷作用下的挠度解答Tab. 1 Classic solution of rectangular plate by local distributed load with four edges clamped

将本文的计算方法退化到均布载荷情况，即取s=a,h=b，选取b=100 mm t=4 mm，E=2.06×105MPa 进行验证，2 种方法计算结果如表2 所示。

表 2 本文方法与经典解答对比Tab. 2 Results comparison between two methods

可以发现2 种方法得到的结果非常一致，本文解答在均布载荷的特殊情况下是正确的。

2.2 局部载荷作用下本文解答与数值解的一致性验证

表 3 本文解答与数值解对比Tab. 3 Comparison with numerical solution

可知，本文解答与数值解差别很小。

3 局部载荷作用下四边固支矩形板弹性解答的应用

3.1 基于2 倍弹性斜率法的快速校核方法

近年来，随着极地环境不断变化，越来越多的目光转向极地船舶的发展。在求解极地船舶所受极限冰载荷作用时，可使用有限元软件进行非线性有限元计算，但冰船遭遇情况复杂，工况繁多，且对网格精度要求高，需要极大的工作量。陈伟杰等[10]提出结合载荷-变形曲线和2 倍弹性斜率准则得到极地船舶强度评估方法，由此快速得到极限冰载荷的初步设计值，现将本文的矩形板弹性解答应用于该校核方法中。

修正后的矩形板载荷-变形公式的塑性阶段解析解[10]