膨胀环受力历史对应力应变关系的影响

利用两类实验装置开展了无氧铜TU1膨胀环实验研究,发现:电磁膨胀环在加载阶段,样品受体力作用,满足均匀变形的假定;而爆炸膨胀环在加载阶段,样品内壁受面力冲击作用,不满足均匀变形的假定。针对这个差异,发展了一种考虑冲击阶段变形不均匀性的新方法,利用回收样品几何变形,将冲击阶段试样环内轴向塑性应变、径向塑性应变纳入等效塑性应变的计算中,通过修正后的方法更准确地获得了材料的应力应变关系。     

微薄梁三点弯曲的尺度效应研究

应用高灵敏度的力传感器以及时间序列电子散斑干涉法，同时测出了不同厚度纯镍薄片三点 弯曲试件的抗力与变形，得到薄梁中心点处的载荷与挠度曲线. 应用Fleck和Hutchinson 的偶应力理论，结合平面应变弯曲模型，建立了薄梁处于弹性状态和弹塑性状态的 控制方程, 应用Runge-Kutta法进行数值求解，并将计算得到的载荷-挠度曲线以及无量纲化弯矩-表面 应变曲线和实验结果进行了比较. 在理论计算过程中，没有拟合任何材料参数，所有的材料 参数均来自实验测量的结果，材料特征尺度也是根据Stolken和Evans的工作给出 的. 结果表明: 应用偶应力理论预测的结果和实验结果符合良好，而经典理论的预测结果与 实验不相符合.     

HRR理论的近似解析解

本文根据文献[3]给出的HRR 解的角函数数值表,利用数值拟合方法,给出了HRR理论的近似解析解.本文结果不仅形式简单,而且和精确解符合得很好.对幂硬化指数n在0和1之间的任何弹塑性材料均适用.给进一步研究裂纹尖端附近的弹塑性应力应变场以及工程安全评定提供了一个方便的解析工具.为了计算方便,我们把应力、应变、位移无量纲化,它们与实际应力、应变、位移的关系为     

近场动力学框架下的键基对应模型

传统键基近场动力学模型存在泊松比限制的问题,为了解决这一问题发展了态基近场动力学模型。其中非常规态的近场动力学模型通过定义非局部的变形梯度将近场力和传统应力关联起来,方便使用传统本构,但是态基近场动力学计算效率低于键基近场动力学。结合态基模型和键基模型的优势,提出键基对应模型,定义了基于键的变形梯度,参考连续介质力学中变形梯度的极分解过程,将键的变形分为转动部分和伸长部分。从而进一步定义了应变,通过物理方程求应力,进而计算键传递的近场力。键基对应模型解决了键基近场动力学的泊松比限制问题,也不需要进行近场动力学微观材料常数的计算。数值算例和理论推导证明了键变形梯度定义以及近场力计算方式的正确性。     

应力作用下NiTi形状记忆合金微结构演化的相场模拟及其本征应变率敏感性

基于Ginzburg-Landau动力学控制方程建立了NiTi形状记忆合金非等温相场模型,实现了对NiTi合金内应力诱导马氏体相变的数值模拟。同时将晶界能密度引入系统局部自由能密度,从而考虑多晶系统中晶界的重要作用。数值计算了单晶和多晶NiTi形状记忆合金在单轴机械载荷作用下微结构的动态演化过程和宏观力学行为,并重点研究了晶粒尺寸为60 nm的NiTi纳米多晶在低应变率下（0.0005～15 s?1）力学行为的本征应变率敏感性。研究结果表明,单晶NiTi合金系统高温拉伸-卸载过程中马氏体相变均匀发生,未形成奥氏体-马氏体界面。而纳米多晶系统在加载阶段出现了马氏体带的形成-扩展现象,在卸载阶段出现了马氏体带的收缩-消失现象。相同外载作用过程中,NiTi单晶系统的宏观应力-应变曲线具有更大的滞回环面积,拥有更优的超弹性变形能力。计算结果显示,在中低应变率下纳米晶NiTi形状记忆合金应力-应变关系表现出较明显的应变率相关性,应变率升高导致材料相变应力提升。这一应变率相关性主要源于相场模型中外加载荷速率与马氏体空间演化速度的相互竞争关系。     

纳米丝应变率效应的分子动力学模拟

采用分子动力学模拟了零温时不同应变率作用下纳米丝的拉伸力学行为．计算结果表明在缺乏热激活软化机制条件下，纳米丝应变率效应呈现出与宏观应变率试验结果相一致的特征．纳米丝在不同的应变率范围具有不同的变形机制．在应变率不敏感区和敏感区，纳米丝主要以位错运动作为塑性变形机制；在应变率突变区，纳米丝通过局部原子混乱区的持续扩展乃至整体结构的非晶化作为塑性变形机制．     

动态拉伸试验中试样应变测试的有效性分析

为了评估将试样通过胶粘连接到加载杆的Hopkinson杆装置所获得试样应变的有效性,对四种强度刚度差异较大的纤维增强复合材料进行了动态拉伸试验。试验时,试样通过环氧胶和杆夹层粘接,试样的应变分别按照Hopkinson杆一维应力波理论计算和试样上应变计直接准确测量得到。结果证明：对小变形碳纤维复合材料,按一维应力波理论计算的应变与试样上直接所测应变值偏差超过100%;对较大变形的GFRP和KFRP层合板,两者偏差小于40%。说明采用Hopkinson杆一维应力波理论计算的试样应变不准确。为修正不准确性,一是通过大量数据分析建立按一维应力波理论计算值与直接测量应变之间的关系式,用此式可使此试验装置获得有效的试样应变;二是借助ABAQUS有限元模拟分析得出粘胶层以及试样过渡弧段的变形,用一维应力波理论计算的应变减去此变形,也可获得有效的试样应变。     

固体力学跨尺度计算若干问题研究

本文展示了固体力学领域跨尺度计算的若干问题和研究概况。（1）建立位错动力学与有限元耦合DDD-FEM的计算模型,实现了能够基于纳米尺度离散位错运动机制计算分析连续介质有限变形晶体塑性问题,提出微纳尺度（200 nm~10 μm）晶体塑性流动应力解析公式,结合试验数据揭示了在无应变梯度下强度和变形的尺寸效应;（2）建立具有微相分离结构的纳米尺度粗粒化分子动力学模型CG-MD,计算获得聚脲材料在时域和频域下的存储模量和损耗模量,通过动态加载分析的DMA试验和超声波试验的数据验证,解决了连续介质尺度下微相分离高分子共聚物的设计难题;（3）通过数据驱动关联高分辨率的微米尺度CT影像和临床低分辨率的毫米尺度CT影像的特征值,建立了围关节松质骨小梁的等效模量和结构张量,为骨组织增材制造点阵结构设计和实现个性化骨缺损重建奠定了基础。     

应用固有应变理论计测对接焊钢管残余应力场的实验研究

本文采用改进后的固有应变法对对接焊管接头残余应力计测进行了理论分析,确定了测量方案,研究果结表明对于可处理成轴对称问题的对接焊管残余应力问题只需不太多的测量值就能获得残余应力全场。本文就一个对接焊钢管试件进行了实验,实验计测结果同其它文献较为一致。它显示了应用改进后的固有应变理论决定对接焊钢管残余场的有效性。     

纳米金属丝拉伸破坏及其应变率效应

建立了适于研究纳米金属快速变形破坏过程的分子动力学模型,并对不同应变率工况下不同截面尺寸单晶镍纳米丝的零温单向拉伸破坏过程进行了分子动力学模拟.模拟得到各种纳米镍丝的应力-应变曲线、屈服应变、屈服强度、断裂强度和初始弹性模量,提出了纳米金属丝快速变形力学性能的应变率效应预测公式并加以验证.计算表明金属纳米丝的屈服应变与尺寸和应变率无关,屈服强度、断裂强度和弹性模量与应变率呈对数关系.     

晶体塑性变形离散滑移模型及有限元分析

基于韧性单晶体实验现象，建立了描述晶体塑性变形的离散滑移模型．该模型的主要特点是：晶体滑移变形在宏观上是不均匀的，滑移带的分布是离散的．利用晶体塑性理论对模型进行了有限变形有限元分析，计算结果揭示了晶体滑移的离散行为，模拟的应力 应变曲线与实验曲线相吻合     

A new finite element method for strain gradient theories and applications to fracture analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J₂ deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.     

基于试验与有限元耦合技术的延性材料全程单轴本构关系获取方法

摘 要: 材料拉伸直至断裂的全程单轴本构关系对材料大变形和断裂机理研究具有重要意义。传统拉伸试验获取的材料真应力-真应变曲线在试样颈缩后不可测。借助可以精确测量三维变形的DIC(Digital image correlate) 技术和有限元分析技术(Finite element analysis),本文提出了基于漏斗试样拉伸试验获取材料全程单轴本构关系的新方法,即TF(Test and FEA)方法。该方法将TF方法获取的材料全程单轴应力应变关系曲线作为有限元软件中的材料本构关系对漏斗试样拉伸变形过程进行模拟,其模拟载荷-位移曲线、漏斗根部直径-位移曲线和漏斗变形轮廓线等均与试验结果吻合良好,试样表面模拟应变也与DIC测试结果吻合, 根据不同半径漏斗试样模拟获得的全程真应力-真应变曲线保持良好一致性。最后,还对试样颈缩断面的力学行为进行了讨论,并给出了304不锈钢、汽轮机叶片材料2Cr12Ni4Mo3VNBN和 1Gr12Ni3Mo2VN、汽轮机转子材料30Cr2Ni4MoV的全程单轴本构关系模型参数、破断应力和破断应变。     

The initiation and growth of adiabatic shear bands

A simple version of thermo/viscoplasticity theory is used to model the formation of adiabatic shear bands in high rate deformation of solids. The one dimensional shearing deformation of a finite slab is considered. For the constitutive assumptions made in this paper, homogeneous shearing produces a stress/strain response curve that always has a maximum when strain and rate hardening, plastic heating, and thermal softening are taken into account. Shear bands form if a perturbation is added to the homogeneous fields just before peak stress is obtained with these new fields being used as initial conditions. The resulting initial/boundary value problem is solved by the finite element method for one set of material parameters. The shear band grows slowly at first, then accelerates sharply, until finally the plastic strain rate in the center reaches a maximum, followed by a slow decline. Stress drops rapidly throughout the slab, and the central temperature increases rapidly as the peak in strain rate develops.     

A multi-scale computational model of crystal plasticity at submicron-to-nanometer scales

Plastic flow in crystal at submicron-to-nanometer scales involves many new interesting problems. In this paper, a unified computational model which directly combines 3D discrete dislocation dynamics (DDD) and continuum mechanics is developed to investigate the plastic behaviors at these scales. In this model, the discrete dislocation plasticity in a finite crystal is solved under a completed continuum mechanics framework: (1) an initial internal stress field is introduced to represent the preexisting stationary dislocations in the crystal; (2) the external boundary condition is handled by finite element method spontaneously; and (3) the constitutive relationship is based on the finite deformation theory of crystal plasticity, but the discrete plastic strains induced by the slip of the newly nucleated or propagating dislocations are calculated by dislocation dynamics methodology instead of phenomenological evolution equations used in conventional crystal plasticity. These discrete plastic strains are then localized to the continuum material points by a Burgers vector density function proposed by us. Various processes, such as loop dislocation evolution, dislocation junction formation etc., are simulated to verify the reliability of this computational model. Specifically, a uniaxial compression test for micro-pillars of Cu is simulated by this model to investigate the ‘dislocation starvation hardening’ observed in the recent experiment.     

高分子材料平面应变拉伸变形局部化

将基于应变软化玻璃状高分子材料微观特征建立的ＢＰＡ８－链分子网络模型引入ＵｐｄａｔｉｎｇＬａｇｒａｎｇｅ有限元方法，建立了适于变形局部化分析的大变形弹塑性有限元驱动应力法．在此基础上，数值模拟了初始各向同性高分子材料平面应变拉伸变形局部化的传播过程．探讨了ＢＰＡ模型对具有加工硬化特性的结晶性高分子材料变形分析的适应性；分析了局部化传播过程中颈缩截面的非均匀应力三轴效应；最后，讨论了网格尺寸以及初始几何不均匀性对颈缩扩散以及应力三轴效应的影响     

平面弹性力学问题的离散元法

根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.     

Strain gradient crystal plasticity analysis of a single crystal containing a cylindrical void

The effects of void size and hardening in a hexagonal close-packed single crystal containing a cylindrical void loaded by a far-field equibiaxial tensile stress under plane strain conditions are studied. The crystal has three in-plane slip systems oriented at the angle 60° with respect to one another. Finite element simulations are performed using a strain gradient crystal plasticity formulation with an intrinsic length scale parameter in a non-local strain gradient constitutive framework. For a vanishing length scale parameter the non-local formulation reduces to a local crystal plasticity formulation. The stress and deformation fields obtained with a local non-hardening constitutive formulation are compared to those obtained from a local hardening formulation and to those from a non-local formulation. Compared to the case of the non-hardening local constitutive formulation, it is shown that a local theory with hardening has only minor effects on the deformation field around the void, whereas a significant difference is obtained with the non-local constitutive relation. Finally, it is shown that the applied stress state required to activate plastic deformation at the void is up to three times higher for smaller void sizes than for larger void sizes in the non-local material.     

Crystal plasticity finite element modeling of mechanically induced martensitic transformation (MIMT) in metastable austenite

A new crystal plasticity model incorporating the mechanically induced martensitic transformation in metastable austenitic steel has been formulated and implemented into the finite element analysis. The kinetics of martensite transformation is modeled by taking into consideration of a nucleation-controlled phenomenon, where each potential martensitic variant based on Kurdjumov–Sachs (KS) relationship has different nucleation probability as a function of the interaction energy between externally applied stress and lattice deformation. Therefore, the transformed volume fractions are determined following selective variants given by the crystallographic orientation of austenitic matrix and applied stress in the frame of the crystal plasticity finite element. The developed finite element program is capable of considering the effect of volume change by the Bain deformation and the lattice-invariant shear during the martensitic transformation by effectively modifying the evolution of plastic deformation gradient of the conventional rate-dependent crystal plasticity finite element. The validation of the proposed model has been carried out by comparing with the experimentally measured data under simple loading conditions. Good agreements with the measurements for the stress–strain responses, transformed martensitic volume fractions and the influence of strain rate on the deformation behavior will enable the model to be promising for the future applications to the real forming process of the TRIP aided steel.     

The near tip stress and strain fields in cracked single crystals

The numerical analyses of stationary mathematically sharp Mode I crack in FCC and BCC crystals with elastic-ideally plastic (EIP) and fast hardening saturation (FHS) law are carried out in the present paper. From the calculated results, it is shown that: for the cases of small strain, EIP crystal cracks, the features of concentrated deformation patterns and the stress state in near-crack tip deformation fields are identical to the earlier analytical solutions, but along the angular sector boundaries, there exist narrow complex stress zones. The overall characteristics of deformation patterns for the cases of EIP and FHS are similar. The behaviours of crack tip opening can be characterized by crack-tip-opening-displacement (CTOD). For the case of FHS, finite deformation BCC crystal crack, our calculations are qualitatively in agreement with recent experimental observations. The project supported by National Natural Science Foundation of China