In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolso...In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.展开更多
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.展开更多
In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
On the basis of Hartmann Shack sensor imaging analysis, a new method is presented with which the wavefront slope can be obtained when the object is incoherent and extended. This method, which is demonstrated by both ...On the basis of Hartmann Shack sensor imaging analysis, a new method is presented with which the wavefront slope can be obtained when the object is incoherent and extended. This method, which is demonstrated by both theoretical interpreting and computer simulation, explains how to measure the wavefront slope difference between two sub apertures through the determination of image displacements on detector plane. It includes a fast and accurate digital algorithm for detecting wavefront disturbance, which is much suitable for realization in such electrical hardwares as digital signal processors.展开更多
This paper presents the flight dynamical behavior of the thrust vectoring aircraft with extended bifurcation and continuation methods. In contrast to the standard bifurcation and continuation methods, the extended met...This paper presents the flight dynamical behavior of the thrust vectoring aircraft with extended bifurcation and continuation methods. In contrast to the standard bifurcation and continuation methods, the extended methods are capable of calculating the continuation curves of the equilibrium points for the particular type of trimming flight. Therefore, these methods can not only give the performance measures of aircraft, but also determine the stability of trimming points. In this paper, the methods are used to verify the effectiveness of the thrust vectoring control law, to define the flight envelope boundary, to analyze the stability and controllability of trimming flight, and to predict the departures of the instable flight. The result shows that the extended methods provide more flight dynamic information and are useful in preliminary design of the thrust vectoring aircraft.展开更多
This paper applies an extended auxiliary equation method to obtain exact solutions of the KdV equation with variable coefficients. As a result, solitary wave solutions, trigonometric function solutions, rational funct...This paper applies an extended auxiliary equation method to obtain exact solutions of the KdV equation with variable coefficients. As a result, solitary wave solutions, trigonometric function solutions, rational function solutions, Jacobi elliptic doubly periodic wave solutions, and nonsymmetrical kink solution are obtained. It is shown that the extended auxiliary equation method, with the help of a computer symbolic computation system, is reliable and effective in finding exact solutions of variable coefficient nonlinear evolution equations in mathematical physics.展开更多
The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that ...The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that remeshing for moving discontinuities can be overcome. The extended finite element method is presented for hydro-mechanical modeling of impermeable discontinuities in rock. The governing equation of XFEM for hydraulic fracture modeling is derived by the virtual work principle of the fracture problem considering the water pressure on crack surface. The coupling relationship between water pressure gradient on crack surface and fracture opening width is obtained by semi-analytical and semi-numerical method. This method simplifies coupling analysis iteration and improves computational precision. Finally, the efficiency of the proposed method for modeling hydraulic fracture problems is verified by two examples and the advantages of the XFEM for hydraulic fracturing analysis are displayed.展开更多
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.展开更多
An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primar...An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primary group with an average geometrical size larger than or in the same order of magnitude of wavelength of a concerned stress wave is defined as 'macro-joints',while the secondary group with a high density and relatively small geometrical size compared to the wavelength is known as 'micro-defects'.The rock mass with micro-defects is modeled as an equivalent viscoelastic medium while the macro-joints in the rock mass are modeled explicitly as physical discontinuities.Viscoelastic properties of a micro-defected sedimentary rock are obtained by longitudinally impacting a cored long sedimentary rod with a pendulum.Wave propagation coefficient and dynamic viscoelastic modulus are measured.The EDDM is then successfully employed to analyze the wave propagation across macro-joint in VRM.The effect of the rock viscosity on the stress wave propagation is evaluated by comparing the results of VRM from the presented EDDM with those of an elastic rock mass (ERM) from the conventional displacement discontinuity method (CDDM).The CDDM is a special case of the EDDM under the condition that the rock viscosity is ignored.Comparison of the reflected and transmitted waves shows that the essential rock viscosity has a significant effect on stress wave attenuation.When a short propagation distance of a stress wave is considered,the results obtained from the CDDM approximate to the EDDM solutions,however,when the propagation distance is sufficiently long relative to the wavelength,the effect of rock viscosity on the stress wave propagation cannot be ignored.展开更多
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.展开更多
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.展开更多
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℃.展开更多
A stuck drill string results in a major non-productive cost in extended reach drilling engineering. The first step is to determine the depth at which the sticking has occurred. Methods of measurement have been proved ...A stuck drill string results in a major non-productive cost in extended reach drilling engineering. The first step is to determine the depth at which the sticking has occurred. Methods of measurement have been proved useful for determining the stuck points, but these operations take considerable time. As a result of the limitation with the current operational practices, calculation methods are still preferred to estimate the stuck point depth. Current analytical methods do not consider friction and are only valid for vertical rather than extended reach wells. The numerical method is established to take full account of down hole friction, tool joint, upset end of drill pipe, combination drill strings and tubular materials so that it is valid to determine the stuck point in extended reach wells. The pull test, torsion test and combined test of rotation and pulling can be used to determine the stuck point. The results show that down hole friction, tool joint, upset end of drill pipe, tubular sizes and materials have significant effects on the pull length and/or the twist angle of the stuck drill string.展开更多
Drilling engineering has great uncertainty and it always involves huge investment and high risk. Risk analysis of extended reach drilling (ERD) is very important to prevent complex failures and to improve drilling e...Drilling engineering has great uncertainty and it always involves huge investment and high risk. Risk analysis of extended reach drilling (ERD) is very important to prevent complex failures and to improve drilling efficiency. Nowadays there are few reports on how to analyze quantitatively the drilling risk for extended reach wells (ERWs). Based on the fuzzy set theory, a comprehensive fuzzy evaluation model for analyzing risks of ERD is proposed in this paper. Well B6ERW07 is a planned 8,000-meter ERW with a high ratio of horizontal displacement (HD) to vertical depth (VD) in the Liuhua Oilfield, the South China Sea, China. On the basis of the evaluation model developed in this study, the risk for drilling Well B6ERW07 was evaluated before drilling. The evaluation result shows that the success rate of drilling this well is predicted to be 51.9%, providing important rational and scientific information for the decisionmakers.展开更多
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 this paper, a new approach of maneuvering target tracking algorithm based on the autoregressive extended Viterbi(AREV) model is proposed. In contrast to weakness of traditional constant velocity(CV) and constant ac...In this paper, a new approach of maneuvering target tracking algorithm based on the autoregressive extended Viterbi(AREV) model is proposed. In contrast to weakness of traditional constant velocity(CV) and constant acceleration(CA) models to noise effect reduction, the autoregressive(AR) part of the new model which changes the structure of state space equations is proposed. Also using a dynamic form of the state transition matrix leads to improving the rate of convergence and decreasing the noise effects. Since AR will impose the load of overmodeling to the computations, the extended Viterbi(EV) method is incorporated to AR in two cases of EV1 and EV2. According to most probable paths in the interacting multiple model(IMM) during nonmaneuvering and maneuvering parts of estimation, EV1 and EV2 respectively can decrease load of overmodeling computations and improve the AR performance. This new method is coupled with proposed detection schemes for maneuver occurrence and termination as well as for switching initializations. Appropriate design parameter values are derived for the detection schemes of maneuver occurrences and terminations. Finally, simulations demonstrate that the performance of the proposed model is better than the other older linear and also nonlinear algorithms in constant velocity motions and also in various types of maneuvers.展开更多
Taking the Konopelchenko-Dubrovsky system as a simple example, some familles of rational formal hyperbolic function solutions, rational formal triangular periodic solutions, and rational solutions are constructed by u...Taking the Konopelchenko-Dubrovsky system as a simple example, some familles of rational formal hyperbolic function solutions, rational formal triangular periodic solutions, and rational solutions are constructed by using the extended Riccati equation rational expansion method presented by us. The method can also be applied to solve more nonlinear partial differential equation or equations.展开更多
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.展开更多
For the appearance of the additive perturbation of controller gain when the controller parameter has minute adjustment at the initial running stage of system,to avoid the adverse effects,this paper investigates the mi...For the appearance of the additive perturbation of controller gain when the controller parameter has minute adjustment at the initial running stage of system,to avoid the adverse effects,this paper investigates the mixed H_2/H_∞ state feedback attitude control problem of microsatellite based on extended LMI method.Firstly,the microsatellite attitude control system is established and transformed into corresponding state space form.Then,without the equivalence restriction of the two Lyapunov variables of H_2 and H∞performance,this paper introduces additional variables to design the mixed H_2/H_∞ control method based on LMI which can also reduce the conservatives.Finally,numerical simulations are analyzed to show that the proposed method can make the satellite stable within 20 s whether there is additive perturbation of the controller gain or not.The comparative analysis of the simulation results between extended LMI method and traditional LMI method also demonstrates the effectiveness and feasibility of the proposed method in this paper.展开更多
基金supported by the Key Laboratory of Road Construction Technology and Equipment(Chang’an University,No.300102253502)the Natural Science Foundation of Shandong Province of China(GrantNo.ZR2022YQ06)the Development Plan of Youth Innovation Team in Colleges and Universities of Shandong Province(Grant No.2022KJ140).
文摘In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.
文摘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.
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.
文摘On the basis of Hartmann Shack sensor imaging analysis, a new method is presented with which the wavefront slope can be obtained when the object is incoherent and extended. This method, which is demonstrated by both theoretical interpreting and computer simulation, explains how to measure the wavefront slope difference between two sub apertures through the determination of image displacements on detector plane. It includes a fast and accurate digital algorithm for detecting wavefront disturbance, which is much suitable for realization in such electrical hardwares as digital signal processors.
文摘This paper presents the flight dynamical behavior of the thrust vectoring aircraft with extended bifurcation and continuation methods. In contrast to the standard bifurcation and continuation methods, the extended methods are capable of calculating the continuation curves of the equilibrium points for the particular type of trimming flight. Therefore, these methods can not only give the performance measures of aircraft, but also determine the stability of trimming points. In this paper, the methods are used to verify the effectiveness of the thrust vectoring control law, to define the flight envelope boundary, to analyze the stability and controllability of trimming flight, and to predict the departures of the instable flight. The result shows that the extended methods provide more flight dynamic information and are useful in preliminary design of the thrust vectoring aircraft.
基金Project supported by the Fundamental Research Funds for the Central Universities (Grant No. 2010B17914) and the National Natural Science Foundation of China (Grant No. 10926162).
文摘This paper applies an extended auxiliary equation method to obtain exact solutions of the KdV equation with variable coefficients. As a result, solitary wave solutions, trigonometric function solutions, rational function solutions, Jacobi elliptic doubly periodic wave solutions, and nonsymmetrical kink solution are obtained. It is shown that the extended auxiliary equation method, with the help of a computer symbolic computation system, is reliable and effective in finding exact solutions of variable coefficient nonlinear evolution equations in mathematical physics.
基金Project(2011CB013505)supported by the National Basic Research Program of ChinaProject(51279100)supported by the National Natural Science Foundation of China
文摘The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that remeshing for moving discontinuities can be overcome. The extended finite element method is presented for hydro-mechanical modeling of impermeable discontinuities in rock. The governing equation of XFEM for hydraulic fracture modeling is derived by the virtual work principle of the fracture problem considering the water pressure on crack surface. The coupling relationship between water pressure gradient on crack surface and fracture opening width is obtained by semi-analytical and semi-numerical method. This method simplifies coupling analysis iteration and improves computational precision. Finally, the efficiency of the proposed method for modeling hydraulic fracture problems is verified by two examples and the advantages of the XFEM for hydraulic fracturing analysis are displayed.
基金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.
文摘An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primary group with an average geometrical size larger than or in the same order of magnitude of wavelength of a concerned stress wave is defined as 'macro-joints',while the secondary group with a high density and relatively small geometrical size compared to the wavelength is known as 'micro-defects'.The rock mass with micro-defects is modeled as an equivalent viscoelastic medium while the macro-joints in the rock mass are modeled explicitly as physical discontinuities.Viscoelastic properties of a micro-defected sedimentary rock are obtained by longitudinally impacting a cored long sedimentary rod with a pendulum.Wave propagation coefficient and dynamic viscoelastic modulus are measured.The EDDM is then successfully employed to analyze the wave propagation across macro-joint in VRM.The effect of the rock viscosity on the stress wave propagation is evaluated by comparing the results of VRM from the presented EDDM with those of an elastic rock mass (ERM) from the conventional displacement discontinuity method (CDDM).The CDDM is a special case of the EDDM under the condition that the rock viscosity is ignored.Comparison of the reflected and transmitted waves shows that the essential rock viscosity has a significant effect on stress wave attenuation.When a short propagation distance of a stress wave is considered,the results obtained from the CDDM approximate to the EDDM solutions,however,when the propagation distance is sufficiently long relative to the wavelength,the effect of rock viscosity on the stress wave propagation cannot be ignored.
基金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.
基金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.
基金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℃.
基金support from the national projects(Grant No.:2011ZX05009-005and2010CB226703)
文摘A stuck drill string results in a major non-productive cost in extended reach drilling engineering. The first step is to determine the depth at which the sticking has occurred. Methods of measurement have been proved useful for determining the stuck points, but these operations take considerable time. As a result of the limitation with the current operational practices, calculation methods are still preferred to estimate the stuck point depth. Current analytical methods do not consider friction and are only valid for vertical rather than extended reach wells. The numerical method is established to take full account of down hole friction, tool joint, upset end of drill pipe, combination drill strings and tubular materials so that it is valid to determine the stuck point in extended reach wells. The pull test, torsion test and combined test of rotation and pulling can be used to determine the stuck point. The results show that down hole friction, tool joint, upset end of drill pipe, tubular sizes and materials have significant effects on the pull length and/or the twist angle of the stuck drill string.
基金support from the project of CNOOC China Limited-Shenzhen (Grant No. Z2007SLSZ-034)the foundation project of the State Key Laboratory of Petroleum Resource and Prospecting (Grant No. PRPDX2008-08) is gratefully acknowledged
文摘Drilling engineering has great uncertainty and it always involves huge investment and high risk. Risk analysis of extended reach drilling (ERD) is very important to prevent complex failures and to improve drilling efficiency. Nowadays there are few reports on how to analyze quantitatively the drilling risk for extended reach wells (ERWs). Based on the fuzzy set theory, a comprehensive fuzzy evaluation model for analyzing risks of ERD is proposed in this paper. Well B6ERW07 is a planned 8,000-meter ERW with a high ratio of horizontal displacement (HD) to vertical depth (VD) in the Liuhua Oilfield, the South China Sea, China. On the basis of the evaluation model developed in this study, the risk for drilling Well B6ERW07 was evaluated before drilling. The evaluation result shows that the success rate of drilling this well is predicted to be 51.9%, providing important rational and scientific information for the decisionmakers.
基金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.
文摘In this paper, a new approach of maneuvering target tracking algorithm based on the autoregressive extended Viterbi(AREV) model is proposed. In contrast to weakness of traditional constant velocity(CV) and constant acceleration(CA) models to noise effect reduction, the autoregressive(AR) part of the new model which changes the structure of state space equations is proposed. Also using a dynamic form of the state transition matrix leads to improving the rate of convergence and decreasing the noise effects. Since AR will impose the load of overmodeling to the computations, the extended Viterbi(EV) method is incorporated to AR in two cases of EV1 and EV2. According to most probable paths in the interacting multiple model(IMM) during nonmaneuvering and maneuvering parts of estimation, EV1 and EV2 respectively can decrease load of overmodeling computations and improve the AR performance. This new method is coupled with proposed detection schemes for maneuver occurrence and termination as well as for switching initializations. Appropriate design parameter values are derived for the detection schemes of maneuver occurrences and terminations. Finally, simulations demonstrate that the performance of the proposed model is better than the other older linear and also nonlinear algorithms in constant velocity motions and also in various types of maneuvers.
基金The project partially supported by the State Key Basic Research Program of China under Grant No. 2004CB318000
文摘Taking the Konopelchenko-Dubrovsky system as a simple example, some familles of rational formal hyperbolic function solutions, rational formal triangular periodic solutions, and rational solutions are constructed by using the extended Riccati equation rational expansion method presented by us. The method can also be applied to solve more nonlinear partial differential equation or equations.
基金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.
文摘For the appearance of the additive perturbation of controller gain when the controller parameter has minute adjustment at the initial running stage of system,to avoid the adverse effects,this paper investigates the mixed H_2/H_∞ state feedback attitude control problem of microsatellite based on extended LMI method.Firstly,the microsatellite attitude control system is established and transformed into corresponding state space form.Then,without the equivalence restriction of the two Lyapunov variables of H_2 and H∞performance,this paper introduces additional variables to design the mixed H_2/H_∞ control method based on LMI which can also reduce the conservatives.Finally,numerical simulations are analyzed to show that the proposed method can make the satellite stable within 20 s whether there is additive perturbation of the controller gain or not.The comparative analysis of the simulation results between extended LMI method and traditional LMI method also demonstrates the effectiveness and feasibility of the proposed method in this paper.