首届全国爆炸与冲击动力学青年学者学术研讨会

<正>爆炸与冲击动力学是一个交叉性的力学分支学科,主要研究爆炸、冲击和能量突然沉积等强动载荷下介质、材料与结构的力学响应、效应及工程技术应用。当前,爆炸与冲击动力学的发展重点和学科前沿主要有非平衡爆轰与爆轰波结构,复杂介质状态方程的本构理论与材料动态力学,复杂结构高速撞击与侵彻动力学,超高速碰撞新原理和新方法,多场耦合模型建立与多尺度高精度计算方法以及在武器装备、航空航天和民用安全等领域的应用研究等。爆炸与冲击动力     

混凝土爆破损伤的SPH-FEM耦合法数值模拟

为了提高计算效率以及更好展现爆炸荷载下混凝土破坏过程,采用SPH-FEM耦合法对混凝土爆破成坑进行模拟。首先结合前人给出的C30混凝土Holmquist-Johnson-Cook（HJC）部分本构参数,通过理论推导等方法确定出剩余的参数;然后代入模型中计算,将数值解与实测数据进行对比;最后以峰值压力和峰值加速度作为考察对象,对HJC模型中21个参数敏感性进行分析。结果表明:SPH-FEM耦合法能直观地模拟爆炸荷载作用下爆坑的发展全过程,且能够较好地处理SPH边界问题;基于所给出的C30混凝土HJC本构参数,采用SPH-FEM耦合法对混凝土爆破破坏进行模拟,计算结果与实测数据吻合度高,表明HJC本构参数的确定具有合理性。此外,还发现HJC本构参数对爆破问题结果的敏感度各不相同,指出对峰值压力和峰值加速度均有较大影响的参数在确定的时候需引起足够的重视。     

利用虚拟激励法结合CFD技术进行复合材料机翼的阵风响应分析

利用随机振动高效算法-虚拟激励法结合CFD技术,分析复合材料机翼模型的大气紊流响应。计算过程中的虚拟激励法和CFD方法均应用自主产权软件平台DDJ-W。虚拟激励法的使用显著提高了计算效率,同时CFD技术的引入可以准确地确定结构的气动参数。利用Dryden模型得到大气紊流频谱激励下机翼的功率谱响应和方差。计算过程表明,虚拟激励法结合CFD技术可以准确高效地进行复合材料机翼的阵风响应分析。     

不同荷载模式下矿用救生舱受力性能的数值模拟

针对某煤矿井下救生舱,利用ANSYS和LS-DYNA软件建立救生舱、瓦斯/瓦斯煤尘和巷道三维有限元模型。分别进行了真实内瓦斯/瓦斯煤尘爆炸作用下救生舱流固耦合数值分析,等效三角波爆炸冲击荷载作用下救生舱动态模拟和静水压力荷载作用下救生舱的极限承载力和极限变形的数值计算,确定救生舱动/静态承载能力变化规律和变形模式,给出瓦斯/瓦斯煤尘爆炸超压与救生舱承载能力和变形模式的变化规律。研究结果表明,等效三角波爆炸冲击荷载作用下该救生舱整体结构的极限超压明显比流固耦合计算分析结果大,流固耦合计算分析结果更接近于实验结果。该救生舱满足爆炸冲击波荷载下的变形要求时,其在静水压荷载作用下的变形也能满足。     

爆炸与冲击问题的大规模高精度计算

爆炸与冲击问题的数值模拟在国防和民用安全领域具有重要的工程实用价值.由于爆炸与冲击问题是一个多物质在高应变率、高温及高压条件下的强非线性的瞬态动力学问题,给数值模拟带来了很多的困难,为此,针对爆炸与冲击问题数值模拟中的一些关键和难点问题开展了研究.提出了三维非线性双曲守恒系统的伪弧长自适应网格算法,分析了算法的实现过程,数值结果表明该算法有效地提高了冲击波强间断处的分辨率.发展了针对气相爆轰数值模拟的附加龙格-库塔方法,对非线性对流项进行显示计算,化学反应源项进行半隐式计算,有效地解决了源项引起的刚性问题,计算结果表明该算法可以准确地捕捉和描述爆轰波的复杂结构和典型特征.针对三维工程实际物理问题中的大规模计算需求,给出了三维多物质流体动力学欧拉数值方法的并行化方法,开发了三维爆炸与冲击问题并行计算程序,并给出了针对该并行程序的测试方法.上述工作有利地解决了爆炸与冲击问题大规模、高精度计算中的一些难题.最后,开展了大口径聚能射流侵彻混凝土靶问题的数值模拟和实验研究,通过典型爆炸与冲击工程问题的计算验证了所研究数值方法的有效性.      

爆室内强爆炸中的流固耦合计算

为评估强爆炸的毁伤效应,在已有的辐射流体力学及结构动力学响应部分工作的基础上,本文采用流固耦合方法对耦合界面进行了处理,并考虑辐射输运引起的能量沉积过程,使得强爆炸产生的辐射波、冲击波和结构动力学响应两部分可以较好地耦合计算,从而能够较完整地模拟爆室内强爆炸的辐射流体动力学过程.计算结果与定性物理分析相符,证明了算法的...     

代理模型在机翼油箱晃动适航审定中的应用

飞机飞行过程中,机翼油箱易受扰动或激励作用而产生低频晃动,引起机翼油箱壁面变形失效.本文将响应面代理模型引入飞机复合材料燃油箱晃动适航审定工作,以晃动频率、晃动幅值和油箱载液比为自变量,以油箱外壁面最大主应变峰值为响应,建立了油箱上下蒙皮响应面统计模型,给出油箱上下蒙皮最大主应变峰值与激励频率、幅值和油箱载液比等因素的显式响应函数.开展物理全尺寸试验和流固耦合数值模拟验证代理模型的有效性.结果 表明,上下蒙皮最大主应变峰值呈非线性趋势,存在最危险载液比使得主应变峰值最大化,且油箱破坏多发生在下蒙皮.     

三维显式有限元程序及炸药冲击起爆应用

冲击和爆炸数值模拟在军事工程和民用领域都应用广泛。基于显式有限元理论,本文介绍了自主研发的三维冲击和起爆流体动力学程序。该程序的功能包含了材料本构、状态方程、炸药反应速率方程和接触碰撞算法,能够模拟复杂条件下结构冲击和炸药起爆问题。由于采用面向对象的编程方法,程序结构相对简单并且容易修改。首先,进行了Taylor杆碰撞数值模拟,并将计算结果与实验结果以及DYNA2D结果进行了比较;接着,基于Lee-Tarver点火增长模型,模拟了冲击加载下裸炸药和带壳炸药的冲击起爆问题。计算结果表明,三维程序计算结果与实验值和商业软件结果都非常吻合,该程序能够较准确地模拟三维结构的冲击和炸药起爆问题。同时,基于面向对象开发的三维计算程序实现了模块化功能,易于后续开发,为冲击和爆炸计算程序的开发提供了参考。     

第二届全国爆炸与冲击动力学青年学者学术研讨会报告综述

简要介绍了第二届全国爆炸与冲击动力学青年学者学术研讨会的情况,概况总结了与会学者所做的报告内容. 报告包括 3 个特邀报告、14 个专题邀请报告和 24 个主题邀请报告,其中主题邀请报告分为爆轰和爆炸动力学、结构动力学和多尺度高性能计算、材料动力学和实验测试技术、复合材料结构的动态力学行为、轻质结构的吸能特性和优化设计等 5 个主题.      

多物理场耦合作用分析的近场动力学理论与方法

广义来说, 近场动力学(peri-dynamics,PD)是假设每个物质点在承受一定范围内的非接触相互作用下,研究整个物理系统演化过程的理论,为涉及非连续和非局部相互作用的问题提供了一个统一的数学框架,具有广泛的适用性.在简要介绍诸多工程对于多物理场模型和数值计算软件的迫切需求后,针对现有商用软件在处理结构非连续演化问题时遇到的瓶颈,引入近场动力学理论和方法. 概述近场动力学固体力学模型,系统阐述近场动力学扩散模型和近场动力学多物理场耦合建模的研究现状和进展,主要涉及电子元器件、电子封装和岩土工程领域的多物理场耦合建模,包括热--力、湿--热--力、热--氧、热--力--氧、力--电、热--电、力--热--电、多孔介质的水--力流固相互作用等非耦合、半耦合与完全耦合模型,强调发展耦合方程数值解法的重要性.最后对扩散问题和多物理场耦合问题的近场动力学理论模型、数值算法和工程应用做进一步展望.      

复合圆柱壳冲击压缩数值模拟及稳定性研究

针对复合圆柱壳在炸药爆轰作用下的动力学响应及在此过程中伴随的失稳问题,研究了其制造工艺中可能出现的缺陷以及圆柱壳中铜线螺旋角和直径对复合圆柱壳稳定性产生的影响。采用SPH-FEM耦合算法,建立了复合圆柱壳二维细节模型,并提出了一种基于圆柱壳内壁粒子速度历史的失稳判据,计算了在不同参数条件下复合圆柱壳的失稳时间及对应的压缩率,对影响复合圆柱壳稳定性的因素进行了评估。分析结果表明,在复合圆柱壳制备过程中存在的折返层缺陷和铜线直径对复合圆柱壳的稳定性有较大影响,而螺旋角度对其稳定性影响不大。     

Elastic waves from localized sources in composite laminates

This paper is concerned with the analysis of elastic waves generated by localized dynamic sources in structural composites. The source can be external, involving acoustic wave loading as in the so called leaky Lamb wave experiment, and low-velocity foreign object impact on the surface of the structure, or internal, as in sudden crack initiation and its rapid growth from existing internal flaws. All three problems are of critical importance in the safe operation of composite structures, due to their vulnerability to hidden delaminations that can occur in composite materials when they are subjected to this type of loads. It is well known that both the dynamic surface loading associated with impact, and the sudden “opening” of an internal crack associated with the extension of a preexisting flaw act as sources of elastic waves in the material of the structure. The research reported here consists of model-based analysis of the guided waves generated by surface loading and microcrack initiation in graphite epoxy composite laminates commonly used in aircraft and aerospace structures. The objective of this study is to develop a mechanics based understanding of the causal relationship between the properties of the source and the characteristics of the waves generated by its initiation and propagation. The results of this research are expected to be useful in developing effective health monitoring systems for new as well as aging aircraft and aerospace structures.     

应用新型动光弹系统研究某设备受冲击载荷作用

结构受到冲击时,其应力分布与在静荷载作用下完全不同.本丈应用新型动光弹实验系统研究了某设备在冲击载荷作用下其内部动态应力分布情况.该系统由部分传统的动光弹实验装置、高速摄影系统和数字光弹图像处理系统组成,并且采用激光光源、落锤冲击载荷加载装置.加工制作了该设备模型的光弹性模型,应用到冲击载荷作用下的实验研究中,得到了该...     

圆筒内旋转细长管与流体界面耦合的数值方法

圆筒内旋转细长管是石油钻采工程中特有结构,细长管不仅与圆筒发生碰撞接触,还与管内流体和管外环空流体耦合,是一个复杂的非线性流固耦合系统。细长管固体域离散成梁单元,采用非线性碰撞接触动力学方程描述;管内外流体离散成六面体单元,采用计算流体动力学方程描述,在耦合界面处用任意拉格朗日欧拉法动网格来处理运动界面。根据梁单元位移...     

民机复合材料结构冲击后压缩设计值的确定方法

冲击后压缩设计许用值是复合材料飞机结构设计的一个非常重要的参数。本文通过分析常用的目视检测方法和飞机结构冲击能量统计结果,从凹坑深度和冲击能量截止值两方面定义了目视勉强可见损伤;根据民机设计兼顾安全性和经济性的特点,作者提出压缩设计许用值中的目视勉强可见损伤应综合考虑,慎重确定。文章还研究了复合材料层压板的抗冲击性能,从工程应用角度给出了结构冲击后压缩设计许用值的确定要点,即用小试样的冲击后压缩试验结果推导的基准系数和环境系数,对元件乃至典型结构件的冲击后压缩试验结果进行修正。     

离散元法在求解三维冲击动力学问题中的应用

提出了三维连结型离散模型,建立了可实现连结型模型(用于连续介质)-接触型模型(用于非连续介质)转化的三维离散元计算程序,用来模拟连续介质转变为非连续介质的力学过程．利用该计算程序对冲击载荷下混凝土块体内(连续体情况下)的应力波传播过程进行了数值模拟．将计算结果的数值与LS-DYNA程序计算的结果进行比较,验证了该计算程序的计算精度．在此基础上,模拟了混凝土块体的动态破坏(连续介质向非连续介质转化)过程．其计算结果可用动画显示,得到的破坏形式与由实验得到的破坏形式相近．两个算例说明该离散元模型及其计算程序是模拟计算伴随有连续介质向非连续介质转变的动态破坏问题的有力工具．     

冰载荷下不锈钢复合板挠度特性分析

以复合板中面的挠度响应作为不锈钢复合板抗冲击性能的评价指标,基于能量法和经典层合板理论,考虑层间结构参数设计,通过横向载荷下的弯曲平衡微分方程,建立冰载荷下不锈钢复合板挠度响应简化解析模型。该分析模型将整个动态响应分析过程分为冰载荷计算分析和动力学方程求解两个阶段。分析了冰载荷模型的面倾角、冲击速度和碰撞位置对冰载荷的影响,确定极端工况参数,汇总接触面的节点力数据;分析了层厚比对挠度响应的影响规律;基于LS-DYNA有限元仿真以及数值算例分析,对比挠度响应仿真结果和解析计算值,验证了本文简化解析模型的准确性,研究结果对不锈钢复合板抗冲击性能分析和评估具有一定的参考价值。     

饱和地基任意形状衬砌内部均布突加荷载作用下的动力响应

地下衬砌结构经常会受到内部动荷载的作用,内荷载引起的衬砌结构的动应力集中备受关注.论文利用波函数展开法和Laplace变换,推导了饱和土中突加荷载作用下衬砌结构和土体的位移、应力、孔压表达式.应用复变函数和保角变换,将任意形状边界映射为圆形边界,利用饱和土和衬砌结构的连续条件和边界条件,求得了任意形状衬砌结构的动力响应...     

Elastoplastic constitutive modeling under the complex loading driven by GRU and small-amount data

In this paper, a data-driven method to model the three-dimensional engineering structure under the cyclic load with the one-dimensional stress-strain data is proposed. In this method, one-dimensional stress-strain data obtained under uniaxial load and different loading history is learned offline by gate recurrent unit (GRU) network. The learned constitutive model is embedded into the general finite element framework through data expansion from one dimension to three dimensions, which can perform stress updates under the three-dimensional setting. The proposed method is then adopted to drive numerical solutions of boundary value problems for engineering structures. Compared with direct numerical simulations using the J2 plasticity model, the stress-strain response of beam structure with elastoplastic materials under forward loading, reverse loading and cyclic loading were predicted accurately. Loading path dependent response of structure was captured and the effectiveness of the proposed method is verified. The shortcomings of the proposed method are also discussed.