Octasection‐based refinement of finite element approximations of tetrahedral meshes that guarantees shape quality

During the last 2 years, a multidomain formalism for structural dynamics based on a multi‐time‐step algorithm and local linear modal reduction was proposed by Gravouil, Combescure, Herry & Faucher. In the first part of this paper, we extended modal reduction to subdomains undergoing finite rigid‐body rotations. Here, we focus on the consequences of local modal projection (either linear or geometrically non‐linear) on the treatment of interface problems between subdomains. In particular, we address the issues of the invertibility and efficiency of the solution process. We illustrate our propositions with specific interpretations of the examples presented in Part 1 and present an additional example to demonstrate the properties of special sets of modes. Copyright © 2004 John Wiley & Sons, Ltd.     

Locking in the incompressible limit for the element‐free Galerkin method

Volumetric locking (locking in the incompressible limit) for linear elastic isotropic materials is studied in the context of the element‐free Galerkin method. The modal analysis developed here shows that the number of non‐physical locking modes is independent of the dilation parameter (support of the interpolation functions). Thus increasing the dilation parameter does not suppress locking. Nevertheless, an increase in the dilation parameter does reduce the energy associated with the non‐physical locking modes; thus, in part, it alleviates the locking phenomena. This is shown for linear and quadratic orders of consistency. Moreover, the biquadratic order of consistency, as in finite elements, improves the locking behaviour. Although more locking modes are present in the element‐free Galerkin method with quadratic consistency than with standard biquadratic finite elements. Finally, numerical examples are shown to validate the modal analysis. In particular, the conclusions of the modal analysis are also confirmed in an elastoplastic example. Copyright © 2001 John Wiley & Sons, Ltd.     

Model reduction for LPV systems based on approximate modal decomposition

The paper presents a novel model order reduction technique for large‐scale linear parameter‐varying (LPV) systems. The approach is based on decoupling the original dynamics into smaller dimensional LPV subsystems that can be independently reduced by parameter‐varying reduction methods. The decomposition starts with the construction of a modal transformation that separates the modal subsystems. Hierarchical clustering is applied then to collect the dynamically similar modal subsystems into larger groups. The resulting parameter‐varying subsystems are then independently reduced. This approach substantially differs from most of the previously proposed LPV model reduction techniques, since it performs manipulations on the LPV model itself, instead of on a set of linear time‐invariant models defined at fixed scheduling parameter values. Therefore, the interpolation, which is often a challenging part in reduction techniques, is inherently solved. The applicability of the developed algorithm is thoroughly investigated and demonstrated by numerical case studies.     

A semi‐local spectral/hp element solver for linear elasticity problems

We develop an efficient semi‐local method for speeding up the solution of linear systems arising in spectral/hp element discretization of the linear elasticity equations. The main idea is to approximate the element‐wise residual distribution with a localization operator we introduce in this paper, and subsequently solve the local linear system. Additionally, we decouple the three directions of displacement in the localization operator, hence enabling the use of an efficient low energy preconditioner for the conjugate gradient solver. This approach is effective for both nodal and modal bases in the spectral/hp element method, but here, we focus on the modal hierarchical basis. In numerical tests, we verify that there is no loss of accuracy in the semi‐local method, and we obtain good parallel scalability and substantial speed‐up compared to the original formulation. In particular, our tests include both structure‐only and fluid‐structure interaction problems, with the latter modeling a 3D patient‐specific brain aneurysm. Copyright © 2014 John Wiley & Sons, Ltd.     

The global modal parameterization for non‐linear model‐order reduction in flexible multibody dynamics

In flexible multibody dynamics, advanced modelling methods lead to high‐order non‐linear differential‐algebraic equations (DAEs). The development of model reduction techniques is motivated by control design problems, for which compact ordinary differential equations (ODEs) in closed‐form are desirable. In a linear framework, reduction techniques classically rely on a projection of the dynamics onto a linear subspace. In flexible multibody dynamics, we propose to project the dynamics onto a submanifold of the configuration space, which allows to eliminate the non‐linear holonomic constraints and to preserve the Lagrangian structure. The construction of this submanifold follows from the definition of a global modal parameterization (GMP): the motion of the assembled mechanism is described in terms of rigid and flexible modes, which are configuration‐dependent. The numerical reduction procedure is presented, and an approximation strategy is also implemented in order to build a closed‐form expression of the reduced model in the configuration space. Numerical and experimental results illustrate the relevance of this approach. Copyright © 2006 John Wiley & Sons, Ltd.     

Multi‐time‐step and two‐scale domain decomposition method for non‐linear structural dynamics

In this paper we propose a method to improve the means of taking into account the specific time‐scale and space‐scale characteristics in time‐dependent non‐linear problems. This method enables the use of arbitrary time steps in each subdomain: these can be coupled by prescribing continuous velocities at the interfaces, which are modelled using a dual Schur formulation. For certain subdomains, in space, we adopt a two‐scale resolution technique inspired by the multigrid methods in order to obtain the part of the solution related to small variation lengths on a refined scale and the part corresponding to large variation lengths on a coarse scale. For non‐linear problems, we propose an algorithm with a single iteration level to deal with both the non‐linear equilibrium and the two space scales thanks to a two‐grid method in which the relaxation steps are performed using a non‐linear, preconditioned conjugate gradient algorithm. Finally, we present an example which demonstrates the feasibility of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

A method for interpolating on manifolds structural dynamics reduced‐order models

A rigorous method for interpolating a set of parameterized linear structural dynamics reduced‐order models (ROMs) is presented. By design, this method does not operate on the underlying set of parameterized full‐order models. Hence, it is amenable to an online real‐time implementation. It is based on mapping appropriately the ROM data onto a tangent space to the manifold of symmetric positive‐definite matrices, interpolating the mapped data in this space and mapping back the result to the aforementioned manifold. Algorithms for computing the forward and backward mappings are offered for the case where the ROMs are derived from a general Galerkin projection method and the case where they are constructed from modal reduction. The proposed interpolation method is illustrated with applications ranging from the fast dynamic characterization of a parameterized structural model to the fast evaluation of its response to a given input. In all cases, good accuracy is demonstrated at real‐time processing speeds. Copyright © 2009 John Wiley & Sons, Ltd.     

A finite element approach combining a reduced‐order system,Padé approximants,and an adaptive frequency windowing for fast multi‐frequency solution of poro‐acoustic problems

In this work, a solution strategy is investigated for the resolution of multi‐frequency structural‐acoustic problems including 3D modeling of poroelastic materials. The finite element method is used, together with a combination of a modal‐based reduction of the poroelastic domain and a Padé‐based reconstruction approach. It thus takes advantage of the reduced‐size of the problem while further improving the computational efficiency by limiting the number of frequency resolutions of the full‐sized problem. An adaptive procedure is proposed for the discretization of the frequency range into frequency intervals of reconstructed solution. The validation is presented on a 3D poro‐acoustic example. Copyright © 2013 John Wiley & Sons, Ltd.     

Optimal modal reduction of vibrating substructures

A structure which consists of a main part and a number of attached substructures is considered. A ‘model reduction’ scheme is developed and applied to each of the discrete substructures. Linear undamped transient vibrational motion of the structure is assumed, with general external forcing and initial conditions. The goal is to replace each discrete substructure by another substructure with a much smaller number of degrees of freedom, while minimizing the effect this reduction has on the dynamic behaviour of the main structure. The approach taken here involves Ritz reduction and the Dirichlet‐to‐Neumann map as analysis tools. The resulting scheme is based on a special form of modal reduction, and is shown to be optimal in a certain sense, for long simulation times. The performance of the scheme is demonstrated via numerical examples, and is compared to that of standard modal reduction. Copyright © 2003 John Wiley & Sons, Ltd.     

Non‐linear finite element modal approach for the large amplitude free vibration of symmetric and unsymmetric composite plates

A theoretical analysis is presented for the large amplitude vibration of symmetric and unsymmetric composite plates using the non‐linear finite element modal reduction method. The problem is first reduced to a set of Duffing‐type modal equations using the finite element modal reduction method. The main advantage of the proposed approach is that no updating of the non‐linear stiffness matrices is needed. Without loss of generality, accurate frequency ratios for the fundamental mode and the higher modes of a composite plate at various values of maximum deflection are then determined by using the Runge–Kutta numerical integration scheme. The procedure for obtaining proper initial conditions for the periodic plate motions is very time consuming. Thus, an alternative scheme (the harmonic balance method) is adopted and assessed, as it was employed to formulate the large amplitude free vibration of beams in a previous study, and the results agreed well with the elliptic solution. The numerical results that are obtained with the harmonic balance method agree reasonably well with those obtained with the Runge–Kutta method. The contribution of each linear mode to the maximum deflection of a plate can also be obtained. The frequency ratios for isotropic and composite plates at various maximum deflections are presented, and convergence of frequencies with the number of finite elements, number of linear modes, and number of harmonic terms is also studied. Copyright © 2005 John Wiley & Sons, Ltd.     

Interface‐reduction for the Craig–Bampton and Rubin method applied to FE–BE coupling with a large fluid–structure interface

Component mode‐based model‐order reduction (MOR) methods like the Craig–Bampton method or the Rubin method are known to be limited to structures with small coupling interfaces. This paper investigates two interface‐reduction methods for application of MOR to systems with large coupling interfaces: for the Craig–Bampton method a direct reduction method based on strain energy considerations is investigated. Additionally, for the Rubin method an iterative reduction scheme is proposed, which incrementally constructs the reduction basis. Hereby, attachment modes are tested if they sufficiently enlarge the spanned subspace of the current reduction basis. If so, the m‐orthogonal part is used to augment the basis. The methods are applied to FE–BE coupled systems in order to predict the vibro‐acoustic behavior of structures, which are partly immersed in water. Hereby, a strong coupling scheme is employed, since for dense fluids the feedback of the acoustic pressure onto the structure is not negligible. For two example structures, the efficiency of the reduction methods with respect to numerical effort, memory consumption and computation time is compared with the exact full‐order solution. Copyright © 2008 John Wiley & Sons, Ltd.     

On the use of modal derivatives for nonlinear model order reduction

Modal derivative is an approach to compute a reduced basis for model order reduction of large‐scale nonlinear systems that typically stem from the discretization of partial differential equations. In this way, a complex nonlinear simulation model can be integrated into an optimization problem or the design of a controller, based on the resulting small‐scale state‐space model. We investigate the approximation properties of modal derivatives analytically and thus lay a theoretical foundation of their use in model order reduction, which has been missing so far. Concentrating on the application field of structural mechanics and structural dynamics, we show that the concept of modal derivatives can also be applied as nonlinear extension of the Craig–Bampton family of methods for substructuring. We furthermore generalize the approach from a pure projection scheme to a novel reduced‐order modeling method that replaces all nonlinear terms by quadratic expressions in the reduced state variables. This complexity reduction leads to a frequency‐preserving nonlinear quadratic state‐space model. Numerical examples with carefully chosen nonlinear model problems and three‐dimensional nonlinear elasticity confirm the analytical properties of the modal derivative reduction and show the potential of the proposed novel complexity reduction methods, along with the current limitations. Copyright © 2016 John Wiley & Sons, Ltd.     

Practical application of a non‐linear earthquake analysis according to DIN EN 1998‐1 / Praktische Anwendung einer nichtlinearen Erdbebenanalyse nach DIN EN 1998‐1

Seismic safety verification can be performed by maintaining constructive rules or by calculation. Verification by calculation can be performed with a linear simplified or linear multi‐modal response spectrum analysis. Alternatively, a non‐linear quasi‐static verification is also possible according to DIN EN 1998‐1, which was not available in DIN 4149. In this article, the non‐linear quasi‐static earthquake verification according to DIN EN 1998‐1 is presented in practice, using the example of a building in Mittenwald/Germany. The verification has been checked and accepted by an independent building supervision report.     

A new methodology to solve non‐linear equation systems using genetic algorithms. Application to combined cyclegas turbine simulation

This paper shows a methodology to sort out the equations of a non‐linear system in order to solve it by the fixed‐point method. The arrangement of the equations is established by a genetic algorithm that deals with a population of possible resolution processes of the system. The method is specially useful in the following situations: first, when the system is very non‐linear and has many variables (where the Newton–Raphson method does not work properly); second, when the number of equations and variables may be altered because the equation system may change in each simulation and, therefore, more than one only solution process is needed if the fixed‐point process is employed. As an example, the methodology has been applied to solve the equation system that models the behaviour of a heat recovery steam generator of a combined cycle power plant at full load and part load conditions. Copyright © 2005 John Wiley & Sons, Ltd.     

A model updating strategy of non‐linear vibrating structures

The objective of this paper is to present a model updating strategy of non‐linear vibrating structures. Because modal analysis is no longer helpful in non‐linear structural dynamics, a special attention is devoted to the features extracted from the proper orthogonal decomposition and one of its non‐linear generalizations based on auto‐associative neural networks. The efficiency of the proposed procedure is illustrated using simulated data from a three‐dimensional portal frame. Copyright © 2004 John Wiley & Sons, Ltd.     

A time‐parallel implicit method for accelerating the solution of non‐linear structural dynamics problems

The parallel implicit time‐integration algorithm (PITA) is among a very limited number of time‐integrators that have been successfully applied to the time‐parallel solution of linear second‐order hyperbolic problems such as those encountered in structural dynamics. Time‐parallelism can be of paramount importance to fast computations, for example, when space‐parallelism is unfeasible as in problems with a relatively small number of degrees of freedom in general, and reduced‐order model applications in particular, or when reaching the fastest possible CPU time is desired and requires the exploitation of both space‐ and time‐parallelisms. This paper extends the previously developed PITA to the non‐linear case. It also demonstrates its application to the reduction of the time‐to‐solution on a Linux cluster of sample non‐linear structural dynamics problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Non‐linear random response of laminated composite shallow shells using finite element modal method

This paper investigates the large‐amplitude multi‐mode random response of thin shallow shells with rectangular planform at elevated temperatures using a finite element non‐linear modal formulation. A thin laminated composite shallow shell element and the system equations of motion are developed. The system equations in structural node degrees‐of‐freedom (DOF) are transformed into modal co‐ordinates, and the non‐linear stiffness matrices are transformed into non‐linear modal stiffness matrices. The number of modal equations is much smaller than the number of equations in structural node DOF. A numerical integration is employed to determine the random response. Thermal buckling deflections are obtained to explain the intermittent snap‐through phenomenon. The natural frequencies of the infinitesimal vibration about the thermally buckled equilibrium positions (BEPs) are studied, and it is found that there is great difference between the frequencies about the primary (positive) and the secondary (negative) BEPs. All three types of motion: (i) linear random vibration about the primary BEP, (ii) intermittent snap‐through between the two BEPs, and (iii) non‐linear large‐amplitude random vibration over the two BEPs, can be predicted. Copyright © 2006 John Wiley & Sons, Ltd.     

Accuracy of one‐step integration schemes for damped/forced linear structural dynamics

This paper proposes an energy‐based measure for the evaluation of the local truncation error of two‐level one‐step integration schemes. The measure applies to multiple degree of freedom systems and does not necessarily require modal reduction to a scalar model; it naturally handles the structural damping and external forcing terms that are generally and mistakenly neglected in error analyses, and it segregates the error associated with the free and forced response components of the problem. To illustrate the approach, two examples associated with the application of the trapezoidal scheme and of a high‐order scheme proposed in the literature are analyzed. The latter reveals the shortcomings of the standard approach that is based on the undamped/unforced linear oscillator and therefore highlights the need for the proposed framework. Indeed, the scheme order of accuracy is below expectation when structural damping or external forcing is considered, in the numerically dissipative setting. Developments on the basis of the time discontinuous Galerkin (TDG) method are then proposed to recover the scheme high‐order accuracy. Additionally, they show the similarity that exists between schemes related to the TDG method and the ones obtained by integration by parts of the equation of motion. Copyright © 2014 John Wiley & Sons, Ltd.     

应用时频分析方法辨识时变系统的模态参数

应用Gabor展开及Hilbert变换进行时变线性系统模态参数的辨识。利用非平稳信号的Gabor变换将结构响应信号表示在时频面上，并通过Gabor展开重构分离模态分量，建立每一阶模态的时变线性模型。对单模态时变线性模型应用Hilbert变换来辨识随时间变化的模态频率和阻尼。通过对刚度和阻尼慢变的两自由度系统模态参数的仿真辨识验证辨识方法的有效性。仿真结果表明：本文方法为时变线性系统的模态参数辨识提供了一条新的途径。