An Approach to Variation Simulation of Final Aircraft Assembly with Presence of Sealant

The article discusses the process of aeronautical structure assembly in the presence of a sealant between the parts to be joined. An attempt to estimate the influence of sealant on assembly quality in terms of variation analysis is presented. The sealant is considered as a highly viscous liquid that is applied to the surfaces of the assembled parts before the start of final assembly. The modeling approach is based on simulation of two-way coupled fluid-structure interaction between fluid sealant and compliant structural parts. Reynolds lubrication approximation is used in the fluid dynamics problem and variational formulation of contact problem combined with static condensation is used in the structural one. The joining of two aircraft panels is used as a numerical test for demonstration of developed approach. Various phenomena connected with the presence of sealant are demonstrated. In particular, the difference in the fastener loosening due to sealant flow between different types of fasteners is investigated. Results of variation simulation show that presence of sealant should be considered among determining factors in the analysis of assembly quality.

Performances of Hybridized-, Embedded-, and Weighted-Interior Penalty Discontinuous Galerkin Methods for Heterogeneous and Anisotropic Diffusion Problems

The present paper discusses families of Interior Penalty Discontinuous Galerkin (IP) methods for solving heterogeneous and anisotropic diffusion problems. Specifically, we focus on distinctive schemes, namely the Hybridized-, Embedded-, and Weighted-IP schemes, leading to final matrixes of different sizes and sparsities. Both the Hybridized- and Embedded-IP schemes are eligible for static condensation, and their degrees of freedom are distributed on the mesh skeleton. In contrast, the unknowns are located inside the mesh elements for the Weighted-IP variant. For a given mesh, it is well-known that the number of degrees of freedom related to the standard Discontinuous Galerkin methods increases more rapidly than those of the skeletal approaches (Hybridized- and Embedded-IP). We then quantify the impact of the static condensation procedure on the computational performances of the different IP classes in terms of robustness, accuracy, and CPU time. To this aim, numerical experiments are investigated by considering strong heterogeneities and anisotropies. We analyze the fixed error tolerance versus the run time and mesh size to guide our performance criterion. We also outlined some relationships between these Interior Penalty schemes. Eventually, we confirm the superiority of the Hybridized- and Embedded-IP schemes, regardless of the mesh, the polynomial degree, and the physical properties (homogeneous, heterogeneous, and/or anisotropic).

p-Multilevel Preconditioners for HHO Discretizations of the Stokes Equations with Static Condensation

We propose a p-multilevel preconditioner for hybrid high-order (HHO) discretizations of the Stokes equation, numerically assess its performance on two variants of the method, and compare with a classical discontinuous Galerkin scheme. An efficient implementation is proposed where coarse level operators are inherited using $$L^2$$ L 2 -orthogonal projections defined over mesh faces and the restriction of the fine grid operators is performed recursively and matrix-free. Both h- and k-dependency are investigated tackling two- and three-dimensional problems on standard meshes and graded meshes. For the two HHO formulations, featuring discontinuous or hybrid pressure, we study how the combination of p-coarsening and static condensation influences the V-cycle iteration. In particular, two different static condensation procedures are considered for the discontinuous pressure HHO variant, resulting in global linear systems with a different number of unknowns and matrix non-zero entries. Interestingly, we show that the efficiency of the solution strategy might be impacted by static condensation options in the case of graded meshes.

The pursuit of a dream, Francisco Javier Sayas and the HDG methods

Franciso Javier Sayas, man of grit and determination, left his hometown of Zaragoza in 2007 in pursuit of a dream, to become a scholar in the USA. I hosted him in Minneapolis, where he spent three long years of an arduous transition before obtaining a permanent position at the University of Delaware. There, he enthusiastically worked on the unfolding of his dream until his life was tragically cut short by cancer, at only 50. In this paper, I try to bring to light the part of his academic life he shared with me. As we both worked on hybridizable discontinuous Galerkin methods, and he wrote a book on the subject, I will tell Javier’s life as it developed around this topic. First, I will show how the ideas of static condensation and hybridization, proposed back in the mid 60s, lead to the introduction of those methods. This background material will allow me to tell the story of the evolution of the hybridizable discontinuous Galerkin methods and describe Javier’s participation in it. Javier faced death with open eyes and poised dignity. I will end with a poem he liked.

Two-field formulations for isogeometric Reissner–Mindlin plates and shells with global and local condensation

In this paper, mixed formulations are presented in the framework of isogeometric Reissner–Mindlin plates and shells with the aim of alleviating membrane and shear locking. The formulations are based on the Hellinger-Reissner functional and use the stress resultants as additional unknowns, which have to be interpolated in appropriate approximation spaces. The additional unknowns can be eliminated by static condensation. In the framework of isogeometric analysis static condensation is performed globally on the patch level, which leads to a high computational cost. Thus, two additional local approaches to the existing continuous method are presented, an approach with discontinuous stress resultant fields at the element boundaries and a reconstructed approach which is blending the local control variables by using weights in order to compute the global ones. Both approaches allow for a static condensation on the element level instead of the patch level. Various numerical examples are investigated in order to verify the accuracy and effectiveness of the different approaches and a comparison to existing elements that include mechanisms against locking is carried out.

Detecting internal resonances during model reduction

Model order reduction of geometrically nonlinear dynamic structures is often achieved via a static condensation procedure, whereby high-frequency modes are assumed to be quasi-statically coupled to a small set of lower frequency modes, which form the reduction basis. This approach is mathematically justifiable for structures characterized by slow/fast dynamics, such as thin plates and slender beams, and has been shown to provide highly accurate results. Nevertheless, selecting the reduction basis without a priori knowledge of the full-order dynamics is a challenging task; retaining redundant modes will lead to computationally suboptimal reduced-order models (ROMs), while omitting dynamically significant modes will lead to inaccurate results, and important features such as internal resonances may not be captured. In this study, we demonstrate how the error associated with static condensation can be efficiently approximated during model reduction. This approximate error can then be used as the basis of a method for predicting when dynamic modal interactions will occur, which will guide the reduction basis selection process. Equivalently, this may serve as a tool for verifying the accuracy of ROMs without the need for full-order simulations. The proposed method is demonstrated using a simple oscillator and a finite element model of a clamped–clamped beam.