Adiabatic shear band localization of inelastic single crystals in symmetric double-slip process

Summary The main objective of the present paper is the development of a viscoplastic regularization procedure valid for an adiabatic dynamic process for multi-slips of single crystals. The next objective is to focus attention on the investigation of instability criteria, and particularly on shear band localization conditions.To achieve this aim, an analysis of acceleration waves is given, and advantage is taken of the notion of the instantaneous adiabatic acoustic tensor. If zero is an eigenvalue of the acoustic tensor, then the associated discontinuity does not propagate, and one speaks of a stationary discontinuity. This situation is referred to as the strain localization condition, and corresponds to a loss of hyperbolicity of the dynamical equations. It has been proved that for an, adiabatic process of rate-dependent (elastic-viscoplastic) crystal, the wave speed of discontinuity surface always remains real and different from zero. It means that for this case the initial-value problem is well-posed. However, for an adiabatic process of rate-independent(elastic-plastic) crystal, the wave speed of discontinuity surface can be equal zero. Then the necessary condition for a localized plastic deformation along the shear band to be formed is as follows: the determinant of the instantaneous adiabatic acoustic tensor is equal to zero. This condition for localization is equivalent to that obtained by using the standard bifurcation method. Based on this idea, the conditions for adiabatic shear band localization of plastic deformation have been investigated for single crystals. Particular attention has been focused on the discussion of the influence of thermal expansion, thermal plastic, softening and spatial covariance effects on shear band localization criteria for a planar model of an f.c.c. crystal undergoing symmetric primary-conjugate double slip. The results obtained have been compared with available experimental observations.Finally, it is noteworthy that the viscoplasticity regularization procedure can be used in the developing of an unconditionally stable numerical integration algorithm for simulation of adiabatic inelastic flow processes in ductile single crystals, cf. [21].The paper has been prepared within research programme sponsored by the Committee of Scientific Research under Grant 3 P404 031 07.     

金属板料成形的快速有限元分析

研究了基于形变理论的金属板形成快速有限元分析的方法－反向方法,并实现了计算程序。通过实例的计算结果和实验以及增量方法进行了比较,表明此方法能够定性地分析成形工件中的变形情况。由于计算速度快、建立分析模模型简单,此方法可用于设计早期估计零件的可成形性,以及部分工艺参数对形的影响。     

The effect of the sheet rigidity on strength of the spot weld tension shear joint

This paper presents some test and analysis results for a spot welded joint subjected to tensile and alternate load. The effect of sheet rigidity on the tensile strength and fatigue life of the spot welded joint is studied by using the stress intensity factorsK _I,K _II,K _III and an effective stress intensity factor K_max calculated by the finite element method for crack around the nugget. The results show that the effective stress intensity factor K_max is an essential parameter for estimating the fatigue life of the spot welded joint.     

Springback simulation and analysis of strong anisotropic sheet metals in U-channel bending process

The springback phenomenon of strong anisotropic sheet metals with U-channel bending as well as deep-drawing is numerically studied in detail by using Updating Lagrange FEM based on virtual work-rate principle, Kirchhoff shell element models and the Barlat-Lian planar anisotropic yield function. Simulation results are compared with a benchamark test. Very good agreement is obtained between numerical and test results. The focus of the present study is on the numerical simulation of the springback characteristics of the strong anisotropic sheet metals after unloading. The effects of the planar anisotropy coefficients and yield function exponent in the B-L yield function on the springback characteristics are discussed in detail. Some conclusions are given. The project supported by the National Natural Science Foundation of China (No. 19832020) and Provincial Natural Science Foundation of Jilin China (No.20000519)     

Analysis of shear localization in porous materials

A set of constitutive equations are derived based on the authors' lower bound yield loci for porous materials. By using these equations, the conditions for shear localization in porous materials are then investigated and the results are compared with those of Gurson's equations and the finite element analysis. The advantages of the present constitutive equations are fully illustrated.     

塑性各向异性板料本构关系及冲压成形的数值模拟

推导了具有一般屈服函数形式的弹塑性速率型本构关系;给出了用于板料成形的Hill塑性各向异性屈服模型下本构关系的具体形式;用有限元动力显式计算程序MSC/DYTRAN模拟了金属板料的冲压成形;通过算例分析,考察了塑性各向异性对凸耳形成和大小以及对成形模拟结果准确性的影响;数值结果和实验结果表明：各向（厚向）异性本构模型比各向同性本构模型更真实地反映了板料的成形性。     

Numerical investigation of the influence of shear band localization on the resonant behavior in the VHCF regime

The monitored resonant behavior of fatigue specimens of metastable austenitic stainless steel (AISI304) is correlated with its damage accumulation in the very high cycle fatigue (VHCF) regime. The resonant behavior is studied experimentally and shows a distinct transient characteristic. Microscopic examinations indicate that during VHCF a localized plastic deformation in shear bands arises on the specimen surface. Hence, this work focuses on the effect of damage accumulation in shear bands on the resonant behavior of AISI304 in the VHCF regime. A microstructural simulation model is proposed that takes into account specific mechanisms in shear bands proven by experimental results. The simulation model is solved numerically using the two-dimensional boundary element method and the resonant behavior is characterized by evaluating the force-displacement hysteresis loop. Simulation of shear bands agrees well with microscopic examinations and plastic deformation in shear bands influences the transient characteristic of the resonant behavior.     

A thin-sheet, combined tension and bending specimen

Validating stress intensity factor solutions for combined tension and bending is an arduous task because the necessary experimental data are not readily available. Toward this end, a tension and bending test specimen was designed to produce controllable levels of both tension stress and bending stress at the crack location. The specimen was made from 2024-T3 clad aluminum, which is commonly used in aircraft structures. The need for testing multiple specimens of various geometries and stress levels prompted the development of an analytical tool for specimen design. An extention of the Schijve line model, based on simple beam theory, was developed to calculate the stress distributions of tension and bending through the length of the specimen. A comparison of measured static strain levels with those predicted by the model showed the model to be accurate to within 5 percent, confirming its efficacy for specimen design. As expected, for the same remote stress (100 MPa), cracks in the tension and bending specimens grew faster than those in middle-cracked tension specimens.     

金属材料平面应变拉伸变形局部化有限元分析

本文对于涉及韧性金属大变形中颈缩与剪切带断裂一类高度非线性变形局部化问题进行了弹塑性有限元数值模拟。采用改进的Ｊ２形变理论微分形式公式与交叉三角形四边形单元有限元网格，详细研究了应变硬化指数及初始表面不均匀特性的平面应变拉伸颈缩和剪切带形成的综合影响，给出此类问题的断裂机制图。     

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

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

非经典大变形弹塑性本构方程及其算法研究

本文介绍了作者提出的非经典大变形弹塑性木构方程的增量形式，发展了相应的有限元分析方法。由于不采用屈服面的概念及相应的流动规律与加载准则，与传统的算法相比，本文提出的算法不仅简化了计算过程，而且改善了收敛速度和计算精度。作为验证和应用的例子，分析了圆柱试件的大变形拉伸颈缩过程和短圆柱的镦粗过程。与实验结果的比较证明了所发展方法的有效性。     

Flange earing and its control on deepdrawing of anisotropy circular sheets

The Hill's quadric anisotropy yield function and the Barlat-Lian anisotropy yield function describing well anisotropy sheet metal with stronger texture are introduced into a quadric-flow corner constitutive theory of elastic-plastic finite deformation suitable for deformation localization analysis. And then, the elastic-plastic large deformation finite element formulation based on the virtual power principle and the discrete Kirchhoff shell element model including the yield functions and the constitutive theory are established. The focus of the present research is on the numerical simulation of the flange earing of the deep-drawing of anisotropy circular sheets, based on the investigated results, the schemes for controlling the flange earing are proposed. Supported by NSFC(No. 19832020) and National Automobile Dynamic Simulation Laboratory of China.     

随机有限元方法在断裂分析中的应用

在幂律非线性随机有限元基础上,以单边裂纹板为例给出计算含量钢继裂参数,J（J积分）,δ（裂纹张开位移）,Δ（由裂纹引起的裂纹板上下底面相对位移）,θ（由裂纹引起的裂纹板上下底在相对转角）及其对基本随机变量变化率的方法和分析算例。     

金属板料拉延二次成形的有限元法模拟

建立了二次拉延成形加工的有限元分析计算模型;采用一次成形-回弹计算-二次成形的连续计算过程模拟了实际加工过程;有限元计算采用动力显式计算程序MSC/DYTRAN;用主从面（master surface-slave surface)模型定义板料和模具的接触,摩擦力用库仑定律计算;利用动力松弛法对成形过程中的回弹进行了计算。模拟结果和实际零件比较,证明模型合理,自满稳定,结果可靠,具有良好的应用价值。     

基于三阶剪切变形理论的悬臂输流管道自由振动

采用三阶剪切变形理论,结合有限元法研究了悬臂输流管道的自由振动问题．利用虚功原理建立了输流管系统的有限元方程,同时将悬臂端弹性支承以势能的形式引入到系统方程中,求解了系统前三阶的复频率．分别探讨了流体速度和弹簧刚度对系统复频率实部和虚部的影响,重点分析了弹簧刚度与前三阶固有频率间的关系．在弹性支承刚度为零的特例下,对比了本文结果与Timoshenko梁理论的结果,证明了本文方法的可靠性．研究发现系统固有频率的实部恒为负值,表明一端带有弹性支承的约束形式有利于提高悬臂输流管道自由振动的稳定性;流体的流动对管道振动起到了阻尼作用,在流动速度足够大的情况下,各阶振动固有频率均趋于零;当弹簧刚度为无穷大,且流体速度足够大时,输流管道将发生失稳．     

Numerical simulation of progressive shear localization and scale effect in cohesionless soil media

Objective and methodsA multi-yield surface model of the total strain formulation, which takes into account the confinement dependency and has numerical consistency with the crack analysis of reinforced concrete, is presented to simulate the progressive failure and the shear band formation in a cohesionless soil media.ResultsFor experimental verification, a comparison of the numerical analysis results with the reality is carried out in view of the progressive failure. The scale effect on the nominal strength is discussed with regard to the shear localization.ConclusionIt is shown that by means of the proposed concept the scale effect of the static bearing capacity of foundation can be well simulated.     

板材冲压成形变形局部化与动边界摩擦约束

提出适于具有强烈弯曲效应的大应变，大位移与界面摩擦三重非线性空间曲壳有限元列式。给出一种动边界摩擦约束的元级代数消元方法，就板材在方盒型冲头作用中压胀形过程，引入几种现代非经典塑性本构理论以及空单元技术，研究了成形过程中的变形局部化与断裂问题。     

Differential quadrature method for bending of orthotropic plates with finite deformation and transverse shear effects

Based on the Reddy ‘s theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature ( DQ ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert ( DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.     

Dynamic fracture of concrete compact tension specimen: Experimental and numerical study

Compared to quasi-static loading concrete loaded by higher loading rates acts in a different way. There is an influence of strain-rate and inertia on resistance, failure mode and crack pattern. With increase of loading rate failure mode changes from mode-I to mixed mode. Moreover, theoretical and numerical investigations indicate that after the crack reaches critical velocity there is progressive increase of resistance and crack branching. These phenomena have recently been demonstrated and discussed by O?bolt et al. (2011) on numerical study of compact tension specimen (CTS) loaded by different loading rates. The aim of the present paper is to experimentally verify the results obtained numerically. Therefore, the tests and additional numerical studies on CTS are carried out. The experiments fully confirm the results of numerical prediction discussed in O?bolt et al. (2011). The same as in the numerical study it is shown that for strain rates lower than approximately 50/s the structural response is controlled by the rate dependent constitutive law, however, for higher strain rates crack branching and progressive increase of resistance is observed. This is attributed to structural inertia and not the rate dependent strength of concrete. Maximum crack velocity of approximately 800 m/s is measured before initiation of crack branching. The comparison between numerical and experimental results shows that relatively simple modeling approach based on continuum mechanics, rate dependent microplane model and standard finite elements is capable to realistically predict complex phenomena related to dynamic fracture of concrete.     

Determining the shear angle,forces, and sizes of shearing elements during metal cutting

A new solution for the shear angle is proposed which is a generalization of the solution Lee-Shaffer solution and allows the determination of the cutting force and the shearing-element size. Merchant’s experimental data are processed taking into account the resistance force at the cutting edge, and it is shown that accounting for this force leads to the need to increase the internal friction angle in the calculated dependences in order to match theory with experiment. It is shown that the obtained theoretical results agree well with experimental results. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 177–186, January–February, 2009.