Dynamic fracture of concrete – compact tension specimen

The behavior of concrete structures is strongly influenced by the loading rate. Compared to quasi-static loading concrete loaded by impact loading acts in a different way. First, there is a strain-rate influence on strength, stiffness, and ductility, and, second, there are inertia forces activated. Both influences are clearly demonstrated in experiments. Moreover, for concrete structures, which exhibit damage and fracture phenomena, the failure mode and cracking pattern depend on loading rate. In general, there is a tendency that with the increase of loading rate the failure mode changes from mode-I to mixed mode. Furthermore, theoretical and experimental investigations indicate that after the crack reaches critical speed of propagation there is crack branching. The present paper focuses on 3D finite-element study of the crack propagation of the concrete compact tension specimen. The rate sensitive microplane model is used as a constitutive law for concrete. The strain-rate influence is captured by the activation energy theory. Inertia forces are implicitly accounted for through dynamic finite element analysis. The results of the study show that the fracture of the specimen strongly depends on the loading rate. For relatively low loading rates there is a single crack due to the mode-I fracture. However, with the increase of loading rate crack branching is observed. Up to certain threshold (critical) loading rate the maximal crack velocity increases with increase of loading rate, however, for higher loading rates maximal velocity of the crack propagation becomes independent of the loading rate. The critical crack velocity at the onset of crack branching is found to be approximately 500 m/s.     

Ductile–brittle transitions in the fracture of plastically-deforming,adhesively-bonded structures. Part I: Experimental studies

Rate effects for adhesively-bonded joints in steel sheets failing by mode-I fracture and plastic deformation were examined. Three types of test geometries were used to provide a range of crack velocities between 0.1 and 5000 mm/s: a DCB geometry under displacement control, a wedge geometry under displacement control, and a wedge geometry loaded under impact conditions. Two fracture modes were observed: quasi-static crack growth and dynamic crack growth. The quasi-static crack growth was associated with a toughened mode of failure; the dynamic crack growth was associated with a more brittle mode of failure. The experiments indicated that the fracture parameters for the quasi-static crack growth were rate independent, and that quasi-static crack growth could occur even at the highest crack velocities. Effects of rate appeared to be limited to the ease with which a transition to dynamic fracture could be triggered. This transition appeared to be stochastic in nature, it did not appear to be associated with the attainment of any critical value for crack velocity or loading rate. While the mode-I quasi-static fracture behavior appeared to be rate independent, an increase in the tendency for dynamic fracture to be triggered as the crack velocity increased did have the effect of decreasing the average energy dissipated during fracture at higher loading rates.     

An improved atomistic simulation based mixed-mode cohesive zone law considering non-planar crack growth

A novel and improved atomistic simulation based cohesive zone law characterizing interfacial debonding is developed which explicitly accounts for the non-planarity of the crack propagation. Group of atoms in the simulation constituting cohesive zones which are used to obtain local stress and crack opening displacement data are determined dynamically during the non-planar crack growth as they cannot be determined apriori. The methodology is used to study the debonding of Σ5 (2 1 0)/[0 0 1] symmetric tilt grain boundary interface in a Cu bicrystal under several mixed mode loading conditions. Simulations show that such bicrystalline specimen exhibits three types of energy dissipative mechanisms – shear coupled GB migration (SCM) away from the crack-tips, change in spacial orientation of GB structural units rendering highly disordered grain boundary near the crack tips and brittle intergranular fracture. Which combination of these three deformation mechanism will be active influencing the degree of non-planarity of the crack propagation at various stages of loading depends on the loading mode-mixity. As the ratio of shear component of the loading parallel to the GB plane and normal to the tilt axis with respect to the normal loading increases (thereby increasing the mode-mixity), overall strain-to-failure also increases and SCM tends to become the dominant deformation mechanism. Through this framework, analytical functional forms and parameters describing cohesive laws for both normal and shear traction as a function of the mode-mixity of the loading and crack opening displacement are predicted.     

Fatigue properties of age-hardened Al alloy 2017-T4 under ultrasonic loading

Fatigue properties of age-hardened Al alloy 2017-T4 under ultrasonic loading frequency (20 kHz) were investigated and compared with the results under conventional loading of rotating bending (50 Hz). The growth of a crack retarded at about 500 µm in surface length under ultrasonic loading, while at about 20 µm under rotating bending. Although striations being a typical fracture mechanism were observed under conventional loading, most of fracture surface was covered with many facets under ultrasonic loading. These facets were also observed under rotating bending in nitrogen gas. The difference in growth mechanism depending on the loading frequency and the retardation of a crack growth under ultrasonic loading may be caused by the environment at the crack tip due to high crack growth rate under ultrasonic loading.     

Lattice orientation effects on void growth and coalescence in fcc single crystals

Void growth and coalescence in fcc single crystals were studied using crystal plasticity under uniaxial and biaxial loading conditions and various orientations of the crystalline lattice. A 2D plane strain unit cell with one and two cylindrical voids was employed using three-dimensional 12 potentially active slip systems. The results were compared to five representative orientations of the tensile axis on the stereographic triangle. For uniaxial tension conditions, the void volume fraction increase under the applied load is strongly dependent on the crystallographic orientation with respect to the tensile axis. For some orientations of the tensile axis, such as [1 0 0] or [1 1 0], the voids exhibited a growth rate twice as fast compared with other orientations ([1 0 0], [2 1 1]). Void growth and coalescence simulations under uniaxial loading indicated that during deformation along some orientations with asymmetry of the slip systems, the voids experienced rotation and shape distortion, due mainly to lattice reorientation. Coalescence effects are shown to diminish the influence of lattice orientation on the void volume fraction increase, but noteworthy differences are still present. Under biaxial loading conditions, practically all differences in the void volume fraction for different orientations of the tensile axes during void growth vanish. These results lead to the conclusion that at microstructural length scales in regions under intense biaxiality/triaxiality conditions, such as crack tip or notched regions, the plastic anisotropy due to the initial lattice orientation has only a minor role in influencing the void growth rate. In such situations, void growth and coalescence are mainly determined by the stress triaxiality, the magnitude of accumulated strain, and the spatial localization of such plastic strains.     

CTS试件中复合型疲劳裂纹扩展

针对复合型循环载荷作用下的金属构件中的裂纹扩展问题进行了实验分析和理论建模. 首先 采用紧凑拉剪试件(CTS)和 Richard研制的复合型载荷加载装置，对承受复合型循环载荷的裂纹进行了实验研究. 实验选择了两种金属材料试件，分别承受3种形式的复合型循环载荷的作用，在裂纹尖端具 有相同的初始应力场强度的条件下考察复合型循环载荷对裂纹扩展规律的影响. 实验结果表明，疲劳裂纹的扩展速率与加载角度有关. 对于同样金属材料的试件，当裂尖处 初始应力场强度相等时，载荷越接近于II型，裂纹增长速率越快. 采用等效应力强度 因子(I型和II型应力强度因子的组合)、裂纹扩展速率及复合强度等参数，以实验数据为 基础，建立了一个疲劳裂纹扩展模型，用来预测裂纹在不同模式疲劳载荷作用下的扩展速率. 为验证其有效性，该模型被应用于钢制试件的数值模拟计算中. 实验结果与模拟计算曲线保 持一致，表明该模型可以用来估算带裂纹金属构件的寿命.     

Periodic group crack problems in an infinite plate

This paper investigates periodic group crack problems in an infinite plate. The periodic group crack is composed of infinite groups with numbering from j = −∞, …, −2, −1, 0, 1, 2, …, to j = ∞, and the groups are placed periodically. The same loading condition and the same geometry are assumed for cracks in all groups. A singular integral equation is used to solve the problems. The singular integral equation is formulated on cracks of the 0th group (or the central group) with the collection of influences from the infinite groups. The influences of many neighboring groups to the central group are evaluated exactly. Meantime, the influences of many remote groups to the central group can be summed up into one term approximately. The stress intensity factors at crack tips can be evaluated from the solution of the singular integral equation. It is found from some sample problems that the obtained results are very accurate. Finally, several numerical examples are presented and interaction among the group cracks is addressed.     

Thermo-mechanical response of nylon 101 under uniaxial and multi-axial loadings: Part I,Experimental results over wide ranges of temperatures and strain rates

A comprehensive study of the thermo-mechanical response of a thermoplastic polymer, nylon 101 is presented. Quasi-static and dynamic compression uniaxial and multi-axial experiments (stress states) were performed at a wide range of strain rates (10⁻⁵ to 5000 s⁻¹) and temperatures (−60 to 177 °C or −76 to 350 °F). The material is found to be non-linearly dependent on strain rate and temperature. The change in volume after plastic deformation is investigated and is found to be negligibly small. The relaxation and creep responses at room temperature are found to be dependent on strain rate and the stress–strain level at which these phenomena are initiated. Total deformation is decomposed into visco-elastic and visco-plastic components; these components have been determined at different levels of deformation. Results from non-proportional uniaxial to biaxial compression, and torsion experiments, are also reported for three different strain rates at room temperature. It is shown that nylon 101 has a response dependent on the hydrostatic pressure.     

Influence of temperature on a carbon–fibre epoxy composite subjected to static and fatigue loading under mode-I delamination

In this paper, interlaminar crack initiation and propagation under mode-I with static and fatigue loading of a composite material are experimentally assessed for different test temperatures. The material under study is made of a 3501-6 epoxy matrix reinforced with AS4 unidirectional carbon fibres, with a symmetric laminate configuration [0°]_16/S. In the experimental programme, DCB specimens were tested under static and fatigue loading. Based on the results obtained from static tests, fatigue tests were programmed to analyse the mode-I fatigue behaviour, so the necessary number of cycles was calculated for initiation and propagation of the crack at the different temperatures. G–N curves were determined under fatigue loading, N being the number of cycles at which delamination begins for a given energy release rate. G_ICmax–a, a–N and da/dN–a curves were also determined for different G_cr rates (90%, 85%, 75%, etc.) and different test temperatures: 90 °C, 50 °C, 20 °C, 0 °C, ?30 °C and ?60 °C.     

What physical phenomena can be neglected when modelling concrete at high temperature? A comparative study. Part 1: Physical phenomena and mathematical model

The paper deals with modelling of hygro-thermal performance and thermo-chemical degradation of concrete exposed to high temperature. Several possible simplifications in modelling of heat and mass transport phenomena in heated concrete are considered and their effect on the results of numerical simulations is analyzed.A mathematical model of concrete at high temperature, already extensively validated with respect to experiments, is used as the reference model. It is based on mechanics of multiphase porous media and considers all important couplings and material nonlinearities, as well as different properties of water above the critical point of water, i.e. 647.3 K (374.15 °C).In this part of the paper, first physical phenomena, as well as heat and mass flux and sources in a concrete element are studied, both during slow and fast heating process, to examine the relative importance of different flux components. Then, the mathematical model of concrete at high temperature, developed by Authors in the last 10 years, is briefly presented and for the first time all the constitutive relationships of the model are summarized and discussed in detail. Finally, the method of numerical solution of the model equations is thoroughly presented.In the companion paper (part II) a brief literature review of the existing mathematical models of concrete at high temperature and a summary of their main features and physical assumptions will be presented. Then, extensive numerical studies will be performed with several simplified models, neglecting chosen physical phenomenon or flux components, to evaluate the difference between the results obtained with the simplified models and with the reference model. The study will concern hygric, thermal and degradation performance of 1-D and 2-D axisymmetric concrete elements during fast and slow heating. The analysis will allow us to indicate which simplifications in modeling of concrete at high temperature are practically and physically possible, without generating excessive errors with respect to the full reference model.     

Thermo-mechanical behaviour of TRIP 1000 steel sheets subjected to low velocity perforation by conical projectiles at different temperatures

This paper presents and analyzes the behaviour of TRIP 1000 steel sheets subjected to low velocity perforation by conical projectiles. The relevance of this material resides in the potential transformation of retained austenite to martensite during impact loading. This process leads to an increase in strength and ductility of the material. However, this transformation takes place only under certain loading conditions strongly dependent on the initial temperature and deformation rate. In order to study the material behaviour under impact loading, perforation tests have been performed using a drop weight tower. Experiments were carried out at two different initial temperatures T₀ = 213 K and T₀ = 288 K, and within the range of impact velocities 2.5 m/s ? V₀ ? 4.5 m/s. The experimental setup enabled the measuring of impact velocity, residual velocity, load-time history and failure mode. In addition, dry and lubricated contacts between the striker and the plate have been investigated. Finally, by using X-ray diffraction it has been shown that no martensitic transformation takes place during the perforation process. The causes involving the none-appearance of martensite are examined.     

Novel multistage solid–liquid circulating fluidized bed: Hydrodynamic characteristics

The present work proposes a novel radially cross-flow multistage solid–liquid circulating fluidized bed (SLCFB). The SLCFB primarily consists of a single multistage column (having an inner diameter of 100 mm and length of 1.40 m), which is divided into two sections wherein both the steps of utilization or loading (e.g., adsorption and catalytic reaction) and regeneration of the solid phase can be carried out simultaneously in continuous mode. The hydrodynamic characteristics were studied using ion exchange resin as the solid phase and water as the fluidizing medium. The loading and flooding states were determined for three particle sizes; i.e., 0.30, 0.42, and 0.61 mm. The effects of the superficial liquid velocity and solid feed rate on the solid hold-up were investigated under loading and flooding conditions. The solid hold-up increases with an increase in the solid feed rate and decreases with an increase in the superficial liquid velocity. An artificial-intelligence formalism, namely the multilayer perceptron neural network (MLPNN), was employed for the prediction of the solid hold-up. The input space of MLPNN-based model consists of four parameters, representing operating and system parameters of the proposed SLCFB. The developed MLPNN-based model has excellent prediction accuracy and generalization capability.     

Thermo-elastic crack branching in general anisotropic media

Thermal fields may exist in addition to mechanical loading, for example, due to short term exposure to fire. In this paper, the branching of cracks in the presence of combined thermal and mechanical loads is investigated for general anisotropic media by employing the theory of Stroh’s dislocation formalism, extended to thermo-elasticity in matrix notation. A general solution to the thermo-elastic crack problem for an anisotropic material under arbitrary loading is obtained in a compact form. Green’s functions are also presented for a thermal dislocation (heat vortex) and a conventional dislocation (or, referred as mechanical dislocation), which are formulated considering the cuts located at an arbitrary angle with respect to the x₁ axis of the coordinate system (x₁, x₂, x₃). Using the derived compact expressions, the interaction between the crack and the dislocation is studied and a closed form solution for this interaction is obtained. The branching portion of the thermo-elastic crack is modelled as a continuous distribution of dislocations. This problem is then converted into a set of singular integral equations. Numerical results are presented to illustrate the possible effects of thermal loading on the propagation of the branched crack.     

Coplanar propagation paths of 3D cracks in infinite bodies loaded in shear

Bower and Ortiz, recently followed by Lazarus, developed a powerful method, based on a theoretical work of Rice, for numerical simulation of planar propagation paths of mode 1 cracks in infinite isotropic elastic bodies. The efficiency of this method arose from the need for the sole 1D meshing of the crack front. This paper presents an extension of Rice’s theoretical work and the associated numerical scheme to mixed-mode (2 + 3) shear loadings. Propagation is supposed to be channeled along some weak planar layer and to remain therefore coplanar, as in the case of a geological fault for instance. The capabilities of the method are illustrated by computing the propagation paths of cracks with various initial contours (circular, elliptic, rectangular, heart-shaped) in both fatigue and brittle fracture. The crack quickly reaches a stable, almost elliptic shape in all cases. An approximate but accurate analytic formula for the ratio of the axes of this stable shape is derived.     

Experimental investigation on strength and failure behavior of pre-cracked marble under conventional triaxial compression

Fractures in natural rocks have an important effect on the strength and failure behavior of rock mass, which are often evaluated in rock engineering practice. The theoretical evaluation of mechanical behavior of fractured rock mass has no satisfactory answer due to the role of confining pressure and crack geometry. Therefore, in this paper, conventional triaxial compression experiments were carried out to study the strength and failure behavior of marble samples with two pre-existing closed cracks in non-overlapping geometry. Based on the experimental results of a number of triaxial compression tests, the effect of crack coalescence on the axial supporting capacity and deformation property were investigated with different confining pressures. The results show that intact samples and flawed samples (marble with pre-existing cracks) have different deformation properties after peak stress, which change from brittleness to plasticity and ductility with the increase of confining pressure. The peak strength and failure mode are found depending not only on the geometry of flaw, but also on the confining pressure. The strength of flawed samples shows distinct non-linear behavior, which is in a better agreement with non-linear Hoek–Brown criterion than linear Mohr–Coulomb criterion. For a kind of rock that has been evaluated as a Hoek–Brown material, a new evaluation criterion is put forward by adopting optimal approximation polynomial theory, which can be used to confirm more precisely the strength parameters (cohesion and internal friction angle) of flawed samples. For intact samples, the marble leads to typical shear failure mode with a single fracture surface under different confining pressures, while for flawed samples, under uniaxial compression and a lower confining pressure (σ₃ = 10 MPa), tests for coarse and medium marble (the coarse and medium refer to the grain size) exhibit three basic failure modes, i.e., tensile mode, shear mode, and mixed mode (tensile and shear). Shear mode is associated with lower strength behavior. However, under higher confining pressures (σ₃ = 30 MPa), for coarse marble, the axial supporting capacity is not related to the geometry of flaw. The friction among crystal grains determines the strength behavior of coarse marble. For medium marble, the failure mode and deformation behavior are dependent on the crack coalescence in the sample. The present research provides increased understanding of the fundamental nature of rock failure under conventional triaxial compression.     

Relating inhomogeneous deformation to local texture in zirconium through grain-scale digital image correlation strain mapping experiments

The effect of local texture on inhomogeneous plastic deformation is studied in zirconium subjected to uniaxial compression. Cross-rolled commercially pure Zr 702 plate that had a strong basal (0 0 0 1) texture through the plate thickness, and a non-basal texture in cross-section, was obtained. At a compressive strain rate of 1 s^?1, samples loaded either in the through-thickness or in-plane directions exhibited significant differences in yield strength, hardening response and failure mechanisms. These macroscopic differences are related to microstructural features by combining information from electron backscattered diffraction with real time in situ imaging and subsequent full-field strain measurements obtained using digital image correlation. Experimental results indicate that the through-thickness loaded zirconium samples, which show a strong basal-texture in the loading direction, do not deform homogeneously – implying the lack of a representative volume element. The detailed surface deformation fields provided by digital image correlation allow for a qualitative and quantitative analysis of the relationship between grain orientation and patterns of deformation bands that form as the precursors to development of an adiabatic shear band in the through-thickness loaded sample. For the in-plane loaded samples, inhomogeneities still exist at the microscale, but the collective behavior of several grains leads to a homogeneous response at the macroscale. It is observed that local texture for hcp polycrystals, which are significantly slip restricted, can directly affect both local and global response, even at low to moderate plastic strains.     

The fracture behaviour of adhesively-bonded composite joints: Effects of rate of test and mode of loading

The present paper discusses the results of an investigation into the effects of test rate and the mode of loading on the fracture energy, G_c, of adhesively-bonded fibre-composite joints. Various carbon-fibre reinforced-polymer (CFRP) matrix composite substrates have been bonded using two different types of automotive structural epoxy-adhesives. They have been tested via loading the bonded joints in mode I (tensile), mode II (in-plane shear) and mixed-mode I/II from slow rates (i.e., of about 10^?5 m/s) up to relatively high rates of test of about 15 m/s. The high-rate tests were photographed using a high-speed digital video camera to record the deformation of the joint and the fracture behaviour. An analysis strategy has been developed for the various modes of loading (i) to account for the observed fracture behaviour, (ii) to circumvent the problems posed by oscillations in the load traces due to the presence of dynamic effects in the faster tests, and (iii) to account for the kinetic energy associated with the moving specimen arms in the faster tests. Based on the analysis strategy developed, the effect of the test rate on the fracture energy, G_c, for the different loading modes for the joints has been ascertained. Furthermore, various different fracture paths were observed in the tests. They were either cohesive, in the adhesive layer, or interlaminar in the composite substrates. The exact fracture path observed was a function of (i) the type of composite substrate, (ii) the type of adhesive, and (iii) the mode of loading employed. However, the nature of the fracture path was found to be quite insensitive to the test rate. Essentially, it was found that joints subjected to mixed-mode I/II loading were more likely to exhibit an interlaminar fracture path in the composite substrates than when loaded in either pure modes I or II. The propensity for a given joint to exhibit such a fracture path via delamination of the composite substrate has been explained by calculating the transverse tensile stresses induced in the loaded composite arms, and comparing this value to the measured transverse tensile strength of the composite. Following this approach, the underlying reasons for the observed fracture path were identified and could be predicted. Also, the proposed scheme provides a route to design against delamination failure occurring in adhesively-bonded fibre-composite test specimens.     

A mode III crack in a piezoelectric semiconductor of crystals with 6 mm symmetry

A finite mode III crack in a piezoelectric semiconductor of 6 mm crystals is analyzed. Fourier transform is employed to reduce the mixed boundary value problem to a pair of dual-integral equations. Numerical solution of these equations yields coupled electromechanical fields, the intensity factor and the energy release rate near the crack tip. Numerical results are presented graphically to show the fracture behavior which is affected by the semiconduction.     

Numerical study on coalescence of two pre-existing coplanar flaws in rock

Crack propagation and coalescence processes are the fundamental mechanisms leading to progressive failure processes in rock masses, in which parallel non-persistent rock joints are commonly involved. The coalescence behavior of the latter, which are represented as pre-existing coplanar flaws (cracks), is numerically investigated in the present study. By using AUTODYN as the numerical tool, the present study systematically simulates the coalescence of two pre-existing coplanar flaws in rock under compression. The cumulative damage failure criterion is adopted in the numerical models to simulate the cumulative damage process in the crack initiation and propagation. The crack types (shear or tensile) are identified by analyzing the mechanics information associated with the crack initiation and propagation processes. The simulation results, which are generally in a good accordance with physical experimental results, indicate that the ligament length and the flaw inclination angle have a great influence on the coalescence pattern. The coalescence pattern is relatively simple for the flaw arrangements with a short ligament length, which becomes more complicated for those with a long ligament length. The coalescence trajectory is composed of shear cracks only when the flaw inclination angle is small (such as β ? 30°). When the pre-existing flaws are steep (such as β ? 75°), the coalescence trajectory is composed of tensile cracks as well as shear cracks. When the inclination angle is close to the failure angle of the corresponding intact rock material, and the ligament length is not long (such as L ? 2a), the direct shear coalescence is the more favorable coalescence pattern. In the special case that the two pre-existing flaws are vertical, the model will have a direct tensile coalescence pattern when the ligament length is short (L ? a), while the coalescence between the two inner flaw tips is not easy to achieve if the ligament length is long (L ? 2a).