Three-dimensional continuum analysis for a unidirectional composite with a broken fiber

The three-dimensional problem of a periodic unidirectional composite with a penny-shaped crack traversing one of the fibers is analyzed by the continuum equations of elasticity. The solution of the crack problem is represented by a superposition of weighted unit normal displacement jump solutions, everyone of which forms a Green’s function. The Green’s functions for the unbounded periodic composite are obtained by the combined use of the representative cell method and the higher-order theory. The representative cell method, based on the triple discrete Fourier transform, allows the reduction of the problem of an infinite domain to a problem of a finite one in the transform space. This problem is solved by the higher-order theory according to which the transformed displacement vector is expressed by a second order expansion in terms of local coordinates, in conjunction with the equilibrium equations and the relevant boundary conditions. The actual elastic field is obtained by a numerical evaluation of the inverse transform. The accuracy of the suggested approach is verified by a comparison with the exact analytical solution for a penny-shaped crack embedded in a homogeneous medium. Results for a unidirectional composite with a broken fiber are given for various fiber volume fractions and fiber-to-matrix stiffness ratios. It is shown that for certain parameter combinations the use of the average stress in the fiber, as it is employed in the framework of the shear lag approach, for the prediction of composite’s strength, leads to an over estimation. To this end, the concept of “point stress concentration factor” is introduced to characterize the strength of the composite with a broken fiber. Several generalizations of the proposed approach are offered.     

The effect of a fiber loss in periodic composites

Three types of analyses are combined to investigate the effect of missing fibers in periodic continuous fiber composites that are subjected to thermomechanical loadings. The representative cell method is employed in the first analysis for the construction of Green’s functions elastic fields for the fiber–matrix interfacial jumps problem. As a result, the infinite domain problem is reduced, in conjunction with the discrete Fourier transform, to a finite domain problem on which Born–von Karman type boundary conditions are applied. In the second analysis, the transformed elastic field is determined by a second-order expansion of the displacement vector in terms of local coordinates, and by imposing the equilibrium equations, the interfacial traction and displacement conditions, and the Born–von Karman type boundary conditions. The actual non-periodic elastic field at any point is obtained from the Fourier-transformed fields by a numerical inversion. In the third one, a micromechanical analysis for periodic continuous fiber composites in which all fibers are perfectly bonded to the surrounding matrix provides the elastic field within the phases. A superposition of the thermoelastic fields obtained from the first and third analysis provides the traction-free boundary conditions at the interface of the missing fibers. The accuracy of the offered approach is verified by comparison with analytical solutions that exist in some special cases. Results show the effect of a missing fiber in boron/epoxy and glass/epoxy composites that are subjected to various types of thermomechanical loadings.     

Heating and cooling effects generated by the dynamic stresses of cracked interfaces and layers in laminated composites

The temperature field generated by the sudden application of a far-field mechanical loading of a periodically layered composite with an interfacial crack or with a cracked layer parallel to the interfaces is determined. As a result of the crack’s existence, the periodicities of the structure and the thermoelastic field are lost. The complexity of the resulting problem is resolved by the combined application of the representative cell method and the full (two-way) dynamic thermomechanical equations. In the former analysis, due to the loss of periodicity the dynamic thermoelastic Green’s functions are generated, in conjunction with the double finite discrete Fourier transform. In the latter one, the transformed displacements and temperature are expressed by second-order expansions and the strong-form of the elastodynamic and energy equations together with the various interfacial and the so called Born–von Karman boundary conditions are imposed in the average sense (in the transform domain). The results exhibit the induced temperature field at any point in the plane of the crack. The generated temperature fields show the cooling and heating zones in both Mode I and Mode II deformations. In addition, the adiabatic assumption (according to which the heat conduction is a priori ignored) is assessed by comparing the computed temperature field with the corresponding one based on the full thermomechanical coupling.     

Constitutive equations for porous plane-strain gradient elasticity obtained by homogenization

This paper addresses the problem of plane-strain gradient elasticity models derived by higher-order homogenization. A microstructure that consists of cylindrical voids surrounded by a linear elastic matrix material is considered. Both plane-stress and plane-strain conditions are assumed and the homogenization is performed by means of a cylindrical representative volume element (RVE) subjected to quadratic boundary displacements. The constitutive equations for the equivalent medium at the macroscale are obtained analytically by means of the Airy’s stress function in conjunction with Fourier series. Furthermore, a failure criterion based on the maximum hoop stress on the void surface is formulated. A mixed finite-element formulation has been implemented into the commercial finite-element program Abaqus. Using the constitutive relations derived, numerical simulations were performed in order to compute the stress concentration at a hole with varying parameters of the constitutive equations. The results predicted by the model are discussed in comparison with the results of the theory of simple materials.     

Higher-order crack tip fields for functionally graded material plate with transverse shear deformation

The crack tip fields are investigated for a cracked functionally graded material (FGM) plate by Reissner’s linear plate theory with the consideration of the transverse shear deformation generated by bending. The elastic modulus and Poisson’s ratio of the functionally graded plates are assumed to vary continuously through the coordinate y, according to a linear law and a constant, respectively. The governing equations, i.e., the 6th-order partial differential equations with variable coefficients, are derived in the polar coordinate system based on Reissner’s plate theory. Furthermore, the generalized displacements are treated in a separation-of-variable form, and the higher-order crack tip fields of the cracked FGM plate are obtained by the eigen-expansion method. It is found that the analytic solutions degenerate to the corresponding fields of the isotropic homogeneous plate with Reissner’s effect when the in-homogeneity parameter approaches zero.     

Multiple cracking of fiber/matrix composites—Analysis of normal extension

This paper is concerned with the fracture of a fiber embedded in a matrix of finite radius. There is a periodic array of cracks in the fiber along the central axis of the medium. The paper accounts for the cases of axial extension and residual temperature change of the composite medium. Fourier and Hankel transforms are used to reduce the problem to the solution of a system of dual integral equations, which are solved by the singular integral equation technique. Rigorous fracture mechanics analysis, which exactly satisfies all boundary conditions of the problem, is conducted. Numerical solutions for the crack tip field and the stress in the fiber are obtained for various values such as crack radius, crack spacing and fiber volume fraction.     

Dynamic response of an elastic medium containing crack

A fundamental solution for an infinite elastic medium containing a penny-shaped crack subjected to dynamic torsional surface tractions is attempted. A double Laplace–Hankel integral transform with respect to time and space is applied both to motion equation and boundary conditions yielding dual integral equations. The solution of the derived dual integral equations is based on an analytic procedure using theorems of Bessel functions and ordinary differential equations. The dynamic displacements’ field is obtained by inversion of the corresponding Laplace–Hankel transformed variable. Results of a representative example for a crack subjected to pulse surface tractions are obtained and discussed.     

Damped vibration of a graphene sheet using a higher-order nonlocal strain-gradient Kirchhoff plate model

A higher-order nonlocal strain-gradient model is presented for the damped vibration analysis of single-layer graphene sheets (SLGSs) in hygrothermal environment. Based on Kirchhoff plate theory in conjunction with a higher-order (bi-Helmholtz) nonlocal strain gradient theory, the equations of motion are obtained using Hamilton's principle. The higher-order nonlocal strain gradient theory has lower- and higher-order nonlocal parameters and a material characteristic parameter. The presented model can reasonably interpret the softening effects of the SLGS, and indicates a reasonably good match with the experimental flexural frequencies. Finally, the roles of viscous and structural damping coefficients, small-scale parameters, hygrothermal environment and elastic foundation on the vibrational responses of SLGSs are studied in detail.     

考虑电场梯度的挠曲电纳米板弯曲性能分析

具有尺度依赖的挠曲电效应在器件的设计中扮演着越来越关键的角色, 研究人员在微纳米尺度多物理场分析中进行了大量工作. 基于考虑挠曲电和电场梯度效应的弹性介电材料非经典理论, 以二维纳米板为例, 通过理论建模, 分析纳米板在弯曲问题中的力?电耦合行为. 根据Mindlin假设给出板的位移场和电势场的一阶截断, 选取板的材料为立方晶体(m3m点群), 将广义三维本构方程代入到高阶应力、高阶偶应力、高阶电位移和高阶电四极矩的表达式中得到相应的二维本构方程, 利用弹性电介质变分原理得到板的控制方程和边界上的线积分等式, 分别将二维本构方程和边界上外法线的方向余弦代入, 得到板的高阶弯曲方程、高阶电势方程以及对应的四边简支边界条件. 利用四边简支矩形板的高阶弯曲方程、高阶电势方程和相应的边界条件, 根据Navier解理论, 求解纳米板的电势场, 重点分析电场梯度对板内一阶电势的影响. 数值计算结果表明: 电场梯度对纳米板中由挠曲电效应产生的一阶电势有削弱作用, 且材料参数g₁₁越大, 一阶电势受到的削弱越大; 同时电场梯度的存在消除了纳米板在受横向集中载荷作用时一阶电势的奇异性. 本文是对具有挠曲电效应和电场梯度效应的纳米板结构分析理论的一个扩展, 为微纳米尺度器件的结构设计提供参考.      

Thermoelastic stress intensity factors for a cracked layer between two different anisotropic solids

Steady-state anisotropic thermoelasticity equations are used to obtain the stress intensity factors for a cracked layer sandwiched between two different anisotropic elastic solids. The anisotropy is assumed to arise from discrete fibers whose orientation could alter with reference to the crack edges. A generalized plane deformation prevails in the dissimilar media domain with a line of discontinuity disturbing a uniform heat flow. The flexibility/stiffness matrix approach is used such that the crack problem reduces to solving two sets of singular integral equations. Numerical values of the crack tip stress-intensity factors are obtained for various crack size, crack location, crack surface insulation, fiber volume fraction and orientation angles. The results are displayed graphically.     

Elastic-plastic fields near crack tip growing step-by-step in a perfectly plastic solid

In this paper, based on the three-dimensional flow theory of plasticity, the fundametal equations for plane strain problem of elastic-perfectly plastic solids are presented. By using these equations the elastic-plastic fields near the crack tip growing step-by-step in an elastic incompressible-perfectly plastic solid are analysed.The first order asymptotic solutions for the stress field and velocity fields near the crack tip are obtained. The solutions show the evolution process of elastic unloading domain and the development process of central fan domain and reveal the possibility of the presence of the secondary plastic domain. The second order asymptotic solution for stress field is also presented.     

Thermo‐hydrodynamic‐mechanical multiphase flow model for CO2 sequestration in fracturing porous media

In this paper, we introduce a fully coupled thermo‐hydrodynamic‐mechanical computational model for multiphase flow in a deformable porous solid, exhibiting crack propagation due to fluid dynamics, with focus on CO₂ geosequestration. The geometry is described by a matrix domain, a fracture domain, and a matrix‐fracture domain. The fluid flow in the matrix domain is governed by Darcy's law and that in the crack is governed by the Navier–Stokes equations. At the matrix‐fracture domain, the fluid flow is governed by a leakage term derived from Darcy's law. Upon crack propagation, the conservation of mass and energy of the crack fluid is constrained by the isentropic process. We utilize the representative elementary volume‐averaging theory to formulate the mathematical model of the porous matrix, and the drift flux model to formulate the fluid dynamics in the fracture. The numerical solution is conducted using a mixed finite element discretization scheme. The standard Galerkin finite element method is utilized to discretize the diffusive dominant field equations, and the extended finite element method is utilized to discretize the crack propagation, and the fluid leakage at the boundaries between layers of different physical properties. A numerical example is given to demonstrate the computational capability of the model. It shows that the model, despite the relatively large number of degrees of freedom of different physical nature per node, is computationally efficient, and geometry and effectively mesh independent. Copyright © 2017 John Wiley & Sons, Ltd.     

Electroelastic analysis of a piezoelectric cylindrical fiber with a penny-shaped crack embedded in a matrix

The electroelastic response of a penny-shaped crack in a piezoelectric cylindrical fiber embedded in an elastic matrix is investigated in this study. Fourier and Hankel transforms are used to reduce the problem to the solution of a pair of dual integral equations. They are then reduced to a Fredholm integral equation of the second kind. Numerical values on the stress intensity factor, energy release rate and energy density factor for piezoelectric composites are obtained to show the influence of applied electric fields.     

On the truth of nanoscale for nanobeams based on nonlocal elastic stress field theory： equilibrium, governing equation and static deflection

This paper has successfully addressed three critical but overlooked issues in nonlocal elastic stress field theory for nanobeams： （i） why does the presence of increasing nonlocal effects induce reduced nanostructural stiffness in many, but not consistently for all, cases of study, i.e., increasing static deflection, decreasing natural frequency and decreasing buckling load, although physical intuition according to the nonlocal elasticity field theory first established by Eringen tells otherwise？ （ii） the intriguing conclusion that nanoscale effects are missing in the solutions in many exemplary cases of study, e.g., bending deflection of a cantilever nanobeam with a point load at its tip; and （iii） the non-existence of additional higher-order boundary conditions for a higher-order governing differential equation. Applying the nonlocal elasticity field theory in nanomechanics and an exact variational principal approach, we derive the new equilibrium conditions, do- main governing differential equation and boundary conditions for bending of nanobeams. These equations and conditions involve essential higher-order differential terms which are opposite in sign with respect to the previously studies in the statics and dynamics of nonlocal nano-structures. The difference in higher-order terms results in reverse trends of nanoscale effects with respect to the conclusion of this paper. Effectively, this paper reports new equilibrium conditions, governing differential equation and boundary condi- tions and the true basic static responses for bending of nanobeams. It is also concluded that the widely accepted equilibrium conditions of nonlocal nanostructures are in fact not in equilibrium, but they can be made perfect should the nonlocal bending moment be replaced by an effective nonlocal bending moment. These conclusions are substantiated, in a general sense, by other approaches in nanostructural models such as strain gradient theory, modified couple stress models and experiments.     

Penny shaped crack in a piezoelectric cylinder surrounded by an elastic medium subjected to combined in-plane mechanical and electrical loads

The fracture problem of a penny shaped crack in a piezoelectric ceramic cylinder surrounded by an infinite elastic medium under in-plane normal mechanical and electrical loads is considered with the electric continuous boundary conditions on the crack surface. By using the potential theory and Hankel transform, a system of dual integral equations is obtained, and expressed to a Fredholm integral equation of the second kind. The mechanical and electrical field equations and all sorts of field intensity factors of mode I are obtained, and the numerical values of various field intensity factors for PZT-6B piezoelectric ceramic surrounded by several different elastic media are graphically shown for a uniform load and a ring-shaped load, respectively. And the effects of the size of the piezoelectric cylinder and the elastic material properties on various field intensity factors are obtained.     

Penny-shaped crack in a fiber-reinforced matrix

Using the slender inclusion model developed earlier the elastostatic interaction problem between a penny-shaped crack and elastic fibers in an elastic matrix is formulated. For a single set and for multiple sets of fibers oriented perpendicularly to the plane of the crack and distributed symmetrically on concentric circles the problem is reduced to a system of singular integral equations. Techniques for the regularization and for the numerical solution of the system are outlined. For various fiber geometries numerical examples are given and distribution of the stress intensity factor along the crack border is obtained. Sample results showing the distribution of the fiber stress and a measure of the fiber-matrix interface shear are also included.     

弱界面复合材料中的基体裂纹

利用二维弹性力学模型研究了纤维增强复合材料中基体裂纹与弱界面的相互作用机理．文中首先导出各向异性弹性多层介质中刃型位错的基本解，然后运用这些基本解建立了弱界面复合材料中典型的Ｈ型缺陷的奇异积分方程组，通过求解这些方程得到外载荷的大小、弱界面的结合强度、界面的残余压力和摩擦系数、纤维与基体的弹性模量比等微结构参量与基体裂纹附近的应力场的关系     

纤维增强韧性基体界面力学行为

分析了纤维增强韧性基体的界面力学行为及其失效机理,按剪滞理论和应变理化规律研究微复合材料的弹塑性变形和应力状态,讨论了幂硬化和线性硬化基体的弹塑性变形和界面应力分布,并给出纤维应力和位移的表达式。按最大剪应力强度理论建立了纤维／基体界面失效准则,推导出弹塑性界面失效的平均剪应力随纤维埋入长度的变化关系。     

Functionally graded penny-shaped cracks under dynamic loading

Dynamic load is applied to a functionally graded material with penny-shaped cracks. The materials are also transversely isotropic depending only on the axial coordinate z. The elastic region may be regarded to consist of many thin layers such that properties are constants within each layer, but they may vary from layer to layer. Laplace and Hankel transform are used in conjunction with the stiffness matrix approach. The Dual integral equations are then obtained by application of appropriate boundary and interface conditions. Stress intensity factors are then determined in the Laplace transform domain. Inversion yields the results in the time domain. Numerical examples show that multiple crack configurations in functionally graded materials can be treated where the continuously varied material properties can be divided into a finite number of layers with different properties.     

多场作用下输流单层碳纳米管的动力学特性

以非局部弹性理论为基础,采用欧拉-伯努利梁模型,考虑管型区域内滑移边界条件以及碳纳米管的小尺度效应,应用哈密顿原理获得了温度场与轴向磁场共同作用下的输流单层固支碳纳米管（SWCNT）的振动控制方程以及边界条件,依靠微分变换法（DTM法）对此高阶偏微分方程进行求解,通过数值计算研究了多场中单层固支输流碳纳米管的振动与失稳问题。结果表明：温度场、轴向磁场强度、Knudsen数及小尺度参数都会对系统振动频率以及失稳临界流速产生影响。