等径通道挤压过程三维有限元模拟

利用三维有限元模型对等径通道挤压过程中变形分布的不均匀性进行了模拟,通过对不同摩擦系数下截面等效塑性应变分布比较发现:沿三个方向截面上的变形分布都不是均匀的,这说明二维有限元模型不能真实模拟等径通道挤压过程中试样中的变形分布。此外,摩擦对等效塑性应变分布及挤压力的影响较大,截面变形不均匀参数和最大等效塑性应变出现的位置随摩擦系数的增大而变化。压力-位移曲线的变化可以结合试样的变形过程来解释。     

花岗岩单轴压缩侧向变形及脆性破坏机制实验研究

为研究花岗岩侧向变形及脆性破坏机制,对花岗岩试件进行单轴压缩实验。利用动态应变采集系统、数字散斑相关方法(DSCM)和显微观测手段,记录并分析花岗岩试件在单轴压缩过程中的宏观侧向应变、局部侧向应变以及破裂面形貌,并与水泥砂浆试件的破坏过程对比,讨论了花岗岩脆性破坏机制。实验与分析结果表明:(1)花岗岩试件在加载初期发生侧向收缩变形,产生并发展于压密阶段,消失于线弹性阶段初期,这主要由于试件内部裂纹闭合造成的;此后,宏观侧向应变持续增长,当侧向应变与轴向应变之比接近0.5时试件破坏;(2)在峰值载荷前很长一段时间内,局部侧向应变在一定范围内波动,临近试件破坏时局部侧向应变最大值和最小值均出现较大幅度的波动,二者差值迅速增大,试件不均匀程度增大,最终导致试件破坏;(3)在峰值载荷前有无塑性屈服阶段是峰值载荷后脆性破坏程度的重要影响因素,而宏观裂纹的贯通程度是峰值载荷后应力降大小的决定因素。     

SHPB试验中岩石试件的端面不平行修正

为研究短圆柱体岩石试件端面不平行对岩石动力学特性测试结果的影响，采用有限元分析软件LS-DYNA对9种端面不平行度和5种杨氏模量的岩石试件开展SHPB(split Hopkinson pressure bar)试验数值模拟，对岩石选用HJC(Holmquist-Johnson-Cook)本构模型。数值模拟结果表明，当端面不平行度在0.40%以内时，端面不平行对动态应力测试结果的影响可忽略不计；但对动态应变测试结果的影响较大。当杨氏模量一定时，平均应变率测试误差和峰值应变测试误差随端面不平行度增大呈线性增大；当端面不平行度一定时，平均应变率测试误差和峰值应变测试误差随杨氏模量增大也呈线性增大。对数值模拟得到的平均应变率测试误差和峰值应变测试误差实施二元线性回归分析，提出了SHPB试验中端面不平行岩石试件平均应变率和峰值应变的修正公式。     

单轴载荷下X80钢的包申格效应研究

本文通过单轴拉伸和压缩试验研究了X80管线钢的包申格效应(BE)。采用正向与反向加载方法研究材料变形历史特性。测定了X80钢的简单拉伸试验曲线,其应力-应变关系表明,该材料具有理想弹塑性特点。为了得到X80钢的BE,在不同预变形下对几个试件分别进行加载,并当给定的预应变值分别达到0.63%,0.67%,0.95%,1.27%和1.55%时就卸载。随后再进行反向加载实验,并记录应力应变曲线。该钢材反向加载时出现加工硬化,且屈服强度比正向加载时要低。正反向加载之间的屈服强度差值随着预应变增加而增大;当预应变超过0.95%时,反向屈服强度达到恒量。实验表明,X80钢的反向加载特性可用Remberg-Osgood关系拟合。最后给出了屈服强度降和预塑性应变之间的经验公式。     

岩石试件端面效应的变形局部化数值模拟研究

试件端面效应直接影响试件内的应力分布、试件破坏形式和剪切带图案。采用FLAC - 3D模拟了端面效应对岩石试件剪切带图案、速度场和位移场的影响。试样端部的缚束作用较强时 ,仅在试样中部出现X形的剪切带 ;反之 ,在试样端部出现了倾斜的剪切带。压头刚性越大 ,剪切带倾角越大。块体具有整体平移的特性 ,剪切带边缘位移梯度却较大 ,反映了块体沿剪切带发生滑动的剪切破坏本质。数值模拟结果可在众多的试验结果中找到佐证。     

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

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

冲击载荷下圆环纵向压缩力学行为研究

选用外径、内径和高度三者比值为6∶3∶2的铝合金圆环试件作为研究对象,在粗磨、细磨和润滑三种端面粗糙条件下,利用SHPB实验方法研究了其纵向动态压缩行为。实验结果表明:随着端面摩擦系数的增大,对应的屈服强度和硬化强度增大明显。同时,圆环纵向压缩应力应变曲线有一个明显的应力下降段。其原因是圆环在塑性压缩过程从内而外发生了动态塑性失稳。由于端面应力不均匀和压缩过程端面摩擦系数不均匀变化,试件前后端面的压缩变形明显不一致。数值模拟分析了不同表面粗糙情况下的圆环纵向压缩过程,得到了和实验相类似的现象。     

钛镍形状记忆合金冲击变形后形状记忆效应的研究

采用SHPB技术和可控速率循环加温条件下变形恢复量测定装置研究了冲击及静载变形后的TiNi形状记忆合金的单程及双程形状记忆特性。发现马氏体状态下的TiNi合金的力学特性显示出明显的应变率强化效应 ,并且高应变率压缩应力应变曲线呈现流动平台。应变率对形状记忆效应的影响具有双重性 ,当外加应力或残余应变较小 ,可逆非弹性变形机制起主导作用时 ,提高应变率可以增加其单程形状记忆效应 ;而随外加应力或残余应变增大 ,当基于位错机制的不可逆非弹性变形机制起主导作用时 ,应变率提高却抑制了其单程形状记忆效应。应变率对TiNi合金双程形状记忆效应的影响视塑性变形的大小而异 ,高应变率动载后的双程形状记忆效应在较小塑性应变时 ,比静载后的要强 ;但在较大塑性应变时两者差别不大。     

基于能量等效原理的应变局部化分析:Ⅱ.有限元解法

以非局部塑性理论为基础,应用状态空间理论,通过局部和非局部两个状态空间的塑性能量耗散率等效原理,提出了一种求解应变局部化问题的新方法,以得到与网格无关的数值解.针对二维问题的屈服函数和流动法则导出了求解非局部内变量的一般方程,并提出了在有限元环境中求解应变局部化问题的应力更新算法.为了验证所提出的方法,对1个一维拉杆和3个二维平面应变加载试件进行了有限元分析.数值结果表明,塑性应变的分布和载荷-位移曲线都随着网格的变小而稳定地收敛,应变局部化区域的尺寸只与材料内尺度有关,而对有限元网格的大小不敏感.对于一维问题,当有限元网格尺寸减小时,数值解收敛于解析解.对于二维剪切带局部化问题,数值解随着网格尺寸的减小而稳定地向唯一解收敛.当网格尺寸减小时,剪切带的宽度和方向基本上没有变化.而且得到的塑性应变分布和网格变形是平滑的.这说明,所提方法可以克服经典连续介质力学模型导致的网格相关性问题,从而获得具有物理意义的客观解.此模型只需要单元之间的位移插值函数具有C~0连续性,因而容易在现有的有限元程序中实现而无需对程序作大的修改.     

轴向尺寸对泡沫铝动静态力学性能的影响

通过材料密度均匀性分析得知,泡沫材料密度分布的均匀性越差,尺寸的影响越明显,对于密度分布区间小、并符合正态分布的试件,尺寸的影响较小.根据Hopkinson杆动态测试的应力均匀性假定,利用石英片监测试件左右端面达到应力均匀所需的时间,试件轴向尺寸越大,惯性效应的影响也越大.选择合适的尺寸,将试件的惯性效应(波动效应)与...     

The frictional sliding response of elasto-plastic materials in contact with a conical indenter

Over the past decade, many computational studies have explored the mechanics of normal indentation. Quantitative relationships have been well established between the load–displacement hysteresis response and material properties. By contrast, very few studies have investigated broad quantitative aspects of the effects of material properties, especially plastic deformation characteristics, on the frictional sliding response of metals and alloys. The response to instrumented, depth-sensing frictional sliding, hereafter referred to as a scratch test, could potentially be used for material characterization. In addition, it could reproduce a basic tribological event, such as asperity contact and deformation, at different length scales for the multi-scale modeling of wear processes. For these reasons, a comprehensive study was undertaken to investigate the effect of elasto-plastic properties, such as flow strength and strain hardening, on the response to steady-state frictional sliding. Dimensional analysis was used to define scaling variables and universal functions. The dependence of these functions on material properties was assessed through a detailed parametric study using the finite element method. The strain hardening exponent was found to have a greater influence on the scratch hardness and the pile-up height during frictional sliding than observed in frictionless normal indentation. When normalized by the penetration depth, the pile-up height can be up to three times larger in frictional sliding than in normal indentation. Furthermore, in contrast to normal indentation, sink-in is not observed during frictional sliding over the wide range of material properties examined. Finally, friction between indenter and indented material was introduced in the finite element model, and quantitative relationships were also established for the limited effects of plastic strain hardening and yield strength on the overall friction coefficient. Aspects of the predictions of computational simulations were compared with experiments on carefully selected metallic systems in which the plastic properties were systematically controlled. The level of accuracy of the predicted frictional response is also assessed by recourse to the finite element method and by comparison with experiment.     

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

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

板材冲压成形变形局部化与动边界摩擦约束

提出适于具有强烈弯曲效应的大应变，大位移与界面摩擦三重非线性空间曲壳有限元列式。给出一种动边界摩擦约束的元级代数消元方法，就板材在方盒型冲头作用中压胀形过程，引入几种现代非经典塑性本构理论以及空单元技术，研究了成形过程中的变形局部化与断裂问题。     

平面幂律塑性Ⅰ型裂纹尖端应力场全解

在航空航天、船舶、石油管道和核电等领域,服役结构或部件在长期极端条件下运行,不可避免地会产生裂纹,因此,为研究含裂纹结构的准静态断裂行为,必须了解裂纹尖端附近区域的应力应变场特点.对于幂律材料裂纹构元,研究平面应变和平面应力条件下Ⅰ型裂纹尖端应力场的解析分布.基于能量密度等效和量纲分析,推导了能量密度中值点代表性体积单元(representative volume element, RVE)的等效应力解析方程,并定义其为应力因子,进而针对有限平面应变和平面应力紧凑拉伸(compact tension, CT)试样和单边裂纹弯曲(single edge bend, SEB)试样,以应力因子作为应力特征量,并构造用于表征裂尖等效应力等值线的蝶翅轮廓式和扇贝轮廓式三角特殊函数,提出描述幂律塑性条件下平面I型裂纹尖端应力场的半解析模型.该半解析模型形式简单,对CT和SEB试样的裂尖应力场的预测结果与有限元分析的结果比较表明,两者之间均密切吻合,模型公式可直接用于预测Ⅰ型裂纹尖端应力分布,方便于断裂安全评价和理论发展.     

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

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

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

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

Development of a solution method for frictional contact problems with complex loading

In this work, solution methods for frictional contact problems are extended to the case of moving punches and to the external loading history-dependent system states. To solve the frictional contact problems in the contact area, an iterative method is developed and implemented. Solutions of two-dimensional problems are constructed using the boundary element method. Numerical analysis is aimed at the quantitative study of effects such as the interaction of contact pressure and friction forces, estimates of the friction force differences due to the differences in the choice of local basis for the calculation of normal pressure and friction forces, and evaluation of the effects of complex loading (rotation of the rigid punch after its preliminary penetration into the solid). We find that, for the same definition of the friction force, different initial approximations lead to the same solution. At the same time, the friction forces defined either as projections onto the common tangent plane or as projections onto the plane tangent to the punch can differ quite substantially. Similar conclusions are derived for the solutions corresponding to single or multiple loading steps. The work relies on the variational principle for the solution of contact problems and numerical algorithms developed for the problems with one-sided constraints. The variational principle was first applied by Signorini [1] to the determination of the stress-strain state in a linearly deformed body in a rigid smooth shell. The modern view of the problem and its generalizations to the frictional problems and some other problems involving unilateral constraints in given in the monograph [2]. Finite difference and finite element methods in application to the problems with unilateral constraints are described in [3]. Analytical solution methods are developed in the monographs [4–6].     

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.     

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

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