非均质结构PBX炸药的动态压缩过程模拟

采用离散元方法,构建了考虑PBX炸药晶体颗粒与黏结剂非均匀性的计算模型,通过在炸药试样两端采用相向飞片加载、将原SHPB实验的应力边界替代为速度边界的方法,开展了PBX炸药的动态压缩数值模拟研究,再现了考虑PBX炸药非均质特性的SHPB实验过程,获得了不同加-卸载路径下炸药应力应变曲线以及相应的损伤破坏图像。模拟结果表明:在动态压缩过程中,虽然PBX炸药处于整体应力平衡,由于PBX高度非均质性,内部应力分布并不均匀,晶体间应力以应力桥形式传递;应力曲线为试样整体平均,试样内局部所受真实应力可能高于曲线峰值应力;出现卸载回滞曲线的试样发生局部损伤破坏,而完全软化曲线的试样发生整体失稳破坏。     

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

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

ANALYSIS OF MECHANICAL CHARACTERISTICS OF FROZEN CLAY AND MODIFIED S-M CREEP MODEL1)

为合理描述人工冻结法施工的矿井土层力学特性,对某矿井土样的重塑土进行单轴抗压强度实验和单轴压缩蠕变实验,得到重塑冻结黏土在不同温度及载荷加载等级下的应力-应变曲线和蠕变曲线。实验结果表明：冻结温度越低,土样的单轴强度越大;相同冻结温度下,土样的蠕变变形随着应力水平的升高而增大。单轴压缩蠕变的等时应力-应变曲线随时间发展向应变轴靠拢;土样经历初始蠕变和等速蠕变两个阶段,在较高应力水平下有进入加速蠕变的趋势;对S-M模型中各参数的意义进行修正并考虑温度的影响,得到人工冻结黏土改进S-M蠕变显式模型,然后采用粒子群算法对人工冻结黏土S-M蠕变模型参数进行优化。改进S-M蠕变显式模型理论计算值与实验值吻合良好,表明改进S-M蠕变显式模型能较好模拟人工冻结黏土的蠕变特性。改进S-M蠕变显式模型为人工冻土蠕变计算提供一种新方法。     

冻结黏土力学特性分析及改进S-M蠕变模型

为合理描述人工冻结法施工的矿井土层力学特性,对某矿井土样的重塑土进行单轴抗压强度实验和单轴压缩蠕变实验,得到重塑冻结黏土在不同温度及载荷加载等级下的应力-应变曲线和蠕变曲线。实验结果表明:冻结温度越低,土样的单轴强度越大;相同冻结温度下,土样的蠕变变形随着应力水平的升高而增大。单轴压缩蠕变的等时应力-应变曲线随时间发展向应变轴靠拢;土样经历初始蠕变和等速蠕变两个阶段,在较高应力水平下有进入加速蠕变的趋势;对S--M模型中各参数的意义进行修正并考虑温度的影响,得到人工冻结黏土改进S--M蠕变显式模型,然后采用粒子群算法对人工冻结黏土S--M蠕变模型参数进行优化。改进S--M蠕变显式模型理论计算值与实验值吻合良好,表明改进S--M蠕变显式模型能较好模拟人工冻结黏土的蠕变特性。改进S--M蠕变显式模型为人工冻土蠕变计算提供一种新方法。     

混凝土HJC本构模型参数的研究

运用非线性有限元动力分析软件LS-DYNA模拟了混凝土SHPB试验.将试验采集的入射波数据作为压杆的端面载荷,模拟条件尽量符合试验条件,确保模拟精度.将数值模拟与试验结果相结合,研究了混凝土HJC模型参数的确定方法,得到了C60混凝土的相关计算参数,模拟结果得到的应力-应变曲线与不同应变率下混凝土SHPB试验应力-应变曲线吻合较好.分析了HJC模型的本构特性,比较了与金属JC模型的差异与共性,讨论了HJC模型存在的问题.研究结果表明,未考虑混凝土弹性响应的应变率效应,横向效应的影响与采用的失效方式密切相关.     

110甲基乙烯基硅橡胶材料参数的确定

应用INSTRON5544对110甲基乙烯基硅橡胶材料进行了单轴大变形压缩的力学特性测试。本文根据已有的Mooney-Rivilin模型和Yeoh模型的本构关系,采用泰勒级数展开的办法获得了基于小变形的本构关系表达式。利用本文所得本构关系,将实验数据进行了非线性拟合处理,得到了相应的Mooney-Rivilin模型和Yeoh模型的材料参数。拟合所得的名义应力-应变曲线与实测的名义应力-应变曲线吻合较好。利用非线性拟合所得的Mooney-Rivilin模型和Yeoh模型的材料参数,应用ABAQUS软件进行有限元分析计算。通过有限元的分析计算,验证了Mooney-Rivilin模型和Yeoh模型材料参数的正确性与可行性。     

传爆药静态压缩力学性能及损伤特性研究

为研究聚黑-14C（JH-14C）传爆药静态压缩力学性能及损伤特性,开展准静态压缩实验,获得了不同应变率下的应力-应变曲线,建立了描述不同应变率下JH-14C力学行为的非线性本构模型;利用扫描电镜（SEM）对回收试样进行细微形貌观测,获得了准静态压缩JH-14C损伤特性的表征。结果表明:JH-14C压缩强度随应变率的升高而提高;实验与计算结果对照验证了本构模型的有效性;准静态压缩实验中,JH-14C主要损伤模式为脱湿和穿晶断裂。     

基于数字图像相关性方法的脆性材料拉氏反分析实验技术

利用分离式Hopkinson压杆作为加载系统,借助超高速数字图像相关性分析方法,发展了长杆试件拉氏反分析实验技术,并用于研究脆性材料在小变形条件下的动态本构特性。通过超高速相机实时拍摄冲击加载下长杆试件变形的散斑图像,再对散斑图像进行数字图像相关性（digital image correlation,DIC）分析,获得长杆试件表面速度场和应变场。随后,以脆性材料PMMA（polymethyl methacrylate）为例,从DIC分析得到的速度场中提取出不同拉格朗日位置上质点速度时程曲线,构建路径线连接整个速度场,再结合零初始条件,数值求解得到了试件中的应力时程曲线,消去时间参数后,获得了脆性材料PMMA的动态应力-应变曲线,并与Hopkinson压杆实验和准静态压缩实验的结果进行了对比,揭示了PMMA材料在小应变条件下的黏弹性本构响应特征。     

准静态压缩应力-应变曲线测量方法的探索

本文介绍了用应变片直接测量材料的准静态应力-应变曲线的试验研究。在MTS810材料试验机上分别对93W、G50、砂浆等几种材料进行了准静态压缩试验。由于仪器的系统误差不能由MTSCod规准确的得到材料的真实的弹性变形。为此在试件的中部贴应变片得到材料的弹性变形,塑性变形仍旧由MTSCod规记录,从而得到试样的真实的应变,准确获得准静态压缩应力-应变曲线。试验结果表明:在试件中部贴应变片的方法能够准确得到该材料的杨氏模量;在试件两端垫块的刀口上安置Cod规可得到材料的应力应变曲线。两者的合理组合即可得到准确而完整的准静态压缩应力-应变曲线。试验中还发现准静态实验中试件的断面加工不平,偏心压缩等都会影响E的准确测量。     

基于BP神经网络与小冲杆试验确定在役管道钢弹塑性性能方法研究

为了能够在不停输油气工况下获得在役管道材料的弹塑性力学性能,提出了一种人工智能BP (backpropagation)神经网络、小冲杆试验与有限元模拟相结合,通过确定材料真应力-应变曲线从而获得材料弹塑性力学性能的方法.首先,通过系统改变Hollomon公式中的参数K, n值,获得457组具有不同弹塑性力学性能的假想材料本构关系,其次,将得到的本构关系代入经试验验证的含有Gurson-Tvergaard-Needleman(GTN)损伤参数的小冲杆试验二维轴对称有限元模型,通过有限元计算得到了与真应力-应变曲线一一对应的457条不同假想材料的载荷-位移曲线,最终将两组数据作为数据库输入BP神经网络进行训练,建立了同种材料小冲杆试验载荷-位移曲线与真应力-应变曲线之间的关联关系.通过此关联关系,可利用试验得到的小冲杆载荷-位移曲线获取在役管道钢的真应力-应变曲线,从而确定其弹塑性力学性能.通过对比BP神经网络得到的X80管道钢真应力-应变曲线与单轴拉伸试验的结果以及引用现有文献中不同材料的试验数据对此关系进行验证,证明了该方法的准确性与广泛适用性.     

Four-Point Bending Test of Determining Stress-Strain Curves Asymmetric between Tension and Compression

A method of determining both uniaxial tension and compression stress-strain curves from the result of a single four-point bending test was demonstrated. Stress-strain curves of magnesium showing tension-compression asymmetry due to twinning deformation and those of an S45C steel due to the Bauschinger effect were calculated. The Mayville-Finnie equation was modified slightly for this calculation. The calculation is sensitive to small change in the slope of bending curve, revealing an aspect of inverse problem.     

人工冻结深部软岩单轴力学性能试验及蠕变模型1）

为了了解深部软岩在冻结条件下的单轴力学性能,以东北地区的原状泥砂岩为试验对象,利用自行研制的WDT-100型人工冻土试验仪器,对其进行不同温度下的人工冻土单轴抗压强度试验和单轴蠕变试验,得到泥砂岩单轴压缩应力-应变关系曲线,各温度下试样的单轴抗压强度以及蠕变曲线.单轴压缩试验结果表明：试样在给定温度和加载速率条件下,单轴压缩应力-应变关系曲线都有较为明显的屈服点,并且都在屈服点后,强度有所提高,出现硬化现象.单轴蠕变试验结果表明：单轴压缩蠕变曲线有非线性特征,单轴压缩蠕变的等时应力-应变曲线向应变轴靠拢;单轴压缩时蠕变模量随时间的增长而降低.最后采用遗传算法优化模型参数,得出泥砂岩蠕变经验方程.与试验结果对比,发现拟合情况较好.     

动态压缩下Zr基非晶合金失效释能机理

为研究Zr基非晶合金动态压缩条件下的失效释能机理，采用力学试验机、霍普金森杆、高速摄影、差示扫描量热分析（differential scanning calorimetry, DSC）、扫描电镜（scanning electron microscope，SEM）等，得到了材料应力应变曲线、高速摄影图像、失效式样微观形貌及DSC曲线，根据实验数据计算了材料的晶化激活能，并拟合了材料的JH-2（Johnson-Holmquist II）模型，对材料动态失效过程进行有限元数值模拟。实验结果表明，压缩条件下材料为脆性断裂，断口处观察到典型的脉状纹样及液滴状结构，材料失效过程伴随着释能现象；数值模拟结果表明，材料裂纹局部的瞬时内能大于材料晶化激活能。动态压缩下材料的失效释能机理即为材料破碎释放储存的弹性势能，并导致材料局部晶化释能，释能强度与应变率成正相关。     

Single-test determination of the friction coefficient and stress-strain curve for application in metal-forming analyses

The compression test of a ring is customarily used to determine the friction coefficient in metal-forming problems. The aim of this paper is to describe a procedure to obtain, for metals in the large deformation range, the stress-strain curve and friction coefficient through the same test. The friction coefficient is evaluated by comparing curves obtained by a finite element code with the experimental curves. The stress-strain curve is obtained by inverting an analytical solution, calculated by a lower-upper-bound approach. Good results in terms of accuracy are obtained.     

AUTOMATIC IDENTIFICATION OF MESOSCOPIC PARAMETERS OF DISCRETE ELEMENTS OF SAND MATERIAL 1)

本文依照四川大学某理论力学教学班实际教学改革的经验,通过具体举例说明计算机技术在理论力学教学中尚未开发的巨大潜力,减小数学运算使学生重点关注理论力学基本原理,批量求解使学生提前接触实际应用案例,普适性求解使学生深化对处理多变量问题的认识。同时本文还探讨了理论力学教学改革的方向和目的,即把教学重点放在普适的教学方法上,借助计算机求解多刚体、多变量问题,以期有助于理论力学教材更好地适应时代发展。     

泡沫铝材料准静态本构关系的理论和实验研究

应用Chen和Lu提出的适用于可压缩弹塑性固体的唯象本构模型框架，建立了泡沫铝的准静态本构模型，推导了三维等比例加载和环向受约束轴向加载下的宏观应力-应变曲线. 对两种泡沫铝材料(开孔和闭孔)进行了4类准静态试验，即单轴压缩、三维静水压缩、三维等比例压缩和侧向受约束轴向压缩实验. 利用单轴压缩和三维静水压缩的实验结果得到了泡沫铝材料的本构参数曲线，并由此预测它在三维等比例压缩和侧向受约束轴向压缩情况下的响应. 理论预测与相应实验结果相比较，三维等比例压缩的结果比较吻合，但与侧向受约束轴向压缩的结果却相差很大. 分析表明，理论预测与侧向受约束轴向压缩实验结果的偏差是由于泡沫铝试件与约束筒之间的摩擦造成的. 研究结果说明, Chen-Lu模型能够很好地描述泡沫铝材料在压缩占主导的应力状态下的响应.     

A non-linear uniaxial integral constitutive equation incorporating rate effects,creep and relaxation

A previously proposed first order non-linear differential equation for uniaxial viscoplasticity, which is non-linear in stress and strain but linear in stress and strain rates, is transformed into an equivalent integral equation. The proposed equation employs total strain only and is symmetric with respect to the origin and applies for tension and compression. The limiting behavior for large strains and large times for monotonic, creep and relaxation loading is investigated and appropriate limits are obtained. When the equation is specialized to an overstress model it is qualitatively shown to reproduce key features of viscoplastic behavior. These include: initial linear elastic or linear viscoelastic response: immediate elastic slope for a large instantaneous change in strain rate normal strain rate sensitivity and non-linear spacing of the stress-strain curves obtained at various strain rates; and primary and secondary creep and relaxation such that the creep (relaxation) curves do not cross. Isochronous creep curves are also considered. Other specializations yield wavy stress-strain curves and inverse strain rate sensitivity. For cyclic loading the model must be modified to account for history dependence in the sense of plasticity.     

双基推进剂的高应变率力学特性及其含损伤ZWT本构

为了解双基推进剂在冲击载荷下的力学特性及本构行为,利用材料试验机和分离式霍普金森压杆（SHPB）对双基推进剂进行了单轴压缩实验,并对实验数据的有效性进行了检验。用二波法对实验数据进行处理,得到了双基推进剂的应力应变曲线。实验结果表明:双基推进剂具有明显的应变率相关性,动态下屈服应力与静态下相比明显提高,且屈服应力表现为应变率对数的双线性关系;双基推进剂屈服应变表现为延脆转换特性,在低应变率下表现为延展性,高应变率下表现为冲击脆化特性。利用含损伤朱王唐（ZWT）本构模型对实验结果进行了拟合,得到了模型中的本构参数,并对损伤因子项进行了分析。通过模型预测曲线与实验曲线的对比,发现含损伤ZWT本构能较好地描述双基推进剂在0～0.14应变范围内的力学特性。