Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

外包钢-砼组合梁纯扭作用下抗扭刚度的实验及数值模拟

为了研究新型外包钢-砼T形截面组合梁在纯扭作用下的变形性能,设计了5根不同配箍率的的足尺悬臂组合梁。通过对5根悬臂梁的抗扭性能的实验研究,得到了组合梁的扭矩-扭率关系曲线。利用有限元分析软件ANSYS,对组合梁的抗扭性能进行了非线性有限元分析,得到了混凝土与外包钢在极限阶段的应力云图。根据实验以及有限元结果分析了组合梁在整个加载过程中扭转刚度的变化。基于现行砼结构设计规范,提出了组合梁从开裂到极限阶段抗扭刚度的计算公式,可供组合梁受扭设计参考。把有限元模型和公式的计算结果与实验结果进行比较,三者吻合较好。     

THE USE OF VISCO-ELASTOPLASTIC DAMAGE CONSTITUTIVE MODEL TO SIMULATE NONLINEAR BEHAVIOR OF CONCRETE

A visco-elastoplastic damage constitutive model is proposed for simulating nonlinear behavior of concrete. Based on traditional plastic theory, the irreversible deformation is simulated in effective stress space. In order to reflect different stiffness degradation mechanism of concrete under tensile and compressive loading conditions, both tensile and compressive damage variables are introduced, and then on the basis of energy release rate, the model is firmly derived within the concept of irreversible thermodynamics. The rate-dependent model is considered by introducing viscous regularization into the inelastic strain and damage variable, and combined with an additional elastic condition. Fully implicit backward-Euler algorithm is used to perform constitutive integration. Results of numerical examples using the proposed model agree well with test results for specimens under uniaxial tension and compression, biaxial loading and triaxial loading. Failure processes of single-edge-notched (SEN) beam and double-edge-notched (DEN) specimen are also simulated to further validate the proposed model.     

高应力高水压下砂岩三轴蠕变特性试验研究

对砂岩进行高围压高水压条件下的三轴压缩蠕变试验。试验表明,在整个加载过程中,孔隙水压力主要起到增强轴向变形和横向变形的作用;但在加载的初始阶段,孔隙水压力在一定程度上抑制了轴向变形。当应力水平高于屈服应力时,横向蠕变速率明显大于轴向蠕变速率,且横向蠕变率先进入加速蠕变阶段。本文提出一个新的非线性黏性元件,并引入一个能判断是否进入加速蠕变阶段的计时器元件,组建一个非线性黏塑性加速蠕变启动模型,将该黏塑性模型与Burgers模型串联,构建一个新的非线性黏弹塑性蠕变模型,推导了该模型在常规三轴应力状态下的本构方程。基于试验结果,通过对非线性优化算法(BFGS)的Matlab编程,实现对本文提出蠕变模型的参数识别,识别效果比较理想。对比试验曲线与拟合曲线,二者相当吻合,验证了新提出的非线性黏弹塑性蠕变模型的正确性。     

A lattice deformation based model of metallic lattice sandwich plates subjected to impulsive loading

A lattice structure deformation mechanism based theoretical model is developed to predict the dynamic response of square lattice sandwich plates under impulsive loading. The analytical model is established on the basis of the three-stage framework proposed by Fleck and Deshpande (2004). In the first stage, the impulse transmitted from air shock loading to the sandwich plates by fluid-structure interaction is analytically calculated. The lattice core suffers non-uniform compression in the second stage due to the clamped boundary conditions. The structure deformation mechanism is introduced in the lattice core compression and the analytical nominal stress–strain curve of core compression accords well with previous experimental results. In the final stage, the sandwich plate is analyzed as a continuum plate with non-uniform thickness deduced by inconsistent deformation of the front and back sheets.The experiment results of square metallic sandwich plates with tetrahedral lattice core are presented and compared with analytical prediction to validate the theoretical model. Good agreements are found between the predicted and testing results for both the impulse transmitted to the sandwich plates and the maximum deflection of the back face sheet.     

复杂加载下混凝土的弹塑性本构模型

混凝土材料在不同应力路径下或复杂加载条件下会表现出差异性显著的应力应变关系,在小幅循环加载条件下,其应力应变关系会表现出类似于弹性变形的滞回曲线.在不同应力水平下,混凝土的应力应变关系以及破坏特性都具有静水压力相关特点,即随着静水压力增大,各向异性强度特性弱化.此外,混凝土受压及受拉破坏机理不同,因而对应于混凝土硬化损伤亦有不同,即可分为受压硬化损伤,受拉硬化损伤及两者的混合硬化损伤类型.基于Hsieh模型,对该模型进行了三点改进.（1）针对小幅循环加载下混凝土无塑性变形的试验规律,而模型中在应力水平较低的循环加载条件下始终存在塑性变形的预测问题,采用在边界面模型框架下,设置了应力空间的弹性域,初始屈服面与后续临界状态屈服面几何相似的假定.（2）基于广义非线性强度准则将原模型采用变换应力方法将其推广为三维弹塑性本构模型,采用变换后模型可合理的考虑不同应力路径对于子午面以及偏平面上静水压力效应形成的影响,并避免了边界面应力点奇异问题.（3）分别对拉压两种加载损伤模式建议了相应的硬化参数表达式,可分别用于描述上述加载中产生的应变软化及强度退化行为.基于多种加载路径模拟表明:所建立的三维弹塑性本构模型可合理地用于描述混凝土的一般应力应变关系特性.      

受拉钢筋混凝土构件破坏过程的数值模拟

采用三维材料破坏过程分析MFPA3D系统,对钢筋混凝土构件轴心受拉条件下的受力、变形与内部裂纹萌生、扩展及最终破坏全过程进行了数值试验研究。数值模型中引入统计分布函数反映了混凝土的非均匀性影响,并采用具有残余强度的弹性损伤本构模型及其破坏单元材料性质退化方法,利用位移加载方式对钢筋混凝土构件实施拉伸加载。通过对钢筋、素混凝土方形体以及钢筋混凝土方形柱体构件在拉伸作用下破坏过程的数值试验,分析了钢筋与混凝土两种材料之间的相互作用、约束机理和破坏机理。数值试验成果对于深入了解钢筋和混凝土的联合受力规律和钢筋在开裂前后对整体钢筋混凝土结构的作用机制有参考价值。     

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

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

T形短肢剪力墙静力性能有限元仿真分析

利用有限元分析软件ANSYS,首先采用其中三维实体单元SOLID65建立了T形短肢剪力墙有限元分析模型,从弹性到混凝土开裂直至破坏的全过程进行了仿真试验分析。分析了影响短肢剪力墙受力的几种因素:混凝土强度等级、配筋率、轴压比、墙肢截面高厚比对短肢剪力墙承载能力、变形能力及延性的影响,剖析了短肢剪力墙破坏过程及其原因。比较真实的反映了短肢剪力墙在轴压力和逐步加载侧向力共同作用下的响应。试验结果表明:增加混凝土等级和轴压比能提高试件的开裂、屈服和极限荷载,但应综合考虑其与变形能力、延性的关系。截面配筋率具有其特殊性,配筋率在1.4%¹.6%之间时试件的承载能力、变形能力及延性较好。墙肢截面高厚比是不稳定因素但在高厚比为6.5⁷.1时,延性及变形能力较强。     

弹体贯穿混凝土数值模拟的改进材料模型

研究混凝土结构在冲击载荷下的力学特性对武器以及防护结构的设计和评估具有重要意义,而合适的材料模型可以更准确地预测混凝土结构的力学行为和破坏模式。因此,本文中提出了一种改进的混凝土塑性损伤材料模型来描述其在冲击载荷下的力学响应。该改进模型考虑了压力-体积应变关系、应变率效应、洛德角效应和塑性损伤累积对混凝土材料力学特性的影响,并引入了一个与损伤相关的硬化/软化函数来描述压缩状态下的应变硬化和软化行为。随后,通过对3个独立的强度面进行线性插值得到了该改进模型的破坏强度面,并采用部分关联流动法则考虑了混凝土材料的体积膨胀特性。最后,开展了单个单元在不同加载条件下和弹体贯穿钢筋混凝土靶的数值模拟,验证了该改进模型的可行性、准确性以及预测性能提升。     

Analytical solutions for the time-dependent behaviour of composite beams with partial interaction

This paper presents a generic modelling for the time-dependent analysis of composite steel–concrete beams with partial shear interaction that occurs due to the deformation of the shear connection. The time effects considered in this modelling are those that arise from shrinkage and creep deformations of the concrete slab, and these effects are modelled using algebraic representations such as those of the age-adjusted effective modulus method (AEMM) and the mean stress method (MS), which are viscoelastic models for concrete deformation that can be stated algebraically. The generic model lends itself to closed form solutions for the analysis of composite beams subjected to a generic applied loading under a variety of end conditions. In this paper, the generic model is applied for the time-dependent analysis of composite beams that are simply supported and encastré, and to a propped cantilever, that are subjected to uniformly distributed loading and shrinkage deformations. Various representations of the structural behaviour of these beams are given in closed form which can also be used to benchmark available modelling techniques, i.e. finite element and finite difference formulations, which require a spatial discretisation to be specified as well as the time discretisation to perform a time analysis.     

An experimental study of inhomogeneous cyclic plastic deformation of 1045 steel under multiaxial cyclic loading

An experimental study was conducted on the inhomogeneous cyclic plastic deformation of 1045 steel under multiaxial cyclic loading. Thin-walled tubular specimens were used and small strain gages were bonded on the specimen surface to characterize the local deformation. The controlled loading paths included cyclic tension–compression, cyclic torsion, proportional axial-torsion, 90°-out-of-phase axial-torsion, and fully reversed torsion with a constant axial stress. The maximum stress in each experiment was lower than the lower yield stress of the material. It was found that the cyclic plastic deformation within the gage section of the specimen under multiaxial stress state followed the three-stage process that was observed from uniaxial loading, namely, incubation, propagation, and saturation. The plastic deformation was significantly inhomogeneous during the propagation stage, and the inhomogeneity continued through the saturation stage. The duration of each stage and the saturated strains were dependent on the cyclic stress amplitude and the loading path. Multiaxial stress state reduced the incubation stage. With identical equivalent stress magnitude, the nonproportional loading path resulted in the shortest incubation and propagation stages, and the saturated equivalent plastic strain magnitude was the smallest. Although the deformation over the gage section was inhomogeneous, the plastic deformation in a given local area was found to be practically isotropic.     

混凝土多尺度应力响应方程及其数值模拟

针对混凝土的多相多尺度材料组成特征及其复杂力学响应问题, 首先, 根据混凝土中各组成材料的几何特征, 将C-S-H凝胶、硬化水泥浆体、砂浆及混凝土细观组成分别视为纳观、微观、亚细观和细观尺度上的复合材料, 并利用颗粒空间堆积方法, 重构了混凝土各尺度复合材料的简化几何模型; 其次, 基于重构的几何模型和等效夹杂理论, 通过等效刚度的升阶计算和应力响应的降阶计算, 建立各尺度复合材料应力响应之间的过渡关系, 推导混凝土多尺度应力响应方程, 并编制相应的计算程序; 最后, 以单轴压缩载荷作用为例, 数值计算载荷作用下混凝土各尺度复合材料中的应力响应, 分析骨料空间位置和相互作用以及水化产物刚度、几何形状和空间取向对其应力响应的影响规律. 结果表明, 单轴压缩载荷作用下, 混凝土细观组成中的应力分布并不均匀; 骨料颗粒之间的距离影响到混凝土中的应力分布, 其有效影响范围约为骨料粒径的6倍; 水泥水化产物的刚度、几何形状和空间取向是影响其应力分布的重要因素, 刚度越大, 所受应力越大, 与载荷作用方向的夹角越小, 长椭球形水化产物沿载荷作用方向的应力越大, 扁椭球形水化产物与之相反.      

Irreversible deformation with subsequent unloading of a spherical viscoelastoplastic layer

Analytic solutions are obtained for a sequence of one-dimensional quasistatic problems describing viscoelastic deformation processes in the material of a hollow ball and the plastic flow nucleation and evolution processes occurring in the ball as the pressure on the outer boundary increases. The unloading process under slow removal of the loading pressure is considered as well. The stress fields and the elastic and plastic strain fields in the spherical layer material, the law of motion of the elastoplastic boundary, and the residual stress level and distribution are computed. It is assumed that at the stage preceding the plastic flow the material obeys the viscoelastic Voigt model and the loading surface is determined by the von Mises plastic flow condition.     

Rate (time)-dependent deformation behavior: an overview of some properties of metals and solid polymers

The inelastic deformation behaviors of metals and polymers are discussed with the aim of finding a common base that would simplify academic and engineering analyses. Only monotonic loading conditions at room temperature are considered. For loading at different rates, nonlinear relations between loading rate and stress level, creep stress level and creep strain, and relaxation rate and stress were common to both type of materials. There are, of course, significant differences in elastic properties, strength levels and the strains involved. Special properties such as relaxation behaviors and creep anomalies can be qualitatively and quantitatively reproduced by the state variable model VBO (viscoplasticity theory based on overstress). Since experimental investigations typically concentrate on one particular aspect of inelastic deformation behavior such as creep or strain-rate dependence, it is often difficult to gather a comprehensive data set for a given material. In spite of this, considerable similitude in the deformation behavior of metals and polymers in various test conditions has nevertheless been established.     

退火和调质42CrMo合金钢单轴循环塑性行为的实验研究

在室温下对退火和调质42CrMo合金钢进行了单轴应变控制和应力控制的系统循环实验,并对它们的应变循环和应力循环特性进行比较.揭示和分析了应变幅值、平均应变及其历史对材料应变循环特性的影响以及应力幅值、平均应力及其历史对棘轮行为的影响.讨论了应变循环和应力循环间的交互作用以及不同热处理工艺下材料循环变形行为间的区别.研究发现材料的热处理工艺、平均应力和应力幅值及其历史对材料的棘轮行为都有很大的影响.得到了一些有助于进行合理本构描述的结果.     

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

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

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

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

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

Coupled viscoplasticity damage constitutive model for concrete materials

A coupled viscoplasticity damage constitutive model for concrete materials is developed within the framework of irreversible thermodynamics.Simultaneously the Helmholtz free energy function and a non-associated flow potential function are given, which include the internal variables of kinematic hardening,isotropic hardening and dam- age.Results from the numerical simulation show that the model presented can describe the deformation properties of the concrete without the formal hypotheses of yield criterion and failure criteria,such as the volume dilatancy under the compression,strain-rate sen- sitivity,stiffness degradation and stress-softening behavior beyond the peak stress which are brought by damages and fractures.Moreover,we could benefit from the application of the finite element method based on this model under complex loading because of not having to choose different constitutive models based on the deformation level.