槽形托輥帶式輸送機(jī)設(shè)計(jì)【30米】【說(shuō)明書(shū)+CAD】
槽形托輥帶式輸送機(jī)設(shè)計(jì)【30米】【說(shuō)明書(shū)+CAD】,30米,說(shuō)明書(shū)+CAD,槽形托輥帶式輸送機(jī)設(shè)計(jì)【30米】【說(shuō)明書(shū)+CAD】,槽形托輥帶式,輸送,設(shè)計(jì),30,說(shuō)明書(shū),仿單,cad
ORIGINAL ARTICLE Pseudo-constructal theory for shape optimization of mechanical structures Jean Luc Marcelin Received: 10 January 2007 /Accepted: 1 May 2007 /Published online: 25 May 2007 # Springer-Verlag London Limited 2007 Abstract This work gives some applications of a pseudo- constructal technique for shape optimization of mechanical structures. In the pseudo-constructal theory developed in this paper, the main objective of optimization is only the minimization of total potential energy. The other objectives usually used in mechanical structures optimization are treated like limitations or optimization constraints. Two applications are presented; the first one deals with the optimization of the shape of a drop of water by using a genetic algorithm with the pseudo-constructal technique, and the second one deals with the optimization of the shape of a hydraulic hammers rear bearing. Keywords Shapeoptimization . Constructal . Geneticalgorithms 1 Introduction This paper introduces a pseudo-constructal approach to shape optimization based on the minimization of the total potential energy. We are going to show that minimizing the total potential energy of a structure to find the optimal shape might be a good idea in some cases. The reference to the constructal theory can be justified in some way for the following reasons. According to Bejan 1, shape and structure spring from the struggle for better performance in both engineering and nature; the objective and constraints principle used in engineering is the same mechanism from which the geometry in natural flow systems emerges. Bejan 1 starts with the design and optimization of engineering systems and discovers a deterministic principle for the generation of geometric form in natural systems. This observation is the basis of the new constructal theory. Optimal distribution of imperfection is destined to remain imperfect. The system works best when its imperfections are spread around so that more and more internal points are stressed as much as the hardest working parts. Seemingly universal geometric forms unite the flow systems of engineering and nature. Bejan 1 advances a new theory in which he unabashedly hints that his law is in the same league as the second law of thermodynamics, because a simple law is purported to predict the geometric form of anything alive on earth. Many applications of the constructal theory were developed in fluids mechanics, in particular for the optimization of flows 210. On the other hand, there exists, to our knowledge, little examples of applications in solids or structures mechanics. So we have at least half of the references to papers in fluid dynamics (most of the same author), because the constructal method was developed first by the same author, Adrian Bejan, with only references to papers in fluid dynamics. The constructal theory rests on the assumption that all creations of nature are overall optimal compared to the laws which control the evolution and the adaptation of the natural systems. The constructal principle consists of distributing the imperfections as well as possible, starting from the smallest scales to the largest. The constructal theory works with the total macroscopic structure starting from the assembly of elementary struc- tures, by complying with the natural rules of optimal distribution of the imperfections. The objective is the research of lower cost. Int J Adv Manuf Technol (2008) 38:16 DOI 10.1007/s00170-007-1080-2 J. L. Marcelin (*) Laboratorie Sols Solides Structures 3S, UMR CNRS C5521, Domaine Universitaire, BP n53, 38041 Grenoble Cedex 9, France e-mail: Jean-Luc.Marcelinujf-grenoble.fr However, a global and macroscopic solution for the optimization of mechanical structures having least cost as the objective can be very close to the constructal theory, from where the term pseudo-constructal comes. The constructal theory is a predictive theory, with only one single principle of optimization from which all rises. The same applies to the pseudo-constructal step which is the subject of this article. The single principle of optimiza- tion of the pseudo-constructal theory is the minimization of total potential energy. Moreover, in our examples presented hereafter, the pseudo-constructal principle will be associated with a genetic algorithm, with the result that our optimization will be very close to the natural laws. The objective of this paper is thus to show how the pseudo-constructal step can apply to the mechanics of the structures, and in particular to the shape optimization of mechanical structures. The basic idea is very simple: a mechanical structure in a balanced state corresponds to a minimal total potential energy. In the same way, an optimal mechanical structure must also correspond to a minimal total potential energy, and it is this objective which must intervene first over all the others. It is this idea which will be developed in this article. Two examples will be presented thereafter. The idea to minimize total potential energy in order to optimize a mechanical structure is not brand new. Many papers already deal with this problem. What is new, is to make this approach systematic. The only objective of optimization becomes the minimization of energy. In Gosling 11, a simple method is proposed for the difficult case of form-finding of cablenet and membrane structures. This method is based upon basic energy concepts. A truncated strain expression is used to define the total potential energy. The final energy form is minimized using the Powell algorithm. In Kanno and Ohsaki 12, the minimum principle of complementary energy is established for cable networks involving only stress components as variables in geometrically nonlinear elasticity. In order to show the strong duality between the minimization problems of total potential energy and complementary energy, the convex formulations of these problems are investigated, which can be embedded into a primal-dual pair of second-order programming problems. In Taroco 13, shape sensitivity analysis of an elastic solid in equilibrium is presented. The domain and boundary integral expressions of the first and second-order shape derivatives of the total potential energy are established. In Warner 14, an optimal design problem is solved for an elastic rod hanging under its own weight. The distribution of the cross- sectional area that minimizes the total potential energy stored in an equilibrium state is found. The companion problem of the design that stores the maximum potential energy under the same constraint conditions is also solved. In Ventura 15, the problem of boundary conditions enforcement in meshless methods is solved. In Ventura 15, the moving least-squares approximation is introduced in the total potential energy functional for the elastic solid problem and an augmented Lagrangian term is added to satisfy essential boundary conditions. The principle of minimization of total potential energy is in addition at the base of the general finite elements formulation, with an aim of finding the unknown optimal nodal factors 16. 2 The methods used In the pseudo-constructal theory developed in this paper, the main objective of optimization is only the minimization of total potential energy. The other objectives usually used in mechanical structures optimization are treated here like limitations or optimization constraints. For example, one may have limitations on the weight, or to not exceed the value of a stress. The idea which will be developed in this paper is thus very simple. A mechanical structure is described by two types of parameters: variables known as discretization variables (for example, degrees of freedom in displacement for finite elements method), and geometrical variables of design (for example parameters which make it possible to describe the mechanical structure shape). Total potential energy depends on an implicit or explicit way of determin- ing discretization and design variables at the same time. One thus will carry out a double optimization of the mechanical structure, compared to the discretization and design variables; the objective being to minimize total potential energy overall. Clearly, the problem of optimiza- tion of a mechanical structure will be addressed by the following approach: Objective: to minimize total potential energy Variables of optimization: concurrently determining discretization variables (in the case of a traditional use of the finite element method in mechanics of struc- tures), and design variables describing the shape of the structure Optimization limitations: Weight or volume Displacements or strains Stresses Frequencies The problem of optimization of a mechanical structure will be solved in the following way, while reiterating on 2 Int J Adv Manuf Technol (2008) 38:16 these stages, if needed (according to the nature of the problem): Stage 1 Minimization of the total potential energy of the mechanical structure compared to the only dis- cretization variables of the structure (degrees of freedom in finite elements). It acts here as an optimization without optimization limitations. The only limitations at this stage are of purely mechanical origin, and relate to the boundary conditions and to the external efforts applied to the structure. In this stage 1, the design variables remain fixed, and one obtains the implicit or explicit expressions of the degrees of freedom according to the design variables (which can be the variables which make it possible to describe the shape, in the case of a shape optimization, for example). One will see in the examples of the following part that these expressions can be explicit or implicit and which is the suitable treatment following the cases. In the case of a finite elements method of calculation, this stage 1 is the basis of finite elements calculation to obtain the degrees of freedom of the mechanical structure. Indeed, in finite elements, displacements with the nodes of the mechanical structure mesh are obtained by minimization of total potential energy 16. Stage 2 The expressions of the degrees of freedom of the mechanical structure according to the design variables obtained previously are then injected into the total potential energy of the mechanical structure (one will see in the second example of thefollowingparthowonetreatsthecasewherethe degrees of freedom are implicit functions of the design variables). One then obtains an expression of the total potential energy which depends only on the design variables (in explicit or implicit form). Stage 3 One then carries out a second and new minimi- zation of the total potential energy obtained in the preceding form, but this time compared to the design variables while respecting the technolog- ical limitations or the optimization constraints of the problem. This method can be applied with more or less facility according to the nature of the problem. It is clear, for example, that if the discretization variables can be expressed in an explicit way according to the design variables, the setting in of stages 2 to 3 is immediate, and without iterations. If the discretization variables cannot be expressed in an explicit way according to the design variables, or if the topology of the structure is not fixed, or if the behavior is not linear, it will be necessary to proceed by successive iterations on stages 1 to 3. It is the case of the examples presented in the following part, and one will see on this occasion which type of strategy one can adopt for these iterations. To summarize, in the pseudo-constructal step, the main objective is only the minimization of total potential energy, the other possible objectives are treated like limitations or optimization constraints. The optimization method used for our examples is GA (genetic algorithm), as described in 17. Examples with similar instructional value can also be found in many books, e.g. in 18. This evolutionary method is very convenient for our pseudo-constructal method. The author has worked extensively in GAs and published in some reputed journals on this topic 1931. As the topic of GAs is still relatively new in the structural mechanics commu- nity, we provide here some details of exactly what is used in this GA. A multiple point crossover is used rather than a single point crossover. The selection scheme used at each generation is entirely stochastic. For our examples, the number of generations is equal to that used for conver- gence. The results provided for our examples were consistently reproduced by using different seeds in the GA. It has been proved that a rather standard genetic algorithm is sufficient for our examples. 3 Examples Even though potential energy may be a good measure for some optimizations, potential energy is not what gives the shape to a water droplet, nor defines the optimal shape for a hammer, which is why potential energy is not the only objective; but the optimization problem is a multiobjective one and the objective functions for the two examples are then clearly formulated. 3.1 Example 1: optimization of the shape of a drop of water The first test example is the optimization of the shape of a drop of water (Fig. 1). This problem is equivalent to an equal resistance tank calculated by the membrane theory. The objective is to see if the pseudo-constructal theory gives the natures optimum design. 3.1.1 The methods used The geometry of the drop of water is defined by the generating line of a thin axisymmetric shell. This line is described by successive straight or circular segments described in a given sense and defined by input data of master point coordinates and radius values. The initial data are a set of nodal points connected by straight segments. Each nodal point is identified by its two cylindrical Int J Adv Manuf Technol (2008) 38:16 3 coordinates (r, z), and a real R which represents the radius of the circle tangent to the two straight segments intersect- ing at the point. The other computer calculations give the coordinates of any boundary point and especially the tangent points necessary to define the circular arc lengths. The design of the drop of water is described by three arcs of circles as indicated in Fig. 1. Analysis is performed by the finite element method with three-node parabolic elements using the classical Love- Kirchoff shell theory. An automatic mesh generator creates the finite element mesh of each straight or circular segment considered as a macro finite element. The objective is to obtain a shape for the drop of water giving rise to a minimum total potential energy (which is the main objective) and an equal resistance tank (which is the only constraint or limitation of the problem). In fact, for the drop of water problem, the goal is a multi- objective one, the two objectives ( f 1 =minimum total potential energy and f 2 =equal resistance) are combined in a multi-objective: f=f 1 +f 2 . The constraint or limitation of the problem is taken into account by a penalization of the total potential energy as indicated in Marcelin et al. 19. 3.1.2 The results The design of the drop of water is described by three arcs of a circle (Fig. 1). Their centers and radius are the design variables. So, there are nine design variables: r1, z1, R1 for circle 1; r2, z2, R2 for circle 2; and r3, z3, R3 for circle 3. In the genetic algorithm, each of these design variables is coded by three binary digits. The tables of coding-decoding will be the following: For r1: For z1: For R1: For r2: For z2: For R2: For r3: For z3: For R3: All these binary digits are put end to end to form a chromosome length of 27 binary digits. GA is run for a population of 30 individuals, a number of generations of 50, a probability of crossing of 0.8, and a probability of mutation of 0.1. The optimal solution corresponds to the chromosome 100 100 011 011 010 011 100 011 101 which gives the solution of Fig. 1, for which: r1=18, z1=17,and R1=0.065 r2=13.75, z2=12.2 and R2=7.7 r3=4.1, z3=21.4 and R3=21 It is very close to the natures optimal solution for the shape of a drop of water. The model of the water drop modelled by three arcs of a circle is imperfect. However, the constructal theory optimizes the imperfections, and 105 20 25 r 5 10 20 z 3 2 1 15 15 Fig. 1 Optimization of the shape of a drop of water 000 001 010 011 100 101 110 111 16 16.5 17 17.5 18 18.5 19 19.5 000 001 010 011 100 101 110 111 15 15.5 16 16.5 17 17.5 18 18.5 000 001 010 011 100 101 110 111 0.050 0.055 0.060 0.065 0.070 0.075 0.080 0.085 000 001 010 011 100 101 110 111 13 13.25 13.5 13.75 14 14.25 14.5 14.75 000 001 010 011 100 101 110 111 12 12.1 12.2 12.3 12.4 12.5 12.6 12.7 000 001 010 011 100 101 110 111 7.4 7.5 7.6 7.7 7.8 7.9 8 8.1 000 001 010 011 100 101 110 111 3.7 3.8 3.9 4 4.1 4.2 4.3 4.4 000 001 010 011 100 101 110 111 21.1 21.2 21.3 21.4 21.5 21.6 21.7 21.8 000 001 010 011 100 101 110 111 18.5 19 19.5 20 20.5 21 21.5 22 4 Int J Adv Manuf Technol (2008) 38:16 finds the nearest solution to that of nature. So, the constructal principle consists of distributing the imperfec- tions as well as possible. 3.2 Example 2: optimization of the shape of an axisymmetric structure In this part, the very localized optimization of the rear bearing of a hydraulic hammer is presented. The bearing in question (Fig. 2) breaks after relatively few cycles of operation. For axisymmetric structures, analysis is performed by the finite element method in which the special character of a GA optimization process has been considered to ease the calculations and to save computer time. First, because just a few parts of the structure must often be modified, the substructure concept is used to separate the “fixed” and the “mobile” parts. The fixed parts are calculated twice: once at the beginning and also at the end of the optimization process. Only the reduced stiffness matrices of these substructures are added to the matrices of the mobile parts. Related to this division, an automatic generator creates the finite element mesh of each substructure considered as a macro finite element. These macro elements are either triangular (six nodes) or quadrilateral (eight nodes). Following a well-known technique, the same subdivision is used in the parent space to obtain the mesh itself, which is obviously made out of the same types of elements. During the optimization process this mesh is controlled and a new discretization can be chosen if necessary. To summarize, the optimization problem is the following: Objective function Minimization of the total potential energy. It is important to note that another important objective (the minimization of the maximum value of the Von Mises equivalent stress along the mobile contour) is taken here as a constraint of the problem. This second objective is necessary to achieve the minimization of the rear bearing of the hydraulic hammer . Design variables The design variables are radius r and width X near the radius (Fig. 2). Constraints The side constraints are established in such a way that only small changes in geometry are allowed. They take into account the technological constraints. They are included in the coding of the design variables. Another important constraint is that the maximum value of the Von Mises equivalent stress along the mobile contour must not exceed a certain value. The constraints are taken into account by a penalization of the total potential energy as indicated in 19. The tables of coding-decoding are the following: For r: For X: All these binary digits are put end to end to form a chromosome length of eight binary digits. The GA is run for a population of 12 individuals, a number of generations of 30, a probability of crossing of 0.5, and a probability of mutation of 0.06. The optimal solution corresponds to the chromosome 1101 1000 which gives the solution of Fig. 2, for which: r=1.95, X=6.0 The automatic optimization of the shape of this product has,
收藏