复杂应力状态下岩土材料非线性本构模型研究

基于统一强度理论,将洛德参数引入统一强度理论,推导出材料的统一强度参数,进而分析中间主应力以及主剪应力系数对材料统一强度参数的影响。在此基础上,对邓肯张双曲线模型进行改进,使得该模型能够反映复杂应力状态下的应力应变关系。采用粘土的平面应变试验进行验证,结果表明:在平面应变情况下,当b=0.4时,模型能够较好地反映粘土的应力应变关系,及其强度参数。分析了洛德参数以及中间主应力系数对模型的影响,进一步说明岩土材料存在主应力效应。该模型能够反映不同材料在复杂应力状态下的应力应变关系,有其更为广泛的适用性。     

压杆/试样表面接触变形对SHPB实验应变测量的影响

采用SHPB（split Hopkinson pressure bar）实验技术测量了3种不同尺寸纯铁试样的动态压缩应力应变关系,根据实验结果提出一个经验模型,定量分析了SHPB实验中压杆/试样表面接触变形对应变测量的影响。分析表明,在轴向应力平衡条件下,表面的接触变形对弹性段的应变测量影响显著;对塑性段应变测量的影响与试样的强化模量和长度相关,当试样强化模量较大而长度较小时,这种影响将不可忽略,可根据影响量的经验分析模型对应变进行修正。     

压缩-剪切复合应力下PMMA材料的静态力学特性研究

通过压缩具有一定倾斜角(0°,10°,15°,20°和25°)试件和双剪切模型试件,实现了单轴压缩、压缩-剪切复合应力以及纯剪切三种应力状态,得到PMMA(聚甲基丙烯酸甲酯)在相应应力状态下的应力-应变曲线,同时对不同应力状态下试件的破坏模式进行了分析。结果表明:在不同受力环境中材料的强度和破坏的机理不同;同单轴压缩状态下相比,材料在压缩-剪切复合应力状态下屈服极限、强度极限以及破坏应变均不同程度的增大,呈现明显的"剪切增强"现象。单轴压缩与压缩-剪切应力状态下试件的破坏模式均为在试件短对角面上出现明显的剪切屈服带,由应力分析得出试件剪应力在短对角面上达到最大,引起在此平面上分子链间滑动从而产生应变软化形成剪切屈服带;双剪切试件的破坏模式为与剪切面呈45°的斜面。     

显式模拟类橡胶材料Mullins效应滞回圈

通过显式、直接的方法提出一个多轴可压缩应变能函数,用来模拟类橡胶材料在加载——卸载作用下,由于Mullins效应而产生的应力——应变滞回圈. 本文的创新点在于将表征能量耗散的变量引入到应变能函数.新的弹性势具有以下两个特点:第一,在加载情况下,新引入的变量不会对弹性势产生任何影响,因此,只要给出合适的形函数显式表达,3个基准实验,包括单轴拉伸和压缩,等双轴拉伸和压缩,以及平面应变,都可精确模拟;第二,新引入的变量在卸载情况下将被激活.在不同的卸载应力下,变量将发生改变,从而影响弹性势,使其最终产生不同的应力——应变关系卸载曲线,与对应的加载曲线共同构成应力——应变滞回圈.通过对Mullins效应实验数据进行分析和研究,得出了卸载形函数在不同卸载应力下变化的规律,并预测不同卸载应力下的应力——应变关系.最后,我们将得到精确匹配实验数据的数值模拟结果,从而证明本文方法不仅可以精确匹配至少3个基准实验,还可以模拟和预测类橡胶材料在加载——卸载作用下由于Mullins效应而产生的滞回圈.      

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

Crack-tip fields in steady crack-growth with linear strain-hardening

Singular stress and strain fields are found at the tip of a crack growing steadily and quasi-statically into an elastic-plastic strain-hardening material. The material is characterized byJ₂ flow theory together with a bilinear effective stress-strain curve. The cases of anti-plane shear, plane stress and plane strain are each considered. Numerical results are given for the order of the singularity, details of the stress and strain-rate fields, and the near-tip regions of plastic loading and elastic unloading.     

Size dependent response of large shell elements under in-plane tensile loading

Large-scale thin-walled structures with a low weight-to-stiffness ratio provide the means for cost and energy efficiency in structural design. However, the design of such structures for crash and impact resistance requires reliable FE simulations. Large shell elements are used in those simulations. Simulations require the knowledge of the true stress–strain response of the material until fracture initiation. Because of the size effects, local material relation determined with experiments is not applicable to large shell elements. Therefore, a numerical method is outlined to determine the effect of element size on the macroscopic response of large structural shell elements until fracture initiation. Macroscopic response is determined by introducing averaging unit into the numerical model over which volume averaged equivalent stress and plastic strain are evaluated. Three different stress states are considered in this investigation: uniaxial, plane strain and equi-biaxial tension. The results demonstrate that fracture strain is highly sensitive to size effects in uniaxial tension whereas in plane strain or equi-biaxial tension size effects are much weaker. In uniaxial and plane strain tension the fracture strain for large shell elements approaches the Swift diffuse necking condition.     

Crack nucleation near stress concentrators in quasi-brittle materials

The results of combined tension (compression) and shear experiments with plexiglass specimens are used to construct a Coulomb-Mohr-type strength curve. Fracture experiments were performed with plexiglass square plates with internal through cuts. The specimens were subjected to compression on a Zwick/Roell testing machine until cracks appeared. In the process of testing, nucleation of symmetric cracks was observed, which propagated as the load was further increased. The fracture character (tensile or shear) could not be decided until a subsequent numerical analysis of the plate stress-strain state was performed by the finite elementmethod. The stress concentration loci were found to coincide with the crack nucleation loci. In the plane (σ, τ ), Mohr’s circles were constructed for the stress states at the stress concentration points. If the point at which Mohr’s curve touches the limit curve is known, then one can determine the plane on which the normal and tangential stresses attain critical values, and hence one can determine the crack propagation direction. The experimental results and numerical solutions were found to be in good agreement.     

双向受压状态下的混凝土率效应试验研究

利用三峡大学自行研制的大型多功能三轴材料试验机,完成了4种侧应力等级和3个数量级加载速率的C30混凝土立方体试件双向受压试验,系统地探讨了不同侧应力和不同加载速率下混凝土的强度特性及变形特性。研究表明,侧应力和应变速率都能提高混凝土的应变峰值,但侧应力对混凝土的应变峰值的影响较应变速率的影响大;随着侧应力的增加,混凝土的率敏感性逐渐降低;在有侧应力的情况下,随应变速率的增加,混凝土弹性模量变化逐渐减小;随应变速率的增加,应力应变曲线上升段更加陡峭,曲线更加饱满,混凝土剩余强度越来越高;双轴受压情况下,混凝土的应力应变曲线下降段趋势比无侧应力状态下较平缓,随侧应力的增加下降段曲线越来越平缓,混凝土的应力应变曲线更加饱满,峰值点无明显尖峰出现,混凝土剩余强度越来越高。     

On the role of strain and stress state in ductile failure

Using plane strain and axisymmetric notched tensile specimens in combination with finite deformation stress analysis, the strain to initiate failure in a range of structural steels is shown to be co-related by the state of stress for both axisymmetric and plane states of strain. The implications of this result for ductile failure terminated by flow localisation is discussed in the light of the theoretical work on localised flow.     

基于PVDF薄膜传感器的猪肝动态压缩力学性能测试

瞬态冲击载荷作用下肝脏的力学响应是损伤生物力学的重要研究内容。本文提出了一种可用于软组织动态压缩力学特性测试的改进SHPB(分离式霍普金森压杆,Split Hopkinson Pressure Bar)方法。该方法采用PVDF(聚偏氟乙烯,Polyvinylidene Fluor)压电薄膜传感器测量实验过程中试件两端面的受力,以此来计算试件的应力,从而无需测量透射杆上的微弱透射信号。猪肝试样前后端面的PVDF压电信号对比表明,加载过程中试样达到了动态应力平衡状态。试样动态压缩中的惯性效应主要在加载的初始阶段对透射应力信号造成较大影响,在大变形阶段惯性效应引起的轴向应力较小。利用此方法对猪肝组织进行三种高应变率(1800s^-1,2500s^-1,3500s^-1)的动态压缩实验,并采用基于真实应变的惯性效应公式对实验数据进行修正计算。结果表明:猪肝组织在准静态与高应变率时的应力应变曲线都呈现出凹向上的非线性特征,即曲线初始阶段应力增长较缓慢,当应变达到15%后应力值则迅速增大;猪肝组织也具有明显的应变率效应,即随着应变率的增加,应力应变曲线的整体应力值也随之增大。最后,采用黏超弹性本构模型描述了猪肝组织的动态应力应变曲线。     

韧性材料颈缩主应力面上应变分布假定的研究

为了深入研究韧性材料主应力面的应变分布情况，讨论Bridgman关于主应力面上应变均布假定的精度．本文首先结合断口学基本理论与在室温条件下拉伸断裂的Q235钢、40Cr钢和45#钢断口电子显微镜扫描结果，从微观角度分析其韧窝分布后发现应变均布假定存在一定的误差；然后基于体积守恒和Aramis三维应变观测系统测得的试验结果，将计算值和试验值进行对比分析发现，三种材料中除含碳量较高的45#钢外，应变均布假定并不适用于分析颈缩主应力面上应变分布．本文的研究成果可为进一步研究韧性材料颈缩大塑性阶段的应力、应变本构关系提供一定的参考．     

Discretization of the plane problem for a cracked body with nonlinear stress–strain diagram under tension

The plane problem for a cracked body with a piecewise-linear stress–strain diagram under tension is reduced by the Fourier transformation to a system of nonlinear algebraic equations. The system is numerically solved for plane strain and stress states of a perfect elastoplastic material to study plastic zones, stress and strain distributions, and displacements of crack faces     

非局部平面应变和平面应力问题界定及其精确性讨论

在非局部弹性理论框架下对平面应变和平面应力状态重新界定.首先,分别在其相应简化假设下推导控制方程,并与经典局部情况进行比较.然后,引入变形协调条件对两类非局部平面问题的精确性进行讨论.其中,对于非局部平面应力状态,通过应变协调方程的Fourier变换形式来进行研究,使问题得以简化.通过以上分析,最终得到一些有价值的结论.     

Q235钢缺口试样拉压低循环疲劳的实验研究

以Q235钢制U型缺口板试样为研究对象,用有限元方法计算其缺口根部等效应变幅对应的试样标距段位移,以此控制试验机进行拉压循环疲劳试验。然后用局部应力应变法对试验测得的寿命结果进行分析。结果表明：无论用有限元还是修正Neuber公式计算缺口根部的应力应变,局部应力应变法的疲劳寿命评估只适用于缺口半径较大的试样;对缺口半径较小试样的估计寿命明显低于实测值,且有限元法比修正Neuber法更保守。进而又对试样缺口区域应变梯度的影响进行了探讨：参照有限元计算的应变梯度,利用Taylor模型估算了缺口根部的屈服应力和流动应力;在此基础上重新计算应变分布并估计试样的疲劳寿命,结果证实考虑应变梯度影响可改善缺口试样的疲劳寿命估计。     

Quasiconvex envelopes of energies for nematic elastomers in the small strain regime and applications

We provide some explicit formulas for the quasiconvex envelope of energy densities for nematic elastomers in the small strain regime and plane strain conditions. We then demonstrate their use as a powerful tool for the interpretation of mechanical experiments. Analytical formulas characterizing the stress-strain response in pure shear are derived, providing an easily testable benchmark for future numerical and experimental investigations on the mechanics of nematic elastomers.     

A damage model for crack initiation and growth

A damage accumulation model is presented for the study of the problem of crack initiation and stable growth in an elastic-plastic material. A centre-cracked specimen subjected to a uniform stress perpendicular to the crack plane is considered. A coupled stress and failure analysis is performed by using a finite element computer program based on J₂-plasticity theory in conjunction with the strain energy density theory. After initial yielding, each material element follows a different equivalent uniaxial stress-strain behavior depending on the amount of energy dissipation by permanent deformation. A host of uniaxial stress-strain curves constituting parts of the same stress-strain curve were assigned to material elements for each increment of loading. The path-dependent nature of the onset of crack initiation and growth was revealed. The proposed model predicts faster crack growth rates than those obtained on the basis of a single uniaxial stress-strain curve and is closer to experimental observation.     

不同加载状态下TA2钛合金绝热剪切破坏响应特性

一般认为绝热剪切现象在宏观上表现为材料动态本构失稳,即热软化大于应变硬化.本文采用帽型受迫剪切试样研究TA2钛合金的动态力学特性和本构失稳过程.首先对剪切区加载应力状态进行理论和数值分析,通过合理设计帽型试样,剪切区变形可近似按剪切状态处理;结合二维数字图像相关法(two-dimensional digital image correlation,DIC-2D)直接测试试样剪切区应变演化,给出帽型受迫剪切实验的等效应力-应变响应曲线.进一步,利用Hopkinson压杆对TA2钛合金开展动态压缩及帽型剪切对比试验研究,比较压缩、剪切试验得到的等效应力-应变曲线,采用"冻结"试样方法分析试样中绝热剪切局域化演化过程,探讨不同加载状态下TA2钛合金的绝热剪切破坏现象及其动态力学响应特性.实验结果表明,在塑性变形初始阶段,动态压缩及剪切加载下的等效应力-应变曲线符合较好,但随塑性损伤发展及绝热剪切带形成,两者出现分离,表明损伤及绝热剪切演化过程与应力状态相关.剪切试样实验得到的本构"软化"特性能够反映绝热剪切带起始、破坏演化过程的力学响应特性,而在动态压缩实验中,即使试样中已出现双锥形的绝热剪切带及局部裂纹分布,其表观等效应力-应变曲线并不出现软化特征,动态压缩实验无法得到关于绝热剪切起始、发展以及破坏的本构软化响应特性.     

Two-dimensional analysis of the split hopkinson pressure bar system

The split Hopkinson pressure bar is widely used to measure the dynamic properties of solid materials. This paper presents the results of the first comprehensive two-dimensional numerical analysis of the technique, and quantitatively describes the effects of realistic friction and of variations in both the specimen geometry and the imposed strain-rate on the validity of the assumptions used in analyzing experimental data. A two-dimensional axisymmetric numerical analysis is used to compute all components of the stress, strain and strain-rate tensors at each mesh point within the specimen and the elastic bars. The calculated response of the pressure bars is used to reconstruct the stress-strain behavior of the specimen and this is compared to both the input stress-strain curve and the actual calculated stress-strain states in the specimen. Thus, the validity of the assumptions and the corrections used in the analysis of the data is determined.Inertia and friction between the specimen and the elastic bars affect the response of the specimen differently for different length-to-diameter ratios. Inertia effects produce stress waves propagating radially and axially in the specimen and may result in an oscillating reconstructed stress-strain curve. If the ends of the specimen are well lubricated and care is taken to minimize the effects of inertia, the reconstructed stress-strain curve agrees with the input. However, serious stress and strain nonuniformity exists when the ends are not lubricated and this results in a reconstructed stress-strain curve where, for any given strain, the stress magnitude is larger than the correct value. A comparison of the calculations with experiment shows excellent agreement for various interface conditions. Finally, the inertia correction of E.D.H. Davies and S.C. Hunter (1963) is found to be reasonable.