核电站反应堆压力容器冷热态密封分析
|
摘要
本论文来源于国防科学技术工业委员会关于先进压水堆核电站先进反应堆设计技术研究项目（批文号：科工委[1999]149号）的子课题“反应堆压力容器瞬态密封程序研制”。本课题研究的目的是为核电站国产化及自主设计提供技术准备，提高设计质量和设计水平，确保先进压水堆核容器密封系统设计的可靠性，确保核安全，使核容器密封设计软件达到世界先进水平。
本文开发了反应堆压力容器瞬态密封有限元分析程序FEARPV的前处理程序，实现了参数化的有限元分析模型的自动生成，建立较接近反应堆压力容器密封结构实际情况的计算模型。
应用反应堆压力容器瞬态密封有限元分析程序对反应堆压力容器模拟体密封结构进行了螺栓预紧、升温升压、恒温恒压、降温降压等全过程分析。所得计算结果与中国核动力运行研究所（武汉）和中国核动力院设计研究所在压力容器模拟体上所测的大量应力、应变、温度、螺栓力、分离量等数据进行了对比，在两种载荷工况下，无论是冷态还是热态，无论是定性规律还是定量数据，分析结果和实验结果都吻合良好。
应用反应堆压力容器瞬态密封有限元分析程序对电站核压力容器密封结构进行了螺栓预紧、升温升压、恒温恒压、降温降压等全过程的三维瞬态密封分析。由本文分析计算可知，虽然电站核压力容器的结构及计算工况较压力容器模拟体复杂得多，然而反应堆压力容器瞬态密封有限元分析程序对塑性迭代和接触迭代仍具有很好的收敛性，数值稳定性好，得到的结果具有良好的规律性，在升温、升压、恒温、恒压、降温和降压等过程中其规律和压力容器模拟体基本相似。计算结果表明，由于密封区斜面的存在，使密封面的接触区增大，接触压力分布更趋合理，并有助于减小内外槽的分离量。
The theme comes from the project The Program Development for Transient Sealing in Reactor Pressure Vessel. This project is the subproject of the research project commissioned by the Technical Industry Committee, Department of Defense, which is about the advanced reactor design in advanced pressurized water reactor nuclear power station. The aim of this research project is to supply the technical preparation for the nuclear power station domestalization and autonomous design, to improve the design quality and level, to ensure the reliability of the sealing system in the advanced pressurized water reactor vessel, to ensure the nuclear safety, and to raise the design software level of the nuclear vessel sealing system to the advanced world class.
In this thesis, the pre-processing program of FEARPV (Finite Element Analysis for Sealing System in Reactor Pressure Vessel) is developed, parameterization finite element automatical modeling is achieved, and computational model which is close to actual sealing system structure in the reactor pressure vessel is built.
The program FEARPV is applied to analysis the sealing structure of the pressure vessel model during the total process, including bolt pretensioning process, heating and compression process, thermostated heating and constant pressure process, and the cooling and decompressing process. Then the computational results are compared with the experimental data including stress, strain, temperature, bolt force and separation, which are obtained by the China Nuclear Power Operation Research Institute (Wuhan) and China Nuclear Power Academy Design Research Institute. Either under cooling or heating cycle, no matter the qualitative regularity or quantitative data, the computational results are in common with the experimental data.
The program FEARPV is applied to 3-D Transient sealing analysis of nuclear pressure vessel during the total process, including bolt pretensioning process, heating and compression process, thermostated heating and constant pressure process, and the cooling and decompressing process. According to the analysis and computational results in this thesis, although the actual structure and work conditions of the pressure vessel in nuclear power station are much more complex than that of the pressure vessel model, the program FEARPV still can get quite good convergence on the plasticity and
contact iteration. And also the program FEARPV has good value stability, and good regularity. During the heating, compression, constant temperature, constant pressure, cooling and decompression process, the regularity is basically in common with that of the pressure vessel model. According to the computational results, due to the existence of sealing zone slope, the contact zone of sealing zone is increased, the distribution of contact stress become more reasonable, and the separation of inner and outer canal is decreased.
本文开发了反应堆压力容器瞬态密封有限元分析程序FEARPV的前处理程序，实现了参数化的有限元分析模型的自动生成，建立较接近反应堆压力容器密封结构实际情况的计算模型。
应用反应堆压力容器瞬态密封有限元分析程序对反应堆压力容器模拟体密封结构进行了螺栓预紧、升温升压、恒温恒压、降温降压等全过程分析。所得计算结果与中国核动力运行研究所（武汉）和中国核动力院设计研究所在压力容器模拟体上所测的大量应力、应变、温度、螺栓力、分离量等数据进行了对比，在两种载荷工况下，无论是冷态还是热态，无论是定性规律还是定量数据，分析结果和实验结果都吻合良好。
应用反应堆压力容器瞬态密封有限元分析程序对电站核压力容器密封结构进行了螺栓预紧、升温升压、恒温恒压、降温降压等全过程的三维瞬态密封分析。由本文分析计算可知，虽然电站核压力容器的结构及计算工况较压力容器模拟体复杂得多，然而反应堆压力容器瞬态密封有限元分析程序对塑性迭代和接触迭代仍具有很好的收敛性，数值稳定性好，得到的结果具有良好的规律性，在升温、升压、恒温、恒压、降温和降压等过程中其规律和压力容器模拟体基本相似。计算结果表明，由于密封区斜面的存在，使密封面的接触区增大，接触压力分布更趋合理，并有助于减小内外槽的分离量。
The theme comes from the project The Program Development for Transient Sealing in Reactor Pressure Vessel. This project is the subproject of the research project commissioned by the Technical Industry Committee, Department of Defense, which is about the advanced reactor design in advanced pressurized water reactor nuclear power station. The aim of this research project is to supply the technical preparation for the nuclear power station domestalization and autonomous design, to improve the design quality and level, to ensure the reliability of the sealing system in the advanced pressurized water reactor vessel, to ensure the nuclear safety, and to raise the design software level of the nuclear vessel sealing system to the advanced world class.
In this thesis, the pre-processing program of FEARPV (Finite Element Analysis for Sealing System in Reactor Pressure Vessel) is developed, parameterization finite element automatical modeling is achieved, and computational model which is close to actual sealing system structure in the reactor pressure vessel is built.
The program FEARPV is applied to analysis the sealing structure of the pressure vessel model during the total process, including bolt pretensioning process, heating and compression process, thermostated heating and constant pressure process, and the cooling and decompressing process. Then the computational results are compared with the experimental data including stress, strain, temperature, bolt force and separation, which are obtained by the China Nuclear Power Operation Research Institute (Wuhan) and China Nuclear Power Academy Design Research Institute. Either under cooling or heating cycle, no matter the qualitative regularity or quantitative data, the computational results are in common with the experimental data.
The program FEARPV is applied to 3-D Transient sealing analysis of nuclear pressure vessel during the total process, including bolt pretensioning process, heating and compression process, thermostated heating and constant pressure process, and the cooling and decompressing process. According to the analysis and computational results in this thesis, although the actual structure and work conditions of the pressure vessel in nuclear power station are much more complex than that of the pressure vessel model, the program FEARPV still can get quite good convergence on the plasticity and
contact iteration. And also the program FEARPV has good value stability, and good regularity. During the heating, compression, constant temperature, constant pressure, cooling and decompression process, the regularity is basically in common with that of the pressure vessel model. According to the computational results, due to the existence of sealing zone slope, the contact zone of sealing zone is increased, the distribution of contact stress become more reasonable, and the separation of inner and outer canal is decreased.
引文
[1] 许连义. 中国核电发展的新篇章. 动力工程. 1997，17(5)：7~10，6
[2] 陈广. 上海核电国产化优势初步分析. 工业技术进步. 1997，(3)：13~15
[3] 吴宗鑫，张作义. 世界核电发展趋势和高温气冷堆. 核科学与工程. 2000，20(3)：211~219，231
[4] 压水堆压力容器法兰密封技术专题报告．中国核动力研究设计情报室，1990
[5] D H Van Campen, P J Deen, D G H Latzko. Deformation of Large-diameter, High-pressure Vessel Flanges. Proc. 1st Int. Conf. on Pressure Vessel Tech., Delft, 1969
[6] C Menken. Influence of Bolt Loading on Deformation of Pressure Vessel Flanges. Proc. 1st Int. Conf. on Pressure Vessel Tech., Delft, 1969
[7] N W Murray, D G Stuart. Proc. Symp. on Pressure Vessel Research towards Better Design. Int. Mech. Eng., London, 1970
[8] D H Van Campen. Some Applications of Engineering Mechanics in Pressure Vessel Analysis Thesis. Delft, 1972
[9] Z Zudans, T C Yen, W H Steigelmann. Thermal stress Techniques in the Nuclear Industry. New York, American Elsevier Publishing Company, 1975
[10] A bakker, H A C M Spass. A Computational Procedure for the Stress and Deformation Analysis of High-pressure Vessel Flanges. Int. J. Pres. Ves. & Piping. 1977, (5)：135~149
[11] H A C M Spass. Deformation Behavior of Large High-pressure Vessel Flanges Development in Stress Analysis for Pressurised Components. R W Nichols (ed.). Elsevier Applied Science Publishers, London, 1977
[12] H Sacher, M Mayk. Structure Analysis of a Reactor Pressure Vessel by Two- and Three- Dimensional Models. Int. J. Pres. Ves. & Piping. 1982, (10)：147~150
[13] D H Cascales, C Militello. An Accurate Simple Model to Evaluate the Integral Flange Rotation. Int. J. Pres. Ves. & Piping. 1987, (30)：54~65
[14] D H Cascales, C Militello, W J Mulhall. Bolted Flanged Joint Leakage Analysis with Non-linear Gasket Behavior. Int. J. Pres. Ves. & Piping. 1987, (30)：78~92
[15] Ceorge F. P.E.Leon. An overview of the US PVRC Research Program on Bolted Flanged Connections, Pressure Vessel Technology Vol. 1 Design & Analysis, Pergramon Press 1988：32~65
[16] Gu Fangyu, Zeng Fanmao. Pretensioning of Pressure Vessel Bolts as a Problem Including
Elastic Contact with Adjustable Gaps. Int. J. Pres. Ves. & Piping. 1985, (19)：12~25
[17] 谷芳毓，曾凡茂，程中易．可调间隙弹性接触问题．核动力工程. 1982, (10)：25~29
[18] 曲家棣．核容器瞬态热分析、冷态密封分析程序系统．上海核工程研究设计院．1985
[19] 李润方，陈如欣，李水安．反应堆压力容器法兰密封系统密封分析程序系统理论文本．重庆大学，中国核动力研究设计院．1991
[20] 李润方，林腾蛟．反应堆压力容器局部三维冷热态密封分析．机械工程学报，1993，29(6):104-108．
[21] 李润方，李水安，李健，陈骏，林腾蛟. 核容器法兰密封系统密封性能数值分析. 核动力工程. 1994，15(1)：12-18．
[22] 曲家棣，盛贤科，王佩珠. GH169度银O形环冷、热态性能测试. 压力容器. 1997，14(3)：183~188
[23] 曲家棣，王佩珠，窦一康，盛贤科，谢世球，金惠琴. 压力容器热态密封实验与分析. 1997，14(4)：1~7
[24] 彭向和，高芝晖，曾祥国. O形密封环的弹塑性大变形和循环松弛分析. 重庆大学学报. 1998，21(2)：46~50
[25] 于溯源，何树延，刘俊杰. HR-200压力容器主法兰密封的分析法设计. 第十一届全国反应堆结构力学会议，成都，2000, 10
[26] 左卫东，于溯源，刘俊杰等. 10MW高温气冷堆压力容器主法兰结构的有限元接触分析. 核动力工程. 2001，22(3):226-231.
[27] R W Clough. The Finite Element Method in Plane Stress Analysis. Proceedings of the Second Conference on Electronic Computation. American Society of Civil Engineers, New York, 1960
[28] C H Lee, S Kobayashi. Elastoplastic Analysis of Plane stain and Axisymmetric Flat Punch Indentation by the Finite Element Method. Int. J. Mech. Sci.. 1970, (12)：25~38
[29] J T Oden, D R Bhandari, O Yagawa, T J Chang. A New Approach to the Finite Element Formulation and Solution of a Class of Problem in Coupled Thermoelasto-viscoplasticity of Crystalline Solids. Nuclear Engineering Design. 1973, (24)：135~152
[30] D R J Owen, E M Salonen. Three Dimensional Elasto-plastic Finite Element Analysis，Int. J. Num. Mech. Engng. 1975, (19)：248~265
[31] Hindon, Owen. Finite Element Programming. Academic Press, 1977：56~78
[32] Owen, Hinton. Finite Element in Plasticity. Pineridge Press Limited, 1980：102~150
[33] B M Iron. Engineering Applications of Numerical Integration in Stiffness Methods. AIAAJ. 1966, 4(11)：378~392
[34] B M Irons, A Frontal Solution Program. Int. J. Num. Mech. Engng. 1970, 2：456~473
[35] E A Wilson, B Parsons. Finite Element Analysis of Elastic Contact Problems using Differential Displacements. Int. J. Num. Mech. Engng. 1970, 2：248~256
[36] S K Chan, I S Tuba. A Finite Element Method for Contact Problems of Solid Bodies, II, Application to Turbine Blade Fastening. Int. J. Mech. Sci. 1971, 13：132~145
[37] S Ohte. Finite Element Analysis of Elastic Contact Problems. Bull. JSME 16, 1973
[38] A Francavilla, O C Zienkiewicz. A Note on Numerical Computation of Elastic Contact Problem. Int. J. Num. Mech. Engng. 1975, 19：913-924
[39] B Fredriksson. Finite Element Solution of Surface Nonlinearities in Structural Mechanics with Special Emphasis to Contact and Fracture Mechanics Problem. Computer and structure. 1976, 6：514~530
[40] 冈本纪明．有限要素法たょゐ非线性接触问题の解析. 日本机械学会论文. 1977, 372~386
[41] 茑纪夫．有限要素法の接触问题への应用．机械の研究. 1978, 2(8)：28~43
[42] J T Stadter, R O Weiss. Analysis of Contact through Finite Element Gaps. Computer and structure. 1979, 10(6)： 867-873
[43] 陈万吉．用有限元混合法分析弹性接触问题．大连工学院学报．1979, (2)：12~15
[44] 陈万吉．耦合热弹性接触问题有限元分析．固体力学学报．1981, (3)：17~21
[45] R F Li. Finite Element Analysis for Elastic-plastic Contact Problems with Friction. C/R/523/85, 1985
[46] 郭小明. 岩土介质弹塑性耦合问题的数学规划模型. 东南大学学报. 1994，24(5)：101~105
[47] 李锡夔. 蠕变-弹塑性-损伤耦合响应的有限元分析. 计算结构力学及其应用. 1995，12(1)：26~38
[48] 张承柱，刘信声. 应变空间表述的混凝土弹塑性耦合本构模型. 清华大学学报. 1996，36(3)：59~64
[49] 倪强，王乘，黄玉盈，张标标. 焊接过程的耦合三维热弹塑性有限元分析. 华中理工大学学报. 1997，25(1)：84~86
[50] 吕安军，曹志远. 地下结构和岩土介质弹塑性耦合分析的摄动半解析法. 应用数学和力学. 1998，19(4)：335~339
[51] 孙秦，徐晓飞. 弹塑性耦合张量的宏细观对偶形式. 应用力学学报. 1998，15(4)：19~21
[52] 张军，吴昌华. 轮轨接触问题的弹塑性分析. 铁道学报. 2000，22(3)：16~21
[53] 栾茂田，黎勇. 非连续变形计算力学模型弹塑性数值分析方法及其应用. 岩土工程学报. 2001，23(5)：530~535
[54] 高华林，蔡嗣经，童光煦. 有限单元法微机前处理程序的开发. 矿业研究与开发. 1994，14(2)：58~60
[55] 胡于进，王坚，陆金桂，周济. 任意形状平面域不完全Delaunay三角化研究. 华中理工大学学报. 1994，22(10)：4~8
[56] 陈大宏，王革，邹家祥. CAD技术在有限元分析前处理中的应用. 1994，16(增刊)：70~74
[57] 张广健，黄由玲. 用CAGD技术进行有限元前处理. 工程力学. 1995，12(4)：124~130
[58] 杜群贵. 三维实体有限元网格自动 Delaunay剖分. 华南理工大学学报(自然科学版). 1996，24(9)：46~50
[59] 安晓卫，郝长忠. 旋转结构的有限元网格自动生成. 机械设计. 1996，(12)：40~42
[60] 安晓卫，张华. 用于不规则几何形体的有限元网格自动生成处理技术. 沈阳工业学院学报. 1997，16(4)：82~85
[61] 赵铁成，沈月芬，朱国桢，许先壁. 锅炉锅筒下降管三通的三维有限元自动剖分. 锅炉技术. 1997，(6)：10~14
[62] 林腾蛟，李润方，陈兵奎，林超. 混合标架下宽斜齿轮三维接触有限元前处理技术. 机械传动. 1997，21(2)：9~12
[63] 吴淑芳，李占国. 二维有模成形自适应数值模拟系统. 塑性工程学报. 1998，5(3)：22~26
[64] 雷震宇，黄洪钟. 160t铁路救援起重机转台有限元建模专家系统数据库的研究. 建筑机械. 1998，(5)：21~24
[65] 饶志伟，赵洪伦. 车体建模CAD系统的设计. 计算机辅助设计与制造. 1998，(2)：28~29
[66] 王琳，周理兵，马志云. 一种适于大型水轮发电机磁场计算前处理的新方法. 1998，18(1)：39~41
[67] 林腾蛟，李润方，陈兵奎，李亦彦. 鼓形齿联轴器轮齿三维有限元模型自动生成及强度计算. 1998，11(3)：18~19
[68] 左正兴. 一种简单高效的有限元网格优化算法及应用. 机械设计. 1998，(8)：21~24
[69] 查德根，沈景凤，黄宗益. 焊接结构件有限元分析中的前处理技术. 上海水产大学学报. 1998，7(4)：340~344
[70] 虞春，周雄辉，张永清. 特征造型与有限元分析的集成化研究. 计算机辅助设计与图形学学报. 1999，11(1)：62~65
[71] 安晓卫，张华，钟震宇. 三通管件的有限元网格的自动剖分. 沈阳工业学院学报. 1999，18(3)：5~10
[72] 段红杰，阎德元. 复合材料整体异形构件的CAD技术研究. 北京航空航天大学学报. 1999，25(5)：585~587
[73] 王璋奇，安利强. 基于SuperSAP的有限元前处理接口程序设计. 华北电力大学学报. 1999，26(1)：74~77
[74] 范峰. 空间网格结构CAD软件ASSAP 1系统研制与开发. 哈尔滨建筑大学学报. 1999，32(5)：16~19
[75] 倪其民，钱瑞明，王水来，阮雪榆. 回转支承智能设计系统研究与开发. 机械科学与技术. 1999，18(5)：709~711
[76] 秦浦雄，王书庆. 桥梁结构的三维建模技术. 计算机辅助工程. 2000，(2)：46~52
[77] 谭也平. 有限元前处理系统的图形集成实现. 计算机辅助工程. 2000，(1)：53~58
[78] 潘子杰，杨文通，郭九生. 任意平面区域有限元三角形网格全自动剖分. 北京工业大学学报. 2001，27(4)：493~496
[79] 陈建梅，关湛湛. 电器磁场计算可视化前处理软件开发研究. 三峡大学学报(自然科学版). 2001，23(1)：35~37
[80] 安利强，孟宪国，王璋奇. 转子热应力有限元前后处理方法. 华北电力大学学报. 2001，28(1):68~71
[81] 王璋奇，安利强，彭震中. 基于参数化建模的转子有限元剖分. 热能动力工程. 2001，16(9)：538~539
[82] 马智慧，胡兆吉，李鸣. 板料成形中的载荷及位移边界条件的处理. 化工装备技术. 2002，23(2)：29~31
[83] 王璋奇，安利强，王孟. 有限元分析结果的可视化处理方法. 工程图学学报. 2002，(1)：72~77
[84] 李建国，刘实，李殿中，戎利建，范寸，李依依. 汽车纵梁冲压成形及回弹过程的计算机模拟. 塑性工程学报. 2001，8(2)：4~7
[85] W. E. Lorence, H. E. Cline. MarchingCube: A high resolution 3D surface construction algorithm. Computer Graphics(ACM Siggraph Proceedings). 1988, 22(4)：51~58
[86] 程世才，林国成，满波. 实用有限元后处理程序的研制及开发. 哈尔滨工业大学学报. 1996，28(6)：45~49
[87] 仲梁维. 三维有限元后处理程序设计. 微型电脑应用. 1997，1：62~64，67
[88] 费少梅. 有限元计算结果的可视化. 内燃机工程. 1997，17(4)：26~28
[89] 王乐天，谢永慧，安宁，孟庆集. 汽轮机叶片三维有限元分析后处理. 西安交通大学学报. 1997，31(12)：33~38
[90] 林鹏，卓家寿，裴安荣. 等值面绘制及其工程应用. 西安理工大学学报. 1998，14(4)：413~416
[91] 虞春，周雄辉，张永清. 基于特征的智能有限元后处理. 计算机辅助工程. 1999，(1)：19~24
[92] 杨晓松，顾云宪，李云鹏，关振群. 三维有限元模型的任意剖切及其等值线与彩色云图生成的方法. 中国图象图形学报. 1999，4(7)：574~578
[93] 刘理，刘土光，张标标，李天匀. 有限元分析后处理程序中的等值面生成算法. 华中理工大学学报. 1999，27(7)：80~82
[94] 闵作兴，陈国华，吴衍平. 发动机有限元可视化应用软件研究. 华中理工大学学报. 1999，
27(11)：57~58
[95] 孙黎明，李世其. 有限元后处理中的一种快速多等值面面绘制算法. 武汉纺织工学院学报. 1999，12(1)：69~73
[96] 杨晓松，原仓周，李云鹏，关振群，顾云宪. 基于scan-buffer三维不规则数据场的直接体绘制技术. 大连理工大学学报. 2000，40(2)：198~202
[97] 苏海东. 有限元后处理图形软件FEPS的开发. 长江科学院院报. 2000，17(3)：19~22
[98] 李昱明，周飞，郭新贵，彭颖红. 八节点六面体有限元网格自动消隐技术. 计算机工程与应用. 2000，(11)：89~91
[99] 罗婷，张颖江，陈雯. 利用OpenGL实现有限元后处理. 湖北工学院学报. 2000，15(4)：7~9
[100] 詹梅，刘郁丽，杨合. 三维有限元模拟等值线生成的新方法. 重型机械. 2000，(1)：36~38
[101] 章勤，李品，苗京. 一种等值线图绘制算法及其并行实现. 华中科技大学学报. 2001，29(9)：4~6
[102] 诸葛良，刘健，王明武. 有限元后处理的一种动态可视化方法. 机械设计与研究. 2001，17(1)：29~30
[103] 李昱明，周飞，郭新贵，彭颖红. 基于六面体八节点单元的三维有限元数值模拟表面等值线的生成算法. 锻压技术. 2001，(2)：10~13
[104] 杨晓松，李云鹏，顾云宪. 基于Slicing的快速有限元网格体绘制算法. 计算机辅助设计与图形学学报. 2001，13(4)：348~352
[105] 李润方，龚剑霞．接触问题数值方法及其在机械设计中的应用．重庆：重庆大学出版社，1990：2~26
[106] 林腾蛟．核压力容器法兰密封分析三维计算机程序系统[学位论文]．重庆大学硕士论文，1992：32~40
[107] 陈如欧，胡忠民．塑性有限元及其在金属成形中的应用．重庆：重庆大学出版社，1989：24~50
[108] 李大潜，李瑞遐，沈玮熙，谭永基，郑宋穆. 有限元素法续讲. 北京：科学出版社. 1979：175~214
[109] 黄振中．求解瞬时温度场的有限元显式算法．应用数学和力学，1987, (6):5~8
[110] K H Huebner. The Finite Element Method for Engineers. John Wiley & Sons Inc., 1982：29~58
[111] 谢贻权，何福保．弹性和塑性力学中的有限单元法．北京：机械工业出版社，1981：102~160
[112] 马书尧，郭占科. 有限元软件中的 VC＋＋和 FORTRAN混合语言编程. 河北工业大学学报. 2001，30(6)：21~26
[113] 聂义勇. 有限元工程化的几项关键技术. 小型微型计算机系统. 1997，18(11)：7~12
[114] Ho-Lek. Finite Element Mesh Generation Methods: a Review and Classification. Computer Aided Design. 1980, (20)：27~38
[115] 黄志超，包忠诩，周天瑞，陈泽中. 有限元网格划分技术研究. 南昌大学学报(工科版). 2001，23(4)：25~31，44
[116] Martins P A F, Barata Marques M J M. Model 3——a Three-Dimensional Mesh Generation. Computer & Structures, 1992, 42(4)：511~529
[117] Patrick Knupp, Atanly Steinberg. Fundamentals Grid Generation. CRC Press,1993：22~103
[118] Berry K J. Parametric 3D Finite Element Mesh Generation. Computers & Structures, 1989, 33(4)：969~976
[119] Cavendish J C, Field D A, Frey W H. An Approach to Automatic Three Dimension Mesh Generation. International Journal for Numerical Methods in Engineering, 1985, 21：329~347
[120] 丁永祥. 约束Delaunay三角剖分与有限元网格自动网格生成. 华中理工大学学报. 1995，(6)：39~43
[121] 杨伟军，包忠诩，柳和生. 三维表面有限元网格生成的修正Delaynay法. 金属成型工艺. 1999，(1)：42~44
[122] Blacker Ted D. Paving: a New Approach to Automated Quadrilateral Mesh Generation. International Journal for Numerical Methods in Engineering, 1991, 32：811~847
[123] 王世军，黄玉美，张广鹏. 一种全四边形有限元网格生成方法——堆砌法. 机械工程学报. 2000，(10)：72~75
[124] Shephard Mark S. Automatic Three Dimensional Mesh Generation by the Finite Octree Technique. International Journal for Numerical Methods in Engineering, 1971, 20：709~749
[125] 魏红宁，周本宽. 基于叉树法的自适应有限元局部网格加密研究. 计算力学学报，1997，14(3)：324~331
[126] Yerry, M A, Shephard, M S. A modified quadtree approachto finite element mesh generation. IEEE Comput. Graph. & Appl. 1983, (2)：39~46
[127] Lai Mingwu, Benzley Steven E, Sjaardema Greg. A Multiple Source and Target Sweeping Method for Genarating All Hexahedral Finite Element Meshes. 5th International Meshing Roundtable. Sandia National Laboratories, 1996：217~228
[128] Robert Schneiders, Rolf Bunten. Automatic Generation of Hexahedral Finite Element Meshes. Computer Aided Geometric Design. 1995, (12)：693~707
[129] 傅沛福，吴淑芳，胡平，李运兴，柳玉启. 全自动自适应网格细化. 计算力学学报. 1997，14(1)：114~118
[130] 吴淑芳，傅沛福. 网格自动重分及细化技术. 塑性工程学报. 1997，(4)：70~73
[131] 丁俊，唐国安. 平面有限元网格生成的线段转换法. 计算机辅助设计与图形学报. 1998，
(10)：41~45
[132] 王璋奇，安利强，王孟. 有限元分析结果的可视化处理方法. 工程图学学报. 2002，(1)：72~77
[133] 武君胜. 有限元分析计算中的可视化技术及实现. 计算机辅助工程. 1995，4(1)：47~51
[134] 谢春，马兰，彭颖红. 有限元计算结果的可视化处理. 计算机辅助设计与图形学报. 2000，12(2)：81~84
[2] 陈广. 上海核电国产化优势初步分析. 工业技术进步. 1997，(3)：13~15
[3] 吴宗鑫，张作义. 世界核电发展趋势和高温气冷堆. 核科学与工程. 2000，20(3)：211~219，231
[4] 压水堆压力容器法兰密封技术专题报告．中国核动力研究设计情报室，1990
[5] D H Van Campen, P J Deen, D G H Latzko. Deformation of Large-diameter, High-pressure Vessel Flanges. Proc. 1st Int. Conf. on Pressure Vessel Tech., Delft, 1969
[6] C Menken. Influence of Bolt Loading on Deformation of Pressure Vessel Flanges. Proc. 1st Int. Conf. on Pressure Vessel Tech., Delft, 1969
[7] N W Murray, D G Stuart. Proc. Symp. on Pressure Vessel Research towards Better Design. Int. Mech. Eng., London, 1970
[8] D H Van Campen. Some Applications of Engineering Mechanics in Pressure Vessel Analysis Thesis. Delft, 1972
[9] Z Zudans, T C Yen, W H Steigelmann. Thermal stress Techniques in the Nuclear Industry. New York, American Elsevier Publishing Company, 1975
[10] A bakker, H A C M Spass. A Computational Procedure for the Stress and Deformation Analysis of High-pressure Vessel Flanges. Int. J. Pres. Ves. & Piping. 1977, (5)：135~149
[11] H A C M Spass. Deformation Behavior of Large High-pressure Vessel Flanges Development in Stress Analysis for Pressurised Components. R W Nichols (ed.). Elsevier Applied Science Publishers, London, 1977
[12] H Sacher, M Mayk. Structure Analysis of a Reactor Pressure Vessel by Two- and Three- Dimensional Models. Int. J. Pres. Ves. & Piping. 1982, (10)：147~150
[13] D H Cascales, C Militello. An Accurate Simple Model to Evaluate the Integral Flange Rotation. Int. J. Pres. Ves. & Piping. 1987, (30)：54~65
[14] D H Cascales, C Militello, W J Mulhall. Bolted Flanged Joint Leakage Analysis with Non-linear Gasket Behavior. Int. J. Pres. Ves. & Piping. 1987, (30)：78~92
[15] Ceorge F. P.E.Leon. An overview of the US PVRC Research Program on Bolted Flanged Connections, Pressure Vessel Technology Vol. 1 Design & Analysis, Pergramon Press 1988：32~65
[16] Gu Fangyu, Zeng Fanmao. Pretensioning of Pressure Vessel Bolts as a Problem Including
Elastic Contact with Adjustable Gaps. Int. J. Pres. Ves. & Piping. 1985, (19)：12~25
[17] 谷芳毓，曾凡茂，程中易．可调间隙弹性接触问题．核动力工程. 1982, (10)：25~29
[18] 曲家棣．核容器瞬态热分析、冷态密封分析程序系统．上海核工程研究设计院．1985
[19] 李润方，陈如欣，李水安．反应堆压力容器法兰密封系统密封分析程序系统理论文本．重庆大学，中国核动力研究设计院．1991
[20] 李润方，林腾蛟．反应堆压力容器局部三维冷热态密封分析．机械工程学报，1993，29(6):104-108．
[21] 李润方，李水安，李健，陈骏，林腾蛟. 核容器法兰密封系统密封性能数值分析. 核动力工程. 1994，15(1)：12-18．
[22] 曲家棣，盛贤科，王佩珠. GH169度银O形环冷、热态性能测试. 压力容器. 1997，14(3)：183~188
[23] 曲家棣，王佩珠，窦一康，盛贤科，谢世球，金惠琴. 压力容器热态密封实验与分析. 1997，14(4)：1~7
[24] 彭向和，高芝晖，曾祥国. O形密封环的弹塑性大变形和循环松弛分析. 重庆大学学报. 1998，21(2)：46~50
[25] 于溯源，何树延，刘俊杰. HR-200压力容器主法兰密封的分析法设计. 第十一届全国反应堆结构力学会议，成都，2000, 10
[26] 左卫东，于溯源，刘俊杰等. 10MW高温气冷堆压力容器主法兰结构的有限元接触分析. 核动力工程. 2001，22(3):226-231.
[27] R W Clough. The Finite Element Method in Plane Stress Analysis. Proceedings of the Second Conference on Electronic Computation. American Society of Civil Engineers, New York, 1960
[28] C H Lee, S Kobayashi. Elastoplastic Analysis of Plane stain and Axisymmetric Flat Punch Indentation by the Finite Element Method. Int. J. Mech. Sci.. 1970, (12)：25~38
[29] J T Oden, D R Bhandari, O Yagawa, T J Chang. A New Approach to the Finite Element Formulation and Solution of a Class of Problem in Coupled Thermoelasto-viscoplasticity of Crystalline Solids. Nuclear Engineering Design. 1973, (24)：135~152
[30] D R J Owen, E M Salonen. Three Dimensional Elasto-plastic Finite Element Analysis，Int. J. Num. Mech. Engng. 1975, (19)：248~265
[31] Hindon, Owen. Finite Element Programming. Academic Press, 1977：56~78
[32] Owen, Hinton. Finite Element in Plasticity. Pineridge Press Limited, 1980：102~150
[33] B M Iron. Engineering Applications of Numerical Integration in Stiffness Methods. AIAAJ. 1966, 4(11)：378~392
[34] B M Irons, A Frontal Solution Program. Int. J. Num. Mech. Engng. 1970, 2：456~473
[35] E A Wilson, B Parsons. Finite Element Analysis of Elastic Contact Problems using Differential Displacements. Int. J. Num. Mech. Engng. 1970, 2：248~256
[36] S K Chan, I S Tuba. A Finite Element Method for Contact Problems of Solid Bodies, II, Application to Turbine Blade Fastening. Int. J. Mech. Sci. 1971, 13：132~145
[37] S Ohte. Finite Element Analysis of Elastic Contact Problems. Bull. JSME 16, 1973
[38] A Francavilla, O C Zienkiewicz. A Note on Numerical Computation of Elastic Contact Problem. Int. J. Num. Mech. Engng. 1975, 19：913-924
[39] B Fredriksson. Finite Element Solution of Surface Nonlinearities in Structural Mechanics with Special Emphasis to Contact and Fracture Mechanics Problem. Computer and structure. 1976, 6：514~530
[40] 冈本纪明．有限要素法たょゐ非线性接触问题の解析. 日本机械学会论文. 1977, 372~386
[41] 茑纪夫．有限要素法の接触问题への应用．机械の研究. 1978, 2(8)：28~43
[42] J T Stadter, R O Weiss. Analysis of Contact through Finite Element Gaps. Computer and structure. 1979, 10(6)： 867-873
[43] 陈万吉．用有限元混合法分析弹性接触问题．大连工学院学报．1979, (2)：12~15
[44] 陈万吉．耦合热弹性接触问题有限元分析．固体力学学报．1981, (3)：17~21
[45] R F Li. Finite Element Analysis for Elastic-plastic Contact Problems with Friction. C/R/523/85, 1985
[46] 郭小明. 岩土介质弹塑性耦合问题的数学规划模型. 东南大学学报. 1994，24(5)：101~105
[47] 李锡夔. 蠕变-弹塑性-损伤耦合响应的有限元分析. 计算结构力学及其应用. 1995，12(1)：26~38
[48] 张承柱，刘信声. 应变空间表述的混凝土弹塑性耦合本构模型. 清华大学学报. 1996，36(3)：59~64
[49] 倪强，王乘，黄玉盈，张标标. 焊接过程的耦合三维热弹塑性有限元分析. 华中理工大学学报. 1997，25(1)：84~86
[50] 吕安军，曹志远. 地下结构和岩土介质弹塑性耦合分析的摄动半解析法. 应用数学和力学. 1998，19(4)：335~339
[51] 孙秦，徐晓飞. 弹塑性耦合张量的宏细观对偶形式. 应用力学学报. 1998，15(4)：19~21
[52] 张军，吴昌华. 轮轨接触问题的弹塑性分析. 铁道学报. 2000，22(3)：16~21
[53] 栾茂田，黎勇. 非连续变形计算力学模型弹塑性数值分析方法及其应用. 岩土工程学报. 2001，23(5)：530~535
[54] 高华林，蔡嗣经，童光煦. 有限单元法微机前处理程序的开发. 矿业研究与开发. 1994，14(2)：58~60
[55] 胡于进，王坚，陆金桂，周济. 任意形状平面域不完全Delaunay三角化研究. 华中理工大学学报. 1994，22(10)：4~8
[56] 陈大宏，王革，邹家祥. CAD技术在有限元分析前处理中的应用. 1994，16(增刊)：70~74
[57] 张广健，黄由玲. 用CAGD技术进行有限元前处理. 工程力学. 1995，12(4)：124~130
[58] 杜群贵. 三维实体有限元网格自动 Delaunay剖分. 华南理工大学学报(自然科学版). 1996，24(9)：46~50
[59] 安晓卫，郝长忠. 旋转结构的有限元网格自动生成. 机械设计. 1996，(12)：40~42
[60] 安晓卫，张华. 用于不规则几何形体的有限元网格自动生成处理技术. 沈阳工业学院学报. 1997，16(4)：82~85
[61] 赵铁成，沈月芬，朱国桢，许先壁. 锅炉锅筒下降管三通的三维有限元自动剖分. 锅炉技术. 1997，(6)：10~14
[62] 林腾蛟，李润方，陈兵奎，林超. 混合标架下宽斜齿轮三维接触有限元前处理技术. 机械传动. 1997，21(2)：9~12
[63] 吴淑芳，李占国. 二维有模成形自适应数值模拟系统. 塑性工程学报. 1998，5(3)：22~26
[64] 雷震宇，黄洪钟. 160t铁路救援起重机转台有限元建模专家系统数据库的研究. 建筑机械. 1998，(5)：21~24
[65] 饶志伟，赵洪伦. 车体建模CAD系统的设计. 计算机辅助设计与制造. 1998，(2)：28~29
[66] 王琳，周理兵，马志云. 一种适于大型水轮发电机磁场计算前处理的新方法. 1998，18(1)：39~41
[67] 林腾蛟，李润方，陈兵奎，李亦彦. 鼓形齿联轴器轮齿三维有限元模型自动生成及强度计算. 1998，11(3)：18~19
[68] 左正兴. 一种简单高效的有限元网格优化算法及应用. 机械设计. 1998，(8)：21~24
[69] 查德根，沈景凤，黄宗益. 焊接结构件有限元分析中的前处理技术. 上海水产大学学报. 1998，7(4)：340~344
[70] 虞春，周雄辉，张永清. 特征造型与有限元分析的集成化研究. 计算机辅助设计与图形学学报. 1999，11(1)：62~65
[71] 安晓卫，张华，钟震宇. 三通管件的有限元网格的自动剖分. 沈阳工业学院学报. 1999，18(3)：5~10
[72] 段红杰，阎德元. 复合材料整体异形构件的CAD技术研究. 北京航空航天大学学报. 1999，25(5)：585~587
[73] 王璋奇，安利强. 基于SuperSAP的有限元前处理接口程序设计. 华北电力大学学报. 1999，26(1)：74~77
[74] 范峰. 空间网格结构CAD软件ASSAP 1系统研制与开发. 哈尔滨建筑大学学报. 1999，32(5)：16~19
[75] 倪其民，钱瑞明，王水来，阮雪榆. 回转支承智能设计系统研究与开发. 机械科学与技术. 1999，18(5)：709~711
[76] 秦浦雄，王书庆. 桥梁结构的三维建模技术. 计算机辅助工程. 2000，(2)：46~52
[77] 谭也平. 有限元前处理系统的图形集成实现. 计算机辅助工程. 2000，(1)：53~58
[78] 潘子杰，杨文通，郭九生. 任意平面区域有限元三角形网格全自动剖分. 北京工业大学学报. 2001，27(4)：493~496
[79] 陈建梅，关湛湛. 电器磁场计算可视化前处理软件开发研究. 三峡大学学报(自然科学版). 2001，23(1)：35~37
[80] 安利强，孟宪国，王璋奇. 转子热应力有限元前后处理方法. 华北电力大学学报. 2001，28(1):68~71
[81] 王璋奇，安利强，彭震中. 基于参数化建模的转子有限元剖分. 热能动力工程. 2001，16(9)：538~539
[82] 马智慧，胡兆吉，李鸣. 板料成形中的载荷及位移边界条件的处理. 化工装备技术. 2002，23(2)：29~31
[83] 王璋奇，安利强，王孟. 有限元分析结果的可视化处理方法. 工程图学学报. 2002，(1)：72~77
[84] 李建国，刘实，李殿中，戎利建，范寸，李依依. 汽车纵梁冲压成形及回弹过程的计算机模拟. 塑性工程学报. 2001，8(2)：4~7
[85] W. E. Lorence, H. E. Cline. MarchingCube: A high resolution 3D surface construction algorithm. Computer Graphics(ACM Siggraph Proceedings). 1988, 22(4)：51~58
[86] 程世才，林国成，满波. 实用有限元后处理程序的研制及开发. 哈尔滨工业大学学报. 1996，28(6)：45~49
[87] 仲梁维. 三维有限元后处理程序设计. 微型电脑应用. 1997，1：62~64，67
[88] 费少梅. 有限元计算结果的可视化. 内燃机工程. 1997，17(4)：26~28
[89] 王乐天，谢永慧，安宁，孟庆集. 汽轮机叶片三维有限元分析后处理. 西安交通大学学报. 1997，31(12)：33~38
[90] 林鹏，卓家寿，裴安荣. 等值面绘制及其工程应用. 西安理工大学学报. 1998，14(4)：413~416
[91] 虞春，周雄辉，张永清. 基于特征的智能有限元后处理. 计算机辅助工程. 1999，(1)：19~24
[92] 杨晓松，顾云宪，李云鹏，关振群. 三维有限元模型的任意剖切及其等值线与彩色云图生成的方法. 中国图象图形学报. 1999，4(7)：574~578
[93] 刘理，刘土光，张标标，李天匀. 有限元分析后处理程序中的等值面生成算法. 华中理工大学学报. 1999，27(7)：80~82
[94] 闵作兴，陈国华，吴衍平. 发动机有限元可视化应用软件研究. 华中理工大学学报. 1999，
27(11)：57~58
[95] 孙黎明，李世其. 有限元后处理中的一种快速多等值面面绘制算法. 武汉纺织工学院学报. 1999，12(1)：69~73
[96] 杨晓松，原仓周，李云鹏，关振群，顾云宪. 基于scan-buffer三维不规则数据场的直接体绘制技术. 大连理工大学学报. 2000，40(2)：198~202
[97] 苏海东. 有限元后处理图形软件FEPS的开发. 长江科学院院报. 2000，17(3)：19~22
[98] 李昱明，周飞，郭新贵，彭颖红. 八节点六面体有限元网格自动消隐技术. 计算机工程与应用. 2000，(11)：89~91
[99] 罗婷，张颖江，陈雯. 利用OpenGL实现有限元后处理. 湖北工学院学报. 2000，15(4)：7~9
[100] 詹梅，刘郁丽，杨合. 三维有限元模拟等值线生成的新方法. 重型机械. 2000，(1)：36~38
[101] 章勤，李品，苗京. 一种等值线图绘制算法及其并行实现. 华中科技大学学报. 2001，29(9)：4~6
[102] 诸葛良，刘健，王明武. 有限元后处理的一种动态可视化方法. 机械设计与研究. 2001，17(1)：29~30
[103] 李昱明，周飞，郭新贵，彭颖红. 基于六面体八节点单元的三维有限元数值模拟表面等值线的生成算法. 锻压技术. 2001，(2)：10~13
[104] 杨晓松，李云鹏，顾云宪. 基于Slicing的快速有限元网格体绘制算法. 计算机辅助设计与图形学学报. 2001，13(4)：348~352
[105] 李润方，龚剑霞．接触问题数值方法及其在机械设计中的应用．重庆：重庆大学出版社，1990：2~26
[106] 林腾蛟．核压力容器法兰密封分析三维计算机程序系统[学位论文]．重庆大学硕士论文，1992：32~40
[107] 陈如欧，胡忠民．塑性有限元及其在金属成形中的应用．重庆：重庆大学出版社，1989：24~50
[108] 李大潜，李瑞遐，沈玮熙，谭永基，郑宋穆. 有限元素法续讲. 北京：科学出版社. 1979：175~214
[109] 黄振中．求解瞬时温度场的有限元显式算法．应用数学和力学，1987, (6):5~8
[110] K H Huebner. The Finite Element Method for Engineers. John Wiley & Sons Inc., 1982：29~58
[111] 谢贻权，何福保．弹性和塑性力学中的有限单元法．北京：机械工业出版社，1981：102~160
[112] 马书尧，郭占科. 有限元软件中的 VC＋＋和 FORTRAN混合语言编程. 河北工业大学学报. 2001，30(6)：21~26
[113] 聂义勇. 有限元工程化的几项关键技术. 小型微型计算机系统. 1997，18(11)：7~12
[114] Ho-Lek. Finite Element Mesh Generation Methods: a Review and Classification. Computer Aided Design. 1980, (20)：27~38
[115] 黄志超，包忠诩，周天瑞，陈泽中. 有限元网格划分技术研究. 南昌大学学报(工科版). 2001，23(4)：25~31，44
[116] Martins P A F, Barata Marques M J M. Model 3——a Three-Dimensional Mesh Generation. Computer & Structures, 1992, 42(4)：511~529
[117] Patrick Knupp, Atanly Steinberg. Fundamentals Grid Generation. CRC Press,1993：22~103
[118] Berry K J. Parametric 3D Finite Element Mesh Generation. Computers & Structures, 1989, 33(4)：969~976
[119] Cavendish J C, Field D A, Frey W H. An Approach to Automatic Three Dimension Mesh Generation. International Journal for Numerical Methods in Engineering, 1985, 21：329~347
[120] 丁永祥. 约束Delaunay三角剖分与有限元网格自动网格生成. 华中理工大学学报. 1995，(6)：39~43
[121] 杨伟军，包忠诩，柳和生. 三维表面有限元网格生成的修正Delaynay法. 金属成型工艺. 1999，(1)：42~44
[122] Blacker Ted D. Paving: a New Approach to Automated Quadrilateral Mesh Generation. International Journal for Numerical Methods in Engineering, 1991, 32：811~847
[123] 王世军，黄玉美，张广鹏. 一种全四边形有限元网格生成方法——堆砌法. 机械工程学报. 2000，(10)：72~75
[124] Shephard Mark S. Automatic Three Dimensional Mesh Generation by the Finite Octree Technique. International Journal for Numerical Methods in Engineering, 1971, 20：709~749
[125] 魏红宁，周本宽. 基于叉树法的自适应有限元局部网格加密研究. 计算力学学报，1997，14(3)：324~331
[126] Yerry, M A, Shephard, M S. A modified quadtree approachto finite element mesh generation. IEEE Comput. Graph. & Appl. 1983, (2)：39~46
[127] Lai Mingwu, Benzley Steven E, Sjaardema Greg. A Multiple Source and Target Sweeping Method for Genarating All Hexahedral Finite Element Meshes. 5th International Meshing Roundtable. Sandia National Laboratories, 1996：217~228
[128] Robert Schneiders, Rolf Bunten. Automatic Generation of Hexahedral Finite Element Meshes. Computer Aided Geometric Design. 1995, (12)：693~707
[129] 傅沛福，吴淑芳，胡平，李运兴，柳玉启. 全自动自适应网格细化. 计算力学学报. 1997，14(1)：114~118
[130] 吴淑芳，傅沛福. 网格自动重分及细化技术. 塑性工程学报. 1997，(4)：70~73
[131] 丁俊，唐国安. 平面有限元网格生成的线段转换法. 计算机辅助设计与图形学报. 1998，
(10)：41~45
[132] 王璋奇，安利强，王孟. 有限元分析结果的可视化处理方法. 工程图学学报. 2002，(1)：72~77
[133] 武君胜. 有限元分析计算中的可视化技术及实现. 计算机辅助工程. 1995，4(1)：47~51
[134] 谢春，马兰，彭颖红. 有限元计算结果的可视化处理. 计算机辅助设计与图形学报. 2000，12(2)：81~84