模拟结构倒塌破坏的有限质点法

有限质点法是以向量式力学为基础的新兴结构分析方法,本文将其应用于冲击荷载作用下的网壳结构的倒塌破坏模拟中。以空间杆单元为例建立了有限质点法的基本方程,推导了求解几何和材料非线性问题的基本公式。为计算断裂问题,建立了空间杆单元的断裂准则和断裂模型,发展了有限质点法进行断裂分析的基本算法。通过对某双层网壳冲击荷载下破坏过程的模拟和分析,验证了该方法在结构倒塌破坏过程模拟中的有效性和适用性。     

冲击荷载下单层球面网壳的失效机理

为研究单层球面网壳在冲击荷载下的失效机理,在ANSYS/LS-DYNA中建立60m跨度K8型单层球面网壳与圆柱形冲击物的数值模型并进行数值分析,总结了网壳结构的4种失效模式。通过对失效全过程的分析,从能量的角度将失效过程分为能量施加、能量传递与损失、能量消耗3个阶段。之后分别从能量传递与杆件破坏形式2方面揭示了网壳的失效机理。能量分析表明:剩余能量(Elf)对结构最终动力响应及失效模式起决定作用,而Elf只是初始冲击能量中除去冲击物穿透损失与网壳局部破坏损失后的剩余部分。通过对杆件破坏形式的分析发现:杆件的破坏可能滞后于冲击荷载的作用,且杆件的破坏形式决定其传递能量的能力,当杆件发生拉伸破坏时,其强度被充分利用,传递的能量最多,Elf值较大,网壳整体破坏严重。杆件的破坏形式与Elf及网壳整体的失效模式间有很好的对应关系。     

火灾下混凝土结构破坏模拟的纤维梁单元模型

为分析和模拟结构构件在火灾下的失效破坏过程,本文基于建筑结构分析中常用的纤维梁单元模型,建立了钢筋混凝土梁、柱构件的火灾破坏数值模型.此模型将构件截面划分成多个纤维,从而可以模拟构件截面的不均匀温度场分布以及高温下混凝土材料的开裂、压碎和钢筋屈服等行为.并根据全拉格朗日描述方法,推导了纤维梁单元的切线刚度矩阵,建立了纤维梁单元的增量求解方程.最后,将本文模型的模拟结果和多个具体试验结果进行了分析比较,进一步验证了模型的可靠性.     

钢筋混凝土柱破坏过程扩展有限元数值模拟

为描述钢筋混凝土柱的破坏过程,分别建立了由三维实体单元和三维杆单元构成的精细化模型以及结合纤维梁单元的纤维梁与实体单元耦合模型。耦合模型中,钢筋混凝土柱下段混凝土和钢筋分别用三维实体单元与三维杆单元建模,通过耦合连接实现精细单元与粗糙单元之间的变形协调。基于混凝土弹塑性本构关系,运用扩展有限元方法模拟钢筋混凝土柱的荷载-位移曲线以及开裂过程,并进行比较。数值模拟与试验结果的比较表明,模拟得到的裂缝位置以及荷载-位移曲线与试验结果吻合较好。扩展有限元具有无需重划网格、无需预设裂纹的优点,能有效地模拟筋混凝土柱的开裂过程。耦合模型具有计算效率高的优势,且能较好地模拟构件的破坏模式。     

基于能量平衡原理的有限单元耦合方法

结构非线性数值计算分析应真实反映局部损伤破坏细节,以作为损伤演化全过程分析的依据。对同类构件,有限单元耦合方法可以解决破坏细节与整体模拟的空间尺度差异问题。基于能量平衡原理,建立了梁与实体单元、梁与壳单元以及壳与实体单元的耦合方程,适用于结构的损伤数值计算。对某RC框架结构原位推覆试验的损伤数值分析表明,有限单元耦合模型能正确反映整体结构的承载力和变形性能,并且能准确反映局部损伤破坏细节。     

用于模拟颗粒增强复合材料破坏过程的梁—颗粒细观模型的实验验证

本文的主要目的是验证梁－颗粒细观模型在模拟混凝土和砂岩类颗粒增强复合材料连续破坏过程的有效性，文中首先介绍了梁－颗粒细数值模型的基本原理，然后给出了由细钢丝粘结成的正方体试验的单轴抗压实验结果，最后用梁－颗粒细观模型对物理实验进行了数值模拟，研究结果表明，物理实验和数值模拟所得到的试件破坏模式和荷载－位移曲线，两者基本一致，从而初步证明了梁－颗粒细观数值模型及模拟颗粒增强复合材料破坏过程，以及解释     

泡沫混凝土弹塑性损伤模型在组合式防护结构中的应用

为了将新型泡沫混凝土动态弹塑性损伤模型应用到防护结构中,首先开展组合式防护结构预制孔装药爆炸试验;随后利用新泡沫混凝土材料模型对试验进行数值模拟验证,并将新模型的模拟结果与LS-DYNA中Soil and Foam模型的模拟结果进行对比;最后,基于验证的数值模型,开展以梯度泡沫混凝土作为分配层的组合式防护结构预制孔装药爆炸的数值模拟,探讨梯度泡沫混凝土层界面层数和排列方式对组合式防护结构抗爆性能的影响。结果表明,新泡沫混凝土材料模型的模拟结果与试验结果吻合良好,与Soil and Foam模型相比,新模型在应力波传播和损伤破坏方面预测更好,泡沫混凝土层界面层数和排列方式对作用在主体结构上的应力以及分配层的损伤破坏情况有一定的影响。     

砌体结构破坏倒塌过程的离散单元模拟

基于传统颗粒离散单元模型（DEM）,开发了一种绑定式离散单元模型（BTDEM）用于砌体结构崩塌机制的研究,并对一个典型砌体结构的破坏倒塌过程进行数值模拟。结果证实BTDEM方法能够详细地模拟出砌体建筑结构的崩塌过程,同时还可以模拟出砌块的破碎,是一种普适性较强的模拟砌体结构建筑的高效数值方法。BTDEM方法具有计算速度...     

梁-颗粒模型在混凝土侵彻问题中的应用

用二维梁-颗粒模型BPM2D（beam-particle model in two dimensions）模拟了混凝土侵彻问题。在模型中用3种类型梁单元形成混凝土数值试样。每种类型梁单元的力学性质均按韦伯分布随机赋值来模拟混凝土细观结构的非均匀性。利用此模型对在钢弹侵彻下混凝土圆板的破坏过程进行了数值模拟。通过将计算结果与实验数据及LS-DYNA程序模拟结果比较,表明梁-颗粒模型可有效地应用于计算和模拟材料从连续介质向非连续介质转变的动态破坏问题。     

爆炸荷载作用下钢筋混凝土梁的动力响应及破坏形态分析

对基于Timoshenko梁理论建立的非线性动力有限元法作了改进。根据压区理论得到混凝土的平均剪应力和平均剪应变关系,建立了能反映箍筋的抗剪作用的材料模型;此外,对结构在爆炸荷载作用下可能出现的各种响应现象进行了描述,以准确地预测梁破坏时不同位置截面上钢筋和混凝土的受力、变形及破坏情况。应用改进的材料模型,对爆炸荷载作用下的五个钢筋混凝土试验梁的动力响应和破坏形态进行了数值模拟,结果表明,该数值方法能较好地模拟钢筋混凝土梁的弯曲、弯剪和剪切等破坏形态。     

Inverse Estimation of Cohesive Fracture Properties of Asphalt Mixtures Using an Optimization Approach

Tensile cracking in asphalt pavements due to vehicular and thermal loads has become an experimental and numerical research focus in the asphalt materials community. Previous studies have used the discrete element method (DEM) to study asphalt concrete fracture. These studies used trial-and-error to obtain local fracture properties such that the DEM models approximate the experimental load-crack mouth opening displacement response. In the current study, we identify the cohesive fracture properties of asphalt mixtures via a nonlinear optimization method. The method encompasses a comparative investigation of displacement fields obtained using both digital image correlation (DIC) and heterogeneous DEM fracture simulations. The proposed method is applied to two standard fracture test geometries: the single-edge notched beam test, SE(B), under three-point bending, and the disk-shaped compact tension test, DC(T). For each test, the Subset Splitting DIC algorithm is used to determine the displacement field in a predefined region near the notch tip. Then, a given number of DEM simulations are performed on the same specimen. The DEM is used to simulate the fracture of asphalt concrete with a linear softening cohesive contact model, where fracture-related properties (e.g., maximum tensile force and maximum crack opening) are varied within a predefined range. The difference between DIC and DEM displacement fields for each set of fracture parameters is then computed and converted to a continuous function via multivariate Lagrange interpolation. Finally, we use a Newton-like optimization technique to minimize Lagrange multinomials, yielding a set of fracture parameters that minimizes the difference between the DEM and DIC displacement fields. The optimized set of fracture parameters from this nonlinear optimization procedure led to DEM results which are consistent with the experimental results for both SE(B) and DC(T) geometries.     

Size dependent response of large shell elements under in-plane tensile loading

Large-scale thin-walled structures with a low weight-to-stiffness ratio provide the means for cost and energy efficiency in structural design. However, the design of such structures for crash and impact resistance requires reliable FE simulations. Large shell elements are used in those simulations. Simulations require the knowledge of the true stress–strain response of the material until fracture initiation. Because of the size effects, local material relation determined with experiments is not applicable to large shell elements. Therefore, a numerical method is outlined to determine the effect of element size on the macroscopic response of large structural shell elements until fracture initiation. Macroscopic response is determined by introducing averaging unit into the numerical model over which volume averaged equivalent stress and plastic strain are evaluated. Three different stress states are considered in this investigation: uniaxial, plane strain and equi-biaxial tension. The results demonstrate that fracture strain is highly sensitive to size effects in uniaxial tension whereas in plane strain or equi-biaxial tension size effects are much weaker. In uniaxial and plane strain tension the fracture strain for large shell elements approaches the Swift diffuse necking condition.     

Cohesive fracture simulation and failure modes of FRP–concrete bonded interfaces

A work-of-fracture method using three-point bend beam (3PBB) specimen, commonly employed to determine the fracture energy of concrete, is adapted to evaluate the mode-I cohesive fracture of fiber reinforced plastic (FRP) composite–concrete adhesively bonded interfaces. In this study, a bilinear damage cohesive zone model (CZM) is used to simulate cohesive fracture of FRP–concrete bonded interfaces. The interface cohesive process damage model is proposed to simulate the adhesive–concrete interface debonding; while a tensile plastic damage model is used to account for the cohesive cracking of concrete near the bond line. The influences of the important interface parameters, such as the interface cohesive strength, concrete tensile strength, critical interface energy, and concrete fracture energy, on the interface failure modes and load-carrying capacity are discussed in detail through a numerical finite element parametric study. The results of numerical simulations indicate that there is a transition of the failure modes controlling the interface fracture process. Three failure modes in the mode-I fracture of FRP–concrete interface bond are identified: (1) complete adhesive–concrete interface debonding (a weak bond), (2) complete concrete cohesive cracking near the bond line (a strong bond), and (3) a combined failure of interface debonding and concrete cohesive cracking. With the change of interface parameters, the transition of failure modes from interface debonding to concrete cohesive cracking is captured, and such a transition cannot be revealed by using a conventional fracture mechanics-based approach, in which only an energy criterion for fracture is employed. The proposed cohesive damage models for the interface and concrete combined with the numerical finite element simulation can be used to analyze the interface fracture process, predict the load-carrying capacity and ductility, and optimize the interface design, and they can further shed new light on the interface failure modes and transition mechanism which emulate the practical application.     

高速冲击问题的离散元法数值模拟

以能量等效为原则,构建了新的求解弹塑性轴对称问题的离散元模型。通过引入断裂准则,使模型不仅可以应用于连续介质问题,还可以应用于从连续介质向非连续介质的动态转化过程。最后,通过模拟钢板受冲击载荷产生层裂的过程,说明模型在高速冲击问题中的适用性及其有效性。     

冲击载荷作用下路面结构的沉降及破坏特征

随着公路在导弹的无依托随机发射中扮演更加重要的角色，研究发射期间其在发射筒冲击载荷作用下的沉降及破坏特征对提高我国武器威慑力具有重要意义。选择低等级道路结构为研究对象，基于河南省、山东省、云南省及广东省的低等级道路统计信息，建立典型低等级道路结构的1/4简化计算模型，采用一种基于有限元和离散元耦合的显式数值分析方法CDEM，并引入塑性-局部化-破裂耦合的本构模型，实现冲击载荷作用下路面沉降的三维全时程模拟，显式地展现裂纹萌生和扩展过程，实现路面结构从连续状态到非连续状态的过渡。通过与落锤弯沉仪试验的监测结果对比验证道路简化模型的正确性及CDEM的计算精度，随后各选取一个沥青混凝土和水泥混凝土道路结构进行数值模拟，并从沉降量、破裂度、破坏特征等方面分析。计算结果表明:沉降量时程曲线与冲击载荷时程曲线变化趋势一致，沉降量在冲击载荷峰值点处达到最大值；破裂面主要产生于冲击载荷急剧增大时期，约占最终全部破裂面积的97%；破裂面大体分为路面结构层内的竖向破裂面和层间的水平破裂面，破坏类型包括拉伸破坏与剪切破坏；路面结构上部易产生环状破裂面，下部易产生径向破裂面。     

含冲击损伤高强中模碳纤维复合材料层压板压缩剩余强度分析与试验验证

论文以碳纤维复合材料层压板为研究对象,发展了一种模拟复合材料层压板冲击及冲击后压缩的一体化数值分析方法．基于Puck 失效准则和粘聚区模型描述层内损伤与层间损伤,分别采用基于断裂能的双线性型、函数型以及直接折减型等不同损伤折减方法构建了层内损伤预测与演化模型;建立了碳纤维复合材料冲击后压缩数值仿真模型,通过开展不同能量冲击后压缩试验,验证了所发展的数值分析方法的有效性;研究结果表明,采用Puck 失效准则和基于断裂能的双线性损伤演化模型预测冲击后压缩强度时具有较高精度．     

基于DEM-FEM耦合方法的海上风机结构冰激振动分析

冰载荷是海上风机在寒区安全运行的重要影响因素之一,由其引发的冰激振动给风机结构带来了严重的危害. 本文通过离散元(discrete element method, DEM)-有限元(finiteelement method, FEM)耦合方法建立了寒区单桩式风机结构的 冰激振动模型.采用具有粘结-破碎性能的球体离散单元描述平整海冰损伤破坏行为,采用梁单元和三角形平板壳单元构造带有抗冰锥体的单桩式风机有限元模型.采用DEM-FEM耦合方法模拟不同冰速、冰厚条件下单桩式风机与平整冰相互作用过程,并且与IEC规范和ISO标准经验公式对比验证该耦合模型计算冰载荷的准确性.对比风机塔筒顶端和基础顶端的位移和加速度响 应时程,定性地给出风机结构不同部位振动响应行为差异性.风机不同部位动力特性差异原因为风机结构独特结构特点:下部为大刚度桩基和上部为高柔度塔筒,使其动力特征表现为主从式结构特性. “主-从式结构”特征使得结构在复杂的冰载荷作用下,风机塔筒(子结构)和桩基(主结构)表现为不同的响应行为,风机不同部位振动周期和加速度谱两者出现差异. 本文研究成果为海上风机抗冰设计和疲劳分析提供了有益参考.      

基于细观建模的电子束预控弹体破裂机理数值研究

为研究电子束预控弹体的破裂机理。提出了基于电子束预控弹体细观几何特性的参数化建模方法,建立了含基体、熔融区、过渡区和空腔区弹体的三维有限元模型,采用LS-DYNA软件对典型弹体的爆炸驱动和破裂过程进行了数值模拟分析。结果表明:电子束预控弹体破裂过程可分为:弹体膨胀后的空腔区在环向拉应力作用下产生拉伸断裂、过渡区产生裂纹扩展和拉伸断裂以及空腔区底部基体在两侧拉应力和底部压应力作用下产生与弹体内壁法线呈45°的剪切破坏3个阶段。数值模拟结果与回收的破片截面形状和破坏模式吻合较好。研究结果对电子束预控弹体破片成型控制具有参考价值。     

Discrete fracture modeling of asphalt concrete

A heterogeneous fracture approach is presented for modeling asphalt concrete that is composed of solid inclusions and a viscous matrix, and is subjected to mode-I loading in the fracture test configuration. A heterogeneous fracture model, based on the discrete element method (DEM), is developed to investigate various fracture toughening mechanisms of asphalt materials using a high-resolution image processing technique. An energy-based bilinear cohesive zone model is used to model the crack initiation and propagation of materials, and is implemented as a user-defined model within the discrete element method. Experimental fracture tests are performed to investigate various fracture behavior of asphalt concrete and obtain material input parameters for numerical models. Also, bulk material properties are necessary for each material phase for heterogeneous numerical models; these properties are determined by uniaxial complex modulus tests and indirect tensile strength tests. The main objective of this study is to integrate the experimental tests and numerical models in order to better understand the fracture mechanisms of asphaltic heterogeneous materials. Experimental results and numerical simulations are compared at different test conditions with excellent agreement. The heterogeneous DEM fracture modeling approach has the potential capability to understand various crack mechanisms of quasi-brittle materials.