小应变几何非线性问题有限元求解.docVIP

小应变几何非线性问题有限元求解 　　摘要： 本文以小应变几何非线性问题为研究对象，以最小势能原理或虚位移原理基础，对小变形几何非线性问题的求解方程作有限元分析；按修正拉氏表述、全量拉氏表述给出小变形几何非线性有限元方法求解位移基本量的迭代公式，这种迭代公式对于小变形几何非线性问题的有限元分析具有非常重要的理论价值及应用价值。 　　Abstract: Taking the geometric nonlinear problems of small strain as research object, the paper makes finite element analysis for the solving equations of small deformation geometrically nonlinear question by the principle of minimum potential energy or imaginary displacement principle basis. According to the fixed Laplace expression, total content Laplace expression, the small deformation geometry nonlinear finite element method to solve the basic amount of iterative formula displacement is given, and this iterative formula has a very important theory value and application value for small deformation geometrically nonlinear finite element analysis. 　　关键词： 有限元方法；几何非线性；小应变；拉氏表述 　　Key words: finite element method；geometry nonlinear；small strain；Laplace expression 　　中图分类号：O34 文献标识码：A 文章编号：1006-4311（2011）２4-0044-02 　　0引言 　　有限元方法[1-3]在实际工程中得到了广泛的应用，针对小应变几何非线性问题有限元求解方面有较多的应用成果[5-8]，如大跨度桥梁几何非线性问题[5]、大跨度斜拉桥非线性静力分析[6]、平面梁杆结构几何非线性分析[7]及平面桁架的几何非线性有限元分析[8]等，同时针对基础理论的公式推导也有相应的成果，如微分刚度计算[9]、两节点曲线索单元精细分析的非线性有限元法[10]等。这些成果都有较好的工程应用前景，笔者采用理论推导方式，在有限元基本求解方程的基础上，结合预应力锚固体系的理论分析及实践研究[4]的课题项目，推导可直接用于课题研究的有限元求解迭代公式。 　　1问题提出 　　结构出现大位移时，就属于几何非线性问题，平衡方程应在变形后的位形上建立，而变形后的位形是比较复杂的，通常用变形前已知坐标来表示，只要位移不符合“小位移”或“小变形”的条件，在应变位移关系中就有非线性项，结构刚度阵就不再是常数，而是位移的函数，致使有限元方程组成为非线性的代数方程组。在大位移问题中，结构各部分应变可能是大的，也可能是小的，但在梁、板、壳等结构中，常常是大位移小应变。如图1所示，悬臂梁中微小段AB位移到A′B′，从局部看A′B′的应变是微小的，从整体看位移是很大的。 　　在工程结构中小应变大位移问题普遍存在，有单独研究必要，同时在小应变问题中，有关应力和应变的定义可沿用线性理论中人们熟悉的概念，比较容易理解，至于大应变问题，需重新定义应力和应变，一般用连续介质力学理论，连续介质力学描述位移有两种方法，以变形前初始坐标为自变量的拉格朗日坐标或物质坐标；以变形后瞬时坐标为自变量的欧拉坐标或空间坐标。 　　２修正拉氏表述 　　大位移的主要影响在于不能忽略位移引起的整体几何形状的变化，如图2所示，单元由AB位移到A′B′，其方向的变化（转动角θ）就是不能忽略的几何形状的变化，它是节点位移uA、vA、uB、vB的函数。 　　 　　迭代过程收敛后，外荷载可增加到另一新水平，再开始新的迭代，找到新的平衡位置，这样就把迭代法与增量法结合起来完成最后的求解。 　　3全量拉氏表述 　　拉氏描述法始终以初始位置为参考系，一切位移、刚度、微分、积分等都相对初始位置为参考，在变形过程中，相对初始位置的位移是位移全量，如图3所示，一端固定铰接的杆单元，长为L截面为A，节点位移仅有u1、v