碳纤维-不锈钢层板热载条件下冲击动态响应及层间损伤仿真研究

为了研究碳纤维不锈钢层板的冲击动态响应以及热载荷条件下的冲击性能,采用ABAQUS/Explicit,编写基于复合材料渐进损伤用户子程序VUMAT;引入Johnson-Cook模型,仿真计算了碳纤维增强环氧树脂基复合材料-SS304不锈钢层板热载条件下冲击动态响应过程;分析了其冲击动态响应及渐进损伤,着重讨论了热载荷条件对碳纤维金属层板的冲击能量吸收、接触力等抗冲击性能及失效模式的影响.结果显示,高速冲击载荷作用下,纤维层的脆性断裂、金属层的塑性变形以及纤维层与金属层之间的脱层是碳纤维不锈钢层板的主要失效形式.热载荷的存在直接影响了冲头的接触力,随环境温度升高,接触力总体上降低,子弹的速度衰减越慢,剩余速度增大.结果表明,热载荷降低了纤维金属板的冲击动能吸收特性,弱化了碳纤维金属板的抗冲击性能.无论是纤维金属层板的整体破坏,还是纤维失效、基体失效和脱层失效,热载荷都产生了重要影响.     

芳纶纤维平纹织布在爆炸载荷下的动态响应与失效行为的数值分析

为了研究柔性纤维织布的抗爆性能,通过数值模拟的方法对芳纶纤维织布在爆炸冲击下的响应与失效行为进行了分析。对国产芳纶纤维织布H1000D-AP220进行了力学性能试验,建立了柔性平纹织布的本构模型和爆炸冲击数值分析模型,对不同厚度不同铺层角度织布进行了爆炸冲击数值分析,获得了织布在不同爆炸载荷下的动态响应和失效模式。结果表明,织布在爆炸冲击载荷下主要表现为中心撕裂破孔和简支边界处拉伸撕裂2种典型的失效模式,并伴有明显褶皱,宽度方向出现织布向内侧收缩翻转现象;相比于中间层织布,迎爆面和背爆面吸能较多;在计算分析中改变了织布的层叠角度,获得了更好的抗爆效果。     

爆炸冲击载荷下航空铝合金平板动态响应数值分析方法

为得到适用于爆炸冲击载荷下航空铝合金平板动态响应的数值分析方法,采用LS-DYNA显式动力学分析软件对爆炸冲击载荷下的铝合金平板进行数值仿真计算.主要研究了不同的任意Lagrange-Euler(拉格朗日-欧拉)网格(ALE)输运步算法、流固耦合方式、流固耦合点数量、网格尺寸、有限元单元类型对计算结果的影响.通过计算结果与实验结果的分析对比,表明采用van Leer+HIS输运步算法、罚函数耦合方式、在流体网格与结构网格之间采用3个耦合点、结构网格尺寸与空气域网格尺寸比例设为2∶1、结构单元采用163号壳单元时可以较为准确地计算航空铝合金平板在爆炸冲击载荷下的动态响应,并且能提高计算效率,节约计算时间.     

爆炸冲击载荷下航空铝合金平板动态响应数值分析方法

为得到适用于爆炸冲击载荷下航空铝合金平板动态响应的数值分析方法,采用LS-DYNA显式动力学分析软件对爆炸冲击载荷下的铝合金平板进行数值仿真计算.主要研究了不同的任意Lagrange-Euler(拉格朗日-欧拉)网格（ALE）输运步算法、流固耦合方式、流固耦合点数量、网格尺寸、有限元单元类型对计算结果的影响.通过计算结果与实验结果的分析对比,表明采用van Leer+HIS输运步算法、罚函数耦合方式、在流体网格与结构网格之间采用3个耦合点、结构网格尺寸与空气域网格尺寸比例设为2∶1、结构单元采用163号壳单元时可以较为准确地计算航空铝合金平板在爆炸冲击载荷下的动态响应,并且能提高计算效率,节约计算时间.     

Blast Damage Simulation With the Discontinuous Galerkin Finite Element Method of Bond⁃Based Peridynamics北大核心CSCD

近场动力学是一种积分型非局部的连续介质力学理论,已广泛应用于固体材料和结构的非连续变形与破坏分析中,其数值求解方法主要采用无网格粒子类的显式动力学方法.近年来,弱形式近场动力学方程的非连续Galerkin有限元法得到发展,该方法不仅可以描述考察体的非局部作用效应和非连续变形特性,还可以充分利用有限单元法高效求解的特点,并继承了有限元法能直接施加局部边界条件的优点,可有效避免近场动力学的表面效应问题.该文阐述了键型近场动力学的非连续Galerkin有限元法的基本原理,导出了计算列式,给出了具体算法流程和细节,计算模拟了脆性玻璃板动态开裂分叉问题,并对爆炸冲击荷载作用下混凝土板的毁伤过程进行了计算分析.研究结果表明,该方法能够再现爆炸冲击荷载作用下结构的复杂破裂模式和毁伤破坏过程,且具有较高的计算效率,是模拟结构爆炸冲击毁伤效应的一种有效方法.     

碳纤维/树脂基复合材料曲壁蜂窝夹芯结构的三点弯曲性能北大核心CSCD

为研究碳纤维复合材料(CFRP)曲壁蜂窝结构在三点弯曲载荷作用下的承载特性与失效模式,对不同芯层高度、面板厚度的结构进行了理论预报、数值模拟及试验.首先,根据夹芯结构的主要失效模式,提出了相应的理论预报公式,并绘制了失效机制图;其次,建立了CFRP曲壁蜂窝夹芯结构的有限元仿真模型,对其在三点弯曲载荷作用下的典型失效行为进行模拟;最后,通过模压成型工艺制备了不同尺寸的CFRP曲壁蜂窝夹芯结构,并将试验结果与理论、模拟结果进行比较.结果表明,蜂窝夹芯结构承载能力与芯层高度、面板厚度密切相关,结构芯层及面板刚度随其尺寸的减小而下降,导致结构失效模式由面芯脱黏失效变为面板压溃失效.     

钢管混凝土柱抗爆性能数值模拟与实验验证

利用LS-DYNA非线性有限元程序,基于多物质流固耦合方法,建立了爆炸荷载作用下钢管混凝土柱的动态响应数值模型.对比分析了模拟结果与足尺构件的爆炸破坏实验结果,验证了数值模型和计算方法的有效性,并运用参数化分析方法,研究了截面形式、比例距离、混凝土强度及钢材等级、截面形状特性等关键参数对钢管混凝土柱抗爆性能的影响.研究结果表明:钢管混凝土柱具有优越的抗爆性能,所建立的数值模型能够有效地分析钢管混凝土柱在爆炸荷载作用下的动态影响及破坏形态;圆形截面钢管混凝土柱的抗爆性能优于方形截面;提高材料等级、减小圆形截面钢管混凝土柱的径厚比、增大矩形截面钢管混凝土柱的长宽比,均有利于提升钢管混凝土柱的抗爆特性.     

求解层间界面反平面剪切破坏的剪切梁模型(Ⅰ)——基本特性

在反平面剪切载荷及侧压力共同作用下引起的裂纹及裂纹扩展导致的层间界面失效,是岩土工程层间界面及砌体结构中界面层上典型的失效方式.运用弹性力学和断裂力学的理论原理,提出了能够反映上述层间界面断裂失效问题力学特性的剪切梁模型.文中采用具有应力软化特性的“粘性裂纹”(内聚力裂纹)模型来表述层间裂纹前方损伤过程区的本构行为.对通过粘性层结合在一起的两个弹性板,在反平面剪切载荷及侧压力共同作用下的力学行为作了解析分析计算,研究了层间界面裂纹扩展规律.     

粘弹性运动带动力响应分析

基于Kelvin粘弹性材料本构模型及带运动方程,建立了运动带非线性动力学分析模型.基于该模型和Lie群分析方法推导了匀速运动及简谐运动带线性问题的解析解;基于该非线性模型的数值仿真讨论了运动带材料参数、带稳态运动速度、扰动速度对系统动态响应的影响.结果表明:1)当带匀速运动时,无论系统是线性还是非线性,运动带横向振动"频率"都随着带运动稳态速度增加而减小.2)随着材料粘性增加,系统耗散能力逐渐增强,动态响应逐渐减小.3)当带运动速度简谐波动时,系统动态响应随扰动速度增大而增大.扰动频率对带横向振动影响较大.     

横观各向同性压电材料中正三角形孔口快速传播裂纹的解析解

利用复变函数方法,通过引入合适的数值保角映射研究了横观各向同性压电材料中正三角形孔口快速传播裂纹的反平面剪切问题,并在电非渗透型与电渗透型两种边界条件下,结合柯西积分,导出了力-电耦合作用下以速度v传播时的Ⅲ型裂纹的动态应力强度因子和电位移强度因子的解析解.最后,考虑面内电载荷和面外机械载荷共同作用,分析了三角形孔尺寸、裂纹尺寸、外载变化对裂尖场强度因子的影响.     

Fluid–structure interaction solver for compressible flows with applications to blast loading on thin elastic structures

In this study, we report the development and application of a fluid–structure interaction (FSI) solver for compressible flows with large-scale flow-induced deformation of the structure. The FSI solver utilizes a partitioned approach to strongly couple a sharp interface immersed boundary method-based flow solver with an open-source finite-element structure dynamics solver. The flow solver is based on a higher-order finite-difference method using a Cartesian grid, where it employs the ghost-cell methodology to impose boundary conditions on the immersed boundary. Higher-order accuracy near the immersed boundary is achieved by combining the ghost-cell approach with a weighted least squares error method based on a higher-order approximate polynomial. We present validations for two-dimensional canonical acoustic wave scattering on a rigid cylinder at a low Mach number and for flow past a circular cylinder at a moderate Mach number. The second order spatial accuracy of the flow solver was established in a grid refinement study. The structural solver was validated according to a canonical elastostatics problem. The FSI solver was validated based on comparisons with published measurements and simulations of the large-scale deformation of a thin elastic steel panel subjected to blast loading in a shock tube. The solver correctly predicted the oscillating behavior of the tip of the panel with reasonable fidelity and the computed shock wave propagation was qualitatively consistent with the published results. In order to demonstrate the fidelity of the solver and to investigate the coupled physics of the shock–structure interaction for a thin elastic plate, we employed the solver to simulate a 6.4 kg TNT blast loading on the thin elastic plate. The initial conditions for the blast were taken from previously reported field tests. Using numerical schlieren, the shock front propagation, Mach reflection, and vortex shedding at the tip of the plate were visualized during the impact of the shock wave on the plate. We discuss the coupling between the nonlinear dynamics of the plate and blast loading. The plate oscillates under the influence of blast loading and the restoration of elastic forces. The time-varying displacement of the tip of the plate is the superimposition of two dominant frequencies, which correspond to the first and second modes of the natural frequency of a vibrating plate. The effects of the material properties and length of the plate on the flow-induced deformation are briefly discussed. The proposed FSI solver is a versatile computational tool for simulating the impact of a blast wave on thin elastic structures and the results presented in this study may facilitate the design of thin structures subjected to realistic blast loadings.     

地下结构物在爆炸冲击波作用下的动力分析

提出了一种求解任意形地下结构物在爆炸冲击波作用下的动应力集中问题的半解析方法．爆炸冲击波以平面SH波的形式入射,并用Fourier变换方法将其转换到频域,不同形状地下结构物的导纳函数由复变函数和保角映射的方法求得．利用Fourier逆变换,进一步合成得到时域中的地下结构的动力响应,最后,对正方形、三角形及马蹄形孔洞附近的动应力集中系数作了数值计算,并给出了具体结果．     

Phase Field Modeling of Brittle and Ductile Fracture

The phase field modeling of brittle fracture was a topic of intense research in the last few years and is now well-established. We refer to the work [1-3], where a thermodynamically consistent framework was developed. The main advantage is that the phase-field-type diffusive crack approach is a smooth continuum formulation which avoids the modeling of discontinuities and can be implemented in a straightforward manner by multi-field finite element methods. Therefore complex crack patterns including branching can be resolved easily. In this paper, we extend the recently outlined phase field model of brittle crack propagation [1-3] towards the analysis of ductile fracture in elastic-plastic solids. In particular, we propose a formulation that is able to predict the brittle-to-ductile failure mode transition under dynamic loading that was first observed in experiments by Kalthoff and Winkler [4]. To this end, we outline a new thermodynamically consistent framework for phase field models of crack propagation in ductile elastic-plastic solids under dynamic loading, develop an incremental variational principle and consider its robust numerical implementation by a multi-field finite element method. The performance of the proposed phase field formulation of fracture is demonstrated by means of the numerical simulation of the classical Kalthoff-Winkler experiment that shows the dynamic failure mode transition. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Local and extended models for high-speed,dynamic loading of metals

The purpose of this work is to compare and contrast two different modeling approaches for high-speed, dynamic loading of metals. A so-called “local” model including hardening, strain-rate hardening as well as temperature and ductile damage softening behaviour is extended by gradient plasticity. In particular, both models are compared in the context of their application to the modeling and simulation of dynamic shear-banding in the alloy Inconel 718. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Phase Field Model for Ductile to Brittle Failure Mode Transition

The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in dynamic problems with complex crack topologies. This can be overcome by a diffusive crack modeling based on the introduction of a crack phase field. We outline a conceptual framework for phase field models of crack propagation in brittle elastic and ductile elastic-plastic solids under dynamic loading and investigate the ductile to brittle failure mode transition observed in the experiment performed by Kalthoff and Winkeler [3]. We develop incremental variational principles and consider their numerical implementations by multi-field finite element methods. To this end, we define energy storage and dissipation functions for the plastic flow including the fracture phase field. The introduction of local history fields that drive the evolution of the crack phase field inspires the construction of robust operator split schemes. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)