The impact of geometric non-linearities on the fatigue analysis of trailing edge bond lines in wind turbine rotor blades

The accurate prediction of stress histories for the fatigue analysis is of utmost importance for the design process of wind turbine rotor blades. As detailed, transient, and geometrically non-linear three-dimensional finite element analyses are computationally weigh too expensive, it is commonly regarded sufficient to calculate the stresses with a geometrically linear analysis and superimpose different stress states in order to obtain the complete stress histories. In order to quantify the error from geometrically linear simulations for the calculation of stress histories and to verify the practical applicability of the superposition principal in fatigue analyses, this paper studies the influence of geometric non-linearity in the example of a trailing edge bond line, as this subcomponent suffers from high strains in span-wise direction. The blade under consideration is that of the IWES IWT-7.5-164 reference wind turbine. From turbine simulations the highest edgewise loading scenario from the fatigue load cases is used as the reference. A 3D finite element model of the blade is created and the bond line fatigue assessment is performed according to the GL certification guidelines in its 2010 edition, and in comparison to the latest DNV GL standard from end of 2015. The results show a significant difference between the geometrically linear and non-linear stress analyses when the bending moments are approximated via a corresponding external loading, especially in case of the 2010 GL certification guidelines. This finding emphasizes the demand to reconsider the application of the superposition principal in fatigue analyses of modern flexible rotor blades, where geometrical nonlinearities become significant. In addition, a new load application methodology is introduced that reduces the geometrically non-linear behaviour of the blade in the finite element analysis.

Caracterização analítica e geométrica da metodologia geral de determinação de distribuições de Weibull para o regime eólico e suas aplicações

O regime eólico de uma região pode ser descrito por distribuição de frequências que fornecem informações e características extremamente necessárias para uma possível implantação de sistemas eólicos de captação de energia na região e consequentes aplicações no meio rural em regiões afastadas. Estas características, tais como a velocidade média anual, a variância das velocidades registradas e a densidade da potência eólica média horária, podem ser obtidas pela frequência de ocorrências de determinada velocidade, que por sua vez deve ser estudada através de expressões analíticas. A função analítica mais adequada para distribuições eólicas é a função de densidade de Weibull, que pode ser determinada por métodos numéricos e regressões lineares. O objetivo deste trabalho é caracterizar analítica e geometricamente todos os procedimentos metodológicos necessários para a realização de uma caracterização completa do regime eólico de uma região e suas aplicações na região de Botucatu - SP, visando a determinar o potencial energético para implementação de turbinas eólicas. Assim, foi possível estabelecer teoremas relacionados com a forma de caracterização do regime eólico, estabelecendo a metodologia concisa analiticamente para a definição dos parâmetros eólicos de qualquer região a ser estudada. Para o desenvolvimento desta pesquisa, utilizou-se um anemômetro da CAMPBELL.

I. Quantal effects in biochemical cooperativity and a proposed mechanism for the differentiation of calcium signaling in synaptic plasticity. II. Evolutionary algorithms for the optimization of methods in computational chemistry

In Part 1 of this thesis, we propose that biochemical cooperativity is a fundamentally non-ideal process. We show quantal effects underlying biochemical cooperativity and highlight apparent ergodic breaking at small volumes. The apparent ergodic breaking manifests itself in a divergence of deterministic and stochastic models. We further predict that this divergence of deterministic and stochastic results is a failure of the deterministic methods rather than an issue of stochastic simulations.

Ergodic breaking at small volumes may allow these molecular complexes to function as switches to a greater degree than has previously been shown. We propose that this ergodic breaking is a phenomenon that the synapse might exploit to differentiate Ca$^{2+}$ signaling that would lead to either the strengthening or weakening of a synapse. Techniques such as lattice-based statistics and rule-based modeling are tools that allow us to directly confront this non-ideality. A natural next step to understanding the chemical physics that underlies these processes is to consider \textit{in silico} specifically atomistic simulation methods that might augment our modeling efforts.

In the second part of this thesis, we use evolutionary algorithms to optimize \textit{in silico} methods that might be used to describe biochemical processes at the subcellular and molecular levels. While we have applied evolutionary algorithms to several methods, this thesis will focus on the optimization of charge equilibration methods. Accurate charges are essential to understanding the electrostatic interactions that are involved in ligand binding, as frequently discussed in the first part of this thesis.

Strain gradient plasticity modelling of high-pressure torsion

Gradient plasticity modelling combining a micro-structure-related constitutive description of the local material behaviour with a particular gradient plasticity frame is presented. The constitutive formulation is based on a phase mixture model in which the dislocation cell walls and the cell interiors are considered as separate 'phases', the respective dislocation densities entering as internal variables. Two distinct physical mechanisms, which give rise to gradient plasticity, are considered. The first one is associated with the occurrence of geometrically necessary dislocations leading to first-order strain gradients; the second one is associated with the reaction stresses due to plastic strain incompatibilities between neighbouring grains, which lead to second-order strain gradients. These two separate variants of gradient plasticity were applied to the case of high-pressure torsion: a process known to result in a fairly uniform, ultrafine grained structure of metals. It is shown that the two complementary variants of gradient plasticity can both account for the experimental results, thus resolving a controversial issue of the occurrence of a uniform micro-structure as a result of an inherently non-uniform process. © 2007 Elsevier Ltd. All rights reserved.

A model of grain fragmentation based on lattice curvature

A new model is proposed that aims to capture within a single modelling frame all the main microstructural features of a severe plastic deformation process. These are: evolution of the grain size distribution, misorientation distribution, crystallographic texture and the strain-hardening of the material. The model is based on the lattice curvature that develops in all deformed grains. The basic assumption is that lattice rotation within an individual grain is impeded near the grain boundaries by the constraining effects of the neighbouring grains, which gives rise to lattice curvature. On that basis, a fragmentation scheme is developed which is integrated in the Taylor viscoplastic polycrystal model. Dislocation density evolution is traced for each grain, which includes the contribution of geometrically necessary dislocations associated with lattice curvature. The model is applied to equal-channel angular pressing. The role of texture development is shown to be an important element in the grain fragmentation process. Results of this modelling give fairly precise predictions of grain size and grain misorientation distribution. The crystallographic textures are well reproduced and the strength of the material is also reliably predicted based on the modelling of dislocation density evolution coupled with texture development.

Above ground biomass estimation in Mediterranean multiple use systems with high spatial resolution satellite images

Forest biomass has been having an increasing importance in the world economy and in the evaluation of the forests development and monitoring. It was identified as a global strategic reserve, due to its applications in bioenergy, bioproduct development and issues related to reducing greenhouse gas emissions. The estimation of above ground biomass is frequently done with allometric functions per species with plot inventory data. An adequate sampling design and intensity for an error threshold is required. The estimation per unit area is done using an extrapolation method. This procedure is labour demanding and costly. The mail goal of this study is the development of allometric functions for the estimation of above ground biomass with ground cover as independent variable, for forest areas of holm aok (Quercus rotundifolia), cork oak (Quercus suber) and umbrella pine (Pinus pinea) in multiple use systems. Ground cover per species was derived from crown horizontal projection obtained by processing high resolution satellite images, orthorectified, geometrically and atmospheric corrected, with multi-resolution segmentation method and object oriented classification. Forest inventory data were used to estimate plot above ground biomass with published allometric functions at tree level. The developed functions were fitted for monospecies stands and for multispecies stands of Quercus rotundifolia and Quercus suber, and Quercus suber and Pinus pinea. The stand composition was considered adding dummy variables to distinguish monospecies from multispecies stands. The models showed a good performance. Noteworthy is that the dummy variables, reflecting the differences between species, originated improvements in the models. Significant differences were found for above ground biomass estimation with the functions with and without the dummy variables. An error threshold of 10% corresponds to stand areas of about 40 ha. This method enables the overall area evaluation, not requiring extrapolation procedures, for the three species, which occur frequently in multispecies stands.

Above ground biomass functions with vegetation indices for multiple use systems of two evergreen oaks

Remote sensing is a promising approach for above ground biomass estimation, as forest parameters can be obtained indirectly. The analysis in space and time is quite straight forward due to the flexibility of the method to determine forest crown parameters with remote sensing. It can be used to evaluate and monitoring for example the development of a forest area in time and the impact of disturbances, such as silvicultural practices or deforestation. The vegetation indices, which condense data in a quantitative numeric manner, have been used to estimate several forest parameters, such as the volume, basal area and above ground biomass. The objective of this study was the development of allometric functions to estimate above ground biomass using vegetation indices as independent variables. The vegetation indices used were the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Simple Ratio (SR) and Soil-Adjusted Vegetation Index (SAVI). QuickBird satellite data, with 0.70 m of spatial resolution, was orthorectified, geometrically and atmospheric corrected, and the digital number were converted to top of atmosphere reflectance (ToA). Forest inventory data and published allometric functions at tree level were used to estimate above ground biomass per plot. Linear functions were fitted for the monospecies and multispecies stands of two evergreen oaks (Quercus suber and Quercus rotundifolia) in multiple use systems, montados. The allometric above ground biomass functions were fitted considering the mean and the median of each vegetation index per grid as independent variable. Species composition as a dummy variable was also considered as an independent variable. The linear functions with better performance are those with mean NDVI or mean SR as independent variable. Noteworthy is that the two better functions for monospecies cork oak stands have median NDVI or median SR as independent variable. When species composition dummy variables are included in the function (with stepwise regression) the best model has median NDVI as independent variable. The vegetation indices with the worse model performance were EVI and SAVI.

Improving magnetic coupling for battery charging through 3D magnetic flux

The spatial distribution of the magnetic field and the coupling between the coils in the Wireless Power Transfer (WPT) systems is an important aspect to consider in the system design and efficiency optimization. The presented study in this paper is based on tests performed on a physical model. The transmitting (primary) equipment, is an electrical three-phase system, capable to be connected in star or delta (both electrically and geometrically). The measured results allow to describe graphically the magnetic field distribution in three dimensions. The analytical formulas aim to help to understand and to quantify the physical phenomena but they cannot be considered a universal approach and the measurement results help to understand better the observable facts. In the WPT, the key issues that will influence the efficiency, are the alignment of the coils, the spatial orientation of the magnetic field, the detachment and the tilt between the windings, all they changing the magnetic coupling between the transmitter and the receiver of energy. This research is directed not only to the magnetic field distribution but finally, to optimize the energy transfer efficiency.

Convergence time to equilibrium distributions of autonomous and periodic non-autonomous graphs

We present some estimates of the time of convergence to the equilibrium distribution in autonomous and periodic non-autonomous graphs, with ergodic stochastic adjacency matrices, using the eigenvalues of these matrices. On this way we generalize previous results from several authors, that only considered reversible matrices.

A finite element framework for the interplay between delamination and buckling of rubber-like bi-material systems and stretchable electronics

In this study, a finite element (FE) framework for the analysis of the interplay between buckling and delamination of thin layers bonded to soft substrates is proposed. The current framework incorporates the following modeling features: (i) geometrically nonlinear solid shell elements, (ii) geometrically nonlinear cohesive interface elements, and (iii) hyperelastic material constitutive response for the bodies that compose the system. A fully implicit Newton–Raphson solution strategy is adopted to deal with the complex simultaneous presence of geometrical and material nonlinearities through the derivation of the consistent FE formulation. Applications to a rubber-like bi-material system under finite bending and to patterned stiff islands resting on soft substrate for stretchable solar cells subjected to tensile loading are proposed. The results obtained are in good agreement with benchmark results available in the literature, confirming the accuracy and the capabilities of the proposed numerical method for the analysis of complex three-dimensional fracture mechanics problems under finite deformations.