基体微尺度效应对弹塑性多孔洞材料本构势及孔洞长大的影响

基于尺度相关的应变梯度塑性SG理论，对含孔洞的球形体胞模型进行了解析分析，得到了在基体梯度塑性环境下球形孔洞的演化规律，给出了弹塑性孔洞材料的宏观屈服面方程．与现有的基于尺度无关塑性理论的Curson模型相比，该模型考虑了基体材料的特征长度l与孔洞半径α之比λ(λ＝l／α)对多孔材料宏观屈服面和孔洞演化规律的影响．当不计基体材料的塑性梯度效应和硬化效应时，该模型能退化到经典的Curson模型．     

单晶有限变形的热力耦合计算模型

以弹性变形梯度作为基本变量,结合热力学理论构造了单晶有限变形的热、力耦合计算模型。该模型考虑了温度、变温速率以及塑性耗散等条件对单晶有限变形的影响,相对于传统的以弹性变形梯度为基本变量的晶体塑性模型,算法能够体现温度效应的影响。采用隐式的积分方法对建立的控制方程进行计算以保证求解过程的稳定。以1100Al单晶为例计算了不同升温、降温速率,以及不同应变率影响下的材料应力-应变的响应。结果表明,模型能较好地反映变温过程中,单晶各向异性性质的演化以及应力、应变之间关系的变化。     

薄壁管预扭冲击拉伸实验装置的研制

材料在复合应力下的拉压或扭转行为常常与单轴拉压或者单轴扭转应力下不同,拉压以及扭转常体现出相互影响的特点。而复合应力下的塑性波可以有效地反映材料在复合应力下的本构行为特点。本文通过对霍普金森压杆进行改造,建立了一套薄壁管预扭冲击拉伸的实验装置,可以在薄壁管内产生拉扭耦合塑性波,并对率相关材料304不锈钢进行了薄壁管预扭拉伸实验研究,得到了该材料的拉扭耦合塑性波。结果表明,304不锈钢薄壁管的拉扭耦合塑性波具有明显的耦合快波和耦合慢波的双波结构,并且快波慢波之间没有恒值区间隔。同时也表明,该装置可以实现预期效果,并且可以有效避免薄壁管屈曲的产生。     

梯度多胞牺牲层的抗爆炸分析

运用一维非线性塑性冲击波模型和细观有限元模型对密度梯度多胞牺牲层的抗爆炸性能进行了分析。基于率无关的刚性-塑性硬化模型,建立了描述冲击波在多胞牺牲层中传播的控制方程,分别给出了正、负密度梯度多胞材料在指数型爆炸载荷作用下的响应特性。研究了可正好吸收爆炸能量的梯度多胞牺牲层的临界厚度与载荷强度、覆盖层质量、多胞材料的密度梯度等参数之间的关系,给出了以临界厚度和支撑端应力峰值为指标的密度梯度设计图。运用二维细观有限元模型验证了基于非线性塑性冲击波模型的抗爆炸分析的有效性。     

金属切屑塑性流动的稳定性

对金属正交切削过程中切屑形成机制和材料塑性流动行为进行实验研究和理论分析. 通过对4 种常用金属材料正交切削过程的实验研究和切屑形貌的微观观察,确定了连续切屑转变成锯齿切屑的临界速度. 结果表明该临界速度与材料性能相关. 在实验观察基础上,提出描述材料正交切削过程的二维分析模型. 该模型假设切屑形成区为包括主剪切区和次剪切区的一个平行四边形. 载荷有主剪切区中的剪应力和次剪切区中的正压力;通过量纲分析得到描述材料正交切削过程的无量纲主控参数和无量纲形式的基本控制方程;应用线性稳定性分析方法建立平面应变状态下评价材料塑性流动稳定性的普遍准则;求得切屑形成区内材料塑性变形的速度和应力近似解. 讨论切屑形成、形貌转变以及相关的塑性失稳机制. 分析结果表明, 表征材料惯性与阻尼之比的无量纲参数— 雷诺数可以作为主控参数描述金属切削过程以及切屑材料塑性流动的稳定性.      

金属切屑塑性流动的稳定性

对金属正交切削过程中切屑形成机制和材料塑性流动行为进行实验研究和理论分析.通过对4种常用金属材料正交切削过程的实验研究和切屑形貌的微观观察,确定了连续切屑转变成锯齿切屑的临界速度.结果表明该临界速度与材料性能相关.在实验观察基础上,提出描述材料正交切削过程的二维分析模型.该模型假设切屑形成区为包括主剪切区和次剪切区的一个平行四边形.载荷有主剪切区中的剪应力和次剪切区中的正压力;通过量纲分析得到描述材料正交切削过程的无量纲主控参数和无量纲形式的基本控制方程;应用线性稳定性分析方法建立平面应变状态下评价材料塑性流动稳定性的普遍准则;求得切屑形成区内材料塑性变形的速度和应力近似解.讨论切屑形成、形貌转变以及相关的塑性失稳机制.分析结果表明,表征材料惯性与阻尼之比的无量纲参数—雷诺数可以作为主控参数描述金属切削过程以及切屑材料塑性流动的稳定性.     

功能梯度材料纯弯曲梁的弹塑性分析

假设功能梯度材料为一理想弹塑性材料,其弹性模量和屈服强度沿梁高度方向按照幂函数变化,在小变形及平截面假设下研究功能梯度材料纯弯曲梁的弹塑性性能.根据Mises屈服准则导出了纯弯曲梁的弹性极限弯矩的解析表达式,建立了梁在弹塑性状态时截面弯矩与截面弹、塑性区分布之间的关系式,给出了梁进入塑性极限状态时中性轴的位置以及塑性极限弯矩的解析计算公式.数值算例的结果表明,功能梯度材料梁的弹塑性性能与均匀材料梁不同,其屈服不一定首先产生于截面最大应力点,而可能有多种不同的屈服模态及相应的塑性扩展.弹性模量及屈服强度的梯度变化对功能梯度材料纯弯曲梁的中性轴位置、截面弹塑性应力分布以及塑性极限弯矩均有较大影响.研究结果可为功能梯度材料梁的弹塑性分析提供一定的参考.     

描述大应变率范围下材料响应的粘塑性本构模型

以位错动力学理论中的Ｏｒｗａｎ和Ｇｉｌｍａｎ关系为基础建立描述率相关材料非弹性响应的基本方程，选择材料准静态实验的单轴响应作为强化演化的规律，并考虑应变率敏感程度随变形产生变化的特性，建立了适用于大应变率范围内率相关材料的统一型粘塑性本构模型。对铝１１００－０在应变率范围１０－５～１０４ｓ－１内产生的有限塑性应变的单轴响应进行了理论预测，与Ｋｈａｎ和Ｈｕａｎｇ［１］的实验数据及模型预测结果进行了比较，结果表明本文模型具有较高的预测精度，在高应变率和较大应变下不容忽视率敏感参数随变形的变化。     

基于OTTOSEN准则的混凝土粘塑性力学模型

研究了混凝土的弹．粘塑性动态本构关系,提出一个应变率相关的非线性混凝土模型．模型改进Bicanic的混凝土塑性间断面变化规律和Ottosen的四参数屈服准则,建立冲击荷载下的混凝土本构方程,可以应用于地震和爆炸作用下混凝土材料响应的研究．根据混凝土动态实验结果对应变率和材料强度的关系提出合理假设,考虑以往被忽略的材料峰值应变和泊松比随应变率的变化．模型包含静水压力参数和动态塑性损伤因子,可以有效地反映混凝土的动态力学行为,为其动态问题的研究提供有益的思路和有效的工具．     

考虑前屈曲耦合变形时功能梯度简支梁的稳定性分析

在某些边界条件下,功能梯度材料（FGM）梁会由于拉弯耦合产生前屈曲耦合变形,而该变形对FGM梁的稳定性有影响。本文假设FGM梁的材料性质只沿厚度方向进行变化,基于经典非线性梁理论和物理中面概念,推导出FGM梁的平衡方程以及包含前屈曲耦合变形影响的屈曲控制方程,并用打靶法进行数值求解。讨论了前屈曲耦合变形、梯度指数以及材料性质的温度依赖等因素对FGM梁非线性变形和稳定性的影响。     

Investigation on gradient-dependent nonlocal constitutive models for elasto-plasticity coupled with damage

IntroductionIntheframeworkofconventionalplasticity ,materialinstabilityisoneoftheprincipalfactorsthatresultinginthestrainlocalizationphenomenon .Byusingtheterminology‘homogenized’ ,itisreferredtothefactthatinitialflawsandboundaryconditionsnecessarilyinduceanon_homogeneousstressstateinaspecimenduringtesting .Inparticularintheprocessofprogressivefailure ,theflawsandlocalstressconcentrationwillcausestronglyinhomogeneousdeformationofthespecimen[1,2 ].Asthedeformation_inducedfracture/damagephenom…     

Linear stability of WENO schemes coupled with explicit Runge–Kutta schemes

This paper focuses on the results of the linear stability analysis of the finite‐difference weighted essentially non‐oscillatory (WENO) schemes with optimal weights. The standard WENO schemes between the third and 11th order, the order‐optimised WENO schemes of the sixth and eighth order and the bandwidth‐optimised WENO schemes of the third and fourth order are considered. Several explicit Runge–Kutta schemes including the recently published strong stability‐preserving explicit Runge–Kutta schemes are considered for time discretisation. The stability limits as well as dissipation and dispersion properties dependent on the Courant–Friedrichs–Lewy number are presented for a hyperbolic model equation. The different combinations of space and time discretisation schemes are compared in terms of their accuracy and efficiency. For a parabolic model equation, the viscous term is discretised with high‐order central differences. The stability limits for the parabolic problem are presented as well. Numerical results of linear test cases are shown. Copyright © 2011 John Wiley & Sons, Ltd.     

边坡稳定的剪切带计算

为了解决边坡稳定分析中剪切带有限元网格的依赖性问题,采用梯度塑性理论,从本构关系中引入特征长度入手,建立计算模型。提出了一种8节点缩减积分的梯度塑性单元,并采用梯度塑性理论推导了Drucker-Prager屈服准则的软化模型的有限元格式,在ABAQUS中进行了二次开发,嵌入了本文提出的8节点单元和本构模型,并用ABAQUS软件进行了边坡剪切带的计算。计算结果表明,本文提出的方法消除了经典有限元计算的网格依赖性问题,可以得到与单元剖分无关的稳定的剪切带宽度。本文所提出的方法可适用于其他场合的剪切带计算。     

Predictions of bifurcation and instabilities during dynamic extension

Dynamic bifurcation and flow instabilities of cylindrical bars, made of an incompressible strain hardening plastic material, are investigated. A Lagrangian linear perturbation analysis is performed to obtain a fourth order partial differential equation which governs the evolution of the perturbation. The analysis shows that inertia slows down the growth of long wavelengths while bidimensional effects conjugated to strain hardening extinct short wavelengths. The present approach is applied successfully to the analysis of bifurcation and instabilities in (i) a rectangular block during plane strain extension, (ii) a circular bar during uniaxial extension. New results are obtained in the case of rate independent materials and a synthetical point of view is obtained for rate dependent behaviors.     

Finite volume scheme for the solution of fluid flow problems on unstructured non‐staggered grids

A recently developed non‐staggered methodology which uses the principle of applying fourth‐order dissipation to the governing pressure‐correction equation is developed so it can be applied to unstructured grids. A finite volume methodology is used for discretization. The fourth‐order dissipation term is found using second‐order gradient operators. This makes it straightforward to incorporate the dissipation term on unstructured grids. The new methodology is compared with solutions from a standard finite volume second‐order flow solver and is also tested for a standard laminar driven‐lid flow problem with grids systems that do not have a uniform structure. Finally, we demonstrate how the new methodology can be used to predict flow over a wavy boundary. Copyright © 2002 John Wiley & Sons, Ltd.     

A higher‐order unsplit 2D direct Eulerian finite volume method for two‐material compressible flows based on the MOOD paradigms

A higher‐order unsplit multi‐dimensional discretization of the diffuse interface model for two‐material compressible flows proposed by R. Saurel, F. Petitpas and R. A. Berry in 2009 is developed. The proposed higher‐order method is based on the concepts of the Multidimensional Optimal Order Detection (MOOD) method introduced in three recent papers for single‐material flows. The first‐order unsplit multi‐dimensional Finite Volume discretization presented by SPB serves as foundation for the development of the higher‐order unlimited schemes. Specific detection criteria along with a novel decrementing algorithm for the MOOD method are designed in order to deal with the complexity of multi‐material flows. Numerically, we compare errors and computational times on several 1D problems (stringent shock tube and cavitation problems) computed on 2D meshes with the second‐ and fourth‐order MOOD methods using a classical MUSCL method as reference. Several simulations of a 2D shocked R22 bubble in the air are also presented on Cartesian and unstructured meshes with the second‐ and fourth‐order MOOD methods, and qualitative comparisons confirm the conclusions obtained with 1D problems. These numerical results demonstrate the robustness of the MOOD approach and the interest of using more than second‐order methods even for locally singular solutions of complex physics models. Copyright © 2014 John Wiley & Sons, Ltd.     

Time dependent behavior of ferroelectric materials undergoing changes in their material properties with electric field and temperature

This paper presents a time dependent polarization constitutive model suitable for predicting nonlinear polarization and electro-mechanical strain responses of ferroelectric materials subject to various histories of electric fields. The constitutive model is derived based on a single integral form with nonlinear (electric field and temperature dependent) integrand. The total polarization consists of the time-dependent and residual components. The residual component of the polarization is due to polarization switching in the ferroelectric materials. We use an ‘internal clock’ concept to incorporate the effect of electric field on the rate of polarization. The corresponding strain response is determined through the use of third order piezoelectric constant and/or fourth order electrostrictive constant that vary with polarization stage. It is assumed that in absence of polarization, both piezoelectric and electrostrictive constants are zero. To incorporate the effect of temperature on the overall polarization behavior all material parameters in the constitutive model are allowed to change with the ambient temperature. We present numerical studies on the effect of time, temperature, and electric field on the response of ferroelectric material followed by verification of the constitutive model. Experimental data on lead zirconate titanate (PZT) materials available in the literature are used to verify the model.     

马氏体相变过程中形核、失稳及多界面微结构的动力学模拟

材料发生相变的过程中会出现失稳、滞后回线及多界面的微结构等复杂现象，而稳定性的丧失使其动力学方程的求解十分困难。对于形状记忆合金中的马氏体相变，相变过程中材料的等效杨氏模量变为负值，使得传统的动力学方程成为病态的，无法直接求解，必须要进行正则化。而相变的滞后回线与微结构的出现也说明经典的弹性理论不再适用，必须要引入新的能量项以能刻画这些现象。本文在非线性弹性理论的框架下，引入应变梯度界面能和位移非均匀能，利用变分原理建立了材料相变的一维动力学模型。高阶项的引入极大地改善了方程的性质，使数值求解成为可能。计算结果表明，该模型确能较有效地描述相变时的失稳与微结构。     

On the interplay between strain rate and strain rate sensitivity on flow localization in the dynamic expansion of ductile rings

In this work a stability analysis on flow localization in the dynamic expansion of ductile rings is conducted. Within a 1-D theoretical framework, the boundary value problem of a radially expanding thin ring is posed. Based on a previous work, the equations governing the stretching process of the expanding ring are derived and solved using a linear perturbation method. Then, three different perfectly plastic material constitutive behaviours are analysed: the rate independent material, the rate dependent material showing constant logarithmic rate sensitivity and the rate dependent material showing non-constant and non-monotonic logarithmic rate sensitivity. The latter allows to investigate the interaction between inertia and strain rate sensitivity on necking formation. The main feature of this work is rationally demonstrate that under certain loading conditions and material behaviours: (1) decreasing rate sensitivity may not lead to more unstable material, (2) increasing loading rate may not lead to more stable material. This finding reveals that the relation between rate sensitivity and loading rate controls the unstable flow growth. Additionally a finite element model of the ring expansion problem is built in ABAQUS/Explicit. The stability analysis properly reflects the results obtained from the numerical simulations. Both procedures, perturbation analysis and numerical simulations, allow for emphasizing the interplay between rate sensitivity and inertia on strain localization.