In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approac...In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731–2755) are discussed also.展开更多
The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties....The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.展开更多
The numerical solutions for uncertain viscoelastic problems have important theo- retical and practical significance. The paper develops a new approach by combining the scaled boundary finite element method （SBFEM） a...The numerical solutions for uncertain viscoelastic problems have important theo- retical and practical significance. The paper develops a new approach by combining the scaled boundary finite element method （SBFEM） and fuzzy arithmetic. For the viscoelastic problems with zero uncertainty, the SBFEM and the temporally piecewise adaptive algorithm is employed in the space domain and the time domain, respectively, in order to provide an accurate semi- analytical boundary-based approach and to ensure the accuracy of discretization in the time domain with different sizes of time step at the same time. The fuzzy arithmetic is used to address the uncertainty analysis of viscoelastic material parameters, and the transformation method is used for computation with the advantages of effectively avoiding overestimation and reducing the computational costs. Numerical examples are provided to test the performance of the proposed method. By comparing with the analytical solutions and the Monte Carlo method, satisfactory results are achieved.展开更多
The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special fe...The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.展开更多
The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this pap...The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.展开更多
The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to pre...The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.展开更多
In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. Thi...The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.展开更多
The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the...The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the model order by scaling the boundary solution onto the inner element.To this end,tri-lateral elements are emanated from a scaling center,followed by the development of a semi-analytical solution along the radial direction and a finite element solution along the circumferential/boundary direction.The discretization is thus limited to the boundaries of the model,and the semi-analytical radial solution is found through the solution of an eigenvalue problem,which restricts the methods’applicability to heterogeneous media.In this research,we first extracted the SBFEM formulation considering the heterogeneity of the media.Then,we replaced the semi-analytical radial solution with the EqBFs and removed the eigenvalue solution step from the SBFEM.The varying coefficients of the partial differential equation(PDE)resulting from the heterogeneity of the media are replaced by a finite series in the radial and circumferential directions of the element.A weighted residual approach is applied to the radial equation.The equilibrated radial solution series is used in the new formulation of the SBFEM.展开更多
The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The co...The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.展开更多
Using a modified 3D random representative volume（RV）finite element model,the effects of model dimensions（impact region and interval between impact and representative regions）,model shapes（rectangular,square,and c...Using a modified 3D random representative volume（RV）finite element model,the effects of model dimensions（impact region and interval between impact and representative regions）,model shapes（rectangular,square,and circular）,and peening-induced thermal softening on resultant critical quantities（residual stress,Almen intensity,coverage,and arc height）after shot peening are systematically examined.A new quantity,i.e.,the interval between impact and representative regions,is introduced and its optimal value is first determined to eliminate any boundary effect on shot peening results.Then,model dimensions are respectively assessed for all model shapes to reflect the actual shot peening process,based on which shape-independent critical shot peening quantities are obtained.Further,it is found that thermal softening of the target material due to shot peening leads to variances of the surface residual stress and arc height,demonstrating the necessity of considering the thermal effect in a constitutive material model of shot peeing.Our study clarifies some of the finite element modeling aspects and lays the ground for accurate modeling of the SP process.展开更多
Repeated Unit Cell（RUC）is a useful tool in micromechanical analysis of composites using Displacement-based Finite Element（DFE）method,and merely applying Periodic Displacement Boundary Conditions（PDBCs）to RUC is ...Repeated Unit Cell（RUC）is a useful tool in micromechanical analysis of composites using Displacement-based Finite Element（DFE）method,and merely applying Periodic Displacement Boundary Conditions（PDBCs）to RUC is almost a standard practice to conduct such analysis.Two basic questions arising from this practice are whether Periodic Traction Boundary Conditions（PTBCs,also known as traction continuity conditions）are guaranteed and whether the solution is independent of selection of RUCs.This paper presents the theoretical aspects to tackle these questions,which unify the strong form,weak form and DFE method of the micromechanical problem together.Specifically,the solution’s independence of selection of RUCs is dealt with on the strong form side,PTBCs are derived from the weak form as natural boundary conditions,and the validity of merely applying PDBCs in micromechanical Finite Element（FE）analysis is proved by referring to its intrinsic connection to the strong form and weak form.Key points in the theoretical aspects are demonstrated by illustrative examples,and the merits of setting micromechanical FE analysis under the background of a clear theoretical framework are highlighted in the efficient selection of RUCs for Uni Directional（UD）fiber-reinforced composites.展开更多
In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equatio...In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equation method. The numerical approximate schemes for both problems on a bounded or unbounded domain in R3 are proposed and their prior error estimates are obtained.展开更多
For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has ...For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has the accuracy O（h^min{2k,k＋4}） The theoretical analysis coincides the reported numerical results.展开更多
Stress intensity factors （SIFs） for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method （SBFEM）. In the SBFEM, the analytical adva...Stress intensity factors （SIFs） for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method （SBFEM）. In the SBFEM, the analytical advantage of the solution in the radial direction allows SIFs to be directly determined from its definition, therefore no special crack-tip treatment is necessary. Furthermore anisotropic material behavior can be treated easily. Different distributions of surface tractions are considered for the center and double-edge-cracked disks. The benchmark examples are modeled and an excellent agreement between the results in the present study and those in published literature is found. It shows that SBFEM is effective and possesses high accuracy. The SIFs of the cracked orthotropic material circular disks subjected to different surface tractions are also evaluated. The technique of substructure is applied to handle the multiple cracks problem.展开更多
The stress intensity factors （SIF） considering arbitrarily distributed surface tractions are evaluated based on the sealed boundary finite element method （SBFEM）. The semi-analytical solving process for the stress...The stress intensity factors （SIF） considering arbitrarily distributed surface tractions are evaluated based on the sealed boundary finite element method （SBFEM）. The semi-analytical solving process for the stress intensity factors including the effects of surface tractions is presented. Provided are the numerical examples for the evaluation of mode I and Ⅱ stress intensity factors with linear and non-linear distributing forces loaded on the crack surfaces. The crack problems of anisotropy and bimaterial interface are also studied and the stress intensity factors of single-edge-cracked orthotropic material and bi-material interface problems with surface tractions are calculated. Comparisons with the analytical solutions show that the proposed approach is effective and possesses high accuracy.展开更多
A numerical model based on the scaled boundary finite element method is devel- oped to simulate the hydraulic fracturing in concrete-like quasi-brittle materials using cohesive interface elements. The shadow domain me...A numerical model based on the scaled boundary finite element method is devel- oped to simulate the hydraulic fracturing in concrete-like quasi-brittle materials using cohesive interface elements. The shadow domain method developed previously （Yang and Deeks in Eng Fract Mech 143（4）：333 354, 2007） is extended to consider crack-width-dependent hydraulic pres- sure and cohesive traction, so that the stress intensity factors caused by both crack-face forces are semi-analytically calculated separately in the same way. The crack propagation is determined by the criterion of KI ≥ 0, and the propagation direction by the linear elastic fracture mechanics criteria. Two examples of concrete structures are modeled, and the results are in good agreement with the experimental data and others numerical results.展开更多
Roll flattening theory is an important part of plate shape control theories for 20-high mill. In order to improve the accuracy of roll flattening calculation for 20-high mill, a new and more accurate roll flattening m...Roll flattening theory is an important part of plate shape control theories for 20-high mill. In order to improve the accuracy of roll flattening calculation for 20-high mill, a new and more accurate roll flattening model was proposed. In this model, the roll barrel was considered as a finite length semi-infinite body. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distributed force was obtained and an accurate roll flattening model was established. Coupled with roll bending model and strip plastic deformation, a new and more accurate plate control model for 20-high mill was established. Moreover, the effects of the first intermediate roll taper angle and taper length were analyzed. The tension distribution calculated by analytical model was consistent with the experimental results.展开更多
In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite el...In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite element method （SBFEM）, which is a novel semi-analytical method with the advantages of combining the finite element method （FEM） with the boundary element method （BEM）. The whole solution domain is divided into one unbounded sub-domain and one bounded sub-domain by the exterior cylinder. By weakening the governing differential equation in the circumferential direction, the SBFEM equations for both domains can be solved analytically in the radial direction. Only the boundary on the circumference of the exterior porous cylinder is discretized with curved surface finite elements. Meanwhile, by introducing a variable porous-effect parameter G, non-homogeneous materials caused by the complex configuration of the exterior cylinder are modeled without additional efforts. Comparisons clearly demonstrate the excellent accuracy and computational efficiency associated with the present SBFEM. The effects of the wide range wave parameters and the structure configuration are examined. This parametric study will help determine the various hydrodynamic effects of the concentric porous cylindrical structure.展开更多
The paper presents the variational formulation and well posedness of the coupling method offinite elements and boundary elements for radiation problem. The convergence and optimal errorestimate for the approximate sol...The paper presents the variational formulation and well posedness of the coupling method offinite elements and boundary elements for radiation problem. The convergence and optimal errorestimate for the approximate solution and numerical experiment are provided.展开更多
基金The project supported by the National Natural Science Foundation of China(10202023 and 10272103)the Key Project of CAS(KJCX2-SW-L2).
文摘In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731–2755) are discussed also.
基金This research wasfinanciallysupported bythe National Natural Science Foundation of China(Grant No.50639030)a Programfor Changjiang ScholarsInnovative Research Teamin Dalian University of Technology(Grant No.IRTO420)
文摘The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.
文摘The numerical solutions for uncertain viscoelastic problems have important theo- retical and practical significance. The paper develops a new approach by combining the scaled boundary finite element method （SBFEM） and fuzzy arithmetic. For the viscoelastic problems with zero uncertainty, the SBFEM and the temporally piecewise adaptive algorithm is employed in the space domain and the time domain, respectively, in order to provide an accurate semi- analytical boundary-based approach and to ensure the accuracy of discretization in the time domain with different sizes of time step at the same time. The fuzzy arithmetic is used to address the uncertainty analysis of viscoelastic material parameters, and the transformation method is used for computation with the advantages of effectively avoiding overestimation and reducing the computational costs. Numerical examples are provided to test the performance of the proposed method. By comparing with the analytical solutions and the Monte Carlo method, satisfactory results are achieved.
基金The project supported by the National Natural Science Foundation of China (50579081)the Australian Research Council (DP0452681)The English text was polished by Keren Wang
文摘The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.
基金support by the National Natural Science Foundation of China(grant No.51779033,51409038)the National Key Research and Development Plan(grant No.2016YFB0201001)the National Natural Science Foundation of China(grant No.51421064)
文摘The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.
基金Supported by the Key Program of National Natural Science Foundation of China(No.51138001)the Science Fund for Creative Research Groups of National Natural Science Foundation of China(No.51121005)+2 种基金the Fundamental Research Funds for the Central Universities(DUT13LK16)the Young Scientists Fund of National Natural Science Foundation of China(No.51109134)China Postdoctoral Science Foundation(No.2011M500814)
文摘The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
基金supported by the National Natural Science Foundation of China(No.11002054)
文摘The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.
文摘The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the model order by scaling the boundary solution onto the inner element.To this end,tri-lateral elements are emanated from a scaling center,followed by the development of a semi-analytical solution along the radial direction and a finite element solution along the circumferential/boundary direction.The discretization is thus limited to the boundaries of the model,and the semi-analytical radial solution is found through the solution of an eigenvalue problem,which restricts the methods’applicability to heterogeneous media.In this research,we first extracted the SBFEM formulation considering the heterogeneity of the media.Then,we replaced the semi-analytical radial solution with the EqBFs and removed the eigenvalue solution step from the SBFEM.The varying coefficients of the partial differential equation(PDE)resulting from the heterogeneity of the media are replaced by a finite series in the radial and circumferential directions of the element.A weighted residual approach is applied to the radial equation.The equilibrated radial solution series is used in the new formulation of the SBFEM.
基金National Natural Science Foundation of China (50608024 and 50538050).
文摘The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.
基金the financial support from China Scholarship Council (CSC) (No. 201406025083)National Natural Science Foundation of China (NSFC) (Nos. 51305012 and 51675024)+3 种基金Aviation Science Fund of China (No. 2014ZB51)financial support from NSFC (No. 51375031)financial support from NSFC (No. 51628101)National Sciences and Engineering Research Council (NSERC) Discovery grant (No. RGPIN 418469-2012)
文摘Using a modified 3D random representative volume（RV）finite element model,the effects of model dimensions（impact region and interval between impact and representative regions）,model shapes（rectangular,square,and circular）,and peening-induced thermal softening on resultant critical quantities（residual stress,Almen intensity,coverage,and arc height）after shot peening are systematically examined.A new quantity,i.e.,the interval between impact and representative regions,is introduced and its optimal value is first determined to eliminate any boundary effect on shot peening results.Then,model dimensions are respectively assessed for all model shapes to reflect the actual shot peening process,based on which shape-independent critical shot peening quantities are obtained.Further,it is found that thermal softening of the target material due to shot peening leads to variances of the surface residual stress and arc height,demonstrating the necessity of considering the thermal effect in a constitutive material model of shot peeing.Our study clarifies some of the finite element modeling aspects and lays the ground for accurate modeling of the SP process.
文摘Repeated Unit Cell（RUC）is a useful tool in micromechanical analysis of composites using Displacement-based Finite Element（DFE）method,and merely applying Periodic Displacement Boundary Conditions（PDBCs）to RUC is almost a standard practice to conduct such analysis.Two basic questions arising from this practice are whether Periodic Traction Boundary Conditions（PTBCs,also known as traction continuity conditions）are guaranteed and whether the solution is independent of selection of RUCs.This paper presents the theoretical aspects to tackle these questions,which unify the strong form,weak form and DFE method of the micromechanical problem together.Specifically,the solution’s independence of selection of RUCs is dealt with on the strong form side,PTBCs are derived from the weak form as natural boundary conditions,and the validity of merely applying PDBCs in micromechanical Finite Element（FE）analysis is proved by referring to its intrinsic connection to the strong form and weak form.Key points in the theoretical aspects are demonstrated by illustrative examples,and the merits of setting micromechanical FE analysis under the background of a clear theoretical framework are highlighted in the efficient selection of RUCs for Uni Directional（UD）fiber-reinforced composites.
基金This research was supported by the National Natural Science Foundation of China
文摘In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equation method. The numerical approximate schemes for both problems on a bounded or unbounded domain in R3 are proposed and their prior error estimates are obtained.
基金Project supported by the National Natural Science Foundation of China (Nos. 10571046, 10371038)
文摘For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has the accuracy O（h^min{2k,k＋4}） The theoretical analysis coincides the reported numerical results.
基金financially supported by the National Natural Youth Foundation of China (Grant Nos. 51109134,51009019, 11102118 and 51208310)the Liaoning Province Education Administration Foundation (Grant No. L2010413)+1 种基金the China Postdoctoral Science Foundation (Grant No. 2011M500557)the Natural Science Foundation of Liaoning Province (Grant No.20102164)
文摘Stress intensity factors （SIFs） for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method （SBFEM）. In the SBFEM, the analytical advantage of the solution in the radial direction allows SIFs to be directly determined from its definition, therefore no special crack-tip treatment is necessary. Furthermore anisotropic material behavior can be treated easily. Different distributions of surface tractions are considered for the center and double-edge-cracked disks. The benchmark examples are modeled and an excellent agreement between the results in the present study and those in published literature is found. It shows that SBFEM is effective and possesses high accuracy. The SIFs of the cracked orthotropic material circular disks subjected to different surface tractions are also evaluated. The technique of substructure is applied to handle the multiple cracks problem.
基金The present research workis financially supported by the National Natural Science Foundation of China (Grant No90510018)China Postdoctorial Science Foundation (Grant No20060390985)
文摘The stress intensity factors （SIF） considering arbitrarily distributed surface tractions are evaluated based on the sealed boundary finite element method （SBFEM）. The semi-analytical solving process for the stress intensity factors including the effects of surface tractions is presented. Provided are the numerical examples for the evaluation of mode I and Ⅱ stress intensity factors with linear and non-linear distributing forces loaded on the crack surfaces. The crack problems of anisotropy and bimaterial interface are also studied and the stress intensity factors of single-edge-cracked orthotropic material and bi-material interface problems with surface tractions are calculated. Comparisons with the analytical solutions show that the proposed approach is effective and possesses high accuracy.
基金This research is funded by the National Natural Science Foundation of China （Nos. 51779222 and 51378461）, Zhejiang Provincial Natural Science Foundation of China （No. LR14E080002）, and the Fundamental Research Funds for the Central Universities （No. 2017QNA4027）.
文摘A numerical model based on the scaled boundary finite element method is devel- oped to simulate the hydraulic fracturing in concrete-like quasi-brittle materials using cohesive interface elements. The shadow domain method developed previously （Yang and Deeks in Eng Fract Mech 143（4）：333 354, 2007） is extended to consider crack-width-dependent hydraulic pres- sure and cohesive traction, so that the stress intensity factors caused by both crack-face forces are semi-analytically calculated separately in the same way. The crack propagation is determined by the criterion of KI ≥ 0, and the propagation direction by the linear elastic fracture mechanics criteria. Two examples of concrete structures are modeled, and the results are in good agreement with the experimental data and others numerical results.
基金Item Sponsored by National Natural Science Foundation of China(51474190)Natural Sceince Foundation of Hebei Province of China(E2015203311)
文摘Roll flattening theory is an important part of plate shape control theories for 20-high mill. In order to improve the accuracy of roll flattening calculation for 20-high mill, a new and more accurate roll flattening model was proposed. In this model, the roll barrel was considered as a finite length semi-infinite body. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distributed force was obtained and an accurate roll flattening model was established. Coupled with roll bending model and strip plastic deformation, a new and more accurate plate control model for 20-high mill was established. Moreover, the effects of the first intermediate roll taper angle and taper length were analyzed. The tension distribution calculated by analytical model was consistent with the experimental results.
基金supported by the State Key Program of the National Natural Science Foundation of China(Grant No.51138001)China-Germany joint research project(Grant No.GZ566)Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering(Grant No.shlhse-2010-C-03)
文摘In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite element method （SBFEM）, which is a novel semi-analytical method with the advantages of combining the finite element method （FEM） with the boundary element method （BEM）. The whole solution domain is divided into one unbounded sub-domain and one bounded sub-domain by the exterior cylinder. By weakening the governing differential equation in the circumferential direction, the SBFEM equations for both domains can be solved analytically in the radial direction. Only the boundary on the circumference of the exterior porous cylinder is discretized with curved surface finite elements. Meanwhile, by introducing a variable porous-effect parameter G, non-homogeneous materials caused by the complex configuration of the exterior cylinder are modeled without additional efforts. Comparisons clearly demonstrate the excellent accuracy and computational efficiency associated with the present SBFEM. The effects of the wide range wave parameters and the structure configuration are examined. This parametric study will help determine the various hydrodynamic effects of the concentric porous cylindrical structure.
基金This research was supported in part by the Institute for Mathematics and its applications with funds provided by NSF, USA
文摘The paper presents the variational formulation and well posedness of the coupling method offinite elements and boundary elements for radiation problem. The convergence and optimal errorestimate for the approximate solution and numerical experiment are provided.
