加筋板的局部总体屈曲行为-中国力学学会.docVIP

加筋板结构局部总体屈曲分析 张连文1 曾广武2 计宏伟1 （1．天津商业大学 机械学院，天津300134） （2．华中科技大学 交通学院，武汉430074） 摘要：提出一套能够较为准确反映板筋耦合作用且可同时计算加筋板局部和总体屈曲行为的理论方法。运用卡门大挠度理论和附连翼板有效宽度的概念，采用铁木辛柯能量法做近似分析。然后，针对典型边界状况，给出了用双三角级数表示位形函数且具有初始缺陷加筋板的总势能及其变分式。最后编制了程序计算系统BBASP。算例不仅表明该方法是可行和有效的，而且还揭示出，对于特定的板筋刚度比，在一定边长比的范围内存在二次总体屈曲现象。 关键词：加筋板，稳定性，屈曲载荷，二次屈曲，大挠度理论，能量法，变分，有效宽度 引 言 加筋板结构广泛应用于舰船、建筑、航天飞行器等工程领域里。局部屈曲和总体屈曲在研究加筋板的稳定性时需要同时考虑。若筋条满足必要刚度，则加筋板在端面受压时局部和总体屈曲会相继发生。 如何提高这类结构的抗屈曲破坏能力，在过去的40年中一直是工程界和力学界所关注和研究的热点。期间许多学者曾做过大量的试验研究和理论分析工作，但进展仍然不快。因为板筋刚度差异很大导致屈曲破坏形式非常复杂，还有板列屈曲及后屈曲特性、板筋屈曲耦合效应、加筋板长宽比、板筋刚度比、边界条件以及初始缺陷等诸多因素都对破坏产生复杂影响，因此增加了问题的分析难度。 目前，有两种处理这类问题的力学模型理论，即离散梁理论和正交异性板理论[1]。当筋条较为稀疏时，一般采用离散梁理论，不过该方法只能计算总体屈曲载荷。当筋条配置较密时，可以按正交异性板处理[2]，但这种方法比较复杂，且当横向加强筋不超过三根时正交异性板理论会招致显著的精度损失[1]。此外，不论离散梁理论还是正交异性板理论，他们都是单方面考虑了板对筋或筋对板的作用。从而较难反映板筋在屈曲过程中的相互作用。 从研究屈曲行为的数学手段上讲，有半解析法（能量法）和数值法。能量法以位形函数直观，且计算量相对较小而被工程界所普遍采用。常见的能量近似法有铁木辛柯（S. P. Timoshenko）法、里茨（W. Ritz）法以及沃耳米勒（А. С. Водьмцр）法 [3, 4, 5]。其中铁木辛柯法是先假定满足边界条件的位形函数（），并通过这些函数计算系统的总势能，然后应用最小势能原理求出的待定参数。里茨法只假定挠屈函数，通过求解应变协调方程得到应力函数，并以和写出系统总势能表达式，再由最小势能原理求出的待定参数。沃耳米勒法与里茨法十分相似，它也是只假定挠度的表达式并利用应变协调方程得出应力函数，但它是直接利用伽辽金（B. Garlerkin）法求出的待定参数。由于应变协调方程的高度非线性性，只对较为特殊的边界问题才能直接求出应力函数，所以铁木辛柯法是求解一般边值问题普遍而有效的方法。关于数值方法，有差分法和有限元法等。一般来讲采用数值方法得到的结果更为准确，但计算量大并且不能给出直观的位形函数。这方面的工作目前也不少，但本文不作探讨。 前人在利用铁木辛柯法研究加筋板结构的屈曲破坏特性时主要存在如下不足。其一是只用小挠度理论[2]。其二是虽用大挠度理论，甚至也考虑了筋条的偏心距，但在研究总体屈曲时，忽略了板筋刚度在屈曲过程中的变化[6]。其三是尽管用大挠度理论考虑板的应变能，但在考虑外力势能时仍用小挠度理论，甚至常被许多学者所忽略[2, 3, 7]。 我们看到为研究加筋板稳定性而提出的各种方法均不另人满意。 本文运用冯·卡门大挠度理论和附连翼板有效宽度的概念，采用铁木辛柯能量法做近似分析，提出一套能够较为准确反映板筋耦合作用且可同时考虑加筋板局部和总体屈曲行为的理论计算方法，将对加筋板结构的稳定性问题进行较为深入系统的研究，并试图揭示一些新的屈曲特性。为此，首先将在第二节给出加筋板结构的总势能表达式。第三节将给出满足典型边界状况，用双三角级数表示的位形函数。最后用编制的程序计算系统BBASP，通过对加筋板长宽比、板筋刚度比的变化试算了一系列例子，同时与以往的方法所得的结果进行比较并给出合理的解释。算例还特别揭示出，对于特定的板筋刚度比，在一定边长比的范围内存在二次总体屈曲现象。 1 总 势 能 及 其 变 分 设有一根中心纵向加筋的板格，取面板的中面为坐标面，轴垂直于面，并按右手系规则，如图1所示。 图1 加筋矩形薄板及其坐标示意图 设板的边长为, 厚度为, 加强筋的腹板高度为, 厚度为。则矩形板的边长比为，板筋刚度比为。其中，为板、筋杨氏弹性模量，为板的抗弯刚度，且，为泊松比。设初始缺陷，在单位长度的端面外载荷的作用下，板中面位形函数为，且压缩荷载。 加筋板的总势能由平板势能