The typical quadrangular and triangular elements for thin plate bending based on Kirchhoff assumptions are the non- conforming elements with low computational accuracy and limitative application range in fmite element...The typical quadrangular and triangular elements for thin plate bending based on Kirchhoff assumptions are the non- conforming elements with low computational accuracy and limitative application range in fmite element method(FEM). Some compatible elements can be developed by the means of supplementing correction functions, increasing nodes in element or on the boundaries, expanding nodal degrees of freedom(DOF), etc, but these elements are inconvenient to apply in practice for the high calculation complexity. In this paper, in order to overcome the defects of thin plate bending finite element, numerical manifold method(NMM) was introduced to solve thin plate bending deformation problem. Rectangular mesh was adopted as mathematical mesh to form f'mite element cover system, and then 16-cover manifold element was proposed. Numerical manifold formulas were constructed on the basis of minimum potential energy principle, displacement boundary conditions are implemented by penalty function method, and all the element matrixes were derived in details. The 16-cover element has a simple calculation process for employing only the transverse displacement cover DOFs as the basic unknown variables, and has been proved to meet the requirements of completeness and full compatibility. As an application, the presented 16-cover element has been used to analyze bending deformation of square thin plate under different loads and boundary conditions, and the results show that numerical manifold method with compatible element, compared with finite element method, can improve computational accuracy and convergence greatly.展开更多
This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress...This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress can be naturally satisfied. The gene-rallzed variational principles with mixed hybrid incompatible finite elements are alsopresented and proved, and they can reduce the computation of incompatible elements indynamics of viscous barotropic flows.展开更多
There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement compo...There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.展开更多
Phylogenetic relationships among the five key angiosperm lineages, Ceratophyllum, Chloranthaceae, eudicots, magnoliids, and monocots, have resisted resolution despite several large-scale analyses sampling taxa and cha...Phylogenetic relationships among the five key angiosperm lineages, Ceratophyllum, Chloranthaceae, eudicots, magnoliids, and monocots, have resisted resolution despite several large-scale analyses sampling taxa and characters extensively and using various analytical methods. Meanwhile, compatibility methods, which were explored together with parsimony and likelihood methods during the early development stage of phylogenetics, have been greatly under-appreciated and not been used to analyze the massive amount of sequence data to recon-struct the basal angiosperm phylogeny. In this study, we used a compatibility method on a data set of eight genes (mitochondrial atp1, matR, and nad5, plastid atpB, matK, rbcL, and rpoC2, and nuclear 18S rDNA) gathered in an earlier study. We selected two sets of characters that are compatible with more of the other characters than a random character would be with at probabilities of pM<0.1 and pM<0.5 respectively. The resulting data matrices were subjected to parsimony and likelihood bootstrap analyses. Our unrooted parsimony analyses showed that Cerato-phyllum was immediately related to eudicots, this larger lineage was immediately related to magnoliids, and monocots were closely related to Chloranthaceae. All these relationships received 76%–96% bootstrap support. A likelihood analysis of the 8 gene pM<0.5 compatible site matrix recovered the same topology but with low support. Likelihood analyses of other compatible site matrices produced different topologies that were all weakly sup-ported. The topology reconstructed in the parsimony analyses agrees with the one recovered in the previous study using both parsimony and likelihood methods when no character was eliminated. Parts of this topology have also been recovered in several earlier studies. Hence, this topology plausibly reflects the true relationships among the five key angiosperm lineages.展开更多
A new flux-based hybrid subcell-remapping algorithm for staggered multimaterial arbitrary Lagrangian-Eulerian (MMALE) methods is presented. This new method is an effective generalization of the original subcell-remapp...A new flux-based hybrid subcell-remapping algorithm for staggered multimaterial arbitrary Lagrangian-Eulerian (MMALE) methods is presented. This new method is an effective generalization of the original subcell-remapping method to the multi-material regime (LOUBERE, R. and SHASHKOV,M. A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods. Journal of Computational Physics, 209, 105–138 (2005)). A complete remapping procedure of all fluid quantities is described detailedly in this paper. In the pure material regions, remapping of mass and internal energy is performed by using the original subcell-remapping method. In the regions near the material interfaces, remapping of mass and internal energy is performed with the intersection-based fluxes where intersections are performed between the swept regions and pure material polygons in the Lagrangian mesh, and an approximate approach is then introduced for constructing the subcell mass fluxes. In remapping of the subcell momentum, the mass fluxes are used to construct the momentum fluxes by multiplying a reconstructed velocity in the swept region. The nodal velocity is then conservatively recovered. Some numerical examples simulated in the full MMALE regime and several purely cyclic remapping examples are presented to prove the properties of the remapping method.展开更多
The interaction between poly(methymethacrylate) (PMMA) and poly(vinyl chloride) (PVC) has been studied indilute urea solutions of dimethylformamide (DMF) at 28℃ using a dilute solution viscometry method. The results ...The interaction between poly(methymethacrylate) (PMMA) and poly(vinyl chloride) (PVC) has been studied indilute urea solutions of dimethylformamide (DMF) at 28℃ using a dilute solution viscometry method. The results show thatthe polymer mixtures are compatible in DMF solution in the absence of urea. The influence of urea addition on the degree ofcompatibility of the polymer mixtures has been studied in terms of the compatibility parameters (△b_m and △[η]_m). It wasfound that the compatibility of the polymer mixtures is decreased with increasing urea addition, passing through a minimumat 0.5 M urea.展开更多
We propose an explicit,single-step discontinuous Galerkin method on moving grids using the arbitrary Lagrangian-Eulerian approach for one-dimensional Euler equations.The grid is moved with the local fluid velocity mod...We propose an explicit,single-step discontinuous Galerkin method on moving grids using the arbitrary Lagrangian-Eulerian approach for one-dimensional Euler equations.The grid is moved with the local fluid velocity modified by some smoothing,which is found to con-siderably reduce the numerical dissipation introduced by Riemann solvers.The scheme preserves constant states for any mesh motion and we also study its positivity preservation property.Local grid refinement and coarsening are performed to maintain the mesh qual-ity and avoid the appearance of very small or large cells.Second,higher order methods are developed and several test cases are provided to demonstrate the accuracy of the proposed scheme.展开更多
In this paper,we present the negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)method solving nonlinear hyperbolic equations with smooth solutions.The smoothness-increasing ac...In this paper,we present the negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)method solving nonlinear hyperbolic equations with smooth solutions.The smoothness-increasing accuracy-conserving(SIAC)filter is a post-processing technique to enhance the accuracy of the discontinuous Galerkin(DG)solutions.This work is the essential step to extend the SIAC filter to the moving mesh for nonlinear problems.By the post-processing theory,the negative norm estimates are vital to get the superconvergence error estimates of the solutions after post-processing in the L2 norm.Although the SIAC filter has been extended to nonuniform mesh,the analysis of fil-tered solutions on the nonuniform mesh is complicated.We prove superconvergence error estimates in the negative norm for the ALE-DG method on moving meshes.The main dif-ficulties of the analysis are the terms in the ALE-DG scheme brought by the grid velocity field,and the time-dependent function space.The mapping from time-dependent cells to reference cells is very crucial in the proof.The numerical results also confirm the theoreti-cal proof.展开更多
We study the convergence and asymptotic compatibility of higher order collocation methods for nonlocal operators inspired by peridynamics,a nonlocal formulation of continuum mechanics.We prove that the methods are opt...We study the convergence and asymptotic compatibility of higher order collocation methods for nonlocal operators inspired by peridynamics,a nonlocal formulation of continuum mechanics.We prove that the methods are optimally convergent with respect to the polynomial degree of the approximation.A numerical method is said to be asymptotically compatible if the sequence of approximate solutions of the nonlocal problem converges to the solution of the corresponding local problem as the horizon and the grid sizes simultaneously approach zero.We carry out a calibration process via Taylor series expansions and a scaling of the nonlocal operator via a strain energy density argument to ensure that the resulting collocation methods are asymptotically compatible.We fnd that,for polynomial degrees greater than or equal to two,there exists a calibration constant independent of the horizon size and the grid size such that the resulting collocation methods for the nonlocal difusion are asymptotically compatible.We verify these fndings through extensive numerical experiments.展开更多
Free-interface dual-compatibility modal synthesis method(compatibility of both force and displacement on interfaces)is introduced to large-scale civil engineering structure to enhance computation efficiency. The basic...Free-interface dual-compatibility modal synthesis method(compatibility of both force and displacement on interfaces)is introduced to large-scale civil engineering structure to enhance computation efficiency. The basic equations of the method are first set up, and then the mode cut-off principle and the dividing principle are proposed. MATLAB is used for simulation in different frame structures. The simulation results demonstrate the applicability of this substructure method to civil engineering structures and the correctness of the proposed mode cut-off principle. Studies are also conducted on how to divide the whole structure for better computation efficiency while maintaining better precision. It is observed that the geometry and material properties should be considered, and the synthesis results would be more precise when the inflection points of the mode shapes are taken into consideration. Furthermore, the simulation performed on a large-scale high-rise connected structure further proves the feasibility and efficiency of this modal synthesis method compared with the traditional global method. It is also concluded from the simulation results that the fewer number of DOFs in each substructure will result in better computation efficiency, but too many substructures will be time-consuming due to the tedious synthesis procedures. Moreover, the substructures with free interface will introduce errors and reduce the precision dramatically, which should be avoided.展开更多
Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-freque...Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-frequency electromagnetic fields on the environment and which are capable of modifying and destroying our ecosystem. Particular interest was given in this article due to a massive influx of population living near high voltage lines. The analysis and simulation of the influence of low frequency electromagnetic fields on living beings in the vicinity of high voltage sources 132 kV and 220 kV in urban areas in DR Congo is the subject of our research with a view to estimating the level of exposure of humans to low frequency electromagnetic fields. To carry out our research, we used the classic method of analyzing the field produced near a high voltage line based on Maxwell’s image theory, the Maxwell-Gauss theorem and Maxwell-Ampère theorem to model and quantify low-frequency electromagnetic fields in the vicinity of high-voltage lines. The 2D FDTD numerical formulation was developed from telegraphers’ equations and allowed us to obtain models of current and voltage induced by electromagnetic fields on living beings below and near HV lines. The different simulations carried out on the proposed models illustrate the effects of the electrical and geometric parameters of the pylons on the distribution of the electromagnetic field in the vicinity of the HV lines. The results obtained were compared to the safety limits recommended by the standards.展开更多
基金supported by National Natural Science Foundation of China (Grant No. 50775044, Grant No. 50975050)Guangdong Provincial and Ministry of Education Industry-University-Research Integration Project of China (Grant No. 2009B090300044)
文摘The typical quadrangular and triangular elements for thin plate bending based on Kirchhoff assumptions are the non- conforming elements with low computational accuracy and limitative application range in fmite element method(FEM). Some compatible elements can be developed by the means of supplementing correction functions, increasing nodes in element or on the boundaries, expanding nodal degrees of freedom(DOF), etc, but these elements are inconvenient to apply in practice for the high calculation complexity. In this paper, in order to overcome the defects of thin plate bending finite element, numerical manifold method(NMM) was introduced to solve thin plate bending deformation problem. Rectangular mesh was adopted as mathematical mesh to form f'mite element cover system, and then 16-cover manifold element was proposed. Numerical manifold formulas were constructed on the basis of minimum potential energy principle, displacement boundary conditions are implemented by penalty function method, and all the element matrixes were derived in details. The 16-cover element has a simple calculation process for employing only the transverse displacement cover DOFs as the basic unknown variables, and has been proved to meet the requirements of completeness and full compatibility. As an application, the presented 16-cover element has been used to analyze bending deformation of square thin plate under different loads and boundary conditions, and the results show that numerical manifold method with compatible element, compared with finite element method, can improve computational accuracy and convergence greatly.
文摘This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress can be naturally satisfied. The gene-rallzed variational principles with mixed hybrid incompatible finite elements are alsopresented and proved, and they can reduce the computation of incompatible elements indynamics of viscous barotropic flows.
文摘There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.
基金supported by an Early Career Award(DEB 0332298)an ATOL grant(DEB 0431239)from the National Science Foundation,USA.
文摘Phylogenetic relationships among the five key angiosperm lineages, Ceratophyllum, Chloranthaceae, eudicots, magnoliids, and monocots, have resisted resolution despite several large-scale analyses sampling taxa and characters extensively and using various analytical methods. Meanwhile, compatibility methods, which were explored together with parsimony and likelihood methods during the early development stage of phylogenetics, have been greatly under-appreciated and not been used to analyze the massive amount of sequence data to recon-struct the basal angiosperm phylogeny. In this study, we used a compatibility method on a data set of eight genes (mitochondrial atp1, matR, and nad5, plastid atpB, matK, rbcL, and rpoC2, and nuclear 18S rDNA) gathered in an earlier study. We selected two sets of characters that are compatible with more of the other characters than a random character would be with at probabilities of pM<0.1 and pM<0.5 respectively. The resulting data matrices were subjected to parsimony and likelihood bootstrap analyses. Our unrooted parsimony analyses showed that Cerato-phyllum was immediately related to eudicots, this larger lineage was immediately related to magnoliids, and monocots were closely related to Chloranthaceae. All these relationships received 76%–96% bootstrap support. A likelihood analysis of the 8 gene pM<0.5 compatible site matrix recovered the same topology but with low support. Likelihood analyses of other compatible site matrices produced different topologies that were all weakly sup-ported. The topology reconstructed in the parsimony analyses agrees with the one recovered in the previous study using both parsimony and likelihood methods when no character was eliminated. Parts of this topology have also been recovered in several earlier studies. Hence, this topology plausibly reflects the true relationships among the five key angiosperm lineages.
基金Project supported by the China Postdoctoral Science Foundation(No.2017M610823)
文摘A new flux-based hybrid subcell-remapping algorithm for staggered multimaterial arbitrary Lagrangian-Eulerian (MMALE) methods is presented. This new method is an effective generalization of the original subcell-remapping method to the multi-material regime (LOUBERE, R. and SHASHKOV,M. A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods. Journal of Computational Physics, 209, 105–138 (2005)). A complete remapping procedure of all fluid quantities is described detailedly in this paper. In the pure material regions, remapping of mass and internal energy is performed by using the original subcell-remapping method. In the regions near the material interfaces, remapping of mass and internal energy is performed with the intersection-based fluxes where intersections are performed between the swept regions and pure material polygons in the Lagrangian mesh, and an approximate approach is then introduced for constructing the subcell mass fluxes. In remapping of the subcell momentum, the mass fluxes are used to construct the momentum fluxes by multiplying a reconstructed velocity in the swept region. The nodal velocity is then conservatively recovered. Some numerical examples simulated in the full MMALE regime and several purely cyclic remapping examples are presented to prove the properties of the remapping method.
文摘The interaction between poly(methymethacrylate) (PMMA) and poly(vinyl chloride) (PVC) has been studied indilute urea solutions of dimethylformamide (DMF) at 28℃ using a dilute solution viscometry method. The results show thatthe polymer mixtures are compatible in DMF solution in the absence of urea. The influence of urea addition on the degree ofcompatibility of the polymer mixtures has been studied in terms of the compatibility parameters (△b_m and △[η]_m). It wasfound that the compatibility of the polymer mixtures is decreased with increasing urea addition, passing through a minimumat 0.5 M urea.
文摘We propose an explicit,single-step discontinuous Galerkin method on moving grids using the arbitrary Lagrangian-Eulerian approach for one-dimensional Euler equations.The grid is moved with the local fluid velocity modified by some smoothing,which is found to con-siderably reduce the numerical dissipation introduced by Riemann solvers.The scheme preserves constant states for any mesh motion and we also study its positivity preservation property.Local grid refinement and coarsening are performed to maintain the mesh qual-ity and avoid the appearance of very small or large cells.Second,higher order methods are developed and several test cases are provided to demonstrate the accuracy of the proposed scheme.
基金the fellowship of China Postdoctoral Science Foundation,no:2020TQ0030.Y.Xu:Research supported by National Numerical Windtunnel Project NNW2019ZT4-B08+1 种基金Science Challenge Project TZZT2019-A2.3NSFC Grants 11722112,12071455.X.Li:Research supported by NSFC Grant 11801062.
文摘In this paper,we present the negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)method solving nonlinear hyperbolic equations with smooth solutions.The smoothness-increasing accuracy-conserving(SIAC)filter is a post-processing technique to enhance the accuracy of the discontinuous Galerkin(DG)solutions.This work is the essential step to extend the SIAC filter to the moving mesh for nonlinear problems.By the post-processing theory,the negative norm estimates are vital to get the superconvergence error estimates of the solutions after post-processing in the L2 norm.Although the SIAC filter has been extended to nonuniform mesh,the analysis of fil-tered solutions on the nonuniform mesh is complicated.We prove superconvergence error estimates in the negative norm for the ALE-DG method on moving meshes.The main dif-ficulties of the analysis are the terms in the ALE-DG scheme brought by the grid velocity field,and the time-dependent function space.The mapping from time-dependent cells to reference cells is very crucial in the proof.The numerical results also confirm the theoreti-cal proof.
文摘We study the convergence and asymptotic compatibility of higher order collocation methods for nonlocal operators inspired by peridynamics,a nonlocal formulation of continuum mechanics.We prove that the methods are optimally convergent with respect to the polynomial degree of the approximation.A numerical method is said to be asymptotically compatible if the sequence of approximate solutions of the nonlocal problem converges to the solution of the corresponding local problem as the horizon and the grid sizes simultaneously approach zero.We carry out a calibration process via Taylor series expansions and a scaling of the nonlocal operator via a strain energy density argument to ensure that the resulting collocation methods are asymptotically compatible.We fnd that,for polynomial degrees greater than or equal to two,there exists a calibration constant independent of the horizon size and the grid size such that the resulting collocation methods for the nonlocal difusion are asymptotically compatible.We verify these fndings through extensive numerical experiments.
基金Supported by the National Natural Science Foundation of China(No.51108089)Doctoral Programs Foundation of Ministry of Education of China(No.20113514120005)the Foundation of the Education Department of Fujian Province(No.JA14057)
文摘Free-interface dual-compatibility modal synthesis method(compatibility of both force and displacement on interfaces)is introduced to large-scale civil engineering structure to enhance computation efficiency. The basic equations of the method are first set up, and then the mode cut-off principle and the dividing principle are proposed. MATLAB is used for simulation in different frame structures. The simulation results demonstrate the applicability of this substructure method to civil engineering structures and the correctness of the proposed mode cut-off principle. Studies are also conducted on how to divide the whole structure for better computation efficiency while maintaining better precision. It is observed that the geometry and material properties should be considered, and the synthesis results would be more precise when the inflection points of the mode shapes are taken into consideration. Furthermore, the simulation performed on a large-scale high-rise connected structure further proves the feasibility and efficiency of this modal synthesis method compared with the traditional global method. It is also concluded from the simulation results that the fewer number of DOFs in each substructure will result in better computation efficiency, but too many substructures will be time-consuming due to the tedious synthesis procedures. Moreover, the substructures with free interface will introduce errors and reduce the precision dramatically, which should be avoided.
文摘Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-frequency electromagnetic fields on the environment and which are capable of modifying and destroying our ecosystem. Particular interest was given in this article due to a massive influx of population living near high voltage lines. The analysis and simulation of the influence of low frequency electromagnetic fields on living beings in the vicinity of high voltage sources 132 kV and 220 kV in urban areas in DR Congo is the subject of our research with a view to estimating the level of exposure of humans to low frequency electromagnetic fields. To carry out our research, we used the classic method of analyzing the field produced near a high voltage line based on Maxwell’s image theory, the Maxwell-Gauss theorem and Maxwell-Ampère theorem to model and quantify low-frequency electromagnetic fields in the vicinity of high-voltage lines. The 2D FDTD numerical formulation was developed from telegraphers’ equations and allowed us to obtain models of current and voltage induced by electromagnetic fields on living beings below and near HV lines. The different simulations carried out on the proposed models illustrate the effects of the electrical and geometric parameters of the pylons on the distribution of the electromagnetic field in the vicinity of the HV lines. The results obtained were compared to the safety limits recommended by the standards.