On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these meth...On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming (bi)linear macroelements or (bi)quadratic elements, and the rotation by conforming (bi)linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.展开更多
By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved hav...By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.展开更多
The solid and finite element model of metal pushing type continuously variable transmission are established at speed ratio of i =0 5 and i=2 0. In order to solve the problem of the complicated of structure,the...The solid and finite element model of metal pushing type continuously variable transmission are established at speed ratio of i =0 5 and i=2 0. In order to solve the problem of the complicated of structure,the node node rod discrete finite element model is put forward and the whole system is simplified and established.The natural frequency and mode shape of system are solved by iterative Lanczos reduce method for sensitivity analysis in finite element model.The new method and the result can be used to improve the smoothness of the variable transmission system and to propose the theory for reducing noise at operation.展开更多
In design optimization of crane metal structures, present approaches are based on simple models and mixed variables, which are difficult to use in practice and usually lead to failure of optimized results for rounding...In design optimization of crane metal structures, present approaches are based on simple models and mixed variables, which are difficult to use in practice and usually lead to failure of optimized results for rounding variables. Crane metal structure optimal design(CMSOD) belongs to a constrained nonlinear optimization problem with discrete variables. A novel algorithm combining ant colony algorithm with a mutation-based local search(ACAM) is developed and used for a real CMSOD for the first time. In the algorithm model, the encoded mode of continuous array elements is introduced. This not only avoids the need to round optimization design variables during mixed variable optimization, but also facilitates the construction of heuristic information, and the storage and update of the ant colony pheromone. Together with the proposed ACAM, a genetic algorithm(GA) and particle swarm optimization(PSO) are used to optimize the metal structure of a crane. The optimization results show that the convergence speed of ACAM is approximately 20% of that of the GA and around 11% of that of the PSO. The objective function value given by ACAM is 22.23% less than the practical design value, a reduction of 16.42% over the GA and 3.27% over the PSO. The developed ACAM is an effective intelligent method for CMSOD and superior to other methods.展开更多
The research presented here is focused on the vibration condition of a small volume solder solder ball,which is placed on the surface of a soldering pad and is exerted a pulse modulated continuous wave laser heat sour...The research presented here is focused on the vibration condition of a small volume solder solder ball,which is placed on the surface of a soldering pad and is exerted a pulse modulated continuous wave laser heat source. Finite element method is applied to analyzed the temperature field in the solder ball, and experi- ment is conducted to test the vibration. the results show that,that, the temperature field flucturates with the same frequency as that of the laser pulse, which in turn causes a forced vibration of the same frequency in the liquid solder ball.展开更多
In this paper,we compute a phase field(diffuse interface)model of CahnHilliard type for moving contact line problems governing the motion of isothermal multiphase incompressible fluids.The generalized Navier boundary ...In this paper,we compute a phase field(diffuse interface)model of CahnHilliard type for moving contact line problems governing the motion of isothermal multiphase incompressible fluids.The generalized Navier boundary condition proposed by Qian et al.[1]is adopted here.We discretize model equations using a continuous finite element method in space and a modified midpoint scheme in time.We apply a penalty formulation to the continuity equation which may increase the stability in the pressure variable.Two kinds of immiscible fluids in a pipe and droplet displacement with a moving contact line under the effect of pressure driven shear flow are studied using a relatively coarse grid.We also derive the discrete energy law for the droplet displacement case,which is slightly different due to the boundary conditions.The accuracy and stability of the scheme are validated by examples,results and estimate order.展开更多
A novel multiscale algorithm based on the higher-order continuum at both micro-and macrostructural level is proposed for the consideration of the quasi-brittle damage response of heterogeneous materials.Herein,the mic...A novel multiscale algorithm based on the higher-order continuum at both micro-and macrostructural level is proposed for the consideration of the quasi-brittle damage response of heterogeneous materials.Herein,the microlevel damage is modelled by the degradation of the homogenized stress and tangent stiffness tensors,which are then upscaled to govern the localization at the macrolevel.The C^1 continuity finite element employing a modified case of Mindlin’s form II strain energy density is derived for the softening analysis.To the authors’knowledge,the finite element discretization based on the strain gradient theory is applied for the modeling of damage evolution at the microstructural level for heterogeneous materials for the first time.The advantage of the novel C1 finite element formulation in comparison with the standard finite element discretization in terms of the regularization efficiency as well as the objectivity has been shown.An isotropic damage law is used for the reduction of the constitutive and nonlocal material behaviour,which is necessary for the physically correct description of the localization formation in quasi-brittle materials.The capabilities of the derived finite element to capture the fully developed localization zones are tested on a random representative volume element(RVE)for several different loading cases.By employing the conventional second-order computational homogenization,the microstructural material constitutive response is averaged over the whole RVE area.In order to model the loss of structural integrity when sharp localization is formed across RVE,the specific conditions which detect a completely formed localization zone are developed.A new failure criterion at the microstructural level has been proposed.The derived finite element formulation,as well as the multiscale damage algorithm,are implemented into the finite element program ABAQUS.The capabilities of the presented multiscale scheme to capture the effects of the deformation localization are demonstrated by few benchmark numerical examples.展开更多
We propose a new high order accurate nodal discontinuous Galerkin(DG)method for the solution of nonlinear hyperbolic systems of partial differential equations(PDE)on unstructured polygonal Voronoi meshes.Rather than u...We propose a new high order accurate nodal discontinuous Galerkin(DG)method for the solution of nonlinear hyperbolic systems of partial differential equations(PDE)on unstructured polygonal Voronoi meshes.Rather than using classical polynomials of degree N inside each element,in our new approach the discrete solution is represented by piecewise continuous polynomials of degree N within each Voronoi element,using a continuous finite element basis defined on a subgrid inside each polygon.We call the resulting subgrid basis an agglomerated finite element(AFE)basis for the DG method on general polygons,since it is obtained by the agglomeration of the finite element basis functions associated with the subgrid triangles.The basis functions on each sub-triangle are defined,as usual,on a universal reference element,hence allowing to compute universal mass,flux and stiffness matrices for the subgrid triangles once and for all in a pre-processing stage for the reference element only.Consequently,the construction of an efficient quadrature-free algorithm is possible,despite the unstructured nature of the computational grid.High order of accuracy in time is achieved thanks to the ADER approach,making use of an element-local space-time Galerkin finite element predictor.The novel schemes are carefully validated against a set of typical benchmark problems for the compressible Euler and Navier-Stokes equations.The numerical results have been checked with reference solutions available in literature and also systematically compared,in terms of computational efficiency and accuracy,with those obtained by the corresponding modal DG version of the scheme.展开更多
基金supported by NSFC(11571266,91430106,11171168,11071132)NSFC-RGC(China-Hong Kong)(11661161017)
文摘On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming (bi)linear macroelements or (bi)quadratic elements, and the rotation by conforming (bi)linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.
基金Project supported by the National Natural Science Foundation of China (No.10471038)
文摘By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.
文摘The solid and finite element model of metal pushing type continuously variable transmission are established at speed ratio of i =0 5 and i=2 0. In order to solve the problem of the complicated of structure,the node node rod discrete finite element model is put forward and the whole system is simplified and established.The natural frequency and mode shape of system are solved by iterative Lanczos reduce method for sensitivity analysis in finite element model.The new method and the result can be used to improve the smoothness of the variable transmission system and to propose the theory for reducing noise at operation.
基金Supported by National Natural Science Foundation of China(Grant No.51275329)the Youth Fund Program of Taiyuan University of Science and Technology,China(Grant No.20113014)
文摘In design optimization of crane metal structures, present approaches are based on simple models and mixed variables, which are difficult to use in practice and usually lead to failure of optimized results for rounding variables. Crane metal structure optimal design(CMSOD) belongs to a constrained nonlinear optimization problem with discrete variables. A novel algorithm combining ant colony algorithm with a mutation-based local search(ACAM) is developed and used for a real CMSOD for the first time. In the algorithm model, the encoded mode of continuous array elements is introduced. This not only avoids the need to round optimization design variables during mixed variable optimization, but also facilitates the construction of heuristic information, and the storage and update of the ant colony pheromone. Together with the proposed ACAM, a genetic algorithm(GA) and particle swarm optimization(PSO) are used to optimize the metal structure of a crane. The optimization results show that the convergence speed of ACAM is approximately 20% of that of the GA and around 11% of that of the PSO. The objective function value given by ACAM is 22.23% less than the practical design value, a reduction of 16.42% over the GA and 3.27% over the PSO. The developed ACAM is an effective intelligent method for CMSOD and superior to other methods.
文摘The research presented here is focused on the vibration condition of a small volume solder solder ball,which is placed on the surface of a soldering pad and is exerted a pulse modulated continuous wave laser heat source. Finite element method is applied to analyzed the temperature field in the solder ball, and experi- ment is conducted to test the vibration. the results show that,that, the temperature field flucturates with the same frequency as that of the laser pulse, which in turn causes a forced vibration of the same frequency in the liquid solder ball.
基金Zhenlin Guo is partially supported by the 150th Anniversary Postdoctoral Mobility Grants(2014-15 awards)and the reference number is PMG14-1509Shuangling Dong is supported by the National Natural Science Foundation of China(No.51406098)And also thank the China Postdoctoral Science Foundation(No.2014M560967).
文摘In this paper,we compute a phase field(diffuse interface)model of CahnHilliard type for moving contact line problems governing the motion of isothermal multiphase incompressible fluids.The generalized Navier boundary condition proposed by Qian et al.[1]is adopted here.We discretize model equations using a continuous finite element method in space and a modified midpoint scheme in time.We apply a penalty formulation to the continuity equation which may increase the stability in the pressure variable.Two kinds of immiscible fluids in a pipe and droplet displacement with a moving contact line under the effect of pressure driven shear flow are studied using a relatively coarse grid.We also derive the discrete energy law for the droplet displacement case,which is slightly different due to the boundary conditions.The accuracy and stability of the scheme are validated by examples,results and estimate order.
基金This work has been fully supported by Croatian Science Foundation under the project“Multiscale Numerical Modelling of Material Deformation Responses from Macro-to Nanolevel”(2516).
文摘A novel multiscale algorithm based on the higher-order continuum at both micro-and macrostructural level is proposed for the consideration of the quasi-brittle damage response of heterogeneous materials.Herein,the microlevel damage is modelled by the degradation of the homogenized stress and tangent stiffness tensors,which are then upscaled to govern the localization at the macrolevel.The C^1 continuity finite element employing a modified case of Mindlin’s form II strain energy density is derived for the softening analysis.To the authors’knowledge,the finite element discretization based on the strain gradient theory is applied for the modeling of damage evolution at the microstructural level for heterogeneous materials for the first time.The advantage of the novel C1 finite element formulation in comparison with the standard finite element discretization in terms of the regularization efficiency as well as the objectivity has been shown.An isotropic damage law is used for the reduction of the constitutive and nonlocal material behaviour,which is necessary for the physically correct description of the localization formation in quasi-brittle materials.The capabilities of the derived finite element to capture the fully developed localization zones are tested on a random representative volume element(RVE)for several different loading cases.By employing the conventional second-order computational homogenization,the microstructural material constitutive response is averaged over the whole RVE area.In order to model the loss of structural integrity when sharp localization is formed across RVE,the specific conditions which detect a completely formed localization zone are developed.A new failure criterion at the microstructural level has been proposed.The derived finite element formulation,as well as the multiscale damage algorithm,are implemented into the finite element program ABAQUS.The capabilities of the presented multiscale scheme to capture the effects of the deformation localization are demonstrated by few benchmark numerical examples.
文摘We propose a new high order accurate nodal discontinuous Galerkin(DG)method for the solution of nonlinear hyperbolic systems of partial differential equations(PDE)on unstructured polygonal Voronoi meshes.Rather than using classical polynomials of degree N inside each element,in our new approach the discrete solution is represented by piecewise continuous polynomials of degree N within each Voronoi element,using a continuous finite element basis defined on a subgrid inside each polygon.We call the resulting subgrid basis an agglomerated finite element(AFE)basis for the DG method on general polygons,since it is obtained by the agglomeration of the finite element basis functions associated with the subgrid triangles.The basis functions on each sub-triangle are defined,as usual,on a universal reference element,hence allowing to compute universal mass,flux and stiffness matrices for the subgrid triangles once and for all in a pre-processing stage for the reference element only.Consequently,the construction of an efficient quadrature-free algorithm is possible,despite the unstructured nature of the computational grid.High order of accuracy in time is achieved thanks to the ADER approach,making use of an element-local space-time Galerkin finite element predictor.The novel schemes are carefully validated against a set of typical benchmark problems for the compressible Euler and Navier-Stokes equations.The numerical results have been checked with reference solutions available in literature and also systematically compared,in terms of computational efficiency and accuracy,with those obtained by the corresponding modal DG version of the scheme.