一个基于塑性的损伤本构模型的参数辨识

研究了一个用于混凝土的基于塑性的损伤模型本构参数辨识问题。把从实验获得的应力-应变曲线与数值计算中获得的应力-应变曲线的差别,作为局部水平上最小二乘法的目标函数。为了求解这个反问题,局部水平上求解损伤弹塑性正问题的子程序被嵌入到本文的反问题的迭代格式之中。灵敏度系数矩阵是通过有限差分方法近似计算得到的。给出的数值计算例子计算了单轴压缩试验结果的参数辨识问题。采用反分析得到的模型参数值,对单向拉伸及三种不同侧压作用下的压缩试验进行了数值模拟。数值结果表明:本文采用的应力反分析计算格式稳定,且具有合理的准确性,数值计算得到的应力-应变结果可以较准确地拟合实验曲线。     

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

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

等p路径十砂土本构关系的归一化特性及数值建模方法

用神经网络对土体进行建模能反映应力路径相关性、反映土的剪胀剪缩以及反映体应力、剪应力对体应变、剪应变的交互影响,因而成为一种比较理想的建模方式.能否在样本有限的情况下获得精度比较高的本构模型正是主要的研究目的.通过研究中密砂在等p路径下的三轴试验曲线,发现其应力-应变关系曲线在常规应力范围内具有归一化特性.选择合适的归一化指标对砂土三轴试验数据进行归一化,以归一化的试验数据为训练样本进行神经网络训练,得到了比较理想的砂土的神经网络本构模型.本构模型仿真值与试验值符合较好,表明所给出的建模方法是合理的.提出的建模方法可以在所有试验数据的基础上自动实现概率寻优,能有效降低噪声信号的干扰、减小试验数据的分散造成的影响.     

STUDY ON THE PARAMETER VALUE DOMAIN OF MACRO-MICRO CONSTITUTIVE MODEL

宏微观耦合本构模型的参数识别往往通过反分析方法进行,为了使参数识别结果具有高的置信度,需要确定合适的参数取值范围. 基于动态再结晶过程的微观机理以及相应本构方程的数学特征,提出一个确定参数取值范围的方法. 首先详细给出考虑动态再结晶的黏塑性本构模型,并根据模型构造物理机理,提出通过6 步确定该模型参数取值范围的方法;其次,对300M 低合金钢进行不同温度、应变速率下的热变形试验,测试宏观的流动应力-应变数据及微观的组织数据;然后应用提出的方法,依据试验数据,确定参数取值范围;最后,基于确定参数取值范围中获得的知识,对模型进行局部修改,使模型模拟结果更接近实验结果.     

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蠕变显式模型为人工冻土蠕变计算提供一种新方法。     

Uniaxial stress-strain curves from a bending test

A method which can be used to obtain both tension and compression stress-strain curves from a single bend test was first proposed in 1910. Curiously, this apparently useful technique has been very little application since that time. For this reason, the present investigation was conducted to examine the merits and shortcomings of the test. Three materials which show different stress-strain curves in tension and compression, beryllium, cast iron and prestrained copper, were studied to test the versatility of the method. A special four-point bend fixture, designed to minimize loading errors, was designed for this study and is described. Generally, the stress-strain curves obtained from bending and conventional tests showed good agreement for strains less than about five percent. Finally, some advantages and disadvantages of the bending tests are pointed out.     

基于“AM-GoogLeNet+BP”联合数据驱动的混凝土细观模型压缩应力-应变曲线预测

本文结合GoogLeNet卷积神经网络和BP神经网络分别在图像数据挖掘和数据分析方面的良好性能,采用“AM-GoogLeNet+BP”联合数据驱动方法,对混凝土细观模型(含砂浆、骨料及孔隙)的单轴压缩应力-应变曲线进行了有效预测.通过引入力学参量对图像数据驱动的训练结果进行优化,从而提升了神经网络的物理可解释性.基于Python语言实现混凝土细观模型在Abaqus中的自动建模及细观图像生成过程,并将生成的细观图像数据库与相应的压缩应力-应变曲线作为训练数据集.在GoogLeNet中分别引入SENet, ECANet和CBAM三种代表性注意力机制并对三种注意力机制的性能进行对比和分析,以自适应方式提升神经网络对混凝土各相组分的分析能力,并以此得到混凝土细观模型的初步应力-应变预测曲线;将骨料体积分数、孔隙率及初步峰值应力等物理参量作为输入引入BP神经网络以改善峰值应力的预测精度,并与将物理参量直接引入卷积神经网络输入层的方法进行了对比,最后定量给出了骨料体积分数和孔隙率对峰值应力的影响权重.结果表明,对于不同骨料体积分数及孔隙率的混凝土细观模型,该方法均展现了较高的预测精度.本文采用的“...     

滑轨导向式静/动态双轴拉伸实验技术

基于液压伺服高速加载系统,发展了一种材料双轴拉伸力学性能测试技术。利用锥面接触导向驱动方法,把加载锤竖直方向的驱动力转化为水平方向的双轴驱动力,从而实现对十字形试样平面双轴加载。借助有限元数值模拟手段优化了锥面接触角和十字形试样尺寸。当接触锥角为45°时,既有较好的水平驱动转化效率,同时又保持较小的接触力,确保水平驱动加载各组件在弹性变形范围内,可多次重复使用。确定了加载臂狭缝个数、狭缝与减薄区边缘长度和标距段厚度等试样设计关键参数,在十字形试样测试标距段内实现了均匀变形。设计了测力夹持一体化导杆和非接触光学全场应变测试系统,准确获得了试样的应力和应变。利用此平面双轴拉伸加载装置,开展2024-T351铝合金板单轴拉伸实验和激光探测同步性验证实验,验证装置设计的可行性;开展铝合金板材在不同加载速率下的双轴拉伸实验,得到在双轴加载下铝合金板材应力应变曲线,并与单轴加载下实验结果进行了对比分析。     

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

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

Experimental study and numerical modeling of brittle fracture of carbonate rock under uniaxial compression

The brittle carbonate rock taken from the Tarim Oilfield is tested in laboratory under uniaxial compression. The acoustic emission (AE) is used to monitor the microcracking activity in rock during the experiment. Moreover, the 3D tomograms of carbonate rock after uniaxial compression are obtained by using CT imaging technology, which indicates that microcracks mutually interconnect and eventually form macroscopic fractures after failure. The PFC2D is used to model the behavior of brittle rock including microcracks propagation. The stress–strain curve and cracks distribution in rock model are obtained from the PFC simulation. The numerical results agree with the experimental test well.     

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

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

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.     

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

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

MACRO AND MICROSCOPIC PARAMETER MAPPING RELATIONSHIP FOR RED SANDSTONE 1)

颗粒流离散元软件PFC2D在岩土工程中的应用十分广泛,存在的主要问题是如何标定其细观参数,目前大多使用"试凑法",此方法的缺点在于工作量大、效率低等,本文提出以反向传播算法(back propagation,BP)神经网络的方式代替此方法,利用PFC2D 内置的FISH 以及Python 语言对其进行二次开发,使之自动运行并获取40组宏细观参数样本。结果表明：BP神经网络可以快速准确地建立宏细观参数映射关系,误差均保持在0.01 以内,模拟得到的应力--应变曲线与室内试验曲线高度吻合,且无需大量的数据样本便可创建网络,效率较高;另外,经验证本文选用的平直节理模型,可以有效地解决平行粘结模型UCS/TS值偏小问题,确定平直节理模型可以更好地模拟岩石。     

应变率历史对应力应变曲线的影响

利用普通SHPB实验系统、双试件SHPB实验系统，对一特种钢材进行了不同应变率历史的动态压缩实验，获得了不同应变率历史所对应的应力应变曲线。通过量化平均应变率相同的情况下不同应变率历史所对应的应力应变曲线的差别，以及量化应变率历史的恒定程度，初步分析了应变率历史对应力应变曲线的影响。研究结果表明:特别是在较高平均应变率下，应变率历史对试件材料的应力应变曲线有明显的影响，在材料动态本构关系研究中应当考虑应变率历史的影响。     

黄土固结蠕变特性试验研究

基于杨凌地区黄土的单向固结试验资料,分别研究了加载条件下原状黄土和扰动黄土的固结蠕变特性,由试验数据整理分析出不同含水率、不同固结压力下的固结蠕变曲线簇,根据蠕变曲线形态特征,推导出单向固结条件下黄土的应变-时间关系并拟合得到该关系中的参数,在分析参数与固结压力之间关系的基础上推导出黄土的一维固结蠕变本构关系,经试验证明,该关系能较好的描述杨凌地区黄土的固结蠕变特性。另外,基于试验资料以及Buisman次固结系数的概念,分析出次固结系数与固结压力之间的定量关系,为确定固结压力对黄土次固结特性的影响提供了一种新途径。     

FAILURE MODE AND CONSTITUTIVE MODEL OF PLAIN HIGH-STRENGTH HIGH-PERFORMANCE CONCRETE UNDER BIAXIAL COMPRESSION AFTER EXPOSURE TO HIGH TEMPERATURES

An orthotropic constitutive relationship with temperature parameters for plain highstrength high-performance concrete （HSHPC） under biaxial compression is developed. It is based on the experiments performed for characterizing the strength and deformation behavior at two strength levels of HSHPC at 7 different stress ratios including a=σs ： σ3=0.00：-1,-0.20：-1,-0.30 ： -1,-0.40：-1,-0.50：-1,-0.75：-1,-1.00：-1, after the exposure to normal and high temperatures of 20, 200, 300, 400, 500 and 600℃, and using a large static-dynamic true triaxial machine. The biaxial tests were performed on 100 mm×100 mm×100 mm cubic specimens, and friction-reducing pads were used consisting of three layers of plastic membrane with glycerine in-between for the compressive loading plane. Based on the experimental results, failure modes of HSHPC specimens were described. The principal static compressive strengths, strains at the peak stress and stress-strain curves were measured; and the influence of the temperature and stress ratios on them was also analyzed. The experimental results showed that the uniaxial compressive strength of plain HSHPC after exposure to high temperatures does not decrease dramatically with the increase of temperature. The ratio of the biaxial to its uniaxial compressive strength depends on the stress ratios and brittleness-stiffness of HSHPC after exposure to different temperature levels. Comparison of the stress-strain results obtained from the theoretical model and the experimental data indicates good agreement.     

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

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

Measurement of the shapes of foil bulge-test samples

Accurate measurements of important tensile properties of thin metal foils are often quite difficult to achieve in uniaxial tests because of sample-preparation difficulties and the tensile instability called necking. Consequently, hydraulic bulge tests have been introduced as a successful means of suppressing these problems through the use of a simplified specimen geometry and biaxial rather than uniaxial tensile-stress states. Considerable effort has been made by various investigators to relate such biaxial stress-strain and ductility data to uniaxial data, generally following the assumption that the bulge is shaped like a spherical cap. The present study evaluates this assumption for foils by measuring actual shapes with unprecedented accuracy and detail using the two-source holographic technique and a polynomial-spline computer analysis of the resulting interferograms. These measurements were made on nine specimens of 0.127-mm-thick annealed rolled copper foil which had been deformed into bulges of varying heights up to rupture. A comparison is made between the measured shapes and the spherical-cap shape generally assumed in the interpretation of bulge-test data. The spherical assumption gives results which are reasonably valid for the later stages of deformation. Indeed, the stress-strain curve obtained from bulge testing corresponds closely with the uniaxial tensile curves for this material. The strain at failure (i.e., elongation) was greater in the biaxial bulge test than in the uniaxial test but not nearly as great as the strain expected from a theoretical model proposed by Hill. However, all the specimens measured exhibited localized areas with larger radii of curvature. The presence of these “flats” may be associated with a mode of failure in the bulge test which corresponds to necking instability in the uniaxial test, and thereby account for the limited strain to failure.     

滞后细观模型在岩石力学中的应用

对以砂岩为代表的所谓``NME材料'的力学行为研究方面的一些新的概念和模型进行了评介.首先介绍了一种基于所谓``滞后单元'的描述滞后现象的物理模型------Preisach-Mayergoyz(P-M)模型,然后详细阐述了P-M模型应用于模拟岩石的非线性滞后应力应变关系的过程和结果.这种唯象模型很好地描述了宏观上的滞后表现和``离散记忆'效应.接着本文对应变能耗散的力学机制进行了简单分析. 最后,介绍了一种描述弹性波在``NME材料'中传播规律的数学方法, 该方法从一般的弹性波传播规律出发,分析了``NME材料'特殊的力学性质给弹性波传播带来的影响,揭示了产生特殊的弹性波传播规律的原因.