中厚板的弯曲问题理论研究进展

1、中厚板的弯曲问题理论研究进展古典的薄板理论，对世纪板进行分析的时候引入了克希霍夫-勒夫假设，假定板的横向变形为零，形式上相当于假定垂直于板中面的各个面内剪切模量无穷大，而板的其他各个方向的弹性常数仍是真是材料的弹性常数。或者说，在古典薄板理论中，用这种各向异性材料代替了真实的各向同性材料。在大量由板壳组成的工程结构中，根据Kirchhoff假设建立的薄板近似理论的计算结果已经满足工程计算的精度，但如果板比较厚，或者即使板很薄，在集中力作用点附近以及薄板边界周围，近似理论不仅不能取得满意的结果，甚至会导致错误的结论。近似理论的这些缺点要求研究者提出一些新的假定，建立新的理论，它既能避免数学方面的

2、困难，又能克服采用Kirchhoff假定忽略横向剪切应变所引起的误差。随着工程技术要求的不断提高，薄板与薄壳理论已无法满足要求，促使国内外力学工作者对厚板壳理论与中厚板壳理论进行了大量的研究。 为解决上述问题，在二十世纪中期，以赖斯纳（Reissner）为代表的一批学者，提出了考虑剪切变形的板的理论，一般称为中等厚度板（中厚板）理论。Reissner关于中厚板的基本假设有：1. 与薄板理论相同，应力沿板厚仍是线性分布的。2. 横剪力引起的变形不能略去。即薄板理论的基本假定不成立。 Reissner 首先采用直线假设， 即变形前垂直于中面的直中法线变形后仍为直线，代替直法线假设，考虑了横向切应力

3、对板变形的影响，该理论同时计及法向应力对应变的影响，采用广义余能的变分原理导出考虑剪切变形的基本方程以及边界条件，最终得到如下一组的基本方程：消去，得到六阶偏微分方程： （1）式中w为壳中面挠度，q为横向荷载，为中面法线转角。(1)式的定解问题在每一条边上需要3个自然边界条件。Reissner用此理论研究了矩形板的扭转和有圆孔的无限大板的弯曲、扭转。之后Green A，Mindlin R D相继发展了Reissner 理论。几乎在Reissner 理论的同时，Hencky H也提出横向剪切的精化理论，假设位移u ，v 是z的线性函数，横向剪切力与z无关。因此实际上是考虑横向切应力平均值对板变形

4、的影响。该假设忽略法向变形z的影响，利用Lagrange 的虚位移原理，导出与( 1) 中第1式相似的挠度方程式及包括转角的六阶方程组。50 年代初期， Kromm A 提出另一种精化理论。他力求严格满足三维弹性体应力平衡方程，以克服Reissner 和Hencky论中的不足之处。其主要特点是假定是z 的函数，由应力、应变关系、应力平衡微分方程和上下两面边界条件， 最终得到基本方程组为：式中函数f(z)为挠曲函数， 即根据变形状态是弯曲还是剪切进行选择。苏联学者Bacob B 在1957年提出考虑切应力对变形影响的平板理论，假定变形前垂直于中面的直法线变形后发生扭曲，且剪切沿厚度方向按抛物线规

5、律分布；平板处于平面应力状态等。利用Lagangian变分方程导出的方程组与Reissner 理论相比，可以发现后者方程中被略去的项和保留的项都属阶的量。1975年Vladimir Panc在他的专著中介绍了各种主要的精化理论，包括由他提出的逐次近似法。该法的基本思想是经典理论解属第1次近似解，第2次近似是利用经典理论的解计算横向切应力x z，yz 对板变形的影响以修正直法线假设引起的误差，第3次近似是考虑与横向切应力相应的齐次应力状态。所导出的方程组与Reissner理论相比，扭矩表达式完全相同，而其余内力因素和平衡基本微分方程式只是在一些次要项上有差别，这是因为逐次近似理论中没有包括法向应

6、变z 的影响。Levinson M提出一种精化理论，与Reissner- Mindlin理论相比新理论允许板变形时法线可以扭曲，因此有可能满足上下表面切应力自由的条件，依此就能导出以挠度w，中面处横截面转角x，y 为变量的3个解偶微分方程，由这组方程得到的各种载荷下周边简支板的结果比Mindlin 理论更加精确。其实该理论的基本假设、位移表达式和基本方程式与上述Bacob理论完全相同，唯一的差别只是中面内线素的旋转角的正向定义与后者正好相反。上述这些精化理论，特别是Reissner理论，目前广泛应用于复合材料板壳理论、板壳断裂和板壳有限元分析中。随着计算机的发展，数值计算越来越得到人们的重视，

7、在此基础上发展了关于厚板理论的一系列数值计算方法，比如有限元法，加权残差法，差分法，变分法，边界元法，有限条法等许多近似方法。本文中将主要简单地介绍下具有代表性的运用比较广泛的加权残差法以及有限元法。1. 加权残数法在工程技术和科学研究中，我们经常会遇到各种各样的定解问题。即在一定的几何边界条件或初始条件下，求解问题的控制方程，而目前关于控制方程中的微分方程或积分方程的研究，仅仅能得到少量问题的精确解。因此，大量的定解问题只能通过近似求解来实现。加权残数法就是诸多近似方法中的一种。加权残数法用于解板的弯曲问题，需要在全域内积分或配点，配点数量及位置都直接影响计算精度。1.1 中厚板弯曲的基本方

8、程现在考虑一中厚板，厚度为h，承受横向分布荷载q(x，y)作用。对于各向同性中厚度矩形板的弯曲，为了简化方程，引入函数F(x，y)和f(x，y)，得到转角和挠度方程为：而函数F(x,y)和f(x,y)满足方程：式中D为板的抗弯刚度，G 为抗剪模量，为泊松比。板的内力方程为：考虑到板的静力平衡, 则得中厚板的基本方程：其中k为截面系数。根据分解刚度法, 得到板的总势能泛函为：其中，式中Wb、Ws 分别为板的弯曲和剪切挠度，b、s 分别为板的弯曲和剪切势能。1.2 加权残数解法弹性力学问题常可归结为求一未知的(列阵)函数，使其在给定域内满足该问题的偏(或常)微分方程及边界条件，即：对于多变量函数的

9、泛函变分问题, 康托洛维奇提出设泛函的自变函数为满足边界条件的函数系列( ，)及待定函数，将变分问题的近似解写成：其一级近似式为将式中右端两函数之一选定为满足边界条件的某种函数，由泛函变分极值条件得到关于另一函数的常微分方程，进而求出此方程的近似解，便得到问题的一级近似解。设有两个自变量和两个自由函数的泛函具有：取u,v的一级近似表达式为：式中， 为满足y 方向边界条件的试函数。将式(13)代入式(12)并沿y积分有：使泛函极值=0实现的两个函数和A应当满足欧拉方程和相应的边界条件：显然，上式为常微分方程。一般来说，求解式(15)的过程较复杂，特别是对于二阶以上近似解就更为复杂，为此，本文选取

10、满足边界条件的试函数：将式(16)代入式(15)，得到包含有待定关系ai 和bi的残值Ra和Rb，取权函数Wai和Wbi使得残值在某种平均意义上等于零。据此，可得到用于求解待定系数ai、bi (i = 1, 2, r ) 的代数方程组。最后得到问题的近似解。现以x方向为例，给出适合于不同边界条件的试函数。(1) 简支简支试函数:(2) 固端固端试函数:(3) 固端简支试函数:(4) 固端自由试函数:(5) 简支自由试函数：1. 有限元法在板壳有限元分析中，位移型非协调元薄板单元得到成功的应用，解决了大量的工程问题。随着各种新型材料的出现和广泛应用，构造计入横向剪切变形的板单元是非常必要的。在位移型单元研究方面，已经提出了一些厚薄板通用的单元。以Reissner- Mindlin理论为基础，构造了一种计入横向剪切变形的矩形板单元；在位移模式中，挠度和转角采用不同阶次的插值函数。数值算例表明，该单元对工程中的薄板、中厚板均能得到令人满意的数值结果。1. 有限元公式基本符号与关系式广义位移向量 广义应变向 曲率和扭率 广义应力向量 几何关系 式中弹性关系