In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
High energy gas fracturing is a simple approach of applying high pressure gas to stimulate wells by gen- erating several radial cracks without creating any other damages to the wells. In this paper, a numerical algori...High energy gas fracturing is a simple approach of applying high pressure gas to stimulate wells by gen- erating several radial cracks without creating any other damages to the wells. In this paper, a numerical algorithm is proposed to quantitatively simulate propagation of these fractures around a pressurized hole as a quasi-static phenomenon. The gas flow through the cracks is assumed as a one-dimensional transient flow, governed by equations of conservation of mass and momentum. The fractured medium is modeled with the extended finite element method, and the stress intensity factor is calculated by the simple, though sufficiently accurate, displacement ex- trapolation method. To evaluate the proposed algorithm, two field tests are simulated and the unknown parameters are determined through calibration. Sensitivity analyses are performed on the main effective parameters. Considering that the level of uncertainty is very high in these types of engineering problems, the results show a good agreement with the experimental data. They are also consistent with the theory that the final crack length is mainly determined by the gas pressure rather than the initial crack length produced by the stress waves.展开更多
The reliability of real-time embedded software directly determines the reliability of the whole real-time embedded sys- tem, and the effective software testing is an important way to ensure software quality and reliab...The reliability of real-time embedded software directly determines the reliability of the whole real-time embedded sys- tem, and the effective software testing is an important way to ensure software quality and reliability. Based on the analysis of the characteristics of real-time embedded software, the formal method is introduced into the real-time embedded software testing field and the real-time extended finite state machine (RT-EFSM) model is studied firstly. Then, the time zone division method of real-time embedded system is presented and the definition and description methods of time-constrained transition equivalence class (timeCTEC) are presented. Furthermore, the approaches of the testing sequence and test case generation are put forward. Finally, the proposed method is applied to a typical avionics real- time embedded software testing practice and the examples of the timeCTEC, testing sequences and test cases are given. With the analysis of the testing result, the application verification shows that the proposed method can effectively describe the real-time embedded software state transition characteristics and real-time requirements and play the advantages of the formal methods in accuracy, effectiveness and the automation supporting. Combined with the testing platform, the real-time, closed loop and automated simulation testing for real-time embedded software can be realized effectively.展开更多
Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Fi...Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Firstly,the single-edge notched beam test was used to analyze the temperature effect and calculate the material parameters.Then,the mechanical responses were studied using numerical analysis.It is concluded that 5℃ can be selected as the critical temperature that affects the fracture properties,and numerical simulations indicate that crack propagation is found to significantly affect the stress state of the epoxy asphalt mixture.The maximum principal stress at the crack surface exhibits different trends at various temperatures.Numerical solution of stress intensity factor can well meet the theoretical solution,especially when the temperature is lower than 5℃.展开更多
Extended finite element method (XFEM) implementation of the interaction integral methodology for evaluating the stress intensity factors (SIF) of the mixed-mode crack problem is presented. A discontinuous function...Extended finite element method (XFEM) implementation of the interaction integral methodology for evaluating the stress intensity factors (SIF) of the mixed-mode crack problem is presented. A discontinuous function and the near-tip asymptotic function are added to the classic finite element approximation to model the crack behavior. Two-state integral by the superposition of actual and auxiliary fields is derived to calculate the SIFs. Applications of the proposed technique to the inclined centre crack plate with inclined angle from 0° to 90° and slant edge crack plate with slant angle 45°, 67.5° and 90° are presented, and comparisons are made with closed form solutions. The results show that the proposed method is convenient, accurate and computationallv efficient.展开更多
In this paper, the extended finite element method (XFEM) is adopted to analyze the interaction between a single macroscopic inclusion and a single macroscopic crack as well as that between multiple macroscopic or micr...In this paper, the extended finite element method (XFEM) is adopted to analyze the interaction between a single macroscopic inclusion and a single macroscopic crack as well as that between multiple macroscopic or microscopic defects under thermal/mechanical load. The effects of different shapes of multiple inclusions on the material thermomechanical response are investigated, and the level set method is coupled with XFEM to analyze the interaction of multiple defects. Further, the discretized extended finite element approximations in relation to thermoelastic problems of multiple defects under displacement or temperature field are given. Also, the interfaces of cracks or materials are represented by level set functions, which allow the mesh assignment not to conform to crack or material interfaces. Moreover, stress intensity factors of cracks are obtained by the interaction integral method or the M-integral method, and the stress/strain/stiffness fields are simulated in the case of multiple cracks or multiple inclusions. Finally, some numerical examples are provided to demonstrate the accuracy of our proposed method.展开更多
In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination...In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.展开更多
Self-deformation cracking is the cracking caused by thermal deformation, autogenous volume deformation or shrinkage deformation. In this paper, an extended finite element calculation method was deduced for concrete cr...Self-deformation cracking is the cracking caused by thermal deformation, autogenous volume deformation or shrinkage deformation. In this paper, an extended finite element calculation method was deduced for concrete crack propagation under a constant hydration and hardening condition during the construction period, and a corresponding programming code was developed. The experimental investigation shows that initial crack propagation caused by self-deformation loads can be analyzed by this program. This improved algorithm was a preliminary application of the XFEM to the problem of the concrete self-deformation cracking during the hydration and hardening period. However, room for improvement exists for this algorithm in terms of matching calculation programs with mass concrete temperature fields containing cooling pipes and the influence of creep or damage on crack propagation.展开更多
This study develops a mesoscopic framework and methodology for the modeling of microcracks in concrete. A new algorithm is first proposed for the generation of random concrete meso-structure including microcracks and ...This study develops a mesoscopic framework and methodology for the modeling of microcracks in concrete. A new algorithm is first proposed for the generation of random concrete meso-structure including microcracks and then coupled with the extended finite element method to simulate the heterogeneities and discontinuities present in the meso-structure of concrete. The proposed procedure is verified and exemplified by a series of numerical simulations. The simulation results show that microcracks can exert considerable impact on the fracture performance of concrete. More broadly, this work provides valuable insight into the initiation and propagation mechanism of microcracks in concrete and helps to foster a better understanding of the micro-mechanical behavior of cementitious materials.展开更多
In this paper,the problem of modeling crack in 2D viscoelastic media is studied using the extended finite element method.The paper focuses on the definition of enrichment functions suitable for cracks assessment in vi...In this paper,the problem of modeling crack in 2D viscoelastic media is studied using the extended finite element method.The paper focuses on the definition of enrichment functions suitable for cracks assessment in viscoelastic media and the generalized domain integrals used in the determination of crack tip parameters.The opening mode and mixed mode solutions of crack tip fracture problems in viscoelastic media are also undertaken.The results obtained by the proposed method show good agreement with the analytical methods and provide reasonable background information to enhance the modeling of crack growth in viscoelastic media.展开更多
In this paper, we study Nitsche extended finite element method (XFEM) for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular famil...In this paper, we study Nitsche extended finite element method (XFEM) for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular family of grids and prove the optimal convergence rate of the scheme with respect to the mesh size. Main efforts are devoted onto classifying the cases of intersection between the elements and the interface and prove a weighted trace inequality for the extended finite element functions needed, and the general framework of analysing XFEM c^n be implemented then.展开更多
This paper analyzes two extended finite element methods(XFEMs)for linear quadratic optimal control problems governed by Poisson equation in non-convex domains.We follow the variational discretization concept to discre...This paper analyzes two extended finite element methods(XFEMs)for linear quadratic optimal control problems governed by Poisson equation in non-convex domains.We follow the variational discretization concept to discretize the continuous problems,and apply an XFEM with a cut-off function and a classic XFEM with a fixed enrichment area to discretize the state and co-state equations.Optimal error estimates are derived for the state,co-state and control.Numerical results confirm our theoretical results.展开更多
The extended finite element method(XFEM) is a numerical method for modeling discontinuities within a classical finite element framework. Based on the algorithm of XFEM, the major factors such as integral domain factor...The extended finite element method(XFEM) is a numerical method for modeling discontinuities within a classical finite element framework. Based on the algorithm of XFEM, the major factors such as integral domain factor and mesh density which all influence the calculation accuracy of stress intensity factor(SIF) are discussed,and the proper parameters to calculate the SIF are given. The results from the case analysis demonstrate that the crack path is the most sensitive to the crack growth increment size, and the crack path is not mesh-sensitive. A reanalysis method for the XFEM has been introduced. The example presented shows that there is a significantly reduced computational cost for each iteration of crack growth achieved by using the reanalysis method and the reanalysis approach has increasing benefits as the mesh density increases or the value of crack growth increments size decreases.展开更多
Fracking is one of the kernel technologies in the remarkable shale gas revolution. The extended finite element method is used in this paper to numerically investigate the interaction between hydraulic and natural frac...Fracking is one of the kernel technologies in the remarkable shale gas revolution. The extended finite element method is used in this paper to numerically investigate the interaction between hydraulic and natural fractures, which is an important issue of the enigmatic fracture network formation in fracking. The criteria which control the opening of natural fracture and crossing of hydraulic fracture are tentatively presented. Influence factors on the interaction process are systematically analyzed, which include the approach angle, anisotropy of in-situ stress and fluid pressure profile.展开更多
The finite element method (FEM) is one of the most popular and efficient methods for computational modeling in scien- tific research and engineering . To expand its application to more complex problems, lots of new ...The finite element method (FEM) is one of the most popular and efficient methods for computational modeling in scien- tific research and engineering . To expand its application to more complex problems, lots of new FEM-based methods have been developed in recent decades.展开更多
The extended finite element method is used to analyze a plate with two parallel edge cracks impacted by a cylindrical projectile. The influence of the impact speed, crack length, plate thickness and notch tip radius o...The extended finite element method is used to analyze a plate with two parallel edge cracks impacted by a cylindrical projectile. The influence of the impact speed, crack length, plate thickness and notch tip radius on the crack initiation and propagation is studied. Dynamics equations are solved by an implicit time integration scheme which is unconditionally stable. Very good agreement is achieved between numerical predictions and experimental results. The critical velocity of the crack initiation under different conditions is examined. The influence of the crack length is greater than that of the impact speed, plate thickness and notch tip radius.展开更多
In this paper,a concurrent simulation framework for fatigue crack growth analysis is proposed using a novel small time scale model for fatigue mechanism analysis and the extended finite element method(X-FEM)for fatigu...In this paper,a concurrent simulation framework for fatigue crack growth analysis is proposed using a novel small time scale model for fatigue mechanism analysis and the extended finite element method(X-FEM)for fatigue crack growth simulation.The proposed small time scale fatigue model does not require the cycle counting as those using the classical fatigue analysis methods and can be performed concurrently with structural/mechanical analysis.The X-FEM greatly facilitates crack growth simulation without remeshing requirements ahead of the crack tip as in the classical finite element method.The basic concept and theory of X-FEM was briefly introduced and numerical predictions of stress intensity factors are verified with reference solutions under both uniaxial and multiaxial loadings.The small time scale fatigue model is integrated into the numerical simulation algorithm for concurrent fatigue crack growth analysis.Model predictions are compared with available experimental observations for model validation.展开更多
An extended multiscale finite element method (EMsFEM) is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multi...An extended multiscale finite element method (EMsFEM) is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multiscale base functions to capture the small-scale features of the coarse elements in the multiscale finite element analysis.On the basis of our existing work for periodic truss materials, the construction methods of the base functions for continuum heterogeneous materials are systematically introduced. Numerical experiments show that the choice of boundary conditions for the construction of the base functions has a big influence on the accuracy of the multiscale solutions, thus,different kinds of boundary conditions are proposed. The efficiency and accuracy of the developed method are validated and the results with different boundary conditions are verified through extensive numerical examples with both periodic and random heterogeneous micro-structures.Also, a consistency test of the method is performed numerically. The results show that the EMsFEM can effectively obtain the macro response of the heterogeneous structures as well as the response in micro-scale,especially under the periodic boundary conditions.展开更多
Horizontal well-stimulation is the key to unconventional resource exploration and development.The development mode of the well plant helps increase the stimulated reservoir volume.Nevertheless,fracture interference be...Horizontal well-stimulation is the key to unconventional resource exploration and development.The development mode of the well plant helps increase the stimulated reservoir volume.Nevertheless,fracture interference between wells reduces the fracturing effect.Here,a 2D hydro-mechanical coupling model describing hydraulic fracture(HF)propagation is established with the extended finite element method,and the effects of several factors on HF propagation during multiple wells fracturing are analyzed.The results show that with an increase in elastic modulus,horizontal principal stress difference and injection fluid displacement,the total fracture area and the reservoir stimulation efficiency are both improved in all three fracturing technologies.After a comparison of the three technologies,the method of improved zipper fracturing is proposed,which avoids mutual interference between HFs,and the reservoir stimulation effect is improved significantly.The study provides guidance for optimizing the fracturing technology of multiple horizontal wells.展开更多
制造业的生产物流方式处于不断变革中,对其建模仿真可为制造系统规划设计、分析及改造提供决策支持。依“人-机-物-环-法”分类给出了智能车间制造系统中实体元素的描述,结合EFSM(extended finite state machine)和组件化建模思想,建立...制造业的生产物流方式处于不断变革中,对其建模仿真可为制造系统规划设计、分析及改造提供决策支持。依“人-机-物-环-法”分类给出了智能车间制造系统中实体元素的描述,结合EFSM(extended finite state machine)和组件化建模思想,建立了生产和物流组件化EFSM模型;阐述了智能车间多作业生产的建模过程以及组件模型实例化方法;通过EFSM-DEVS(discrete event system specification)模型自动转换及DEVS引擎完成了仿真运行。仿真结果表明:该方法所建立的模型更符合车间实际状况,适用性更广;组件化建模思想能构造更具扩展性的软件;建模及仿真运行的3D可视化使软件直观性更好,其仿真结果与AnyLogic保持一致。展开更多
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金support of Iran National Science Foundation is also gratefully appreciated
文摘High energy gas fracturing is a simple approach of applying high pressure gas to stimulate wells by gen- erating several radial cracks without creating any other damages to the wells. In this paper, a numerical algorithm is proposed to quantitatively simulate propagation of these fractures around a pressurized hole as a quasi-static phenomenon. The gas flow through the cracks is assumed as a one-dimensional transient flow, governed by equations of conservation of mass and momentum. The fractured medium is modeled with the extended finite element method, and the stress intensity factor is calculated by the simple, though sufficiently accurate, displacement ex- trapolation method. To evaluate the proposed algorithm, two field tests are simulated and the unknown parameters are determined through calibration. Sensitivity analyses are performed on the main effective parameters. Considering that the level of uncertainty is very high in these types of engineering problems, the results show a good agreement with the experimental data. They are also consistent with the theory that the final crack length is mainly determined by the gas pressure rather than the initial crack length produced by the stress waves.
基金supported by the Aviation Science Foundation of China
文摘The reliability of real-time embedded software directly determines the reliability of the whole real-time embedded sys- tem, and the effective software testing is an important way to ensure software quality and reliability. Based on the analysis of the characteristics of real-time embedded software, the formal method is introduced into the real-time embedded software testing field and the real-time extended finite state machine (RT-EFSM) model is studied firstly. Then, the time zone division method of real-time embedded system is presented and the definition and description methods of time-constrained transition equivalence class (timeCTEC) are presented. Furthermore, the approaches of the testing sequence and test case generation are put forward. Finally, the proposed method is applied to a typical avionics real- time embedded software testing practice and the examples of the timeCTEC, testing sequences and test cases are given. With the analysis of the testing result, the application verification shows that the proposed method can effectively describe the real-time embedded software state transition characteristics and real-time requirements and play the advantages of the formal methods in accuracy, effectiveness and the automation supporting. Combined with the testing platform, the real-time, closed loop and automated simulation testing for real-time embedded software can be realized effectively.
基金Project(50578038)supported by the National Natural Science Foundation of China
文摘Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Firstly,the single-edge notched beam test was used to analyze the temperature effect and calculate the material parameters.Then,the mechanical responses were studied using numerical analysis.It is concluded that 5℃ can be selected as the critical temperature that affects the fracture properties,and numerical simulations indicate that crack propagation is found to significantly affect the stress state of the epoxy asphalt mixture.The maximum principal stress at the crack surface exhibits different trends at various temperatures.Numerical solution of stress intensity factor can well meet the theoretical solution,especially when the temperature is lower than 5℃.
基金Projects(41172244,41072224) supported by the National Natural Science Foundation of ChinaProject(2009GGJS-037) supported by the Foundation of Youths Key Teacher by the Henan Educational Committee,China
文摘Extended finite element method (XFEM) implementation of the interaction integral methodology for evaluating the stress intensity factors (SIF) of the mixed-mode crack problem is presented. A discontinuous function and the near-tip asymptotic function are added to the classic finite element approximation to model the crack behavior. Two-state integral by the superposition of actual and auxiliary fields is derived to calculate the SIFs. Applications of the proposed technique to the inclined centre crack plate with inclined angle from 0° to 90° and slant edge crack plate with slant angle 45°, 67.5° and 90° are presented, and comparisons are made with closed form solutions. The results show that the proposed method is convenient, accurate and computationallv efficient.
基金supported by the National Natural Science Foundation of China (Grants 11471262, 50976003, 51136005)
文摘In this paper, the extended finite element method (XFEM) is adopted to analyze the interaction between a single macroscopic inclusion and a single macroscopic crack as well as that between multiple macroscopic or microscopic defects under thermal/mechanical load. The effects of different shapes of multiple inclusions on the material thermomechanical response are investigated, and the level set method is coupled with XFEM to analyze the interaction of multiple defects. Further, the discretized extended finite element approximations in relation to thermoelastic problems of multiple defects under displacement or temperature field are given. Also, the interfaces of cracks or materials are represented by level set functions, which allow the mesh assignment not to conform to crack or material interfaces. Moreover, stress intensity factors of cracks are obtained by the interaction integral method or the M-integral method, and the stress/strain/stiffness fields are simulated in the case of multiple cracks or multiple inclusions. Finally, some numerical examples are provided to demonstrate the accuracy of our proposed method.
基金Project(2017YFC0404802)supported by the National Key R&D Program of ChinaProjects(U1965206,51979143)supported by the National Natural Science Foundation of China。
文摘In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.
基金Funded by the National Natural Science Foundation of China(Nos.51409264,51509020,51209219)Zhejiang Provincial Natural Science Foundation(No.LQ15E090003)
文摘Self-deformation cracking is the cracking caused by thermal deformation, autogenous volume deformation or shrinkage deformation. In this paper, an extended finite element calculation method was deduced for concrete crack propagation under a constant hydration and hardening condition during the construction period, and a corresponding programming code was developed. The experimental investigation shows that initial crack propagation caused by self-deformation loads can be analyzed by this program. This improved algorithm was a preliminary application of the XFEM to the problem of the concrete self-deformation cracking during the hydration and hardening period. However, room for improvement exists for this algorithm in terms of matching calculation programs with mass concrete temperature fields containing cooling pipes and the influence of creep or damage on crack propagation.
基金supported by the National Basic Research Program of China(2014CB046904)the Hubei Provincial Key Laboratory of Safety for Geotechnical and Structural Engineering at Wuhan University(HBKLCIV201207)the China Postdoctoral Science Foundation(2013M540604)
文摘This study develops a mesoscopic framework and methodology for the modeling of microcracks in concrete. A new algorithm is first proposed for the generation of random concrete meso-structure including microcracks and then coupled with the extended finite element method to simulate the heterogeneities and discontinuities present in the meso-structure of concrete. The proposed procedure is verified and exemplified by a series of numerical simulations. The simulation results show that microcracks can exert considerable impact on the fracture performance of concrete. More broadly, this work provides valuable insight into the initiation and propagation mechanism of microcracks in concrete and helps to foster a better understanding of the micro-mechanical behavior of cementitious materials.
基金supported by the National Basic Research Program of China ("973" Program) (Grant No 2007CB714104)
文摘In this paper,the problem of modeling crack in 2D viscoelastic media is studied using the extended finite element method.The paper focuses on the definition of enrichment functions suitable for cracks assessment in viscoelastic media and the generalized domain integrals used in the determination of crack tip parameters.The opening mode and mixed mode solutions of crack tip fracture problems in viscoelastic media are also undertaken.The results obtained by the proposed method show good agreement with the analytical methods and provide reasonable background information to enhance the modeling of crack growth in viscoelastic media.
基金The first author is partially supported by the U.S. Department of Energy, Office of Science, Office of Advanced Scientific Computing Research as part of the Collaboratory on Mathematics for Mesoscopic Modeling of Materials under Award Number DE-SC-0009249, and the Key Program of National Natural Science Foundation of China with Grant No. 91430215. The second author is supported by State Key Laboratory of Scientific and Engineering Computing (LSEC), National Center for Mathematics and Interdisciplinary Sciences of Chinese Academy of Sciences (NCMIS), and National Natural Science Foundation of China with Grant No. 11471026 he is thankful to the Center for Computational Mathematics and Applications, the Pennsylvania State University, where he worked on this manuscript as a visiting scholar. The authors are grateful to Professor Jinchao Xu, Dr. Yuanming Xiao and Dr. Maximilian Metti for their valuable suggestions and discussions, to Professor Haijun Wu for his valuable help on preparing the numerical example, and to the anonymous referee for the valuable comments and suggestion which lead to improvements of the paper.
文摘In this paper, we study Nitsche extended finite element method (XFEM) for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular family of grids and prove the optimal convergence rate of the scheme with respect to the mesh size. Main efforts are devoted onto classifying the cases of intersection between the elements and the interface and prove a weighted trace inequality for the extended finite element functions needed, and the general framework of analysing XFEM c^n be implemented then.
基金supported by National Natural Science Foundation of China(Grant No.11771312)。
文摘This paper analyzes two extended finite element methods(XFEMs)for linear quadratic optimal control problems governed by Poisson equation in non-convex domains.We follow the variational discretization concept to discretize the continuous problems,and apply an XFEM with a cut-off function and a classic XFEM with a fixed enrichment area to discretize the state and co-state equations.Optimal error estimates are derived for the state,co-state and control.Numerical results confirm our theoretical results.
基金the National Basic Research Program(973) of China(No.2011CB013505)the National Natural Science Foundation of China(No.51279100)
文摘The extended finite element method(XFEM) is a numerical method for modeling discontinuities within a classical finite element framework. Based on the algorithm of XFEM, the major factors such as integral domain factor and mesh density which all influence the calculation accuracy of stress intensity factor(SIF) are discussed,and the proper parameters to calculate the SIF are given. The results from the case analysis demonstrate that the crack path is the most sensitive to the crack growth increment size, and the crack path is not mesh-sensitive. A reanalysis method for the XFEM has been introduced. The example presented shows that there is a significantly reduced computational cost for each iteration of crack growth achieved by using the reanalysis method and the reanalysis approach has increasing benefits as the mesh density increases or the value of crack growth increments size decreases.
基金supported by the National Natural Science Foundation of China (Grant No. 11372157)the Special Research Grant for Doctor Discipline by Ministry of Education of China (Grant No. 20120002110075)the Foundation for the Author of National Excellent Doctoral Dissertation of China (FANEDD) (Grant No. 201326)
文摘Fracking is one of the kernel technologies in the remarkable shale gas revolution. The extended finite element method is used in this paper to numerically investigate the interaction between hydraulic and natural fractures, which is an important issue of the enigmatic fracture network formation in fracking. The criteria which control the opening of natural fracture and crossing of hydraulic fracture are tentatively presented. Influence factors on the interaction process are systematically analyzed, which include the approach angle, anisotropy of in-situ stress and fluid pressure profile.
文摘The finite element method (FEM) is one of the most popular and efficient methods for computational modeling in scien- tific research and engineering . To expand its application to more complex problems, lots of new FEM-based methods have been developed in recent decades.
基金supported by the National Natural Science Foundation of China(Nos.11272096 and 11472086)the Research Fund for the Doctoral Program of Higher Education of China(No.20112304110015)
文摘The extended finite element method is used to analyze a plate with two parallel edge cracks impacted by a cylindrical projectile. The influence of the impact speed, crack length, plate thickness and notch tip radius on the crack initiation and propagation is studied. Dynamics equations are solved by an implicit time integration scheme which is unconditionally stable. Very good agreement is achieved between numerical predictions and experimental results. The critical velocity of the crack initiation under different conditions is examined. The influence of the crack length is greater than that of the impact speed, plate thickness and notch tip radius.
基金support is greatly appreciated.The research reported in this paper was supported by funds from NSF(No.CMMI-0900111,Project Manager:Dr.Mahendra Singh).
文摘In this paper,a concurrent simulation framework for fatigue crack growth analysis is proposed using a novel small time scale model for fatigue mechanism analysis and the extended finite element method(X-FEM)for fatigue crack growth simulation.The proposed small time scale fatigue model does not require the cycle counting as those using the classical fatigue analysis methods and can be performed concurrently with structural/mechanical analysis.The X-FEM greatly facilitates crack growth simulation without remeshing requirements ahead of the crack tip as in the classical finite element method.The basic concept and theory of X-FEM was briefly introduced and numerical predictions of stress intensity factors are verified with reference solutions under both uniaxial and multiaxial loadings.The small time scale fatigue model is integrated into the numerical simulation algorithm for concurrent fatigue crack growth analysis.Model predictions are compared with available experimental observations for model validation.
基金supported by the National Natural Science Foundation(10721062,11072051,90715037,10728205,91015003, 51021140004)the Program of Introducing Talents of Discipline to Universities(B08014)the National Key Basic Research Special Foundation of China(2010CB832704).
文摘An extended multiscale finite element method (EMsFEM) is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multiscale base functions to capture the small-scale features of the coarse elements in the multiscale finite element analysis.On the basis of our existing work for periodic truss materials, the construction methods of the base functions for continuum heterogeneous materials are systematically introduced. Numerical experiments show that the choice of boundary conditions for the construction of the base functions has a big influence on the accuracy of the multiscale solutions, thus,different kinds of boundary conditions are proposed. The efficiency and accuracy of the developed method are validated and the results with different boundary conditions are verified through extensive numerical examples with both periodic and random heterogeneous micro-structures.Also, a consistency test of the method is performed numerically. The results show that the EMsFEM can effectively obtain the macro response of the heterogeneous structures as well as the response in micro-scale,especially under the periodic boundary conditions.
基金funded by Shaanxi Natural Science Basic Research Program Project Study on Liquid Propellant High Energy Gas Fracturing Mechanism in Radial Well Based on Phase Field Method(No.2019JQ-824)NSFC Projects Evolution Mechanism and Effectiveness Evaluation of Fracture Network Produced by Volume Fracturing with Tighter Clusters in Continental Shale Oil Reservoir(No.52274040)+1 种基金Study on Thermal Secondary Pore Evolution and Salt Precipitation Regulation Mechanism in Fire Flooding Reservoirs Based on Multi-field Coupling of Thermal-Flow-Solid-Chemical(No.52274039)Xi’an Shiyou University Youth Scientific Research and Innovation Team Operation Funds in 2018 Flow Mechanism of Complex Reservoirs and High Efficiency Development and Oil Production Technology(No.115080020).
文摘Horizontal well-stimulation is the key to unconventional resource exploration and development.The development mode of the well plant helps increase the stimulated reservoir volume.Nevertheless,fracture interference between wells reduces the fracturing effect.Here,a 2D hydro-mechanical coupling model describing hydraulic fracture(HF)propagation is established with the extended finite element method,and the effects of several factors on HF propagation during multiple wells fracturing are analyzed.The results show that with an increase in elastic modulus,horizontal principal stress difference and injection fluid displacement,the total fracture area and the reservoir stimulation efficiency are both improved in all three fracturing technologies.After a comparison of the three technologies,the method of improved zipper fracturing is proposed,which avoids mutual interference between HFs,and the reservoir stimulation effect is improved significantly.The study provides guidance for optimizing the fracturing technology of multiple horizontal wells.
文摘制造业的生产物流方式处于不断变革中,对其建模仿真可为制造系统规划设计、分析及改造提供决策支持。依“人-机-物-环-法”分类给出了智能车间制造系统中实体元素的描述,结合EFSM(extended finite state machine)和组件化建模思想,建立了生产和物流组件化EFSM模型;阐述了智能车间多作业生产的建模过程以及组件模型实例化方法;通过EFSM-DEVS(discrete event system specification)模型自动转换及DEVS引擎完成了仿真运行。仿真结果表明:该方法所建立的模型更符合车间实际状况,适用性更广;组件化建模思想能构造更具扩展性的软件;建模及仿真运行的3D可视化使软件直观性更好,其仿真结果与AnyLogic保持一致。