桥梁结构有限元分析中的数值计算.pdfVIP

第 25 卷 　第 3 期 重 　庆 　交 　通 　学 　院 　学 　报 2006 年 6 月 Vo125 No3 JOURNAL OF CHONGQINGJ IAOTONG UNIVERSITY June 2006 桥梁结构有限元分析中的数值计算 陈庆军 , 　顾安邦 (重庆交通学院 ,重庆 400074) 摘要 :有限元作为一种数值计算方法 ,其核心在于对总刚度矩阵的各种处理. 本文通过对其中数值方法的讨论 ,指 出以混合编程的方式 ,利用 FORTRAN 上和数值计算资源及其在数值计算上的优势 ,能够很好的解决这个问题 ,提 高程序的速度与效率. 并给出了几个有代表性的函数的例子. 关 　键 　词 :有限元 ;数值计算 ;混合编程 ( ) 中图分类号:TP3111 　　文献标识码 :A 　　文章编号 :1001716X 2006 ( 　　随着桥梁发展的需求 ,能够更真实模拟桥梁结 界面以外 ,其核心程序 求解平衡方程、稀疏矩阵的 ) 构的受力行为 ,更精确计算桥梁各部位内力、应力变 运算等 的计算精度和速度也是至关重要的. 计算精 得越来越重要. 因此随着计算机在桥梁工程中应用 度关系到我们所得到的结果的正确性 ,而计算速度 的普及 ,有限元这种数值方法得到了广泛的应用. 则体现了我们计算成本的大小 ,与桥梁建设的经济 有限元为一种数值方法 ,认为一个结构可以看 性有着直接的关系. 本文就有限元数值计算的精度 作是由有限个力学小单元互相连结而组成的集合 和速度上做了一些探讨. 体 ;表征单元力学特性的刚度矩阵可以比喻作建筑 1 　有限元计算中数值计算方法的选取 物中的砖瓦 ,装配在一起就能提供整个结构的力学 在大型结构分析中 ,总刚度矩阵往往呈现对称 特性. 其主要的步骤有 : 稀疏的特点 ,对其进行约束处理和结点编号后 ,总刚 ) 1 离散结构. 将要分析的结构物分割成有限个 度矩阵成为一个带状稀疏并且正定的矩阵. 对其进 单元体 ,并在单元的指定点设置结点 ,使相邻单元的 行运算耗时较多. 目前对这些数值运算通常有以下 有关参数具有一定的连续性 ,构成一个集合体来代 两种方式 : 替原来的结构. 1) 研究现有的数值方法 , 自己编写程序解决问 ) 2 选择位移模式 [N ]. 用结点的位移去描述结 题. 构中非结点上的位移. 如求解平衡方程组的改进平方根法 ,其算法针 3) 分析单元的力学特性. 导出单元的刚度矩阵. 对对称、正定的方程组 ,较容易编写程序 ;求解矩阵 4) 集合所有单元的平衡方程 ,建立整个结构的 特征值的雅可比方法 ;矩阵的加、减、乘等运算. 这些 平衡方程. (