两种弹塑性算法的比较

分析了弹塑性一致切线模量算法与径向返回算法，发现一致切线模量算法并不如有关文献报道的那样，比径向返回算法有快得多的收敛速度，两个算例表明了这一点。这使对一致切线模量算法有比较全面的看法。     

有限元实现含应变非强化区的双面循环塑性模型

在有限单元方法的框架下,由广义径向返回法来实现含有应变非强化区的循双面环塑性本构模型.在小变形的假设下,本文详细描述了有限元方法实现率无关的双面循环塑性本构模型的完整算法.采用向后Euler隐式算法进行离散积分,其最终可化简为一个非线性方程,并且由Peg-asus方法来求解.对径向返回法进行线性化,可以得到切向刚度算子.最后,文中通过两个算例表征了离散积分的精确性和材料模型的正确性.     

软岩三维非线性统一弹粘塑性本构模型有限元分析

在作者已建立的软岩三维非线性统一弹粘塑性软化本构模型基础上,讨论了软岩三维非线性统一弹粘塑性本构模型有限元分析过程,即推导了软岩三维非线性统一弹粘塑性屈服面(破坏面)的流动矢量表达式,便于编程与通用的有限元程序接口,进行二次开发;利用张量求导的原则,导出了流动矢量对应力求导公式中关键矩阵的具体表达式,并结合平面应力、平面应变和轴对称问题的特点进行了简化,得出了适合于分析二维问题的关键矩阵的表达式;采用流动矢量奇异性的统一处理方法,有效地处理了各种角点奇异性问题;最后给出了有限元分析的流程图,便于编制有限元程序.     

粘塑性薄壁管中复合应力波的传播特性研究

以本构关系一般理论为基础,导出了计及材料功硬化效应和应变率硬化效应的粘塑性薄壁管的本构关系及管中复合应力波的控制方程,应用有限差分方法研究了在压扭复合冲击载荷作用下粘塑性薄壁管中复合应力波的传播特性与演化规律,分析了复合应力波的耦合效应以及薄壁管中粘塑性参数和功硬化效应对复合应力波传播与演化规律的影响,并对有关现象进行了分析和解释。     

弹塑性板壳结构非线性有限元分析

本文根据塑性流动理论的基本公式，由隐式积分导出了与路径无关的变量更新算法和一致切线模量。采用单元广义应力应变直接离散塑性流动定律，构造了杂交应力单元一致切线刚度矩阵的显式表达式，编制了结构有限元程序ＳＡＦＥ，数值算例表明：本文的计算方法和计算程序是正确可靠的，可用于弹塑性板壳结构的非线性分析，计算结果屈曲临界载荷和极限承载能力。     

三维非线性有限元与弹性边界元耦合数值方法

本文系统地讨论了以下三个问题:(1) 有限元与边界元耦合中的几个数值问题,其中包括:边界积分方程的凝聚、等效刚度矩阵的对称化及面力不连续的处理;(2) 弹塑性有限元与弹性边界元的耦合;(3) 弹粘塑性有限元与弹性边界元的耦合及数值计算稳定性条件。     

Bammann-Chiesa-Johnson粘塑性本构模型材料参数的一种识别方法

粘塑性本构模型能否成功模拟金属高应变率大应变变形过程依赖于材料参数识别结果的好坏。由于BCJ模型考虑了应变率、温度与材料硬化之间的耦合效应以及应变率、温度历史效应,同时模型中包含了多个材料参数,因此很难通过试验直接识别模型的材料参数。本文针对BCJ模型中的耦合效应和历史效应,基于对模型中材料参数物理涵义的界定,给出了一种对材料参数解耦、分离并进行估计的方法,获得了模型材料参数估计公式,估计了材料参数的取值范围。在此基础上,编制了BCJ模型应力积分径向返回算法和粒子群优化算法的计算程序,应用重新设计了BCJ模型耦合效应和历史效应的反分析方法,在参数取值范围内对材料参数进行了优化识别。以OFHC Cu为例,应用提出的识别方法对BCJ模型的材料参数进行了识别,计算结果和试验结果符合较好。     

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

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

无限粘弹介质中圆孔时变问题的解析分析

本文推导粘弹介质中圆孔孔径时变时的应力和位移.由粘弹解与弹性解的对应关系得到粘弹时变应力解.用直接解方程法求径向位移,最终归结为求解关于待定函数的l阶非齐次微分方程.将半径时变函数泰勒展开,用幂级数解法得到一般情况下的解.在寻找定解条件时,采用了对待定函数的光滑化处理,认为在t=0的微小邻域内函数仍满足微分方程,通过积分得到与待定系数数目相同的定解条件,从而获得本问题径向位移解析解.对Maxwell粘弹模型的求解证明了该法的可靠性.文中解适用于任意粘弹模型和孔径任意时变的情况.     

耦合蠕变损伤的粘塑性模型用于疲劳-蠕变交互作用的研究

该文对550℃下30CrMoNiV5-11转子钢进行了轴向应力控制的带保持时间的疲劳实验.按照Le-maitre有效应力的概念,采用Kachanov损伤演化方程,推导了耦合损伤的粘塑性流动方程和随动硬化方程,使用返回映射和向后欧拉方法实现应力更新,并在有限元软件ANSYS的接口程序UMAT中嵌入用户程序.对试样疲劳-蠕变交互行为进行模拟,结果表明,耦合蠕变损伤变量后,模拟结果较不含损伤变量模型的模拟结果精度大为提高.     

Calculation for path-domain independentJ integral with elasto-viscoplastic consistent tangent operator concept-based boundary element methods

This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals. The project supported by National Natural Science Foundation of China (9713008) and Zhejiang Natural Science Foundation Special Funds No. RC.9601     

Implicit BEM formulations for usual and sensitivity problems in elasto-plasticity using the consistent tangent operator concept

This paper presents boundary element method (BEM) formulations for usual and sensitivity problems in (small strain) elasto-plasticity using the concept of the local consistent tangent operator (CTO). “Usual” problems here refer to analysis of nonlinear problems in structural and solid continua, for which Simo and Taylor first proposed the use of the CTO within the context of the finite element method (FEM). A new implicit BEM scheme for such problems, using the CTO, is presented first. A formulation for sensitivity analysis follows. It is shown that the sensitivity of the strain increment, associated with an infinitesimal variation of some design parameter, solves a linear problem which is governed by the (converged value of the) same global CTO as the one that appears in the usual problem. Numerical results for both usual and sensitivity problems are shown for a one-dimensional example. They demonstrate the effectiveness of the present approach. In particular, accurate sensitivities with respect to material parameters (e.g., exponent of the power-type hardening law) are obtained even with few integration cells and for large load increments.     

Design sensitivity analysis for parameters affecting geometry,elastic–viscoplastic material constant and boundary condition by consistent tangent operator-based boundary element method

This paper presents a design sensitivity analysis method by the consistent tangent operator concept-based boundary element implicit algorithm. The design variables for sensitivity analysis include geometry parameters, elastic–viscoplastic material parameters and boundary condition parameters. Based on small strain theory, Perzyna’s elastic–viscoplastic material constitutive relation with a mixed hardening model and two flow functions is considered in the sensitivity analysis. The related elastic–viscoplastic radial return algorithm and the formula of elastic–viscoplastic consistent tangent operator are derived and discussed. Based on the direct differentiation approach, the incremental boundary integral equations and related algorithms for both geometric and elastic–viscoplastic sensitivity analysis are developed. A 2D boundary element program for geometry sensitivity, elastic–viscoplastic material constant sensitivity and boundary condition sensitivity has been developed. Comparison and discussion with the results of this paper, analytical solution and finite element code ANSYS for four plane strain numerical examples are presented finally.     

Shape sensitivity analysis of elastic structures by multi-domain boundary element method

Based on the direct differentiation of boundary element equations to design variables, an algorithm of the boundary element method which is suitable for the shape sensitivity analysis of serial multi-domain elastic structures is built. By the algorithm, it is not necessary to form and solve the equations for the structure as a whole. Each domain can be solved separately. Two computational examples are presented to demonstrate the accuracy and versatility of the algorithm.This work has been partially supported by National Natural Science Foundation of China     

Analytic Sensitivity Analysis for Shape Optimization

ＩｎｔｒｏｄｕｃｔｉｏｎＳｈａｐｅｏｐｔｉｍｉｚａｔｉｏｎｆｏｒｅｌａｓｔｏｍｅｒｉｓａｎｉｍｐｏｒｔａｎｔｍｅａｓｕｒｅｔｏｒｅｌｉｅｖｅｓｔｒｅｓｓｃｏｎｔｒｉｔｉｏｎ ,ｐｒｅｖｅｎｔｂｒｅａｋｉｎｇ ,ｉｍｐｒｏｖｅｌｏａｄ＿ｂｅａｒｉｎｇｃａｐａｃｉｔｙ .Ｓｔｒｕｃｔｕｒａｌｒｅｐｅｔｉｔｉｏｎａｎａｌｙｓｉｓ,ｓｅｎｓｉｔｉｖｉｔｙａｎａｌｙｓｉｓａｎｄｏｐｔｉｍａｌａｌｇｏｒｉｔｈｍａｒｅｔｈｒｅｅｉｍｐｏｒｔａｎｔｌｉｎｋｓｏｆｓｈａｐｅｏｐｔｉｍｉｚａｔｉｏｎ .Ｉｎｔｈｅｃｏｕｒｓ…     

基于等几何边界元法的声学敏感度分析

基于非均匀有理B样条(NURBS)曲面建模技术,边界物理量同样用NURBS基函数插值,推导出三维声场等几何边界积分方程。进一步以控制点为设计变量,用直接微分法推导出等几何敏感度边界积分方程,给出声场声压对形状参量的敏感度。针对边界积分方程中的超奇异积分,使用奇异相消技术并结合Cauchy主值积分和Hadamard有限部分积分处理,给出了超奇异积分的NURBS插值半解析表达式。数值算例验证了本文算法求解声学结构形状敏感度的有效性,为声学结构的整体形状优化打下基础。     

Structural–acoustic sensitivity analysis of radiated sound power using a finite element/ discontinuous fast multipole boundary element scheme

The complete interaction between the structural domain and the acoustic domain needs to be considered in many engineering problems, especially for the acoustic analysis concerning thin structures immersed in water. This study employs the finite element method to model the structural parts and the fast multipole boundary element method to model the exterior acoustic domain. Discontinuous higher‐order boundary elements are developed for the acoustic domain to achieve higher accuracy in the coupling analysis. Structural–acoustic design sensitivity analysis can provide insights into the effects of design variables on radiated acoustic performance and thus is important to the structural–acoustic design and optimization processes. This study is the first to formulate equations for sound power sensitivity on structural surfaces based on an adjoint operator approach and equations for sound power sensitivity on arbitrary closed surfaces around the radiator based on the direct differentiation approach. The design variables include fluid density, structural density, Poisson's ratio, Young's modulus, and structural shape/size. A numerical example is presented to demonstrate the accuracy and validity of the proposed algorithm. Different types of coupled continuous and discontinuous boundary elements with finite elements are used for the numerical solution, and the performances of the different types of finite element/continuous and discontinuous boundary element coupling are presented and compared in detail. Copyright © 2016 John Wiley & Sons, Ltd.     

Design of microstructure-sensitive properties in elasto-viscoplastic polycrystals using multi-scale homogenization

Evolution of properties during processing of materials depends on the underlying material microstructure. A finite element homogenization approach is presented for calculating the evolution of macro-scale properties during processing of microstructures. A mathematically rigorous sensitivity analysis of homogenization is presented that is used to identify optimal forging rates in processes that would lead to a desired microstructure response. Macro-scale parameters such as forging rates are linked with microstructure deformation using boundary conditions drawn from the theory of multi-scale homogenization. Homogenized stresses at the macro-scale are obtained through volume-averaging laws. A constitutive framework for thermo-elastic–viscoplastic response of single crystals is utilized along with a fully-implicit Lagrangian finite element algorithm for modelling microstructure evolution. The continuum sensitivity method (CSM) used for designing processes involves differentiation of the governing field equations of homogenization with respect to the processing parameters and development of the weak forms for the corresponding sensitivity equations that are solved using finite element analysis. The sensitivity of the deformation field within the microstructure is exactly defined and an averaging principle is developed to compute the sensitivity of homogenized stresses at the macro-scale due to perturbations in the process parameters. Computed sensitivities are used within a gradient-based optimization framework for controlling the response of the microstructure. Development of texture and stress–strain response in 2D and 3D FCC aluminum polycrystalline aggregates using the homogenization algorithm is compared with both Taylor-based simulations and published experimental results. Processing parameters that would lead to a desired equivalent stress–strain curve in a sample poly-crystalline microstructure are identified for single and two-stage loading using the design algorithm.     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 2D ELASTICITY

针对边界元法中高阶单元中几乎奇异积分计算难题,解剖了二维边界元法高阶单元的几何特征,定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近,从而分离出积分核中主导的奇异函数部分,其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分,再对奇异核函数积分导出解析公式,从而建立了一种新的半解析法,用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例,采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力,并与线性边界元正则化算法结果作了比较,结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构,采用正则化算法的线性边界元分析比有限元有显著优势,而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。     

An implicit algorithm for hypoelasto-plastic and hypoelasto-viscoplastic endochronic theory in finite strain isotropic–kinematic-hardening model

This paper is concerned with objective stress update algorithm for elasto-plastic and elasto-viscoplastic endochronic theory within the framework of additive plasticity. The elastic response is stated in terms of hypoelastic model and endochronic constitutive equations are stated in unrotated frame of reference. A trivially incrementally objective integration scheme for rate constitutive equations is established. Algorithmic modulus consistent with numerical integration algorithm of constitutive equations is extracted. The implementation is validated by means of a set of simple deformation paths (simple shear, extension and rotation), two benchmark test in nonlinear mechanics (the necking of a circular bar and expansion of a thick-walled cylinder), a test which demonstrates the capabilities of the proposed model in simulation of cyclic loading and ratcheting in finite strain case (cyclically loaded notched bar) and finally, the analysis of a tensile test, which presents a shear band with a finite thickness independent of the finite element mesh using endochronic viscoplastic constitutive model.