Distributed multilevel optimization for complex structures

Optimization problems concerning complex structures with many design variables may entail an unacceptable computational cost. This problem can be reduced considerably with a multilevel approach: A structure consisting of several components is optimized as a whole (global) as well as on the component level. In this paper, an optimization method is discussed with applications in the assessment of the impact of new design considerations in the development of a structure. A strategy based on fully stressed design is applied for optimization problems in linear statics. A global model is used to calculate the interactions (e.g., loads) for each of the components. These components are then optimized using the prescribed interactions, followed by a new global calculation to update the interactions. Mixed discrete and continuous design variables as well as different design configurations are possible. An application of this strategy is presented in the form of the full optimization of a vertical tail plane center box of a generic large passenger aircraft. In linear dynamics, the parametrization of the component interactions is problematic due to the frequency dependence. Hence, a modified method is presented in which the speed of component mode synthesis is used to avoid this parametrization. This method is applied to a simple test case that originates from noise control. \u

Design optimization applied in structural dynamics

This paper introduces the design optimization strategies, especially for structures which have dynamic constraints. Design optimization involves first the modeling and then the optimization of the problem. Utilizing the Finite Element (FE) model of a structure directly in an optimization process requires a long computation time. Therefore the Backpropagation Neural Networks (NNs) are introduced as a so called surrogate model for the FE model. Optimization techniques mentioned in this study cover the Genetic Algorithm (GA) and the Sequential Quadratic Programming (SQP) methods. For the applications of the introduced techniques, a multisegment cantilever beam problem under the constraints of its first and second natural frequency has been selected and solved using four different approaches

Coupled Simulation of Transient Heat Flow and Electric Currents in Thin Wires: Application to Bond Wires in Microelectronic Chip Packaging

This work addresses the simulation of heat flow and electric currents in thin wires. An important application is the use of bond wires in microelectronic chip packaging. The heat distribution is modeled by an electrothermal coupled problem, which poses numerical challenges due to the presence of different geometric scales. The necessity of very fine grids is relaxed by solving and embedding a 1D sub-problem along the wire into the surrounding 3D geometry. The arising singularities are described using de Rham currents. It is shown that the problem is related to fluid flow in porous 3D media with 1D fractures [C. D'Angelo, SIAM Journal on Numerical Analysis 50.1, pp. 194-215, 2012]. A careful formulation of the 1D-3D coupling condition is essential to obtain a stable scheme that yields a physical solution. Elliptic model problems are used to investigate the numerical errors and the corresponding convergence rates. Additionally, the transient electrothermal simulation of a simplified microelectronic chip package as used in industrial applications is presented.Comment: all numerical results can be reproduced by the Matlab code openly available at https://github.com/tc88/ETwireSi

A global-local optimization method for problems in structural dynamics

The optimization of complex structures involving many design variables and constraints can be performed using a multi-level approach: a structure consisting of several components is optimized as a whole (global) and on the component level (local). Earlier work [1], [2], [3], described a multilevel technique developed for the optimization the Airbus A380 vertical tail plane. In this application, a global model is used to calculate the loads on each of the components. These components are then optimized using the prescribed loads, followed by a new global calculation to update the loads. The component optimization strategy is based on Neural Networks (NN) and Genetic Algorithms (GA). This paper describes a strategy that makes this global-local optimization method possible for problems in structural dynamics. It is established that a parametrization of the component interactions (e.g. component loads) is problematic due to frequency dependence. Hence, a modified method is proposed in which the speed of Component Mode Synthesis (CMS) is used to avoid this parametrization. The effectiveness of this method is demonstrated in a test case concerning the placement of sensor and actuator locations in Active Structural Acoustic Control (ASAC). Special attention is paid to the behavior of the optimization strategy

Explicit mixed strain–displacement finite elements for compressible and quasi-incompressible elasticity and plasticity

The final publication is available at Springer via http://dx.doi.org/ 10.1007/s00466-016-1305-zThis paper presents an explicit mixed finite element formulation to address compressible and quasi-incompressible problems in elasticity and plasticity. This implies that the numerical solution only involves diagonal systems of equations. The formulation uses independent and equal interpolation of displacements and strains, stabilized by variational subscales. A displacement sub-scale is introduced in order to stabilize the mean-stress field. Compared to the standard irreducible formulation, the proposed mixed formulation yields improved strain and stress fields. The paper investigates the effect of this enhancement on the accuracy in problems involving strain softening and localization leading to failure, using low order finite elements with linear continuous strain and displacement fields (P1P1 triangles in 2D and tetrahedra in 3D) in conjunction with associative frictional Mohr–Coulomb and Drucker–Prager plastic models. The performance of the strain/displacement formulation under compressible and nearly incompressible deformation patterns is assessed and compared to analytical solutions for plane stress and plane strain situations. Benchmark numerical examples show the capacity of the mixed formulation to predict correctly failure mechanisms with localized patterns of strain, virtually free from any dependence of the mesh directional bias. No auxiliary crack tracking technique is necessary.Peer ReviewedPostprint (author's final draft

Variational Calculation on A=3 and 4 Nuclei with Non-Local Potentials

The application of the hyperspherical harmonic approach to the case of non-local two-body potentials is described. Given the properties of the hyperspherical harmonic functions, there are no difficulties in considering the approach in both coordinate and momentum space. The binding energies and other ground state properties of A=3 and 4 nuclei are calculated using the CD Bonn 2000 and N3LO two-body potentials. The results are shown to be in excellent agreement with corresponding ones obtained by other accurate techniques.Comment: 12 pages, 6 tables, RevTex

First-principles determined charge and orbital interactions in Fe3_3O4_4

The interactions between charge and orbitally ordered $d$-electrons are important in many transition metal oxides. We propose an effective energy model for such interactions, parameterized with DFT+U calculations, so that energy contributions of both electronic and lattice origin can be simultaneously accounted for. The model is applied to the low-temperature phase of magnetite, for which we propose a new ground state structure. The effective interactions on the B-lattice of Fe$_3$O$_4$ can be interpreted in terms of electrostatics and short-range Kugel-Khomskii exchange coupling. The frustration between optimal charge and orbital orderings leads to a complex energy landscape whereby the supercell for the charge ordering, orbital ordering and ionic displacements can all be different.Comment: 15 pages, 4 figure

Layer-resolved magnetic exchange interactions of surfaces of late 3d elements: effects of electronic correlations

We present the results of an ab initio study of magnetic properties of Fe, Co and Ni surfaces. In particular, we discuss their electronic structure and magnetic exchange interactions (Jij), as obtained by means of a combination of density functional theory and dynamical mean-field theory. All studied systems have a pronounced tendency to ferromagnetism both for bulk and surface atoms. The presence of narrow-band surface states is shown to enhance the magnetic moment as well as the exchange couplings. The most interesting results were obtained for the Fe surface where the atoms have a tendency to couple antiferromagnetically with each other. This interaction is relatively small, when compared to interlayer ferromagnetic interaction, and strongly depends on the lattice parameter. Local correlation effects are shown to lead to strong changes of the overall shape of the spectral functions. However, they seem to not play a decisive role on the overall picture of the magnetic couplings studied here. We have also investigated the influence of correlations on the spin and orbital moments of the bulk-like and surface atoms. We found that dynamical correlations in general lead to enhanced values of the orbital moment.Comment: 13 pages, 12 figure