-
摘要: 为了研究复杂载荷作用下薄板结构的受力, 基于离散Kirchhoff原理, 推导了三角形弯曲薄板单元的应变矩阵、刚度矩阵中的显式面积坐标积分方法, 根据变分原理构造了任意四边形弯曲单元DKQ, 引入对称及反对称畸变模式进行网格畸变敏感性分析, 将DKQ弯曲单元与平面应力单元组合, 得到了用于薄板分析的四边形平板单元, 并编制了有限元分析程序。计算结果表明: 相对偏移量从-0.18变化至0.18, 反对称畸变模式下挠度最大误差为1.83%, 而对称畸变模式下挠度最大误差为0.99%;对某地铁车体结构, 计算结果与ANSYS结果误差在3.5%之内, 这说明构造的DKQ弯曲单元对网格畸变不敏感, 具备良好的位移解收敛性和计算精度。Abstract: In order to analyze the states of thin-walled structures under complicated loads by discrete Kirchhoff principle, the strain matrix of triangle bending element and explicit area coordinate integral method were discussed in detail.According to variant principle, DKQ (Discrete Kirchhoff Quadrilateral) bending element was constructed.By introducing symmetric and antisymmetric mesh distortion modes, mesh distortion sensitivity analysis was carried out.Through combining DKQ element with planar stress element, quadrilateral plate element used for general thin-walled structures analysis was implemented by finite element program.Numerical test result indicates that when relative offset value varies from-0.18 to 0.18, the deflection maximum errors of antisymmetric and symmetric modes are 1.83% and 0.99% respectively.For a car-body of some subway vehicle, the differences between the analysis values and those obtained by ANSYS are less than 3.5%.The result shows that DKQ bending element is insensitive to mesh distortion and has good convergence capability and accuracy for displacement field.
-
Key words:
- vehicle engineering /
- structure analysis /
- finite element /
- qudrilateral bending element
-
表 1 C的数值解
Table 1. Numeric Solution of C
理论解[9] 简支: 0.011 60 固支: 0.005 61 N 数值解 数值解 2 0.012 063 0.005 979 3 0.011 833 0.005 822 4 0.011 744 0.005 745 5 0.011 698 0.005 704 
下载: 导出CSV
表 2 结果比较
Table 2. Comparison of Analysis Results
当量应力/MPa X/mm Y/mm Z/mm ANSYS 135.595 -0.679 21 3.920 9 -4.021 8 本文方法 131.421 -0.655 59 3.863 8 -4.008 3
下载: 导出CSV
-
[1] 姚曙光, 许平. 高速磁浮列车国产化[J]. 交通运输工程学报, 2004, 4 (2): 40-44. http://transport.chd.edu.cn/article/id/200402010YAO Shu-guang, XU Ping. Chinese Domestic Components of High-Speed Maglev Train Body[J]. Journal of Traffic and Transportation Engineering, 2004, 4 (2): 40-44. (in Chinese) http://transport.chd.edu.cn/article/id/200402010 [2] 郜永杰, 翟婉明. 轨道结构强度有限元分析[J]. 交通运输工程学报, 2004, 4 (2): 36-39. http://transport.chd.edu.cn/article/id/200402009GAO Yong-jie, ZHAI Wan-ming. Finite Element Analysis of Track Structure Strength[J]. Journal of Traffic and Transportation Engineering, 2004, 4 (2): 36-39. (in Chinese) http://transport.chd.edu.cn/article/id/200402009 [3] BATOZ J L, BATHE K J, HO L W. AStudy of Three-Node Triangular Plate Bending Elements[J]. International Journal for Numerical Methodsin Engineering, 1980, 15 (12): 1 771-1 812. doi: 10.1002/nme.1620151205 [4] BATOZJ L, TAHAR MB. Evaluation of a New Quadrilateral Thin Plate Bending Element[J]. International Journal for Numerical Methods in Engineering, 1982, 18 (11): 1 655-1 677. doi: 10.1002/nme.1620181106 [5] BERNADOU M. On the Convergence of a Discrete Kirchhoff Triangle Method Valid for Shells of Arbitrary Shape[J]. Computer Methods in Applied Mechanics and Engineering, 1994, 118 (4): 373-391. [6] PIMPINELLI G. An Assumed Strain Quadrilateral Element with Drilling Degrees of Freedom[J]. Finite Elements in Analysis and Design, 2004, 41 (3): 267-283. [7] CHEN Wan-ji. Refined Non-Conforming Triangular Elementsfor Analysis of Shell Structures[J]. International Journal for Numerical Methods in Engineering, 1999, 46 (3): 433-455. [8] 彭国伦. Fortran 95程序设计[M]. 北京: 中国电力出版社, 2002. [9] 龙志飞, 岑松. 有限元法新论原理. 程序. 进展[M]. 北京: 中国水利水电出版社, 2001. [10] BATHE KJ. Finite Element Procedures[M]. NJ: Prentice Hall, 1996. -
点击查看大图
图(6) / 表(2)
计量
- 文章访问数: 301
- HTML全文浏览量: 115
- PDF下载量: 884
- 被引次数: 0
