Lower bound limit analysis using non‐linear programming

This paper describes a new formulation, based on linear finite elements and non‐linear programming, for computing rigorous lower bounds in 1, 2 and 3 dimensions. The resulting optimization problem is typically very large and highly sparse and is solved using a fast quasi‐Newton method whose iteration count is largely independent of the mesh refinement. For two‐dimensional applications, the new formulation is shown to be vastly superior to an equivalent formulation that is based on a linearized yield surface and linear programming. Although it has been developed primarily for geotechnical applications, the method can be used for a wide range of plasticity problems including those with inhomogeneous materials, complex loading, and complicated geometry. Copyright © 2002 John Wiley & Sons, Ltd.     

Micromechanical modelling of ductile crack growth in the binder phase of WC–Co

This study focuses on numerical simulation of ductile failure in the Co binder phase of WC–Co hardmetal. The growth of edge cracks under mode I loading is considered. A computational micromechanics approach is taken where the Co binder ligaments are explicitly represented in finite element models. An embedding technique is employed. Crystal plasticity theory is used to represent plastic deformation in the Co ligaments. Crack propagation in the binder is simulated using an element removal technique based on a modified Rice and Tracey model for ductile void growth, and fracture resistance curves are generated. Parameter studies are performed for variations in microstructrual parameters such as numbers of Co ligaments ahead of the crack tip and local Co volume fraction. The importance of thermal residual stresses and finite element mesh density are also investigated.     

Efficient and accurate approach for powder compaction problems

 In this paper, a new approach for powder cold compaction simulations is presented. A density-dependent plastic model within the framework of finite strain multiplicative hyperelastoplasticity is used to describe the highly nonlinear material behaviour; the Coulomb dry friction model is used to capture friction effects at die-powder contact; and an Arbitrary Lagrangian–Eulerian (ALE) formulation is used to avoid the (usual) excessive distortion of Lagrangian meshes caused by large mass fluxes. Several representative examples, involving structured and unstructured meshes are simulated. The results obtained agree with the experimental data and other numerical results reported in the literature. It is shown that, contrary to other Lagrangian and adaptive h-remeshing approaches recently reported for this type of problems, the present approach verifies the mass conservation principle with very low relative errors (less than 1% in all ALE examples and exactly in the pure Lagrangian examples). Moreover, thanks to the use of an ALE formulation and in contrast with other simulations, the presented density distributions do not present spurious oscillations. Received: 20 March 2002 / Accepted: 15 October 2002 The partial financial support of the Ministerio de Ciencia y Tecnología (grant number DPI 2001-2204) is gratefully acknowledged.     

铁基形状记忆合金0Cr12Ni6Mn16Si5的研制

探讨了0Cr12Ni6Mn16Si5形状记忆合金的工业化生产工艺,0Cr12Ni6Mn16Si5可以采用EF1+AOD+EF2进行冶炼.合金具有强的热敏感性和高温变形抗力,高温可加工热塑性温度区间比较窄,温度超过1*!100℃,塑性直线下降.铸锭的开锻温度为1*!080℃,终锻温度为900℃;锻造管坯只有采用挤压方法实施管材成型.挤压成型的毛管可以采用冷轧与冷拔相接合方法进一步加工成薄壁管.     

A mesomodel for the numerical simulation of the multiphasic behavior of materials under anisothermal loading (application to two low-carbon steels)

The determination of residual stresses induced by welding or heat treatment operations requires the use of complex models taking into account thermal, metallurgical and mechanical phenomena. In this paper, we propose a mechanical model in which each phase can follow its own constitutive law. This model also takes into account phase transformation plasticity, which is treated independently of the behavior of each phase. This model has been implemented into the French FEM code Castem 2000. The interest of the proposed method is that it allows one to mix any type of nonlinear behavior using Taylor homogenization hypothesis. There is no need to develop a theory to get the equations of the homogenized material law. Two numerical examples demonstrate the efficiency and the flexibility of this approach. The results obtained are compared to experimental values for a typical welding situation and a high-temperature response. This comparison seems to indicate that viscous effects in the materials have a significative influence on the residual stresses produced by welding.     

Analysis of crack tip plasticity for microstructurally small cracks using crystal plasticity theory

Crack tip plastic zone sizes and crack tip opening displacements (CTOD) for stationary microstructurally small cracks are calculated using the finite element method. To simulate the plastic deformation occurring at the crack tip, a two-dimensional small strain constitutive relationship from single crystal plasticity theory is implemented in the finite element code ANSYS as a user-defined plasticity subroutine. Small cracks are modeled in both single grains and multiple grains, and different crystallographic conditions are considered. The computed plastic zone sizes and CTOD are compared with those found using conventional isotropic plasticity theory, and significant differences are observed.     

轴对称空间问题的塑性解法

本文从理想材料轴对称弹塑性问题的所有基本方程出发,研究了全塑性时轴对称空间问题的塑性应变分布规律和应力场,并推导了控制应力分布的定解方程,给出了全塑性轴对称问题的通用解法。     

Efficient return algorithms for associated plasticity with multiple yield planes

A new return method for implicit integration of linear isotropic yield criteria is presented. The basic idea is to perform all the manipulations in the principal stress space and thereby achieve very simple formulae for calculating the plastic corrector stresses, based on the constant gradient of such criteria. The return formulae are in closed form and no iteration is required. The method accounts for three types of stress return: return to a single yield plane, to a discontinuity line at the intersection of two yield planes and to a discontinuity point at the intersection between three or more yield planes. The infinitesimal and the consistent elastoplastic constitutive matrix are calculated for each type of stress return, as are the conditions to ascertain which type of return is required. The method is exemplified with the Mohr–Coulomb yield criterion. Copyright © 2005 John Wiley & Sons, Ltd.     

An efficient return algorithm for non-associated plasticity with linear yield criteria in principal stress space

An efficient return algorithm for stress update in numerical plasticity computations is presented. The yield criterion must be linear in principal stress space and can be composed of any number of yield planes. Each of these yield planes may have an associated or non-associated flow rule. The stress return and the formation of the constitutive matrix is carried out in principal stress space. Here the manipulations simplify and rely on geometrical arguments. The singularities arising at the intersection of yield planes are dealt with in a straightforward way also based on geometrical considerations. The method is exemplified on non-associated Mohr–Coulomb plasticity throughout the paper.     

Numerical Study of Scale Effects of Ultra-Thin Nickel Beams

Recent experiments have shown that metallic materials display significant size effects when the characteristic length scale of non-uniform plastic deformation is close to a micron. Couple stress plasticity has been developed to explain such phenomena by Fleck and Hutchinson. The mechanical behaviors of ultra-thin nickel beams in different boundary conditions were studied with the hybrid element developed for couple stress plasticity before. Strong scale effects are found when the beam's thickness is close to the material characteristic length scale. Such phenomena will disappear if the beam' s thickness is greatly larger than the material characteristic length scale. The scale effect is the beams inherent property and it does not change with the change of support conditions.