变形—损伤耦合的超塑性恒压轴对称充模胀形的有限元模拟

采用大变形刚粘塑性有限元法模拟超塑性恒压轴对称充模胀形过程，分析了模具几何参数及材料参数对胀形过程中材料的流变行为、胀形制件厚度分布和成形时间的影响规律，给出了质点的流动轨迹，不同时刻制件的剖面形状及应力、应变分布；基于修正的Ｇｕｒｓｏｎ粘塑性势推导了内部空洞体积分数累积增大模型并据此进行了变形－损伤耦合计算。     

轴对称件超塑约束胀形中的摩擦效应

尽管超塑胀形作为一种金属成形方法正日渐受到世界各国的普遍关注，但对超塑约束胀形理论的研究报道却还很少，尤其对胀形过程的有限元模拟研究就更为罕见。针对这种情况，采用大变形刚粘塑性有限元法模拟了轴对称零件向圆筒形凹模内超塑约束胀形的变形过程，着重研究了工具工件之界面摩擦对胀形件厚度分布不均匀性和胀形板料向凹模角部充填性的影响．结果表明，随着摩擦的降低，胀形件的侧向较厚部分能有所减薄，可以改善整个胀形件的厚度均匀性，但当摩擦因子Ａｍ≤０．２（相当于摩擦系数μ≤０．１２）时，胀形件极顶部分的减薄过大；摩擦较小时，胀形板料向凹模角部的充填性较好；在考虑到极点附近厚度适度减薄和胀形板料对凹模角部充填性好的前提下，工艺上应当适当减小摩擦，其最佳状态是μ值约为０．３．为了检验所用刚粘塑性有限元法模拟的可靠性，将计算结果与试验结果作了对比，发现两者相当吻合。     

材料参数对板材胀形过程综合影响的数值研究

本文将Hosford 与Hill 各向异性屈服函数应用于刚粘塑性有限元方法,分析了圆形薄板液压胀形过程.研究了材料性能参数:硬化指数n、速率敏感指数m、厚向各向异性参数R、屈服函数非多项式指数M 对液压胀形过程的综合影响.并通过数值分析,找出了临界厚向断裂应变-ε_3~(ov·)与材料参数关系的经验方程式.屈服表面形状对极限厚度应变的影响,可以用Barlat 提出的包含了R 和M 影响的p 值表示出来.     

弹性/粘塑性线性硬化模型及有限元算法

1.粘塑性力学的本构方程粘塑性力学的近日研究说明它不仅可用于解决那种对变形率敏感的材料的力学问题,例如材料在高速载荷作用下的时滞现象问题,而且还适用于解决一般塑性力学、应力松弛、蠕变等问题.按照文的定义:弹性-粘塑性,材料在其弹性及塑性阶段均出现粘滞现象;弹性/粘塑性,材料仅在其塑     

超塑性恒压充模胀形过程各力学量场的影响因素及影响规律

采用刚粘塑性有限元法模拟了超塑性恒压充模胀形过程，分析了应变速率敏感性指数ｍ、厚向异性系数Ｒ、初始板厚、界面摩擦等因素对应力、应变和应变速率等力学量场的影响；通过对铝合金ＬＹ１２ＣＺ板料的相应实验对模拟结果进行了验证     

双相钢板料的单向拉伸断裂失效研究(Ⅱ)——弧长法非线性有限元分析

为了获得docol800DP双相钢板料的等效塑性应变和应力三轴度在单向拉伸断裂失效过程中的变化历程,对其拉伸过程进行了弧长法非线性有限元分析。根据试验得到的材料参数建立了双相钢单向拉伸试件的有限元模型,考虑大变形引起的几何非线性和塑性强化引起的材料非线性,模拟了试件拉伸变形的全过程。计算结果表明:弧长法可较好地模拟试件变形的全程,计算得到的典型变形阶段、集中性失稳带的分布及方向与试验结果吻合很好;发生失稳变形时的应变与理论分析结果一致。最后,给出了起裂点处失效参数的变化历程,为定量化研究双相钢材料的断裂失效模型提供准确的数据支持。     

摆辗成形后的摩托车启动齿轮三维塑性应变场的光塑性分析

本文叙述了采用聚碳酸脂材料模拟低碳钢和低碳合金钢一类金属的摩托车启动齿轮(下面陶称为端面齿轮)摆辗成形过程。运用光塑性法计算分析了氐形后的光塑性模型内部三维塑性应雯分布;对采用螺旋形轨迹和玫瑰线轨迹成形后的光塑性模型内部的三维塑性应变分布进行了分析比较。该方法已为摆辗最佳工艺参数的选择提供了有用的依据。     

材料的弹粘塑性本构关系和应力波的演化

试验表明，大多数工程材料在冲击载荷作用之下的变形一般都同时包含有可恢复的瞬态性弹性变形和不可恢复的粘滞性塑性变形，即其本构关系可以用弹粘塑性模型来描述。本文从内变量理论出发，探讨了时率相关材料的弹粘塑性本构关系的一般特性，建立了增量型的弹粘塑性本构关系的一般理论框架和普适的表达式，并且对两种最常用的本构模型——Bodner-Partom模型和Johnson-Cook模型给出了在一维应变条件下的具体形式。通过计算和讨论一维应变粘塑性靶板中冲击波的衰减机制和应力波的演化规律，特别是考察各种粘塑性本构模型中的材料参数对冲击波的衰减和应力波的演化的影响，得出了一些可以直接应用或具有一定借鉴价值的结果，为研究应力波的其他衰减机制以及在人防工程中智能防护层设计时新材料的选取奠定了基础。     

等径通道挤压过程三维有限元模拟

利用三维有限元模型对等径通道挤压过程中变形分布的不均匀性进行了模拟,通过对不同摩擦系数下截面等效塑性应变分布比较发现:沿三个方向截面上的变形分布都不是均匀的,这说明二维有限元模型不能真实模拟等径通道挤压过程中试样中的变形分布。此外,摩擦对等效塑性应变分布及挤压力的影响较大,截面变形不均匀参数和最大等效塑性应变出现的位置随摩擦系数的增大而变化。压力-位移曲线的变化可以结合试样的变形过程来解释。     

粗晶超塑性单轴拉伸颈缩模拟

利用粘塑性本构模型模拟粗晶超塑性单轴拉伸。数值结果表明，颈缩的位置及发展过程受拉伸应变速率的影响。不只一处分散不均匀变形相互牵制与协调，使材料得以在接近均匀的状态下经受大的变形，模拟得到的局部应变速率演化曲线，可以预测变形局部化发展的情况。     

Experimental mechanics at velocity extremes —Very high strain rates

The early experimental work of Clark and Wood with regard to von Kármán's theory on the effect of material flow and fracture at high strain rates has led to many controversial issues on these effects. Interest has been greatly revived in recent years because of the increased emphasis on such high-velocity forming processes as explosive and capacitor discharge. Considerable new work has been performed by Ling-Temco-Vought, Inc., for the Air Force, the results of which are presented in this paper. Data have been accumulated on tensile and compression specimens, spherical bulging and cylindrical bulging for a wide variety of materials. This high-speed information has been generated with the use of a special projectile impact machine and special presses utilizing various combinations of explosive and capacitor-discharge energy, with strain rates to 10¹/sec. The effect of velocity on ductility is discussed for total strain distribution, uniform strain, double necking and critical impact velocity. The modes of failure for various part shapes are presented and related to the forming velocity.     

板材多点成形过程的有限元分析

多点成形过程采用静力隐式格式进行数值模拟是比较合适的。本文建立了用于多点成形过程分析的静力隐式弹塑性大变形有限元方法 ,给出了对稳定迭代收敛过程效果较好的板壳有限单元模型、处理多点不连续接触边界的接触单元方法以及增量变形过程中应力及塑性应变计算的多步回映计算方法。基于这些方法编制了计算软件 ,应用该软件进行了矩形板的液压胀形过程及球形模具拉伸成形过程的有限元分析 ,数值计算结果与典型的实验结果及计算结果吻合很好。最后给出了球形、圆柱形目标形状的实际多点成形过程的数值模拟结果。     

A multi-length scale sensitivity analysis for the control of texture-dependent properties in deformation processing

Material property evolution during processing is governed by the evolution of the underlying microstructure. We present an efficient technique for tailoring texture development and thus, optimizing properties in forming processes involving polycrystalline materials. The deformation process simulator allows simulation of texture formation using a continuum representation of the orientation distribution function. An efficient multi-scale sensitivity analysis technique is then introduced that allows computation of the sensitivity of microstructure field variables such as slip resistances and texture with respect to perturbations in macro-scale forming parameters such as forging rates, die shapes and preform shapes. These sensitivities are used within a gradient-based optimization framework for computational design of material property distribution during metal forming processes. Effectiveness of the developed computational scheme is demonstrated through computationally intensive examples that address control of properties such as Young’s modulus, strength and magnetic hysteresis loss in finished products.     

SECTIONAL FINITE ELEMENT ANALYSIS OF COUPLED DEFORMATION BETWEEN ELASTOPLASTIC SHEET METAL AND VISCO-ELASTOPLASTIC BODY

The present paper is devoted to developing a new numerical simulation method for the analysis of viscous pressure forming (VPF), which is a sheet flexible-die forming (FDF) process. The pressure-carrying medium used in VPF is one kind of semisolid, flowable and viscous material and its deformation behavior can be described by the visco-elastoplastic constitutive model. A sectional finite element model for the coupled deformation analysis between the visco- elastoplastic pressure-carrying medium and the elastoplastic sheet metal is proposed. The resolution of the Updated Lagrangian (UL) formulation is based on a static explicit approach. The frictional contact between sheet metal and visco-elastoplastic pressure-carrying medium is treated by the penalty function method. Coupled deformation between sheet metal and visco-elastoplastic pressure-carrying medium with large slip is analyzed to validate the developed algorithm. Finally, the viscous pressure bulging (VPB) process of DC06 sheet metal is simulated. Good agreement between numerical simulation results and experimental measurements shows the validity of the developed algorithm.     

Analysis of tube hydroforming in a square cross-sectional die

A mathematical model considering sliding friction between the tube and die is proposed to explore the plastic deformation behavior of the tubes during a hydraulic expansion process in a square cross-sectional die. This model is used to predict the forming pressure needed to hydroform a circular tube into a square cross-section and the thickness distribution of the product. FE simulations on tube expansion have also been carried out to compare the analytical results with those by the proposed mathematical model. The effects of the friction coefficient between the die and tube, upon the forming pressures needed and thickness distributions after the expansion process are discussed. Experiments of tube expansion in a square die using a self-designed apparatus are also conducted. The analytical results of forming pressure and thickness distributions of the formed parts are compared with the experimental measurements to verify the validity of this newly proposed model.     

A thermal-elastic analysis on laser-induced reverse bulging and plugging in circular brass foil

A new failure mode is observed in circular brass foils induced by laser beam. The new failure is based on the following experimental facts : (1) the peripheries of the circular brass foils are fixed and the surfaces of the foils are radiated by laser beam ; (2) the laser beam used is considered to be non-Gaussian spatially, actually an approximately uniform distribution limited in a certain size spot ; (3) the pulse on time of laser beam should be 250 μs, i.e. so called long duration pulse laser. The failure process consists of three stages ; i.e. thermal bulging, localized shear deformation and perforation by plugging. The word reverse in reverse bulging and plugging mode means that bulging and plugging occur in the direction of incident laser beam. To study the newly-discovered type of failure quantitatively, analytical solutions for the axisymmetric temperature field and deflection curve are derived. The calculated results show that the newly discovered failure mode is attributed to the spatial structure effect of laser beam indeed.     

Finite element modeling of plastic anisotropy induced by texture and strain-path change

Consideration of plastic anisotropy is essential in accurate simulations of metal forming processes. In this study, finite element (FE) simulations have been performed to predict the plastic anisotropy of sheet metals using a texture- and microstructure-based constitutive model. The effect of crystallographic texture is incorporated through the use of an anisotropic plastic potential in strain-rate space, which gives the shape of the yield locus. The effect of dislocation is captured by use of a hardening model with four internal variables, which characterize the position and the size of the yield locus. Two applications are presented to evaluate the accuracy and the efficiency of the model: a cup drawing test and a two-stage pseudo-orthogonal sequential test (biaxial stretching in hydraulic bulging followed by uniaxial tension), using an interstitial-free steel sheet. The experimental results of earing behavior in the cup drawing test, maximum pressure and strain distribution in bulging, and transient hardening in the sequential test are compared against the FE predictions. It is shown that the current model is capable of predicting the plastic anisotropy induced by both the texture and the strain-path change. The relative significance of texture and strain-path change in the predictions is discussed.     

Post-bifurcation analysis of a thin-walled hyperelastic tube under inflation

We consider the problem of bulging, or necking, of an infinite thin-walled hyperelastic tube that is inflated by an internal pressure, with the axial stretch at infinity maintained at unity. We present a simple procedure that can be used to derive the bifurcation condition and to determine the near-critical behaviour analytically. It is shown that there is a bifurcation with zero mode number and that the associated axial variation of near-critical bifurcated configurations is governed by a first-order differential equation that admits a locally bulging or necking solution. This result suggests that the corresponding bifurcation pressure can be identified with the so-called initiation pressure which featured in recent experimental studies. This is supported by good agreement between our theoretical predictions and one set of experimental data. It is also shown that the Gent material model can support both bulging and necking solutions whereas the Varga and Ogden material models can only support bulging solutions. Relevance of the present method to the study of non-linear wave propagation in a fluid-filled distensible tube is also discussed.     

<Emphasis Type="Italic">Ex Vivo</Emphasis> Intervertebral Disc Bulging Measurement Using a Fibre Bragg Grating Sensor

Advances using optical fibres as sensors may represent an important contribution for development of minimally invasive techniques in biomedical and biomechanical applications. Concerning spine injuries, intervertebral disc (IVD) degeneration is a major clinical issue since it represents gross structural disruption and it is irreversible. Measuring biomechanical parameters of the IVD should contribute for better understanding on its mechanical response to external applied forces. The purpose of this study was to explore the potential of a Fibre Bragg Grating (FBG) sensor to measure strain caused by bulging of the intervertebral disc under axial compression. Disc bulging is a consequence of IVD compression and a technique to register this behaviour is addressed in this study. Needle-mounted sensors were already used to measure IVD pressure in cadaveric material. In this study we also explored the possibility of using needles only for sensor guiding and positioning leaving sensor directly in contact with the IVD material. An ex vivo porcine dorsal functional spinal unit was instrumented with a FBG sensor and submitted to axial compression. Results suggest the sensor’s ability to measure strain response to load. Bulging of the annulus fibrosus as a consequence of axial compression was confirmed using the FBG sensor. Hysteresis and viscoelastic behaviour were observable suggesting that energy is dissipated by the deformation of the annulus and that unloading time was insufficient for disc recovery. Nevertheless the relatively low strain sensitivity of the sensor as well as signal artefacts caused by transverse loading may constitute a problem in the analysis and interpretation of strain data. The technique may not be suitable for measurement of physiologic bulging being more indicative of the radial force exerted by the annulus.