首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
An interaction (energy) integral is derived for the computation of mixed-mode stress intensity factors (SIFs) in nonhomogeneous materials with continuous or discontinuous properties. This method is based on a conservation integral that relies on two admissible mechanical states (actual and auxiliary fields). In general, the interaction energy contour integral is converted into an equivalent domain integral in numerical computations. It can be seen from the equivalent domain integral, the integrand does not involve any derivatives of material properties. Moreover, the formulation can be proved valid even when the integral domain contains material interfaces. Therefore, it is not necessary to limit the material properties to be continuous for the present method. Due to these advantages the application range of the interaction integral method can be greatly enlarged. The numerical implementation of the derived expression is combined with the extended finite element method (XFEM). Using this method, the influences of material properties on the mixed-mode SIFs are investigated for four types of material properties selected in this work. Numerical results show that the mechanical properties and their first-order derivatives can affect mode I and II SIFs greatly, while the higher-order derivatives affect the SIFs very slightly.  相似文献   

2.
Magneto-electro-elastic (MEE) materials usually consist of piezoelectric (PE) and piezomagnetic (PM) phases. Between different constituent phases, there exist lots of interfaces with discontinuous MEE properties. Complex interface distribution brings a great difficulty to the fracture analysis of MEE materials since the present fracture mechanics methods can hardly solve the fracture parameters efficiently of a crack surrounded by complex interfaces. This paper develops a new domain formulation of the interaction integral for the computation of the fracture parameters including stress intensity factors (SIFs), electric displacement intensity factor (EDIF) and magnetic induction intensity factor (MIIF) for linear MEE materials. The formulation derived here does not involve any derivatives of material properties and moreover, it can be proved that an arbitrary interface in the integral domain does not affect the validity and the value of the interaction integral. Namely, the interaction integral is domain-independent for material interfaces and thus, its application does not require material parameters to be continuous. Due to this advantage, the interaction integral becomes an effective approach for extracting the fracture parameters of MEE materials with complex interfaces. Combined with the extended finite element method (XFEM), the interaction integral is employed to solve several representative problems to verify its accuracy and domain-independence. Good results show the effectiveness of the present method in the fracture analysis of MEE materials with continuous and discontinuous properties. Finally, the particulate MEE composites composed of PE and PM phases are considered and four schemes of different property-homogenization level are proposed for comparing their effectiveness.  相似文献   

3.
A plane crack problem of nonhomogeneous materials with interfaces subjected to static thermal loading is investigated. A modified interaction energy integral method (IEIM) is developed to obtain the mixed-mode thermal stress intensity factors (TSIFs). Compared with the previous IEIM, the original point of this paper is: the domain-independence of the modified IEIM still stands in nonhomogeneous materials with interfaces under thermal loading. Therefore, the modified IEIM can still be applied to obtain the TSIFs of nonhomogeneous material even if the integral domain includes interfaces. The modified IEIM is combined with the extended finite element method (XFEM) to solve several thermal fracture problems of nonhomogeneous materials. Good agreement can be obtained compared with the analytic solutions and the domain-independence of the IEIM is verified. Therefore, the present method is effective to study the TSIFs of nonhomogeneous materials even when the materials contain interfaces. The influence of the discontinuity of the material properties (thermal expansion coefficient, thermal conductivity and Young’s modulus) on the TSIFs is investigated. The results show that the discontinuity of both thermal expansion coefficient and Young’s modulus affects the TSIFs greatly, while the discontinuity of thermal conductivity does not arouse obvious change of the TSIFs.  相似文献   

4.
Piezoelectric materials and structures contain more or less electromechanical interfaces in engineering applications. It is difficult to obtain the fracture parameters efficiently of the piezoelectric materials with complex interfaces. This paper presents a domain-independent interaction integral for material nonhomogeneity and discontinuity which can be used for solving the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) of piezoelectric materials with complex interfaces efficiently. The interaction integral is based on the J-integral by superimposition of two admissible states and the present formulation does not involve any derivatives of mechanical and electric properties. Moreover, it is proved that the interface in the integral domain does not affect the value of the interaction integral and thus, the present method does not require electromechanical parameters of piezoelectric materials to be continuous. The interaction integral method combined with the extended finite element method (XFEM) is used to investigate the influences of material continuity on the SIF and the EDIF and the results show that the material parameters and their first-order derivatives affect both the SIF and the EDIF greatly, while the higher-order derivatives affect both of them slightly.  相似文献   

5.
This paper presents a modified interaction energy integral method to analyze the thermal stress intensity factors (TSIFs) and electric displacement intensity factor (EDIF) in nonhomogeneous piezoelectric materials under thermal loading. This modified method is demonstrated to be domain-independent, even when the nonhomogeneous piezoelectric materials contain interfaces with thermo-electro-mechanical properties. As a result, the method is shown to be convenient for determining the TSIFs and EDIF in nonhomogeneous piezoelectric materials with interfaces. Several examples are shown, and they successfully verify the domain-independence of the present interaction energy integral. The study results also show that the mismatch of material properties can significantly influence the TSIFs and EDIF, particularly when the crack tip is close to the interface. Crack angles and temperature boundary conditions are also shown to significantly influence the TSIFs and EDIF.  相似文献   

6.
A new model, piecewise-exponential model (PE model), is developed to investigate the crack problem of the functionally graded materials (FGMs) with arbitrary properties. In the PE model, the functionally graded material is divided into some nonhomogeneous layers along the gradient direction of the properties, with each layer’s properties varying exponentially. By this way, the real material properties can be approached by a series of exponential functions. Since the real material properties are used on both surfaces of each nonhomogeneous layer, the nature of continuously varying properties of FGMs can be approached accurately. The influences of the local nonhomogeneity on the crack-tip fields can be fully considered. By using the new model, the fracture problem of a functionally graded strip with arbitrary properties and a crack vertical to the free surfaces is studied. The integral transform method, the theory of residues and the theory of singular integral equation are applied. Some representative samples with different kinds of nonhomogeneous properties are analyzed and the corresponding stress intensity factors (SIFs) are presented. It is shown that the PE mode is effective for investigating the crack problems of the FGMs with arbitrary properties.  相似文献   

7.
A strategy for the mixed-dimensional coupling of finite shell elements and 3D boundary elements is presented. The stiffness formulation for the boundary element domain is generated by the 3D symmetric Galerkin boundary element method and thus can be assembled to the global finite element formulation. Based on the equality of work at the coupling interface, coupling equations in an integral sense are derived for curved coupling interfaces and formulated as multipoint constraints in terms of kinematic quantities. Several examples show the highly accurate results compared to a strict kinematic coupling technique.  相似文献   

8.
In the present work, graded finite element and boundary element methods capable of modeling behaviors of structures made of nonhomogeneous functionally graded materials (FGMs) composed of two constituent phases are presented. A numerical implementation of Somigliana’s identity in two-dimensional displacement fields of the isotropic nonhomogeneous problems is presented using the graded elements. Based on the constitutive and governing equations and the weighted residual technique, effective boundary element formulations are implemented for elastic nonhomogeneous isotropic solid models. Results of the finite element method are derived based on a Rayleigh–Ritz energy formulation. The heterogeneous structures are made of combined ceramic–metal materials, in which the material properties vary continuously along the in-plane or thickness directions according to a power law. To verify the present work, three numerical examples are provided in the paper.  相似文献   

9.
Thermal effects are well known to manifest themselves as additional volume integral terms in the direct formulation of the boundary integral equation (BIE) for linear elastic solids when using the boundary element method (BEM). This domain integral has been successfully transformed in an exact manner to surface ones only in isotropy and in 2D anisotropy, thereby restoring the BEM as a truly boundary solution technique. The difficulties with extending it to 3D general anisotropic solids lie in the mathematical complexity of the Green’s function and its derivatives for such materials. These quantities are required items in the BEM formulation. In this paper, the exact, analytical transformation of the volume integral associated with thermal effects to surface ones is achieved for a transversely isotropic material using a similar approach which the authors have previously employed for the same task in BEM for 2D general anisotropy. A numerical scheme, however, needs to be employed to evaluate some of the new terms introduced in the surface integrals that arise from this process here. The mathematical soundness of the formulation is demonstrated by a few examples; the numerical results obtained are checked by alternative means, including those obtained from the commercial FEM code, ANSYS.  相似文献   

10.
采用一种新型的杂交元模型和一种单胞模型来解决周期分布多边形夹杂角部的奇异性应力相互干涉的问题。新型杂交元模型是基于广义Hellinger-Reissner变分原理建立的,其中奇异性应力场分量和位移场分量是采用有限元特征分析法的数值特征解得到的。使用当前的新型杂交元模型,只需要在夹杂角部邻域的周界上划分一维单元,避免了像传统有限元模型那样需要划分高密度二维单元。文中给出了代表奇异性应力场强度的夹杂角部广义应力强度因子数值解,并考虑材料属性、夹杂尺寸和夹杂位置关系的影响。算例中,考虑了夹杂和基体完全接合的情况,并给出了考核例。结果表明:当前模型能得到高精度数值解,且收敛性好;与传统有限元法和积分方程方法相比,该模型更具有通用性,为非均质材料的细观力学分析打下了基础。  相似文献   

11.
Anti-plane dynamic fracture analysis is presented for functionally graded materials (FGM) with arbitrary spatial variations of material properties. The FGM with the material properties varying continuously in an arbitrary manner is modeled as a multi-layered medium with the elastic modulus and mass density varying linearly in each sub-layer and continuous at the interfaces between two adjacent sub-layers. With this linearly inhomogeneous multi-layered model, the problem of a crack in a graded interfacial zone bonded to two homogeneous half-spaces or in a coating bonded to a homogeneous half-space subjected to the anti-plane shear impact load is investigated. Laplace and Fourier transforms and transfer matrix are applied to reduce the associated mixed boundary value problem to a Cauchy singular integral equation which is solved numerically in the Laplace transformed domain. The dynamic stress intensity factors (DSIF) are obtained by using the numerical technique of Laplace inversion.  相似文献   

12.
梯度材料中矩形裂纹的对偶边界元方法分析   总被引:2,自引:0,他引:2  
肖洪天  岳中琦 《力学学报》2008,40(6):840-848
采用对偶边界元方法分析了梯度材料中的矩形裂纹. 该方法基于层状材料基本解,以非裂纹边界的位移和面力以及裂纹面的间断位移作为未知量. 位移边界积分方程的源点配置在非裂纹边界上,面力边界积分方程的源点配置在裂纹面上. 发展了边界积分方程中不同类型奇异积分的数值方法. 借助层状材料基本解,采用分层方法逼近梯度材料夹层沿厚度方向力学参数的变化. 与均匀介质中矩形裂纹的数值解对比,建议方法可以获得高精度的计算结果. 最后,分析了梯度材料中均匀张应力作用下矩形裂纹的应力强度因子,讨论了梯度材料非均匀参数、夹层厚度和裂纹与夹层之间相对位置对应力强度因子的影响.   相似文献   

13.
Finite elements with graded properties are used to simulate elastic wave propagation in functionally graded materials. The graded elements are formulated with continuously nonhomogeneous material property fields and compared to conventionally formulated homogeneous elements. An example problem is solved for the two-dimensional case to show the potential benefits of the graded formulation. It is observed that the conventional elements give a discontinuous stress field in the direction perpendicular to the material property gradient, while the graded elements give a continuous distribution. In the one-dimensional case, the solutions are compared to the analytical solutions by Chiu and Erdogan [J. Sound Vib. 222 (3) (1999) 453]. The results show that for identical levels of mesh refinement, the graded formulation produces similar spatial resolution, and temporal resolution for the range of boundary value problems studied. Observations are explained and their implications for numerical modeling are discussed.  相似文献   

14.
The interaction integral is an accurate and robust scheme for evaluating mixed-mode stress intensity factors. This paper extends the concept to orthotropic functionally graded materials and addresses fracture mechanics problems with arbitrarily oriented straight and/or curved cracks. The gradation of orthotropic material properties are smooth functions of spatial coordinates, which are integrated into the element stiffness matrix using the so-called “generalized isoparametric formulation”. The types of orthotropic material gradation considered include exponential, radial, and hyperbolic-tangent functions. Stress intensity factors for mode I and mixed-mode two-dimensional problems are evaluated by means of the interaction integral and the finite element method. Extensive computational experiments have been performed to validate the proposed formulation. The accuracy of numerical results is discussed by comparison with available analytical, semi-analytical, or numerical solutions.  相似文献   

15.
功能梯度材料结构的热应力边界元分析   总被引:2,自引:1,他引:1  
高效伟  杨恺 《力学学报》2011,43(1):136-143
导出了一种新的可对功能梯度材料结构进行二维和三维热应力分析的积分方程, 利用该方程并结合多区域边界元三步求解技术, 可对由任意多种介质组成的变物性参数复合结构进行热应力分析; 采用与弹性模量无关的Kelvin解作为问题的基本解, 导出的积分方程含有由材料的非均质性以及温度变化引起的域积分; 使用径向积分法将所有的域积分转换成等价的边界积分, 从而建立起只需要边界离散的无内部网格边界元算法; 最后对两种典型的飞行器防热结构进行二维和三维热应力算例分析, 并通过与有限元计算结构对比验证结果的正确性.   相似文献   

16.
Some composite materials are constructed of two dissimilar half-planes bonded by a nonhomogeneous elastic layer. In the present study, a crack is situated at the interface between the upper half-plane and the bonding layer of such a material, and another crack is located at the interface between the lower half-plane and the bonding layer. The material properties of the bonding layer vary continuously from those of the lower half-plane to those of the upper half-plane. Incoming shock stress waves impinge upon the two interface cracks normal to their surfaces. Fourier transformations were used to reduce the boundary conditions for the cracks to two pairs of dual integral equations in the Laplace domain. To solve these equations, the differences in the crack surface displacements were expanded in a series of functions that are zero-valued outside the cracks. The unknown coefficients in the series were solved using the Schmidt method so as to satisfy the conditions inside the cracks. The stress intensity factors were defined in the Laplace domain and were inverted numerically to physical space. Dynamic stress intensity factors were calculated numerically for selected crack configurations.  相似文献   

17.
A new integral equation formulation of two-dimensional infinite isotropic medium (matrix) with various inclusions and cracks is presented in this paper. The proposed integral formulation only contains the unknown displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks. In order to solve the inclusion–crack problems, the displacement integral equation is used when the source points are acting on the inclusion–matrix interfaces, whilst the stress integral equation is adopted when the source points are being on the crack surfaces. Thus, the resulting system of equations can be formulated so that the displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks can be obtained. Based on one point formulation, the stress intensity factors at the crack tips can be achieved. Numerical results from the present method are in excellent agreement with those from the conventional boundary element method.  相似文献   

18.
功能梯度材料动态断裂力学的径向积分边界元法   总被引:1,自引:0,他引:1  
高效伟  郑保敬  刘健 《力学学报》2015,47(5):868-873
采用径向积分边界元法分析功能梯度材料动态断裂力学问题. 该方法使用与弹性模量无关的弹性静力学开尔文基本解作为问题的基本解,在导出的边界-域积分方程中含有由材料的非均质性和惯性项引起的域积分,通过径向积分法将域积分转化为等效的边界积分,得到只含边界积分的纯边界积分方程;从而建立只需边界离散的无内部网格边界元算法. 采用候博特方法求解关于时间二阶导数的系统离散的常微分方程组. 最后通过数值算例验证本文方法的精度和有效性.   相似文献   

19.
黏弹性体界面裂纹的冲击响应   总被引:3,自引:0,他引:3  
研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。  相似文献   

20.
The problem of thermoelastic contact mechanics for the coating/substrate system with functionally graded properties is investigated, where the rigid flat punch is assumed to slide over the surface of the coating involving frictional heat generation. With the coefficient of friction being constant, the inertia effects are neglected and the solution is obtained within the framework of steady-state plane thermoelasticity. The graded material exists as a nonhomogeneous interlayer between dissimilar, homogeneous phases of the coating/substrate system or as a nonhomogeneous coating deposited on the substrate. The material nonhomogeneity is represented by spatially varying thermoelastic moduli expressed in terms of exponential functions. The Fourier integral transform method is employed and the formulation of the current thermoelastic contact problem is reduced to a Cauchy-type singular integral equation of the second kind for the unknown contact pressure. Numerical results include the distributions of the contact pressure and the in-plane component of the surface stress under the prescribed thermoelastic environment for various combinations of geometric, loading, and material parameters of the coated medium. Moreover, in order to quantify and characterize the singular behavior of contact pressure distributions at the edges of the flat punch, the stress intensity factors are defined and evaluated in terms of the solution to the governing integral equation.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号