An experimental and numerical study of the dynamic response of composites under impact at low temperatures

An analytical method was applied to investigate the effects of strain rate, low temperature, and impact direction on the energy absorption of composite specimens. The performance of an analytical model was found to be adequate for modeling composites before failure initiation. Results indicated that impact performance of composites is affected over the range of temperatures. Failure mechanism was changed from matrix cracking at room temperature to delamination and fiber breakage at low temperatures. Moreover, it was shown that about 70% of total energy absorbed by a specimen was used to destroy the composite under different types of failure mechanisms.     

A wave propagation model for the high velocity impact response of a composite sandwich panel

A solution methodology to predict the residual velocity of a hemispherical-nose cylindrical projectile impacting a composite sandwich panel at high velocity is presented. The term high velocity impact is used to describe impact scenarios where the projectile perforates the panel and exits with a residual velocity. The solution is derived from a wave propagation model involving deformation and failure of facesheets, through-thickness propagation of shock waves in the core, and through-thickness core shear failure. Equations of motion for the projectile and effective masses of the facesheets and core as the shock waves travel through sandwich panel are derived using Lagrangian mechanics. The analytical approach is mechanistic involving no detail account of progressive damage due to delamination and debonding but changes in the load-bearing resistance of the sandwich panel due to failure and complete loss of resistance from the facesheets and core during projectile penetration. The predicted transient deflection and velocity of the projectile and sandwich panel compared fairly well with results from finite element analysis. Analytical predictions of the projectile residual velocities were also found to be in good agreement with experimental data.     

Preliminary study of counter impact with two-stage light gas gun using electrothermal–chemical gun technology

Recently, hypervelocity impacts between space structures and space debris have been brought to attention with the advance of space development. The impact velocity at low earth orbit (LEO) is up to 15 km/s. Such impact velocity cannot be attained by the current two-stage light gas gun (TSLGG). Therefore we have a plan of counter impacts that used two sets of TSLGGs in order to raise impact velocity. There are three problems in order to realize the counter impacts. They are launch timing, recovery of target and alignment of gun axis. In this research, we considered launch timing of TSLGG among them. It is important to control variation of delay time from propellant ignition to projectile launch in the counter impacts. In a conventional TSLGG, combustion rate changes every time and impact position will vary widely. In order to stabilize combustion rate, we developed a new TSLGG that used the technology of electrothermal–chemical (ETC) gun for its 1st stage. In a conventional TSLGG, the variation of delay time is 20–30 ms. On the other hand, in a new TSLGG, we have succeeded in suppressing the delay time variation less than 200 μs. It is concluded that the TSLGG based on the technology of ETC gun is effective to control the delay time in this paper.     

冲击荷载下单层球面网壳动力响应分析与试验研究

为分析带下部支承结构的Kiewitt-6型单层球面网壳在冲击荷载下动力响应,在ANSYS/LS-DYNA中建立钢管柱支承的单层球面网壳数值分析模型,据冲击响应特点,总结四种响应模式,研究冲击能量、冲击位置、环梁刚度对上部网壳动力响应影响。进行钢管柱支承的单层球壳模型冲击试验,测量、分析结构的动应力、动位移及加速度,研究冲击柱破坏模式。结果表明,四种响应模式以冲击柱破坏模式(轻微损伤、局部凹陷、压弯破坏、剪切破坏)为典型特征;除响应模式4,网壳动力响应随冲击能量增大而增大;柱中为最不利冲击位置;环梁刚度增大,网壳动力响应减小;有限元分析结果与实测结果吻合较好,验证数值计算方法的正确性。     

Analytical and numerical investigation of polyurea layered aluminium plates subjected to high velocity projectile impact

This paper discusses on the penetration of high velocity projectiles through aluminium–polyurea composite layered plate systems. An analytical model has been proposed to predict the residual velocity of aluminium–polyurea composite plates, and validated with both experimental and numerical investigations. Full metal jacket (FMJ) projectiles (5.56 mm × 45 mm), corresponding to NATO standard SS109, were fired at the aluminium–polyurea composite layered plate systems from a distance of 10.0 m at a fixed velocity of 945 m/s. Four different composite plate configurations were used with thicknesses varying from 16 to 34 mm. Each configuration consisted of six specimens. Residual velocities for each individual test were recorded. Numerical simulations of the penetration process have been performed using advanced finite element code LS-DYNA®. The well-established Johnson–Cook and Mooney–Rivlin material models were used to represent the stress–strain behaviour of aluminium and polyurea in the numerical analysis. The analytical and numerical models provided good approximations for the residual velocities measured during the experimental tests. Polyurea layers contributed positively towards the reduction of residual velocity of the projectile in composite plate systems. In addition, ballistic limit curves for different composite systems have been established based on the validated models. As the results showed that polyurea contributes positively towards the reduction of residual velocity of projectiles, the findings of this study can be effectively used for the similar applications in future armour industry.     

钢管混凝土梁在侧向冲击荷载作用下动力响应的试验研究和数值模拟

为研究不同轴压比对钢管混凝土柱侧向冲击性能的影响,该文对已有的冲击试验进行了验证。在验证模型基础上建立了不同轴压比的钢管混凝土柱在侧向冲击作用下的动力有限元模型,计算得到连续变化轴压比下钢管混凝土的抗冲击性能参数。结果表明：轴压比对钢管混凝土的侧向抗冲击性能有显著影响。在轴压比小于0.7时,轴力对钢管混凝土构件侧向抗冲击能力有提高作用,当轴压比大于0.7时,轴力对钢管混凝土侧向抗冲击能力有削弱作用。     

Formation of amorphous and nanocrystalline phases in high velocity oxy-fuel thermally sprayed a Fe–Cr–Si–B–Mn alloy

High velocity oxy-fuel (HVOF) thermal spray was used to deposit a Fe–Cr–Si–B alloy coating onto stainless steel (1Cr18Ni9Ti) substrate. Microstructures of the powder and the coating were investigated by X-ray diffraction (XRD), scanning election microscopy (SEM), transmission election microscopy (TEM) and differential scanning calorimeter (DSC). The coating had layered morphologies due to the deposition and solidification of successive molten or half-molten splats. The microstructures of the coating consisted of a Fe–Cr-rich matrix and several kinds of borides. The Fe–Cr-rich matrix contained both amorphous phase and nanocrystalline grains with a size of 10–50 nm. The crystallization temperature of the amorphous phase was about 605 °C. The formation of the amorphous phase was attributed to the high cooling rates of molten droplets and the proper powder compositions by effective addition of Cr, Mn, Si and B. The nanocrystalline grains could result from crystallization in amorphous region or interface of the amorphous phase and borides by homogeneous and heterogeneous nucleation.     

Crystallization of an amorphous B–C–N precursor with a Li–B–N catalyst at high pressures and temperatures

An orthorhombic B–C–N compound was synthesized using an amorphous B–C–N precursor and a Li–B–N catalyst at 6 GPa and 1773 K. The results of energy dispersive spectrometry and electronic energy loss spectrometry suggest a stoichiometry of B:C:N = 1:3.3:1. In addition, the Li–B–N catalyst improves the crystallizations of the B–C–N compound, graphite and BN and therefore might be a profitable catalyst in ultrahigh pressure experiments.     

An experimental investigation on the low velocity impact response of composite plates repaired by VARIM and hand lay-up processes

This paper presents results of an experimental investigation on the impact response of repaired and unrepaired glass/epoxy composite plates. Repaired samples were prepared by two different manufacturing methods; vacuum assisted resin infusion process and hand lay-up technique. In order to compare impact response of the repaired and unrepaired samples a number of single impact tests were performed under various impact energies. Damage process of the samples is analyzed from cross-examining load–deflection curves and damaged specimens. From the visual inspection, for the impacted side of the samples, it is noted that the main damage modes for repaired samples are matrix and fiber cracks around point of impact and delaminations while severe matrix cracks expanded through fiber directions are the dominant damage mode for unrepaired samples. At the back surfaces, delaminations and fiber–matrix debonding oriented in the fiber directions are observed for unrepaired samples. However, for repaired samples the fiber fractures through repair line as well as the delaminations become dominant modes. For a reasoning justification in discussing impact test results, interlaminar fracture toughness (Mode I and Mode II) and flexural tests for repaired and non-repaired samples were also conducted.     

Energy absorption capacity of a sustainable Ultra-High Performance Fibre Reinforced Concrete (UHPFRC) in quasi-static mode and under high velocity projectile impact

This paper investigates the energy absorption capacity of a sustainable Ultra-High Performance Fibre Reinforced Concrete (UHPFRC) in quasi-static mode and under high velocity projectile impact. The design of the sustainable concrete mixtures aims on achieving a densely compacted cementitious matrix with a relatively low binder amount, employing the modified Andreasen & Andersen particle packing model. The experiments on UHPFRC are performed using a 4-point bending test and high velocity projectile impact tests. The obtained results show that although the utilization of hybrid steel fibre enhances the mechanical properties of the developed UHPFRC, the application of fibres with hooked ends is crucial in improving the energy absorption capacity of the sustainable UHPFRC in quasi-static mode. However, under high velocity projectile impact, the UHPFRC mixture with hybrid fibres shows a much better energy absorption capacity than the one with hooked steel fibres only, particularly in resisting the scabbing at the rear surface. The intrinsic mechanisms for the energy absorption capacity of the sustainable UHPFRC in quasi-static mode and under high velocity projectile impact are studied and analysed.     

Effect of the crab waist and of the micro-beta on the beam–beam instability

We calculate the luminosity and the strengths of the beam–beam resonances for colliders with large horizontal crossing angles. Achievable luminosities of such colliders can reach high values provided that the number of particles in colliding beams can be increased while the vertical β-function can be decreased till the mini-beta range. The crab waist option of the optics in the interaction region decreases (or, even vanishes) the strengths of two-dimensional betatron weak–strong beam–beam resonances and of their synchro-betatron satellites provided that β-functions at the interaction point can be decreased till the micro-beta range. This can help to increase the achievable value of the collider luminosity.     

An implicit boundary element solution with consistent linearization for free surface flows and non‐linear fluid–structure interaction of floating bodies

In this work, a new comprehensive method has been developed which enables the solution of large, non‐linear motions of rigid bodies in a fluid with a free surface. The application of the modern Eulerian–Lagrangian approach has been translated into an implicit time‐integration formulation, a development which enables the use of larger time steps (where accuracy requirements allow it). Novel features of this project include: (1) an implicit formulation of the rigid‐body motion in a fluid with a free surface valid for both two or three dimensions and several moving bodies; (2) a complete formulation and solution of the initial conditions; (3) a fully consistent (exact) linearization for free surface flows valid for any boundary elements such that optimal convergence properties are obtained when using a Newton–Raphson solver. The proposed framework has been completed with details on implementation issues referring mainly to the computation of the complete initial conditions and the consistent linearization of the formulation for free surface flows. The second part of the paper demonstrates the mathematical and numerical formulation through numerical results simulating large free surface flows and non‐linear fluid structure interaction. The implicit formulation using a fully consistent linearization based on the boundary element method and the generalized trapezoidal rule has been applied to the solution of free surface flows for the evolution of a triangular wave, the generation of tsunamis and the propagation of a wave up to overturning. Fluid–structure interaction examples include the free and forced motion of a circular cylinder and the sway, heave and roll motion of a U‐shaped body in a tank with a flap wave generator. The presented examples demonstrate the applicability and performance of the implicit scheme with consistent linearization. Copyright © 2001 John Wiley & Sons. Ltd.     

An experimental and numerical study of the dynamic behaviour of a counter-flow evaporator

The dynamic behaviour of a counter-flow, water-heated evaporator is studied experimentally and numerically. The frequency distribution of the random oscillations of the mixture-vapour transition point and the superheat temperature at the exit of the evaporator is obtained for steady operation of the system. These oscillations are well correlated. The transition point movement is found to cause fluctuations in the refrigerant temperature over 1 m downstream of its range of motion. Step changes in the refrigerant flow rate and the heating water flow rate demonstrate the non-linear characteristics of the evaporator where the time constants for step increases and step decreases of the same magnitude differ significantly. The distributed model predicts the variation of the superheat temperature and the evaporator pressure following step changes in the inputs with good accuracy.     

A Chebyshev collocation multidomain method to solve the Reissner–Mindlin equations for the transient response of an anisotropic plate subjected to impact

The transient response of an anisotropic rectangular plate subjected to impact is described through a Chebyshev collocation multidomain discretization of the Reissner–Mindlin plate equations. The trapezoidal rule is used for time-integration. The spatial collocation derivative operators are represented by matrices, and the subdomains are patched by natural and essential conditions. At each time level the resulting governing matrix equation is reduced by two consecutive block Gaussian eliminations, so that an equation for the variables at the subdomain corners has to be solved. Back-substitution gives the variables at all other collocation points. The time history as represented by computed contour plots has been compared with analytical results and with photos produced by holographic interferometry. The agreements are satisfactory. © 1997 John Wiley & Sons, Ltd.     

An experimental and theoretical study on the prediction of forming limit diagrams using new BBC yield criteria and M–K analysis

In the present paper, a comprehensive study on the prediction of forming limit diagrams (FLDs) for an AA3003-O aluminium alloy is developed theoretically and experimentally. For obtaining the experimental FLDs, an out-of-plane formability test was performed based on the technique proposed by Ozturk and Lee [F. Ozturk, D. Lee, J. Mater. Process. Technol. 170 (2005) 247–253]. The classical Marciniak–Kuczynski (M–K) model and some new yield criteria are utilized to simulate the necking phenomenon and calculate the limit strains theoretically. The employed yield functions are: the BBC2000, BBC2002, and BBC2003 yield criteria proposed by Banabic et al. [D. Banabic, S.D. Comsa, T. Balan, in: Proceedings of the Cold Metal Forming 2000 Conference, Cluj-Napoca, 2000, p. 217; D. Banabic, T. Kuwabara, T. Balan, D.S. Comsa, D. Julean, Int. J. Mech. Sci. 45 (2003) 797–811; D. Banabic, H. Aretz, D.S. Comsa, L. Paraianu, Int. J. Plast. 21 (2005) 493–512]. To calibrate and determine each particular coefficients of performed yield functions an appropriate error-function is defined and minimized by a Newton algorithm. To compare the calculated yield stresses and r-values with experimental data a relative root mean square deviation method presented by Leacock [Alan G. Leacock, J. Mech. Phys. Solids 54 (2006) 425–444] is used. Work-hardening effects on the FLD are analyzed by using Swift and Voce hardening laws. The effect of yield surface on the prediction of numerical FLDs and the number of experimental anisotropy parameters on the accuracy of yield functions are also studied.     

An experimental study on the shattering behavior of a high strength armour steel under blast and long rod penetrator impact

Armour steel plates with drilled holes are filled with explosive and subjected to single or multiple blasts to induce shattering. The critical diameter for shattering under explosive detonation, found from the experiments matches closely with that of the long rod penetrator impacted plates. The damage pattern and fracture surface of the tested samples under blast effect have been compared with plates impacted with long rod penetrators for matching the observed shattering behavior. Difference in the shattering behavior of the armour steel subjected to single and multiple blasts has also been presented.     

Application of the higher‐order Cauchy–Born rule in mesh‐free continuum and multiscale simulation of carbon nanotubes

This paper investigates the application of a recently proposed higher‐order Cauchy–Born rule in the continuum simulation and multiscale analysis of carbon nanotubes (CNTs). A mesh‐free computational framework is developed to implement the numerical computation of the hyper‐elastic constitutive model that is derived from the higher‐order Cauchy–Born rule. The numerical computation reveals that the buckling pattern of a single‐walled carbon nanotube (SWCNT) can be accurately displayed by taking into consideration the second‐order deformation gradient, and fewer mesh‐free nodes can provide a good simulation of homogeneous deformation. The bridging domain method is employed to couple the developed mesh‐free method and the atomistic simulation. The coupling method is used to simulate the bending buckling of an SWCNT and the tensile failure of an SWCNT with a single‐atom vacancy defect, and good computational results are obtained. Copyright © 2008 John Wiley & Sons, Ltd.     

Investigation of fluid–structure interactions using a velocity‐linked P2/P1 finite element method and the generalized‐ α method

A velocity‐linked algorithm for solving unsteady fluid–structure interaction (FSI) problems in a fully coupled manner is developed using the arbitrary Lagrangian–Eulerian method. The P2/P1 finite element is used to spatially discretize the incompressible Navier–Stokes equations and structural equations, and the generalized‐ α method is adopted for temporal discretization. Common velocity variables are employed at the fluid–structure interface for the strong coupling of both equations. Because of the velocity‐linked formulation, kinematic compatibility is automatically satisfied and forcing terms do not need to be calculated explicitly. Both the numerical stability and the convergence characteristics of an iterative solver for the coupled algorithm are investigated by solving the FSI problem of flexible tube flows. It is noteworthy that the generalized‐ α method with small damping is free from unstable velocity fields. However, the convergence characteristics of the coupled system deteriorate greatly for certain Poisson's ratios so that direct solvers are essential for these cases. Furthermore, the proposed method is shown to clearly display the advantage of considering FSI in the simulation of flexible tube flows, while enabling much larger time‐steps than those adopted in some previous studies. This is possible through the strong coupling of the fluid and structural equations by employing common primitive variables. Copyright © 2012 John Wiley & Sons, Ltd.     

A molecular dynamics study on the impact of defects and functionalization on the Young modulus of boron–nitride nanotubes

In this article, we examine the Young modulus of (6,m) boron–nitride nanotubes with vacancy and functionalization defects. We employ molecular dynamics simulations using the Parrinello–Rahman approach. To this end, all systems are modeled with a reactive many-body bond order Tersoff potential with parameters due to Matsunaga et al. [K. Matsunaga, Y. Iwamoto, Molecular dynamics study of atomic structure and diffusion behavior in amorphous silicon nitride containing boron, Journal of American Ceramics Society 84(10) (2001) 2213–2219], which is able to accurately describe covalent bonding. We apply external stress to periodically repeated tubes in vacuum and derive stress–strain curves for various tensile loads at standard temperature and pressure. In addition to the Young modulus, we study visualized stress-per-atom snapshots of the simulation runs. Our results show that the decrease in Young modulus with increasing defect concentration is independent of the chirality of the tube for vacancy defects. Also, we observe that functionalization does not weaken the tube. There is even indication of a relative strengthening for armchair types.     

Theoretical and experimental study on characteristics of a novel silica gel–water chiller under the conditions of variable heat source temperature

In this paper, a transient model of a silica gel–water adsorption chiller, which is developed in Shanghai Jiao Tong University (SJTU), is developed in order to simulate the evaporating, condensing, and adsorption temperature. Furthermore, this model is verified by a series of experiments. The theoretical studies and experimental data show that the coefficient of performance (COP) is influenced significantly by the variation rates of the heat source temperatures. The results also show that when this chiller is driven by solar energy, a buffer tank should be adopted in the system in order to get better performance when solar insolation is low, and should not be utilized when solar insolation is high, otherwise low COP will be gotten for the reason of the consumption of high electric energy.