This paper presents the dimension split element-free Galerkin (DSEFG) method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-d...This paper presents the dimension split element-free Galerkin (DSEFG) method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares (IMLS) approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.展开更多
Splitting modulus variational principle in linear theory of solid mechanics was introduced, the principle for thin plate was derived, and splitting modulus finite element method of thin plate was established too. The ...Splitting modulus variational principle in linear theory of solid mechanics was introduced, the principle for thin plate was derived, and splitting modulus finite element method of thin plate was established too. The distinctive feature of the splitting model is that its functional contains one or more arbitrary additional parameters, called splitting factors, so stiffness of the model can be adjusted by properly selecting the splitting factors. Examples show that splitting modulus method has high precision and the ability to conquer some ill-conditioned problems in usual finite elements. The cause why the new method could transform the ill-conditioned problems into well-conditioned problem, is analyzed finally.展开更多
In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split i...In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.展开更多
The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain (CBTD) method. A governing equation in the cylindrical coordinate system is formulat...The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain (CBTD) method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,展开更多
A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal...A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.展开更多
In this paper,we consider the initial-boundary value problem(IBVP)for the micropolar Naviers-Stokes equations(MNSE)and analyze a first order fully discrete mixed finite element scheme.We first establish some regularit...In this paper,we consider the initial-boundary value problem(IBVP)for the micropolar Naviers-Stokes equations(MNSE)and analyze a first order fully discrete mixed finite element scheme.We first establish some regularity results for the solution of MNSE,which seem to be not available in the literature.Next,we study a semi-implicit time-discrete scheme for the MNSE and prove L2-H1 error estimates for the time discrete solution.Furthermore,certain regularity results for the time discrete solution are establishes rigorously.Based on these regularity results,we prove the unconditional L2-H1 error estimates for the finite element solution of MNSE.Finally,some numerical examples are carried out to demonstrate both accuracy and efficiency of the fully discrete finite element scheme.展开更多
Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to e...Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.展开更多
In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bil...In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ||·||div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.展开更多
A splitting positive definite mixed finite element method is proposed for second-order viscoelasticity wave equation. The proposed procedure can be split into three independent symmetric positive definite integro-diff...A splitting positive definite mixed finite element method is proposed for second-order viscoelasticity wave equation. The proposed procedure can be split into three independent symmetric positive definite integro-differential sub-system and does not need to solve a coupled system of equations. Error estimates are derived for both semidiscrete and fully discrete schemes. The existence and uniqueness for semidiscrete scheme are proved. Finally, a numerical example is provided to illustrate the efficiency of the method.展开更多
By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IE...By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IEFG) method is used for 2D transient heat conduction problems, and the finite difference method is applied in the splitting direction. The discretized system equation is obtained based on the Galerkin weak form of 2D problem; the essential boundary conditions are imposed with the penalty method; and the finite difference method is employed in the time domain. Four exemplary problems are chosen to verify the efficiency of the DSEFG method. The numerical solutions show that the efficiency and precision of the DSEFG method are greater than ones of the IEFG method for 3D problems.展开更多
Split Hopkinson pressure bar(SHPB) apparatus, usually used for testing behavior of material in median and high strain-rate, is now widely used in the study of rock dynamic constitutive relation, damage evolvement mech...Split Hopkinson pressure bar(SHPB) apparatus, usually used for testing behavior of material in median and high strain-rate, is now widely used in the study of rock dynamic constitutive relation, damage evolvement mechanism and energy consumption. However, the possible reasons of sampling disturbance, machining error and so on often lead to the scattering of test results, and bring ultimate difficulty for forming general test conclusion. Based on the stochastic finite element method, the uncertain parameters of specimen density ρs, specimen radius Rs, specimen elastic modulus Es and specimen length Ls in the data processing of SHPB test were considered, and the correlation between the parameters and the test results was analyzed. The results show that the specimen radius Rs has direct correlation with the test result, improving the accuracy in preparing and measuring of specimen is an effective way to improve the accuracy of test and minish the scattering of results for SHPB test.展开更多
Based on the thermal kinetic and mass conservation, a series of controlling equations for the finite element are derived and related programs are developed to simulate the damage microcrack healing process controlled ...Based on the thermal kinetic and mass conservation, a series of controlling equations for the finite element are derived and related programs are developed to simulate the damage microcrack healing process controlled by surface diffusion. Two kinds of typical models for microcrack splitting are proposed, i.e., the grain boundary energy existing on the crack surface and residual stresses applying on the crack surface. And the conditions of microcrack splitting in the two models are given as a function of the microcrack aspect ratio. The microcrack with traction-free surfaces will directly evolve into a spheroid.展开更多
In this study, the 3-dimensional discrete element method is firstly introduced to explain the fracturing damage process of the dynamic split experiment of a special brittle glass ZnS. The corresponding dynamic split e...In this study, the 3-dimensional discrete element method is firstly introduced to explain the fracturing damage process of the dynamic split experiment of a special brittle glass ZnS. The corresponding dynamic split experiment is also performed by using the split Hopkinson pressure bar. Then the numerical results correspond closely to those obtained by experiments, and the fracturing damage mode shows that the sample under high strain rate loading would crack along vertical diameter in the band region between two loading edges, which differs from the static damage mode. Furthermore, by comparing a group of contrast numerical tests, the numerical results prove that loading area upon the top side of samples would influence the fracture mode of dynamic split experiments, which indicates that the narrow loading plane is better.展开更多
基金supported by the National Natural Science Foundation of China (Grants 11571223, 51404160)Shanxi Province Science Foundation for Youths (Grant 2014021025-1)
文摘This paper presents the dimension split element-free Galerkin (DSEFG) method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares (IMLS) approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
文摘Splitting modulus variational principle in linear theory of solid mechanics was introduced, the principle for thin plate was derived, and splitting modulus finite element method of thin plate was established too. The distinctive feature of the splitting model is that its functional contains one or more arbitrary additional parameters, called splitting factors, so stiffness of the model can be adjusted by properly selecting the splitting factors. Examples show that splitting modulus method has high precision and the ability to conquer some ill-conditioned problems in usual finite elements. The cause why the new method could transform the ill-conditioned problems into well-conditioned problem, is analyzed finally.
基金financially supported by the National Natural Science Foundation of China(Grant No.51349011)the Foundation of Si’chuan Educational Committee(Grant No.17ZB0452)+1 种基金the Innovation Team Project of Si’chuan Educational Committee(Grant No.18TD0019)the Longshan Academic Talent Research Support Program of the Southwest of Science and Technology(Grant Nos.18LZX715 and 18LZX410)
文摘In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.
文摘The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain (CBTD) method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,
文摘A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.
基金supported by the National Natural Science Foundation of China(Grant Nos.11871467,11471329).
文摘In this paper,we consider the initial-boundary value problem(IBVP)for the micropolar Naviers-Stokes equations(MNSE)and analyze a first order fully discrete mixed finite element scheme.We first establish some regularity results for the solution of MNSE,which seem to be not available in the literature.Next,we study a semi-implicit time-discrete scheme for the MNSE and prove L2-H1 error estimates for the time discrete solution.Furthermore,certain regularity results for the time discrete solution are establishes rigorously.Based on these regularity results,we prove the unconditional L2-H1 error estimates for the finite element solution of MNSE.Finally,some numerical examples are carried out to demonstrate both accuracy and efficiency of the fully discrete finite element scheme.
基金Supported by the National Natural Science Foundation of China(51379125,51490675,11432009,51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of the Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for Doctoral Students
文摘Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.
基金Supported by the National Natural Science Foundation of China(No.10971203,11271340,11101384)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20094101110006)
文摘In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ||·||div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.
基金This work was supported in part by the National Natural Science Foundation of China (Grant No, 11061021), the Program of Higher-level Talents of Inner Mongolia University (No. Z200901004), and the Scientific Research Projection of Higher Schools of Inner Mongolia (Nos. N J10006, N J10016, NJZZ12011).
文摘A splitting positive definite mixed finite element method is proposed for second-order viscoelasticity wave equation. The proposed procedure can be split into three independent symmetric positive definite integro-differential sub-system and does not need to solve a coupled system of equations. Error estimates are derived for both semidiscrete and fully discrete schemes. The existence and uniqueness for semidiscrete scheme are proved. Finally, a numerical example is provided to illustrate the efficiency of the method.
基金supported by the National Natural Science Foundation of China(Grant Nos.11571223,and 51404160)the Science and Technology Innovation Foundation of Higher Education of Shanxi Province(Grant No.2016163)
文摘By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IEFG) method is used for 2D transient heat conduction problems, and the finite difference method is applied in the splitting direction. The discretized system equation is obtained based on the Galerkin weak form of 2D problem; the essential boundary conditions are imposed with the penalty method; and the finite difference method is employed in the time domain. Four exemplary problems are chosen to verify the efficiency of the DSEFG method. The numerical solutions show that the efficiency and precision of the DSEFG method are greater than ones of the IEFG method for 3D problems.
基金Projects(50490274, 50534030) supported by the National Natural Science Foundation of ChinaProject supported by the Natural Science Foundatin of Hunan Province, China
文摘Split Hopkinson pressure bar(SHPB) apparatus, usually used for testing behavior of material in median and high strain-rate, is now widely used in the study of rock dynamic constitutive relation, damage evolvement mechanism and energy consumption. However, the possible reasons of sampling disturbance, machining error and so on often lead to the scattering of test results, and bring ultimate difficulty for forming general test conclusion. Based on the stochastic finite element method, the uncertain parameters of specimen density ρs, specimen radius Rs, specimen elastic modulus Es and specimen length Ls in the data processing of SHPB test were considered, and the correlation between the parameters and the test results was analyzed. The results show that the specimen radius Rs has direct correlation with the test result, improving the accuracy in preparing and measuring of specimen is an effective way to improve the accuracy of test and minish the scattering of results for SHPB test.
基金The project supported by the National Natural Science Foundation of China(19972053 and 59889101)the National Outstanding Young Scientist Fund of China(59925104)
文摘Based on the thermal kinetic and mass conservation, a series of controlling equations for the finite element are derived and related programs are developed to simulate the damage microcrack healing process controlled by surface diffusion. Two kinds of typical models for microcrack splitting are proposed, i.e., the grain boundary energy existing on the crack surface and residual stresses applying on the crack surface. And the conditions of microcrack splitting in the two models are given as a function of the microcrack aspect ratio. The microcrack with traction-free surfaces will directly evolve into a spheroid.
基金supported by the National Natural Science Foundation of China (10732010,10972010 and 11028206)
文摘In this study, the 3-dimensional discrete element method is firstly introduced to explain the fracturing damage process of the dynamic split experiment of a special brittle glass ZnS. The corresponding dynamic split experiment is also performed by using the split Hopkinson pressure bar. Then the numerical results correspond closely to those obtained by experiments, and the fracturing damage mode shows that the sample under high strain rate loading would crack along vertical diameter in the band region between two loading edges, which differs from the static damage mode. Furthermore, by comparing a group of contrast numerical tests, the numerical results prove that loading area upon the top side of samples would influence the fracture mode of dynamic split experiments, which indicates that the narrow loading plane is better.