DYNAMIC CRACK MODELS ON PROBLEM OF BRIDGING FIBER PULL-OUT OF COMPOSITE MATERIALS

An elastic analysis of an internal central crack with bridging fibers parallel to the free surface in an infinite orthotropic anisotropic elastic plane was performed. A dynamic model of bridging fiber pull-out of composite materials was presented. Resultingly the fiber failure is governed by maximum tensile stress, the fiber breaks and hence the crack extension should occur in self-similar fashion. By the methods of complex functions, the problem studied can be transformed into the dynamic model to the Reimann-Hilbert mixed boundary value problem, and a straightforward and easy analytical solution is presented. Analytical study on the crack propagation subjected to a ladder load and an instantaneous pulse loading is obtained respectively for orthotropic anisotropic body. By utilizing the solution, the concrete solutions of this model are attained by ways of superposition.     

Asymmetrical dynamic fracture model of bridging fiber pull-out of unidirectional composite materials

An elastic analysis of an internal crack with bridging fibers parallel to the free surface in an infinite orthotropic anisotropic elastic plane is studied, and asymmetrical dynamic fracture model of bridging fiber pull-out of unidirectional composite materials is presented for analyzing the distributions of stress and displacement with the internal asymmetrical crack under the loading conditions of an applied non-homogenous stress and the traction forces on crack faces yielded by the bridging fiber pull-out model. Thus the fiber failure is ascertained by maximum tensile stress, the fiber ruptures and hence the crack propagation should also appear in the modality of self-similarity. The formulation involves the development of a Riemann-Hilbert problem. Analytical solution of an asymmetrical propagation crack of unidirectional composite materials under the conditions of two increasing loads given is obtained, respectively. In terms of correlative material properties, the variable rule of dynamic stress intensity factor was depicted very well. After those analytical solutions were utilized by superposition theorem, the solutions of arbitrary complex problems could be gained.     

An asymmetrical dynamic model for bridging fiber pull-out of unidirectional composite materials

An elastic analysis of an internal crack with bridging fibers parallel to the free surface in an infinite orthotropic elastic plane is studied. An asymmetrical dynamic model for bridging fiber pull-out of unidirectional composite materials is presented for analyzing the distributions of stress and displacement with the internal asymmetrical crack under the loading conditions of an applied non-homogenous stress and the traction forces on crack faces yielded by the bridging fiber pull-out model. Thus the fiber failure is determined by maximum tensile stress, resulting in fiber rupture and hence the crack propagation would occur in a self-similarity manner. The formulation involves the development of a Riemann-Hilbert problem. Analytical solution of an asymmetrical propagation crack of unidirectional composite materials under the conditions of two moving loads given is obtained, respectively. After those analytical solutions were utilized by superposition theorem, the solutions of arbitrary complex problems could be obtained.     

复合材料桥纤维拔出的动态裂纹模型

在无限大正交各向异性体弹性平面上对复合材料桥纤维平行自由表面的内部中央裂纹提出了桥纤维拔出的动态裂纹模型。通过复变函数将其转化为Reimann-Hilbert混合边界值问题。求得了裂纹在坐标原点受载荷Px/t、Px2/t作用的解析解。利用这一解析解可通过迭加原理求得任意复杂问题的解。     

Dynamic crack tip fields and dynamic crack propagation characteristics of anisotropic material

Based on mechanics of anisotropic material, the dynamic crack propagation problem of I/II mixed mode crack in an infinite anisotropic body is investigated. Expressions of dynamic stress intensity factors for modes I and II crack are obtained. Components of dynamic stress and dynamic displacements around the crack tip are derived. The strain energy density theory is used to predict the dynamic crack extension angle. The critical strain energy density is determined by the strength parameters of anisotropic materials. The obtained dynamic crack tip fields are unified and applicable to the analysis of the crack tip fields of anisotropic material, orthotropic material and isotropic material under dynamic or static load. The obtained results show Crack propagation characteristics are represented by the mechanical properties of anisotropic material, i.e., crack propagation velocity M and fiber direction α. In particular, the fiber direction α and the crack propagation velocity M give greater influence on the variations of the stress fields and displacement fields. Fracture angle is found to depend not only on the crack propagation but also on the anisotropic character of the material.     

Debonding and fiber pull-out in reinforced composites

In order to evaluate the strength of fiber-reinforced composites, there is first the need to investigate the interfacial debonding and the pull-out of fibers in a fractured composite with intact fibers. This type of problem in crack bridging has been investigated by several authors based on different models and assumptions [1–7]. In this study, we will consider a three-dimensional model of a single fiber of finite length bonded by a finite cylindrical matrix with an initial crack existing in a portion of the interface. In the model, one end of the cylinder is so constrained that the axial component of displacement vanishes. A tensile stress is applied to the fiber at the other end. The aim is to determine the pull-out of the fiber and the critical condition for interfacial debonding. Both the fiber and the matrix are treated as elastic materials. Analysis is made based on a method using Papkovich-Neuber displacement potential functions for the problem of an elastic solid subjected to axisymmetrical boundary conditions. Solutions are found by means of the technique of trigonometrical series. Effects of initial misfit strains and frictional sliding between the fiber and the matrix over the interfacial crack are also included in the study.     

Residual strength of a fiber reinforced metal matrix composite with a crack

The residual strength of a cracked unidirectional fiver reinforced metal matrix composite is studied. We propose a bridging model based on the Dugdale strip yielding zones in the matrix ahead of the crack tips that accounts for ductile deformations of the matrix and fiber debonding and pull-out in the strip yielding zone. The bridging model is used to study the fracture of an anisotropic material and its residual strength is calculated numerically. The predicted results for a SiC/titanium composite agree well with the existing experimental data. It is found that a higher fiber bridging stress and a larger fiber pull-out length significantly contribute to the composite's residual strength. The composite's strength may be more notch-insensitive than the corresponding matrix material's strength depending on several factors such as fiber-matrix interface properties and the ratio of the matrix modulus to an ‘effective modulus’ of the composite.     

复合材料桥连的断裂动力学模型

复合材料产生裂纹后，其纤维处形成“桥连”，这是一个不可避免的现象。由于桥连问题很复杂．在数学方法的处理上有很大困难，至今人们研究大多是桥连的静力学问题．而对其动力学问题研究得很少。为了便于分析复合材料的问题，将桥连处用载荷代替，当裂纹高速扩展时．其纤维也连续地断裂。只有建立复合材料的桥连动力学模型，才能更好地研究复合材料的断裂动力学问题。通过复变函数论的方法，将所讨论的问题转化为Riemann—Hilbert问题。利用建立的动态模型和自相似方法，得到了正交异性体中扩展裂纹受运动的集中力P及阶跃载荷作用下位移、应力和动态应力强度因子的解析解，并通过叠加原理，最终求得了该模型的解。     

复合材料的裂纹动力学问题的模型

复合材料产生裂纹后，其纤维处形成“桥连”，这是一个不可避免的现象。由于桥连问题很复杂，在数学方法的处理上有很大困难，至今人们研究的大多是桥连的静力学问题，而对其动力学问题研究得很少。只有建立复合材料的桥连动力学模型，才能更好地研究复合材料的断裂动力学问题。为了便于分析复合材料的问题，将桥连处用载荷代替，当裂纹高速扩展时，其纤维也连续地断裂。通过复变函数论的方法，将所讨论的问题转化为Riemann-Hilbet问题。利用建立的动态模型和自相似方法，得到了正交异性体中扩展裂纹受运动的集中力Px/t及均布载荷作用下位移、应力和动态应力强度因子的解析解，并通过迭加原理，最终求得了该模型的解。     

Effects of curvilinear anisotropy on radially symmetric stresses in anisotropic linearly elastic solids

It has been known for some time that certain radial anisotropies in some linear elasticity problems can give rise to stress singularities which are absent in the corresponding isotropic problems. Recently related issues were examined by other authors in the context of plane strain axisymmetric deformations of a hollow circular cylindrically anisotropic linearly elastic cylinder under uniform external pressure, an anisotropic analog of the classic isotropic Lamé problem. In the isotropic case, as the external radius increases, the stresses rapidly approach those for a traction-free cavity in an infinite medium under remotely applied uniform compression. However, it has been shown that this does not occur when the cylinder is even slightly anisotropic. In this paper, we provide further elaboration on these issues. For the externally pressurized hollow cylinder (or disk), it is shown that for radially orthotropic materials, the maximum hoop stress occurs always on the inner boundary (as in the isotropic case) but that the stress concentration factor is infinite. For circumferentially orthotropic materials, if the tube is sufficiently thin, the maximum hoop stress always occurs on the inner boundary whereas for sufficiently thick tubes, the maximum hoop stress occurs at the outer boundary. For the case of an internally pressurized tube, the anisotropic problem does not give rise to such radical differences in stress behavior from the isotropic problem. Such differences do, however, arise in the problem of an anisotropic disk, in plane stress, rotating at a constant angular velocity about its center, as well as in the three-dimensional problem governing radially symmetric deformations of anisotropic externally pressurized hollow spheres. The anisotropies of concern here do arise in technological applications such as the processing of fiber composites as well as the casting of metals.     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Dynamic behavior of rectangular crack in three-dimensional orthotropic elastic medium by means of non-local theory

The dynamic behavior of a rectangular crack in a three-dimensional(3D)orthotropic elastic medium is investigated under a harmonic stress wave based on the non-local theory. The two-dimensional(2D) Fourier transform is applied, and the mixedboundary value problems are converted into three pairs of dual integral equations with the unknown variables being the displacement jumps across the crack surfaces. The effects of the geometric shape of the rectangular crack, the circular frequency of the incident waves,and the lattice parameter of the orthotropic elastic medium on the dynamic stress field near the crack edges are analyzed. The present solution exhibits no stress singularity at the rectangular crack edges, and the dynamic stress field near the rectangular crack edges is finite.     

Helical Fiber Pull-out in Biological Materials

Many biological materials, such as wood and bone, possess helicoid microstructures at microscale, which can serve as reinforcing elements to transfer stress between crack surfaces and improve the fracture toughness of their composites. Failure processes, such as fiber/matrix interface debonding and sliding associated with pull-out of helical fibers, are responsible mainly for the high energy dissipation needed for the fracture toughness enhancement. Here we present systemic analyses of the pull-out behavior of a helical fiber from an elastic matrix via the finite element method(FEM) simulation, with implications regarding the underlying toughening mechanism of helicoid microstructures. We find that, through their uniform curvature and torsion, helical fibers can provide high pull-out force and large interface areas, resulting in high energy dissipation that accounts, to a large extent, for the high toughness of biological materials. The helicity of fiber shape in terms of the helical angle has significant effects on the force-displacement relationships as well as the corresponding energy dissipation during fiber pull-out.     

An elastic–plastic crack bridging model for brittle-matrix fibrous composite beams under cyclic loading

A fibrous composite beam with an edge crack is submitted to a cyclic bending moment and the crack bridging actions due to the fibers. Assuming a general elastic-linearly hardening crack bridging model for the fibers and a linear-elastic law for the matrix, the statically indeterminate bridging actions are obtained from compatibility conditions. The elastic and plastic shake-down phenomena are examined in terms of generalised cross-sectional quantities and, by employing a fatigue crack growth law, the mechanical behaviour up to failure is captured. Within the framework of the proposed fracture mechanics-based model, the cyclic crack bridging due to debonding at fiber–matrix interface of short fibers is analysed in depth. By means of some simplifying assumptions, such a phenomenon can be described by a linear isotropic tensile softening/compressive hardening law. Finally, numerical examples are presented for fibrous composite beams with randomly distributed short fibers.     

Fracture dynamics problems of mode I semi-infinite crack for anisotropic orthotropic body

By means of the theory of complex functions, fracture dynamics problems of mode I semi- infinite crack for anisotropic orthotropic body were researched. Analytical solutions of stress, displacement, and dynamic stress intensity factor under the action of moving increasing loads Px ³/t ³, Pt ⁴/x ³, respectively, are very easily obtained utilizing the approaches of self-similar functions. In the light of relevant material’s coefficients, the alterable rule of dynamic stress intensity factor was depicted very well. The correlative closed solutions are attained based on the Riemann–Hilbert problems. After those analytical solutions were applied by the superposition principle, the solutions of discretional complex problems could be attained.     

An anisotropic elastoplastic constitutive formulation generalised for orthotropic materials

This paper presents a finite strain constitutive model to predict a complex elastoplastic deformation behaviour that involves very high pressures and shockwaves in orthotropic materials using an anisotropic Hill’s yield criterion by means of the evolving structural tensors. The yield surface of this hyperelastic–plastic constitutive model is aligned uniquely within the principal stress space due to the combination of Mandel stress tensor and a new generalised orthotropic pressure. The formulation is developed in the isoclinic configuration and allows for a unique treatment for elastic and plastic orthotropy. An isotropic hardening is adopted to define the evolution of plastic orthotropy. The important feature of the proposed hyperelastic–plastic constitutive model is the introduction of anisotropic effect in the Mie–Gruneisen equation of state (EOS). The formulation is further combined with Grady spall failure model to predict spall failure in the materials. The proposed constitutive model is implemented as a new material model in the Lawrence Livermore National Laboratory (LLNL)-DYNA3D code of UTHM’s version, named Material Type 92 (Mat92). The combination of the proposed stress tensor decomposition and the Mie–Gruneisen EOS requires some modifications in the code to reflect the formulation of the generalised orthotropic pressure. The validation approach is also presented in this paper for guidance purpose. The \({\varvec{\psi }}\) tensor used to define the alignment of the adopted yield surface is first validated. This is continued with an internal validation related to elastic isotropic, elastic orthotropic and elastic–plastic orthotropic of the proposed formulation before a comparison against range of plate impact test data at 234, 450 and \({\mathrm {895\,ms}}^{\mathrm {-1}}\) impact velocities is performed. A good agreement is obtained in each test.     

Concentrated impact loading on one face of a semi-infinite crack in an anisotropic elastic material

The response of an unbounded anisotropic elastic body containing a semi-infinite crack subjected to a concentrated impact force on one of the crack faces is studied. An exact solution of the dynamic stress intensity factors is obtained from a linear superposition of the solution of Lamb’s problem and a solution of a dislocation emitting from the crack tip. The stress intensity factors exhibit square-root singularity upon the arrival of the Rayleigh wave at the crack tip. As the Rayleigh wave passes through the crack tip, the stress intensity factors either instantaneously assume the static values or gradually approach to zero. Several numerical examples are given for isotropic, cubic and orthotropic materials.     

Weight functions for interface cracks in dissimilar anisotropic materials

Bueckner‘s work conjugate integral customarily adopted for linear elastic materials is established for an interface crack in dissimilar anisotropic materials. The difficulties in separating Stroh‘s six complex arguments involved in the integral for the dissimilar materials are overcome and then the explicit function representations of the integral are given and studied in detail. It is found that the pseudo-orthogonal properties of the eigenfunction expansion form (EEF) for a crack presented previously in isotropic elastic cases, in isotopic bimaterial cases, and in orthotropic cases are also valid in the present dissimilar arbitrary anisotropic cases. The relation between Bueckner‘s work conjugate integral and the J-integral in these cases is obtained by introducing a complementary stressdisplacement state. Finally, some useful path-independent integrals and weight functions are proposed for calculating the crack tip parameters such as the stress intensity factors.     

含孔正交各向异性复合材料板的动应力集中

运用一种改进的非结构化四边形格子法,对含孔正交各向异性板条受面内冲击拉伸时弹性应力波的传播过程和孔边的动应力集中进行了研究．非结构化格子法采用与有限元类似的网格剖分方法,并基于围绕每个节点的积分平衡方程,并自然满足复杂边界的自由边界条件．计算中不需存储刚度矩阵,因而计算速度快、效率高、节省内存,在解决应力波传播问题中具有显著的优越性．通过对多种工况进行数值模拟,分析了材料的各向异性性质、纤维方向、孔径比、加载脉冲周期等参数对孔边动应力的影响,得到了一些规律性的结果．并与现有实验结果进行对比,验证了该方法的有效性．