考虑后屈曲行为的单层球面网壳稳定承载力分析-analysis of stability bearing capacity of single-layer spherical reticulated shell considering post-buckling behavior.docxVIP

摘要屈曲是造成连续薄壳与单层网壳结构失效的重要原因之一。研究发现经典 屈曲理论对薄壳屈曲荷载的估计值偏高，因此许多学者尝试对屈曲荷载进行更 准确可靠地预测。本文考虑单层球面网壳的后屈曲行为进行研究，以期对其屈 曲荷载进行准确估计。薄壳结构有限元分析本文使用 ABAQUS 建立有限元模型，考虑几何缺陷，对壳体的后屈曲行为 进行非线性分析，求解采用了 Riks 弧长法跟踪荷载位移曲线。分析结果表明， 当初始几何缺陷较小时，初始屈曲荷载与经典屈曲理论计算值非常接近，然而 随着初始缺陷的增加，模拟的初始屈曲荷载会显著下降，与经典理论值差异很 大。对于球面薄壳，不同几何缺陷的荷载-位移曲线在后期趋向于重合，形成一 “平台”，表明初始屈曲荷载对几何缺陷非常敏感，而后屈曲荷载对几何缺陷不 敏感。这是经典屈曲理论高估屈曲荷载的一个主要原因。 薄壳结构屈曲荷载的经验预测实际工程中，几何缺陷主要来自建造误差。经典屈曲理论并未考虑结构的 非线性这一重要因素，而这一因素一般被认为是经典理论高估屈曲荷载的另一 主要原因。本文使用有限元法模拟时考虑了几何非线性，计算出了后屈曲荷载。 以往对轴心受压圆柱壳的研究表明，后屈曲荷载值为柱壳屈曲试验所得承载能 力的平均值，薄壳的屈曲应力与厚度的 1.5 次方成正比，由此预测屈曲荷载与 厚度的 2.5 次方(t2.5)成正比，该预测得到了有限元分析验证。因此采用该荷载 值作为后屈曲荷载可避开几何缺陷的影响。采用“平台”处的荷载作为屈曲荷 载更具现实意义。有限元分析表明，薄壳后屈曲后变形的一种形式是“局部凹 陷”，而应变能集中在“局部凹陷”的周边。因此本文利用功能原理对“局部凹 陷”进行分析，推导出了球面薄壳的后屈曲承载力经验公式，并依据拟壳化处 理，将这一承载力计算式扩展到单层球面网壳中。利用拟壳法预测单层球面网壳屈曲荷载利用拟壳法可将单层球面网壳可以看成连续壳进行分析。对单层球面网壳 有限元分析发现其后屈曲行为与球面薄壳类似。初始屈曲荷载对初始几何缺陷 敏感，而后屈曲荷载对初始几何缺陷并不敏感。因此，根据后屈曲荷载，将单 层球面网壳等效为球面薄壳，并确定了球面薄壳的等效厚度teq 。依据薄膜刚度 以及弯曲刚度推导了等效薄膜厚度tm以及等效弯曲厚度tb。不改变材料弹性模 量，通过teq建立了三向网格单层球面网壳与球面的等效关系。 单层球面网壳屈曲试验为了验证经验公式以及有限元模拟结果，本文进行了三向网格单层球面网 壳屈曲试验。试验采用模拟均布荷载加载，记录了网壳顶点的变形，形成了荷 载-位移曲线。除其中一个模型形成局部凹陷后承载力可以继续上升外，其余模 型均在屈曲后垮塌。试验结果、有限元模拟结果以及经验公式计算结果吻合良 好，从而验证了本文的可靠性和有效性。关键词：等效厚度，有限元，非线性，后屈曲，球面网壳，球壳AbstractBuckling of thin continuous shells and single-layer reticulated shells has attracted interest of many researchers because it is the main cause of failure for such structures. Many attempts have been made to predict the buckling load as it has been found out that the classical buckling theory overestimates this load. The purpose of this study is to predict a more realistic buckling load for thin spherical shells and single-layer reticulated spherical shells considering their post-buckling behavior.Finite Element (FE) modeling of thin shell structuresThe post-buckling behavior of thin spherical shells and single-layer reticulated shells was studied by conducting nonlinear FE modeling using ABAQUS, where initial geometrical imperfections were introduced and Riks method used to trace the load deflectio