结构与土接触面的弹黏塑性本构模型

基于Goodman接触面单元,在Clough--Duncan双曲线关系式的基础上改进了三元件宾厄姆体,建立了结构与土接触面的弹黏塑本构模型,使用ABAQUS的FRIC子程序对该模型进行了二次开发.通过改装的应力控制式结构-土接触面剪切试验仪进行了接触面长期剪切试验,验证了该接触面本构模型.研究结果表明:结构与土接触面的弹黏塑本构模型能较好地模拟接触面的力学特性和时间效应;正应力对结构与土接触面的力学行为影响较大,正应力越大,接触面的时间效应越显著.     

弹塑性微凸体侧向接触相互作用能耗

传统的结合面研究多基于光滑刚性平面与等效粗糙表面接触假设,忽略了结合面上微凸体侧向接触及相邻微凸体之间的相互作用,这导致理论模型与实际结合面存在较大出入.针对承受法向静、动态力的机械结合面,从微观上研究了微凸体侧向接触及相互作用的接触能耗.将法向静、动态力分解为法向分力和切向分力,获取弹性/弹塑性/塑性阶段考虑微凸体侧接触及相互作用的加、卸载法向分力-变形和切向分力-位移的关系.通过力的合成定理,从而获取加、卸载法向合力与总变形之间的关系,由于法向分力产生的塑性变形及切向分力产生的摩擦,导致加载、卸载法向合力-总变形曲线存在迟滞回线.通过对一个加、卸载周期内的法向合力-总变形曲线积分,获得一个周期的微凸体接触能耗,包括应变能耗及摩擦能耗.仿真分析表明:微凸体在3个阶段的能耗均随变形的增大而非线性增大.微凸体侧向接触角度越大,能耗越大,且在弹性阶段最为明显.在弹性阶段,仅存在侧向的摩擦能耗,故结合面在低载荷作用下必须采用双粗糙表面假设.在塑性阶段,由于微凸体接触能耗为应变能耗,且接触角对其能耗影响甚微,故结合面在大载荷作用下可采用单平面假设对其进行研究.相对于KE和Etsion模型,本文提出的模型与Bartier的实验结果更吻合.     

钙质砂与结构接触面的本构模型研究

基于接触面的宏、细观物理特征,建立了单调加载条件下钙质砂与结构接触面的弹塑性增量本构关系。从接触面的宏观条件上考虑,该模型将弹性模量取为法向压力的指数函数,采用非关联流动法则和Mohr-columb屈服函数,以及切向塑性功为硬化参量,适用于多数接触摩擦问题。在细观上将滑动面抽象为锯齿面,同时将摩擦系数取为塑性功的双曲线函数,以考虑钙质砂颗粒破碎对接触面力学特性的影响。模型概念简单、参数较少,通过理论计算与钙质砂拉拔试验结果比较,显示了模型的合理性。     

Duncan—Chang土体非线性弹性模型向弹塑性模型的转化

本文中从传统的弹塑性概念出发，将Ｄｕｎｃａｎ－Ｃｈａｎｇ土体双曲线应力－应变的非线性关系转化为一种新型且简易的弹塑性模型，以１９８号油罐工程为算例阐明其应用，并验证了本模型的可靠性。     

A work softening joint element used in dynamic analysis of soil-structure

This paper introduces an elasto-plastic joint element characterised by strain hardening and softening in the analysis of dynamic soil-structure interaction. The phenomena of separation and sliding on the contact surface between soil and structure can be better simulated and the process can also be described. The interaction problems in a typical soil-structure system are analyzed in terms of elasto-plastic joint element as well as elastic ones. The results show that the elasto-plastic joint element is much better than the elastic one in modelling, especially in that the relative displacements accross the joint element can be much greater than that of the elastic case. Separation and sliding are not only related to the coefficient of friction and cohesion but also to their changes with plastic volumetric strain. The project is supported by the Sciences Foundation of the Chinese Academy of Sciences and the Disaster Provention Institute of Kyoto University (Japan).     

石墨烯与柔性基底界面载荷传递的拉曼光谱实验

本文研究了CVD制备的大尺寸石墨烯与柔性PET基底在拉伸变形过程中切向界面载荷传递的问题,采用原位拉曼光谱实验给出了加载过程中石墨烯的正应变、正应力以及界面切应力的分布曲线。分析表明,石墨烯与PET基底间的载荷传递存在四个阶段,分别是初始阶段、粘附阶段、滑移阶段和界面脱粘破坏阶段。在此基础上,本文对50μm、140μm、270μm和600μm四种尺寸石墨烯试件的界面力学性能进行测量,得到了不同尺寸石墨烯试件的界面力学性能参数,并初步给出了基底变形引起的石墨烯切向界面粘接能的变化,同时分析了试件尺寸对石墨烯界面力学性能的影响。实验结果表明,石墨烯材料和柔性基底最大切应力与临界脱粘切向界面粘接能等界面强度指标受到尺寸的显著影响,尺寸越小切向界面强度越高,反之,尺寸越大则越低。     

Constitutive modeling of soil-structure interface through the concept of critical state soil mechanics

The behavior of soil-structure interface can be crucial to the overall response of a soil-structure system. The numerical simulation of soil-structure interaction problem requires proper modeling of the interface. The similarity between the behavior of soil and interface is first analyzed in the present paper. With this similarity, the concept of critical state soil mechanics (CSSM), which has been successfully used in the modeling of soil behavior, is used to develop a constitutive model for soil-structure interface in the framework of generalized plasticity. The model is capable of modeling strain hardening, softening, normal dilatancy and stress-path dependency of interface between sandy soil and structures during shearing. The effects of normal pressure as well as density of sand are captured in the model. The performance of the model is verified with various experimental results. The unified modeling of the behavior of interfaces with different roughness, different density of soil and different normal pressures using the concept of CSSM is also successfully attempted.     

Numerical simulation of elasto-plastic deformation of composites: evolution of stress microfields and implications for homogenization models

The deformation of a composite made up of a random and homogeneous dispersion of elastic spheres in an elasto-plastic matrix was simulated by the finite element analysis of three-dimensional multiparticle cubic cells with periodic boundary conditions. “Exact” results (to a few percent) in tension and shear were determined by averaging 12 stress-strain curves obtained from cells containing 30 spheres, and they were compared with the predictions of secant homogenization models. In addition, the numerical simulations supplied detailed information of the stress microfields, which was used to ascertain the accuracy and the limitations of the homogenization models to include the nonlinear deformation of the matrix. It was found that secant approximations based on the volume-averaged second-order moment of the matrix stress tensor, combined with a highly accurate linear homogenization model, provided excellent predictions of the composite response when the matrix strain hardening rate was high. This was not the case, however, in composites which exhibited marked plastic strain localization in the matrix. The analysis of the evolution of the matrix stresses revealed that better predictions of the composite behavior can be obtained with new homogenization models which capture the essential differences in the stress carried by the elastic and plastic regions in the matrix at the onset of plastic deformation.     

由Bezier曲线构造的聚合物熔体剪切黏度模型

为提高聚合物熔体剪切黏度模型的描述精度,提出了一个基于二次Bezier 曲线的黏度模型. 模型采用分段函数描述,在对数坐标系中,低剪切速率时的牛顿区和高剪切速率时的幂律区采用线性函数,介于二者之间的过渡区采用二次Bezier 曲线. 通过牛顿区和幂律区的直线延长线构造Bezier 曲线的控制多边形,从而保证三段曲线的光滑过渡. 模型可以明确给出任意温度下低剪切速率时牛顿区的结束点,以及高剪切速率时幂律区的开始点. 拟合算例表明,所提出模型的拟合精度明显高于Cross-Arrhenius 黏度模型.      

老化作用下沥青集料界面力学行为研究

为了探究老化对沥青混合料界面力学性能的影响,将不同温度老化处理的沥青制成的三明治试件进行剪切试验以模拟实际路面发生的破坏情况。通过分析应力-剪切位移的关系与界面破坏情况,研究在不同程度老化作用下沥青集料界面剪切破坏与沥青内部剪切破坏机理,并根据本试验的情况在已有的模型基础上完善出相对符合本试验的理论模型。试验结果表明:每组试验均产生单波峰与双波峰曲线,通过开尔文模型可知,之所以出现双波峰情况,是沥青弹性与粘性的共同作用导致的。随着沥青老化程度的加深,沥青流变性能降低,沥青集料界面的粘附性呈下降趋势,沥青的凝聚性呈先下降后上升趋势,老化加深到一定程度时,双波峰迅速减少,单波峰迅速增加,原因是随着老化程度的加深,沥青凝聚性反向增强,界面破坏只发生在沥青集料界面。     

机械结合面切向接触阻尼计算模型

针对两粗糙表面在法向力和切向力共同作用下相互接触时结合面切向阻尼的问题进行了研究。首先,根据KE模型对单个微凸体在弹性、弹塑性、塑性变形阶段的切向接触行为进行了分析,获得了微凸体在三个变形阶段的黏滑特性;然后,基于GW统计模型建立了一种在微凸体法向弹性、弹塑性和塑性变形机制基础上,考虑微凸体黏滑摩擦行为的机械结合面切向接触阻尼统计模型;最后,分别讨论了机械结合面的法向预载荷、切向激振频率和切向动态位移幅值对机械结合面切向阻尼的影响。研究表明:结合面切向接触阻尼系数随着结合面法向载荷的增大而增大,随着切向激振频率和切向动态位移幅值的增大而减小;在高频率、大幅值下,结合面切向接触阻尼系数几乎与动态位移幅值和激振频率无关。为了验证模型的准确性,构建了动态切向力作用下的结合面切向阻尼试验,其试验结果与理论仿真变化规律与量级基本一致,从而证明了本文所提出的切向阻尼模型的有效性。      

Critical state shear behavior of the soil-structure interface determined by discrete element modeling

The interface between soil and structure can be referred to as a soil-structure system, and its behavior plays an important role in many geotechnical engineering practices. In this study, results are presented from a series of monotonic direct shear tests performed on a sand-structure interface under constant normal stiffness using the discrete element method (DEM). Strain localization and dilatancy behavior of the interface is carefully examined at both macroscopic and microscopic scales. The effects of soil initial relative density and normal stress on the interface shear behavior are also investigated. The results show that a shear band progressively develops along the structural surface as shear displacement increases. At large shear displacement a unique relationship between stress ratio and void ratio is reached in the shear band for a certain normal stress, indicating that a critical state exists in the shear band. It is also found that the thickness and void ratio of the shear band at the critical state decreases with increasing normal stress. Comparison of the DEM simulation results with experimental results provides insight into the shear behavior of a sand-structure interface and offers a means for quantitative modeling of such interfaces based on the critical state soil mechanics.     

Interfacial stress analysis and prediction of debonding for a thin plate bonded to a curved substrate

This paper focuses on the analytical and numerical modeling of the interface between a rigid substrate with simple constant curvature and a thin bonded plate. The interfacial behavior is modeled by independent cohesive laws in the normal and tangential directions, coupled with a mixed-mode fracture criterion. The newly developed analytical model determines the interfacial shear and normal stress distributions as functions of the substrate curvature, during the various behavioral stages of the interface prior to the initiation of debonding. The model is also able to predict the debonding load and the effective bond length. In the numerical model the interface is modeled by zero-thickness node-to-segment contact elements, in which both the geometrical relationships between the nodes of the discretized problem and the interface constitutive laws are suitably defined. Numerical results and comparisons between the predictions of the two models are presented.     

Micromechanics of elasto-plastic materials reinforced with ellipsoidal inclusions

The effective mechanical behavior of an elasto-plastic matrix reinforced with a random and homogeneous distribution of aligned elastic ellipsoids was obtained by the finite element simulation of a representative volume element (RVE) of the microstructure and by homogenization methods. In the latter, the composite behavior was modeled by linearization of the local behavior through the use of the tangent or secant stiffness tensors of the phases. “Quasi-exact” results for the tensile deformation were attained by averaging of the stress-strain curves coming from the numerical simulation of RVEs containing a few dozens of ellipsoids. These results were used as benchmarks to assess the accuracy of the homogenization models. The best approximations to the reference numerical results were provided by the incremental and the second-order secant methods, while the classical or first-order secant approach overestimated the composite flow stress, particularly when the composite was deformed in the longitudinal direction. The discrepancies among the homogenization models and the numerical results were assessed from the analysis of the stress and strain microfields provided by the numerical simulations, which demonstrated the dominant effect of the localization of the plastic strain in the matrix on the accuracy of the homogenization models.     

A Finite Strain Analysis of Cavity Formation at a Rigid Inhomogeneity

The formation of a cavity by inclusion-matrix interfacial separation is examined by analyzing the response of a plane rigid inclusion embedded in an unbounded incompressible matrix subject to remote equibiaxial dead load traction. A vanishingly thin interfacial cohesive zone, characterized by normal and tangential interface force-separation constitutive relations, is assumed to govern separation behavior. Rotationally symmetric cavity shapes (circles) are shown to be solutions of an interfacial integral equation depending on the strain energy density of the matrix, the interface force constitutive relation and the remote loading. Nonsymmetrical cavity formation, under rotationally symmetric conditions of geometry and loading, is treated within the theory of infinitesimal strain superimposed on a given finite strain state. Rotationally symmetric and nonsymmetric bifurcations are analyzed and detailed results, for the Mooney–Rivlin strain energy density and for an exponential interface force-separation law, are presented. For the nonsymmetric rigid body displacement mode, a simple formula for the critical load is presented. The effect on bifurcation behavior of interfacial shear stiffness and other interface parameters is treated as well. In particular we demonstrate that (i) for the smooth interface nonsymmetric bifurcation always precedes rotationally symmetric bifurcation, (ii) unlike rotationally symmetric bifurcation, there is no threshold value of interface parameter for which nonsymmetric bifurcation will not occur and (iii) interfacial shear may significantly delay the onset of nonsymmetric bifurcation. Also discussed is the range of validity of a nonlinear infinitesimal strain theory previously presented by the author (Levy [1]). This revised version was published online in July 2006 with corrections to the Cover Date.     

Continuous analysis of stresses from arbitrary surface loads on a half space

A new form of elemental surface load on a half space is introduced, presuming a quasi-pyramidal variation of load which is doubly linear in each of four rectangular parts of a surface rectangle. Approximations of arbitrary load distributions by sums of such elements are continuous, piecewise linear in two directions and well adaptable. The loads may be normal or tangential. The explicit solutions obtained for all stress and displacement components due to each elemental load involve only elementary functions, are free of the discontinuities which arise with stepwise elements, and are suitable for computing. Some illustrative stress distributions are presented for elemental loads and for multiple pyramidal loads involving both normal and tangential loads. The value of the load continuity in the more complicated analyses of surface cracks is also illustrated.     

Prediction of wheel performance by analysis of normal and tengential stress distributions under the wheel-soil interface

Prediction of wheel performance by analysis of normal stress distribution under the wheel-soil interface was reported by one of our research members. In this study analysis of both normal and tangential stress distributions are included for the prediction of wheel performance. A visco-elastic soil model based on a three-element Maxwell model is used to evaluate normal stress distribution under a wheel running on soft ground. The values of the parameters characterizing the visco-elastic behavior of the soil can be derived from plate penetration tests. A rigid wheel-soil interface model is used to evaluate the tangential stress distribution under the wheel-soil interface. Shear deformation modulus, cohesion and angle of internal shearing resistance of the soil are derived from shear-displacement tests. Test results indicate that both maximum normal and shear stress occur in front of the wheel axle, and the location of peak normal stress shifts backwards towards the wheel axle while that of tangential stress shifts forwards when slippage is increased from a low value. Increasing slippage also causes a decrease in normal stress and an increase in tangential stress. Coefficients of traction and tractive efficiency are low at low slippage, increase with an increase in slippage, and level off at higher slippage.     

Shear band formation analysis in soils by the subloading surface model with tangential stress rate effect

The generalized elastoplastic constitutive equation for soils is proposed based on the subloading surface model extended so as to describe the dependence of both the magnitude and the direction of inelastic stretching on the stress rate tangential to the subloading surface [Int J Plasticity 17 (2001) 117]. It would be applicable to the analysis of deformation of soils in both normal-yield and subyield states for not only lower but also higher stress ratio than that in the critical state. Then, the shear band formation in the rectangular specimen subjected to the biaxial compression under the undrained plane strain condition is analyzed by the generalized equation, and thus the condition for shear band formation and the shear band inclination are discussed in relation to material properties and the state of stress, i.e. the stress-ratio and the normal-yield ratio. These results reveal that the tangential stretching term makes easy to fulfill the necessary condition of shear band formation for not only normal-yield but also subyield states, and further the formation is affected by the material parameter prescribing the approaching degree of the stress to the normal-yield state.     

On attainable lower boundary of the three-dimensional Coulomb-Tresca invariant

An attainable lower boundary of the three-dimensional Coulomb-Tresca invariant is constructed to facilitate the search of plasticity conditions for isotropic bodies which, just as the Tresca-Saint Venant plasticity condition, ensure the hyperbolic analytic type of three-dimensional equations of the mathematical theories of plasticity based on the generalized associate flow law. The construction is performed by using the system of three “two-dimensional” tangential stresses related to the given three-dimensional stress state. It is proved that the Coulomb-Tresca invariant attains its lower bound in any plane strain state where the out-of-plane principal normal stress is intermediate (or median).     

柔性梁斜碰撞问题的非线性动力学建模和实验研究

本文研究柔性梁点面斜碰撞问题。用Hertz接触模型处理法向撞击力,分别用Hertz切向接触模型和Coulomb摩擦力模型处理粘滞状态和滑动状态的摩擦力。从精确的应变与位移的关系出发,用绝对节点坐标法建立了柔性梁的动力学方程。为了准确地处理斜碰撞切向运动的复杂状态,提出滑动-粘滞切换的准则,在此基础上,设计了斜碰撞实验,数值对比了法向撞击力和法向速度的时间历程的仿真计算结果与实验结果,验证了Hertz理论在斜碰撞情况下的正确性。另一方面切向速度的实验与理论的结果对照表明滑动-粘滞切换准则的有效性。