双非线性四边形平面应力单元分析石拱桥的极限承载力胡柏学蔡.docVIP

双非线性四边形平面应力单元 分析石拱桥的极限承载力 胡柏学1 蔡松柏2 邓继华3 （1．湖南省交通科研所 湖南长沙，410000 2．长沙理工大学 湖南长沙 410000 3．湖南师范大学 资源与环境学院 湖南长沙410081） 摘 要：本文采用共旋坐标法导出了四边形平面应力单元在大转动、小应变条件下的几何非线性和材料非线性的单元切线刚度矩阵；提出了在n+1维荷载-位移空间中求解非线性有限元方程组的增量法迭代格式，该格式保留了迭代刚度矩阵的部分带状性能；在此基础上编制相应的有限元程序NSTSAP。通过对1973年交通部公路桥涵标准图中的23个单跨石拱桥的计算与分析表明本文提出的共旋坐标法四边形平面应力单元的列式、程序和求解器具有较强的非线性分析能力，能较好地分析石拱桥受力全过程的非线性性能。本文计算结果提供了对小跨度石拱桥非线性性能的更多的了解，可供工程设计人员参考。 关键词： 石拱桥; 共旋坐标法; 有限元; 位移增量法; 几何非线性; 材料非线性 Geometrical and Material Nonlinear Quadrilateral Plane Stress Element for Ultimate Bearing Capacity of Stone Bridges 1, 1, Cai Songbai2 and Deng Jihua3 (1.Communication Institute of Hunan Province, Changsha, 410082, P.R. China 2. College of Resources and Environmental Science, Hunan Normal University, Changsha, 410081, P.R. China 3. University of Changsha Communication, Changsha, 410082, P.R. China) Abstract：A tangent stiffness matrix for quadrilateral plane stress element under large rotation with small strain is proposed in this paper. In collaboration with the element, a particular incremental iteration scheme of nonlinear FEM equations in the (n+1)-dimensional load-displacement space is adopted in this nonlinear FEM. In doing so, the method partially retained the band characteristic of the structural stiffness matrix. Based on the above mentioned procedure and its corresponding formulation, a FORTRAN computer program, named NPSSAP, has been developed. Computations for 23 typical stone arch bridges of the Standard Design of Highway Bridges and Culverts of the Communication Ministry have verified that the element’s formulation, the program and the solver in this work have strong nonlinear computation ability, and are cable of analyzing nonlinear behavior of stone arch bridge. Computations presented reflect preferably a comprehensive understanding of nonlinear characteristics of stone arch bridge of small span, could be beneficial for the engineer designers. Keywords：Stone arch bridge, co-rotational procedure, FEM, displacemen