长持时爆炸冲击波荷载作用下梁板组合结构的动力响应

为研究钢筋混凝土梁板组合结构在长持时远爆冲击波荷载作用下的动力响应及毁伤形态,通过实验获得了梁板组合结构的破坏形态和背爆面中心点位移变化。利用有限元软件对钢筋混凝土梁板组合结构的动态响应过程进行数值模拟研究,模拟得到的结构破坏现象与实验吻合较好。在此基础上,分析了梁板组合结构在相同冲量、不同峰值爆炸荷载作用下组合结构的动态响应和破坏过程,并结合挠跨比与破坏形态划分破坏模式。研究结果表明,相同冲量作用下,随着爆炸荷载峰值强度增加,梁板组合构件的破坏程度逐渐增加,破坏模式从弯曲破坏向弯剪联合破坏转换,最后呈现冲切破坏模式;组合构件中板部分发生破坏的时间早于交叉梁部分、破坏程度大于交叉梁。     

爆炸荷载作用下RC梁单边直剪破坏数值分析

本文通过建立数值模型,选择相应的参数信息得出数值计算结果对既有的爆炸荷载作用下RC简支梁的动力响应试验结果进行对比分析,认为所选数值模型及计算方法能够很好地与试验相匹配.据此,文中通过建立RC梁在一端作用有0.2kg的TNT炸药爆炸模型分析该梁的破坏模式、提取梁中不同部位的位移时程曲线得出RC梁破坏的具体形态及由局部破坏导致整根构件失效等一系列过程.爆压分布的测定对RC梁发生该种破坏类型时荷载的简化具有指导意义.通过截取梁两端支反力可知受到单边直剪破坏作用的RC梁端的剪切作用力远高于另一端.对另外三种近似工况建立数值模型,通过数值计算分析出RC梁构件发生单边直剪破坏风险的诸条件及要素.与此同时,确定了单边直剪破坏荷载类型简化形式.     

爆炸荷载作用下影响RC梁破坏形态的主要因素分析

在冲击、爆炸等强荷载作用下，钢筋混凝土结构通常是发生弯曲破坏。室内外试验结果表明，当动荷载峰值较大、作用时间较短（脉冲荷载）时，钢筋混凝土结构易发生脆性的剪切破坏，其主要原因是脉冲荷载的高频成份丰富和加载速率高等特点容易激发结构构件中的剪切变形和剪 应力。本文基于Timoshenko梁理论，提出了一种爆炸荷载作用下钢筋混凝土梁破坏形态的有限差分预报方法，利用该方法分析了荷载加载速率、截面高度、混凝土强度、配筋率等对钢筋混凝土破坏形态的影响规律，并指出，随着荷载升压段加载速率、截面高度、主筋配筋率的增加或混凝土强度的降低，梁的破坏形态从弯曲破坏转变为剪切破坏。     

钢筋混凝土拱的水下抗爆性能

为探究水下爆炸荷载作用下钢筋混凝土拱的动力响应特性和破坏特征,制作了两个钢筋混凝土拱试件,并开展了水下爆炸试验。试验分为拱外爆炸和拱内爆炸两组,采用10 g乳化炸药,试验时爆源距结构面最小距离为10 cm（起爆点位于拱结构正上方和正下方）,通过传感器记录爆炸试验中钢筋混凝土拱典型断面处的水压力及加速度时程曲线。基于Arbitrary Lagrange-Euler (ALE)算法,建立了空气-水-炸药-钢筋混凝土拱等多介质动态耦合作用模型,将数值模拟结果与试验结果对比,验证了数值方法的可靠性。采用验证后的数值模型进一步研究了拱外及拱内爆炸荷载作用下钢筋混凝土拱的动力响应差异。结果表明：相同炸药当量下,内部爆炸有更多的能量作用于混凝土拱,使结构的动力响应更强烈;外部爆炸下,拱顶、拱腰处产生较大裂缝;内部爆炸时,迎爆面裂缝数量明显增多,拱肩位置出现裂缝。钢筋混凝土拱形结构抵抗外部爆炸荷载的能力明显强于内部爆炸荷载。     

钢筋混凝土梁受弯破坏过程的细观数值模拟研究

在细观层次上将混凝土视为由粗骨料、硬化水泥砂浆及界面粘结带组成的三相非均质复合材料,能够较好地模拟混凝土加载时的裂缝扩展过程和分布规律。本文采用细观刚体弹簧元法模拟了钢筋混凝土梁的弯曲受力性能和破坏过程,加载方式采用两点对称加载。数值计算得到了钢筋混凝土梁的破坏形态和荷载一变形曲线,并分析了受力过程中纵向钢筋的应力变化。与试验结果对比表明,细观刚体弹簧元法可以有效模拟钢筋混凝土梁的裂缝开展过程、破坏形态和荷载一变形响应。     

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

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

ABAQUS混凝土损伤塑性模型中损伤因子的率相关性及实现方法

ABAQUS程序中最常用的混凝土损伤塑性（concrete damage plasticity, CDP）模型无法实现损伤因子与应变率相关。为了准确描述混凝土材料在高应变率下的损伤特性,基于CDP模型定义了新的应变率场变量,编制了VUSDFLD用户子程序,开发了能够考虑损伤因子率相关性的改进的CDP（modified CDP,MCDP）模型。MCDP模型采用能量法求解混凝土拉压损伤因子,主求解程序能够随着应变率场变量的变化而自动更新不同应变率对应的损伤参数,计算得到的混凝土单轴静态加载结果与CDP模型吻合较好。MCDP模型对高应变率下动态压缩性能的模拟结果表明:混凝土材料在不同应变率下的拉压损伤对其动态力学性能有显著影响,编制的VUSDFLD子程序和MCDP模型能够有效地解决损伤应变率相关的模拟难题,可以准确地模拟爆炸荷载作用下钢筋混凝土梁的动态响应,为预测爆炸冲击等强动载作用下混凝土结构的响应和破坏提供了更可靠的技术途径。     

冲击荷载与火灾联合作用下SFRC梁的力学行为

为了探究冲击荷载与火灾联合作用下钢纤维混凝土(steel fiber reinforced concrete, SFRC)梁的力学性能，联合应用高性能落锤试验系统、四点弯曲实验装置与装配式电炉开展了4根SFRC梁的冲击实验与高温恒载实验，观察了其破坏模式并记录了跨中位移和钢筋应变的时程曲线，探讨了冲击损伤SFRC梁的抗火性能。此外，在实验研究的基础上，考虑材料的应变率强化效应及温度软化效应，建立数值模型，首先对梁进行冲击加载模拟，并以冲击模拟结果为初始状态，采用热-力“顺序”耦合方法，对冲击加载与高温恒载联合作用下SFRC梁的力学行为进行了三维宏观有限元数值模拟。同时，考虑混凝土内部结构非均质性的影响，采用类似步骤，开展了细观模拟。宏/细观模拟结果与实验结果的良好吻合验证了本文数值方法的合理性与有效性，并体现了细观方法的优越性。研究发现，冲击能量较小时，SFRC梁在冲击荷载作用下，尽管局部混凝土开裂，梁整体残余变形较小，抗火性能有一定程度的下降；随着钢纤维掺量增大，混凝土基体抗剪强度增大，SFRC梁在冲击荷载作用下的开裂形态由弯剪裂缝并存向以弯曲裂缝为主转变；冲击损伤SFRC梁在高温恒载作用下裂缝分布较为集中，且发生脆性破坏。     

爆炸荷载下基于细观建模的素/钢筋混凝土板破坏模式

为了获得爆炸荷载下细观结构对素/钢筋混凝土板的影响，采用随机骨料投放建立了素/钢筋混凝土板细观模型。利用LS-DYNA对基于细观建模的钢筋混凝土板进行爆炸荷载作用下的数值模拟，通过与实验以及均质建模方法进行比较，验证了细观建模方法的准确性。进而研究了基于细观建模的素/钢筋混凝土板在不同爆炸荷载下的结构响应，获得了素/钢筋混凝土板的响应过程和破坏模式。结果表明:在低药量（1、2 kg）爆炸荷载下，细观结构对素/钢筋混凝土板的影响较小，其破坏模式以纵横塑性铰线破坏为主，药量越大，铰线越多；在高药量（5、10和15 kg）爆炸荷载下，细观结构对素/钢筋混凝土板的影响较大，与均质模型相比存在较大差异，细观素/钢筋混凝土板以爆坑为中心，产生环向与径向裂纹，药量越大，圆坑越大，裂纹越多，板局部破坏越严重。     

钢板-混凝土组合加固钢筋混凝土梁的非线性有限元分析

以10根采用钢板-混凝土组合加固技术加固的钢筋混凝土梁的试验为基础,建立了钢板-混凝土组合加固钢筋混凝土梁的非线性有限元模型.沿用了钢板-混凝土组合梁有限元分析中栓钉的荷载-滑移曲线及断裂模型,提出了新老混凝土界面的剪切-滑移模型,可以较好地模拟界面受力性能.采用有限元软件ABAQUS,模拟了一次加固、不同损伤加固及持载加固等加载情况下加固构件受力全过程和受弯破坏以及新老混凝土剥离破坏等破坏形态下加固试件的受力性能.在验证了有限元模型的准确性后,分析了钢板厚度、加固部分高度及长度等参数对加固构件承载能力及刚度的影响.最后提出了几点对钢板-混凝土组合加固RC梁的设计建议.     

GEOMETRICALLY NONLINEAR FINITE ELEMENT MODEL OF SPATIAL THIN-WALLED BEAMS WITH GENERAL OPEN CROSS SECTION

Based on the theory of Timoshenko and thin-walled beams, a new finite element model of spatial thin-walled beams with general open cross sections is presented in the paper, in which several factors are included such as lateral shear deformation, warp generated by nonuni- form torsion and second-order shear stress, coupling of flexure and torsion, and large displacement with small strain. With an additional internal node in the element, the element stiffness matrix is deduced by incremental virtual work in updated Lagrangian （UL） formulation. Numerical examples demonstrate that the presented model well describes the geometrically nonlinear property of spatial thin-walled beams.     

燃气爆炸作用下蒸压加气混凝土砌体墙的加固性能

为研究燃气爆炸作用下蒸压加气混凝土砌体墙的加固性能，基于有限元软件LS-DYNA，建立了砌体墙简化数值模型，分析了GB 50779-2012 石油化工控制室抗爆设计规范中建议的荷载作用下砌体墙高度和厚度的影响，对比了玄武岩纤维(basalt fiber reinforced plastic, BFRP)布与喷涂式聚脲对蒸压加气混凝土单向砌体墙的加固效果，并以防止砌体墙倒塌为设计目标，给出了加固建议。研究表明，本文中建立的简化数值模型能较好地模拟燃气爆炸作用下蒸压加气混凝土砌体墙的变形和破坏模式，计算结果与试验吻合良好；《规范》建议荷载作用下，未加固砌体墙以弯曲破坏为主，随着墙体高度增加，破坏模式由弯曲破坏向剪切破坏转变；BFRP布条加固可以有效提高墙体抗弯刚度和压拱效应，而聚脲涂层加固对抗弯刚度提高有限但墙体拉拱效应明显，二者均能显著提高墙体抗爆性能；加固墙体均发生弯曲破坏，BFRP布条材料的断裂一般发生在墙体位移最大处，而聚脲涂层材料的断裂发生在跨端边界处。     

Shear stress analysis in engineering beams using deplanation field of special 2-D finite elements

This paper deals with the 2-D finite element shear stress analysis in beams, loaded by bending with shear and St. Venant’s torsion. The properties of these finite elements, like stiffness matrices as well as load vectors, are derived on the basis of their axial nodal displacements, e.g. by warping field. Proposed finite elements enable stress analysis independently of both cross-sectional member shape and material properties. Stiffness matrices and load vectors are derived for several finite element types. Material is assumed to be isotropic and linear elastic. For justification of the proposed stress analysis procedure, some examples are presented.     

爆炸荷载下钢筋混凝土梁的变形和破坏

利用爆炸压力模拟器进行钢筋混凝土简支梁爆炸冲击实验,详细介绍了实验设计,通过实验系统分析了钢筋混凝土(RC)梁变形破坏特征以及钢筋作用机理和对变形破坏的影响,并建立了RC梁的分离式有限元模型,利用LS-DYNA分析了实验过程,对计算结果与实验结果进行了比较,分析了误差产生的原因,得到了爆炸冲击荷载作用下RC梁的损伤破坏特征和机理,可为毁伤评估和结构抗爆设计提供参考。     

An analysis of the shear failure of rigid-linear hardening beams under impulsive loading

A theoretical rigid-plastic analysis for the dynamic shear failure of beams under impulsive loading is presented when using a travelling plastic shear hinge model which takes into account material strain hardening. The maximum dynamic shear strain and shear strain-rate can be predicted in addition to the permanent transverse deflections and other parameters. The conditions for the three modes of shear failure, i.e., excess deflection failure, excess shear strain failure and adiabatic shear failure are analyzed. The special case of an infinitesimally small plastic zone is discussed and compared with Nonaka's solution for a rigid, perfectly plastic material. The results can also be generalized to examine the dynamic response of fibre-reinforced beams.     

Novel Technique for Static and Dynamic Shear Testing of Ti6Al4V Sheet

Few shear test techniques exist that cover the range of strain rates from static to dynamic. In this work, a novel specimen geometry is presented that can be used for the characterisation of the shear behaviour of sheet metals over a wide range of strain rates using traditional tensile test devices. The main objectives during the development of the shear specimen have been 1) obtaining a homogeneous stress state with low stress triaxiality in the zone of the specimen subjected to shear and 2) appropriateness for dynamic testing. Additionally, avoiding premature specimen failure due to edge effects was aimed at. Most dimensional and practical constraints arose from the dynamic test in which the specimen is loaded by mechanical waves in a split Hopkinson tensile bar device. Design of the specimen geometry is based on finite element simulations using ABAQUS/Explicit. The behaviour of the specimen is compared with the more commonly used simple shear specimen with clamped grips. Advantages of the new technique are shown. The technique is applied to Ti6Al4V sheet. During the high strain rate experiments high speed photography and digital image correlation are used to obtain the local shear strain in the specimen. Comparison of experimental and numerical results shows good correspondence.     

Elastodynamic Green’s function for reinforced concrete beams

A semi-analytical solution procedure for three dimensional wave propagation in reinforced concrete (RC) beams has been presented in this paper. Elastodynamic Green’s function has been derived by employing the compatibility conditions and utilizing the symmetry conditions at the loaded cross section. Numerical procedure developed for the Green’s function has been validated using results available in the literature for an infinite laminated composite plate. Three-dimensional wave propagation analysis has been performed for reinforced concrete beam sections of T and L shapes which are common forms of structural elements. Steel reinforcement has been modeled in the finite element mesh. Effect of corrosion has also been included in the finite element model. Green’s function for reinforced concrete sections affected by corrosion of steel unit normalized frequency has been evaluated for illustration. Accuracy of the solution technique has been evaluated in terms of the percentage error in energy balance between the input energy of the applied unit load and the output energy carried by the propagating wave modes. The percentage error has been found to be negligible in all the cases considered here. A simple and accurate numerical method has been presented here as a tool to evaluate Green’s function for RC beams and can be used to detect corrosion.     

方形钢筋混凝土板的近场抗爆性能

为了探讨钢筋混凝土板在爆炸载荷作用下的抗爆性能,对方形钢筋混凝土板在单向支撑条件下进行了近场爆炸加载实验,实验中采取TNT装药对钢筋混凝土板进行加载.并利用AUTODYN软件采用流固耦合算法,建立了混凝土和钢筋三维分离式实体模型,对钢筋混凝土板的动态响应过程进行数值模拟,且考虑应变率对钢筋和混凝土材料的动态本构特性的影...     

Problems in determining the elastic strain energy function for rubber

Lightly crosslinked natural rubber can be stretched by 600% or more, and recovers almost completely. It is often regarded as a model highly elastic material and characterized by a strain energy function to describe its stress-strain behavior under various types of deformation. A number of such functions have been proposed; some of them appear in current finite element programs. They are usually validated by comparison with measured stress-strain relations by Treloar [7] [L.R.G. Treloar, Stress-strain data for vulcanized rubber under various types of deformation, Trans. Faraday Soc. 40 (1944) 59-70] and Jones and Treloar [15] [D.F. Jones, L.R.G. Treloar, The properties of rubber in pure homogeneous strain, J. Phys. D Appl. Phys. 8 (1975) 1285-1304]. But Treloar pointed out that the relations at high strains became markedly irreversible, and he did not assign a strain energy function for strains greater than about 300%. Rivlin's universal relation between torsional stiffness and tensile stress [14] [R.S. Rivlin, Large elastic deformations of isotropic materials. Part V1: further results in the theory of torsion, shear and flexure, Philos. Trans. R. Soc. A 243 (1949) 251-288] is applied here to show that a typical elastic solid cannot be described by any strain energy function at strains greater than about 300%. Elastic strain energy functions for higher strains, or for other rubbery materials, are thus of doubtful value unless evidence for reversibility of stress-strain relations is adduced or the applicability of a strain energy function is demonstrated.     

40Cr材料动态起裂韧性KId（）的实验测试

描述了利用Ｈｏｐｋｉｎｓｏｎ压杆技术加载三点弯曲试样测试４０Ｃｒ，材料动态起裂韧性ＫＩｄ（）的试验方法。试样上的动态载荷历程由Ｈｏｐｋｉｎｓｏｎ杆直接测得，并分别代入动态有限元程序及近似公式求得动态应力强度因子历史；由贴在试样裂尖附近的应变片确定起裂时间，最终确定起裂时的动态应力强度因子值，即动态起裂韧性ＫＩｄ（）。试验结果表明:利用Ｈｏｐｋｉｎｓｏｎ压杆技术加载三点弯曲试样测试材料动态起裂韧性的方法是可行的，起裂时，动态有限元的位移法、应力法及近似公式法求得的动态应力强度因子值比较吻合；在本文的载荷速率下，４０Ｃｒ材料动态起裂韧性ＫＩｄ（）与准静态裂韧性ＫＩｄ（）相比，降低了约２８％。