钢管穿轧过程的三维热力耦合模拟

描述了使用有限元方法(FEM)模拟钢管在曼内斯曼式穿轧机上的穿孔过程,建立了三维有限元模型,并考虑了在金属成形过程中出现的热力学现象。模拟的结果动态地显示了在工件内部等效应变、等效应变率和温度的分布。基于这些参数,分析了钢管的穿孔过程,同时通过实验测试与有限元模拟的结果进行对比。结果表明,测试的结果与有限元模拟的结果比较吻合。     

基于响应面法橡皮成形侧压块参数优化及验证

橡皮成形是飞机钣金零件制造的一种重要的成形工艺,为提高典型零件橡皮成形的效率,针对橡皮成形过程中产生的起皱问题,基于响应面法对侧压块参数进行优化.首先对新淬火状态下的铝合金板料进行材料试验,分析了其成形性能指标参数,为有限元模拟提供材料本构关系与模型;然后以起皱量化为优化目标,将凸弯边的最大增厚率作为侧压块优化的目标函数,为得到典型零件橡皮成形侧压块精确的几何参数,采用Abaqus有限元软件对带侧压块成形过程进行正交模拟试验,得出侧压块几何参数对起皱指标影响程度的主次关系,并得到最优的参数水平组合.基于响应面建模方法建立了侧压块最优几何参数和起皱指标之间的多项式关系,并用遗传算法进行参数优化,利用优化后的工艺参数进行有限元模拟,阐述了侧压块防皱原理;最后通过实际成形试验,验证了有限元模拟的可行性和优化方法的可靠性.     

Localized remeshing techniques for three‐dimensional metal forming simulations with linear tetrahedral elements

The localized remeshing technique for three‐dimensional metal forming simulations is proposed based on a mixed finite element formulation with linear tetrahedral elements in the present study. The numerical algorithm to generate linear tetrahedral elements is developed for finite element analyses using the advancing front technique with local optimization method which keeps the advancing fronts smooth. The surface mesh generation using mesh manipulations of the boundary elements of the old mesh system was made to improve mesh quality of the boundary surface elements, resulting in reduction of volume change in forming simulations. The mesh quality generated was compared with that obtained from the commercial CAD package for the complex geometry like lumbar. The simulation results of backward extrusion and bevel gear and spider forgings indicate that the currently developed simulation technique with the localized remeshing can be used effectively to simulate the three‐dimensional forming processes with a reduced computation time. Copyright © 2006 John Wiley & Sons, Ltd.     

Role of Friction in Cold Ring Rolling

1. Introduction Cold ring rolling is a main technology used to manu- facture various precise seamless ring shape parts. It has been increasingly used in many industrial fields such as bearing, machine, automobile, petrochemicals, aeronau- tics, astronautics and atomic energy because of its many technical superiorities such as considerable saving in en- ergy and material cost, high quality, high efficiency, and low noise, etc. To research and develop advanced precise cold ring rolling technolog…     

Improvement of a frictional contact algorithm for strongly curved contact problems

One of the challenges in contact problems is the prediction of the actual contact surface and the kind of contact that is established in each region. In numerical simulation of deep drawing problems the contact conditions change continuously during the forming process, increasing the importance of a correct evaluation of these parameters at each load step. In this work a new contact search algorithm devoted to contact between a deformable and a rigid body is presented. The rigid body is modelled by parametric Bézier surfaces, whereas the deformable body is discretized with finite elements. The numerical schemes followed rely on a frictional contact algorithm that operates directly on the parametric Bézier surfaces. The algorithm is implemented in the deep drawing implicit finite element code DD3IMP. This code uses a mechanical model that takes into account the large elastoplastic strains and rotations. The Coulomb classical law models the frictional contact problem, which is treated with an augmented Lagrangian approach. A fully implicit algorithm of Newton–Raphson type is used to solve within a single iterative loop the non‐linearities related with the frictional contact problem and the elastoplastic behaviour of the deformable body. The numerical simulations presented demonstrate the performance of the contact search algorithm in an example with complex tools geometry. Copyright © 2003 John Wiley & Sons, Ltd.     

A new boundary element approach for contact problems with friction

A new boundary element solution algorithm for two-dimensional and axisymmetric contact problems with friction, based on an independent discretization of the contacting surfaces and under static and proportional loading conditions, is presented. The solution procedure uses the element shape functions to distribute the geometry, tractions and displacements on each contact element. The contact constraints are then applied between each contacting node and the opposite contact segment. The overall boundary element matrix equations for the contacting bodies are coupled using the contact conditions at the interface without introducing any additional variables into the solution matrix. The algorithm is applied to several two-dimensional and axisymmetric frictional contact examples and the results obtained are in very good agreement with finite element and analytical solutions.     

Shape optimization of the initial blank in the sheet metal forming process using equivalent static loads

In the sheet metal forming process, forming the final desired shape is difficult to obtain due to wrinkling, tearing, failure of material, etc. Various conditions of the forming process should be controlled for the desired shape. These conditions are the velocity of the punch, the friction factor, the blank holding force, the initial shape of the blank and others. Many researchers have conducted studies to predetermine the initial blank shape. The structural optimization technique is one of them. Non‐linear response structural optimization is required because non‐linearities are involved in the analysis of the metal forming process. When the conventional method is utilized, the cost is extremely high due to repeated non‐linear analysis for function and sensitivity calculation. In this paper, the equivalent static loads (ESLs) method is used to determine the blank shape which leads to the final desired shape and reduced wrinkling. The ESLs method is a structural optimization method where non‐linear dynamic loads are transformed into ESLs, and these ESLs are utilized as external loads in linear response optimization. The design is updated in linear response optimization. Non‐linear analysis is performed with the updated design and the process proceeds in a cyclic manner. An optimization formulation is defined for the examples, the formulated problems are solved to verify the proposed method and the results are discussed. Non‐linear analysis is performed using the commercial software LS‐DYNA, NASTRAN is used for calculating the ESLs and linear response optimization, and an interface program for LS‐DYNA and NASTRAN is developed. Copyright © 2010 John Wiley & Sons, Ltd.     

Improvement of stability conditions, accuracy and uniqueness of penalty approach in contact modeling

The main objective of this paper is to improve stability conditions, uniqueness and convergence of numerical analysis of metal forming processes with contact constraints enforced by the penalty method. A commonly known drawback of this approach is the choice of penalty factor values. When assumed too low, they result in inaccurate fulfillment of the constraints while when assumed too high, they lead to ill-conditioning of the equations system which affects stability and uniqueness of the solution. The proposed modification of the penalty algorithm consists in adaptive estimation of the penalty factor values for the particular system of finite element equations and for the assumed allowed inaccuracy in fulfillment of the contact constraints. The algorithm is tested on realistic examples of sheet metal forming. The finite element code based on flow approach formulation (for rigid-plastic and rigid-viscoplastic material model) has been used.     

板料V形弯曲回弹的动力烛式有限元分析

板料成形后的回弹对精度影响较大，在数值模拟时对回弹进行精确预测显得非常重要。基于连续介质力学及有限变形理论，建立了适合于三给板料成形分析的显式算法的有限元数学模型，采取集中质量矩阵，用动力显式积分的方法，使位移计算显式化，避免了由材料、几何、边界条件等高度非线性因素引起的计算收敛问题。根据该模型开发了动力显式算法的板料成形过程模拟的有限元分析程序DESSFORM3D，应用该软件模拟了包括回弹在内的整个板料V形弯曲的成形过程。通过3个不同凸模行程时计算与实验的板料几何形状对比以及计算结果与实验结果对比，验证了软件计算结果的准确性。     

Coupled thermo-mechanical analysis of metal-forming processes through a combined finite element–boundary element approach

This paper presents a new approach for obtaining the distribution of temperature in the dies during thermo-mechanical numerical analysis of metal forming problems. The proposed approach is based on a solution resulting from the combination of the finite element method with the boundary element method. The finite element method is used to perform the numerical modelling of the thermo-mechanical deformation of the workpiece, taking into account the geometrical and material non-linearities as well as the influence of the temperature distribution on the mechanical behaviour of the material. The boundary element method is applied for computing the distribution of temperatures in the dies. The combination of the two numerical methods is made using the finite element solution of the heat flow exchanged across the die–workpiece interface to define the boundary conditions to be applied on the thermal analysis of the dies. A numerical example of compression under plane-strain conditions is included to show the applicability of the proposed approach. © 1998 John Wiley & Sons, Ltd.     

高强度钢板U形件热冲压凹模结构拓扑优化

为降低热冲压凹模的生产成本和使用成本,基于板料热冲压数值模拟对凹模结构进行了拓扑优化设计.运用有限元软件ABAQUS建立热力耦合有限元模型,对高强钢板U形件的热冲压成形和淬火过程进行了数值仿真.提取凹模与板料间关键工况下接触应力作为凹模拓扑优化的外在载荷,建立约束凹模结构关键区域节点位移的体积最小化拓扑模型,对热冲压凹模结构进行拓扑优化设计,最终实现结构减重20%,且优化后凹模的变形和应力与优化前的结果相差甚微.研究内容对热冲压过程数值模拟和模具结构拓扑优化研究具有一定参考价值.     

Adaptive 2D finite element simulation of metal forming processes

The development and integration of available current methods and the development of new methods for an adaptive finite element analysis of metal forming processes are presented. The analysis includes large-strain, elastic–plastic, and thermal effects. Many numerical methods such as mesh generation, simulation of the contact between the workpiece and tool and die, and optimization of the finite element mesh are integrated and incorporated. In addition, an algorithm is developed which can detect certain severely distorted elements where the area of integration is approaching zero. The advantage of correcting these regions of locally distorted elements is demonstrated. These numerical methods are implemented in a finite element program developed for simulating metal forming processes, with the emphasis on automating the analysis. Examples include an axisymmetric stress simulation of a coldheading process, a plane strain simulation of an extrusion process and a plane strain simulation of orthogonal metal cutting, all with noticeable thermal effects. The orthogonal cutting forces and feed forces calculated are compared with two sets of experimental data, with good agreement.     

板料冲压成形CAD数字化分析基本理论

板料成形过程是涉及几何非线性、材料非线性和边界条件非线性等的多重非线性问题.介绍了板料冲压成形CAD数字化分析所涉及到的基本理论--非线性弹塑性材料的本构关系、板壳成形单元模型、接触问题、有限元方程及求解.     

基于神经网络遗传算法的深腔型零件拉深工艺参数优化

目的 解决冲压成形中工艺参数优化难的问题.方法 以一种深腔型零件的冲压成形为例.首先,借助灰度关联分析法对有限元中的工艺参数进行分析,获取该零件冲压成形中影响成形质量的2个主要因素——冲压速度和压边力.其次,借助拉丁超立方抽样法对上述2个因素进行随机取样,并借助DYNAFORM软件对其进行逐一模拟.再次,将冲压速度和压边力作为输入,最大减薄作为输出,训练在MATLAB中建立的BP神经网络,并借助遗传算法对其进行寻优.结果 最优成形压边力为1.372 MN,最优加载速度为1.5366 m/s.结论 与神经网络遗传算法预测相比,有限元结果的相对误差小于2％,零件试制结果的相对误差小于6％,该方法有较高的预测精度.     

铝合金汽车顶盖充液成形的数值模拟

目的研究铝合金汽车顶盖拉延工序的充液成形工艺。方法基于有限元分析软件Dynaform,利用带局部刚性凹模整形的被动式充液成形工艺,通过建立有限元分析模型,优化成形过程中的关键工艺参数,分析变形规律并进行质量控制。结果成形过程中的液室压力加载路径、压边力、拉延筋,以及坯料形状等工艺参数对成形影响较大。液室压力不宜过早加载。液室压力过大或压边力过小不利于顶部产生充分塑性变形。压边力过大极易造成顶盖圆角处的破裂。结论该成形工艺可行,且数值模拟的准确性及适用性较高,采用该成形工艺可得到表面质量良好,未出现起皱、破裂缺陷的合格零件。     

A simple finite element‐based framework for the analysis of elastic pseudo‐rigid bodies

This article advocates a general procedure for the numerical investigation of pseudo‐rigid bodies. The equations of motion for pseudo‐rigid bodies are shown to be mathematically equivalent to those corresponding to certain constant‐strain finite element approximations for general deformable continua. A straightforward algorithmic implementation is achieved in a classical finite element framework. Also, a penalty formulation is suggested for modelling contact between pseudo‐rigid bodies. Representative planar simulations using a non‐linear elastic model demonstrate the predictive capacity of the pseudo‐rigid theory, as well as the robustness of the proposed computational procedure. Copyright © 1999 John Wiley & Sons, Ltd.     

A combined finite‐discrete element method for simulating pharmaceutical powder tableting

The pharmaceutical powder and tableting process is simulated using a combined finite‐discrete element method and contact dynamics for irregular‐shaped particles. The particle‐scale formulation and two‐stage contact detection algorithm which has been developed for the proposed method enhances the overall calculation efficiency for particle interaction characteristics. The irregular particle shapes and random sizes are represented as a pseudo‐particle assembly having a scaled up geometry but based on the variations of real powder particles. Our simulations show that particle size, shapes and material properties have a significant influence on the behaviour of compaction and deformation. Copyright © 2005 John Wiley & Sons, Ltd.     

运输支板成形仿真及模具设计

目的运输支板的曲线边较长,在钣金翻边过程中易开裂,为提高模具设计的合理性,避免重复试模、返修以及设计失败的风险,开展有限元仿真研究。方法使用仿真软件Pam-Stamp,并与模具实物、产品试验件相结合,确定单动压边、双动压边翻边两次成形工艺,并优化参数、设计新模具。结果新设计的模具经过产品试验件生产,无开裂问题,目前已用于正式产品生产。结论仿真结果的成形趋势与实际生产基本符合。在选择工艺、工装设计方案时,仿真结果具有较高的参考价值。     

球面端坯料热滚切成形仿真及机理分析

坯料端部预加工成球面状是改善楔横轧件凹心缺陷的有效方法.为了避免传统机加工制坯带来的材料损失,本文提出了一种基于金属塑性成形的球面端部热滚切成形方法.建立了球面端坯料热滚切成形的有限元模型,通过仿真分析了成形过程中滚切区金属的位移场变化规律和应变场分布特征,阐述了球面端坯料的成形机理以及在挡块作用下堆料的抑制机理,并对...