基于ABAQUS复合材料薄壁圆筒的屈曲分析.docxVIP

基于ABAQUS复合材料薄壁圆筒的屈曲分析 由于玻璃钢复合材料的薄壁圆筒结构具有强度高、重量轻、刚度大、耐腐蚀，电绝缘及透微波等优点，目前已广泛应用于航空航天和民用领域中。工程中广泛使用的这些薄壁圆筒，当它们受压缩、剪切、弯曲和扭转等荷载作用时，最常见的失效模式为屈曲。因此，为了保证结构的安全，需要进行屈曲分析。对结构进行屈曲分析，涉及到较复杂的弹(塑)性理论和数学计算，要通过求解高阶偏微分方程组，才能求解失稳临界荷载，而且只有少数简单结构才能求得精确的解析解。因此，只能采用能量法、数值方法和有限元方法等近似的分析方法进行分析。近20年来，随着计算机和有限元方法的迅猛发展，形成了许多的实用分析程序，提高了对复杂结构进行屈曲分析的能力和设计水平。abaqus就是其中的杰出代表。 1.屈曲有限元理论 有限元方法中，对结构的屈曲失速问题的分析方法主要存有两类：一类就是通过特征值分析排序屈曲载荷，另一类就是利用融合newton―raphson运算的弧长法去确认读取方向，跟踪失速路径的几何非线性分析方法，能够有效率分析高度非线性屈曲和后屈曲问题。 1.1线性屈曲 假设结构受的外载荷模式为p0。，幅值大小为λ，结构内力为q，则静力均衡方程应属 λp0=λq 进一步实地考察结构在(λ+△λ)p0载荷促进作用下的均衡方程，获得 ke+kss+λ△s+kgu+λu△u=△λp0 由于结构达至保持稳定的临界载荷时存有△λ，代入上式些 ke+λks△σ+kg△u△u=0 该方程对应的特征值问题为 detke+λks△σ+kg△u=0 如果忽略几何刚度增量的影响，屈曲分析的方程又可以进一步精简为 detke+λks△σ=0 该方程即为为解线性屈曲的特征值方程。λ为屈曲失速载荷因子，△u为结构失速形态的特征向量。 1.2非线性屈曲 非线性屈曲分析方法多使用弧长法展开分步运算排序，在增量非线性有限元分析中，沿着均衡路径运算加速度增量的大小(也叫做弧长)和方向，确认载荷增量的自动读取方案，可以用作高度非线性的屈曲失速问题。与抽取特征值的线性屈曲分析较之，弧长法不仅考量刚度奇特的失速点附近的均衡，而且通过跟踪整个失速过程中实际的载荷、加速度关系，赢得结构失速前后的全部信息，适合于高度非线性的屈曲失速问题。 2.abaqus的线性屈曲分析 abaqus中提供更多两种分析方法去确认结构的临界荷载和结构出现屈曲积极响应的特征形状：线性屈曲分析(特征值屈曲分析)、非线性屈曲分析。 线性屈曲分析用于预测一个理想的弹性结构的理论屈曲强度。它是预期的线性屈曲荷载的上限，可以作为非线性屈曲分析的给定荷载，在渐进加载达到此荷载前，非线性求解必然发散；它还可以作为施加初始缺陷或扰动荷载的依据。所以预先进行特征值屈曲分析有助于非线性屈曲分析，进行特征值屈曲分析是必要的。 3.算例 3.1问题概述 图3-1实例模型 如图所示两端开口的复合材料薄壁圆筒，底端固支，顶端作用有均匀分布的轴压边载。半径r=152mm，高度300mm，厚度t=0.804mm，对称铺层[±45,0]s， 单层厚度为0.134mm。复合材料圆筒的材料参数如下表中： 表1as4/3501-6石墨/环氧的弹性参数 弹性模量e1/gpae2/gpag12/gpag13/gpag23/gpaμ 参数1429.763.63.60.3 3.2建模、分割网格及分析过程 3.2.1abaqus进行前处理建立模型 首先先绘制一个半径为152mm的三维轮廓图，然后展开弯曲300mm建立一个具备平面壳体单元的三维变形体。 其次进入property模块生成材料，参考上表编辑其弹性参数。编辑截面进行复合层板的铺层完成材料属性、截面属性的定义。 图3-2铺层 接着生成装配件，定义屈曲分析步。对所建模型施加边界条件和施加载荷。将圆筒的一端固支，另一端施加单位1的均布轴边压载。 最后展开网格分割。考虑到在复合材料层合壳模型中剪切屋盖的影响，泡果 用厚壳单元s4r来模拟它。所得到的模型如图所示： 图3-3模型网格图 3.2.2计算结果分析 在jod模块中创建屈曲分析模块展开分析，可以获得薄壁圆筒的六阶屈曲失速载荷因子。 表2各阶模态的屈曲载荷因子 模态123456 的关系： pcr=λ?p 屈曲载荷因子 48.015 48.01548.9848.98150.18850.188