Finite Element Solution for Torsion Stress Function With Arbitrary Multi-Connected Section
-
摘要: 杆件扭转问题的求解,主要有基于扭转理论翘曲函数的边界元法和有限元法、基于薄壁杆件理论的数值解法和基于扭转理论应力函数的有限元法.根据任意多连通截面直杆扭转问题的应力函数理论,讨论并改进了与微分方程及定解条件等效的泛函,在此基础上推导了求解多连通截面扭转应力函数的有限元列式,将扭转问题的翘曲位移单值条件转化为边界节点上的集中荷载.采用主从节点法满足孔洞边界上应力函数的同值条件,实现了任意多连通复杂截面扭转应力函数的有限元直接求解,通过应力函数积分获得截面的扭转常数.算例验证了方法的可行性和有效性.Abstract: The major three methods can be used to solute the torsion bars’torsion problem. One is the boundary element method and the finite element method that is based on the warping function of torsion theory, the others are numerical solution based on the thinwall theory and the finite element method based on the torsion stress function of torsion theory. According to stress function theory of torsion bars with arbitrary cross section, a functional equivalent to the torsion’s differential equation and definite condition was discussed and improved, finite element formulas were deduced to solute the torsion stress function for multi-connected section, the boundary condition of single warping-displacement value was changed to concentrated force loaded on boundary nodes. The condition that the stress function must be constant value on each hole boundary was satisfied by using master-slave node method, so the torsion stress function with arbitrary multi-connected complex section could be obtained directly by finite element method, and the torsion constant was solved by integrating from the torsion stress function. Examples verified the feasibility and validity of this method.
-
Key words:
- multi-connected section /
- torsion /
- stress function /
- finite element
-
[1] Wilson H B, Farrior D S. Computation of geometrical and inertial properties for general areas and volumes of revolution[J]. Computer Aided Design, 1976,8(4):257-263. [2] Miles R G, Touh J G. A method for the computation of inertial properties for general areas[J]. Computer Aided Design, 1983,15(4):196-200. [3] Eberh A R, Allen C, Williard G H. Calculating precise crosssectional properties for complex geometries[J].Computers in Mechanical Engineering, 1987,6(2):32-37. [4] Friedman Z, Kosmatka J B. Torsion and flexure of a prismatic isotropic beam using the boundary element method[J]. Computers and Structures, 2000,74(4):479-494. [5] 周凌远, 李乔, 张士中. 采用有限单元法计算梁任意形状截面特性[J]. 计算力学学报, 2008,25(5):634639.(ZHOU Ling-yuan, LI Qiao, QIANG Shi-zhong.Evaluation of arbitrary cross sectional properties with finite element method[J].Chinese Journal of Computational Mechanics, 2008,25(5):634639.(in Chinese)) [6] 刘悦藏, 范慕辉. 凸多边形截面扭转杆问题的数值解法[J]. 河北工业大学学报, 1999,28(6):7375.(LIU Yue-cang, FAN Mu-hui.The numerical method for solving the torsional shafts with convex polygon cross sections[J].Journal of Hebei University of Technology, 1999,28(6):73-75. (in Chinese)) [7] 周宗宪, 李宁军.任意多边形截面直杆扭转的加权残值解[J]. 西安公路交通大学学报, 2000,20(2):115-118.(ZHOU Zong-xian, LI Ning-jun.Solution of torsion problems of the prismatic bar with arbitrary section applying weighted residual method[J].Journal of Xi’an Highway University, 2000,20(2):115-118. (in Chinese)) [8] 崔振山, 刘国燕, 刘才.多连通任意截面杆件扭转问题的有限元法[J]. 燕山大学学报, 2001,25(4):297-300.(CUI Zhen-shan, LIU Guo-yan, LIU Cai.Finite element method for torsional bar with arbitrary shaped and multiple-connected cross-sections[J].Journal of Yanshan University, 2001,25(4):297300.(in Chinese)) [9] 李永奇, 张卫红.模拟求解任意截面杆扭转的化复为单有限元法[J]. 兰州理工大学学报, 2006,32(2):158161.(LI Yong-qi, ZHANG Wei-hong.A special simulation for torsional bars with arbitrary cross-section using finite element method[J]. Journal of Lanzhou University Technology, 2006,32(2):158-161.(in Chinese)) [10] 何琳, 王家林, 张德立.任意多连通截面扭转应力函数的温度场比拟解法[J]. 重庆交通大学学报(自然科学版), 2012,31(5):927-929(934).(HE Lin, WANG Jia-lin, ZHANG De-li.The temperature field analogical solution for torsion stress function with arbitrary multi-connected section[J].Journal of Chongqing Jiaotong University(Natural Science), 2012,31(5):927-929(934).(in Chinese)) [11] Gurujee C S, Shah K R. A computer program for thin walled frame analysis[J]. Adv Eng Softerwear, 1989,11(2):58-70. [12] Yoo Chai Hong, Acra V S. Cross-sectional properties of thinwalled multi-cellular section[J].Computers and Structures, 1986,22(1):53-61. [13] Aleksandar P. Computer program for determination of geometrical properties of thinwalled beams with openedclosed section[J]. Computers and Structures, 2000,74(6):705-715. [14] Lubarda V A. On the torsion constant of multicell profiles and its maximization with respect to spar position[J]. Thin-Walled Structures, 2009,47(6/7):798-806. [15] 杜柏松, 葛耀君, 周峥.复杂薄壁断面的自由扭转惯性矩的数值分析[J]. 同济大学学报(自然科学版), 2006,34(9):1170-1174.(DU Bai-song, GE Yao-jun, ZHOU Zheng. Numerical analysis of free torsion constant for complicated thin-walled cross section[J].Journal of Tongji University(Natural Science), 2006,34(9):1170-1174.(in Chinese)) [16] 钱伟长, 林鸿荪, 胡海昌, 叶开沅.弹性柱体的扭转理论[M].北京:科学出版社,1956: 187-206.(CHIEN Wei-zang, LIN Hong-sun, HU Hai-chang, YE Kai-yuan. Torsion Theory of Elastic Cylinder[M].Beijing: Science Press, 1956: 187-206. (in Chinese)) [17] 王敏中, 王炜, 武际可. 弹性力学教程(修订版)[M].北京:北京大学出版社,2011: 127-174.(WANG Min-zhong, WANG Wei, WU Ji-ke.Course of Elastic Mechanics(Revised Edition)[M]. Beijing: Beijing University Press, 2011: 127-174. (in Chinese)) [18] 钱伟长, 叶开沅.弹性力学[M]. 北京:科学出版社,1980: 182183.(CHIEN Wei-zang, YE Kai-yuan. Elastic Mechanics[M]. Beijing: Science Press, 1980: 182-183.(in Chinese)) [19] 吴家龙.弹性力学[M]. 北京:高等教育出版社,2001: 215244.(WU Jia-long. Elasticity[M]. Beijing:Higher Education Press, 2001: 215244.(in Chinese)) [20] 杜柏松.考虑非线性影响的分体双箱梁悬索桥动力特性研究[D]. 博士论文.上海: 同济大学, 2006.(DU Bai-song. Study on dynamic characteristics of suspension bridges with twin-box stiffening girders considering geometric nonlinearity[D]. Ph D Dissertions. Shanghai:Tongji University, 2006. (in Chinese))
点击查看大图
计量
- 文章访问数: 1729
- HTML全文浏览量: 115
- PDF下载量: 1518
- 被引次数: 0