Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

In vivo evaluation of the translations of the gleno-humeral joint using Magnetic resonance imaging

Gleno-humeral joint (GHJ) is the most mobile joint of the human body. This is related to theincongr uence between the large humeral head articulating with the much smaller glenoid (ratio 3:1). The GHJ laxity is the ability of the humeral head to be passively translated on the glenoid fossa and, when physiological, it guarantees the normal range of motion of the joint. Three-dimensional GHJ linear displacements have been measured, both in vivo and in vitro by means of different instrumental techniques. In vivo gleno-humeral displacements have been assessed by means of stereophotogrammetry, electromagnetic tracking sensors, and bio-imaging techniques. Both stereophotogrammetric systems and electromagnetic tracking devices, due to the deformation of the soft tissues surrounding the bones, are not capable to accurately assess small displacements, such as gleno-humeral joint translations. The bio-imaging techniques can ensure for an accurate joint kinematic (linear and angular displacement) description, but, due to the radiation exposure, most of these techniques, such as computer tomography or fluoroscopy, are invasive for patients. Among the bioimaging techniques, an alternative which could provide an acceptable level of accuracy and that is innocuous for patients is represented by magnetic resonance imaging (MRI). Unfortunately, only few studies have been conducted for three-dimensional analysis and very limited data is available in situations where preset loads are being applied. The general aim of this doctoral thesis is to develop a non-invasive methodology based on open-MRI for in-vivo evaluation of the gleno-humeral translation components in healthy subjects under the application of external loads.

Application of Submerged Grouted Anchors in Sheet Pile Quay Walls

Sheet pile walls are one of the oldest earth retention systems utilized in civil engineering projects. They are used for various purposes; such as excavation support system, cofferdams, cut-off walls under dams, slope stabilization, waterfront structures, and flood walls. Sheet pile walls are one of the most common types of quay walls used in port construction. The worldwide increases in utilization of large ships for transportation have created an urgent need of deepening the seabed within port areas and consequently the rehabilitation of its wharfs. Several methods can be used to increase the load-carrying capacity of sheet-piling walls. The use of additional anchored tie rods grouted into the backfill soil and arranged along the exposed wall height is one of the most practical and appropriate solutions adopted for stabilization and rehabilitation of the existing quay wall. The Ravenna Port Authority initiated a project to deepen the harbor bottom at selected wharves. An extensive parametric study through the finite element program, PLAXIS 2D, version 2012 was carried out to investigate the enhancement of using submerged grouted anchors technique on the load response of sheet-piling quay wall. The influence of grout-ties area, length of grouted body, anchor inclination and anchor location were considered and evaluated due to the effect of different system parameters. Also a comparative study was conducted by Plaxis 2D and 3D program to investigate the behavior of these sheet pile quay walls in terms of horizontal displacements induced along the sheet pile wall and ground surface settlements as well as the anchor force and calculated factor of safety. Finally, a comprehensive study was carried out by using different constitutive models to simulate the mechanical behavior of the soil to investigate the effect of these two models (Mohr-Coulomb and Hardening Soil) on the behavior of these sheet pile quay walls.

The reconstruction of skeletal movement: the soft tissue artefact issue

In 3D human movement analysis performed using stereophotogrammetric systems and skin markers, bone pose can only be estimated in an indirect fashion. During a movement, soft tissue deformations make the markers move with respect to the underlying bone generating soft tissue artefact (STA). STA has devastating effects on bone pose estimation and its compensation remains an open question. The aim of this PhD thesis was to contribute to the solution of this crucial issue. Modelling STA using measurable trial-specific variables is a fundamental prerequisite for its removal from marker trajectories. Two STA model architectures are proposed. Initially, a thigh marker-level artefact model is presented. STA was modelled as a linear combination of joint angles involved in the movement. This model was calibrated using ex-vivo and in-vivo STA invasive measures. The considerable number of model parameters led to defining STA approximations. Three definitions were proposed to represent STA as a series of modes: individual marker displacements, marker-cluster geometrical transformations (MCGT), and skin envelope shape variations. Modes were selected using two criteria: one based on modal energy and another on the selection of modes chosen a priori. The MCGT allows to select either rigid or non-rigid STA components. It was also empirically demonstrated that only the rigid component affects joint kinematics, regardless of the non-rigid amplitude. Therefore, a model of thigh and shank STA rigid component at cluster-level was then defined. An acceptable trade-off between STA compensation effectiveness and number of parameters can be obtained, improving joint kinematics accuracy. The obtained results lead to two main potential applications: the proposed models can generate realistic STAs for simulation purposes to compare different skeletal kinematics estimators; and, more importantly, focusing only on the STA rigid component, the model attains a satisfactory STA reconstruction with less parameters, facilitating its incorporation in an pose estimator.

Spatial and temporal characterisation of ground deformation recorded by geodetic techniques

A critical point in the analysis of ground displacements time series is the development of data driven methods that allow the different sources that generate the observed displacements to be discerned and characterised. A widely used multivariate statistical technique is the Principal Component Analysis (PCA), which allows reducing the dimensionality of the data space maintaining most of the variance of the dataset explained. Anyway, PCA does not perform well in finding the solution to the so-called Blind Source Separation (BSS) problem, i.e. in recovering and separating the original sources that generated the observed data. This is mainly due to the assumptions on which PCA relies: it looks for a new Euclidean space where the projected data are uncorrelated. The Independent Component Analysis (ICA) is a popular technique adopted to approach this problem. However, the independence condition is not easy to impose, and it is often necessary to introduce some approximations. To work around this problem, I use a variational bayesian ICA (vbICA) method, which models the probability density function (pdf) of each source signal using a mix of Gaussian distributions. This technique allows for more flexibility in the description of the pdf of the sources, giving a more reliable estimate of them. Here I present the application of the vbICA technique to GPS position time series. First, I use vbICA on synthetic data that simulate a seismic cycle (interseismic + coseismic + postseismic + seasonal + noise) and a volcanic source, and I study the ability of the algorithm to recover the original (known) sources of deformation. Secondly, I apply vbICA to different tectonically active scenarios, such as the 2009 L'Aquila (central Italy) earthquake, the 2012 Emilia (northern Italy) seismic sequence, and the 2006 Guerrero (Mexico) Slow Slip Event (SSE).

Lateral spreading and associated slope processes in fractured rock slabs

Landslides of the lateral spreading type, involving brittle geological units overlying ductile terrains, are a common occurrence in the sandstone and limestone plateaux of the northern Apennines of Italy. These instability phenomena can become particularly risky, when historical towns and cultural heritage sites built on the top of them are endangered. Neverthless, the mechanisms controlling the developing of related instabilities, i.e. toppling and rock falls, at the edges of rock plateaux are not fully understood yet. In addition, the groundwater flow path developing at the contact between the more permeable units, i.e. the jointed rock slab, and the relatively impermeable clay-rich units have not been already studied in details, even if they may play a role in this kind of instability processes, acting as eventual predisposing and/or triggering factors. Field survey, Terrestrial Laser Scanner and Close Range Photogrammetry techniques, laboratory tests on the involved materials, hydrogeological monitoring and modelling, displacements evaluation and stability analysis through continuum and discontinuum numerical codes have been performed on the San Leo case study, with the aim to bring further insights for the understanding and the assessment of the slope processes taking place in this geological context. The current research permitted to relate the aquifer behaviour of the rocky slab to slope instability processes. The aquifer hosted in the fractured slab leads to the development of perennial and ephemeral springs at the contact between the two units. The related piping erosion phenomena, together with slope processes in the clay-shales led to the progressive undermining of the slab. The cliff becomes progressively unstable due to undermining and undergoes large-scale landslides due to fall or topple.

Prediction and measurement of forced response of a composite blade undergoing nonlinear vibration

The main objective of this project is to experimentally demonstrate geometrical nonlinear phenomena due to large displacements during resonant vibration of composite materials and to explain the problem associated with fatigue prediction at resonant conditions. Three different composite blades to be tested were designed and manufactured, being their difference in the composite layup (i.e. unidirectional, cross-ply, and angle-ply layups). Manual envelope bagging technique is explained as applied to the actual manufacturing of the components; problems encountered and their solutions are detailed. Forced response tests of the first flexural, first torsional, and second flexural modes were performed by means of a uniquely contactless excitation system which induced vibration by using a pulsed airflow. Vibration intensity was acquired by means of Polytec LDV system. The first flexural mode is found to be completely linear irrespective of the vibration amplitude. The first torsional mode exhibits a general nonlinear softening behaviour which is interestingly coupled with a hardening behaviour for the unidirectional layup. The second flexural mode has a hardening nonlinear behaviour for either the unidirectional and angle-ply blade, whereas it is slightly softening for the cross-ply layup. By using the same equipment as that used for forced response analyses, free decay tests were performed at different airflow intensities. Discrete Fourier Trasform over the entire decay and Sliding DFT were computed so as to visualise the presence of nonlinear superharmonics in the decay signal and when they were damped out from the vibration over the decay time. Linear modes exhibit an exponential decay, while nonlinearities are associated with a dry-friction damping phenomenon which tends to increase with increasing amplitude. Damping ratio is derived from logarithmic decrement for the exponential branch of the decay.

Hydrodynamic correlation functions of a driven granular fluid in steady state

We study a homogeneously driven granular fluid of hard spheres at intermediate volume fractions and focus on time-delayed correlation functions in the stationary state. Inelastic collisions are modeled by incomplete normal restitution, allowing for efficient simulations with an event-driven algorithm. The incoherent scattering function Fincoh(q,t ) is seen to follow time-density superposition with a relaxation time that increases significantly as the volume fraction increases. The statistics of particle displacements is approximately Gaussian. For the coherent scattering function S(q,ω), we compare our results to the predictions of generalized fluctuating hydrodynamics, which takes into account that temperature fluctuations decay either diffusively or with a finite relaxation rate, depending on wave number and inelasticity. For sufficiently small wave number q we observe sound waves in the coherent scattering function S(q,ω) and the longitudinal current correlation function Cl(q,ω). We determine the speed of sound and the transport coefficients and compare them to the results of kinetic theory.

The Estimation of the Ecological Condition of the City Environment by Methods of Phytoindication (as Applied to the Vasileostrovsky District of St. Petersburg)

The project aimed to use results of contamination of city vegetation with heavy metals and sulphur compounds as the basis for analysing the integral response of trees and shrubs to contamination, through a complex method of phytoindication. The results were used to draw up recommendations on pollution reduction in the city and to develop the method of phytoindication as a means of monitoring environmental pollution in St. Petersburg and other large cities. Field investigations were carried out in August 1996, and 66 descriptions of green areas were made in order to estimate the functional state of plants in the Vasileostrovsky district. Investigations of the spectrum reflecting properties of plants showed considerable variation of albedo meanings of leaves under the influence of various internal and external factors. The results indicated that lime trees most closely reflect the condition of the environment. Practically all the green areas studied were in poor condition, the only exceptions being areas of ash trees, which are more resistant to environmental pollution, and one lime-tree alley in a comparatively unpolluted street. The study identified those types of trees which are more or less resistant to complex environmental pollution and Ms. Terekhina recommends that the species in the present green areas be changed to include a higher number of the more resistant species. The turbidimetric analysis of tree barks for sulphates gave an indication of the level and spatial distribution of each pollutant, and the results also confirmed other findings that electric conductivity is a significant feature in determining the extent of sulphate pollution. In testing for various metals, the lime tree showed the highest contents for all elements except magnesium, copper, zinc, cadmium and strontium, again confirming the species' vulnerability to pollution. Medium rates of concentration in the city and environs showed that city plants concentrate 3 times as many different elements and 10 times more chromium, copper and lead than do those in the suburbs. The second stage of the study was based on the concept of phytoindication, which presupposes that changes in the relation of chemical elements in regional biological circulation under the influence of technogenesis provide a criterion for predicting displacements in people's health. There are certain basic factors in this concept. The first is that all living beings are related ecologically as well as by their evolutionary origin, and that the lower an organism is on the evolutionary scale, the less adaptational reserve it has. The second is that smaller concentrations of chemical elements are needed for toxicological influence on plants than on people and so the former's reactions to geochemical factors are easier to characterise. Visual indicational features of urban plants are well defined and can form the basis of a complex "environment - public health" analysis. Specific plant reactions reflecting atmospheric pollution and other components of urbogeosystems make it possible to determine indication criteria for predicting possible disturbances in the general state of health of the population. Thirdly the results of phytoindication investigations must be taken together with information about public health in the area. It only proved possibly to analyse general indexes of public health based on statistical data from the late 1980s and early 1990s as the data of later years were greatly influenced by social factors. These data show that the rates of illness in St. Petersburg (especially for children) are higher than in Russia as a whole, for most classes of diseases, indicating that the population there is more sensitive to the ecological state of the urban environment. The Vasileostrovsky district had the second highest sick rate for adullts, while the rate of infant mortality in the first year of life was highest there. Ms. Terekhina recommends further studies to more precisely assess the effectiveness of the methods she tested, but has drawn up a proposed map of environmental hazard for the population, taking into account prevailing wind directions.

Impact of Holocene climate changes on alpine and treeline vegetation at Sanetsch Pass, Bernese Alps, Switzerland

In order to infer reactions of treeline and alpine vegetation to climatic change, past vegetation changes are reconstructed on the basis of pollen, macrofossil and charcoal analysis. The sampled sediment cores originate from the small pond Emines, located at the Sanetsch Pass (connecting the Valais and Bern, Switzerland) at an altitude of 2288 m a.s.l. Today's treeline is at ca. 2200 m a.s.l. in the area, though due to special pass (saddle) conditions it is locally depressed to ca. 2060 m a.s.l. Our results reveal that the area around Emines was covered by treeless alpine vegetation during most of the past 12,000 years. Single individuals of Betula, Larix decidua and possibly Pinus cembra occurred during the Holocene. Major centennial to millennial-scale responses of treeline vegetation to climatic changes are evident. However, alpine vegetation composition remained rather stable between 11,500 and 6000 cal. BP, showing that Holocene climatic changes of +/− 1 °C hardly influenced the local vegetation at Emines. The rapid warming of 3–4 °C at the Late Glacial/Holocene transition (11,600 cal. BP) caused significant altitudinal displacements of alpine species that were additionally affected by the rapid upward movement of trees and shrubs. Since the beginning of the Neolithic, vegetation changes at Sanetsch Pass resulted from a combination of climate change and human impact. Anthropogenic fire increase and land-use change combined with a natural change from subcontinental to more oceanic climate during the second half of the Holocene led to the disappearance of P. cembra in the study area, but favoured the occurrence of Picea abies and Alnus viridis. The mid- to late-Holocene decline of Abies alba was primarily a consequence of human impact, since this mesic species should have benefitted from a shift to more oceanic conditions. Future alpine vegetation changes will be a function of the amplitude and rapidity of global warming as well as human land use. Our results imply that alpine vegetation at our treeline pass site was never replaced by forests since the last ice-age. This may change in the future if anticipated climate change will induce upslope migration of trees. The results of this study emphasise the necessity of climate change mitigation in order to prevent biodiversity losses as a consequence of unprecedented community and species displacement in response to climatic change.

Daily organ tracking in intensity-modulated radiotherapy of prostate cancer using an electronic portal imaging device with a dose saving acquisition mode

BACKGROUND AND PURPOSE: Daily use of conventional electronic portal imaging devices (EPID) for organ tracking is limited due to the relatively high dose required for high quality image acquisition. We studied the use of a novel dose saving acquisition mode (RadMode) allowing to take images with one monitor unit per image in prostate cancer patients undergoing intensity-modulated radiotherapy (IMRT) and tracking of implanted fiducial gold markers. PATIENTS AND METHODS: Twenty five patients underwent implantation of three fiducial gold markers prior to the planning CT. Before each treatment of a course of 37 fractions, orthogonal localization images from the antero-posterior and from the lateral direction were acquired. Portal images of both the setup procedure and the five IMRT treatment beams were analyzed. RESULTS: On average, four localization images were needed for a correct patient setup, resulting in four monitor units extra dose per fraction. The mean extra dose delivered to the patient was thereby increased by 1.2%. The procedure was precise enough to reduce the mean displacements prior to treatment to < o =0.3 mm. CONCLUSIONS: The use of a new dose saving acquisition mode enables to perform daily EPID-based prostate tracking with a cumulative extra dose of below 1 Gy. This concept is efficiently used in IMRT-treated patients, where separation of setup beams from treatment beams is mandatory.

Differential GPS as a monitoring tool on Volcano Santa Ana (Illamatepec) and the Coatepeque Caldera, El Salvador

We used differential GPS measurements from a 13 station GPS network spanning the Santa Ana Volcano and Coatepeque Caldera to characterize the inter-eruptive activity and tectonic movements near these two active and potentially hazardous features. Caldera-forming events occurred from 70-40 ka and at Santa Ana/Izalco volcanoes eruptive activity occurred as recently as 2005. Twelve differential stations were surveyed for 1 to 2 hours on a monthly basis from February through September 2009 and tied to a centrally located continuous GPS station, which serves as the reference site for this volcanic network. Repeatabilities of the averages from 20-minute sessions taken over 20 hours or longer range from 2-11 mm in the horizontal (north and east) components of the inter-station baselines, suggesting a lower detection limit for the horizontal components of any short-term tectonic or volcanic deformation. Repeatabilities of the vertical baseline component range from 12-34 mm. Analysis of the precipitable water vapor in the troposphere suggests that tropospheric decorrelation as a function of baseline lengths and variable site elevations are the most likely sources of vertical error. Differential motions of the 12 sites relative to the continuous reference site reveal inflation from February through July at several sites surrounding the caldera with vertical displacements that range from 61 mm to 139 mm followed by a lower magnitude deflation event on 1.8-7.4 km-long baselines. Uplift rates for the inflationary period reach 300 mm/yr with 1σ uncertainties of +/- 26 – 119 mm. Only one other station outside the caldera exhibits a similar deformation trend, suggesting a localized source. The results suggest that the use of differential GPS measurements from short duration occupations over short baselines can be a useful monitoring tool at sub-tropical volcanoes and calderas.

Energy harvesting from body motion using rotational micro-generation

Autonomous system applications are typically limited by the power supply operational lifetime when battery replacement is difficult or costly. A trade-off between battery size and battery life is usually calculated to determine the device capability and lifespan. As a result, energy harvesting research has gained importance as society searches for alternative energy sources for power generation. For instance, energy harvesting has been a proven alternative for powering solar-based calculators and self-winding wristwatches. Thus, the use of energy harvesting technology can make it possible to assist or replace batteries for portable, wearable, or surgically-implantable autonomous systems. Applications such as cardiac pacemakers or electrical stimulation applications can benefit from this approach since the number of surgeries for battery replacement can be reduced or eliminated. Research on energy scavenging from body motion has been investigated to evaluate the feasibility of powering wearable or implantable systems. Energy from walking has been previously extracted using generators placed on shoes, backpacks, and knee braces while producing power levels ranging from milliwatts to watts. The research presented in this paper examines the available power from walking and running at several body locations. The ankle, knee, hip, chest, wrist, elbow, upper arm, side of the head, and back of the head were the chosen target localizations. Joints were preferred since they experience the most drastic acceleration changes. For this, a motor-driven treadmill test was performed on 11 healthy individuals at several walking (1-4 mph) and running (2-5 mph) speeds. The treadmill test provided the acceleration magnitudes from the listed body locations. Power can be estimated from the treadmill evaluation since it is proportional to the acceleration and frequency of occurrence. Available power output from walking was determined to be greater than 1mW/cm³ for most body locations while being over 10mW/cm³ at the foot and ankle locations. Available power from running was found to be almost 10 times higher than that from walking. Most energy harvester topologies use linear generator approaches that are well suited to fixed-frequency vibrations with sub-millimeter amplitude oscillations. In contrast, body motion is characterized with a wide frequency spectrum and larger amplitudes. A generator prototype based on self-winding wristwatches is deemed to be appropriate for harvesting body motion since it is not limited to operate at fixed-frequencies or restricted displacements. Electromagnetic generation is typically favored because of its slightly higher power output per unit volume. Then, a nonharmonic oscillating rotational energy scavenger prototype is proposed to harness body motion. The electromagnetic generator follows the approach from small wind turbine designs that overcome the lack of a gearbox by using a larger number of coil and magnets arrangements. The device presented here is composed of a rotor with multiple-pole permanent magnets having an eccentric weight and a stator composed of stacked planar coils. The rotor oscillations induce a voltage on the planar coil due to the eccentric mass unbalance produced by body motion. A meso-scale prototype device was then built and evaluated for energy generation. The meso-scale casing and rotor were constructed on PMMA with the help of a CNC mill machine. Commercially available discrete magnets were encased in a 25mm rotor. Commercial copper-coated polyimide film was employed to manufacture the planar coils using MEMS fabrication processes. Jewel bearings were used to finalize the arrangement. The prototypes were also tested at the listed body locations. A meso-scale generator with a 2-layer coil was capable to extract up to 234 µW of power at the ankle while walking at 3mph with a 2cm³ prototype for a power density of 117 µW/cm³. This dissertation presents the analysis of available power from walking and running at different speeds and the development of an unobtrusive miniature energy harvesting generator for body motion. Power generation indicates the possibility of powering devices by extracting energy from body motion.

The use of intramedullary helix wire for the treatment of proximal humerus fractures

OBJECTIVE: To present the functional and radiographic outcome 1 and 6 years after application of a new intramedullary fixation device for proximal humerus fractures. DESIGN: Retrospective case series. SETTING: Level II orthopaedic surgery hospital. PATIENTS: Twenty-six consecutive patients (average age 68.9 years) with 2-, 3- and 4-part fractures of the proximal humerus were operated at a single institution. Follow-up was performed after 1 year (26 patients) and 6 years (16 patients). INTERVENTION: All patients were treated with closed reduction and intramedullary helix wires. MAIN OUTCOME MEASUREMENTS: The Constant-Murley score and the University of California Los Angeles (UCLA) score. Clinical complications and radiological posttraumatic arthritis were recorded. RESULTS: The average Constant-Murley score was 70.3 (points) and 70.7 after 1 and 6 years, respectively; the average UCLA score was 27.2 and 31.5 after 1 and 6 years, respectively. Major complications were 4 revisions for 3 secondary fragment displacements and 1 nonunion with partial avascular osteonecrosis in the first postoperative year. Complications were found predominantly in 4-part fractures (3/5, 60%). There were no further complications or progressive posttraumatic arthritis up to 6 years following surgery. CONCLUSION: The helix wire is well suited for displaced or unstable 2- and 3-part proximal humerus fractures. Adequate functional outcome, a low number of implant displacements, a low number of application morbidity, and infrequent implant removals were recorded. The use of this device is not recommended for 4-part fractures.

Structural response and dynamics of fluid-structure-control interaction in wind turbine blades

Reducing the uncertainties related to blade dynamics by the improvement of the quality of numerical simulations of the fluid structure interaction process is a key for a breakthrough in wind-turbine technology. A fundamental step in that direction is the implementation of aeroelastic models capable of capturing the complex features of innovative prototype blades, so they can be tested at realistic full-scale conditions with a reasonable computational cost. We make use of a code based on a combination of two advanced numerical models implemented in a parallel HPC supercomputer platform: First, a model of the structural response of heterogeneous composite blades, based on a variation of the dimensional reduction technique proposed by Hodges and Yu. This technique has the capacity of reducing the geometrical complexity of the blade section into a stiffness matrix for an equivalent beam. The reduced 1-D strain energy is equivalent to the actual 3-D strain energy in an asymptotic sense, allowing accurate modeling of the blade structure as a 1-D finite-element problem. This substantially reduces the computational effort required to model the structural dynamics at each time step. Second, a novel aerodynamic model based on an advanced implementation of the BEM(Blade ElementMomentum) Theory; where all velocities and forces are re-projected through orthogonal matrices into the instantaneous deformed configuration to fully include the effects of large displacements and rotation of the airfoil sections into the computation of aerodynamic forces. This allows the aerodynamic model to take into account the effects of the complex flexo-torsional deformation that can be captured by the more sophisticated structural model mentioned above. In this thesis we have successfully developed a powerful computational tool for the aeroelastic analysis of wind-turbine blades. Due to the particular features mentioned above in terms of a full representation of the combined modes of deformation of the blade as a complex structural part and their effects on the aerodynamic loads, it constitutes a substantial advancement ahead the state-of-the-art aeroelastic models currently available, like the FAST-Aerodyn suite. In this thesis, we also include the results of several experiments on the NREL-5MW blade, which is widely accepted today as a benchmark blade, together with some modifications intended to explore the capacities of the new code in terms of capturing features on blade-dynamic behavior, which are normally overlooked by the existing aeroelastic models.