冷拔钢管与焊管成型过程的非线性有限元分析

1、分类号：密级：天津理工大学研究生学位论文冷拔钢管与焊管成型过程的非线性有限元分析（申请硕士学位）学科专业：机械设计及理论研究方向：机械强度分析及现代设计理论作者姓名：胡建英指导教师：叶金铎年月摘要摘要为了获得不同模具，不同材料模型，不同强化规律下冷拔钢管成型过程中的应力应变分布规律和对成型精度的影响规律，以及单半径和双半径焊管成型过程的应力应变分布规律，本文利用弹塑性大变形和接触非线性有限元法，应用大型多物理场有限元分析软件对空拔管、芯拔管以及单半径和双半径焊管的第一架成型的成型过程进行了有限元仿真，其中考虑了模具的变形。得到的主要结论如下：（）通过对锥模与弧形模下空拔管成型过程的数值模拟，发

2、现在成型直径上弧通过研究不同材料模型和强化规律对空拔管拔制力和尺寸精度的影响，发现通过对锥模与弧形模下芯拔管成型过程的数值模拟，发现弧形模成型比锥模提出了一种数值模拟和实验杂交的实验方法，测量了钢管内表面的残余应形模比锥模成型精度高，在成型壁厚上锥模比弧形模成型精度高；（）真实材料模型下利用随动强化理论进行的仿真总体效果要好一些；（）成型的拔制力与芯棒后拉力大，锥形模成型脱模后钢管的直径与壁厚更接近设计值。（）力。得到了残余应力对液压油缸缸筒承载能力的影响规律。（）通过对单半径和双半径焊管第一架成型过程的数值模拟，发现二者成型过程中的应力应变分布很不均匀，可以考虑多半径成型方式，最后对多辊下的

3、焊管成型过程的仿真提出了研究方向。关键词：非线性有限元法冷拔钢管高精度模具残余应力焊管成型辊杂交实验方法，鹳啦蛐畦啦斌妣吣姚筒呶谳僦啦：，（）（），（），仃（），：独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作和取得的研究成果，除了文中特别加以标注和致谢之处外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含为获得墨洼堡墨盘至或其他教育机构的学位或证书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示了谢意。学位论文作者签名：钢遣寝签字日期：如年月日学位论文版权使用授权书本学位论文作者完全了解墨盗墨兰盘至有关保留、使用学位论文可以

4、将学位论文的全部或部分内容编入的规定。特授权墨盗墨墨盘鲎有关数据库进行检索，并采用影印、缩印或扫描等复制手段保存、汇编，以供查阅和借阅。同意学校向国家有关部门或机构送交论文的复本和电子文件。学位论文作者签名：撼签字日期：如年月户日（保密的学位论文在解密后适用本授权说明）导师基名：叶令碍签字日期：一占年月厂日第一章绪论第一章绪论课题的提出钢管是经济建设中重要的原材料之一，素有工业“血管”之称。精密钢管广泛地被用作机械工业用管（主要是液压与气动缸体管、液压输油管、汽车传动轴管和半轴套管、轴承钢管、皮带传送托辊管、印染辊管等）、石油地质工业用管（油泵管、钻探管、油管、套管、钻杆等）、化工用管（石油裂

5、化管、化肥用高压管、化工设备及管道用管）、电站锅炉和热交换器管等。随着我国经济建设的不断发展和技术的不断进步、精密钢管的应用范围会不断扩大、需求量也会越来越大。高精度精密冷拔钢管广泛地应用于国内液压、气动缸，煤炭井下支架（支柱），石油泵管，千斤顶等制造领域。高精度精密冷拔钢管的推广应用对节约钢材，提高加工工效，节约能源，减少液压缸、气缸加工设备投资有重要意义。焊管的使用领域涉及到石油钻采中所用的套管，海洋石油工业中用于石油、天然气输送的海底油气输送管干线用油气输送管，配气管等。日本、美国、德国在石油焊管技术的应用方面处于领先地位，已有十几年的应用经验，年代首钢等单位虽然开展了焊管应用于油套管的

6、研究，但由于工艺水平差和设备落后，一直未能形成大批量生产能力，目前，我国仍在大量进口油套管。因此建立我国自有的油井焊管生产线，不仅能解决油套管钢管的国产化问题，同时可降低油套管的成本。根据国外资料，焊管比无缝管销售价格低。冷拔钢管的生产与研究现状冷拔是继热轧、挤压或焊接之后对管材进行二次加工以生产精密、薄壁和高机械性能产品的主要方法。冷拔包括空拔和芯拔，这种成型方式具有工艺设备简单，成型效果好的优点，在军工、机械和石油等领域都得到了广泛的应用。用冷拔的方法可以生产外径为毫米、壁厚毫米的钢管和外径毫米的有色金属管。冷加工和相应的热处理相结合可使管材获得高的综合机械性能，除了一般圆管以外还可以生产

7、多种异形管和变断面管。冷拔钢管的生产现状冷拔是在外加拉力作用下，迫使金属坯料通过模孔，以获得相应形状与尺寸制品的塑性加工方法。钢管冷拔的基本方式有：（）空拔。拉拔时管坯内部不放芯头，通过模孔后外径减小，壁厚一般会略有变化。）长芯棒拉拔。将管坯自由地套在表面抛光的芯杆上，使芯杆与管坯一起拉过模孔，以实现减径和减壁，芯杆的长度应略大于管坯的长度，在拉拔一道次之后，需要用脱管法或滚轧使之扩径的方法取出芯杆。（固定短芯头拉拔。拉拔时将带有芯头的芯杼固定，管坯通过模孔实现减径和减第一章绪论壁。（）游动芯头拉拔。在拉拔过程中，芯头靠本身所特有的外形建立起来的力平衡被稳定在模孔中。与固定短芯头拉拔相比，游动

8、芯头拉拔难度大，工艺条件和技术要求较高，配模有一定限制，故不可能完全取代固定短芯头拉拔。（）扩径拉拔。管坯通过扩径后，直径增大，壁厚和长度减小。这种方法主要是在设备能力受到限制而不能生产大直径的管材时采用。近十年来，在研究许多新的拉拔方法的同时，展开了高速拉拔的研究，成功的制造了许多高速连续拉拔机、多线链式拉拔机和圆盘拉拔机。高速拉线机的拉拔速度达到；圆盘拉拔机可生产毋以下的管材，最大圆盘直径为，拉拔速度可达，最大管长为以上；多线链式拉拔机一般可以自动供料、自动穿模、自动套芯杆、自动咬料和挂钩、管材自动下落以及自动调整中心。另外还有管棒材成品连续拉拔矫直机列，在该机列上实现了拉拔、矫直、抛光、

9、切断、退火以及探伤等。冷拔钢管的国内外研究现状冷拔钢管成型过程十分复杂，属于弹塑性大变形和接触非线性相互耦合的多重非线性问题，长期以来对冷拔钢管成型过程的研究还落后与生产实际，模具工艺参数的设计和冷拔工艺的设计主要采用经验设计方法，设计周期长，费用高，产品尺寸精度低，不能满足对精密冷拔钢管需要。国内外已有的研究成果多数是在各种假设的基础上设计模具和工艺，研究工作多数针对单一因素，如拔制力、钢管成型直径或壁厚等，对冷拔成型过程的研究工作相对较少。在冷拔钢管成型过程的数值模拟方面，考虑弹塑性大变形的工作较多，考虑接触非线性的工作较少，考虑模具变形及获取管坯在成型过程的应变及应力分布及脱模后管材残余

10、应力的工作更少。年，用弹塑性有限元分析固定楔型芯棒冷拔钢管的拔制过程，建立了应力分布和残余应力模型；年，使用刚塑性有限元模型分析了空拔管生产过程中管壁的变化与应变分布；年，采用刚塑性有限元模型和实验数据对空拔管过程的应变状态、应变分布与应力分析进行评价；年，颜鹏等人使用二维弹塑性有限元法，分析了空拔管纵裂与残余应力的关系，并且指出机械加工是一种消除残余应力的有效方法；年，杜风山等人使用三维弹塑性有限元法模拟了钢管空拔，研究了管材稳态与非稳态变形；年，黄成江等人使用二维弹塑性有限元法分析了无芯棒拉拔的均壁作用；年，黄成江等人使用三维弹塑性有限元法比较了空拔时，单道次拉拔和多道次拉拔的拔制力、残余

11、应力，并且分析了拉拔过程的热影响：年，叶金铎等人使用二维弹塑性有限元法模拟了空拔管成型，分析了钢管的变形过程、残余应力。已有的研究工作主要分为两类：一是采用各种假设导出拔制力和拔制后钢管直径的公式，也采用数表的方法确定钢管的直径；二是采用有限元的方法通过研究钢管的成型过程，获取拔制力和钢管的直径。在数值模拟方面已有的研究工作，主要集中在对锥模半锥角、摩擦系数等参数的研究。在非线性因素的考虑方面，多数的研究工作没有考虑接触非线性和模具的变形。所用模具的孔型只有苏式模和中式模。第一章绪论焊管的生产与研究现状年代以来焊管机组的建设一直很活跃，由于不断地采用新技术，焊管的焊缝质量大大提高，焊接区的强度

12、、韧性等指标甚至可以达到或超过母材。正因为如此，焊管在继续占据着水、煤气输送、建筑结构、电线套管等较低强度的钢管市场外，在锅炉、石油油井、化工设备、机械构件、汽车工业以及高压流体输送等领域广泛地与无缝钢管进行市场竞争，其市场份额越来越大。焊管用钢也不只是原来的碳素钢，而向合金钢和不锈钢等高质钢发展焊管的生产现状（）小口径焊管生产机组小口径焊管机组（咖以下机组）的成型方法众多，最早是导板成型法和对辊成型法（含边缘弯曲法、中心弯蛆法和圆周弯曲法），年代末引入了排辊成型法。进入年代，为了提高产品质量、扩大产品规格范围，先后开发出众多新的成型法，主要有：可生产薄壁管的立辊成型法（法）和直缘式排辊成型法

13、；生产厚壁管的带有边部成型机架的组装式成型法；生产极厚壁管的复合成型法；既可生产薄壁管又可生产厚壁管的弯曲辊成型法和柔性成型法但法）；生产极薄壁管的无辊自然成型法以及无鼓肚变形辊成型法（法）等。近年来除导板成型法外，所有的成型法都在采用，生产的最薄焊管为，最厚焊管为，生产钢种为普碳钢、合金钢及不锈钢。另外，为提高焊管机组的轧辊共用性、提高带钢边部成型性能及质量稳定性，新日铁公司设计和开发出把上下轧辊以垂直轴为中心进行旋转和交叉，以适应各种不同弯曲率和壁厚间隙形状的轧辊交叉边部成型法（法。新成型法开发的主要目的就是使成型过程中的带钢变形自然流畅、合理均匀，边缘受力条件良好。薄壁管成型主要解决边缘

14、延伸问题，厚壁管成型则主要要求带钢边缘达到平齐，并充满孔型。（）中径焊管生产机组中径焊管机±（最早也是采用对辊成型法，但自年代出现了排辊成型法后，排辊成型法己成为中径焊管机组的主流。年代后除台湾耕安公司在中国新建的套机组采用对辊成型法和日本为美园和欧洲提供的套机组采用成型法外，几乎无例外地采用排辊成型法。中径焊管机组的排辊成型法主要有三种工艺：全排辊成型法、半排辊成型法和直缘式排辊成型法。其制造厂家主要有法国、奥地利、德虱公司，最近年内世界新建了套机组即分别由这三家公司制造。这些机组生产的主要品种为管线管、结构管及套管等。（大口径焊管生产机组生产油气输送管道用大口径直缝埋弧焊管（）生

15、产机组的生产方式主要有法、法、成型法、法、瑚法以及全排辊成型法等。各种成型法中以法生产的钢管质量最好、可靠性最高，但由于机组建设投资巨大，进人年代以后已不再新建。而原来主要用于生产大口径结构管（外径最大！）的辊弯成型工艺，经过一系列改造，生产的直缝埋弧焊管在使用性能和可靠性方面有较大的提高，因而被用来生产油气输送用管。由于其生产灵活性大，设备投资？，第一章绪论机组低许多倍，在一些发展中国家得到推广应用，这就是目前已被广泛认可的法。据统计，年代以来国外建设的大口径直缝埋弧焊管生产机组的套，其中有套为机组。另外，近年来德国公司开发了即成型法用于生产大口径直缝埋弧焊管，并成功地应用于家工厂之中（其中

16、之一可以生产大口径厚壁不锈钢管）。大口径螺旋埋弧焊管生产机组近年来仍在一些国家建设，不但用于生产结构用管，而且也用于生产油气输送用管线管。生产的油气输送管最大钢级已达到。焊管的国内外研究现状焊管属于辊式冷弯成型，与无缝钢管相比，具有设备投入小，钢管尺寸精度高的优点，是一种用途广泛的高效型材。在辊式成型过程中，带钢经一系列轧辊（平辊、立辊）连续弯曲经焊接成为钢管。在成型过程中，带钢要经过非常复杂的弹塑性大挠度弯曲变形，同时伴有中性面的挤压，是横向变形和纵向变形综合的物理、几何非线性问题，因此其理论解析比较困难目前还没有统一的数学模型，也没有完整的数值模拟研究工作。目前，成型过程的分析方法主要分为

17、五类，即整体几何分析法，实验法，能量法”，有限元法和法。整体分析法。整体分析法是不考虑成型辊对带钢成型的局部影响，以成型区中的管坯整体为研究对象。它通过比较各种成型方法中管坯边缘的总体延伸来鉴别这些方法的优劣。实验研究法。实验研究法就是对于所要求解的问题不同设置不同的实验装置。日本有九个工厂用厚度为，毫米的带钢焊接直径毫米的钢管，下辊的调整分下山、上山和水平三种，并采用不同的孔型系统，不同的成型辊布景形式和不同的机架间距，由各厂分别进行试验。通过测量和计算，求出了钢管外表面的带钢变形量【。电阻应变片法在实验法测定中起到了很重要的作用。西德狄盖斯用电测法研究了焊管成型过程中的带钢变形，原苏联以及

18、其他国家，普遍采用电阻应变片来测量成型过程中的纵向变形。电阻应变片法能精确地测量成型过程中的变形，具有灵敏度高、惯性小的优点，尤其是，这种方法能远距离测量多个点的变形。天津理工大学的温殿英等用应变片的方法在液压万能试验机上研究了成型区内金属的应变。但是电测法不能测量管坯与轧辊接触以后钢管的变形。此外，加拿大等人还特意设计了模拟实验装置来研究电焊管成型的轴向应变分布。能量法。能量法是指通过建立含有若干待定系数的全能或全功率范函，然后对该能量或范函取极小化来确定这些参数，进而确定物体变形的分析方法。能量法最早是在年由伽古恩应用到物体变形分析中去的，他用能量方法分析了角材的连续弯曲。在确定板材成型变

19、形状态方面，最初由斯米尔诺夫弗可拉耶夫和古恩提出来的【”。古恩在分析中作了一些简化假定：（）管坯的横截面始终为平面并垂直于轧制轴线；（）垂直于轧制轴线的每个断面上的中性面宽度不变并等于带钢的初始宽度。通过联立拉格朗日坐标和欧拉坐标方程来确定主变形及变形强度、应力强度，进而计算内力的全部功，利用拉格朗日中值定理，并设外力功恒定，确定未知函数，最后求出内力功的极小值。文献的管坯变形分析工作是：（）假设管坯在变形区中的构形是己知的，它由成型工具的孔型所第一章绪论决定的横断面形状按设计规律均匀变化。（）建立拉格朗日坐标系和欧拉坐标系，用拉格朗日坐标系与带钢中性面联系起来，中性面上的每个点便与成型过程中

20、不变的两个数（，）相对应。将变形区的中性面计入欧拉坐标系（，）。引入两个第一类二次型曲面之差，即与成形区中性面相应的和与带钢中性面相应的之差可求出应变张量分量，叫，。这种分析方法属于全量型能量法，不能考虑变形历史的影响。目前日本学者木内学等人建立的数学模型应用比较广泛【】。现将其思想介绍如下：（）采用平断面假设并假设相邻两个平断面可以熬体的沿纵向伸缩，以使断面上纵向力的合力为零。相邻两机架间管坯的构形借助于形状函数（）（）描述：（）（）一（）“式中，一相邻机架间距；待定幂指数；（）描述了带材上一点由号轧辊轧槽上的点（，）向号轧辊轧槽上的点（：，：）流动时的空间轨迹。（一）（）。（一。）（）式中

21、（）是依赖于形状函数（）的非独立形状函数。（）（）将成型前的带材沿长度方向分割成具有初始长度的宽度带状单元，沿带状要素的宽度方向再分成若干等份。假设中性层与中央层一致，沿厚度方向均布且。：，。使用增量法可得到每个要素的应力和变形功表达式，从而得到整段坯料的总变形功的数学表达式，即：一击摹，【（、（。）“，。¨厚度方向的分层（），。塑性变形功增量（。）“，。弹性变形功增量通过对求极小化确定，这时的应力应变就是所求的解，由这时的节点坐标给出的变形曲面即为最佳变形曲面，认为该曲面是真实变形曲面。木内学等人开发的数学模型提供了直观、简明、实用的分析方法，该理论属于增量型能量法，能反映变形的历

22、史，更能反映塑性变形的本质，适合于计算机辅助设计。法。随着科学技术的不断进步，计算机辅助设计也得到了广泛的应用，同样计算机辅助设计也给我们对焊管的研究提供了一种有效的方法，这就是焊管成型法。在这方面，各国学者也做了很多尝试并取得了可喜的成果。奥纳艾特沃提出了用三次多项式表示的边缘变形高度口”。日本的通过收集几个公司的轧辊图和数据，并研究了管截面形状因子和轧辊孔型之间的关系，开发第一章绪论可用三次曲线示出：，（），（）任意阶段车辊的弯曲角；由下面方程给出：（一）【（）一（）】（）引入浮动因子后。使轧辊设计者进行人机对话：一（一）【（）一（）】（）通过这一方程形成一个简单通道截面，将其用于管截面，

23、这种方法属于数据库方法。德国的慕尼黑大学开发了软件系统【】，该系统提供有友好界面，可用模拟技术分析和优化辊弯成形工艺。他们采用一个与实际变形非常相近的几何形状作为模拟的模型，为优化不同成型区域孔型的数学近似几何模型，采用了有限元的分析计算，但该文并未涉及数学模型。美国的采用日本学者木内学的形状函数方法并加以完善，开发了限）软件。该软件在木内学理论膜应变的基础上加入了弯曲校正”】。有限元法。关于方法，目前很多国家的学者们正在作这方面的工作，并取得了一定的进展。蔡松庆等人首次尝试用板壳弹塑性有限元法分析直缝焊管成型过程，他们将成型区中的管坯看作是由一块平板在横向同时逐渐弯曲的结果。刘才、周瑛、韩志

24、武等用有限条法研究了冷弯型钢的成型过程”。徐树成通过对冷弯成型进行分析，用弹塑性大变形有限元理论进行了模拟，通过改写有限元程序，并对现场进行连续辊弯的一、二次成型实验。叶金铎等人将轧辊简化为平面曲线，采用样条有限条法研究了双半径焊管第一架的成型过程，与实验结果进行了比较，研究中假设钢管的材料为理想塑性。日本的小野田等人长期研究焊管成型的有限元模拟问题。他们使用刚塑性有限元分别对焊管管坯在开口孔型和闭口孔型中的变形进行研究，分别得到了管坯处于开口孔型和闭口孔型时的网格划分及边界条件。美国公司的，英国大学的以及加拿大大学的等人，他们使用商用大变形有限元软件包对冷弯成型进行分析，并将其应用于冷弯工具

25、的设计。等人使用非线性有限元软件对直缝焊管的成型应力，应变进行了分析，通过计算建立了预报焊管成型过程的边缘失稳模型，并得到调整轧机参数的依据。目前的对焊管成型过程的分析方法存在的问题有：整体分析法由于它没有考虑局部影响，因此不能提供各个点的应变、应力等全方位的信息。实验研究法只能建立成型过程不完整的近似物理模型。古恩将平断面假说应用于连续弯曲的成型情况是不恰当的，按圣维南原理，平断面假设仅适用于远离外力作用、接触变形区内的应力波及不到的地方，其次，平断面导致了成型区不准确的纵向变形图，这种第一章绪论图是单值的，除带钢中央纤维外，全部纤维处于拉紧状态，拉伸变形从中央向边缘增加。木内学等人开发的数

26、学模型采用一个参数来定变形曲面，误差会加大。蔡松庆等人采用的有限元分析方法从实际结果看，管坯中的应变和管坯的构形比较接近实际情况，但应力及外载与实际情况相差甚远【。刘才、周瑛、韩志武等用的有限条法仍属于有限元的范畴，是为了降低计算量而开发的，计算量降低了，其精度就会受到影响，他们在整个过程中采用了一些假设：单元的板壳划分；忽略了坯料沿轧制方向的运动，使整个过程改变成成型辊将坯料一步一步压下，直到成型为止；忽略了坯料内的剪应变；将坯料作为线性强化的弹塑性体；忽略了中性面的移动：忽略了接触区以外的变形；这样不能真实反映物体的成型过程，仅能够近似模拟整个变形过程。由于直缝焊管成型过程理论上属于弹塑性

27、大变形和接触非线性多重非线性耦合问题，虽然弹塑性大变形和接触非线性有限元算法和软件已经基本完善，但如果精确地刻画轧辊和带钢，模拟整个轧制过程会使单元数量和计算量急剧增加。焊管成型属于多道轧制，如果分道模拟轧制过程，虽然可以降低计算规模，但是因为前道和后道的变形不仅相互依赖相互影响，还伴随弹性释放，因此分道模拟和连续模拟结果还有区别。总之，在现有的计算机软硬件条件下，对焊管成型过程的模拟不仅计算量大，计算难度也很大，还需要进行深入的研究。有限元法的发展许多工程分析问题，如固体力学中的位移场和应力场分析、电磁学中的电磁场分析、振动特性分析、传热学中的温度场分析、流体力学中的流场分析等，都可以归结为

28、在给定边界条件下求解其控制方程（常微分方程或者偏微分方程）的问题。但是，能够用解析方法求出精确解的只是方程性质比较简单，且几何边界相当规则的少数问题。对于大多数的工程技术问题，由于研究对象的几何形状复杂或者它的某些特性是非线性的，则很少有精确解。这类问题的解决通常有两种途径：其一，引入简化假设，将方程和边界简化为能够处理的问题，从而得到它在简化状态下的解。这种方法只是在非常有限的情况下可行，因为过多的简化将可能导致不正确的甚至错误的解。其二，数值模拟。目前，在工程技术领域内常用的数值模拟方法有：有限单元法、边界元法、离散单元法和有限差分法。但就其实用性和应用的广泛性而言，主要还是有限单元法。有

29、限元方法是用于求解工程中各类问题的数值方法。应力分析中稳态的、瞬态的、线性或非线性的问题以及热传导、流体流动和电磁学中的问题都可以用有限元方法进行分析解决。有限单元法的基本思想是：将问题的求解域划分为一系列单元，单元之间仅靠节点联结。单元内部点的待求量可由单元节点量通过选定的函数关系插值得到。由于单元形状简单，易于由平衡关系或者能量关系建立节点之间的方程式，然后将各个单元方程“组集”在一起形成总体代数方程组，加入边界条件后即可对方程求解。单元划分越细，计算结果就越精确。现代有限元方法的起源可以追溯世纪的早期，当时一些研究者应用离散的等价第一章绪论杆拟合模态的弹性体。然而，人们公认（）是应用有限

30、元方法的第一人。使用分段多边形插值法而不是三角剖分法来研究扭转问题。世纪年代，公司采用三角元对机翼进行建模，大大推动了有限元方法的应用。其他不少人接着采用了这一方法。然而，直到世纪年代，人们才广为接受“有限元”这一术语，【在这方面做了不少工作。在世纪年代，研究者开始将有限元方法应用到解决工程中的其他领域，例如热传导和地下渗流中的问题。和（）在年撰写了第一本有限元的专著。年来，有限单元法的应用已经由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳态问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、传热学等连续介质力学领

31、域。在工程分析中的作用已经从分析和校核扩展到优化设计并和计算机辅助设计技术相结合。可以预计，随着现代力学、计算数学和计算机技术等学科的发展，有限单元法作为一个具有巩固理论基础和广泛应用效力的数值分析工具，必将在国民经济建设和科学技术发展中发挥更大作用。仿真工具介绍大型工程结构庞大，计算模型复杂，因而，它的分析软件一般均由前后处理器和分析求解器两个软件集成，并在两者之间常设有功能强大的软件接口，目的是保持这两部分的相互独立性并方便于与其它软件的连接。图中给出了大型有限元分析软件的整体结构图，同时也反映了大型有限元分析软件的操作过程和主要功能构成。其中：前处理器的主要功能是：形成并提交求解器能识别

32、的有限元模型数据文件，包括数据文件的输入输出功能、自动网格划分功能、生成几何图形功能和图形显示及输出功能。因此，前处理器要处理如下数据文件：读写描述有限元模型的数据库文件；生成求解器可识别的输入数据文件；读写并记录各种操作的命令；编辑过程文件等。求解器的主要功能是：对有限元模型形成的整体方程组进行计算和分析，包括带宽优化、方程组求解方法的选取、确定收敛判据等等。因此，它是有限元的理论方法、计算方法和软件分析方法的结晶。后处理的主要功能是：处理求解器分析后的结果，包括识别求解器生成的格式化及非格式化结果文件；将结果写成各种常用的图形文件；还要以图形、动画、曲线、表格和文件形式分别显示求解器分析后

33、的结果等。软件接口的主要功能是：读写可被其它的软件识别的各种文件，实现与其它软件间的数据转换功能，使得几何建模可以由第三方软件来提供，如此增强了软件本身的健壮性。目前我国已引进的主要有限元分析软件有：、以及等等。本文主要使用软件进行仿真分析，下面将加以介绍说明。第一章绪论前处理器求解器后处理几何建模自动划分网形成有限格兀模型一添加载荷和约束设置邕兀特性确定材料特性形成数据文件递交数据文件分析与求解形成结果文件读取结果文件结果后处理软件接口软件接口软件接口图软件整体结构框图软件简介公司由博士创立于年，有限元程序是该公司主要产品，它是第一个通过质量认证的分析设计类软件。由于软件具有建模简单、快速、

34、方便的特点，因而成为大型通用有限元程序的代表。软件是集结构、热、流体、电磁、声学于一体的大型通用有限元分析软件，可广泛地应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电等一般工业及科学研究。软件的主要功能包括建立模型、结构分析、非线性分析、电磁分析、计算流体力学分析、接触分析、压电分析、结构优化等。的建模功能有限元分析软件的前处理器功能很强大，具有强大的建模功能。建模时，需要先建立结构的几何模型，给出材料参数和单元类型，最后划分网格，形成结构的有限元模型。提供了三种创建模型的方法：实体建模方法、直接建模方法和

35、输入在其它计算机辅助设计系统中创建的模型。直接建模的方法就是在的前处理程序中直接定义每个节点的位置以及单元的大小、形状和连通性来创建有限元模型。节点用来定义单元在空间的位置，单元定义了模型的连接性。直接建模的方法适用于线模型和较简单的有规则几何结构。可以自己控制每个单元和节点的编号。但是直接建模的方法往往需要处理大量的数据，也不能进行自适应网格划分，改进网格划分十分困难。对大而复杂的结构，应采用实体建模的方法。第一章绪论提供了两种方法进行实体建模，即自底向上的建模方法和自上向下的建模方法。自底向上的建模方法是先创建关键点，然后依次创建相关的线、面和体等图元。自上而下的建模方法是可以直接创建最高

36、级的图元，如球、棱柱等三维实体，通常称之为几何体素。当定义了一个体素时，程序会自动定义相关的面、线和关键点。可以利用这些高级图元直接构造几何造型。在建模过程中，自上而下的建模方式和自底向上的建模方式可以自由组合使用，使模型的创建更加方便。实体建模的优点是：对于庞大或复杂的模型，特别是对三维实体模型更合适；相对而言需处理的数据少一些；容许对节点和单元不能进行的几何操作（如拖拉和旋转）；支持使用面和体的体素及布尔运算以顺利建立模型；便于施加载荷之后能进行所要求的局部网格细化；便于几何上的改进；便于改变单元类型，不受分析模型的限制。无论是使用自底向上还是自上而下的方法构造实体模型，均由关键点、线、面

37、和体组成。这些图元的层次关系是：体以面为边界，面以线为边界，线以关键点为端点。无论是使用自上而下还是自底向上方法建模，均能使用布尔运算来组合数据集，形成一个实体模型。程序提供了完整的布尔运算，诸如相加、相减、相交、分割、粘结和重叠。在创建复杂实体模型时，对线、面、体的布尔运算操作能减少相当可观的建模工作量。程序还提供了拖拉、延伸、旋转、移动、拷贝、蒙皮等的功能，可以大大减少建模时间。辅助工具如选择和组元、拾取与利用工作平面，为建模提供了极大的方便。可以在其它系统中建立模型并把它输入到中进行分析。它有如下优点：一是可利用已有模型，避免重复工作；二是可利用熟悉的工具去建模。但是，从系统中输入的模型

38、，如果不适于网格划分则需要大量的修补工作。的分析功能软件的分析功能包括结构分析、非线性分析、热分析、电磁场分析、电场分析、流体流动分析、耦合场分析。结构分析包括静力学分析、模态分析、谐波分析、瞬态分析、响应谱分析、随机振动分析和屈曲分析。非线性分析包括结构非线性、材料非线性、几何非线性分析和单元非线性分析。热分析包括稳态热分析、瞬态热分析、相变分析和热结构耦合分析。电磁场分析包括静态电磁场、低频时变电磁场和高频时变电磁场。电场分析包括电流传导、静电分析和电路分析。流体流动分析包括计算流体力学分析层流、湍流和热流体）和声学分析。耦合场分析适合于热应力、磁热、磁结构、流体流动热、流体流动结构、热电

39、、电路耦合电磁场、电磁及压电耦合分析。软件还具有将部分单元等效为一个独立单元的子结构功能、将模型中的某一部分与其余部分分开重新细划网格的子模型功能和参数设计语言（）可以进行结构的优化。的接触分析功能是优于其它有限元软件的，近年来，开发了一系列的接触单元，接触单元通常用来描述两物体相互接触或滑动的界面。刚开始中有节点对节点单元和，接着有节点对地单元，然后有节点对面单元年。虽然能求解一些接触问题，但有相当的局限性。最近几年，中引入一类面对面接触单元：目标单元、和接触单元、，达到真正的面面接触，才具有真正意义上的接触分析功第一章绪论能。接触问题是一种高度非线性行为，需要较大的计算资源，为了进行有效的

40、计算，理解问题的特性和建立合理的模型是很重要的。接触问题存在两个较大的难点：其一，在求解问题之前，不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定：其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供挑选，它们都是非线性的，摩擦使问题的收敛性变得更困难。接触问题分为两种基本类型：刚体柔体的接触，柔体柔体的接触，在刚体柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度，一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体柔体的接触，许多金属成形问题归为此类接触。另一类，柔体柔体的接触，是一种更普遍的类型

41、，在这种情况下，两个接触体都是变形体（有近似的刚度）。支持三种接触方式：点点、点面、面面，每种接触方式使用的接触单元适用于不同的问题。为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对应组元是一个节点。如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触对，接触单元是覆盖在分析模型接触面之上的一层单元。点点接触单元主要是用于模拟点点的接触行为，为了使用点点的接触单元，需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下）

42、。点面接触单元主要用于给点面的接触行为建模，例如两根梁的相互接触。如果通过一组节点来定义接触面，生成多个单元，那么可以通过点面的接触单元来模拟面面的接触问题，面既可以是刚性体也可以是柔性体，这类接触问题的一个典型例子是插头插到插座里。使用这类接触单元，不需要预先知道确切的接触位置，接触面之间也不需要保持一致的网格，并且允许有大的变形和大的相对滑动。支持刚体柔体的面面的接触单元，刚性面被当作“目标面”，分别用和来模拟和的“目标面”，柔性体的表面被当作“接触面”，用，来模拟。一个目标单元和一个接触单元叫做一个“接触对”。与点面接触单元相比，面面接触单元有以下几项优点：支持低阶和高阶单元；支持有大滑

43、动和摩擦的大变形，协调刚度阵计算，单元不对称刚度阵的选项：提供工程中采用的更好的接触结果；没有刚体表面形状的限制，刚体表面的光滑性不是必须的，允许有自然的或网格离散引起的表面不连续。第一章绪论本文的主要研究内容及创新点本文的主要研究内容本文的土要研究内容包括：（）采用弹塑性大变形和接触非线性有限元法研究了牢拔管，芯棒拔管的成型过程，考虑了模具的变形。模拟过程中参数研究包括了模具的形状、尺寸、材料模型和强化规律。（）发明了一种冷拔钢管高精度预应力复合模具。（）研究了冷拔钢管工程应用中的相关问题。（）建立了较为符合实际的焊管成型过程的三维有限元模型，包括轧辊（水平辊和立辊）和管坯，其中建立轧辊的三

44、维模型是为了研究轧辊对管坯的真实作用，建立管坯的三维模型是为了研究有压靠力作用时钢管的力学行为以及铡管脱离轧辊后的回弹。（）采用弹塑性大变形和接触非线性有限元法模拟了焊管第一架和第二架的成型过程，获取了钢管和轧辊在成型过程中的位移、应力和应变场的分布。成型方式考虑了单半径和双半径两种成型方式。本文研究的主要创新点本文研究的主要创新点包括：（）采用三维弹塑性大变形和接触非线性有限元法研究了空拔管、芯棒拔管的成型过程，考虑了模具的变形，初步完成了冷拔钢管和焊管的专用计算软件，方便了参数研究，提高了计算效率，为工程师使用大型分析软件解决了模具与工艺设计问题奠定了基础。（）对工艺参数、材料模型和强化规

45、律对成型过程的影响进行了深入的研究，发现了锥模和弧型模对直径和壁厚精度的不同影响，发展了一种高精度预应力复合模具，解决了采用苏式模和中式模钢管尺寸精度低以及模具磨损快和寿命低的问题。（）采用三维有限元法模拟了单半径和双半径焊管的成型过程，考虑了压靠力对成型过程的影响，发现了单半径和双半径成型方法的不足，为合理设计焊管成型工艺和焊管轧辊奠定了基础。（）发展了一种测量钢管内表面残余应力的杂交实验方法，测量了钢管内表面的残余应力，研究了残余应力对液压缸缸筒承载能力的影响。残余应力，研究了残余应力对液压缸缸筒承载能力的影响。第二章空拔钢管成型过程的有限元分析第二章空拔钢管成型过程的有限元分析本章中的数

46、值模拟工作都是采用大型有限元仿真软件完成的。在本章中主要包括以下几部分内容：数值模拟方法的验证。将已有的实验数据同通过数值模拟方法得到的数据进行对照，以验证数值模拟方法的正确性与可行性。对现实生产中采用不同孔形模具进行钢管的空拔成型过程进行数值模拟。大口径钢管成型过程的数值模拟通过数值模拟分析不同材料模型和不同强化规律对钢管拔制力和成型精度的影响。通过对已有的各种模拟结果的分析，提出新型模具的设计方法。下面就这五方面内容进行详细阐述。数值模拟方法的验证为了说明数值模拟方法的正确性与可行性，在研究中首先对数值模拟方法的正确性与可行性进行了验证。验证的实验数据来源于钢管技术，年第期，钟锡汉等人的无芯棒拔管时管子的壁厚变化一文中提供的实验数据。通过将实验数据与数值模拟得到的数据进行对照发现，数值模拟得到的数据与实验数据在误差范围内基本吻合，从而验证了数值模拟方法的正确性与可行性。表中列出了验证的两组实验数据，下面就第二组实验数据的验证过程介绍如下：模具与单元划分模具和钢管尺寸如图所示。钢管拔制前尺寸为，成品管直径为，模孔直径取为成品管径，模具的外径为，长度为，定径带尺寸为。取长为的一段管材，考虑到对称性，本文取模型，管材与模具均采用节点三维实体单元，单元数量总计为个，其中钢管个单元，模具媳个单元，模具内表面划分目标单元个，钢管外表面划分接触单元个，管材与模具的有限元网格如图所示，图中