Semihyper-Reduction for Finite Element Structures With Nonlinear Surface Loads on the Basis of Stress Modes

Author(s):  
Lukas Koller ◽  
Wolfgang Witteveen ◽  
Florian Pichler ◽  
Peter Fischer

Abstract Model reduction via projection is a common method to accelerate time integration of finite element (FE) structures by reducing the number of degrees-of-freedom (DOFs). However, nonlinear state-dependent surface loads are usually computed based on the nonreduced DOFs of the FE model. When a considerably high number of DOFs are involved in the nonlinear surface loads, their computation becomes a bottleneck. This paper presents a general approach for reduced time integration and reduced force computation for FE models. The required force trial vectors can be computed easily and systematically out of deformation trial vectors, commonly called “modes.” Those force trial vectors, which we call “stress modes,” can be determined a priori so that a nonlinear computation of the full system is not necessary. The new idea in this contribution is that stress recovery is used to decrease the number of equations for the force computation. A general framework for semihyper-reduction (SHR) is developed and its practical implementation is discussed. The term SHR is introduced because it is an intermediate approach between the straight-forward method of using the FE DOFs and pure hyper-reduction (HR) where the FE DOFs are omitted for computing state-depended surface loads. In order to demonstrate the proposed SHR approach practically, a numerical example of a planar crank drive is given, where a hydrodynamic lubrication film separates piston and cylinder. Thereby, very good result quality has been observed in comparison to a finite difference reference solution.

2020 ◽  
Vol 1 (1) ◽  
pp. 93-102
Author(s):  
Carsten Strzalka ◽  
◽  
Manfred Zehn ◽  

For the analysis of structural components, the finite element method (FEM) has become the most widely applied tool for numerical stress- and subsequent durability analyses. In industrial application advanced FE-models result in high numbers of degrees of freedom, making dynamic analyses time-consuming and expensive. As detailed finite element models are necessary for accurate stress results, the resulting data and connected numerical effort from dynamic stress analysis can be high. For the reduction of that effort, sophisticated methods have been developed to limit numerical calculations and processing of data to only small fractions of the global model. Therefore, detailed knowledge of the position of a component’s highly stressed areas is of great advantage for any present or subsequent analysis steps. In this paper an efficient method for the a priori detection of highly stressed areas of force-excited components is presented, based on modal stress superposition. As the component’s dynamic response and corresponding stress is always a function of its excitation, special attention is paid to the influence of the loading position. Based on the frequency domain solution of the modally decoupled equations of motion, a coefficient for a priori weighted superposition of modal von Mises stress fields is developed and validated on a simply supported cantilever beam structure with variable loading positions. The proposed approach is then applied to a simplified industrial model of a twist beam rear axle.


Author(s):  
B Ashby ◽  
C Bortolozo ◽  
A Lukyanov ◽  
T Pryer

Summary In this article, we present a goal-oriented adaptive finite element method for a class of subsurface flow problems in porous media, which exhibit seepage faces. We focus on a representative case of the steady state flows governed by a nonlinear Darcy–Buckingham law with physical constraints on subsurface-atmosphere boundaries. This leads to the formulation of the problem as a variational inequality. The solutions to this problem are investigated using an adaptive finite element method based on a dual-weighted a posteriori error estimate, derived with the aim of reducing error in a specific target quantity. The quantity of interest is chosen as volumetric water flux across the seepage face, and therefore depends on an a priori unknown free boundary. We apply our method to challenging numerical examples as well as specific case studies, from which this research originates, illustrating the major difficulties that arise in practical situations. We summarise extensive numerical results that clearly demonstrate the designed method produces rapid error reduction measured against the number of degrees of freedom.


Author(s):  
Rinat Galiautdinov

In this article, the author considers the possibility of applying modern IT technologies to implement information processing algorithms in UAV motion control system. Filtration of coordinates and motion parameters of objects under a priori uncertainty is carried out using nonlinear adaptive filters: Kalman and Bayesian filters. The author considers numerical methods for digital implementation of nonlinear filters based on the convolution of functions, the possibilities of neural networks and fuzzy logic for solving the problems of tracking UAV objects (or missiles), the math model of dynamics, the features of the practical implementation of state estimation algorithms in the frame of added additional degrees of freedom. The considered algorithms are oriented on solving the problems in real time using parallel and cloud computing.


2004 ◽  
Vol 01 (01) ◽  
pp. 1-15 ◽  
Author(s):  
TED BELYTSCHKO ◽  
HAO CHEN

An enrichment technique for accurately modeling two dimensional crack propagation within the framework of the finite element method is presented. The technique uses an enriched basis that spans the asymptotic dynamic crack-tip solution. The enrichment functions and their spatial derivatives are able to exactly reproduce the asymptotic displacement field and strain field for a moving crack. The stress intensity factors for Mode I and Mode II are taken as additional degrees of freedom. An explicit time integration scheme is used to solve the resulting discrete equations. Numerical simulations of linear elastodynamic problems are reported to demonstrate the accuracy and potential of the technique.


Author(s):  
Lassaad Ben Fekih ◽  
Georges Kouroussis ◽  
David Wattiaux ◽  
Olivier Verlinden ◽  
Christophe De Fruytier

An approach is proposed to identify the modal properties of a subsystem made up of an arbitrary chosen inner module of embedded space equipment. An experimental modal analysis was carried out along the equipment transverse direction with references taken onto its outer housing. In parallel, a numerical model using the finite element (FE) method was developed to correlate with the measured results. A static Guyan reduction has led to a set of master degrees of freedom in which the experimental mode shapes were expanded. An updating technique consisting in minimizing the dynamic residual induced by the FE model and the measurements has been investigated. A last verification has consisted in solving the numeric model composed of the new mass and stiffness matrices obtained by means of a minimization of the error in the constitutive equation method.


Author(s):  
Tamer Wasfy

A new technique for modeling rigid bodies undergoing spatial motion using an explicit time-integration finite element code is presented. The key elements of the technique are: (a) use of the total rotation matrix relative to the inertial frame to measure the rotation of the rigid bodies; (b) time-integration of the rotational equations of motion in a body fixed (material) frame, with the resulting incremental rotations added to the total rotation matrix; (c) penalty formulation for creating connection points (virtual nodes which do not add extra degrees of freedom) on the rigid-body where joints can be placed. The use of the rotation matrix along with incremental rotation updates circumvents the problem of singularities associated with other types of three and four parameter rotation measures. Benchmark rigid multibody dynamics problems are solved to demonstrate the accuracy of the present technique.


2020 ◽  
Vol 54 (4) ◽  
pp. 1309-1337
Author(s):  
Zhaonan Dong ◽  
Emmanuil H. Georgoulis ◽  
Tristan Pryer

Recovered Finite Element Methods (R-FEM) have been recently introduced in Georgoulis and Pryer [Comput. Methods Appl. Mech. Eng. 332 (2018) 303–324]. for meshes consisting of simplicial and/or box-type elements. Here, utilising the flexibility of the R-FEM framework, we extend their definition to polygonal and polyhedral meshes in two and three spatial dimensions, respectively. An attractive feature of this framework is its ability to produce arbitrary order polynomial conforming discretizations, yet involving only as many degrees of freedom as discontinuous Galerkin methods over general polygonal/polyhedral meshes with potentially many faces per element. A priori error bounds are shown for general linear, possibly degenerate, second order advection-diffusion-reaction boundary value problems. A series of numerical experiments highlight the good practical performance of the proposed numerical framework.


2019 ◽  
Vol 9 (21) ◽  
pp. 4495 ◽  
Author(s):  
Mucha

Hybrid simulation is a technique for testing mechanical systems. It applies to structures with elements hard or impossible to model numerically. These elements are tested experimentally by straining them by means of actuators, while the rest of the system is simulated numerically using a finite element method (FEM). Data is interchanged between experiment and simulation. The simulation is performed in real-time in order to accurately recreate the dynamic behavior in the experiment. FEM is very computationally demanding, and for systems with a great number of degrees of freedom (DOFs), real-time simulation with small-time steps (ensuring high accuracy) may require powerful computing hardware or may even be impossible. The author proposed to swap the finite element (FE) model with an artificial neural network (ANN) to significantly lower the computational cost of the real-time algorithm. The presented examples proved that the computational cost could be reduced by at least one number of magnitude while maintaining high accuracy of the simulation; however, obtaining appropriate ANN was not trivial and might require many attempts.


Author(s):  
Tong Y. Yi ◽  
Parviz E. Nikravesh

Abstract This paper presents a method for identifying the free-free modes of a structure by utilizing the vibration data of the same structure with boundary conditions. In modal formulations for flexible body dynamics, modal data are primary known quantities that are obtained either experimentally or analytically. The vibration measurements may be obtained for a flexible body that is constrained differently than its boundary conditions in a multibody system. For a flexible body model in a multibody system, depending upon the formulation used, we may need a set of free-free modal data or a set of constrained modal data. If a finite element model of the flexible body is available, its vibration data can be obtained analytically under any desired boundary conditions. However, if a finite element model is not available, the vibration data may be determined experimentally. Since experimentally measured vibration data are obtained for a flexible body supported by some form of boundary conditions, we may need to determine its free-free vibration data. The aim of this study is to extract, based on experimentally obtained vibration data, the necessary free-free frequencies and the associated modes for flexible bodies to be used in multibody formulations. The available vibration data may be obtained for a structure supported either by springs or by fixed boundary conditions. Furthermore, the available modes may be either a complete set; i.e., as many modes as the number of degrees of freedom of the associated FE model is available, or it can be an incomplete set.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Yan Shang ◽  
Huanpu Wu

PurposeThe paper aims to propose two new 8-node quadrilateral membrane elements with good distortion tolerance for the modified couple stress elasticity based on the unsymmetric finite element method (FEM).Design/methodology/approachThe nodal rotation degrees of freedom (DOFs) are introduced into the virtual work principle and constrained by the penalty function for approximating the test functions of the physical rotation and curvature. Therefore, only the C0 continuity instead of C1 continuity is required for the displacement during the element construction. The first unsymmetric element assumes the test functions of the displacement and strain using the standard 8-node isoparametric interpolations, while these test functions in the second model are further enhanced by the nodal rotation DOFs. Besides, the trial functions in these two elements are constructed based on the stress functions that can a priori satisfy related governing equations.FindingsThe benchmark tests show that both the two elements can efficiently simulate the size-dependent plane problems, exhibiting good numerical accuracies and high distortion tolerances. In particular, they can still exactly reproduce the constant couple stress state when the element shape deteriorates severely into the degenerated triangle. Moreover, it can also be observed that the second element model, in which the linked interpolation technique is used, has better performance than the first one, especially in capturing the steep gradients of the physical rotations.Originality/valueAs the proposed new elements use only three DOFs per node, they can be readily incorporated into the existing finite element (FE) programs. Thus, they are of great benefit to analysis of size-dependent membrane behaviors of micro/nano structures.


Sign in / Sign up

Export Citation Format

Share Document