Procedures of preparation of numerical analysis,consisting in a simulation of cooperation of three different media: steel,liquid and gas undergoes dynamic load were discussed.Modelling of the initial static load of th...Procedures of preparation of numerical analysis,consisting in a simulation of cooperation of three different media: steel,liquid and gas undergoes dynamic load were discussed.Modelling of the initial static load of the mechanical system was presented.By using the MSC.Software products the following exemplary computer simulations were made: dynamic load impact on the hydraulic leg as well as effectiveness of the hydraulic leg protection against overload with help of gas accumulator.展开更多
In this paper, the lubrication problem in nmmerical sindation of forming processes is considered.After enumerationg the difficulties encountered when trying to solve such a problem with the finite element method, a ...In this paper, the lubrication problem in nmmerical sindation of forming processes is considered.After enumerationg the difficulties encountered when trying to solve such a problem with the finite element method, a generalization of the formaulation of Liu[4-6] for the thin film hydrodynamic lubrication re- gime is presented. This method is then aplied to a strip rolling simulation,using the Arbitrary La- grangian eulerian (ALE) formalism.展开更多
It is presented a buckling analysis ofmultilayered plates by the finite elements method with mixed unknowns (displacement and rotations, and transverse interlaminar stresses). In the mixed model, each layer is analy...It is presented a buckling analysis ofmultilayered plates by the finite elements method with mixed unknowns (displacement and rotations, and transverse interlaminar stresses). In the mixed model, each layer is analyzed as a single plate, where the continuity of displacement is ensured by Lagrange multipliers which represent static variables. This procedure is easy to be implemented from the computational point of view. A methodology to solve the eigen problem is presented based on the inverse iteration method. The model is verified successfully with results obtained by other authors.展开更多
The paper is devoted to formulations of decay and mapped elastodynamic infinite elements, based on modified Bessel shape functions. These elements are appropriate for Soil-Structure Interaction (SSI) problems, solve...The paper is devoted to formulations of decay and mapped elastodynamic infinite elements, based on modified Bessel shape functions. These elements are appropriate for Soil-Structure Interaction (SSI) problems, solved in time or frequency domain and can be treated as a new form of the recently proposed Elastodynamic Infinite Elements with United Shape Functions (EIEUSF) infinite elements. The formulation of 2D Horizontal type Infinite Elements (HIE) is demonstrated here, but by similar techniques 2D Vertical (VIE) and 2D Comer (CIE) Infinite Elements can also be formulated. Using elastodynamic infinite elements is the easier and appropriate way to achieve an adequate simulation including basic aspects of Soil-Structure Interaction. Continuity along the artificial boundary (the line between finite and infinite elements) is discussed as well and the application of the proposed elastodynamic infinite elements in the Finite Element Method (FEM) is explained in brief. Finally, a numerical example shows the computational efficiency of the proposed infinite elements.展开更多
Stress concentration and large displacements are usual problems in the components of the structure of agricultural machinery such harvesters coffee, and that finite element method (FEM) can be a tool to minimize its e...Stress concentration and large displacements are usual problems in the components of the structure of agricultural machinery such harvesters coffee, and that finite element method (FEM) can be a tool to minimize its effects. The goal of this paper is to get results of stresses and displacements of a coffee harvester structure by using FEM for static simulation. The main parts of the coffee harvester analyzed were: engine frame, body right and left sides, front and rear end, main beam, coffee reservoir, wheels and fuel tank. Two different design concepts of a coffee harvester machine were analyzed (structure with rear wheels aligned and misaligned) and the results were compared. It was observed that the model with rear wheels misaligned showed maximum displacement lower than the model with rear wheels aligned. Although higher stress was found in the rear wheels misaligned, it was observed that average stresses for the misaligned wheels design were lower in most structural components analyzed. Based on FEM results, the coffee harvester machine with misaligned rear wheels was built and subjected to operational tests without showing any structural failure.展开更多
On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node fi...On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node finite elements to compute the compression distortion of the prestressed anchorage segment. The results indicate that when the prestressed force acts on the anchorage head and segment, the stresses and the displacements in the rock around the anchorage head and segment concentrate on the zone center with the anchor axis, and they decrease with exponential forms. Therefore,the stresses and the displacement spindles are formed. The calculating results of the infinite element are close to the theoretical results. This indicates the method is right. This article introduces a new way to study the mechanism of prestressed anchors. The obtained results have an important role in the research of the anchor mechanism and engineering application.展开更多
A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements f...A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.展开更多
Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular character...Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials' auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.展开更多
This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families l...This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families like Daubechies, Interpolets possess rational filter coefficients, are smooth, symmetric and therefore more suitable for use in numerical methods. Expressions for stiffness and mass matrices are developed based on connection coefficients, which are inner products of basis functions and their derivatives. An example in 1-D was formulated using Central Difference and Newmark schemes for time differentiation. Encouraging results were obtained even for large time steps. Results obtained in 2-D are compared with the standard Finite Difference Method for validation.展开更多
This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations. When k is even, the averaging numerical flux (the average of left an...This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations. When k is even, the averaging numerical flux (the average of left and right limits for the discontinuous finite element at nodes) has the optimal-order ultraconvergence 2k + 2. For nanlinear Hamiltonian systems (e.g., SchrSdinger equation and Kepler system) with momentum conservation, the discontinuous finite element methods preserve momentum at nodes. These properties are confirmed by numerical experiments.展开更多
The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can ov...The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.展开更多
In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discr...In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discrete and full-discrete mixed finite elements for the equations, and obtain the optimal L-2 error estimates.展开更多
Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 n...Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O(h2) and O(h3), respectively. Numerical tests confirm the theoretical analysis.展开更多
The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conf...The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.展开更多
The study of damage in rotating machineries is of fundamental interest in the fields of machine and structure design. A rotating system, supported by bearings and under some dynamic conditions, can generate a variety ...The study of damage in rotating machineries is of fundamental interest in the fields of machine and structure design. A rotating system, supported by bearings and under some dynamic conditions, can generate a variety of problems that are encountered in many different types of rotating machines. One of these problems is the unbalance due to non-homogeneous mass distribution along the shaft. One of the techniques which are widespread today is the identification of parameters and excitation forces that may well followed by monitoring the evolution and change of possible variations of these parameters. Although several methods for the identification of unbalance excitation force are available in the literature, none of them can be considered unrestricted to be applied for all rotating systems. In this study, two methodologies to identify unknown excitations, such as unbalance, have been proposed. This project refers to the analysis of unbalanced forces from displacement parameters and speed by using methods of identification by Fourier series and Legendre polynomials together with the finite element method, state observers in reasons of the problem of absence of signs of rotational displacement, bandpass filter were used to noise suppression of the data collected from the experimental part, Quasi-Newton method to minimize a function in which the bearing stiffness and its damping are unknowns, and also the experimental verification of the methodology, using for this system owned by a rotary mechanical vibrations of the Department of Mechanical Engineering of Faculty of Engineering, campus of llha Solteira.展开更多
Natural frequencies for multilayer plates are calculated by mixed finite element method. The main object of this paper is to use the mixed model for multilayer plates, analyzing each layer as an isolated plate, where ...Natural frequencies for multilayer plates are calculated by mixed finite element method. The main object of this paper is to use the mixed model for multilayer plates, analyzing each layer as an isolated plate, where the continuity of displacements is achieved by Lagrange multipliers (representing static variables). This procedure allows us to work with any model for single plate (so as to ensure the proper behavior of each layer), and the complexity of the multilayer system is avoided by ensuring the condition of displacements by the Lagrange multipliers (static variables). The plate is discretized by finite element modeling based on a primary hybrid model, where the domain is divided by quadrilateral, both for the displacement field and static variables. This mixed element for plates was implemented and several examples of vibrations have been verified successfully by the results obtained by other methods in the literature.展开更多
In the processes of manufacturing, MT (machine tools) plays an important role in the manufacture of work pieces with complex and high dimensional and geometric accuracy. Much of the errors of a machine tool are thos...In the processes of manufacturing, MT (machine tools) plays an important role in the manufacture of work pieces with complex and high dimensional and geometric accuracy. Much of the errors of a machine tool are those which are thermally induced which are from internal and external heat sources acting on the machine. In this paper, a methodology for determining and analyzing the thermal deformation of machine tools using FEM (finite element method) and ANN (artificial neural networks) is presented. After modeling the machine using FEM is defined the location of the heat sources, it is possible to obtain the temperature gradient and the corresponding thermal deformation at predetermined periods. Results obtained with simulations using the software NX.7.5 showed that this methodology is an effective tool in determining the thermal deformation of the machine, correlating the temperature reading at strategic points with volumetric deformation at the tool tip. Therefore, the thermal analysis of the errors in the pair tool part can be established. After training and validation process, the network will be able to make the prediction of thermal errors just stating the temperature values of specific points of each heat source, providing a way for compensation of thermally induced errors.展开更多
Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolatio...Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolation error of anisotropic elements based on Newton's formula of polynomial interpolation as well as its applications.展开更多
Infinitesimal-rotation finite elements allow creating a linear problem that can be exploited to systematically reduce the number of coordinates and obtain efficient solutions for a wide range of applications,including...Infinitesimal-rotation finite elements allow creating a linear problem that can be exploited to systematically reduce the number of coordinates and obtain efficient solutions for a wide range of applications,including those governed by nonlinear equations.This paper discusses the limitations of conventional infinitesimal-rotation finite elements(FE)in capturing correctly the initial stress-free reference-configuration geometry,and explains the effect of these limitations on the definition of the inertia used in the motion description.An alternative to conventional infinitesimal-rotation finite elements is a new class of elements that allow developing inertia expressions written explicitly in terms of constant coefficients that define accurately the reference-configuration geometry.It is shown that using a geometrically inconsistent(GI)approach that introduces the infinitesimal-rotation coordinates from the outset to replace the interpolation-polynomial coefficients is the main source of the failure to capture correctly the reference-configuration geometry.On the other hand,by using a geometrically consistent(GC)approach that employs the position gradients of the absolute nodal coordinate formulation(ANCF)to define the infinitesimal-rotation coordinates,the reference-configuration geometry can be preserved.Two simple examples of straight and tapered beams are used to demonstrate the basic differences between the two fundamentally different approaches used to introduce the infinitesimal-rotation coordinates.The analysis presented in this study sheds light on the differences between the incremental co-rotational solution procedure,widely used in computational structural mechanics,and the non-incremental floating frame of reference formulation(FFR),widely used in multibody system(MBS)dynamics.展开更多
We prove that the error estimates of a large class of nonconforming finite elements are dominated by their approximation errors, which means that the well-known Cea’s lemma is still valid for these nonconforming fini...We prove that the error estimates of a large class of nonconforming finite elements are dominated by their approximation errors, which means that the well-known Cea’s lemma is still valid for these nonconforming finite element methods. Furthermore, we derive the error estimates in both energy and L2 norms under the regularity assumption u ∈ H1+s(Ω) with any s 】 0. The extensions to other related problems are possible.展开更多
文摘Procedures of preparation of numerical analysis,consisting in a simulation of cooperation of three different media: steel,liquid and gas undergoes dynamic load were discussed.Modelling of the initial static load of the mechanical system was presented.By using the MSC.Software products the following exemplary computer simulations were made: dynamic load impact on the hydraulic leg as well as effectiveness of the hydraulic leg protection against overload with help of gas accumulator.
文摘In this paper, the lubrication problem in nmmerical sindation of forming processes is considered.After enumerationg the difficulties encountered when trying to solve such a problem with the finite element method, a generalization of the formaulation of Liu[4-6] for the thin film hydrodynamic lubrication re- gime is presented. This method is then aplied to a strip rolling simulation,using the Arbitrary La- grangian eulerian (ALE) formalism.
文摘It is presented a buckling analysis ofmultilayered plates by the finite elements method with mixed unknowns (displacement and rotations, and transverse interlaminar stresses). In the mixed model, each layer is analyzed as a single plate, where the continuity of displacement is ensured by Lagrange multipliers which represent static variables. This procedure is easy to be implemented from the computational point of view. A methodology to solve the eigen problem is presented based on the inverse iteration method. The model is verified successfully with results obtained by other authors.
文摘The paper is devoted to formulations of decay and mapped elastodynamic infinite elements, based on modified Bessel shape functions. These elements are appropriate for Soil-Structure Interaction (SSI) problems, solved in time or frequency domain and can be treated as a new form of the recently proposed Elastodynamic Infinite Elements with United Shape Functions (EIEUSF) infinite elements. The formulation of 2D Horizontal type Infinite Elements (HIE) is demonstrated here, but by similar techniques 2D Vertical (VIE) and 2D Comer (CIE) Infinite Elements can also be formulated. Using elastodynamic infinite elements is the easier and appropriate way to achieve an adequate simulation including basic aspects of Soil-Structure Interaction. Continuity along the artificial boundary (the line between finite and infinite elements) is discussed as well and the application of the proposed elastodynamic infinite elements in the Finite Element Method (FEM) is explained in brief. Finally, a numerical example shows the computational efficiency of the proposed infinite elements.
文摘Stress concentration and large displacements are usual problems in the components of the structure of agricultural machinery such harvesters coffee, and that finite element method (FEM) can be a tool to minimize its effects. The goal of this paper is to get results of stresses and displacements of a coffee harvester structure by using FEM for static simulation. The main parts of the coffee harvester analyzed were: engine frame, body right and left sides, front and rear end, main beam, coffee reservoir, wheels and fuel tank. Two different design concepts of a coffee harvester machine were analyzed (structure with rear wheels aligned and misaligned) and the results were compared. It was observed that the model with rear wheels misaligned showed maximum displacement lower than the model with rear wheels aligned. Although higher stress was found in the rear wheels misaligned, it was observed that average stresses for the misaligned wheels design were lower in most structural components analyzed. Based on FEM results, the coffee harvester machine with misaligned rear wheels was built and subjected to operational tests without showing any structural failure.
文摘On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node finite elements to compute the compression distortion of the prestressed anchorage segment. The results indicate that when the prestressed force acts on the anchorage head and segment, the stresses and the displacements in the rock around the anchorage head and segment concentrate on the zone center with the anchor axis, and they decrease with exponential forms. Therefore,the stresses and the displacement spindles are formed. The calculating results of the infinite element are close to the theoretical results. This indicates the method is right. This article introduces a new way to study the mechanism of prestressed anchors. The obtained results have an important role in the research of the anchor mechanism and engineering application.
基金supported by the National Natural Science Foundation of China(Nos.11271273 and 11271298)
文摘A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.
文摘Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials' auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.
文摘This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families like Daubechies, Interpolets possess rational filter coefficients, are smooth, symmetric and therefore more suitable for use in numerical methods. Expressions for stiffness and mass matrices are developed based on connection coefficients, which are inner products of basis functions and their derivatives. An example in 1-D was formulated using Central Difference and Newmark schemes for time differentiation. Encouraging results were obtained even for large time steps. Results obtained in 2-D are compared with the standard Finite Difference Method for validation.
基金Project supported by the National Natural Science Foundation of China (No. 10771063)
文摘This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations. When k is even, the averaging numerical flux (the average of left and right limits for the discontinuous finite element at nodes) has the optimal-order ultraconvergence 2k + 2. For nanlinear Hamiltonian systems (e.g., SchrSdinger equation and Kepler system) with momentum conservation, the discontinuous finite element methods preserve momentum at nodes. These properties are confirmed by numerical experiments.
文摘The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.
文摘In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discrete and full-discrete mixed finite elements for the equations, and obtain the optimal L-2 error estimates.
基金Project supported by the National Natural Science Foundation of China (Nos. 10771198 and 11071226)the Foundation of International Science and Technology Cooperation of Henan Province
文摘Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O(h2) and O(h3), respectively. Numerical tests confirm the theoretical analysis.
文摘The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.
文摘The study of damage in rotating machineries is of fundamental interest in the fields of machine and structure design. A rotating system, supported by bearings and under some dynamic conditions, can generate a variety of problems that are encountered in many different types of rotating machines. One of these problems is the unbalance due to non-homogeneous mass distribution along the shaft. One of the techniques which are widespread today is the identification of parameters and excitation forces that may well followed by monitoring the evolution and change of possible variations of these parameters. Although several methods for the identification of unbalance excitation force are available in the literature, none of them can be considered unrestricted to be applied for all rotating systems. In this study, two methodologies to identify unknown excitations, such as unbalance, have been proposed. This project refers to the analysis of unbalanced forces from displacement parameters and speed by using methods of identification by Fourier series and Legendre polynomials together with the finite element method, state observers in reasons of the problem of absence of signs of rotational displacement, bandpass filter were used to noise suppression of the data collected from the experimental part, Quasi-Newton method to minimize a function in which the bearing stiffness and its damping are unknowns, and also the experimental verification of the methodology, using for this system owned by a rotary mechanical vibrations of the Department of Mechanical Engineering of Faculty of Engineering, campus of llha Solteira.
文摘Natural frequencies for multilayer plates are calculated by mixed finite element method. The main object of this paper is to use the mixed model for multilayer plates, analyzing each layer as an isolated plate, where the continuity of displacements is achieved by Lagrange multipliers (representing static variables). This procedure allows us to work with any model for single plate (so as to ensure the proper behavior of each layer), and the complexity of the multilayer system is avoided by ensuring the condition of displacements by the Lagrange multipliers (static variables). The plate is discretized by finite element modeling based on a primary hybrid model, where the domain is divided by quadrilateral, both for the displacement field and static variables. This mixed element for plates was implemented and several examples of vibrations have been verified successfully by the results obtained by other methods in the literature.
文摘In the processes of manufacturing, MT (machine tools) plays an important role in the manufacture of work pieces with complex and high dimensional and geometric accuracy. Much of the errors of a machine tool are those which are thermally induced which are from internal and external heat sources acting on the machine. In this paper, a methodology for determining and analyzing the thermal deformation of machine tools using FEM (finite element method) and ANN (artificial neural networks) is presented. After modeling the machine using FEM is defined the location of the heat sources, it is possible to obtain the temperature gradient and the corresponding thermal deformation at predetermined periods. Results obtained with simulations using the software NX.7.5 showed that this methodology is an effective tool in determining the thermal deformation of the machine, correlating the temperature reading at strategic points with volumetric deformation at the tool tip. Therefore, the thermal analysis of the errors in the pair tool part can be established. After training and validation process, the network will be able to make the prediction of thermal errors just stating the temperature values of specific points of each heat source, providing a way for compensation of thermally induced errors.
基金the National Nutural Science Foundation of China(Grant Nos.10771198,10590353)
文摘Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolation error of anisotropic elements based on Newton's formula of polynomial interpolation as well as its applications.
基金supported,in part,by the National Science Foundation(Grant 1852510).
文摘Infinitesimal-rotation finite elements allow creating a linear problem that can be exploited to systematically reduce the number of coordinates and obtain efficient solutions for a wide range of applications,including those governed by nonlinear equations.This paper discusses the limitations of conventional infinitesimal-rotation finite elements(FE)in capturing correctly the initial stress-free reference-configuration geometry,and explains the effect of these limitations on the definition of the inertia used in the motion description.An alternative to conventional infinitesimal-rotation finite elements is a new class of elements that allow developing inertia expressions written explicitly in terms of constant coefficients that define accurately the reference-configuration geometry.It is shown that using a geometrically inconsistent(GI)approach that introduces the infinitesimal-rotation coordinates from the outset to replace the interpolation-polynomial coefficients is the main source of the failure to capture correctly the reference-configuration geometry.On the other hand,by using a geometrically consistent(GC)approach that employs the position gradients of the absolute nodal coordinate formulation(ANCF)to define the infinitesimal-rotation coordinates,the reference-configuration geometry can be preserved.Two simple examples of straight and tapered beams are used to demonstrate the basic differences between the two fundamentally different approaches used to introduce the infinitesimal-rotation coordinates.The analysis presented in this study sheds light on the differences between the incremental co-rotational solution procedure,widely used in computational structural mechanics,and the non-incremental floating frame of reference formulation(FFR),widely used in multibody system(MBS)dynamics.
基金supported by the Special Funds for Major State Basic Research Project(Grant No.2005CB321701)
文摘We prove that the error estimates of a large class of nonconforming finite elements are dominated by their approximation errors, which means that the well-known Cea’s lemma is still valid for these nonconforming finite element methods. Furthermore, we derive the error estimates in both energy and L2 norms under the regularity assumption u ∈ H1+s(Ω) with any s 】 0. The extensions to other related problems are possible.