高层建筑有限元计算软件主要环节改进研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
随着建筑结构的不断发展,结构的计算分析工作越来越复杂,面对计算量巨大、设计效率要求又高的现状,高性能的有限元计算分析软件的研究越来越重要。本文针对建筑结构有限元计算中的几个关键技术问题进行了系统研究,得出一些切实可行的改进方案,部分已经应用于建筑结构分析有限元软件的开发之中。论文的主要研究内容如下:
     1)研究了有限元计算分析软件的发展历史与进展,概括归纳了当今国内外主流的分析与设计软件的特点、使用领域。在阅读国内外相关文献后,归纳总结了建筑结构有限元分析计算时遇到的一些常见问题,并进行了详细评述。
     2)研究了网格质量的评价指标。分析与总结了通用有限元软件ANSYS与SAP2000中的单个网格剖分评价指标,在现有单个网格评价指标计算的基础上,提出了一种新的网格评价指标,该指标能综合考虑单元的长细度、平行偏离程度等;不仅能用于单个网格的评价,而且能用于多个网格组成的有限元整体计算模型,数值计算结果表明该评价指标是高效而又可靠的。
     3)研究了框架梁与剪力墙连接有限元计算的处理方法。在ALLMAN型带转角自由度膜单元的基础上,通过罚函数法,将独立插值的转角位移场与协调位移场变形梯度的反对称部分建立相等关系并引入单元的泛函当中,构造了改进的带有转角自由度的膜单元。分别采用本文改进单元与梁式约束等方法对双肢剪力墙结构进行对比计算,数值结果表明,本文改进单元具有精度高、数据准备简单,罚参数物理意义明确等优点。
     4)通过ANSYS提供的二次开发UPFs功能,将改进后单元植入标准分析模块,对单元进行分片检验等测试。测试结果表明,单元能通过各种小片检验,无多余零能模式。各种测题中的结果均显示改进后的单元是一种可靠、性能优越的单元。
     5)研究了建筑结构中偏心构件的连接处理方法。经过数值算例的试算与对比,建议以基于主从约束的“偏心刚域”方式处理,无论构件长短、粗细,均能准确的计算出构件的位移与内力,从而避免了“刚性”杆件法计算结果与杆件“刚性”相关导致计算内力失真的问题。
     6)研究了有限元误差估计方法。采用二分加密与自适应计算两种方式,将基于能量误差的误差估计方法应用于简支梁、门形结构等的精度评价,都得出当能量误差小于10%时,位移、内力即能趋于稳定,并针对洞角存在应力集中问题,建议将洞口角点处小范围的面积排除分析域后,采用区域自适应方式进行分析,从而将自适应分析方法推广到高层建筑结构分析之中。并指出在计算资源受限的情况下,对连梁进行加密是一种最为行之有效的提高精度方式,其次是最底层柱的根部。
     综上,本文针对建筑结构有限元计算中的计算效率和计算精度问题进行了研究,提出了可行的改进方法,为推动高效、可靠的建筑结构有限元计算软件发展起到了一定的作用。
With the continuous development of structures, structural analysis becomes increasingly complex. Facing the huge amount of computation, high efficiency design requirements, the research of high-performance structure finite element analysis software(SFEAS) is increasingly important.In this paper, several key issues of SFEAS are systematicly studied and some practical improvements are given, parts of them have been applied to the software develop.The major research are as follows:
     1) Research on the history and progress of finite element analysis software development.The features of today's domestic and international mainstream analysis and design software are summed up.After reading the relevant literature, the problems encountered with structural finite element analysis are summarized, and the detailed commentary is made.
     2) Research on the evaluation index of mesh quality.After the analysis and summary of the single evaluation index of the commercial finite element software ANSYS and SAP2000, a new evaluation index, which could be used by both the single grid evaluation and multi-grid finite element model, is recommended, The results show that the evaluation index is efficient and reliable.
     3) Research on the finite element solution way of wall-beam connection.With the penalty function method, after the constraint relation between the independent rotation field interpolation and the deviatoric strain tensor part of displacement field was established, a high-precision membrane with drilling-dof was constructed.
     4)Using the user programming featrures(UPFs) interface of ANSYS, the membrane was implanted. The results show that it can pass most of patch test with no extra zero energy mode.A variety of patch tests show it is a reliable, superior element.In addition, it has a simple data preparation and penalty parameter has a reasonable definition.
     5)Research on the finite element solution way to connect eccentric components.Through numerical experiments, the "eccentric rigid domain", which based on the master-slave constraint,was proposed. Regardless of component length, thickness, it can accurately calculate the displacement and force, which avoids the wrong results by the "rigid" bars method because the internal forces is dependent on the stiffness of the "rigid" bars.
     6) Research on the finite element error estimation methods.To improve accuracy,with the use of dichotomy and adaptive mesh method, the energy error estimation method is applied to the accuracy of evaluation of the simply supported beams, door-shaped structure calculation. The result shows when the energy error is less than 10%,the displacement and internal force will stabilize. In order to solve the stress concentration,a small regeion near the hole conner is excluded from the adaptive analysis domain,by which can extends the adaptive method to the analysis of high-rise building. When the calculation capabiltiy is limited,in order to get higher acurrancy,a refined mesh is proposed to use in the connect wall-beams.
     In summary, the paper has studied some key issues of sturcture finite element analysis which focus on the efficiency and accuracy problems. Some suggustions are given, which may promote the efficient and reliable structural finite element analysis software.
引文
[1-1]中国建筑科学研究院.第二十一届全国高层建筑结构学术交流会论文集.南京:2010年.
    [1-2]手册编写组.建筑结构静力计算手册(第二版).北京:中国建筑工业出版社,1998年.
    [1-3]中国建筑科学研究院主编.高层建筑混凝土结构技术规程(JGJ3-2002).北京:中国建筑工业出版社,2002年.
    [1-4]R W Clough. The Finite Element method in plane stress analysis.Proc.2nd Conference on Electronic Computation. Pittsburg:ASCE, Sept.8-9,1960.
    [1-5]庄茁译.连续体和结构的非线性有限元.北京:清华大学出版社,2002年.
    [1-6]石亦平,周玉蓉ABAQUS有限元分析实例详解.北京:机械工业出版社2007年.
    [1-7]庄茁译ABAQUS有限元软件6.4版入门指南.北京:清华大学出版社,2004年.
    [1-8]陈火红,杨剑,薛小香,王朋波.新编Marc有限元实例教程.北京:机械工业出版社,2007年.
    [1-9]尚晓江,苏建宇LS-DYNA动力分析方法与工程实例.北京:中国水利水电出版社,2006年.
    [1-10]王呼佳,陈洪军等Ansys工程分析进阶实例.北京:中国水利水电出版社,2006年.
    [1-11]王新敏著ANSYS工程结构数值分析.北京:人民交通出版社,2007年.
    [1-12]张吉.基于概率灵敏度技术的钢结构螺栓节点性能研究.硕士学位论文.2007年.
    [1-13]江见鲸,陆新征,叶列平.混凝土结构有限元分析.北京:清华大学出版社,2004年.
    [1-14]岳戈ADINA流体与流固耦合功能的高级应用.北京:人民交通出版社,2010年.
    [1-15]北京金土木软件技术有限公司.SAP2000中文版使用指南.北京:人民交通出版社,2006年.
    [1-16]赵伟封等.工程结构数值分析及Algor实现.北京:中国电力出版社.2003年.
    [1-17]北京金土木软件技术有限公司.中国建筑标准设计研究院编著ETABS中文版使用指南.北京:中国建筑工业出版社.2004年
    [1-18]中仿科技主页http://www.cntech.com.cn/
    [1-19]MSC中国:http://www.mscsoftware.com.cn
    [1-20]安世亚太网址
    [1-21]GB50011-2001.建筑抗震设计规范.北京:中国建筑工业出版社,2001.
    [2-1]李忠学.梁板壳有限单元中的各种闭锁现象及解决方法.浙江大学学报(工学版),2007(7):1119-1125
    [2-2]E.L.Wilson, R.L.Taylor, W.P.Doherty, and J.Ghaboussi, Incompatible displacement models, Numerical and Computer Methods in Structural Mechanics(Edited by S.J.Fenves, N.Perrone, A.R. Robinson and W. C.Schnobrich). Academic Press, New York pp:43-57. (1973)
    [2-3]E.L. Wilson, A.Ibrahimbegovic. Use of incompatible displacement modes for the calculation of element stiffness or stresses. Finite Elements in A nalysis and Design,1990(7):229-241.
    [2-4]G. P. Bazeley, Y. K. Cheung, B. M. Irons and O. C. Zienkiewicz, Triangular elements in plate bending. Conforming and nonconforming solutions, Proc.1st Conf. Matrix Methods in Structural Mechanics,AFFDL-TR-CC-80, Wright Patterson A. F. Base, Ohio,1966, pp:547-576.
    [2-5]B. M. Irons and A. Razzaque, Experience with the patch test for convergence of finite element methods, Math. Foundations of the Finite Element Method (Ed. A. K. Aziz), Academic Press,1972, pp:557-587.
    [2-6]G. Strang, Variational crimes and the finite element method, Proc. Foundations of the Finite Element Method(Ed. A. K. Aziz), Academic Press,1972, pp:689-710.
    [2-7]E. R. de Arantes e Oliveira, The patch test and the general convergence criteria of the finite element method, Int. J. Solids Struct.,1977 (13):159-178.
    [2-8]R.H.MacNeal. Finite Elements:Their Design and Performance,Marcel Dekker,New York,1994.
    [2-9]B. Fraeijs de Veubeke, Variational principles and the patch test, International Journal for Numerical Methods in Engineering,1974(8):783-801.
    [2-10]R.L.Taylor, J.C.Simo, O.C.Zienkiewicz etc. The Patch Test—A Condition for Assessing FEM Convergence. International Journal for Numerical Methods in Engineering, 1986(22), pp:39-62
    [2-11]李文华,赵彦革,张吉等.基于ANSYS的钢筋混凝土浅梁单调加载变形的数值模拟.土木建筑工程信息技术2010(2):83-90.
    [2-12]中华人民共和国水利部.水工混凝土结构设计规范SL/T191-96.1997年.
    [2-13]R.D. Cook. Concepts and Applications of Finite Element Analysis.Wiley and Sons。 New York.2000
    [2-14]吴晓涵.面向对象结构分析程序设计.北京:科学出版社.2002年.
    [2-15]Cardoso RPR, Yoon JW, Valente RAF. A new approach to reduce membrane and transverse shear locking for one-point quadrature shell elements:linear formulation. International Journal for Numerical Methods in Engineering 2006(2):214-249
    [2-16]Huges T J, Cohen M, Haroun M. Reduced and selective integration technology in finite element analysis of plates. International Journal for Numerical Methods in Engineering 1978(46):203-222.
    [2-17]吴学敏.高层建筑剪力墙中连梁设计的探讨[J].建筑结构学报.1995.16(1):77-78
    [3-1]ANSYS help system-meshing.ANSYS Corporation.
    [3-2]J. A. Stricklin, W. S. Ho, E. Q. Richardson and W. E. Haisler, On isoparametric vs linear strain triangular elements,International Journal for Numerical Methods in Engineering,1977 (11):1041-1043.
    [3-3]L.N. Gifford, More on distorted isoparametric elements, International Journal for Numerical Methods in Engineering,1979 (14):290-291.
    [3-4]J. Robinson,Some new distortion measures for quadrilaterals. Finite Elements in Analysis and Design,1987 (3):183-197.
    [3-5]NAM-SUA LEE, KLAUS-JURGEN BATHE. Effects of element distortions on the performance of isoparametric elements. International journal for numerical methods in engineering.1993(36):3553-3576.
    [3-6]曾攀.有限元分析及应用.北京:清华大学出版社.2004年.
    [3-7]关振群,宋超等.有限元网格生成方法研究的新进展.计算机辅助设计与图形学报,2003(1):1-14.
    [3-8]K Hole. Finite element mesh generation methods:a review and classification. Computer Aided Design,1988(20):27-38.
    [3-9]S.TIMOSHENKO J.N.GOODIER. Theory of elasticity. New York, Toronto, London: McGraw-Hill book company, Inc.,(1951).
    [3-10]R. H. MacNeal, R. L. Harder, A Proposed Standard Set of Problems to Test Finite Element Accuracy, Proceedings,25th SDM Finite Element Validation Forum,(1984).
    [3-1l]中国建筑科学研究院PKPM CAD工程部PMSAP复杂多、高层建筑结构分析与设计软件.用户手册及技术条件.
    [3-12]广东省建筑设计研究院.建筑结构通用分析与设计软件GSSAP说明书,2010年.
    [3-13]王勖成,邵敏.有限单元法基本原理和数值方法.北京:清华大学出版社,1997.
    [3-14]黄晓东,杜群贵,叶邦彦.二维几何特征自适应有限元网格生成(一)-几何特征识别.计算机辅助设计与图形学学报.2004(7):923-928.
    [3-15]杜群贵,刘邵宏,黄晓东.二维几何特征自适应有限元网格生成(二)-算法描述,计算机辅助设计与图形学学报.2004(7):923-928.
    [3-16]傅沛福,吴淑芳,胡平等.全自动自适应网格细化.计算力学学报.1997(2):114-118.
    [3-17]李梅娥,周进雄.任意平面区域的高质量有限元网格自动剖分.应用力学学报.2004(9):125-128.
    [3-18]Ted D. Blacker, Michael B.Stephenson.Paving:A new approach to automated quadrilateral mesh generation. International Journal for Numerical Methods in Engineering,1991(32):811-847.
    [3-19]K Hole. Finite element mesh generation methods:a review and classification. Computer Aided Design,1988(20):27-38.
    [3-20]李华,李笑天,程耿东等.一种全四边形网格生成方法-改进模板法.计算力学学报.2002(2):16-19.
    [3-21]W.A. Cook and W.R. Oakes. Mapping methods for generating three-dimensional meshes.Comp. Mech. Eng.,1983(1):67-72.
    [3-22]李庆杨,关治,白峰山.数值计算原理.北京:清华大学出版社,2000年.
    [3-23]ANSYS Help System. Verification Manual for the Mechanical APDL Application-Benchmark D3.
    [4-1]陈万吉,李勇东.带旋转自由度的精化非协调平面四边形等参元.计算结构力学及其应用,1993(2):22-29.
    [4-2]胡宇峰,戚震华.带旋转自由度的膜单元的研究进展.上海铁道大学学报,1996(6):78-83.
    [4-3]刘红,唐立民等,带转角自由度的拟协调平面元,计算结构力学及其应用.1990(7):23-31.
    [4-4]J. Robinson.Four-node quadrilateral stress membrane element with rotational stiffness. International Journal for Numerical Methods in Engineering,2000(47):907-920.
    [4-5]Allman D J. A compatible triangular element including vertex rotations for plane elasticity analysis. Computers & structures.1984,Vol.19,1-8.
    [4-6]Cook R D. On the Allman triangular and a related quadrilateral element. Computers and structures,1986(22):1065-1067.
    [4-7]中国建筑科学研究院PKPM CAD工程部.TAT多层及高层建筑结构三维分析与设计软件(薄壁柱模型).用户手册及技术条件.2008年.
    [4-8]中国建筑科学研究院高层建筑技术开发部.多层及高层建筑结构空间分析程序TBSA计算分析手册.2002年.
    [4-9]中国建筑科学研究院建筑工程软件研究所SATWE用户手册及技术条件.2010年.
    [4-10]焦柯,吴文勇等.高层结构剪力墙计算若干问题探讨.建筑结构,2007(9):134-137.
    [4-11]王勖成.有限单元法基本原理和数值方法(第3版).北京:清华大学出版社,1997年.
    [4-12]关正西,强洪夫译.有限元分析的概念与应用.西安:西安交通大学出版社.2007年.
    [4-13]Edward L.Wilson.Three-Dimensional Static and Dynamic Analysis of Structures.
    [4-14]龙驭球,辛克贵.广义协调元,土木工程学报,1987(1):1-14.
    [4-15]龙志飞,须寅,龙驭球.具有旋转自由度的广义协调矩形膜元,工程力学,1992(5):68-75.
    [4-16]龙驭球,须寅.构造几何不敏感四边形膜元的广义协调方法.力学学报,1997(11):692-700.
    [4-17]须寅,龙驭球.采用广义协调条件构造具有旋转自由度的四边形膜元.工程力学,1993(3):27-35.
    [4-18]Kanok-Nukulchai. A Simple and Efficient Finite Element for General Shell Analysis. International Journal for Numerical Methods in Engineering.1979(14):179-200.
    [4-19]K. J. BATHE, LEE-WING HO.A simple and effective element for analysis of general shell structures.Computers and structures,1981 (13):673-681.
    [4-20]徐芝纶.弹性力学(第4版).北京:高等教育出版社,2006年.
    [4-21]A. Ibrahimbegovic and F. Frey, Membrane quadrilateral finite elements with rotational degrees of freedom, Engineering Fracture Mechanics,1992(43):13-24.
    [4-22]E. Reissner, A note on variational principles in elasticity, International Journal of Solids and Structures,1965(1):93-95.
    [4-23]朱伯芳.有限单元法原理和应用(第2版).北京:中国水利电力出版社,1998年.
    [4-24]F. Frey, Shell finite elements with six degrees of freedom per node. In Analytical and Computational Models for Shells. ASME, New York (1989):291-317.
    [4-25]R. L. Taylor and J. C. Simo, Bending and membrane elements for analysis of thick and thin shells. Proceedings NUMETA 85 (Edited by J. Middelton and G N. Pande). Pineridge Press, Swansea (1985):587-591.
    [4-26]N. Carpenter, H. Stolarski and T. Belytschko, A flat triangular shell element with improved membrane interpolation. Commun. Appli. Numer. Method. (1985)(1):61-168.
    [4-27]Jetteur and F. Frey, A four node Marguerre element for nonlinear shell analysis. Engng Compur.1986(3):276-282.
    [4-28]T.J.R Hughes.The Finite Element Method:Linear Static and Dynamic Analysis, Printice-HallEnglewood Cliffs,NJ,1987.
    [4-29]袁明武,孙树立.多层与高层建筑中墙的模型.工程力学增刊,1996:63-70.
    [4-30]北京金土木软件技术有限公司,中国建筑标准设计研究院译,结构静力与动力分析-强调地震工程学的物理方法[M].(第四版),北京:中国建筑工业出版社.2006年.
    [4-31]Long Yuqiu, Xu Yin. Generalized conforming quadrilateral membrane element with vertex rigid rotational freedoms. Computers and Structures,1994(52):74-75.
    [4-32]Yuan K.Y.,Huang Y S, Pian THH. New strategy for assumed stresses for 4-node hybrid stress membrane elements. International Journal for Numerical Methods in Engineering, 1992(36):1747-1763.
    [4-33]K.Y. Sze. On immunizing five-beta hybrid-stress element models form'trapezoidal locking' in practical analyses. International Journal for Numerical Methods in Engineering. 2000(47):907-920
    [4-34]A. Ibrahimbegovic, R. L. Taylor and E. L.Wilson, A robust membrane quadrilateral element with drilling degrees of freedom. International Journal for Numerical Methods in Engineering.1990(30):445-457.
    [4-35]D. J. Allman, A quadrilateral finite element including vertex rotations for plane elasticity problems, International Journal for Numerical Methods in Engineering,1988 (26):717-739.
    [4-36]T.H.H. Pian and W. Changchun. General formulation of incompatible shape function and an incompatible isoparametric element. In:Proc. of the Invitational China-American Workshop on FEM. Chengda:1986.
    [4-37]Long YQ, Huang MF. A generalized conforming isoparametric element. Applied Mathematics and Mechanics-English Edition 1988; 9(10):929-936.
    [4-38]岑松,龙志飞,龙驭球.对转角和剪应变场进行合理插值的厚薄板通用四边形单元.工程力学,1999(8):1-15.
    [5-1]朱恒伟,杨光松,方赤峰.组合结构罚有限元分析与罚因子的影响.固体力学学报,1996(3)77-80.
    [5-2]袁明武,孙树立等.一种新的墙单元.计算结构力学及其应用,1996(2):17-24.
    [5-3]钱伟长.变分法及有限元.北京:科学出版社.1980年.
    [5-4]M. M. Denn. Optimization by Variational Methods. McGraw-Hill. New York.1969.
    [5-5]朱合华,谬圆冰等.组合结构有限元计算存在的问题和处理方法.岩土力学.2005(9):1437-1441
    [5-6]马文华,吴新炳.约束罚单元及其在接触问题中的应用.上海力学,1984(4):403-415
    [5-7]GallagherR.H.andZienkiewiezO.C.,optimum StrueturalDesign—Theory and APPlieations, John Wiley and SonsLtd.1973.
    [5-8]陈瑜,龚炳年等译.高层建筑结构分析与设计.北京:地震出版社,1993年.
    [5-9]梁启智.高层建筑结构分析与设计.广州:华南理工大学出版社.1991年
    [5-10]曲圣年.组合结构有限元分析的罚单元法.周体力学学报,1982(4):555—559.
    [5-11]J.F.Abel and M.S.Shephard,An Algorithm for Multipoint Constraints in Finite Element Analysis, International Journal for Numerical Methods in Engineering,1979 (14):464-467.
    [5-12]S.C.Chang and T.W.Lin, Constraint Relation Implementation for Finite Element Analysis from an Element Basis.Advances in Engineering Software,1988 (10):191-194.
    [6-1]K.J. Bathe.Finite Element Procedures.Prentice Hall.Upper Saddle River,New Jersey.1996
    [6-2]Oliverira E R A. Optimization of finite element solutions. Proceedings of the 3rd conference on matrix methods in structural mechanics, Wright-Patterson, Ohio,1971.
    [6-3]郭书祥.自适应有限元方法及其工程应用.力学进展.1997(11):479-488.
    [6-4]杨强,吴浩,周维垣.h-型自适应有限元法分析中的大坝应力取值标准.水利学报,2005(3):1-10.
    [6-5]中国建筑科学研究院建筑工程软件研究所PKPM多高层结构计算软件应用指南.北京:中国建筑工业出版社.2010年.
    [6-6]Z. Uthman and H. Askes. R-adaptive and h-adaptive mesh optimisation based on error assessment,13th UK National ACME conference, ed. R.S. Crouch, University of Sheffield, (2005):127-130.
    [6-7]王辉明,马希龄,赵文.自适应有限元方法及其在ANSYS软件中的应用.新疆大学学报,2003(2):84-87.
    [6-8]Zienkiewicz O C, Zhu J Z. A simple error estimator and adaptive procedure for practical engineering analysis. International Journal for Numerical Methods in Engineering,1987(24): 337-357.
    [6-9]ANSYS help system-Advanced Analysis Techniques.
    [6-10]于爱民,王勖成.有限元计算中疏密网格间过渡单元的构造.清华大学学报,1999(8):100-103.
    [6-11]王建华,杨磊,沈为平.有限元后验误差估计方法的研究进展.力学进展.2000(5):175-190.
    [6-12]J. E. Akin.Finite Element Analysis with Error Estimators An Introduction to the FEM and Adaptive Error Analysis forEngineering Students. Elsevier Butterworth-Heinemann Linacre House, Jordan Hill, Oxford OX28DP30 Corporate Drive, Burlington, MA 01803.2005.
    [6-13]C. Carstensen, M. Maischak, E. Stephan, A posteriori error estimate and h-adaptive algo-rithm on surfaces for Symm's integral equation, Numer. Math.90 (2001):197-213.
    [6-14]C. Carstensen, M. Maischak, D. Praetorius, E. Stephan, Residual-based a posteriori errorestimate for hypersingular equation on surfaces, Numer. Math.97 (2004):397-425.
    [6-15]C. Carstensen, D. Praetorius, Averaging techniques for the effective numerical solution of Symm's integral equation of the first kind, SIAM J. Sci. Comput.27 (2006):1226-1260.
    [6-16]P. Binev, W. Dahmen, R. DeVore, Adaptive finite element methods with convergence rates,Numer. Math.97 (2004):219-268.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700