ISEL - Matemática
Permanent URI for this community
Browse
Browsing ISEL - Matemática by Issue Date
Now showing 1 - 10 of 262
Results Per Page
Sort Options
- Multi-objective optimization of structures topology by genetic algorithmsPublication . Madeira, JFA; Rodrigues, H.; Pina, HeitorThis work develops a computational model for topology optimization of linear elastic structures for situations where more than one objective function is required, each one of them with a different optimal solution. The method is thus developed for multi-objective optimization problems and is based on Genetic Algorithms. Its purpose is to evolve an evenly distributed group of solutions (population) to obtain the optimum Pareto set for the given problem. To reduce computational effort, optimal solutions of each of the single-objective problems are introduced in the initial population. Two numerical examples are presented and discussed to assess the method.
- Multiobjective topology optimization of structures using genetic algorithms with chromosome repairingPublication . Madeira, JFA; Rodrigues, H. C.; Pina, H.In this work, a genetic algorithm (GA) for multiobjective topology optimization of linear elastic structures is developed. Its purpose is to evolve an evenly distributed group of solutions to determine the optimum Pareto set for a given problem. The GA determines a set of solutions to be sorted by its domination properties and a filter is defined to retain the Pareto solutions. As an equality constraint on volume has to be enforced, all chromosomes used in the genetic GA must generate individuals with the same volume value; in the coding adopted, this means that they must preserve the same number of “ones” and, implicitly, the same number of “zeros” along the evolutionary process. It is thus necessary: (1) to define chromosomes satisfying this propriety and (2) to create corresponding crossover and mutation operators which preserve volume. Optimal solutions of each of the single-objective problems are introduced in the initial population to reduce computational effort and a repairing mechanism is developed to increase the number of admissible structures in the populations. Also, as the work of the external loads can be calculated independently for each individual, parallel processing was used in its evaluation. Numerical applications involving two and three objective functions in 2D and two objective functions in3Dare employed as tests for the computational model developed. Moreover, results obtained with and without chromosome repairing are compared.
- Dynamical behaviour on the parameter space: new populational growth models proportional to beta densitiesPublication . Aleixo, Sandra; Rocha, J. Leonel; Pestana, Dinis D.We present new populational growth models, generalized logistic models which are proportional to beta densities with shape parameters p and 2, where p > 1, with Malthusian parameter r. The complex dynamical behaviour of these models is investigated in the parameter space (r, p), in terms of topological entropy, using explicit methods, when the Malthusian parameter r increases. This parameter space is split into different regions, according to the chaotic behaviour of the models.
- Measuring and Controlling the Chaotic Motion of ProfitsPublication . Januário, Cristina; Grácio, Clara; Mendes, Diana A.; Duarte, JorgeThe study of economic systems has generated deep interest in exploring the complexity of chaotic motions in economy. Due to important developments in nonlinear dynamics, the last two decades have witnessed strong revival of interest in nonlinear endogenous business chaotic models. The inability to predict the behavior of dynamical systems in the presence of chaos suggests the application of chaos control methods, when we are more interested in obtaining regular behavior. In the present article, we study a specific economic model from the literature. More precisely, a system of three ordinary differential equations gather the variables of profits, reinvestments and financial flow of borrowings in the structure of a firm. Firstly, using results of symbolic dynamics, we characterize the topological entropy and the parameter space ordering of kneading sequences, associated with one-dimensional maps that reproduce significant aspects of the model dynamics. The analysis of the variation of this numerical invariant, in some realistic system parameter region, allows us to quantify and to distinguish different chaotic regimes. Finally, we show that complicated behavior arising from the chaotic firm model can be controlled without changing its original properties and the dynamics can be turned into the desired attracting time periodic motion (a stable steady state or into a regular cycle). The orbit stabilization is illustrated by the application of a feedback control technique initially developed by Romeiras et al. [1992]. This work provides another illustration of how our understanding of economic models can be enhanced by the theoretical and numerical investigation of nonlinear dynamical systems modeled by ordinary differential equations.
- Toplological Entropy in the Syncronization of Piecewise Linear and Monotone Maps. Coupled Duffing OscillatorsPublication . Caneco, Acilina; Rocha, J. Leonel; Grácio, ClaraIn this paper is presented a relationship between the synchronization and the topological entropy. We obtain the values for the coupling parameter, in terms of the topological entropy, to achieve synchronization of two unidirectional and bidirectional coupled piecewise linear maps. In addition, we prove a result that relates the synchronizability of two m-modal maps with the synchronizability of two conjugated piecewise linear maps. An application to the unidirectional and bidirectional coupled identical chaotic Duffing equations is given. We discuss the complete synchronization of two identical double-well Duffing oscillators, from the point of view of symbolic dynamics. Working with Poincare cross-sections and the return maps associated, the synchronization of the two oscillators, in terms of the coupling strength, is characterized.
- Chaos and crises in a model for cooperative hunting: A symbolic dynamics approachPublication . Duarte, Jorge; Januário, Cristina; Martins, Nuno; Sardanyes, JosepIn this work we investigate the population dynamics of cooperative hunting extending the McCann and Yodzis model for a three-species food chain system with a predator, a prey, and a resource species. The new model considers that a given fraction sigma of predators cooperates in prey's hunting, while the rest of the population 1-sigma hunts without cooperation. We use the theory of symbolic dynamics to study the topological entropy and the parameter space ordering of the kneading sequences associated with one-dimensional maps that reproduce significant aspects of the dynamics of the species under several degrees of cooperative hunting. Our model also allows us to investigate the so-called deterministic extinction via chaotic crisis and transient chaos in the framework of cooperative hunting. The symbolic sequences allow us to identify a critical boundary in the parameter spaces (K, C-0) and (K, sigma) which separates two scenarios: (i) all-species coexistence and (ii) predator's extinction via chaotic crisis. We show that the crisis value of the carrying capacity K-c decreases at increasing sigma, indicating that predator's populations with high degree of cooperative hunting are more sensitive to the chaotic crises. We also show that the control method of Dhamala and Lai [Phys. Rev. E 59, 1646 (1999)] can sustain the chaotic behavior after the crisis for systems with cooperative hunting. We finally analyze and quantify the inner structure of the target regions obtained with this control method for wider parameter values beyond the crisis, showing a power law dependence of the extinction transients on such critical parameters.
- Surgical correction of scoliosis: numerical analysis and optimization of the procedurePublication . Madeira, JFA; Pina, H. L.; Pires, E. B.; Monteiro, JoaquimA previously developed model is used to numerically simulate real clinical cases of the surgical correction of scoliosis. This model consists of one-dimensional finite elements with spatial deformation in which (i) the column is represented by its axis; (ii) the vertebrae are assumed to be rigid; and (iii) the deformability of the column is concentrated in springs that connect the successive rigid elements. The metallic rods used for the surgical correction are modeled by beam elements with linear elastic behavior. To obtain the forces at the connections between the metallic rods and the vertebrae geometrically, non-linear finite element analyses are performed. The tightening sequence determines the magnitude of the forces applied to the patient column, and it is desirable to keep those forces as small as possible. In this study, a Genetic Algorithm optimization is applied to this model in order to determine the sequence that minimizes the corrective forces applied during the surgery. This amounts to find the optimal permutation of integers 1,..., n, n being the number of vertebrae involved. As such, we are faced with a combinatorial optimization problem isomorph to the Traveling Salesman Problem. The fitness evaluation requires one computing intensive Finite Element Analysis per candidate solution and, thus, a parallel implementation of the Genetic Algorithm is developed. Copyright q 2010 John Wiley & Sons, Ltd.
- GA topology optimization using random keys for tree encoding of structuresPublication . Madeira, JFA; Pina, H. L.; Rodrigues, H. C.Topology optimization consists in finding the spatial distribution of a given total volume of material for the resulting structure to have some optimal property, for instance, maximization of structural stiffness or maximization of the fundamental eigenfrequency. In this paper a Genetic Algorithm (GA) employing a representation method based on trees is developed to generate initial feasible individuals that remain feasible upon crossover and mutation and as such do not require any repairing operator to ensure feasibility. Several application examples are studied involving the topology optimization of structures where the objective functions is the maximization of the stiffness and the maximization of the first and the second eigenfrequencies of a plate, all cases having a prescribed material volume constraint.
- Completion problems for real matrices, IIPublication . Matos, Isabel T.; Silva, Fernando C.Thirty years ago, G.N. de Oliveira has proposed the following completion problems: Describe the possible characteristic polynomials of [C-ij], i,j is an element of {1, 2}, where C-1,C-1 and C-2,C-2 are square submatrices, when some of the blocks C-ij are fixed and the others vary. Several of these problems remain unsolved. This paper gives the solution, over the field of real numbers, of Oliveira's problem where the blocks C-1,C-1, C-2,C-2 are fixed and the others vary.
- Solutions of second-order and fourth-order ODEs on the half-linePublication . Enguiça, Ricardo Roque; Gavioli, Andrea; Sanchez, LuisWe start by studying the existence of positive solutions for the differential equation u '' = a(x)u - g(u), with u ''(0) = u(+infinity) = 0, where a is a positive function, and g is a power or a bounded function. In other words, we are concerned with even positive homoclinics of the differential equation. The main motivation is to check that some well-known results concerning the existence of homoclinics for the autonomous case (where a is constant) are also true for the non-autonomous equation. This also motivates us to study the analogous fourth-order boundary value problem {u((4)) - cu '' + a(x)u = vertical bar u vertical bar(p-1)u u'(0) = u'''(0) = 0, u(+infinity) = u'(+infinity) = 0 for which we also find nontrivial (and, in some instances, positive) solutions.