Shape Optimization of Potentially Slender Structures

Author(s):  
Kavous Jorabchi ◽  
Joshua Danczyk ◽  
Krishnan Suresh

Shape optimization lies at the heart of modern engineering design. Through shape optimization, computers can, in theory, ‘synthesize’ engineering artifacts in a fully automated fashion. However, a serious limitation today is that the evolving geometry (during optimization) may become slender, i.e., beam or plate-like. Under such circumstances, modern 3-D computational methods, such as finite element analysis (FEA), will fail miserably, and so will the shape optimization process. Indeed, the recommended method for analyzing slender artifacts is to replace them with 1-D beams/ 2-D plates, prior to discretization and computational analysis, a process referred to as geometric dimensional reduction. Unfortunately explicit geometric reduction is impractical and hard to automate during optimization since one cannot predict a priori when an artifact will become slender. In this paper, we develop an implicit dimensional reduction method where the reduction is achieved through an algebraic process. The proposed method of reduction is computationally equivalent to explicit geometric reduction for comparable computational cost. However, the proposed method can be easily automated and integrated within a shape optimization process, and standard off-the-shelf 3-D finite element packages can be used to implement the proposed methodology.

2015 ◽  
Author(s):  
Γεώργιος Μπαλωμένος

Είναι ευρέως αποδεκτό ότι αβεβαιότητα μπορεί να υπάρξει σε πολλά προβλήματα μηχανικής, όπως στις εισαγόμενες παραμέτρους (φορτίο, γεωμετρία, ιδιότητες υλικών, κλπ.), στις εξαγόμενες παραμέτρους (μετατοπίσεις, τάσεις, κλπ.) και στην μεταξύ τους συσχέτιση. Η ανάλυση αξιοπιστίας είναι σε θέση να αντιμετωπίσει όλες αυτές τις αβεβαιότητες παρέχοντας στους μηχανικούς ακριβείς προβλέψεις σχετικά με την πιθανότητα μιας κατασκευής να συμπεριφέρεται επαρκώς (και με ασφάλεια) στις εξωτερικές φορτίσεις κατά τη διάρκεια της ζωής της.Στην πιθανοτική ανάλυση με τη χρήση των πεπερασμένων στοιχείων (FEA), χρησιμοποιούνται προσεγγιστικές μέθοδοι (όπως οι μέθοδοι σειράς Taylor) προκειμένου να υπολογιστεί ο μέσος όρος και η διακύμανση της δομικής απόκρισης, ενώ η πιθανοτική κατανομή της απόκρισης συνήθως προσεγγίζεται με βάση τη μέθοδο προσομοίωσης Monte Carlo (MCS). Αυτή η διατριβή προάγει την πιθανοτική ανάλυση με τη χρήση των πεπερασμένων στοιχείων FEA, δηλαδή μέσω μη γραμμικής στατικής και δυναμικής ανάλυσης της κατασκευής με τη χρήση αριθμητικών μεθόδων ανάλυσης, συνδυάζοντάς την με την πολλαπλασιαστική μορφή της μεθόδου μείωσης των διαστάσεων (M-DRM). Αυτός ο συνδυασμός επιτρέπει αρκετά ακριβείς εκτιμήσεις τόσο των στατιστικών παραμέτρων (μέση τιμή, τυπική απόκλιση, κλπ) όσο και την εκτίμηση της πιθανοτικής κατανομής της δομικής απόκρισης (τάση, μετατόπιση, κλπ.).Η προτεινόμενη προσέγγιση ενσωματώνεται σε δύο ευρέως γνωστά λογισμικά ανάλυσης πεπερασμένων στοιχείων, δηλαδή στο λογισμικό OpenSees χρησιμοποιώντας γλώσσα προγραμματισμού Tcl και στο λογισμικό ABAQUS χρησιμοποιώντας γλώσσα προγραμματισμού Python. Εν συνεχεία, το OpenSees χρησιμοποιείται για την ανάλυση δομών υπό σεισμική φόρτιση, όπου εκτελούνται ανελαστική στατική ανάλυση (pushover) και δυναμική ανάλυση χρησιμοποιώντας επεξεργασμένα επιταχυνσιογραφήματα από προηγούμενους σεισμούς και λαμβάνοντας υπόψη και επαναλαμβανόμενους σεισμούς. Επίσης, το ABAQUS χρησιμοποιείται για την ανάλυση δομών υπό μη γραμμική στατική φόρτωση, για συνδέσεις πλακών-υποστυλώματος από οπλισμένο σκυρόδεμα και για τοιχία πυρηνικών αντιδραστήρων από προεντεταμένο σκυρόδεμα. Αυτή η έρευνα δείχνει ότι η προτεινόμενη μέθοδος, η οποία βασίζεται σε μικρό αριθμό αναλύσεων πεπερασμένων στοιχείων, είναι αποτελεσματική, υπολογιστικά εφικτή και εύκολα εφαρμόσιμη, παρέχοντας μια εφικτή εναλλακτική λύση για την ανάλυση αξιοπιστίας (και ανάλυση ευαισθησίας) με την χρήση των πεπερασμένων στοιχείων για πραγματικά κατασκευαστικά προβλήματα (κατασκευές τρισδιάστατες και μεγάλης κλίμακας). Τα αποτελέσματα μιας τέτοιας εργασίας έχουν σημασία σε μελλοντικές μελέτες για την εκτίμηση της πιθανότητας της δομικής απόκρισης να υπερβεί ένα όριο ασφαλείας και για τον καθορισμό παραγόντων ασφάλειας που σχετίζονται με αποδεκτές πιθανότητες δομικών βλαβών.


Author(s):  
S. Azarm ◽  
S. M. Bhandarkar ◽  
A. J. Durelli

Abstract This paper describes and compares an experimental and a numerical method for shape optimization of continuum structures. The experimental method is based on a systematic application of photoelasticity. The numerical method is based on coupling of a finite element analysis and an optimizer. The paper demonstrate several means used for interactive communication between the designer and the (experimental or numerical) shape optimization process. The discussion is made taking shape optimization of a hole in a tall beam as an example.


2020 ◽  
Vol 26 ◽  
pp. 78
Author(s):  
Thirupathi Gudi ◽  
Ramesh Ch. Sau

We study an energy space-based approach for the Dirichlet boundary optimal control problem governed by the Laplace equation with control constraints. The optimality system results in a simplified Signorini type problem for control which is coupled with boundary value problems for state and costate variables. We propose a finite element based numerical method using the linear Lagrange finite element spaces with discrete control constraints at the Lagrange nodes. The analysis is presented in a combination for both the gradient and the L2 cost functional. A priori error estimates of optimal order in the energy norm is derived up to the regularity of the solution for both the cases. Theoretical results are illustrated by some numerical experiments.


1996 ◽  
Vol 3 (4) ◽  
pp. 259-268 ◽  
Author(s):  
M.S. Yao

The large number of unknown variables in a finite element idealization for dynamic structural analysis is represented by a very small number of generalized variables, each associating with a generalized Ritz vector known as a basis vector. The large system of equations of motion is thereby reduced to a very small set by this transformation and computational cost of the analysis can be greatly reduced. In this article nonlinear equations of motion and their transformation are formulated in detail. A convenient way of selection of the generalized basis vector and its limitations are described. Some illustrative examples are given to demonstrate the speed and validity of the method. The method, within its limitations, may be applied to dynamic problems where the response is global in nature with finite amplitude.


Author(s):  
Sedat Ozcanan ◽  
Ali Osman Atahan

For guardrail designers, it is essential to achieve a crashworthy and optimal system design. One of the most critical parameters for an optimal road restraint system is the post embedment depth or the post-to-soil interaction. This study aims to assess the optimum post embedment depth values of three different guardrail posts embedded in soil with varying density. Posts were subjected to dynamic impact loads in the field while a detailed finite element study was performed to construct accurate models for the post–soil interaction. It is well-known that experimental tests and simulations are costly and time-consuming. Therefore, to reduce the computational cost of optimization, radial basis function–based metamodeling methodology was employed to create surrogate models that were used to replace the expensive three-dimensional finite element models. In order to establish the radial basis function model, samples were derived using the full factorial design. Afterward, radial basis function–based metamodels were generated from the derived data and objective functions performed using finite element analysis. The accuracy of the metamodels were validated by k-fold cross-validation, then optimized using multi-objective genetic algorithm. After optimum embedment depths were obtained, finite element simulations of the results were compared with full-scale crash test results. In comparison with the actual post embedment depths, optimal post embedment depths provided significant economic advantages without compromising safety and crashworthiness. It is concluded that the optimum post embedment depths provide an economic advantage of up to 17.89%, 36.75%, and 43.09% for C, S, and H types of post, respectively, when compared to actual post embedment depths.


2019 ◽  
Vol 53 (4) ◽  
pp. 1191-1222 ◽  
Author(s):  
Seungil Kim

In this paper, we study finite element approximate solutions to the Helmholtz equation in waveguides by using a perfectly matched layer (PML). The PML is defined in terms of a piecewise linear coordinate stretching function with two parameters for absorbing propagating and evanescent components respectively, and truncated with a Neumann condition on an artificial boundary rather than a Dirichlet condition for cutoff modes that waveguides may allow. In the finite element analysis for the PML problem, we have to deal with two difficulties arising from the lack of full regularity of PML solutions and the anisotropic nature of the PML problem with, in particular, large PML damping parameters. Anisotropic finite element meshes in the PML regions depending on the damping parameters are used to handle anisotropy of the PML problem. As a main goal, we establish quasi-optimal a priori error estimates, that does not depend on anisotropy of the PML problem (when no cutoff mode is involved), including the exponentially convergent PML error with respect to the width and the strength of PML. The numerical experiments that confirm the convergence analysis will be presented.


2010 ◽  
Vol 670 ◽  
pp. 284-290 ◽  
Author(s):  
Themistoklis D. Kefalas ◽  
George Loizos ◽  
Antonios G. Kladas

Even though, the flux distribution at joints of stacked type transformer cores has been investigated thoroughly many issues remain unclear in the case of wound transformer cores. The paper addresses this lack of information by longitudinal and normal flux measurements at step-lap joints of Si-Fe wound cores. Flux measurements are verified by an original finite element analysis where the necessary excitation is performed by means of a pseudo-source. The advantage of the proposed technique is the accurate estimation of the flux distribution at step-lap joints, with a two dimensional model of simple geometry and low computational cost, by using any commercial finite element code.


2004 ◽  
Vol 127 (2) ◽  
pp. 135-140 ◽  
Author(s):  
Tong Hong Wang ◽  
Yi-Shao Lai

In a finite element analysis, when localized behavior of a large model is of particular concern, generally one would refine the mesh until it captures the local solution adequately. Submodeling is an alternative way for solving this kind of problem. It provides a relatively accurate solution at a modest computational cost. For a valid submodeling analysis, the boundaries of the submodel should be sufficiently far away from local features so that St. Venant’s principle holds. Moreover, special treatments are required for solving problems that involve path-dependent characteristics. This paper presents a general procedure to perform submodeling analyses for path-dependent thermomechanical problems without a priori assumptions on the structural response. The procedure was benchmarked using a bimaterial strip and demonstrated through analyses on a bump chip carrier package assembly. The procedure is conducive to the numerical assessment of fatigue lives of electronic packages.


Author(s):  
Goldy Kumar ◽  
Vadim Shapiro

Laminate composites are widely used in automotive, aerospace, medical, and increasingly in consumer industries, due to their reduced weight, superior structural properties and cost-effectiveness. However, structural analysis of complex laminate structures remains challenging. 2D finite element methods based on plate and shell theories may be accurate and efficient, but they generally do not apply to the whole structure, and require identification and preprocessing (dimensional reduction) of the regions where the underlying assumptions hold. Differences in and limitations of theories for thin/thick plates and shells further complicate modeling and simulation of composites. Fully automated structural analysis using 3D elements with sufficiently high order basis functions is possible in principle, but is rarely practiced due to the significant increase in computational integration cost in the presence of a large number of laminate plies. We propose to replace the actual layup of the laminate structure by a simplified material model, allowing for a substantial reduction of the computational cost of 3D FEA. The reduced model, under the usual assumptions made in lamination theory, has the same constitutive relationship as the corresponding 2D plate model of the original laminate, but requires only a small fraction of computational integration costs in 3D FEA. We describe implementation of 3D FEA using the reduced material model in a meshfree system using second order B-spline basis functions. Finally, we demonstrate its validity by showing agreement between computed and known results for standard problems.


2013 ◽  
Vol 740 ◽  
pp. 319-322 ◽  
Author(s):  
Young Choon Lee ◽  
Nam Jin Jeon ◽  
Cheol Kim ◽  
Seo Yeon Ahn ◽  
Myung Jae Cho

Finite element analysis was accomplished for a steering knuckle component of a small bus to see whether the static and fatigue strength requirements were satisfied or not. The knuckle was modeled with ANSYS 10-node quadratic elements. The cyclic fatigue load was applied and Soderberg criteria were applied to check the fatigue life. The knuckle structure has an infinite life (10-6 cycle) judging from the fatigue analyses. Shape optimization based on the gradient based method has been performed in order to find out the knuckle shape that has a minimum weight and satisfies the static and fatigue strength requirements. As a result of shape optimization, the weight of the steering knuckle was reduced 8%.


Sign in / Sign up

Export Citation Format

Share Document