一个反映土体复杂加卸荷路径的弹黏塑性模型

基于Perzyna 过应力理论,构造了一个土体各向异性弹黏塑性本构模型. 模型引入Wheeler 旋转硬化法则,能够较好的描述土体的初始各向异性和应力诱发各向异性. 借助ABAQUS 软件中UMAT 子程序接口,模型采用隐式积分算法——图形返回算法实现. 通过对一组复杂加卸荷路径的三轴不排水剪切试验(HKMD) 及一组分级加载三轴不排水蠕变试验(Sachville clay) 的模拟,表明本模型能够合理反映土体的率效应、应力松弛、蠕变及土体各向异性现象,验证了模型的有效性.     

An anisotropic elastic–viscoplastic model for soft clays

Experimental evidences have shown deficiencies of the existing overstress and creep models for viscous behaviour of natural soft clay. The purpose of this paper is to develop a modelling method for viscous behaviour of soft clays without these deficiencies. A new anisotropic elastic–viscoplastic model is extended from overstress theory of Perzyna. A scaling function based on the experimental results of constant strain-rate oedometer tests is adopted, which allows viscoplastic strain-rate occurring whether the stress state is inside or outside of the yielding surface. The inherent and induced anisotropy is modelled using the formulations of yield surface with kinematic hardening and rotation (S-CLAY1). The parameter determination is straightforward and no additional experimental test is needed, compared to the Modified Cam Clay model. Parameters determined from two types of tests (i.e., the constant strain-rate oedometer test and the 24 h standard oedometer test) are examined. Experimental verifications are carried out using the constant strain-rate and creep tests on St. Herblain clay. All comparisons between predicted and measured results demonstrate that the proposed model can successfully reproduce the anisotropic and viscous behaviours of natural soft clays under different loading conditions.     

动力UH模型及其有限元应用

饱和砂土在循环载荷下具有复杂的应力应变关系, 通常表现出液化过程中的大变形以及往返活动性现象. 为简单有效地模拟上述特性, 在超固结UH模型的基础上, 将其扩展为可考虑砂土动力加载下的本构模型. 具体做法有3点: (1)改变屈服面椭圆长短轴之比, 将比值定义为反映应力诱导各向异性转轴斜率的函数; (2)引入旋转硬化规则, 用来反映应力诱导各向异性; (3)建立一个与旋转硬化规则以及临界状态特性相协调的统一硬化参数. 模型预测结果表明, 所提动力模型可简单、有效地用于砂土在动力载荷下应力应变关系的模拟. 最后将该动力UH模型嵌入到有限元软件中, 三维地基的动力加载模拟结果表明, 动力UH模型可方便地应用于岩土工程实践中.      

Elastoplastic model describing cyclic creeping behaviour of soft clays

The elastoplastic model in the deviatoric stress space is constructed to describe the cyclic undrained creeping behaviour of soft clays under cyclic stress by using Mises yield criterion and the concept of a field of hardening moduli. Furthermore, the effect of model parameters on the deformation is studied, and a method is given to determine quantitatively model parameters from results of cyclic triaxial tests of the saturated soft clay. The project supported by National Natural Science Foundation of China and Tianjin Twenty-one Century Youth Science Foundation     

Micromechanical analysis of the behavior of stiff clay

Cementations formed in geological timescale are observed in various stiff clays.A micromechanical stress strain model is developed for modeling the effect of cementation on the deformation behavior of stiff clay.The proposed approach considers explicitly cementations at intercluster contacts,which is different from conventional model.The concept of inter-cluster bonding is introduced to account for an additional cohesion in shear sliding and a higher yield stress in normal compression.A damage law for inter-cluster bonding is proposed at cluster contacts for the debonding process during mechanical loading.The model is used to simulate numerous stress-path tests on Vallericca stiff clay.The applicability of the present model is evaluated through comparisons between the predicted and the measured results.In order to explain the stress-induced anisotropy arising from externally applied load,the evolution of local stresses and local strains at inter-cluster planes are discussed.     

恒定主应力偏转角下黏土不排水剪切特性分析

天然土的结构性、各向异性对其强度有影响.进行固结压力100kPa、主应力轴恒定偏转角为0°,15°,30°,45°,60°的恒定主应力轴偏转角的固结不排水剪切试验,分析了不同恒定主应力轴偏转角对饱和软黏土固结不排水剪切强度影响;推导了总应力控制恒定主应力轴偏转角剪切试验总应力路径,编制了总应力控制剪切三维点积分程序,对比恒定主应力轴偏转角饱和软黏土固结不排水剪切强度计算值与试验结果,表明程序能很好地反应各向异性对土剪切强度影响.     

正常固结饱和黏性土三剪有限变形弹塑性本构模型研究

基于三剪统一强度准则,利用等量代换法和坐标平移法分别推导出正常固结饱和黏性土的三剪破坏应力比,并将其与修正剑桥模型相结合,得到三剪统一屈服面方程。为反映饱和黏性土的变形非线性及大变形特性,基于有限变形理论,建立了正常固结饱和黏性土的三剪统一有限变形等量代换法和坐标平移法2种弹塑性本构模型。为验证模型的适用性,取江西红黏土制备三种不同压实度的试样,在不同围压下进行常规三轴固结不排水和固结排水压缩试验。将试验数据与本文所提2种有限变形模型及相应的2种小变形模型计算结果进行了对比分析,结果表明,随着变形的发展,有限变形模型相对于小变形模型更接近试验结果,能较好地反映黏性土因高孔隙率（低初始压实度、小围压及压缩变形前期）而产生的大变形特性,虽然小变形模型的计算偏差会随着土样初始压实度和所施加围压的增大而减小,但有限变形模型对不同压实度和围压的计算偏差均相对较小,其中,等量代换法有限变形模型在初始压实度较大或者围压较高时所得计算偏差相对最小。对所提本构模型所做真三轴分析表明,中间主应力影响系数b和初始压实度对土体的强度和变形特性具有一定的影响。主应力差、孔隙水压力和体应变与b值呈正相关性,主应力差与初始压实度呈正相关性,孔隙水压力和体应变与初始压实度呈负相关性。     

Numerical simulation of plastic deformation and damage accumulation in structural materials under various low-cycle loading conditions

This paper presents the results of a finite-element study of elastic-plastic deformation and damage accumulation in structural materials under various cyclic loading conditions. Material behavior is described by the relations of damage mechanics using thermoplastic model which takes into account the plastic deformation of material under cyclic loading and the kinetic equations of the energy theory of damage accumulation. The basic laws of plastic deformation and development of damage in materials under hard, soft, symmetric, and asymmetric low-cycle loading are established.     

The general rheological model of paddy soils in south China

Analyses and experiments show that the Burgers model can be accepted as a general rheological model of paddy soils in South China, among which most are clays. Orthogonal tests and variance analyses show that the effect of soil conditions (mechanical composition, water content, bulk density) on rheological parameters is particularly significant, while that of loading conditions (shape and size of bearing surface, load) is practically not significant.     

A RHEOLOGICAL CONSOLIDATION MODEL FOR SOFT CLAY WITH CONSIDERATION OF TEMPERATURE AND THE ANALYTICAL

基于四元件流变模型,考虑温度影响,建立了宁波软黏土流变固结模型,并利用拉普拉斯变换得到瞬时加载条件下的考虑温度 影响的饱和软黏土流变固结解析解;利用不同温度下软黏土的流变固结试验结果,拟合获得了四元件流变模型参数,进行编程得到理论值, 并与试验值进行对比. 结果表明：该模型能较好反映宁波软黏土流变固结特性,计算结果与试验结果较为吻合.温度的升高导致渗透系数增大,在相同的时间内孔压消散 越快, 固结 越快.     

基于近场动力学的单晶冰弹性各向异性数值模拟方法

天然冰材料在变形与破坏行为上的各向异性特征是冰与结构相互作用中产生复杂载荷过程的关键诱因, 而天然冰各向异性的根源则在于单晶冰的各向异性. 目前, 学术界针对单晶冰各向异性的数值模拟方法研究仍较为缺乏. 为了准确再现天然冰材料的特殊力学性质, 本文基于近场动力学理论, 提出了一种单晶冰弹性各向异性的数值模拟方法. 该方法的核心思想是将单晶冰杨氏模量沿不同加载方向的变化规律引入到近场动力学力密度向量的影响函数中. 以前人实验测试得到的杨氏模量值为参考, 通过开展与C轴呈0°, 45°和90°三个加载方向的单晶冰单轴压缩数值模拟实验, 提出了针对该影响函数的修正和辅助参数标定方法, 最终在15°, 30°, 60°和75°等其他四个加载方向进行了验证. 结果表明: 本文提出的针对影响函数的修正与参数标定方法, 能够较为便捷地找到数值模型杨氏模量与参考杨氏模量相一致的影响函数最优解, 即本文提出的基于影响函数的近场动力学数值模拟方法, 能够合理、准确地模拟单晶冰的弹性各向异性行为. 本文研究成果可为后续多晶冰各向异性数值模拟方法的建立提供基础性参考.      

Dynamic stress equilibration in split Hopkinson pressure bar tests on soft materials

The condition of dynamic stress equilibrium is not satisfied automatically when a split Hopkinson pressure bar (SHPB) is employed to determine the dynamic properties of soft materials. In order to develop guidelines for the proper design of SHPB experiments under valid testing conditions, an integrated experimental/analytical study has been conducted to examine the process of dynamic stress equilibrium in a soft rubber specimen. Dynamic compressive experiments on a RTV 630 and an ethylene-propylene-diene monomer rubber with a SHPB modified for soft material testing were conducted to determine the effects of specimen thickness and loading rate on the stress equilibrating process. An analytical model was employed to analyze the equilibrating processes observed in experiments. It is found that the incident loading rate dominates the initial non-equilibrium stress state, and the specimen thickness mainly affects the dynamic stress equilibrium after the initial stage.     

A finite-strain model for anisotropic viscoplastic porous media: I – Theory

In this work, we propose an approximate homogenization-based constitutive model for estimating the effective response and associated microstructure evolution in viscoplastic (including ideally-plastic) porous media subjected to finite-strain loading conditions. The proposed model is based on the “second-order” nonlinear homogenization method, and is constructed in such a way as to reproduce exactly the behavior of a “composite-sphere assemblage” in the limit of hydrostatic loading and isotropic microstructure. However, the model is designed to hold for completely general three-dimensional loading conditions, leading to deformation-induced anisotropy, whose development in time is handled through evolution laws for the internal variables characterizing the instantaneous “ellipsoidal” state of the microstructure. In Part II of this study, results will be given for the instantaneous response and microstructure evolution in porous media for several representative loading conditions and microstructural configurations.     

A robust integration algorithm for implementing rate dependent crystal plasticity into explicit finite element method

Severe numerical instability in the integration of rate dependent crystal plasticity (RDCP) model is one of the main problems for implementing RDCP into finite element method (FEM), especially for simulating dynamic/transient forming process containing complicated contact conditions under large step length, large strain and high strain rate. In order to overcome the problem, an implicit model is deduced with the primary unknowns of shear strain increments of slip systems under the corotational coordinate system in the paper. The homotopy auto-changing continuation method combined with the Newton–Raphson (N–R) iteration is adopted. The subroutine VUMAT is developed for implementing RDCP model in ABAQUS/Explicit. Simulation results show that the algorithm is stable and accurate in 3D FE simulations on both dynamic simple loading and complicated loading process containing nonlinear contacts under the conditions of the maximal step length of 3.5 × 10⁻⁶ s, the maximal strain of 1.05, the maximal loading speed of 120 mm s⁻¹, and the minimal material rate sensitivity coefficient of 0.01. The predictions of the model on crystal behaviors of anisotropy, rate sensitivity and elasticity, as well as ear profiles in deep cup drawing are in agreement with experiments.     

A simple plasticity model for prediction of non-coaxial flow of sand

A bounding surface plasticity model for non-coaxiality, another aspect of anisotropic behavior of sands under rotation of principal stress axes; is developed in the critical state framework. Numerous experimental evidences exist that corroborate dependence of plastic shear strain rate direction on inherent fabric anisotropy. At first, general expressions for plastic strain rate with respect to possible emerge of non-coaxial flow are obtained. Consequently, using an anisotropy state parameter that is specially developed for this model and accounts for the interaction between imposed loading and soil fabric; effect of anisotropy on plastic flow direction is taken into account. Besides, novel circumstances are proposed for plastic modulus and dilatancy under rotation of principal stress axes. Finally, it is shown that the model is able to simulate successfully the non-coaxial behavior of sands subjected to principal stress axes rotation.     

An anisotropic fibre-network model for mechano-sorptive creep in paper

In this paper a simplified network model for mechano-sorptive creep is presented, which is a further development of an earlier paper [Strömbro, J., Gudmundson, P., 2008. Mechano-sorptive creep under compressive loading – a micromechanical model. International Journal of Solids and Structures 45 (9), 2420–2450.]. It is assumed that the anisotropic hygroexpansion of the fibres leads to large stresses at the fibre bonds when the moisture content changes. The resulting stress state will accelerate creep if the fibre material obeys a constitutive law that is non-linear. Fibre kinks are included in order to capture experimental observations of larger mechano-sorptive effects in compression than in tension. Moisture dependent material parameters and anisotropy in the fibre distribution have been introduced. Theoretical predictions based on the model are compared to experimental results for an anisotropic paper both under tensile and compressive loading at varying moisture content and it is found that the important features in the experiments are captured by the model. Different kinds of drying conditions have also been examined.     

晶粒数量对多晶集合体初始各向异性的影响

Taylor类多晶晶体粘塑性模型被用于研究晶粒数量对随机分布多晶体拉伸塑性各向异性的影响。分别沿包含不同晶粒数量的多晶集合体的各方向进行单向拉伸数值模拟实验，得到多晶集合体各方向在一定等效应变下的等效应力，并用云图和等高线表示在多晶体的参考球面上。定义了描述多晶集合体各向异性程度的参考指标。讨论了三种确定晶体随机取向的方法。计算结果表明：晶粒数量有限的多晶集合体的应力应变响应仍有一定的各向异性，且随着晶粒数量增多，多晶集合体的各向异性程度降低；就所包含晶粒数相同的多晶集合体来说，在确定晶粒随机取向时，选取不同的方法对它的各向异性程度也有一定的影响。     

Modeling multi-axial deformation of planar anisotropic elasto-plastic materials,part I: Theory

In sheet metal forming processes local material points can experience multi-axial and multi-path loadings. Under such loading conditions, conventional phenomenological material formulations are not capable to predict the deformation behavior within satisfying accuracy. While micro-mechanical models have significantly improved the understanding of the deformation processes under such conditions, these models require large sets of material data to describe the micromechanical evolution and quite enormous computation expenses for industrial applications. To reduce the drawbacks of phenomenological material models under the multi-path loadings a new anisotropic elasto-plastic material formulation is suggested. The model enables the anisotropic yield surface to grow (isotropic hardening), translate (kinematic hardening) and rotate (rotation of the anisotropy axes) with respect to the deformation, while the shape of the yield surface remains essentially unchanged.Essentially, the model is formulated on the basis of an Armstrong–Frederick type kinematic hardening, the plastic spin theory for the reorientation of the symmetry axes of the anisotropic yield function, and additional terms coupling these expressions. The capability of the model is illustrated with multi-path loading simulations in ‘tension-shear’ and ‘reverse-shear’ to assess its performance with ‘cross’ hardening and ‘Bauschinger’ effects.     

多种冲击载荷条件下的人体肋骨骨折有限元分析

旨在研究不同碰撞载荷条件下基于不同失效模型的人体肋骨骨折机理. 为此采用已验证的人体有限元胸部模型来分析人体肋骨骨折现象. 该模型基于人体解剖学结构,包含了人体胸椎、腰椎、肋骨、胸骨、肋间软骨、胸腹部器官和其他的软组织,定义的生物材料参数都基于已有的文献记载. 基于人体在损伤生物力学领域内一些较为典型的肋骨骨折失效模型,根据已发表文献中的人体标本实验载荷条件模拟了人体肋骨结构在不同冲击载荷下的骨折现象,并与这些实验结果进行了对比分析. 所引用的实验结果包括单根肋骨强度结构实验和人体胸部正面碰撞块冲击实验. 从文中有限元仿真分析的结果来看,针对不同的载荷条件,不同肋骨骨折失效模型的适用性各不相同. 该人体胸部有限元模型可用于车辆交通事故中冲击载荷条件下的人体肋骨损伤生物力学研究.      

基于ANSYS ACP的复合材料螺旋桨流固耦合计算方法

针对复合材料螺旋桨的流固耦合水弹性行为,为形成双向耦合的集成平台,对复合材料螺旋桨进行了流固耦合计算方法的研究。在ANSYS Workbench平台中,采用ACP模块对复合材料螺旋桨的建模进行了分析,介绍了复合材料螺旋桨的建模流程、流体计算模型和固体计算模型。通过静态加载试验对复合材料螺旋桨建模方法的可靠性进行了验证,计算值与试验值的偏差在5.1%以内,且反映出材料的各向异性。通过空泡水洞模型试验对复合材料螺旋桨流固耦合敞水性能计算方法的精确度进行了验证,在转速为12.7r/s和25r/s两个工况下计算的推力和扭矩与试验值吻合较好。验证了基于ANSYS ACP的复合材料螺旋桨稳态流固耦合计算方法合理可行,为复合材料螺旋桨的流固耦合振动特性分析提供了参考。