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

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

功能梯度变曲率曲梁的几何非线性模型及其数值解

基于弹性曲梁平面问题的精确几何非线性理论,建立了功能梯度变曲率曲梁在机械和热载荷共同作用下的无量纲控制方程和边界条件,其中基本未知量均被表示为变形前的轴线坐标的函数。以椭圆弧曲梁为例,采用打靶法求解非线性常微分方程的两点边值问题,获得了两端固定功能梯度椭圆弧曲梁在横向非均匀升温下的热弯曲变形数值解,分析了材料梯度指数、温度参数、结构几何参数等对曲梁受力及变形的影响。     

功能梯度材料扁锥壳的热弹性大变形分析

研究了功能梯度材料扁薄锥壳在横向非均匀升温场中的几何非线性大变形问题.基于von Kármán几何非线性理论推导出了以中面位移为基本未知量的功能梯度扁薄锥壳在横向非均匀热载荷作用下的轴对称大挠度控制方程.采用打靶法数值求解所得非线性常微分方程边值问题,得到了锥壳的大挠度弯曲变形数值解.给出了锥壳的变形与其形状参数、载荷和材料参数等变化的特征关系曲线,分析和讨论了温度参数和材料梯度变化参数对变形的影响.     

功能梯度半空间体的热弹性问题研究

本文研究半空间非均匀各向同性功能梯度弹性体的热弹性问题.根据功能梯度材料热弹性体的运动方程,热传导方程以及本构方程,利用状态空间法推导出功能梯度材料的热力耦合微分方程,进一步利用特征值分析法和拉普拉斯逆变换进行求解,得到功能梯度热弹性体在分别只受温度和应力荷载作用时,时域内位移,温度,应力三个物理量的解析解.通过图表,数值分析了半空间功能梯度材料弹性体在热-力荷载作用下各物理场的响应,数值曲线的变化趋势反映了热-力荷载之间的耦合效应.这可以很好的利用在材料设计领域,通过控制材料的梯度参数来控制物理场的极值,从而为未来该材料的工厂加工设计提供理论支持.     

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

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

基于物理中面FGM简支梁的弯曲行为

基于一阶非线性梁理论,利用物理中面概念导出了FGM梁的基本方程,分析了热载荷作用下简支FGM梁的弯曲行为.当坐标面置于功能梯度材料(FGM)梁的物理中面上时,其本构方程中,面内力与弯矩并不耦合,使得问题的控制方程以及边界条件得以简化.分析中假设功能梯度材料性质只沿梁厚度方向、并按成分含量的幂指数形式变化;利用打靶法数值地求解了所得方程.数值结果表明:热载荷作用下,夹紧FGM梁发生过屈曲变形,而简支梁则发生较为复杂的热弯曲变形;在同一热载荷作用下,简支FGM梁将会产生三种构形问题;剪切变形对夹紧FGM梁的热变形影响比简支梁更明显.     

功能梯度截顶圆锥壳的热弹性弯曲精确解

研究了功能梯度材料截顶圆锥壳在横向机械载荷与非均匀热载荷同时作用下的变形问题. 基于经典线性壳体理论推导出了以横向剪力和中面转角为基本未知量的功能梯度薄圆锥壳轴对称变形的混合型控制方程. 假设功能梯度圆锥壳的材料性质为沿厚度方向按照幂函数连续变化的形式. 然后采用解析方法求解,得到了问题的精确解. 分别就两端简支和两端固支边界条件,给出了圆锥壳的变形随其载荷、材料参数等变化的特征关系曲线,重点分析和讨论了载荷参数与材料梯度变化参数对变形的影响.      

热环境下贴压电层的功能梯度材料板的自由振动和动力响应

基于Reddy高阶剪切变形理论和广义Kármán型方程,用双重Fourier级数展开法求得了热环境下带压电层的功能梯度复合材料混合层合板的自由振动及动力响应的解析解,分析中考虑了材料热物参数对温度变化的依赖性,讨论了环境温度和控制电压对固有频率及动力响应的影响。     

热荷载作用下Timoshenko功能梯度夹层梁的静态响应

在精确考虑轴线伸长和一阶横向剪切变形的基础上建立了Timoshenko功能梯度夹层梁在热载荷作用下的几何非线性控制方程.采用打靶法数值求解所得强非线性边值问题,获得了两端固支功能梯度夹层梁在横向非均匀升温作用下的静态热过屈曲和热弯曲变形数值解.分析了功能梯度材料参数变化、不同表层厚度和升温参数对夹层梁弯曲变形、拉-弯耦...     

径向弹性约束下功能梯度圆板的热过屈曲

研究了周边具有面内径向弹性约束功能梯度圆板在横向非均匀升温下的热过屈曲行为.基于von Karman薄板理论,推导出了横向非均匀加热功能梯度圆板在径向弹性约束作用下的位移形式的轴对称热过屈曲控制方程.假设功能梯度材料性质沿厚度方向按幂函数连续变化,采用打靶法求解得到非线性常微分方程边值问题,获得了周边简支和夹紧条件下功能梯度圆板的热过屈曲响应.定量分析了径向弹性约束对圆板的临界屈曲温度载荷以及热过屈曲变形的影响,给出了不同弹性约束刚度功能梯度圆板的热过屈曲平衡路径和平衡构形.数值结果表明,径向弹性约束对圆板的热过屈曲平衡路径的影响显著,随着约束刚度的减小,临界屈曲温度载荷增大.     

A new integration algorithm for finite deformation of thermo-elasto-viscoplastic single crystals

An algorithm for integrating the constitutive equations in thermal framework is presented, in which the plastic deformation gradient is chosen as the integration variable. Compared with the classic algorithm, a key feature of this new approach is that it can describe the finite deformation of crystals under thermal conditions. The obtained plastic deformation gradient contains not only plastic defor- mation but also thermal effects. The governing equation for the plastic deformation gradient is obtained based on ther- mal multiplicative decomposition of the total deformation gradient. An implicit method is used to integrate this evo- lution equation to ensure stability. Single crystal 1 100 aluminum is investigated to demonstrate practical applications of the model. The effects of anisotropic properties, time step, strain rate and temperature are calculated using this integration model.     

NUMERICAL SIMULATION OF DEFORMATION BEHAVIOR OF 22MnB5 BORON STEEL AT ELEVATED TEMP ERATURES AND EXPERIMENTAL VERIFICATION

The effects of strain, strain rate and temperature on the mechanical behavior of 22MnB5 boron steel deformed isothermally under uniaxial tension tests and the experimental characterization of 22MnB5 boron steel in the austenitic region have been investigated. Based on the crystal plasticity theory and thermal kinematics, an improved integration model is presented. In this model, the elastic deformation gradient is the integration variable of the governing equation, which contains not only the elastic deformation but also the thermal effects. In the coupled thermo- mechanical process, this model can reveal the evolution of microstructures such as the rotation of a single crystal and the slip systems in each of them. The plastic behavior of the boron steel can be well described by the presented model.     

延性单晶非均匀变形的三维有限元模拟

对延性单晶在拉伸载荷作用下的应变局域化和颈缩等非均匀变形过程进行了三维有限元数值模拟。将相关晶体塑性本构模型及一种新的数值积分方法补充到ABAQUS6.1商用有限元软件中。该方法的特点是，利用晶体塑性的动力学方程，获得一个关于晶体弹性变形梯度的演化方程，采用半隐式积分方案进行求解。本文推导出一种新的应力变本构矩阵。按此方式更新本构矩阵，计算速度和计算稳定性大大提高。加载方式，边界条件和变形程度等因素影响着滑移系的启动状况，这是平面模型所不能预测的。本文利用三维有限元方法模拟了不同取向下滑移系的启动状况，全面地考虑了FCC单晶材料12个可能滑移系在变形过程中的启动状况，合理地模拟了FCC面心立方单晶沿不同取向加载时晶轴旋转导致的应变局域化和颈缩等非均匀变形过程。     

有限变形下多晶晶体塑性模型算法及应用

用Sanna和Zacharia^[1]所提出的延性单晶本构模型的积分算法和Taylor多晶模型假设研究了时间步长和硬化模型的选取对多晶集合体的应力应变响应和织构演化的影响。该算法是利用变形梯度乘法分解获得弹性变形梯度演化方程，用增量迭代法积分该方程，显式更新各滑移上的临界分切剪应力。算例的结果表明该算法具有时间步大，计算效率高的特点，另外，不同硬化模型的选取对多晶集合体应力应变响应的预测有明显的影响但对织构演化的预测影响不大。     

考虑晶体滑移面分解正应力的细观损伤模型

材料内部的解理、滑移面剥离等细观损伤是引起宏观失效的根源, 从细观尺度研究损伤的发生和发展有助于深入认识材料的变形和失效过程. 本文基于晶体塑性理论, 从滑移系的受力和变形出发研究材料的细观损伤, 建立了考虑滑移面分解正应力的细观损伤模型, 为晶体材料解理断裂的分析提供了新方法. 首先, 在晶体弹塑性变形构型的基础上引入损伤变形梯度张量的概念, 从变形运动学着手建立了考虑损伤能量耗散的本构方程, 并推导了塑性流动方程与损伤演化方程; 然后, 建立了相应的数值计算方法, 给出了应力与状态变量的更新算法, 推导了Jacobian矩阵的表达式; 接着, 以$[100]$取向的单晶铜材料为例, 通过有限元计算与试验结果的对比, 并采用粒子群优化算法标定了11个材料细观参数; 最后, 将所提细观损伤模型应用于RVE单轴拉伸过程的模拟, 得到了考虑损伤影响的应力应变曲线, 并分析了材料的塑性流动与损伤演化过程. 结果表明, 本文所提模型能够计算材料在受载过程中的损伤累积效应, 合理反映晶体材料的细观损伤机理.      

Deformation and failure mechanisms of single crystal superalloy under temperature gradient

The crystallographic constitutive model under temperature gradient is developed and introduced to study the deformation and failure mechanisms of single crystal superalloy. Tensile tests of thin-walled pipe specimen at different temperatures without cooled air flow were carried out. Based on the experimental results, the temperature dependence of constitutive model was studied and the basic parameters of constitutive model were obtained. To investigate the deformation and failure mechanisms, the thin-walled pipe specimen with cooled air flow under temperature gradient were tested. Considered the fluid-solid interface (FSI), a finite element method (FEM) was proposed to simulate the process of tension. In FEM, the activation rate of slip system was defined as the failure law of specimen under temperature gradient. The simulation result was in good agreement with the experiment result. Furthermore, the fracture surface of the specimen was observed by the scanning electron microscopy (SEM). The microstructure revealed that the slip deformation belonged to {1 1 1} crystalplane is a principal failure mechanism of single crystal superalloy under temperature gradient. The results of the SEM also implied that the proposed FEM method can be used to systemically study the deformation and failure behavior of single crystal superalloy cooled blade.     

Effects of microscopic boundary conditions on plastic deformations of small-sized single crystals

The finite deformation version of the higher-order gradient crystal plasticity model proposed by the authors is applied to solve plane strain boundary value problems, in order to obtain an understanding of the effect of the higher-order boundary conditions. Numerical solutions are carried out for uniaxial plane strain compression of a single crystal block and for uniform pure bending of a single crystal foil. The compressed block has loading surfaces that are penetrable or impenetrable to dislocations. This allows for a study of the two types of higher-order boundaries available, and a significant effect of higher-order boundary conditions on the overall deformation mode of the block is observed. The bent foil has free surfaces through which dislocations can go out of the material, and we observe a strong size-dependent mechanical response resulting from the surface condition assumed.     

基于晶体塑性理论研究铝材料高压高应变率下的强度特性

为了了解金属材料在极端加载下复杂动态响应过程中的多种机制和效应，重点针对Al材料在高压、高应变率加载下的塑性变形机制，在经典晶体塑性模型的基础上，对其中的非线性弹性、位错动力学和硬化形式进行改进，建立适用于高压、高应变率加载下的热弹-黏塑性晶体塑性模型。该模型可以较好地描述单晶铝和多晶铝材料屈服强度随压力的变化过程，相比宏观模型，用该模型还获得了多晶Al材料在冲击加载下的织构演化规律，揭示了织构择优取向行为和压力的关系。     

Mechanism-based strain gradient crystal plasticity—II. Analysis

In part I of this series (Mechanism-based strain gradient crystal plasticity—I. Theory. J. Mech. Phys. Sol. (2005), accepted for publication), we have proposed a theory of mechanism-based strain gradient crystal plasticity (MSG-CP) to model the effect of inherent anisotropy of a crystal lattice on size-dependent non-uniform plastic deformation at micron and submicron length scales. In the present paper, several example problems are investigated to show how crystal anisotropy is reflected by the MSG-CP theory.