冲击载荷下炸药装药动态响应的有限元分析及热点形成机理的数值模拟

对稳态温度场中受冲击载荷作用的炸药药柱进行了弹粘塑性分析。在Ｐｅｒｚｙｎａ本 构模型的基础上，作了适当的补充和修正，将流动参数γ、弹性模量Ｅ均视为温度的函数，动态 有限元计算结果表明，计算曲线和实验曲线有很好的近似。为模拟材料中不均匀性的影响，在 药柱中心引入一孔洞，有限元计算结果给出含孔洞药柱的粘塑性动态响应、药柱网格变形图以 及药柱等温线，可以清楚看出在孔洞附近区域有局部高温产生。本文的本构模型和计算方法对 于研究冲击载荷下炸药装药的力学响应以及炸药装药中热点形成机理的数值模拟提供了良好 的基础。     

应变控制的热机械疲劳行为的数值模拟

根据高温合金材料的力学性能，以弹粘塑性本构模型为基础，用数值模拟方法研究材料的热机械疲劳循环特性，模型将应变分为弹性应变、温度应变和粘塑性应变三部分，认为材料在高温循环载荷下呈现明显的弹粘塑性特征，根据虚位移原理建立轴对称体的弹粘塑性计算有限元格式，对于循环机械载荷和循环温度载荷，程序中采用了增量法迭代求解，在非线性项中不仅考虑了机械载荷增量的影响，同时也考虑了温度增量的影响，根据应变控制热机械疲劳的特点，发展了应变增量法的有限元计算方法、通过数值模拟，得到材料在各种循环载荷下的应力—应变响应，数值模拟较好地反映了粘塑性变形过程以及温度变化的效应，所描述的不可逆系统在某一时刻的状态完全由当时的状态参数、内变量、承载时间及塑性应变累积量决定，对带缺口试件的模拟结果显示了程序对复杂轴对称结构进行热机械疲劳计算的有效性。     

混凝土材料冲击特性的研究

基于混凝土材料强冲击加载下的试验研究，提出了两种损伤型动态本构模型： 损伤型黏弹性本构模型和损伤与塑性耦合的本构模型. 通过模型计算结果与冲击试验结果的 比较可发现，随着冲击速度的提高，混凝土材料内部产生了显著的塑性变形，由此 损伤型黏弹性本构模型的应用就存在一些不足. 而损伤与塑性耦合的本构模型由于考虑了裂 纹扩展引起的材料强度和刚度的弱化，以及微空洞缺陷塌陷引起的塑性变形，因而能更好地 用于模拟强冲击载荷作用下混凝土材料的冲击响应特性.     

脆性固体损伤和局部软化的有限元分析

本文研究混凝土、岩石一类工程中常用的应变软化材料的有限元分析方法。在作者以往有关粘塑性损伤本构模型的工作基础上，给出了一组便于有限元计算的本构方程表达式。包括损伤弹性矩阵和局部损伤软化矩阵，分别运用于计算硬化和软化阶段的有限元刚度矩阵；对所提出的本构方程的实验验证计算和一些算例的有限元数值分析，表明文中给出的本构方程是可行的，相应的有限元算法能较好地对损伤固体的局部软化效应进行数值分析，并可成功地追踪应力应变响应的软化曲线     

单向纤维加强复合材料的总体弹粘塑性响应

应用GMC方法计算了单向纤维加强复合材料在承受不同载荷之后的总体松驰响应。Bodner和Partom的弹－粘塑性统一模型用于描述非弹性相的本构关系，该模型不假设存在屈服条件，也就不必指定加载和卸载条件，可以在加载和卸载的任何时刻使用相同的公式，计算表明，与有限元方法相比，利用GMC方法和统一模型计算单向纤维加强复合材料的总体非弹性响应简单快捷。     

含裂纹梁刚塑性动态断裂的非线性线弹簧模型分析

用非线性线弹簧模型分析了带裂纹梁的刚塑性动态断裂问题.在塑性势理论基础上,建立了全塑性状态下的弹簧本构关系,并用此关系导出带裂纹梁刚塑性动态断裂分析的基本方程,计算了在冲击载荷作用下,裂纹梁的动态断裂响应.     

材料的弹粘塑性本构关系和应力波的演化

试验表明，大多数工程材料在冲击载荷作用之下的变形一般都同时包含有可恢复的瞬态性弹性变形和不可恢复的粘滞性塑性变形，即其本构关系可以用弹粘塑性模型来描述。本文从内变量理论出发，探讨了时率相关材料的弹粘塑性本构关系的一般特性，建立了增量型的弹粘塑性本构关系的一般理论框架和普适的表达式，并且对两种最常用的本构模型——Bodner-Partom模型和Johnson-Cook模型给出了在一维应变条件下的具体形式。通过计算和讨论一维应变粘塑性靶板中冲击波的衰减机制和应力波的演化规律，特别是考察各种粘塑性本构模型中的材料参数对冲击波的衰减和应力波的演化的影响，得出了一些可以直接应用或具有一定借鉴价值的结果，为研究应力波的其他衰减机制以及在人防工程中智能防护层设计时新材料的选取奠定了基础。     

炸药冲击响应的二维细观离散元模拟

用离散元法在颗粒尺度模拟了塑料粘接炸药的冲击响应。炸药的细观结构基于Voronoi拼图建成。以含孔洞和不含孔洞的炸药模型为例,说明细观结构对热点的形成和分布有重要影响。在这些算例中,粘塑性变形是温度局域化的主要原因,而摩擦和体积粘性等机制可以忽略。     

基于J-C本构模型的Comp.B炸药落锤冲击数值模拟

为研究Comp.B炸药在惯性冲击下的力学响应特性, 基于Johnson-Cook本构模型, 对 点火前的Comp.B炸药大落锤(400kg)冲击实验进行了数值模拟. 在不同应力加载条件下模 拟计算得到了用状态方程和不用状态方程的$\sigma$-$t$曲线, 对比应力峰值、应 力上升时间等主要力学参数, 结果表明考虑了状态方程的数值模 拟结果与实验吻合最好, 可以很好地模拟炸药惯性冲击下的力学响应.     

构型弹体跌落冲击载荷及结构响应特性

为深入认识跌落冲击条件下构型弹体内部的载荷传递规律及结构响应特征，促进战斗部装药安定性评估和结构设计，结合数值模拟和应力波分析手段，研究了构型弹体在跌落过程中的冲击响应特征，主要关注内部药柱的变形和损伤特性，并讨论跌落姿态、装药构型和跌落高度等因素的影响。结果表明，在跌落冲击条件下，构型弹体装药的变形并非由药柱同壳体的直接撞击作用主控，而主要受到弹体内部应力波传播的影响。装药结构最大变形和损伤区域并不位于药柱外侧同壳体相接触的位置，而位于内部区域。冲击应力波在壳体和药柱之间的透射特征、在壳体和装药内部的反射和叠加特性等决定了药柱的主要变形区域及其变形程度。跌落姿态对药柱的响应特征和变形形貌具有重要影响，导致装药安定性风险从高到低排序的跌落姿态依次为尾部向下垂直跌落、水平跌落、头部向下垂直跌落和倾斜跌落。药柱构型也具有重要作用，其中药柱分段界面容易使得药柱变形程度增大，但对装药过载以及变形分布特征的影响相对较小；隔板结构则容易增大装药过载，同时导致药柱的局域变形位置和变形程度均发生改变。跌落高度对药柱变形区域分布特征的影响较小，对载荷幅值、变形程度和分布范围大小等则具有重要作用，随跌落高...     

TC4动态力学性能研究

材料的本构模型参量是结构冲击动力响应数值模拟的基础。本文运用静态实验机和SHB装置,对TC4在常温~750C°、应变率10-4~103s-1下的力学行为进行了研究,得到了相应的塑性本构模型参量。通过Taylor撞击实验及其数值模拟对本构模型参量进行了验证,表明所得参量可较好描述冲击载荷作用下材料的塑性力学行为。     

强动载荷下金属材料塑性变形本构模型评述

爆炸、高速冲击等强动载荷作用下金属材料塑性流动的本构模型研究是冲击动力学领域备受关注的课题,近几十年来的研究成果极大地推动了冲击波物理、冲击动力学、材料科学等领域的发展. 本文对金属材料塑性变形动力学数值计算中重要的、得到普遍应用的本构模型进行了回顾和评述,进一步明确其适用范围、优点和不足,为冲击载荷下材料响应的数值模拟对本构模型的选取提供参考. 其中最新的、适用范围最广的PTW本构模型应给以特别地关注,其研究方法及应用前景都值得我们借鉴和期待.      

描述非线性动力响应的砂质土本构模型

本文提出描述非线性响应的砂擀土本构模型，由本构模型计算所得到的理论曲线与实验曲线发接近此模型能广泛地用于土工问题的弹一塑性有限元分析中。     

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为。计算采用二维轴对称应变条件下的绝热模型。钨合金的本构方程采用热粘塑性形式的Johnson Cook模型。为了得到不同尺度的变形信息 ,计算中用了两种网格 ;先用粗糙网格分析试件变形局部化的概貌 ;接着 ,用细密网格 (在变形局部化区域 ,网格尺寸达到 10 m)分析绝热剪切带的形成和发展。有限元模拟得到的绝热剪切带位置和方向与实验一致。计算结果表明 ,绝热剪切带的形成和发展与试件的应力状态密切相关。     

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模拟

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为.计算采用二维轴对称应变条件下的绝热模型.钨合金的本构方程采用热粘塑性形式的Johnson-Cook模型.为了得到不同尺度的变形信息,计算中用了两种网格;先用粗糙网格分析试件变形局部化的概貌;接着,用细密网格(在变形局部化区域,网格尺寸达到10μm)分析绝热剪切带的形成和发展.有限元模拟得到的绝热剪切带位置和方向与实验一致.计算结果表明,绝热剪切带的形成和发展与试件的应力状态密切相关.     

聚奥-9C装药的传爆管殉爆

为了研究冲击波作用下引信传爆装置的响应规律,进行了以主发炸药为RDX-8701、被发装置为聚奥-9C（JO-9C）装药的传爆管（含导爆药柱）的殉爆实验。通过观测残留传爆药、壳体和见证块变形,判断传爆管的爆炸程度,分析了殉爆过程中JO-9C爆轰波的成长历程及传播规律。采用AUTODYN软件建立了殉爆实验有限元模型,计算模型中主要考虑了主发炸药产生的爆炸冲击波对传爆管的冲击响应。基于流固耦合方法,通过调整距离模拟计算得到了传爆管的临界殉爆距离和殉爆安全距离。结果表明,传爆管上端的侧角受到爆炸冲击后产生的爆轰波沿斜下方传播,使传爆药柱完全爆轰,并引起导爆药柱发生殉爆;数值模拟结果显示,JO-9C装药的传爆管临界殉爆距离为5.7 mm,殉爆安全距离为8.8 mm。     

考虑运动硬化的DD3镍基合金力学行为研究

详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,引入了运动硬化规律,考虑了温度和应变率对晶体塑性变形的影响,通过针对每个滑移系考虑屈服准则和流动规律建立了晶体塑性模型. 对积分过程进行了推导,通过编写ABAQUS材料用户子程序(UMAT), 实现本构模型的有限元积分算法. 在此基础上模拟了DD3镍基单晶合金在单轴拉伸和循环载荷下的响应,并与实验数据进行了对比. 利用该模型可以很好地模拟镍基单晶所具有的各向异性特性,体现了镍基单晶在循环载荷作用下的拉-压不对称性.      

含损伤材料的热粘塑性本构关系和柱壳破裂研究

以含内变量的本构关系理论为基础 ,结合材料损伤演化方程 ,并考虑了温度和损伤对材料参数的影响 ,得到了增量形式的热粘塑性本构关系的普适显式表达式。然后使用Bodner幂函数型粘塑性模型 ,具体推导了其增量形式的热粘塑性本构方程。接着结合在实践中有重要意义的内部爆炸载荷作用下的柱壳破裂问题 ,建立了含损伤热粘塑性柱壳破裂问题的完备方程组 ,使用有限差分方法 ,完成了对问题的数值模拟 ,并对结果进行了分析。计算结果与实验结果符合良好。     

不同升温速率下复合药柱烤燃实验与数值模拟研究

为了研究不同升温速率下复合炸药JO-9159/JB-9014烤燃实验的热反应规律,建立了复合炸药的烤燃模型,利用有限元程序LS-DYNA3D对不同结构的复合药柱在烤燃过程中的热响应情况进行了数值模拟,并利用实验进行了验证,结果显示模拟结果可信。利用已建立的模型对5 K/h、3 K/min和10 K/min等3种不同升温速率下复合药柱烤燃过程进行了数值模拟,结果表明:升温速率和装药结构的不同对复合药柱的点火时间和位置有较大影响,随着升温速率的增大,点火时间变短,点火位置由药柱的中心处逐渐移至药柱的两端边缘,升温速率较小时,复合药柱的热安定性取决于内部高能炸药的特性,升温速率较大时,复合药柱的热安定性与单一钝感药柱性能近似。因此,只有在较大的升温速率下,钝感炸药内部嵌入高能炸药才能既提高整体药柱的威力,又保证其具有较好的热安定性。     

混凝土WW三参数率相关动态本构模型

根据一致粘塑性模型理论，结合所做的混凝土单轴动力特性试验，对常用的混凝土本构模型Willam—Warnke三参数模型进行改进，引入应变率的影响，推导了Willam—Warnke材料的动态率相关本构模型；并与试验的应力应变曲线进行比较，结果表明本文模型能够较好地反应混凝土的动力特性；最后，通过对一受冲击荷载的简支梁进行计算，并与线弹性及率无关塑性结果进行比较，结果表明考虑了应变率的影响后梁的动力特性发生了较大的变化。