scholarly journals Impact of the unit cell choice on the efficiency of dispersion curve calculations using GBMS

2021 ◽  
pp. 1-28
Author(s):  
Vanessa Cool ◽  
Lucas Van Belle ◽  
Claus Claeys ◽  
Elke Deckers ◽  
Wim Desmet

Abstract Recently, the potential of metamaterials and phononic crystals to cope with conflicting requirements of obtaining lightweight structures with desirable noise and vibration properties has been demonstrated. These, often periodic, structures are commonly studied based on their representative unit cell of which the vibro-acoustic performance is examined by means of their wave propagation, visualized by dispersion curves. Typically, the unit cell is discretized using a finite element technique to capture the possibly complex geometry. This leads to a high computation cost for the dispersion curve calculation which can be strongly reduced by applying modal based model order reduction techniques such as the (generalized) Bloch mode synthesis. In this paper, the choice of the unit cell is shown to have an impact on the dispersion curve calculation time. Moreover, the efficiency of (generalized) Bloch mode synthesis strongly depends on the unit cell choice. The highest reduction in computation time is accomplished when the number of boundary degrees-of-freedom is limited.


2021 ◽  
Vol 263 (4) ◽  
pp. 2102-2113
Author(s):  
Vanessa Cool ◽  
Lucas Van Belle ◽  
Claus Claeys ◽  
Elke Deckers ◽  
Wim Desmet

Metamaterials, i.e. artificial structures with unconventional properties, have shown to be highly potential lightweight and compact solutions for the attenuation of noise and vibrations in targeted frequency ranges, called stop bands. In order to analyze the performance of these metamaterials, their stop band behavior is typically predicted by means of dispersion curves, which describe the wave propagation in the corresponding infinite periodic structure. The input for these calculations is usually a finite element model of the corresponding unit cell. Most common in literature are 2D plane metamaterials, which often consist of a plate host structure with periodically added masses or resonators. In recent literature, however, full 3D metamaterials are encountered which are periodic in all three directions and which enable complete, omnidirectional stop bands. Although these 3D metamaterials have favorable vibro-acoustic characteristics, the computational cost to analyze them quickly increases with unit cell model size. Model order reduction techniques are important enablers to overcome this problem. In this work, the Bloch Mode Synthesis (BMS) and generalized BMS (GBMS) reduction techniques are extended from 2D to 3D periodic structures. Through several verifications, it is demonstrated that dispersion curve calculation times can be strongly reduced, while accurate stop band predictions are maintained.



Materials ◽  
2021 ◽  
Vol 14 (2) ◽  
pp. 271
Author(s):  
Jun-Jun Zhai ◽  
Xiang-Xia Kong ◽  
Lu-Chen Wang

A homogenization-based five-step multi-scale finite element (FsMsFE) simulation framework is developed to describe the time-temperature-dependent viscoelastic behavior of 3D braided four-directional composites. The current analysis was performed via three-scale finite element models, the fiber/matrix (microscopic) representative unit cell (RUC) model, the yarn/matrix (mesoscopic) representative unit cell model, and the macroscopic solid model with homogeneous property. Coupling the time-temperature equivalence principle, multi-phase finite element approach, Laplace transformation and Prony series fitting technology, the character of the stress relaxation behaviors at three scales subject to variation in temperature is investigated, and the equivalent time-dependent thermal expansion coefficients (TTEC), the equivalent time-dependent thermal relaxation modulus (TTRM) under micro-scale and meso-scale were predicted. Furthermore, the impacts of temperature, structural parameters and relaxation time on the time-dependent thermo-viscoelastic properties of 3D braided four-directional composites were studied.



2021 ◽  
Vol 263 (1) ◽  
pp. 5301-5309
Author(s):  
Luca Alimonti ◽  
Abderrazak Mejdi ◽  
Andrea Parrinello

Statistical Energy Analysis (SEA) often relies on simplified analytical models to compute the parameters required to build the power balance equations of a coupled vibro-acoustic system. However, the vibro-acoustic of modern structural components, such as thick sandwich composites, ribbed panels, isogrids and metamaterials, is often too complex to be amenable to analytical developments without introducing further approximations. To overcome this limitation, a more general numerical approach is considered. It was shown in previous publications that, under the assumption that the structure is made of repetitions of a representative unit cell, a detailed Finite Element (FE) model of the unit cell can be used within a general and accurate numerical SEA framework. In this work, such framework is extended to account for structural-acoustic coupling. Resonant as well as non-resonant acoustic and structural paths are formulated. The effect of any acoustic treatment applied to coupling areas is considered by means of a Generalized Transfer Matrix (TM) approach. Moreover, the formulation employs a definition of pressure loads based on the wavenumber-frequency spectrum, hence allowing for general sources to be fully represented without simplifications. Validations cases are presented to show the effectiveness and generality of the approach.



Author(s):  
Pavel Karban ◽  
David Pánek ◽  
Ivo Doležel

Purpose A novel technique for control of complex physical processes based on the solution of their sufficiently accurate models is presented. The technique works with the model order reduction (MOR), which significantly accelerates the solution at a still acceptable uncertainty. Its advantages are illustrated with an example of induction brazing. Design/methodology/approach The complete mathematical model of the above heat treatment process is presented. Considering all relevant nonlinearities, the numerical model is reduced using the orthogonal decomposition and solved by the finite element method (FEM). It is cheap compared with classical FEM. Findings The proposed technique is applicable in a wide variety of linear and weakly nonlinear problems and exhibits a good degree of robustness and reliability. Research limitations/implications The quality of obtained results strongly depends on the temperature dependencies of material properties and degree of nonlinearities involved. In case of multiphysics problems characterized by low nonlinearities, the results of solved problems differ only negligibly from those solved on the full model, but the computation time is lower by two and more orders. Yet, however, application of the technique in problems with stronger nonlinearities was not fully evaluated. Practical implications The presented model and methodology of its solution may represent a basis for design of complex technologies connected with induction-based heat treatment of metal materials. Originality/value Proposal of a sophisticated methodology for solution of complex multiphysics problems established the MOR technology that significantly accelerates their solution at still acceptable errors.



2018 ◽  
Vol 140 (7) ◽  
Author(s):  
Mazher Iqbal Mohammed ◽  
Ian Gibson

Highly organized, porous architectures leverage the true potential of additive manufacturing (AM) as they can simply not be manufactured by any other means. However, their mainstream usage is being hindered by the traditional methodologies of design which are heavily mathematically orientated and do not allow ease of controlling geometrical attributes. In this study, we aim to address these limitations through a more design-driven approach and demonstrate how complex mathematical surfaces, such as triply periodic structures, can be used to generate unit cells and be applied to design scaffold structures in both regular and irregular volumes in addition to hybrid formats. We examine the conversion of several triply periodic mathematical surfaces into unit cell structures and use these to design scaffolds, which are subsequently manufactured using fused filament fabrication (FFF) additive manufacturing. We present techniques to convert these functions from a two-dimensional surface to three-dimensional (3D) unit cell, fine tune the porosity and surface area, and examine the nuances behind conversion into a scaffold structure suitable for 3D printing. It was found that there are constraints in the final size of unit cell that can be suitably translated through a wider structure while still allowing for repeatable printing, which ultimately restricts the attainable porosities and smallest printed feature size. We found this limit to be approximately three times the stated precision of the 3D printer used this study. Ultimately, this work provides guidance to designers/engineers creating porous structures, and findings could be useful in applications such as tissue engineering and product light-weighting.



2021 ◽  
Author(s):  
K. BUKENYA ◽  
M. N. OLAYA ◽  
E. J. PINEDA ◽  
M. MAIARU

Woven polymer matrix composites (PMCs) are leveraged in aerospace applications for their desirable specific properties, yet they are vulnerable to high residual stresses during manufacturing and their complex geometry makes experimental results difficult to observe. Process modeling is needed to characterize the effects of the curing and predict end stress states. Finite element software can be used to model woven architectures, however accurate representation of processing conditions remains a challenge when it comes to selecting boundary conditions. The effect of BCs on process-induced stress within woven PMCs is studied. The commercial Finite Element Analysis (FEA) software Abaqus is coupled with user-written subroutines in a process modeling framework. A two-dimensionally (2D) woven PMC repeating unit cell (RUC) is modeled with TexGen and Abaqus. Virtual curing is imposed on the bulk matrix. The BC study is conducted with Free, Periodic, Flat, and Flat-Free configurations. Results show that the end stress state is sensitive to the boundary condition assumptions. Flat BC results show great agreement with Periodic BCs. Residual stress results from process modeling are then compared with a linear-elastic thermal cooldown analysis in Abaqus. Cooldown results indicate an overestimation in matrix stresses compared with process modeling.



2018 ◽  
Vol 226 ◽  
pp. 04036
Author(s):  
Yuriy M. Manatskov ◽  
Torsten Bertram ◽  
Danil V. Shaykhutdinov ◽  
Nikolay I. Gorbatenko

Complex dynamic linear systems of equations are solved by numerical iterative methods, which need much computation and are timeconsuming ones, and the optimization stage requires repeated solution of these equation systems that increases the time on development. To shorten the computation time, various methods can be applied, among them preliminary (estimated) calculation or oversimple models calculation, however, while testing and optimizing the full model is used. Reduced order models are very popular in solving this problem. The main idea of a reduced order model is to find a simplified model that may reflect the required properties of the original model as accurately as possible. There are many methods for the model order reduction, which have their advantages and disadvantages. In this article, a method based on Krylov subspaces and SVD methods is considered. A numerical experiments is given.



2002 ◽  
Vol 256 (4) ◽  
pp. 605-627 ◽  
Author(s):  
M.T. BAH ◽  
A. BHASKAR ◽  
A.J. KEANE


1994 ◽  
Vol 9 (4) ◽  
pp. 272-279 ◽  
Author(s):  
H. Toraya ◽  
T. Ochiai

The accuracy of the unit-cell parameters refined by using the whole-powder-pattern decomposition method is discussed. Powders of W, ZnO, TiO2, BaTiO3 Mg2SiO4, Al2SiO5 (+α-SiO2), and monoclinic ZrO2 were used as test samples. Two internal standard reference materials of Si and CeO2 and two types of powder diffractometers were used for data collections. The systematic peak-shift was corrected by determining the unit-cell parameters and the error function simultaneously during the whole-pattern-fitting. The estimated standard deviations for sample means ranged from <10 ppm (10−6) in cubic symmetry to 20∼50 ppm in monoclinic symmetry. These analyses could be carried out almost automatically in a computation time of less than l min for each sample on a workstation. The use of symmetric experimental profiles, obtained by the suppression of axial divergence, is very effective and of essential importance for improving the accuracy of unit-cell parameters.



1997 ◽  
Vol 119 (3) ◽  
pp. 439-447 ◽  
Author(s):  
M. P. Castanier ◽  
G. O´ttarsson ◽  
C. Pierre

The analysis of the response statistics of mistuned turbomachinery rotors requires an expensive Monte Carlo simulation approach. Simple lumped parameter models capture basic localization effects but do not represent well actual engineering structures without a difficult parameter identification. Current component mode analysis techniques generally require a minimum number of degrees of freedom which is too large for running Monte Carlo simulations at a reasonable cost. In the present work, an order reduction method is introduced which is capable of generating reasonably accurate, very low order models of tuned or mistuned bladed disks. This technique is based on component modes of vibration found from a finite element analysis of a single disk-blade sector. It is shown that the phenomenon of mode localization is well captured by the reduced order modeling technique.



Sign in / Sign up

Export Citation Format

Share Document