Multi-scale ballistic material modeling of cross-plied compliant composites

The open-literature material properties for fiber and polymeric matrix, unit-cell microstructural characteristics, atomic-level simulations and unit-cell based finite-element analyses are all used to construct a new continuum-type ballistic material model for 0°/90° cross-plied highly-oriented polyethylene fiber-based armor-grade composite laminates. The material model is formulated in such a way that it can be readily implemented into commercial finite-element programs like ANSYS/Autodyn [ANSYS/Autodyn version 11.0, User Documentation, Century Dynamics Inc. a subsidiary of ANSYS Inc. (2007)] and ABAQUS/Explicit [ABAQUS version 6.7, User Documentation, Dessault Systems, 2007] as a User Material Subroutine. Model validation included a series of transient non-linear dynamics simulations of the transverse impact of armor-grade composite laminates with two types of projectiles, which are next compared with their experimental counterparts. This comparison revealed that a reasonably good agreement is obtained between the experimental and the computational analyses with respect to: (a) the composite laminates’ capability, at different areal densities, to defeat the bullets with different impact velocities; (b) post-mortem spatial distribution of damage within the laminates; (c) the temporal evolution of composite armor laminate back-face bulging and delamination; and (d) the existence of three distinct penetration stages (i.e. an initial filament shearing/cutting dominated stage, an intermediate stage characterized by pronounced filament/matrix de-bonding/decohesion and the final stage associated with the extensive back-face delamination and bulging of the armor panel).     

Identification of CTOA and fracture process parameters by drop weight test and finite element simulation

This paper presents a new technique that is able to predict ductile fracture propagation occurrences in large metallic structures, by means of an appropriate application of the finite element modelling. This technique takes account of a cohesive zone in the vicinity of the crack tip, where a nodal release technique is implemented. Two parameters, governing the process zone of the material under investigation, have to be determined: the process zone dimension (named “Δ distance”) and the critical value of crack tip opening angle (CTOA). CTOA_C can be determined through an experimental laboratory procedure two specimen CTOA test (TSCT) that is already known and used by researchers who study fracture propagation on pipelines [Demofonti G, et al. Step-by-step procedure for the two specimen CTOA test. In: Proceedings of the Second International Conference on Pipeline Technology, Ostend, vol. II. 1995]. The second parameter required, Δ distance, is determined minimizing the differences of Finite Element results towards experimental data of an instrumented impact test (drop weight tear test). Some interesting improvements, concerning distinction between the initiation energy and the propagation energy accounted in TSCT procedure, are also discussed, in order to successfully extend its use to both high strength and high toughness steels.     

Role of Defects on Domain Wall Propagation in Magnetically Bistable Glass-Covered Microwires

We are trying to reveal the contribution of local defects on DW propagation in amorphous microwires. Consequently, we present a comparative study of single domain wall dynamics and local nucleation fields in Fe- and Co-rich amorphous glass-coated microwires. For this we modified experimental set-up based on the classical Sixtus?CTonks approach introducing additional pick-up coil in order evaluate defects contribution in domain wall propagation. Below some critical magnetic field, H _N , determined by the microwires inhomogeneities, an almost linear DW velocity,?v, dependence on magnetic field, H, is found. Quite fast DW propagation (v?till 3000?m/s at?H about 65?A/m) has been observed. When the applied magnetic field exceeds H _N , new reverse domains can be nucleated and consequently tandem remagnetization mechanism can be realized. The role of defects existing in magnetically bistable microwires is related with nucleation of new reversed domains.     

Failure and penetration response of borosilicate glass during multiple short-rod impact

In Anderson Jr CE, Orphal DL, Behner T, Templeton, DW [Failure and penetration response of borosilicate glass during short-rod impact. Int J Impact Eng 2009, doi:10.1016/ j.ijimpeng.2008.12.002.] it was demonstrated that the failure front (FF) produced by the penetration of a borosilicate glass target by a gold rod ceased to propagate a short time after the rod was fully eroded. This strongly suggests that progression of the FF is not described by a wave equation. Here it is shown that propagation of the FF is reinitiated if a second co-axial rod, spaced a distance from the first, impacts the glass at the bottom of the penetration channel. The experiments were performed in reverse ballistic mode with two short rods spaced apart. In some experiments both rods were gold; in other experiments, one rod was copper and the other gold. FF propagation was measured using high-speed photography; rod penetration was measured using multiple, independent flash X-rays. Much of the observed phenomenology can be modeled assuming that the rod, either first or second, “communicates” with the FF at a speed corresponding to the bulk sound speed of the undamaged glass.     

Modeling and characterization of dynamic failure of borosilicate glass under compression/shear loading

In this paper, we study the impact-induced dynamic failure of a borosilicate glass block using an integrated experimental/analytical approach. Previous experimental studies on dynamic failure of borosilicate glass have been reported by Nie et al. [Nie X, Chen WW, Sun X, Templeton DW. Dynamic failure of borosilicate glass under compression/shear loading – experiments. J Am Ceram Soc, in press.] using the split Hopkinson pressure bar (SHPB) technique. The damage growth patterns and stress histories have been reported for various glass specimen designs. In this study, we propose to use a continuum damage mechanics (CDM)-based constitutive model to describe the initial failure and subsequent stiffness reduction of glass. Explicit finite element analyses are used to simulate the glass specimen impact event. A maximum shear stress-based damage evolution law is used in describing the glass damage process under combined compression/shear loading. The impact test results are used in quantifying the critical shear stress for the borosilicate glass under examination. It is shown that with only two modeling parameters, reasonably good comparisons between the predicted and the experimentally measured failure maps can be obtained for various glass sample geometries. Comparisons between the predicted stress histories for different sample designs are also used as model validations.     

High velocity impact response of Kevlar-29/epoxy and 6061-T6 aluminum laminated panels

The high velocity impact response of composite laminated plates has been experimentally investigated using a nitrogen gas gun. Tests were undertaken on sandwich structures based on Kevlar-29 fiber/epoxy resin with different stacking sequence of 6061-T6 Al plates. Impact testing was conducted using cylindrical shape of 7.62 mm diameter steel projectile at a range of velocities (180–400 m/s) were investigated to achieve complete perforation of the target. The numerical parametric study of ballistic impact caused by same conditions in experimental work is undertaken to predict the ballistic limit velocity, energy absorbed by the target and comparison between simulation by using ANSYS Autodyn 3D v.12 software and experimental work and study the effects of shape of the projectile with different (4, 8 and 12 mm) thicknesses on ballistic limit velocity. The sequence of Al plate position (front, middle and back) inside laminate plates of composite specimen was also studied. The Al back stacking sequence plate for overall results obtained was the optimum structure to resist the impact loading.The results obtained hereby are in good agreement with the experimental (maximum error of 3.64%) data where it has been shown that these novel sandwich structures exhibit excellent energy absorbing characteristics under high velocity impact loading conditions. Hence it is considered suitable for applications of armor system.     

FEM analysis of localized deformation behaviour in pseudoelastic NiTi shape memory alloys

This paper takes into account the localized deformation behaviour of a pseudoelastic NiTi shape memory alloy with finite element method. A three‐dimensional micromechanical model has been developed, in which the difference between the elastic properties of austenite and martensite is considered. The model is implemented as User MATerial subroutine (UMAT) into ABAQUS. Then, a polycrystalline NiTi shape memory alloy with [1 1 1] texture under tensile loading is simulated. The main features of the propagation of the deformation band reported in literatures are captured in the simulation. It is also shown that the initiation and propagation of the deformation bands are strongly affected by the geometry of the specimens.     

Homogenised non-linear soft inclusion for simulation of impact damage in composite structures

This paper presents the development of a homogenised non-linear soft inclusion which captures the geometric and material non-linearity of impact damage zone loaded in tension and compression. The homogenised non-linear soft inclusion can present a conservative worst case damage zone or use experimental data to mimic the behaviour of a particular damage zone in a simple and computationally efficient way that can be used as a structural design tool for composite structures subjected to impact. The development of the non-linear soft inclusion, implemented in an ABAQUS/Explicit VUMAT, is presented at element and coupon level. The non-linear soft inclusion is validated against experimental coupon data and produces a conservative worst case estimate in all cases investigated.     

Dynamic buckling of elastic–plastic square tubes under axial impact—II: structural response

Dynamic elastic–plastic buckling of thin-walled square tubes is studied from the viewpoint of elastic–plastic stress wave propagation, which originates from an axial impact loading. The influence of the impact velocity and the striking mass on the development of the buckling shape is discussed when considering the transient deformation process. It is shown that the maximum load, which results from a high velocity impact load and occurs at t=0, is a function of the impact velocity and is related to the speed of the elastic–plastic stress waves propagating along the tube. The predictions for the initiation of buckling based on a numerical simulation of the axial impact of strain rate insensitive square tubes using the FE code ABAQUS show good agreement with the results from experiments on aluminium alloy tubes impacted at various initial velocities. A comparison between the buckling initiation in square tubes and geometrically equivalent circular tubes reveals differences in the response, which are attributed to the stress wave propagation phenomena and to the structural differences between the two structures.     

新型复合柱抗冲击试验研究及有限元分析

为改善结构柱抗冲击性能,提出最优截面形式,对不同截面的新型复合柱进行了抗冲击试验研究,并基于3种通用动力非线性有限元软件MSC.Marc、ABAQUS和ANSYS/LS-DYNA对撞击试验进行了数值模拟。试验结果表明新型复合柱的整体抗冲击性能较好;在相同外径和用钢量下,实心夹层钢管混凝土抗冲击性能最佳,空心夹层钢管混凝土抗冲击性能次之,内插双H型钢钢管混凝土稍差。ABAQUS建立的有限元模型和试验结果吻合最好,MSC.Marc次之,ANSYS/LS-DYNA稍差。研究结果可为新建大型铁路站房柱抗冲击选择最优截面提供依据,并为选取适合的有限元软件进行抗冲击分析提供参考。     

Ensuring primary realization of pressure unit in the vacuum range without typically utilized static expansion system

After the break-up of Czechoslovakia in 1993 the primary metrological standards authority almost fully remained in Slovakia. So we had to build our primary vacuum metrology anew. But we could do it using uptodate techniques that are more effective and they enable us to fulfil our task more effectively and in a shorter time.The steady demand to widen the lower pressure range and a steady demand for maximal accuracy collide with the financial possibilities open to a small country. So it is impossible to follow the classical pattern for the vacuum standards: liquid columns, static expansion system and dynamic expansion system. But it is possible to utilize contemporary progress in construction of piston manometers. This enables us to discard liquid column and static expansion techniques in favour of piston technique [Tesar J, Prazak D. Vacuum 2002; 67: 311, Tesar J, Prazak D, Krajicek Z. In:Proceedings of international symposium on Pressure and Vacuum, 2003. p. 169, Tesar J, Repa P, Prazak D, Krajicek Z, Peksa L. Vacuum 2004; 76: 491.]. Such system fulfils all the needs of a little industrial country and is economically acceptable.     

外包钢管加固RC柱抗冲击试验研究

为改善结构柱抗冲击性能,设计了RC柱和外包钢管加固RC柱在不同轴压比下的抗冲击试验,比较了两种类型柱抗冲击性能差异。试验通过DHR-9401落锤式冲击实验机进行,运用力传感器、高速摄影机和应变片完整记录了冲击力时程曲线、轴力时程曲线、跨中挠度时程曲线和试件特定位置的应变时程曲线。同时运用三种通用动力非线性有限元软件MSC.Marc、ABAQUS和ANSYS/LS-DYNA对冲击试验进行了数值模拟。研究结果表明外包钢管加固RC柱抗冲击性能要明显强于RC柱,证明外包钢管加固RC柱是一种可行的抗冲击加固措施。ABAQUS建立的有限元模型和试验结果吻合最好,MSC.Marc次之,ANSYS/LS-DYNA稍差。但ANSYS/LS-DYNA可较好的模拟RC柱冲击破坏全过程,其他两种软件则较难实现。研究结果可为大型铁路站房抗脱轨列车撞击提供加固措施,并为选取适合的有限元软件进行抗冲击分析提供参考。     

基于材料全程本构关系对延性断裂韧性的数值模拟方法与应用

利用ANSYS APDL编程语言,采用FAT (Finite-element-analysis Aided Testing)方法获取26NiCrMoV11-5汽轮机转子材料、A508-III压力容器钢和304不锈钢材料全程单轴拉伸本构关系曲线。借助ABAQUS软件中GTN (Gurson-Tvergaard-Needleman)损伤模型,采用有限元逆向法对3种材料漏斗试样单轴拉伸试验过程进行数值仿真,标定相应的GTN模型参数,获得了较好的结果。在此基础上将3种材料GTN模型参数用以模拟紧凑拉伸(Compact Tension,CT)试样的准静态裂纹扩展,基于载荷分离法实现了对26NiCrMoV11-5、A508-III和304不锈钢3种材料J阻力曲线的良好预测。     

On damage accumulations in the cyclic cohesive zone model for XFEM analysis of mixed-mode fatigue crack growth

Predicting mixed-mode fatigue crack propagation is an important and troublesome issue in structure assessment for decades. In the present paper an extended finite element method (XFEM) combined with a new cyclic cohesive zone model (CCZM) is introduced for simulating fatigue crack propagation under mixed-mode loading conditions, which has been implemented in the commercial general purpose software ABAQUS. The algorithm allows introducing a new crack surface at arbitrary locations and directions in a finite element mesh, without re-meshing. The cyclic cohesive zone model is based on the known S–N curves and Goodman diagram for metallic materials and validated by uniaxial tension results. Furthermore, the sensitivity of the model parameter is investigated for mixed-mode fatigue. The virtual crack closure technique has been extended to the cohesive zone model and proposed to calculate the energy release rate for the generalized Paris’ law. Finally, the crack propagation rate and direction under mixed-mode fatigue loading conditions are studied.     

Three-dimensional numerical simulations of dynamic fracture in silicon carbide reinforced aluminum

The three-point bending test by Kolsky-bar apparatus is a convenient technique to test the dynamic fracture properties of materials. This paper presents detailed three-dimensional finite element simulations of a silicon particle reinforced aluminum (SiC_p/Al) experiment (Li et al., [Proceedings of the US Army Symposium on Solid Mechanics]. In the simulations, the interaction between the input bar and the specimen is modeled by coupled boundary conditions. The material model includes large plastic deformations, strain-hardening and strain-rate hardening mechanisms. Furthermore, crack initiation and propagation processes are simulated by a cohesive element model. The simulation results quantitatively agree with the experimental measurements on three fronts: (1) the structural response of the specimen, (2) the time of unstable crack propagation, and (3) the local deformations at the crack-tip zone. The simulations reveal crack propagation characteristics, including crack-tip plastic deformation, crack front curving, and crack velocity profile. The effectiveness of Kolsky-bar type fracture tests is verified. It is shown that a rate-independent cohesive model can describe the complicated dynamic elastic–plastic fracture process in the SiC_p/Al material.     

一种修正的GFRP本构模型及其双层板结构抗冲击特性预测

玻璃纤维增强聚合物基材料(GFRP)是一种典型的各向异性复合材料,在冲击载荷作用下具有复杂的力学响应.Autodyn软件中现有的GFRP材料模型所采用的状态方程仅考虑了体应变的单一影响,忽略了各向异性材料中应力分量与应变分量之间的耦合效应.首先,在已有的解耦理论基础上,对现有模型中的状态方程和强化模型进行了修正,计算得...     

Characteristics of domain wall pinning and depinning in a three-terminal magnetic Y-junction

The characteristics of domain wall (DW) pinning and propagation in a three-terminal magnetic Y-junction were investigated, where the junction consisted of two input and one output wires. The output switching depends strongly on the junction angle (α). Junctions with high angles of α>9.5° lead to DW pinning at the junction, whereas junctions with low angles of α<9.5° have no DW pinning effect. At the critical angle of α = 9.5°, the Y-junction showed a multimode DW propagation, which was ascribed to a moderate transverse field effect.     

Experimental and numerical studies on dynamic crack growth in layered slate rock under wedge impact loads: part I—plane strain problem

The experimental and numerical investigations presented in this paper were carried out to determine the splitting forces and crack propagation scenarios of naturally bedded layered slate rock. Splitting loads were determined by impact splitting of regular‐sized slate blocks under plane strain test loading conditions, using a hydraulic actuator with a wedge‐shaped indenter. The mechanical properties of slate blocks required for numerical analyses were obtained from detailed experimental testing. The velocity of dynamic crack propagation in slate blocks under indenting wedge impact loading was determined using a series of strain gauge sensors. Numerical studies were carried out using ABAQUS, a general purpose, finite element analysis (FEA) program. Mode I dynamic crack propagation was simulated numerically by the gradual releasing of the restrained node on the symmetric plane of the specimens. Mode I stress intensity factors were computed for different crack lengths and the results were compared with the plane strain material fracture toughness obtained from earlier experiments/FEA. Very good agreement was obtained between analysis results and the measured fracture toughness value of slate, for the applied impact splitting load. Using the equation derived from a parametric study, of results obtained from the numerical analysis of different sizes of slate blocks, the maximum theoretical impact splitting force was determined using the plane strain fracture toughness value obtained from FEA. The difference between the loads obtained from the experimental studies and the derived empirical equation, varied between + 4.96% and −32.34%.     

Deterioration of fracture toughness of concrete under acid rain environment

Deterioration of fracture toughness of concrete caused by acid rain attack was investigated through laboratory-accelerated corrosion tests. Fracture toughness of the corroded concrete with different corrosion period was obtained by three point bending tests, and meanwhile the P-COMD curves and elastic modulus E of the corroded specimens with different corrosion degrees were also obtained. The experimental results showed that the deterioration of concrete surface is large and the corrosion depth increases with corrosion duration. The observation by scanning electron microscopy showed that the microstructures of crack tip are porous and loose due to acid corrosion, and many crisscross cracks appear at corrosion layer. Based on the experimental results, equivalent crack length in the fracture process zone and the effective fracture toughness K_IC were also obtained by theoretical analysis. Under the impact of acid corrosion, the concrete material physical and chemical properties changed significantly, which made the corresponding resistance stresses to crack propagation decrease. The curves of fracture toughness versus corrosion depth were obtained in this study.     

一种改进的内聚力损伤模型在复合材料层合板低速冲击损伤模拟中的应用

针对传统内聚力损伤模型(CZM)无法考虑层内裂纹对界面分层影响的缺点,提出了一种改进的适用于复合材料层合板低速冲击损伤模拟的CZM。通过对界面单元内聚力本构模型中的损伤起始准则进行修正,考虑了界面层相邻铺层内基体、纤维的损伤状态及应力分布对层间强度和分层扩展的影响。基于ABAQUS用户子程序VUMAT,结合本文模型及层合板失效判据,建立了模拟复合材料层合板在低速冲击作用下的渐进损伤过程的有限元模型,计算了不同铺层角度和材料属性的层合板在低速冲击作用下的损伤状态。通过数值模拟与试验结果的对比,验证了本文方法的精度及合理性。