Analysis of the static and dynamic behaviour of hydraulic fills

En la actualidad existe un gran conocimiento en la caracterización de rellenos hidráulicos, tanto en su caracterización estática, como dinámica. Sin embargo, son escasos en la literatura estudios más generales y globales de estos materiales, muy relacionados con sus usos y principales problemáticas en obras portuarias y mineras. Los procedimientos semi‐empíricos para la evaluación del efecto silo en las celdas de cajones portuarios, así como para el potencial de licuefacción de estos suelos durantes cargas instantáneas y terremotos, se basan en estudios donde la influencia de los parámetros que los rigen no se conocen en gran medida, dando lugar a resultados con considerable dispersión. Este es el caso, por ejemplo, de los daños notificados por el grupo de investigación del Puerto de Barcelona, la rotura de los cajones portuarios en el Puerto de Barcelona en 2007. Por estos motivos y otros, se ha decidido desarrollar un análisis para la evaluación de estos problemas mediante la propuesta de una metodología teórico‐numérica y empírica. El enfoque teórico‐numérico desarrollado en el presente estudio se centra en la determinación del marco teórico y las herramientas numéricas capaces de solventar los retos que presentan estos problemas. La complejidad del problema procede de varios aspectos fundamentales: el comportamiento no lineal de los suelos poco confinados o flojos en procesos de consolidación por preso propio; su alto potencial de licuefacción; la caracterización hidromecánica de los contactos entre estructuras y suelo (camino preferencial para el flujo de agua y consolidación lateral); el punto de partida de los problemas con un estado de tensiones efectivas prácticamente nulo. En cuanto al enfoque experimental, se ha propuesto una metodología de laboratorio muy sencilla para la caracterización hidromecánica del suelo y las interfaces, sin la necesidad de usar complejos aparatos de laboratorio o procedimientos excesivamente complicados. Este trabajo incluye por tanto un breve repaso a los aspectos relacionados con la ejecución de los rellenos hidráulicos, sus usos principales y los fenómenos relacionados, con el fin de establecer un punto de partida para el presente estudio. Este repaso abarca desde la evolución de las ecuaciones de consolidación tradicionales (Terzaghi, 1943), (Gibson, English & Hussey, 1967) y las metodologías de cálculo (Townsend & McVay, 1990) (Fredlund, Donaldson and Gitirana, 2009) hasta las contribuciones en relación al efecto silo (Ranssen, 1985) (Ravenet, 1977) y sobre el fenómeno de la licuefacción (Casagrande, 1936) (Castro, 1969) (Been & Jefferies, 1985) (Pastor & Zienkiewicz, 1986). Con motivo de este estudio se ha desarrollado exclusivamente un código basado en el método de los elementos finitos (MEF) empleando el programa MATLAB. Para ello, se ha esablecido un marco teórico (Biot, 1941) (Zienkiewicz & Shiomi, 1984) (Segura & Caron, 2004) y numérico (Zienkiewicz & Taylor, 1989) (Huerta & Rodríguez, 1992) (Segura & Carol, 2008) para resolver problemas de consolidación multidimensional con condiciones de contorno friccionales, y los correspondientes modelos constitutivos (Pastor & Zienkiewicz, 1986) (Fiu & Liu, 2011). Asimismo, se ha desarrollado una metodología experimental a través de una serie de ensayos de laboratorio para la calibración de los modelos constitutivos y de la caracterización de parámetros índice y de flujo (Castro, 1969) (Bahda 1997) (Been & Jefferies, 2006). Para ello se han empleado arenas de Hostun como material (relleno hidráulico) de referencia. Como principal aportación se incluyen una serie de nuevos ensayos de corte directo para la caracterización hidromecánica de la interfaz suelo – estructura de hormigón, para diferentes tipos de encofrados y rugosidades. Finalmente, se han diseñado una serie de algoritmos específicos para la resolución del set de ecuaciones diferenciales de gobierno que definen este problema. Estos algoritmos son de gran importancia en este problema para tratar el procesamiento transitorio de la consolidación de los rellenos hidráulicos, y de otros efectos relacionados con su implementación en celdas de cajones, como el efecto silo y la licuefacciones autoinducida. Para ello, se ha establecido un modelo 2D axisimétrico, con formulación acoplada u‐p para elementos continuos y elementos interfaz (de espesor cero), que tratan de simular las condiciones de estos rellenos hidráulicos cuando se colocan en las celdas portuarias. Este caso de estudio hace referencia clara a materiales granulares en estado inicial muy suelto y con escasas tensiones efectivas, es decir, con prácticamente todas las sobrepresiones ocasionadas por el proceso de autoconsolidación (por peso propio). Por todo ello se requiere de algoritmos numéricos específicos, así como de modelos constitutivos particulares, para los elementos del continuo y para los elementos interfaz. En el caso de la simulación de diferentes procedimientos de puesta en obra de los rellenos se ha requerido la modificacion de los algoritmos empleados para poder así representar numéricamente la puesta en obra de estos materiales, además de poder realizar una comparativa de los resultados para los distintos procedimientos. La constante actualización de los parámetros del suelo, hace también de este algoritmo una potente herramienta que permite establecer un interesante juego de perfiles de variables, tales como la densidad, el índice de huecos, la fracción de sólidos, el exceso de presiones, y tensiones y deformaciones. En definitiva, el modelo otorga un mejor entendimiento del efecto silo, término comúnmente usado para definir el fenómeno transitorio del gradiente de presiones laterales en las estructuras de contención en forma de silo. Finalmente se incluyen una serie de comparativas entre los resultados del modelo y de diferentes estudios de la literatura técnica, tanto para el fenómeno de las consolidaciones por preso propio (Fredlund, Donaldson & Gitirana, 2009) como para el estudio del efecto silo (Puertos del Estado, 2006, EuroCódigo (2006), Japan Tech, Stands. (2009), etc.). Para concluir, se propone el diseño de un prototipo de columna de decantación con paredes friccionales, como principal propuesta de futura línea de investigación. Wide research is nowadays available on the characterization of hydraulic fills in terms of either static or dynamic behavior. However, reported comprehensive analyses of these soils when meant for port or mining works are scarce. Moreover, the semi‐empirical procedures for assessing the silo effect on cells in floating caissons, and the liquefaction potential of these soils during sudden loads or earthquakes are based on studies where the underlying influence parameters are not well known, yielding results with significant scatter. This is the case, for instance, of hazards reported by the Barcelona Liquefaction working group, with the failure of harbor walls in 2007. By virtue of this, a complex approach has been undertaken to evaluate the problem by a proposal of numerical and laboratory methodology. Within a theoretical and numerical scope, the study is focused on the numerical tools capable to face the different challenges of this problem. The complexity is manifold; the highly non‐linear behavior of consolidating soft soils; their potentially liquefactable nature, the significance of the hydromechanics of the soil‐structure contact, the discontinuities as preferential paths for water flow, setting “negligible” effective stresses as initial conditions. Within an experimental scope, a straightforward laboratory methodology is introduced for the hydromechanical characterization of the soil and the interface without the need of complex laboratory devices or cumbersome procedures. Therefore, this study includes a brief overview of the hydraulic filling execution, main uses (land reclamation, filled cells, tailing dams, etc.) and the underlying phenomena (self‐weight consolidation, silo effect, liquefaction, etc.). It comprises from the evolution of the traditional consolidation equations (Terzaghi, 1943), (Gibson, English, & Hussey, 1967) and solving methodologies (Townsend & McVay, 1990) (Fredlund, Donaldson and Gitirana, 2009) to the contributions in terms of silo effect (Ranssen, 1895) (Ravenet, 1977) and liquefaction phenomena (Casagrande, 1936) (Castro, 1969) (Been & Jefferies, 1985) (Pastor & Zienkiewicz, 1986). The novelty of the study lies on the development of a Finite Element Method (FEM) code, exclusively formulated for this problem. Subsequently, a theoretical (Biot, 1941) (Zienkiewicz and Shiomi, 1984) (Segura and Carol, 2004) and numerical approach (Zienkiewicz and Taylor, 1989) (Huerta, A. & Rodriguez, A., 1992) (Segura, J.M. & Carol, I., 2008) is introduced for multidimensional consolidation problems with frictional contacts and the corresponding constitutive models (Pastor & Zienkiewicz, 1986) (Fu & Liu, 2011). An experimental methodology is presented for the laboratory test and material characterization (Castro 1969) (Bahda 1997) (Been & Jefferies 2006) using Hostun sands as reference hydraulic fill. A series of singular interaction shear tests for the interface calibration is included. Finally, a specific model algorithm for the solution of the set of differential equations governing the problem is presented. The process of consolidation and settlements involves a comprehensive simulation of the transient process of decantation and the build‐up of the silo effect in cells and certain phenomena related to self‐compaction and liquefaction. For this, an implementation of a 2D axi‐syimmetric coupled model with continuum and interface elements, aimed at simulating conditions and self‐weight consolidation of hydraulic fills once placed into floating caisson cells or close to retaining structures. This basically concerns a loose granular soil with a negligible initial effective stress level at the onset of the process. The implementation requires a specific numerical algorithm as well as specific constitutive models for both the continuum and the interface elements. The simulation of implementation procedures for the fills has required the modification of the algorithm so that a numerical representation of these procedures is carried out. A comparison of the results for the different procedures is interesting for the global analysis. Furthermore, the continuous updating of the model provides an insightful logging of variable profiles such as density, void ratio and solid fraction profiles, total and excess pore pressure, stresses and strains. This will lead to a better understanding of complex phenomena such as the transient gradient in lateral pressures due to silo effect in saturated soils. Interesting model and literature comparisons for the self‐weight consolidation (Fredlund, Donaldson, & Gitirana, 2009) and the silo effect results (Puertos del Estado (2006), EuroCode (2006), Japan Tech, Stands. (2009)). This study closes with the design of a decantation column prototype with frictional walls as the main future line of research.

Design of a Semi-analytical Method to Describe Plasma-Wall Interactions

Electric probes are objects immersed in the plasma with sharp boundaries which collect of emit charged particles. Consequently, the nearby plasma evolves under abrupt imposed and/or naturally emerging conditions. There could be localized currents, different time scales for plasma species evolution, charge separation and absorbing-emitting walls. The traditional numerical schemes based on differences often transform these disparate boundary conditions into computational singularities. This is the case of models using advection-diffusion differential equations with source-sink terms (also called Fokker-Planck equations). These equations are used in both, fluid and kinetic descriptions, to obtain the distribution functions or the density for each plasma species close to the boundaries. We present a resolution method grounded on an integral advancing scheme by using approximate Green's functions, also called short-time propagators. All the integrals, as a path integration process, are numerically calculated, what states a robust grid-free computational integral method, which is unconditionally stable for any time step. Hence, the sharp boundary conditions, as the current emission from a wall, can be treated during the short-time regime providing solutions that works as if they were known for each time step analytically. The form of the propagator (typically a multivariate Gaussian) is not unique and it can be adjusted during the advancing scheme to preserve the conserved quantities of the problem. The effects of the electric or magnetic fields can be incorporated into the iterative algorithm. The method allows smooth transitions of the evolving solutions even when abrupt discontinuities are present. In this work it is proposed a procedure to incorporate, for the very first time, the boundary conditions in the numerical integral scheme. This numerical scheme is applied to model the plasma bulk interaction with a charge-emitting electrode, dealing with fluid diffusion equations combined with Poisson equation self-consistently. It has been checked the stability of this computational method under any number of iterations, even for advancing in time electrons and ions having different time scales. This work establishes the basis to deal in future work with problems related to plasma thrusters or emissive probes in electromagnetic fields.

Nonlinear analysis of a simple model of temperature evolution in a satellite

We analyze a simple model of the heat transfer to and from a small satellite orbiting round a solar system planet. Our approach considers the satellite isothermal, with external heat input from the environment and from internal energy dissipation, and output to the environment as black-body radiation. The resulting nonlinear ordinary differential equation for the satellite’s temperature is analyzed by qualitative, perturbation and numerical methods, which prove that the temperature approaches a periodic pattern (attracting limit cycle). This approach can occur in two ways, according to the values of the parameters: (i) a slow decay towards the limit cycle over a time longer than the period, or (ii) a fast decay towards the limit cycle over a time shorter than the period. In the first case, an exactly soluble average equation is valid. We discuss the consequences of our model for the thermal stability of satellites.

Differential patterns of acquired virulence genes distinguish Salmonella strains

Analysis of several Salmonella typhimurium in vivo-induced genes located in regions of atypical base composition has uncovered acquired genetic elements that cumulatively engender pathogenicity. Many of these regions are associated with mobile elements, encode predicted adhesin and invasin-like functions, and are required for full virulence. Some of these regions distinguish broad host range from host-adapted Salmonella serovars and may contribute to inherent differences in host specificity, tissue tropism, and disease manifestation. Maintenance of this archipelago of acquired sequence by selection in specific hosts reveals a fossil record of the evolution of pathogenic species.

Accelerated evolution as a consequence of transitions to mutualism

Differential rates of nucleotide substitutions among taxa are a common observation in molecular phylogenetic studies, yet links between rates of DNA evolution and traits or behaviors of organisms have proved elusive. Likelihood ratio testing is used here for the first time to evaluate specific hypotheses that account for the induction of shifts in rates of DNA evolution. A molecular phylogenetic investigation of mutualist (lichen-forming fungi and fungi associated with liverworts) and nonmutualist fungi revealed four independent transitions to mutualism. We demonstrate a highly significant association between mutualism and increased rates of nucleotide substitutions in nuclear ribosomal DNA, and we demonstrate that a transition to mutualism preceded the rate acceleration of nuclear ribosomal DNA in these lineages. Our results suggest that the increased rate of evolution after the adoption of a mutualist lifestyle is generalized across the genome of these mutualist fungi.

Lamprey Dlx genes and early vertebrate evolution

Gnathostome vertebrates have multiple members of the Dlx family of transcription factors that are expressed during the development of several tissues considered to be vertebrate synapomorphies, including the forebrain, cranial neural crest, placodes, and pharyngeal arches. The Dlx gene family thus presents an ideal system in which to examine the relationship between gene duplication and morphological innovation during vertebrate evolution. Toward this end, we have cloned Dlx genes from the lamprey Petromyzon marinus, an agnathan vertebrate that occupies a critical phylogenetic position between cephalochordates and gnathostomes. We have identified four Dlx genes in P. marinus, whose orthology with gnathostome Dlx genes provides a model for how this gene family evolved in the vertebrate lineage. Differential expression of these lamprey Dlx genes in the forebrain, cranial neural crest, pharyngeal arches, and sensory placodes of lamprey embryos provides insight into the developmental evolution of these structures as well as a model of regulatory evolution after Dlx gene duplication events.

Declines of biomes and biotas and the future of evolution

Although panel discussants disagreed whether the biodiversity crisis constitutes a mass extinction event, all agreed that current extinction rates are 50–500 times background and are increasing and that the consequences for the future evolution of life are serious. In response to the on-going rapid decline of biomes and homogenization of biotas, the panelists predicted changes in species geographic ranges, genetic risks of extinction, genetic assimilation, natural selection, mutation rates, the shortening of food chains, the increase in nutrient-enriched niches permitting the ascendancy of microbes, and the differential survival of ecological generalists. Rates of evolutionary processes will change in different groups, and speciation in the larger vertebrates is essentially over. Action taken over the next few decades will determine how impoverished the biosphere will be in 1,000 years when many species will suffer reduced evolvability and require interventionist genetic and ecological management. Whether the biota will continue to provide the dependable ecological services humans take for granted is less clear. The discussants offered recommendations, including two of paramount importance (concerning human populations and education), seven identifying specific scientific activities to better equip us for stewardship of the processes of evolution, and one suggesting that such stewardship is now our responsibility. The ultimate test of evolutionary biology as a science is not whether it solves the riddles of the past but rather whether it enables us to manage the future of the biosphere. Our inability to make clearer predictions about the future of evolution has serious consequences for both biodiversity and humanity.

Differential Control of Xanthophylls and Light-Induced Stress Proteins, as Opposed to Light-Harvesting Chlorophyll a/b Proteins, during Photosynthetic Acclimation of Barley Leaves to Light Irradiance

Barley (Hordeum vulgare L.) plants were grown at different photon flux densities ranging from 100 to 1800 μmol m−2 s−1 in air and/or in atmospheres with reduced levels of O2 and CO2. Low O2 and CO2 partial pressures allowed plants to grow under high photosystem II (PSII) excitation pressure, estimated in vivo by chlorophyll fluorescence measurements, at moderate photon flux densities. The xanthophyll-cycle pigments, the early light-inducible proteins, and their mRNA accumulated with increasing PSII excitation pressure irrespective of the way high excitation pressure was obtained (high-light irradiance or decreased CO2 and O2 availability). These findings indicate that the reduction state of electron transport chain components could be involved in light sensing for the regulation of nuclear-encoded chloroplast gene expression. In contrast, no correlation was found between the reduction state of PSII and various indicators of the PSII light-harvesting system, such as the chlorophyll a-to-b ratio, the abundance of the major pigment-protein complex of PSII (LHCII), the mRNA level of LHCII, the light-saturation curve of O2 evolution, and the induced chlorophyll-fluorescence rise. We conclude that the chlorophyll antenna size of PSII is not governed by the redox state of PSII in higher plants and, consequently, regulation of early light-inducible protein synthesis is different from that of LHCII.

Differential expression of the ToxR regulon in classical and E1 Tor biotypes of Vibrio cholerae is due to biotype-specific control over toxT expression.

The two major disease-causing biotypes of Vibrio cholerae, classical and El Tor, exhibit differences in their epidemic nature. Their behavior in the laboratory also differs in that El Tor strains produce two major virulence factors, cholera toxin (CT) and the toxin coregulated pilus (TCP), only under very restricted growth conditions, whereas classical strains do so in standard laboratory medium. Expression of toxin and TCP is controlled by two activator proteins, ToxR and ToxT, that operate in cascade fashion with ToxR controlling the synthesis of ToxT. Both biotypes express equivalent levels of ToxR, but only classical strains appear to express ToxT when grown in standard medium. In this report we show that restrictive expression of CT and TCP can be overcome in El Tor strains by expressing ToxT independently of ToxR. An El Tor strain lacking functional ToxT does not express CT or TCP, ruling out existence of a cryptic pathway for virulence regulation in this biotype. These results may have implications for understanding the evolution of El Tor strains toward reduced virulence with respect to classical strains.

Myc and Max: molecular evolution of a family of proto-oncogene products and their dimerization partner.

The myc gene family encodes a group of transcription factors that regulate cell proliferation and differentiation. These genes are widely studied because of their importance as proto-oncogenes. Phylogenetic analyses are described here for 45 Myc protein sequences representing c-, N-, L-, S-, and B-myc genes. A gene duplication early in vertebrate evolution produced the c-myc lineage and another lineage that later gave rise to the N- and L-myc lineages by another gene duplication. Evolutionary divergence in the myc gene family corresponds closely to the known branching order of the major vertebrate groups. The patterns of sequence evolution are described for five separate highly conserved regions, and these analyses show that differential rates of sequence divergence (= mosaic evolution) have occurred among conserved motifs. Further, the closely related dimerization partner protein Max exhibits significantly less sequence variability than Myc. It is suggested that the reduced variability in max stems from natural selection acting to preserve dimerization capability with products of myc and related genes.

Identification of the human prostatic carcinoma oncogene PTI-1 by rapid expression cloning and differential RNA display.

Elucidating the relevant genomic changes mediating development and evolution of prostate cancer is paramount for effective diagnosis and therapy. A putative dominant-acting nude mouse prostatic carcinoma tumor-inducing gene, PTI-1, has been cloned that is expressed in patient-derived human prostatic carcinomas but not in benign prostatic hypertrophy or normal prostate tissue. PTI-1 was detected by cotransfecting human prostate carcinoma DNA into CREF-Trans 6 cells, inducing tumors in nude mice, and isolating genes displaying increased expression in tumor-derived cells by using differential RNA display (DD). Screening a human prostatic carcinoma (LNCaP) cDNA library with a 214-bp DNA fragment found by DD permitted the cloning of a full-length 2.0-kb PTI-1 cDNA. Sequence analysis indicates that PTI-1 is a gene containing a 630-bp 5' sequence and a 3' sequence homologous to a truncated and mutated form of human elongation factor 1 alpha. In vitro translation demonstrates that the PTI-1 cDNA encodes a predominant approximately 46-kDa protein. Probing Northern blots with a DNA fragment corresponding to the 5' region of PTI-1 identifies multiple PTI-1 transcripts in RNAs from human carcinoma cell lines derived from the prostate, lung, breast, and colon. In contrast, PTI-1 RNA is not detected in human melanoma, neuroblastoma, osteosarcoma, normal cerebellum, or glioblastoma multiforme cell lines. By using a pair of primers recognizing a 280-bp region within the 630-bp 5' PTI-1 sequence, reverse transcription-PCR detects PTI-1 expression in patient-derived prostate carcinomas but not in normal prostate or benign hypertrophic prostate tissue. In contrast, reverse transcription-PCR detects prostate-specific antigen expression in all of the prostate tissues. These results indicate that PTI-1 may be a member of a class of oncogenes that could affect protein translation and contribute to carcinoma development in human prostate and other tissues. The approaches used, rapid expression cloning with the CREF-Trans 6 system and the DD strategy, should prove widely applicable for identifying and cloning additional human oncogenes.

Three-dimensional inverse modelling of magnetic anomaly sources based on a genetic algorithm.

We present a modelling method to estimate the 3-D geometry and location of homogeneously magnetized sources from magnetic anomaly data. As input information, the procedure needs the parameters defining the magnetization vector (intensity, inclination and declination) and the Earth's magnetic field direction. When these two vectors are expected to be different in direction, we propose to estimate the magnetization direction from the magnetic map. Then, using this information, we apply an inversion approach based on a genetic algorithm which finds the geometry of the sources by seeking the optimum solution from an initial population of models in successive iterations through an evolutionary process. The evolution consists of three genetic operators (selection, crossover and mutation), which act on each generation, and a smoothing operator, which looks for the best fit to the observed data and a solution consisting of plausible compact sources. The method allows the use of non-gridded, non-planar and inaccurate anomaly data and non-regular subsurface partitions. In addition, neither constraints for the depth to the top of the sources nor an initial model are necessary, although previous models can be incorporated into the process. We show the results of a test using two complex synthetic anomalies to demonstrate the efficiency of our inversion method. The application to real data is illustrated with aeromagnetic data of the volcanic island of Gran Canaria (Canary Islands).

Evolution of the energy consumed by street lighting in Spain estimated with DMSP-OLS data

We present the results of the analysis of satellite imagery to study light pollution in Spain. Both calibrated and non-calibrated DMSP-OLS images were used. We describe the method to scale the non-calibrated DMSP-OLS images which allows us to use differential photometry techniques in order to study the evolution of the light pollution. Population data and DMSP-OLS satellite calibrated images for the year 2006 were compared to test the reliability of official statistics in public lighting consumption. We found a relationship between the population and the energy consumption which is valid for several regions. Finally the true evolution of the electricity consumption for street lighting in Spain from 1992 to 2010 was derived; it has been doubled in the last 18 years in most of the provinces. (C) 2013 Elsevier Ltd. All rights reserved,

Forensic analysis of buildings affected by mining subsidence based on Differential Interferometry (Part III)

This work presents a forensic analysis of buildings affected by mining subsidence, which is based on deformation data obtained by Differential Interferometry (DInSAR). The proposed test site is La Union village (Murcia, SE Spain) where subsidence was triggered in an industrial area due to the collapse of abandoned underground mining labours occurred in 1998. In the first part of this work the study area was introduced, describing the spatial and temporal evolution of ground subsidence, through the elaboration of a cracks map on the buildings located within the affected area. In the second part, the evolution of the most significant cracks found in the most damaged buildings was monitored using biaxial extensometric units and inclinometers. This article describes the work performed in the third part, where DInSAR processing of satellite radar data, available between 1998 and 2008, has permitted to determine the spatial and temporal evolution of the deformation of all the buildings of the study area in a period when no continuous in situ instrumental data is available. Additionally, the comparison of these results with the forensic data gathered in the 2005–2008 period, reveal that there is a coincidence between damaged buildings, buildings where extensometers register significant movements of cracks, and buildings deformation estimated from radar data. As a result, it has been demonstrated that the integration of DInSAR data into forensic analysis methodologies contributes to improve significantly the assessment of the damages of buildings affected by mining subsidence.

A memetic algorithm for the delineation of local labour markets

Given a territory composed of basic geographical units, the delineation of local labour market areas (LLMAs) can be seen as a problem in which those units are grouped subject to multiple constraints. In previous research, standard genetic algorithms were not able to find valid solutions, and a specific evolutionary algorithm was developed. The inclusion of multiple ad hoc operators allowed the algorithm to find better solutions than those of a widely-used greedy method. However, the percentage of invalid solutions was still very high. In this paper we improve that evolutionary algorithm through the inclusion of (i) a reparation process, that allows every invalid individual to fulfil the constraints and contribute to the evolution, and (ii) a hillclimbing optimisation procedure for each generated individual by means of an appropriate reassignment of some of its constituent units. We compare the results of both techniques against the previous results and a greedy method.