基于二次开发的PC桥梁三维仿真分析关键技术研究
|
摘要
目前,桥梁结构的整体分析一般都是建立在平面或空间杆系有限元理论基础上,部分采用三维实体单元分析但尚不成熟。对于普通钢筋混凝土和简单预应力混凝土桥梁,结构计算采用杆系有限元加局部三维分析基本够用。但对于特殊造型桥梁、多向预应力体系桥梁及弯、宽桥,这些以杆系有限元理论为基础的分析方法已不能完全满足设计者的需要。随着大跨度和复杂的桥型不断涌现、计算机技术的高速发展、有限元理论的逐步完善,全桥结构三维仿真分析(SSEB)已进入实质性研发阶段。通过三维实体结构仿真分析,模拟桥梁构件的空间状态,真实反映边界条件和荷载状况,可得到更详尽、精确和可靠的分析结果,以指导设计和施工。
本文基于大型通用有限元软件ANSYS和MATLAB,探索三维实体模型下进行PC桥梁结构整体仿真分析的途径,研究其关键理论及专业模块开发技术。重点研究了预应力筋的模拟、预应力效应的等效转换和桥梁截面切片技术。研究成果可广泛应用于PC桥梁的三维仿真领域,具有较大的理论意义和工程实用价值。
本文在如下几个方面进行了理论、数值研究,并取得了一些研究成果:
1、根据PC桥梁设计与实际筋施工中预应力布置的特点,成功引入CAGD中的非均匀B-Spline模拟空间曲线型预应力筋,成为课题研究的理论前提。
2、研究了桥梁结构中预应力损失产生机理,推导出任意空间曲线型预应力筋的摩擦损失和考虑反摩阻力影响的锚具回缩损失的数值公式;
3、根据预应力等代荷载原理,实现了预应力效应,即锚固端集中力和体内预应力径向力在实体单元中的等效转换,该方法对各类实体单元具有通用性;
4、鉴于桥梁截面多样化特点,提出一套约束Delaunay三角网格划分方法,并针对实体模型,研究了截面几何特性、截面内力和正截面抗弯承载能力计算理论。
5、实现了实体模型下的截面切片可视化以及相关配套技术,为桥梁精细化分析和处理大容量结果分析数据提供了一条可选途径。
6、基于WINDOWS消息处理机制,实现了对ANSYS软件的无缝控制。采用Fortran语言和C++面向对象编程,开发出基于ANSYS与MATLAB平台的PC桥梁分析模块,实现了实体模型下的曲线PC箱梁桥的静力分析。
由于以ANSYS和MATLAB平台为依托,上述成果的实现只编写了少量程序代码,却能利用ANSYS前、后处理、求解器和MATLAB C++数学库等强大功能,说明本文尝试的软件开发途径是有效的。
本文的研究对于提高PC桥梁结构设计水平和效率,普及三维实体有限元技术的应用,拓展专业模块软件的开发途径具有重要的理论意义和实用价值。
Currently, most of the analyses on the structural simulation for entire bridge are based on the theory of 2D or 3D beam Finite Element Method (FEM), whereas solid element method is seldom to be based on due to its immaturity. It is feasible to use beam FEM combined with partial solid FEM for RC bridges or simple PC bridges simulation, but when it comes to the special shape bridge and milt-directional curved prestressed bridge, such theories cannot fully meet the needs of designers. With the rapid increase in large span and complicate bridge forms, as well as the development in computer technology and FEM theory, the Structural 3D Simulation for Entire Bridge (SSEB) has stepped onto its substantial stage of research and development. By means of SSEB, boundary condition and load situation can be real reflected, and more accurate and reliable analysis results can be obtained to assist design and construction.
Based on FEM software ANSYS and MATLAB, the simulation analysis approach for 3D entity model of PC bridge structures and relevant key theories as well as program developing strategies were studied in this paper. The research is focused on the line shape prestressing tendon simulation, equal translation of prestressed effect and cross-section slice technology. The research results, which can be widely applied in 3D simulation of PC bridges, will have its theoretical and practical significance.
Theoretical and numerical researches in this paper are carried out in the following aspects, with some achievements accomplished:
1. Considering the characteristics of prestressing tendon layout in bridge design and construction, non-uniform B-Spline in CAGD was successfully introduced into simulating 3D curved prestressing tendon, and it became the theoretical premise of this research as such.
2. A numerical formula to calculate the prestressed losses due to friction and the pipeline local defection of 3D curved prestressed tendons was deduced. Also, with consideration to unfriction, another formula to calculate the anchorage losses was deduced in this paper.
3. According to prestress equipollent load theory, prestressed effect was realized. That means the concentrated force on the fixed end of anchorage and the internal prestressing radial force could be transformed equivalently on entity elements. This method can be applied in various solid elements.
4. In view to the diversity of bridge sections, a restraint Delaunay triangle net division method was brought up. Theory of section geometrical feature, section endogen force and right section bending resistance capacity calculation were studied to solid modal.
5. The visualization technique of cross-section slice and matching auxiliary techniques for entity model were put forward and realized by programming, and provides a method for accurate analysis and analysis with high-capacity results data.
6. Based on the message-handling mechanism of Windows, the command line of ANSYS was successfully intercepted. By applying FORTRAN language and C++ object-oriented programming, a professional analysis modular with PC structure handling ability based on ANSYS and MATLAB was explored. Therefore, the static analysis for PC box girder on entity model could be carried out.
Only a few codes were compiled for accomplishing the above achievements, but powerful capabilities of ANSYS such as pre/post processors, solvers, and MATLAB C++ mathematical function library could be used fully. Since the task was built upon the ANSYS and MATLAB platform, it means the attempt of software development in this paper was rather effective.
The paper has profound theoretical and practical significance in improving the design level and efficiency of PC bridges, popularizing the application of finite element method (FEM), and promoting the development of professional software.
本文基于大型通用有限元软件ANSYS和MATLAB,探索三维实体模型下进行PC桥梁结构整体仿真分析的途径,研究其关键理论及专业模块开发技术。重点研究了预应力筋的模拟、预应力效应的等效转换和桥梁截面切片技术。研究成果可广泛应用于PC桥梁的三维仿真领域,具有较大的理论意义和工程实用价值。
本文在如下几个方面进行了理论、数值研究,并取得了一些研究成果:
1、根据PC桥梁设计与实际筋施工中预应力布置的特点,成功引入CAGD中的非均匀B-Spline模拟空间曲线型预应力筋,成为课题研究的理论前提。
2、研究了桥梁结构中预应力损失产生机理,推导出任意空间曲线型预应力筋的摩擦损失和考虑反摩阻力影响的锚具回缩损失的数值公式;
3、根据预应力等代荷载原理,实现了预应力效应,即锚固端集中力和体内预应力径向力在实体单元中的等效转换,该方法对各类实体单元具有通用性;
4、鉴于桥梁截面多样化特点,提出一套约束Delaunay三角网格划分方法,并针对实体模型,研究了截面几何特性、截面内力和正截面抗弯承载能力计算理论。
5、实现了实体模型下的截面切片可视化以及相关配套技术,为桥梁精细化分析和处理大容量结果分析数据提供了一条可选途径。
6、基于WINDOWS消息处理机制,实现了对ANSYS软件的无缝控制。采用Fortran语言和C++面向对象编程,开发出基于ANSYS与MATLAB平台的PC桥梁分析模块,实现了实体模型下的曲线PC箱梁桥的静力分析。
由于以ANSYS和MATLAB平台为依托,上述成果的实现只编写了少量程序代码,却能利用ANSYS前、后处理、求解器和MATLAB C++数学库等强大功能,说明本文尝试的软件开发途径是有效的。
本文的研究对于提高PC桥梁结构设计水平和效率,普及三维实体有限元技术的应用,拓展专业模块软件的开发途径具有重要的理论意义和实用价值。
Currently, most of the analyses on the structural simulation for entire bridge are based on the theory of 2D or 3D beam Finite Element Method (FEM), whereas solid element method is seldom to be based on due to its immaturity. It is feasible to use beam FEM combined with partial solid FEM for RC bridges or simple PC bridges simulation, but when it comes to the special shape bridge and milt-directional curved prestressed bridge, such theories cannot fully meet the needs of designers. With the rapid increase in large span and complicate bridge forms, as well as the development in computer technology and FEM theory, the Structural 3D Simulation for Entire Bridge (SSEB) has stepped onto its substantial stage of research and development. By means of SSEB, boundary condition and load situation can be real reflected, and more accurate and reliable analysis results can be obtained to assist design and construction.
Based on FEM software ANSYS and MATLAB, the simulation analysis approach for 3D entity model of PC bridge structures and relevant key theories as well as program developing strategies were studied in this paper. The research is focused on the line shape prestressing tendon simulation, equal translation of prestressed effect and cross-section slice technology. The research results, which can be widely applied in 3D simulation of PC bridges, will have its theoretical and practical significance.
Theoretical and numerical researches in this paper are carried out in the following aspects, with some achievements accomplished:
1. Considering the characteristics of prestressing tendon layout in bridge design and construction, non-uniform B-Spline in CAGD was successfully introduced into simulating 3D curved prestressing tendon, and it became the theoretical premise of this research as such.
2. A numerical formula to calculate the prestressed losses due to friction and the pipeline local defection of 3D curved prestressed tendons was deduced. Also, with consideration to unfriction, another formula to calculate the anchorage losses was deduced in this paper.
3. According to prestress equipollent load theory, prestressed effect was realized. That means the concentrated force on the fixed end of anchorage and the internal prestressing radial force could be transformed equivalently on entity elements. This method can be applied in various solid elements.
4. In view to the diversity of bridge sections, a restraint Delaunay triangle net division method was brought up. Theory of section geometrical feature, section endogen force and right section bending resistance capacity calculation were studied to solid modal.
5. The visualization technique of cross-section slice and matching auxiliary techniques for entity model were put forward and realized by programming, and provides a method for accurate analysis and analysis with high-capacity results data.
6. Based on the message-handling mechanism of Windows, the command line of ANSYS was successfully intercepted. By applying FORTRAN language and C++ object-oriented programming, a professional analysis modular with PC structure handling ability based on ANSYS and MATLAB was explored. Therefore, the static analysis for PC box girder on entity model could be carried out.
Only a few codes were compiled for accomplishing the above achievements, but powerful capabilities of ANSYS such as pre/post processors, solvers, and MATLAB C++ mathematical function library could be used fully. Since the task was built upon the ANSYS and MATLAB platform, it means the attempt of software development in this paper was rather effective.
The paper has profound theoretical and practical significance in improving the design level and efficiency of PC bridges, popularizing the application of finite element method (FEM), and promoting the development of professional software.
引文
[1]毕桂平,魏红一,范立础.ANSYS在桥梁工程中的应用前景[A].中国土木工程学会桥梁及结构工程学会第十四届年会论文集[C].上海:同济大学出版社,2000:410-412.
[2]毕桂平,魏红一,范立础.ANSYS在桥梁工程中的应用前景[A].中国土木工程学会桥梁及结构工程学会第十四届年会论文集[C].上海:同济大学出版社,2000:410-412.
[3]江见鲸,贺小岗.工程结构计算机仿真分析[M].北京:清华大学出版社,1995.
[4]周小勇,金文成.洛河大桥预应力连续梁V型0#节段结构构造及受力特性分析[J].公路交通科技(应用技术版),2007(4):145-149.
[5]顾祥林,孙飞.混凝土结构的计算机仿真[M].上海:同济大学出版社,2002.
[6]程纬.基于ANSYS二次开发的钢筋混凝土结构倒塌仿真分析及其程序设计[博士后报告].上海:同济大学,2002.
[7]喻定新.以ANSYS为分析器的结构优化软件研制—杆板组合结构优化[硕士学位论文].南宁:广西大学,2004.
[8]段志东.ANSYS二次开发途径研究[硕士学位论文].兰州:兰州铁道学院,2002.
[9]魏志刚,汤文成,刘德仿.有限元分析在结构设计中的应用[J].设计与研究,2005,32(8):12-14.
[10]李正,徐立功.ANSYS二次开发在桥梁设计中的应用[J].江苏交通工程,2003,115(6):22-25.
[11]刘艳萍,杨新华,杨文兵.预应力钢筋混凝土局部有限元分析的ANSYS二次开发[J].华中科技大学学报(城市科学版),2005,22(增刊):87-90.
[12]杨建明.预应力混凝土结构的当量荷载计算[J].建筑结构学报,1997,18(5):2-11.
[13]ANSYS theory reference Release 5.7[M].ANSYS Inc.1994.
[14]Analysis of Geotechnical Problems with ABAQUS[M].ABAQUS inc.2001.
[15]ADINA Theory and Modeling Guide[M].ADINA inc.2002.
[16]About ANSYS/CivilFEM[EB/OL].http://www.civilfem.com/about_us/about_us.php,2008.
[17]桥梁三维预应力空间分析系统软件帮助文件(bridgeKF1.0)[M].中华人民共和国交通部科研所,2007-2.
[18]Guide to ANSYS user programmable features[M].ANSYS Inc.1998.
[19]Duane D,Bruce L.Mastering Matlab 7[M].London:Prentice Hall,2004-10.
[20]Compaq Visual FORTRAN Programmer's Guide[M].Houston:Compaq Computer Corporation,2001-8.
[21]Beck Z.Microsoft Visual C++ 6.0 Programmer's Guide[M].北京:北京希望电子出版社,1998-12.
[22]MSC.Nastran简介[EB/OL].http://www.baisi.net/viewthread.php,2008.
[23]John O H.LS_DYNA Theoretical Manual[M].Livermore Software Cor.1997.
[24]CGAL User and Reference Manual Table of Contents[EB/OL].http://www.cgal.org/Manual/3.3/doc_html/cgal_manual/contents.html.2007.
[25]Spatial公司主页(中文版)[EB/OL].http://www.spatial.com.cn/community/case.html,2008.
[26]史忠植,罗杰文,王茂光,等.互联网智能平台AGrIP构建及其应用[A].信息技术快报,2005,3(5):24-30.
[27]王永岩.有限元法基础与SAP5程序应用[M].北京:煤炭工业出版社,1995.
[28]钟万勰.有限元分析与优化设计软件JIFEX研究开发和应用[EB/OL].http://sail.dlut.edu.cn/chengguo/04.htm,2003.
[29]通用有限元软件Strat用户手册[M].上海交通大学土木系,上海佳构软件科技有限公司,2003-5.
[30]施法中.计算机辅助几何设计与非均匀有理B样条[M].北京:高等教育出版社,2001.
[31]Boor C D.On calculating with B-Splines[J].J Approx.Theory,1972,6(1):50-62.
[32]金文成,周小勇.20节点等参单元的预应力等效节点力计算[J].华中科技大学(城市科学版),2008,25(1):5-7.
[33]张立明.Algor＼Ansys在桥梁工程中的应用方法与实例[M].北京:人民交通出版社,2003.
[34]吕勇刚.CAGD自由曲线曲面造型中均匀样条的研究[博士学位论文][D].杭州:浙江大学,2002.
[35]Sablonniere P.Spline and Bezier polygons associated with a polynomial spline curve[J].Computer Aided Design,1978,10(4):257-261.
[36]Barnhill R E,Riesenfield R F.Computer Aided Geometric Design[M],New York:Academic Press,1974.
[37]Cox M G.The Numerical elvaluation of B-Splines[J].J.Inst.Math.Applications,1972.
[38]Yamaguchi F,Curves and surfaces in computer aided geometric design[M].Berlin:Springer-Verlag,1988.
[39]Ferguson J.Multivariable curve interpolation[J].ACM,1964,11(2):221-228.
[40]Coons S.Rational bi-cubic surface patches[R].Boston:Massachusetts Institute of Technology Project MAC,1968.
[41]Schoenberg I.Spline functions and the problem of graduation[J].Proceedings of the National Academy of Sciences of USA,1964(52):947-950.
[42]Bezier P E.Numerical control-Mathematics and applications[M].London:Forrest Tran,1972:16-53.
[43]Boor D C.On Calculation with B-Spline[J].Approx Theory,1972,6(1):50-62.
[44]Boehm W.Inserting new knots into B-Spline curves.Computer-Aided Design,1980,12(4):199-201.
[45]Forrest A R.Interactive interpolation and approximation by Bezier polynomials[J].The Computer,1972,15(1):71-79.Reprinted in CAD 1990,22(9):527-537.
[46]ISO 10303,Standard for the Exchange of Product Model Data(STEP)[S].
[47]BS EN 29592-2:1991,Computer graphics.Programmer's Hierarchical Interactive Graphics System(PHIGS)[S].
[48]李翠华,郑南宁,张永平.基于样条修匀公式的图像边缘检测[J].电子学报,1999,27(1):2-5.
[49]Prautzsch H.Degree elevation of B-Spline curves[J].Computer Aided Geometric Design,1(12):193-198,1984.
[50]Peng Q.An algorithm for finding the intersection lines between two B-Spline surfaces[J].Computer Aided Design,1984,16(4):191-196.
[51]戴公连,李德建.桥梁结构空间分析设计方法与应用[M].北京:人民交通出版社,2001.10
[52]Riesenfeld R F.Nonuniform B-Spline curves[C].2nd USA-JAPAN Confence proc.1975:551-555.
[53]Hartley D J,Judd C J.Parametrization and shape of B-Spline curves[J].Computer Aided Design,1980,12(5):235-238.
[54]Prautzseh H,Piper B.A fast algorithm to raise the degree of B-Spline curves[J].Computer Aided Geometric Design,1991,8(4):253-266.
[55]李翠华,郑南宁,张永平.基于样条修匀公式的图像边缘检测[J].电子学报,1999,27(1):2-5.
[56]JTG D62-2004,公路钢筋混凝土及预应力混凝土桥涵设计规范[S].北京:人民交通出版社,2004.
[57]李铁军,杨红梅.计算机辅助几何造型设计中曲面拼接的处理[J].辽宁工学院学报,1999,19(4):23-26.
[58]Farin G.Curvature continuity and offsets for piecewise conics,ACM Transactions on Graphics,1989,8(2):89
[59]李卫国,徐步青,胡克辉.空间多曲线钢索预应力摩擦损失试验研究[J].2000,13(1):14-17.
[60]南开大学数学系.空间解析几何引论[M].北京:人民教育出版社,1978-9.
[61]Boor C D.A Practical Guide to Spline[M].Berlin:Springer-Verlag,1978.
[62]李红.数值分析[M].武汉:华中科技大学出版社,2003-11.
[63]李厚源.空间解析几何[M].济南:山东科学技术出版社1983-11.
[64]MathWorks Inc.MATLAB C++ Math Library 2.0 Reference[EB/OL].http://www.thphys.may.ie/CompPhysics/matlab/help/extern/cppmath/cpp_cref.html,2008.
[65]Prautzsch H.Degree elevation of B-Spline curves[J].Computer Aided Geometric Design,1984,12(1):193-198.
[66]施法中.论参数曲线曲面的几何不变性[J].工程图学学报,1991(1):79-86.
[67]王天军.曲面造型中几何连续曲面的研究与应用[博士学位论文][D].北京:北京航天航空大学,1991.
[68]Sapidis N,Farin G.Automatic fairing algorithm for B-Spline curves[J].Computer Aided Design,1990,22(2):121-129.
[69]Stam J.On subdivision schemes generalizing uniform B-Spline surfaces of arbitrary degree[J].Computer Aided Geometric Design,2001,18(5):383-396.
[70]李律,钱济章,李敦.基于ANSYS的预应力筋数值模拟[J].公路工程,2007,32(4):178-179.
[71]金文成,张小飞,周小勇.斜拉桥塔索锚固区空间应力分析[J].公路交通科技,2007,24(9):69-73.
[72]MIDAS/Civil V6.7.1在线帮助手册[Z].北京:MIDAS Information Technology China Co.,Ltd,2005.
[73]周宗泽.桥梁博士系统V3.1版[M].同济大学桥梁工程系,2005-8.
[74]李国平.预应力混凝土结构设计原理[M].北京:人民交通出版社,2000:30-43.
[75]宋玉普,车轶,马德有,等.空间多曲线型预应力钢索的预应力摩擦损失研究[J].土木工程学报,2002,35(6):105-108.
[76]刘庆文,苏小卒,空间曲线型预应力摩擦损失计算[J].四川建筑科学研究,2007,33(1):42-44.
[77]胡狄,陈政清.考虑反向摩阻的后张法PC构件锚固损失的计算[J].中国公路学报,2004,17(1):34-38.
[78]Hung T.Anchorage take-up loss in post- tensioned members[J].Journal of the Prestressed Concrete Institute,1969,14(4):30-35.
[79]程万年.圆弧预应力筋锚固损失基本方程和解[J].力学与实践,1992,14(3):33-35.
[80]胡光祥,陈东杰,赵勇,等.空间曲线预应力束摩擦损失测试和摩擦系数反演[J].建筑科学,2004,20(5):26-29.
[81]叶建.下白石特大桥空间曲线孔道摩擦损失分析及测试[J].福建建筑,2005,93(3):115-116.
[82]刘志文,宋一凡,赵小星,等.空间曲线预应力束摩阻损失参数[J].西安公路交通大学学报,2001,21(3):42-44.
[83]李卫国,徐步青,胡克辉.空间多曲线钢索预应力摩擦损失试验研究[J].石家庄铁道学院学报,2000,13(1):14-17.
[84]周小勇,金文成.独塔单索面斜拉桥-施州大桥设计和施工[J].中国市政工程,2006,124(6):23-25,88.
[85]刘庆文,苏小卒.后张法预应力空间多曲线钢筋锚固损失简化计算[J].铁道建筑,2006(2):10-12.
[86]杨建明.预应力混凝土结构当量荷载计算[J].建筑结构学报,1997,18(5):2211.
[87]赵勇,黄鼎业,李云贵,等.预应力混凝土板的等效荷载计算[J].建筑科学,2002,18(4):8-11.
[88]张义,张燎军.有限元等效应力基于ANSYS的计算[J].水利科技与经济,2006,12(4):266-267.
[89]王磊,刘寒冰,吴斌暄,等.六面体单元下预应力等代荷载计算[J].计算力学学报,2005,22(6):785-787.
[90]Frenet-Serret formulas[EB/OL].http://www.answers.com/topic/frenet-serret-formulas,2008.
[91]陈省身,陈维桓.微分几何讲义(第二版)[M].北京:北京大学出版社,2003-1.
[92]王勋成.有限单元法[M].北京:清华大学出版社,2003-7-1.
[93]Newton-Iteration Method[EB/OL].http://mathworld.wolfram.com/NewtonsIteration.html,2008.
[94]Jacobian[EB/OL].http://mathworld.wolfram.com/Jacobian.html,2008.
[95]“数学手册”编写组.数学手册[M].北京:人民教育出版社,1979.
[96]赵先琼,杨晓红.ANSYS有限元分析与20节点块单元[J].岳阳师范学院学报(自然科学版),2001,14(3):29-31.
[97]Guide to Interface with ANSYS[M].ANSYS Inc.1998.
[98]张义,张燎军.有限元等效应力基于ANSYS的计算[J].水利科技与经济,2006,12(4):266-267.
[99]Filipiak M.Mesh Generation[M].Edinburgh:Edinburgh Parallel Computing Centre,1996:25,31-35.
[100]郭仁忠.空间分析[M].高等教育出版社.2001.
[101]Thompson J,Wasi.Z U A,Mastin C W.Numerical Grid Generation[M].Elsevier Science Publishing Co.inc.1997.
[102]熊英,胡亍进,赵建军.基于映射法和Delaunay方法的曲面三角网格划分算法[J].计算机辅助设计与图形学学报,2002,14(1):56-60.
[103]周培德.计算几何算法分析与设计[M].北京:清华大学出版社,2000:38-130.
[104]骆冠勇,曹洪.一种网格和节点同步生成的二维Delaunay网格划分算法[J].计算机辅助设计与图形学学报,2007,19(5):604-608,615.
[105]Cignoni P,Montani C,Scopigno R.DeWall:A fast divide and conquer Delaunay triangulation algorithm in Ed[J].Computer-Aided Design,1998,30(5):333-341.
[106]Lattuada R,Raper J.Applications of 3D Delaunay triangulation algorithms in geoscientific modeling[EB/OL],http://geog.bbk.ac.uk,1997.
[107]Jonathan R S.Tetrahedral Mesh Generation by Delaunay Refinement[EB/OL],http://www.cs.cmu.edu,1998.
[108]周素琴.旋转矩阵及其性质[J].上海师范大学学报(自然科学版),2001,30(1):89-91.
[109]Wolfram M.Rotation matrix[EB/OL].http://mathworld.wolfram.com/RotationMatrix.html,2003-8-6.
[110]Wikipedia.Rotation matrix[EB/OL].http://en.wikipedia.org/wiki/Rotation_matrix.html,2008-5-22.
[111]郭思旭.高等数学(第四版)-上下册[M].北京:高等教育出版社.1997,7
[112]李娜,周小勇.任意截面几何特性研究[J].山西建筑,2007,33(21):58-60.
[113]Kianoosh H,Richard J B.Development and verification of a numerical model for the analysis of geosynthetic-reinforced soil segmental walls under working stress conditions[J].Can.Geotech.J.2005,42(4):1067-1085.
[114]Irons B M.Quadrature Ruls for Brick Based on Finite Elements[J].Num.Meth Eng.,1971(3):293-294.
[115]王磊,刘寒冰,王延华,等.用曲线积分计算桥梁正截面承载能力[J].哈尔滨工业大学学报.2004,36(10):1314-1317.
[116]王磊.钢筋混凝土构件正截面内力改进条带法计算[J].建筑科学.2004,20(2):5-8
[117]李律,钱济章,李敦.基于ANSYS的预应力筋数值模拟[J].公路工程,2007,32(4):178-179.
[118]陆新征.ANSYS软件简介及工程应用[M].清华大学土木工程系,2004-3.
[119]蒲俊,吉家锋.Matlab6.0数学手册[M].上海:浦东电子出版社,2002.
[120]Guide to Interface with ANSYS[M].ANSYS Inc.1998.
[121]The APDL Prograrmner's Guide[M].ANSYS Inc.1998.
[122]The UIDL Programmer's Guide[M].ANSYS Inc.1998.
[123]段志东,张克华.ANSYS图形用户界面二次开发[J].兰州铁道学院学报(自然科学版),2002,21(1):44-46.
[124]李新坡,袁文忠.ANSYS结构计算二次开发技术研究及应用[J].岩土工程技术,2004,18(2):100-102.
[125]刘维.精通Matlab与C/C++混合程序设计[M].北京:北京航天航空大学出版社,2005:145-212.
[126]谢刚,王小林.MATLAB与ANSYS数据接口的开发研究及应用[J].机械工程与自动化,2005,129(2):52-54.
[127]姜长元.VC++与Matlab混合编程的方法及实现[J].电脑知识与技术;2004,14(4):21-24
[128]李天昀,葛临东.综述MATLAB与VC++的交互编程[J].计算机仿真,2004,21(9):193-196.
[129]使用Matlab Eengine实现和Matlab的进程间通信[J].现代计算机,2002,85(2):31-35.
[130]MIDEVA,MATCOM & Visual MATCOM User's Guide[M].MathTools Ltd.1999-6.
[2]毕桂平,魏红一,范立础.ANSYS在桥梁工程中的应用前景[A].中国土木工程学会桥梁及结构工程学会第十四届年会论文集[C].上海:同济大学出版社,2000:410-412.
[3]江见鲸,贺小岗.工程结构计算机仿真分析[M].北京:清华大学出版社,1995.
[4]周小勇,金文成.洛河大桥预应力连续梁V型0#节段结构构造及受力特性分析[J].公路交通科技(应用技术版),2007(4):145-149.
[5]顾祥林,孙飞.混凝土结构的计算机仿真[M].上海:同济大学出版社,2002.
[6]程纬.基于ANSYS二次开发的钢筋混凝土结构倒塌仿真分析及其程序设计[博士后报告].上海:同济大学,2002.
[7]喻定新.以ANSYS为分析器的结构优化软件研制—杆板组合结构优化[硕士学位论文].南宁:广西大学,2004.
[8]段志东.ANSYS二次开发途径研究[硕士学位论文].兰州:兰州铁道学院,2002.
[9]魏志刚,汤文成,刘德仿.有限元分析在结构设计中的应用[J].设计与研究,2005,32(8):12-14.
[10]李正,徐立功.ANSYS二次开发在桥梁设计中的应用[J].江苏交通工程,2003,115(6):22-25.
[11]刘艳萍,杨新华,杨文兵.预应力钢筋混凝土局部有限元分析的ANSYS二次开发[J].华中科技大学学报(城市科学版),2005,22(增刊):87-90.
[12]杨建明.预应力混凝土结构的当量荷载计算[J].建筑结构学报,1997,18(5):2-11.
[13]ANSYS theory reference Release 5.7[M].ANSYS Inc.1994.
[14]Analysis of Geotechnical Problems with ABAQUS[M].ABAQUS inc.2001.
[15]ADINA Theory and Modeling Guide[M].ADINA inc.2002.
[16]About ANSYS/CivilFEM[EB/OL].http://www.civilfem.com/about_us/about_us.php,2008.
[17]桥梁三维预应力空间分析系统软件帮助文件(bridgeKF1.0)[M].中华人民共和国交通部科研所,2007-2.
[18]Guide to ANSYS user programmable features[M].ANSYS Inc.1998.
[19]Duane D,Bruce L.Mastering Matlab 7[M].London:Prentice Hall,2004-10.
[20]Compaq Visual FORTRAN Programmer's Guide[M].Houston:Compaq Computer Corporation,2001-8.
[21]Beck Z.Microsoft Visual C++ 6.0 Programmer's Guide[M].北京:北京希望电子出版社,1998-12.
[22]MSC.Nastran简介[EB/OL].http://www.baisi.net/viewthread.php,2008.
[23]John O H.LS_DYNA Theoretical Manual[M].Livermore Software Cor.1997.
[24]CGAL User and Reference Manual Table of Contents[EB/OL].http://www.cgal.org/Manual/3.3/doc_html/cgal_manual/contents.html.2007.
[25]Spatial公司主页(中文版)[EB/OL].http://www.spatial.com.cn/community/case.html,2008.
[26]史忠植,罗杰文,王茂光,等.互联网智能平台AGrIP构建及其应用[A].信息技术快报,2005,3(5):24-30.
[27]王永岩.有限元法基础与SAP5程序应用[M].北京:煤炭工业出版社,1995.
[28]钟万勰.有限元分析与优化设计软件JIFEX研究开发和应用[EB/OL].http://sail.dlut.edu.cn/chengguo/04.htm,2003.
[29]通用有限元软件Strat用户手册[M].上海交通大学土木系,上海佳构软件科技有限公司,2003-5.
[30]施法中.计算机辅助几何设计与非均匀有理B样条[M].北京:高等教育出版社,2001.
[31]Boor C D.On calculating with B-Splines[J].J Approx.Theory,1972,6(1):50-62.
[32]金文成,周小勇.20节点等参单元的预应力等效节点力计算[J].华中科技大学(城市科学版),2008,25(1):5-7.
[33]张立明.Algor＼Ansys在桥梁工程中的应用方法与实例[M].北京:人民交通出版社,2003.
[34]吕勇刚.CAGD自由曲线曲面造型中均匀样条的研究[博士学位论文][D].杭州:浙江大学,2002.
[35]Sablonniere P.Spline and Bezier polygons associated with a polynomial spline curve[J].Computer Aided Design,1978,10(4):257-261.
[36]Barnhill R E,Riesenfield R F.Computer Aided Geometric Design[M],New York:Academic Press,1974.
[37]Cox M G.The Numerical elvaluation of B-Splines[J].J.Inst.Math.Applications,1972.
[38]Yamaguchi F,Curves and surfaces in computer aided geometric design[M].Berlin:Springer-Verlag,1988.
[39]Ferguson J.Multivariable curve interpolation[J].ACM,1964,11(2):221-228.
[40]Coons S.Rational bi-cubic surface patches[R].Boston:Massachusetts Institute of Technology Project MAC,1968.
[41]Schoenberg I.Spline functions and the problem of graduation[J].Proceedings of the National Academy of Sciences of USA,1964(52):947-950.
[42]Bezier P E.Numerical control-Mathematics and applications[M].London:Forrest Tran,1972:16-53.
[43]Boor D C.On Calculation with B-Spline[J].Approx Theory,1972,6(1):50-62.
[44]Boehm W.Inserting new knots into B-Spline curves.Computer-Aided Design,1980,12(4):199-201.
[45]Forrest A R.Interactive interpolation and approximation by Bezier polynomials[J].The Computer,1972,15(1):71-79.Reprinted in CAD 1990,22(9):527-537.
[46]ISO 10303,Standard for the Exchange of Product Model Data(STEP)[S].
[47]BS EN 29592-2:1991,Computer graphics.Programmer's Hierarchical Interactive Graphics System(PHIGS)[S].
[48]李翠华,郑南宁,张永平.基于样条修匀公式的图像边缘检测[J].电子学报,1999,27(1):2-5.
[49]Prautzsch H.Degree elevation of B-Spline curves[J].Computer Aided Geometric Design,1(12):193-198,1984.
[50]Peng Q.An algorithm for finding the intersection lines between two B-Spline surfaces[J].Computer Aided Design,1984,16(4):191-196.
[51]戴公连,李德建.桥梁结构空间分析设计方法与应用[M].北京:人民交通出版社,2001.10
[52]Riesenfeld R F.Nonuniform B-Spline curves[C].2nd USA-JAPAN Confence proc.1975:551-555.
[53]Hartley D J,Judd C J.Parametrization and shape of B-Spline curves[J].Computer Aided Design,1980,12(5):235-238.
[54]Prautzseh H,Piper B.A fast algorithm to raise the degree of B-Spline curves[J].Computer Aided Geometric Design,1991,8(4):253-266.
[55]李翠华,郑南宁,张永平.基于样条修匀公式的图像边缘检测[J].电子学报,1999,27(1):2-5.
[56]JTG D62-2004,公路钢筋混凝土及预应力混凝土桥涵设计规范[S].北京:人民交通出版社,2004.
[57]李铁军,杨红梅.计算机辅助几何造型设计中曲面拼接的处理[J].辽宁工学院学报,1999,19(4):23-26.
[58]Farin G.Curvature continuity and offsets for piecewise conics,ACM Transactions on Graphics,1989,8(2):89
[59]李卫国,徐步青,胡克辉.空间多曲线钢索预应力摩擦损失试验研究[J].2000,13(1):14-17.
[60]南开大学数学系.空间解析几何引论[M].北京:人民教育出版社,1978-9.
[61]Boor C D.A Practical Guide to Spline[M].Berlin:Springer-Verlag,1978.
[62]李红.数值分析[M].武汉:华中科技大学出版社,2003-11.
[63]李厚源.空间解析几何[M].济南:山东科学技术出版社1983-11.
[64]MathWorks Inc.MATLAB C++ Math Library 2.0 Reference[EB/OL].http://www.thphys.may.ie/CompPhysics/matlab/help/extern/cppmath/cpp_cref.html,2008.
[65]Prautzsch H.Degree elevation of B-Spline curves[J].Computer Aided Geometric Design,1984,12(1):193-198.
[66]施法中.论参数曲线曲面的几何不变性[J].工程图学学报,1991(1):79-86.
[67]王天军.曲面造型中几何连续曲面的研究与应用[博士学位论文][D].北京:北京航天航空大学,1991.
[68]Sapidis N,Farin G.Automatic fairing algorithm for B-Spline curves[J].Computer Aided Design,1990,22(2):121-129.
[69]Stam J.On subdivision schemes generalizing uniform B-Spline surfaces of arbitrary degree[J].Computer Aided Geometric Design,2001,18(5):383-396.
[70]李律,钱济章,李敦.基于ANSYS的预应力筋数值模拟[J].公路工程,2007,32(4):178-179.
[71]金文成,张小飞,周小勇.斜拉桥塔索锚固区空间应力分析[J].公路交通科技,2007,24(9):69-73.
[72]MIDAS/Civil V6.7.1在线帮助手册[Z].北京:MIDAS Information Technology China Co.,Ltd,2005.
[73]周宗泽.桥梁博士系统V3.1版[M].同济大学桥梁工程系,2005-8.
[74]李国平.预应力混凝土结构设计原理[M].北京:人民交通出版社,2000:30-43.
[75]宋玉普,车轶,马德有,等.空间多曲线型预应力钢索的预应力摩擦损失研究[J].土木工程学报,2002,35(6):105-108.
[76]刘庆文,苏小卒,空间曲线型预应力摩擦损失计算[J].四川建筑科学研究,2007,33(1):42-44.
[77]胡狄,陈政清.考虑反向摩阻的后张法PC构件锚固损失的计算[J].中国公路学报,2004,17(1):34-38.
[78]Hung T.Anchorage take-up loss in post- tensioned members[J].Journal of the Prestressed Concrete Institute,1969,14(4):30-35.
[79]程万年.圆弧预应力筋锚固损失基本方程和解[J].力学与实践,1992,14(3):33-35.
[80]胡光祥,陈东杰,赵勇,等.空间曲线预应力束摩擦损失测试和摩擦系数反演[J].建筑科学,2004,20(5):26-29.
[81]叶建.下白石特大桥空间曲线孔道摩擦损失分析及测试[J].福建建筑,2005,93(3):115-116.
[82]刘志文,宋一凡,赵小星,等.空间曲线预应力束摩阻损失参数[J].西安公路交通大学学报,2001,21(3):42-44.
[83]李卫国,徐步青,胡克辉.空间多曲线钢索预应力摩擦损失试验研究[J].石家庄铁道学院学报,2000,13(1):14-17.
[84]周小勇,金文成.独塔单索面斜拉桥-施州大桥设计和施工[J].中国市政工程,2006,124(6):23-25,88.
[85]刘庆文,苏小卒.后张法预应力空间多曲线钢筋锚固损失简化计算[J].铁道建筑,2006(2):10-12.
[86]杨建明.预应力混凝土结构当量荷载计算[J].建筑结构学报,1997,18(5):2211.
[87]赵勇,黄鼎业,李云贵,等.预应力混凝土板的等效荷载计算[J].建筑科学,2002,18(4):8-11.
[88]张义,张燎军.有限元等效应力基于ANSYS的计算[J].水利科技与经济,2006,12(4):266-267.
[89]王磊,刘寒冰,吴斌暄,等.六面体单元下预应力等代荷载计算[J].计算力学学报,2005,22(6):785-787.
[90]Frenet-Serret formulas[EB/OL].http://www.answers.com/topic/frenet-serret-formulas,2008.
[91]陈省身,陈维桓.微分几何讲义(第二版)[M].北京:北京大学出版社,2003-1.
[92]王勋成.有限单元法[M].北京:清华大学出版社,2003-7-1.
[93]Newton-Iteration Method[EB/OL].http://mathworld.wolfram.com/NewtonsIteration.html,2008.
[94]Jacobian[EB/OL].http://mathworld.wolfram.com/Jacobian.html,2008.
[95]“数学手册”编写组.数学手册[M].北京:人民教育出版社,1979.
[96]赵先琼,杨晓红.ANSYS有限元分析与20节点块单元[J].岳阳师范学院学报(自然科学版),2001,14(3):29-31.
[97]Guide to Interface with ANSYS[M].ANSYS Inc.1998.
[98]张义,张燎军.有限元等效应力基于ANSYS的计算[J].水利科技与经济,2006,12(4):266-267.
[99]Filipiak M.Mesh Generation[M].Edinburgh:Edinburgh Parallel Computing Centre,1996:25,31-35.
[100]郭仁忠.空间分析[M].高等教育出版社.2001.
[101]Thompson J,Wasi.Z U A,Mastin C W.Numerical Grid Generation[M].Elsevier Science Publishing Co.inc.1997.
[102]熊英,胡亍进,赵建军.基于映射法和Delaunay方法的曲面三角网格划分算法[J].计算机辅助设计与图形学学报,2002,14(1):56-60.
[103]周培德.计算几何算法分析与设计[M].北京:清华大学出版社,2000:38-130.
[104]骆冠勇,曹洪.一种网格和节点同步生成的二维Delaunay网格划分算法[J].计算机辅助设计与图形学学报,2007,19(5):604-608,615.
[105]Cignoni P,Montani C,Scopigno R.DeWall:A fast divide and conquer Delaunay triangulation algorithm in Ed[J].Computer-Aided Design,1998,30(5):333-341.
[106]Lattuada R,Raper J.Applications of 3D Delaunay triangulation algorithms in geoscientific modeling[EB/OL],http://geog.bbk.ac.uk,1997.
[107]Jonathan R S.Tetrahedral Mesh Generation by Delaunay Refinement[EB/OL],http://www.cs.cmu.edu,1998.
[108]周素琴.旋转矩阵及其性质[J].上海师范大学学报(自然科学版),2001,30(1):89-91.
[109]Wolfram M.Rotation matrix[EB/OL].http://mathworld.wolfram.com/RotationMatrix.html,2003-8-6.
[110]Wikipedia.Rotation matrix[EB/OL].http://en.wikipedia.org/wiki/Rotation_matrix.html,2008-5-22.
[111]郭思旭.高等数学(第四版)-上下册[M].北京:高等教育出版社.1997,7
[112]李娜,周小勇.任意截面几何特性研究[J].山西建筑,2007,33(21):58-60.
[113]Kianoosh H,Richard J B.Development and verification of a numerical model for the analysis of geosynthetic-reinforced soil segmental walls under working stress conditions[J].Can.Geotech.J.2005,42(4):1067-1085.
[114]Irons B M.Quadrature Ruls for Brick Based on Finite Elements[J].Num.Meth Eng.,1971(3):293-294.
[115]王磊,刘寒冰,王延华,等.用曲线积分计算桥梁正截面承载能力[J].哈尔滨工业大学学报.2004,36(10):1314-1317.
[116]王磊.钢筋混凝土构件正截面内力改进条带法计算[J].建筑科学.2004,20(2):5-8
[117]李律,钱济章,李敦.基于ANSYS的预应力筋数值模拟[J].公路工程,2007,32(4):178-179.
[118]陆新征.ANSYS软件简介及工程应用[M].清华大学土木工程系,2004-3.
[119]蒲俊,吉家锋.Matlab6.0数学手册[M].上海:浦东电子出版社,2002.
[120]Guide to Interface with ANSYS[M].ANSYS Inc.1998.
[121]The APDL Prograrmner's Guide[M].ANSYS Inc.1998.
[122]The UIDL Programmer's Guide[M].ANSYS Inc.1998.
[123]段志东,张克华.ANSYS图形用户界面二次开发[J].兰州铁道学院学报(自然科学版),2002,21(1):44-46.
[124]李新坡,袁文忠.ANSYS结构计算二次开发技术研究及应用[J].岩土工程技术,2004,18(2):100-102.
[125]刘维.精通Matlab与C/C++混合程序设计[M].北京:北京航天航空大学出版社,2005:145-212.
[126]谢刚,王小林.MATLAB与ANSYS数据接口的开发研究及应用[J].机械工程与自动化,2005,129(2):52-54.
[127]姜长元.VC++与Matlab混合编程的方法及实现[J].电脑知识与技术;2004,14(4):21-24
[128]李天昀,葛临东.综述MATLAB与VC++的交互编程[J].计算机仿真,2004,21(9):193-196.
[129]使用Matlab Eengine实现和Matlab的进程间通信[J].现代计算机,2002,85(2):31-35.
[130]MIDEVA,MATCOM & Visual MATCOM User's Guide[M].MathTools Ltd.1999-6.