Clinical versus pre-clinical FE models for vertebral body strength predictions

The finite element analysis is an accepted method to predict vertebral body compressive strength. This study compares measurements obtained from in vitro tests with the ones from two different simulation models: clinical quantitative computer tomography (QCT) based homogenized finite element (hFE) models and pre-clinical high-resolution peripheral QCT-based (HR-pQCT) hFE models. About 37 vertebral body sections were prepared by removing end-plates and posterior elements, scanned with QCT (390/450μm voxel size) as well as HR-pQCT (82μm voxel size), and tested in compression up to failure. Non-linear viscous damage hFE models were created from QCT/HT-pQCT images and compared to experimental results based on stiffness and ultimate load. As expected, the predictability of QCT/HR-pQCT-based hFE models for both apparent stiffness (r2=0.685/0.801r2=0.685/0.801) and strength (r2=0.774/0.924r2=0.774/0.924) increased if a better image resolution was used. An analysis of the damage distribution showed similar damage locations for all cases. In conclusion, HR-pQCT-based hFE models increased the predictability considerably and do not need any tuning of input parameters. In contrast, QCT-based hFE models usually need some tuning but are clinically the only possible choice at the moment.

A numerical study of the influence of microvoids in the transverse mechanical response of unidirectional composites

The effect of porosity on the transverse mechanical properties of unidirectional fiber-reinforced composites is studied by means of computational micromechanics. The composite behavior is simulated by the finite element analysis of a representative volume element of the composite microstructure in which the random distribution of fibers and the voids are explicitly included. Two types of voids – interfiber voids and matrix voids – were included in the microstructure and the actual damage mechanisms in the composite, namely matrix and interface failure, were accounted for. It was found that porosity (in the range 1–5%) led to a large reduction in the transverse strength and the influence of both types of voids in the onset and propagation of damage throughout the microstructure was studied under transverse tension and compression. Finally, the failure locus of the composite lamina under transverse tension/compression and out-of-plane shear was obtained by means of computational micromechanics and compared with the predictions of Puck’s model and with experimental data available in the literature. The results show that the strength of composites is significantly reduced by the presence of voids

Método computacional para la aceleración de análisis paramétricos de modificaciones locales en estructuras complejas sometidas a cargas dinámicas

Dentro del análisis y diseño estructural surgen frecuentemente problemas de ingeniería donde se requiere el análisis dinámico de grandes modelos de elementos finitos que llegan a millones de grados de libertad y emplean volúmenes de datos de gran tamaño. La complejidad y dimensión de los análisis se dispara cuando se requiere realizar análisis paramétricos. Este problema se ha abordado tradicionalmente desde diversas perspectivas: en primer lugar, aumentando la capacidad tanto de cálculo como de memoria de los sistemas informáticos empleados en los análisis. En segundo lugar, se pueden simplificar los análisis paramétricos reduciendo su número o detalle y por último se puede recurrir a métodos complementarios a los elementos .nitos para la reducción de sus variables y la simplificación de su ejecución manteniendo los resultados obtenidos próximos al comportamiento real de la estructura. Se propone el empleo de un método de reducción que encaja en la tercera de las opciones y consiste en un análisis simplificado que proporciona una solución para la respuesta dinámica de una estructura en el subespacio modal complejo empleando un volumen de datos muy reducido. De este modo se pueden realizar análisis paramétricos variando múltiples parámetros, para obtener una solución muy aproximada al objetivo buscado. Se propone no solo la variación de propiedades locales de masa, rigidez y amortiguamiento sino la adición de grados de libertad a la estructura original para el cálculo de la respuesta tanto permanente como transitoria. Adicionalmente, su facilidad de implementación permite un control exhaustivo sobre las variables del problema y la implementación de mejoras como diferentes formas de obtención de los autovalores o la eliminación de las limitaciones de amortiguamiento en la estructura original. El objetivo del método se puede considerar similar a los que se obtienen al aplicar el método de Guyan u otras técnicas de reducción de modelos empleados en dinámica estructural. Sin embargo, aunque el método permite ser empleado en conjunción con otros para obtener las ventajas de ambos, el presente procedimiento no realiza la condensación del sistema de ecuaciones, sino que emplea la información del sistema de ecuaciones completa estudiando tan solo la respuesta en las variables apropiadas de los puntos de interés para el analista. Dicho interés puede surgir de la necesidad de obtener la respuesta de las grandes estructuras en unos puntos determinados o de la necesidad de modificar la estructura en zonas determinadas para cambiar su comportamiento (respuesta en aceleraciones, velocidades o desplazamientos) ante cargas dinámicas. Por lo tanto, el procedimiento está particularmente indicado para la selección del valor óptimo de varios parámetros en grandes estructuras (del orden de cientos de miles de modos) como pueden ser la localización de elementos introducidos, rigideces, masas o valores de amortiguamientos viscosos en estudios previos en los que diversas soluciones son planteadas y optimizadas, y que en el caso de grandes estructuras, pueden conllevar un número de simulaciones extremadamente elevado para alcanzar la solución óptima. Tras plantear las herramientas necesarias y desarrollar el procedimiento, se propone un caso de estudio para su aplicación al modelo de elementos .nitos del UAV MILANO desarrollado por el Instituto Nacional de Técnica Aeroespacial. A dicha estructura se le imponen ciertos requisitos al incorporar un equipo en aceleraciones en punta de ala izquierda y desplazamientos en punta de ala derecha en presencia de la sustentación producida por una ráfaga continua de viento de forma sinusoidal. La modificación propuesta consiste en la adición de un equipo en la punta de ala izquierda, bien mediante un anclaje rígido, bien unido mediante un sistema de reducción de la respuesta dinámica con propiedades de masa, rigidez y amortiguamiento variables. El estudio de los resultados obtenidos permite determinar la optimización de los parámetros del sistema de atenuación por medio de múltiples análisis dinámicos de forma que se cumplan de la mejor forma posible los requisitos impuestos con la modificación. Se comparan los resultados con los obtenidos mediante el uso de un programa comercial de análisis por el método de los elementos .nitos lográndose soluciones muy aproximadas entre el modelo completo y el reducido. La influencia de diversos factores como son el amortiguamiento modal de la estructura original, el número de modos retenidos en la truncatura o la precisión proporcionada por el barrido en frecuencia se analiza en detalle para, por último, señalar la eficiencia en términos de tiempo y volumen de datos de computación que ofrece el método propuesto en comparación con otras aproximaciones. Por lo tanto, puede concluirse que el método propuesto se considera una opción útil y eficiente para el análisis paramétrico de modificaciones locales en grandes estructuras. ABSTRACT When developing structural design and analysis some projects require dynamic analysis of large finite element models with millions of degrees of freedom which use large size data .les. The analysis complexity and size grow if a parametric analysis is required. This problem has been approached traditionally in several ways: one way is increasing the power and the storage capacity of computer systems involved in the analysis. Other obvious way is reducing the total amount of analyses and their details. Finally, complementary methods to finite element analysis can also be employed in order to limit the number of variables and to reduce the execution time keeping the results as close as possible to the actual behaviour of the structure. Following this third option, we propose a model reduction method that is based in a simplified analysis that supplies a solution for the dynamic response of the structure in the complex modal space using few data. Thereby, parametric analysis can be done varying multiple parameters so as to obtain a solution which complies with the desired objetive. We propose not only mass, stiffness and damping variations, but also addition of degrees of freedom to the original structure in order to calculate the transient and steady-state response. Additionally, the simple implementation of the procedure allows an in-depth control of the problem variables. Furthermore, improvements such as different ways to obtain eigenvectors or to remove damping limitations of the original structure are also possible. The purpose of the procedure is similar to that of using the Guyan or similar model order reduction techniques. However, in our method we do not perform a true model order reduction in the traditional sense. Furthermore, additional gains, which we do not explore herein, can be obtained through the combination of this method with traditional model-order reduction procedures. In our procedure we use the information of the whole system of equations is used but only those nodes of interest to the analyst are processed. That interest comes from the need to obtain the response of the structure at specific locations or from the need to modify the structure at some suitable positions in order to change its behaviour (acceleration, velocity or displacement response) under dynamic loads. Therefore, the procedure is particularly suitable for parametric optimization in large structures with >100000 normal modes such as position of new elements, stiffness, mass and viscous dampings in previous studies where different solutions are devised and optimized, and in the case of large structures, can carry an extremely high number of simulations to get the optimum solution. After the introduction of the required tools and the development of the procedure, a study case is proposed with use the finite element model (FEM) of the MILANO UAV developed by Instituto Nacional de Técnica Aeroespacial. Due to an equipment addition, certain acceleration and displacement requirements on left wing tip and right wing tip, respectively, are imposed. The structure is under a continuous sinusoidal wind gust which produces lift. The proposed modification consists of the addition of an equipment in left wing tip clamped through a rigid attachment or through a dynamic response reduction system with variable properties of mass, stiffness and damping. The analysis of the obtained results allows us to determine the optimized parametric by means of multiple dynamic analyses in a way such that the imposed requirements have been accomplished in the best possible way. The results achieved are compared with results from a commercial finite element analysis software, showing a good correlation. Influence of several factors such as the modal damping of the original structure, the number of modes kept in the modal truncation or the precission given by the frequency sweep is analyzed. Finally, the efficiency of the proposed method is addressed in tems of computational time and data size compared with other approaches. From the analyses performed, we can conclude that the proposed method is a useful and efficient option to perform parametric analysis of possible local modifications in large structures.

Ultimate strength of continuous composite beams in combined bending and shear

The contributions of the concrete slab and composite action to the vertical shear strength of continuous steel-concrete composite beams are ignored in current design codes, which result in conservative designs. This paper investigates the ultimate strength of continuous composite beams in combined bending and shear by using the finite element analysis method. A three-dimensional finite element model has been developed to account for the geometric and material nonlinear behaviour of continuous composite beams. The finite element model is verified by experimental results and then used to study the effects of the concrete slab and shear connection on the vertical shear strength. The moment-shear interaction strength of continuous composite beams is also investigated by varying the moment/ shear ratio. It is shown that the concrete slab and composite action significantly increase the ultimate strength of continuous composite beams. Based on numerical results, design models are proposed for the vertical shear strength and moment-shear interaction of continuous composite beams. The proposed design models, which incorporates the effects of the concrete slab, composite action, stud pullout failure and web shear buckling, are compared with experimental results with good agreement. (C) 2003 Elsevier Ltd. All rights reserved.

Comparison of beam theory and finite-element analysis with in vivo bone strain data from the alligator cranium

The mechanical behavior of the vertebrate skull is often modeled using free-body analysis of simple geometric structures and, more recently, finite-element (FE) analysis. In this study, we compare experimentally collected in vivo bone strain orientations and magnitudes from the cranium of the American alligator with those extrapolated from a beam model and extracted from an FE model. The strain magnitudes predicted from beam and FE skull models bear little similarity to relative and absolute strain magnitudes recorded during in vivo biting experiments. However, quantitative differences between principal strain orientations extracted from the FE skull model and recorded during the in vivo experiments were smaller, and both generally matched expectations from the beam model. The differences in strain magnitude between the data sets may be attributable to the level of resolution of the models, the material properties used in the FE model, and the loading conditions (i.e., external forces and constraints). This study indicates that FE models and modeling of skulls as simple engineering structures may give a preliminary idea of how these structures are loaded, but whenever possible, modeling results should be verified with either in vitro or preferably in vivo testing, especially if precise knowledge of strain magnitudes is desired. (c) 2005 Wiley-Liss, Inc.

Biomechanics of the rostrum in crocodilians: A comparative analysis using finite-element modeling

This article reports the use of simple beam and finite-element models to investigate the relationship between rostral shape and biomechanical performance in living crocodilians under a range of loading conditions. Load cases corresponded to simple biting, lateral head shaking, and twist feeding behaviors. The six specimens were chosen to reflect, as far as possible, the full range of rostral shape in living crocodilians: a juvenile Caiman crocodilus, subadult Alligator mississippiensis and Crocodylus johnstoni, and adult Caiman crocodilus, Melanosuchus niger, and Paleosuchus palpebrosus. The simple beam models were generated using morphometric landmarks from each specimen. Three of the finite-element models, the A. mississippiensis, juvenile Caiman crocodilus, and the Crocodylus johnstoni, were based on CT scan data from respective specimens, but these data were not available for the other models and so these-the adult Caiman crocodilus, M. niger, and P. palpebrosus-were generated by morphing the juvenile Caiman crocodilus mesh with reference to three-dimensional linear distance measured from specimens. Comparison of the mechanical performance of the six finite-element models essentially matched results of the simple beam models: relatively tall skulls performed best under vertical loading and tall and wide skulls performed best under torsional loading. The widely held assumption that the platyrostral (dorsoventrally flattened) crocodilian skull is optimized for torsional loading was not supported by either simple beam theory models or finite-element modeling. Rather than being purely optimized against loads encountered while subduing and processing food, the shape of the crocodilian rostrum may be significantly affected by the hydrodynamic constraints of catching agile aquatic prey. This observation has important implications for our understanding of biomechanics in crocodilians and other aquatic reptiles.

Sviluppo e applicazione di nuovi criteri per la verifica a fatica multiassiale di componenti strutturali in leghe di acciaio e alluminio

La presente Tesi ha per oggetto lo sviluppo e la validazione di nuovi criteri per la verifica a fatica multiassiale di componenti strutturali metallici . In particolare, i nuovi criteri formulati risultano applicabili a componenti metallici, soggetti ad un’ampia gamma di configurazioni di carico: carichi multiassiali variabili nel tempo, in modo ciclico e random, per alto e basso/medio numero di cicli di carico. Tali criteri costituiscono un utile strumento nell’ambito della valutazione della resistenza/vita a fatica di elementi strutturali metallici, essendo di semplice implementazione, e richiedendo tempi di calcolo piuttosto modesti. Nel primo Capitolo vengono presentate le problematiche relative alla fatica multiassiale, introducendo alcuni aspetti teorici utili a descrivere il meccanismo di danneggiamento a fatica (propagazione della fessura e frattura finale) di componenti strutturali metallici soggetti a carichi variabili nel tempo. Vengono poi presentati i diversi approcci disponibili in letteratura per la verifica a fatica multiassiale di tali componenti, con particolare attenzione all'approccio del piano critico. Infine, vengono definite le grandezze ingegneristiche correlate al piano critico, utilizzate nella progettazione a fatica in presenza di carichi multiassiali ciclici per alto e basso/medio numero di cicli di carico. Il secondo Capitolo è dedicato allo sviluppo di un nuovo criterio per la valutazione della resistenza a fatica di elementi strutturali metallici soggetti a carichi multiassiali ciclici e alto numero di cicli. Il criterio risulta basato sull'approccio del piano critico ed è formulato in termini di tensioni. Lo sviluppo del criterio viene affrontato intervenendo in modo significativo su una precedente formulazione proposta da Carpinteri e collaboratori nel 2011. In particolare, il primo intervento riguarda la determinazione della giacitura del piano critico: nuove espressioni dell'angolo che lega la giacitura del piano critico a quella del piano di frattura vengono implementate nell'algoritmo del criterio. Il secondo intervento è relativo alla definizione dell'ampiezza della tensione tangenziale e un nuovo metodo, noto come Prismatic Hull (PH) method (di Araújo e collaboratori), viene implementato nell'algoritmo. L'affidabilità del criterio viene poi verificata impiegando numerosi dati di prove sperimentali disponibili in letteratura. Nel terzo Capitolo viene proposto un criterio di nuova formulazione per la valutazione della vita a fatica di elementi strutturali metallici soggetti a carichi multiassiali ciclici e basso/medio numero di cicli. Il criterio risulta basato sull'approccio del piano critico, ed è formulato in termini di deformazioni. In particolare, la formulazione proposta trae spunto, come impostazione generale, dal criterio di fatica multiassiale in regime di alto numero di cicli discusso nel secondo Capitolo. Poiché in presenza di deformazioni plastiche significative (come quelle caratterizzanti la fatica per basso/medio numero di cicli di carico) è necessario conoscere il valore del coefficiente efficace di Poisson del materiale, vengono impiegate tre differenti strategie. In particolare, tale coefficiente viene calcolato sia per via analitica, che per via numerica, che impiegando un valore costante frequentemente adottato in letteratura. Successivamente, per validarne l'affidabilità vengono impiegati numerosi dati di prove sperimentali disponibili in letteratura; i risultati numerici sono ottenuti al variare del valore del coefficiente efficace di Poisson. Inoltre, al fine di considerare i significativi gradienti tensionali che si verificano in presenza di discontinuità geometriche, come gli intagli, il criterio viene anche esteso al caso dei componenti strutturali intagliati. Il criterio, riformulato implementando il concetto del volume di controllo proposto da Lazzarin e collaboratori, viene utilizzato per stimare la vita a fatica di provini con un severo intaglio a V, realizzati in lega di titanio grado 5. Il quarto Capitolo è rivolto allo sviluppo di un nuovo criterio per la valutazione del danno a fatica di elementi strutturali metallici soggetti a carichi multiassiali random e alto numero di cicli. Il criterio risulta basato sull'approccio del piano critico ed è formulato nel dominio della frequenza. Lo sviluppo del criterio viene affrontato intervenendo in modo significativo su una precedente formulazione proposta da Carpinteri e collaboratori nel 2014. In particolare, l’intervento riguarda la determinazione della giacitura del piano critico, e nuove espressioni dell'angolo che lega la giacitura del piano critico con quella del piano di frattura vengono implementate nell'algoritmo del criterio. Infine, l’affidabilità del criterio viene verificata impiegando numerosi dati di prove sperimentali disponibili in letteratura.

Application of the finite element method to problems in fracture mechanics

Numerical techniques have been finding increasing use in all aspects of fracture mechanics, and often provide the only means for analyzing fracture problems. The work presented here, is concerned with the application of the finite element method to cracked structures. The present work was directed towards the establishment of a comprehensive two-dimensional finite element, linear elastic, fracture analysis package. Significant progress has been made to this end, and features which can now be studied include multi-crack tip mixed-mode problems, involving partial crack closure. The crack tip core element was refined and special local crack tip elements were employed to reduce the element density in the neighbourhood of the core region. The work builds upon experience gained by previous research workers and, as part of the general development, the program was modified to incorporate the eight-node isoparametric quadrilateral element. Also. a more flexible solving routine was developed, and provided a very compact method of solving large sets of simultaneous equations, stored in a segmented form. To complement the finite element analysis programs, an automatic mesh generation program has been developed, which enables complex problems. involving fine element detail, to be investigated with a minimum of input data. The scheme has proven to be versati Ie and reasonably easy to implement. Numerous examples are given to demonstrate the accuracy and flexibility of the finite element technique.

Modelling the human accommodation system using finite element analysis

The human accommodation system has been extensively examined for over a century, with a particular focus on trying to understand the mechanisms that lead to the loss of accommodative ability with age (Presbyopia). The accommodative process, along with the potential causes of presbyopia, are disputed; hindering efforts to develop methods of restoring accommodation in the presbyopic eye. One method that can be used to provide insight into this complex area is Finite Element Analysis (FEA). The effectiveness of FEA in modelling the accommodative process has been illustrated by a number of accommodative FEA models developed to date. However, there have been limitations to these previous models; principally due to the variation in data on the geometry of the accommodative components, combined with sparse measurements of their material properties. Despite advances in available data, continued oversimplification has occurred in the modelling of the crystalline lens structure and the zonular fibres that surround the lens. A new accommodation model was proposed by the author that aims to eliminate these limitations. A novel representation of the zonular structure was developed, combined with updated lens and capsule modelling methods. The model has been designed to be adaptable so that a range of different age accommodation systems can be modelled, allowing the age related changes that occur to be simulated. The new modelling methods were validated by comparing the changes induced within the model to available in vivo data, leading to the definition of three different age models. These were used in an extended sensitivity study on age related changes, where individual parameters were altered to investigate their effect on the accommodative process. The material properties were found to have the largest impact on the decline in accommodative ability, in particular compared to changes in ciliary body movement or zonular structure. Novel data on the importance of the capsule stiffness and thickness was also established. The new model detailed within this thesis provides further insight into the accommodation mechanism, as well as a foundation for future, more detailed investigations into accommodation, presbyopia and accommodative restoration techniques.

Alternatives to steel grid bridge decks

Most of the moveable bridges use open grid steel decks, because these are factory assembled, light-weight, and easy to install. Open grid steel decks, however, are not as skid resistant as solid decks. Costly maintenance, high noise levels, poor riding comfort and susceptibility to vibrations are among the other disadvantages of these decks. The major objective of this research was to develop alternative deck systems which weigh no more than 25 lb/ft2, have solid riding surface, are no more than 4–5 in. thick and are able to withstand prescribed loading. Three deck systems were considered in this study: ultra-high performance concrete (UHPC) deck, aluminum deck and UHPC-fiber reinforced polymer (FRP) tube deck. UHPC deck was the first alternative system developed as a part of this project. Due to its ultra high strength, this type of concrete results in thinner sections, which helps satisfy the strict self-weight limit. A comprehensive experimental and analytical evaluation of the system was carried out to establish its suitability. Both single and multi-unit specimens with one or two spans were tested for static and dynamic loading. Finite element models were developed to predict the deck behavior. The study led to the conclusion that the UHPC bridge deck is a feasible alternative to open grid steel deck. Aluminum deck was the second alternative system studied in this project. A detailed experimental and analytical evaluation of the system was carried out. The experimental work included static and dynamic loading on the deck panels and connections. Analytical work included detailed finite element modeling. Based on the in-depth experimental and analytical evaluations, it was concluded that aluminum deck was a suitable alternative to open grid steel decks and is ready for implementation. UHPC-FRP tube deck was the third system developed in this research. Prestressed hollow core decks are commonly used, but the proposed type of steel-free deck is quite novel. Preliminary experimental evaluations of two simple-span specimens, one with uniform section and the other with tapered section were carried out. The system was shown to have good promise to replace the conventional open grid decks. Additional work, however, is needed before the system is recommended for field application.

Development of thresholds for surface preparation in repair systems with wet lay -up fiber reinforced polymers

Bonded repair of concrete structures with fiber reinforced polymer (FRP) systems is increasingly being accepted as a cost-efficient and structurally viable method of rapid rehabilitation of concrete structures. However, the relationships between long-term performance attributes, service-life, and details of the installation process are not easy to quantify. Accordingly, there is currently a lack of generally accepted construction specifications, making it difficult for the field engineer to certify the adequacy of the construction process. ^ The objective of the present study, as part of the National Cooperative Highway Research Program (NCHRP) Project 10-59B, was to investigate the effect of surface preparation on the behavior of wet lay-up FRP repair systems and consequently develop rational thresholds that provide sufficient performance. ^ The research program was comprised of both experimental and analytical work for wet lay-up FRP applications. The experimental work included flexure testing of sixty-seven (67) reinforced concrete beams and bond testing of ten (10) reinforced concrete blocks. Four different parameters were studied: surface roughness, surface flatness, surface voids and bug holes, and surface cracks/cuts. The findings were analyzed from various aspects and compared with the data available in the literature. As part of the analytical work, finite element models of the flexural specimens with surface flaws were developed using ANSYS. The purpose of this part was to extend the parametric study on the effects of concrete surface flaws and verify the experimental results based on nonlinear finite element analysis. ^ Test results showed that surface roughness does not appear to have a significant influence on the overall performance of the wet lay-up FRP systems with or without adequate anchorage, and whether failure was by debonding or rupture of FRP. Both experimental and analytical results for surface flatness proved that peaks on concrete surface, in the range studied, do not have a significant effect on the performance of wet lay-up FRP systems. However, valleys of particular size could reduce the strength of wet lay-up FRP systems. Test results regarding surface voids and surface cracks/cuts revealed that previously suggested thresholds for these flaws appear to be conservative, as also confirmed by analytical study. ^

Criação de zonas de alta permeabilidade em fraturas: modelagem numérica do Lajedo de Soledade, Formação Jandaíra, Bacia Potiguar

The fracturing in carbonate rocks has been attracting increasingly attention due to new oil discoveries in carbonate reservoirs. This study investigates how the fractures (faults and joints) behave when subjected to different stress fields and how their behavior may be associated with the generation of karst and consequently to increased secondary porosity in these rocks. In this study I used satellite imagery and unmanned aerial vehicle UAV images and field data to identify and map faults and joints in a carbonate outcrop, which I consider a good analogue of carbonate reservoir. The outcrop comprises rocks of the Jandaíra Formation, Potiguar Basin. Field data were modeled using the TECTOS software, which uses finite element analysis for 2D fracture modeling. I identified three sets of fractures were identified: NS, EW and NW-SE. They correspond to faults that reactivate joint sets. The Ratio of Failure by Stress (RFS) represents stress concentration and how close the rock is to failure and reach the Mohr-Coulomb envelopment. The results indicate that the tectonic stresses are concentrated in preferred structural zones, which are ideal places for carbonate dissolution. Dissolution was observed along sedimentary bedding and fractures throughout the outcrop. However, I observed that the highest values of RFS occur in fracture intersections and terminations. These are site of karst concentration. I finally suggest that there is a relationship between stress concentration and location of karst dissolution in carbonate rocks.

Criação de zonas de alta permeabilidade em fraturas: modelagem numérica do Lajedo de Soledade, Formação Jandaíra, Bacia Potiguar

The fracturing in carbonate rocks has been attracting increasingly attention due to new oil discoveries in carbonate reservoirs. This study investigates how the fractures (faults and joints) behave when subjected to different stress fields and how their behavior may be associated with the generation of karst and consequently to increased secondary porosity in these rocks. In this study I used satellite imagery and unmanned aerial vehicle UAV images and field data to identify and map faults and joints in a carbonate outcrop, which I consider a good analogue of carbonate reservoir. The outcrop comprises rocks of the Jandaíra Formation, Potiguar Basin. Field data were modeled using the TECTOS software, which uses finite element analysis for 2D fracture modeling. I identified three sets of fractures were identified: NS, EW and NW-SE. They correspond to faults that reactivate joint sets. The Ratio of Failure by Stress (RFS) represents stress concentration and how close the rock is to failure and reach the Mohr-Coulomb envelopment. The results indicate that the tectonic stresses are concentrated in preferred structural zones, which are ideal places for carbonate dissolution. Dissolution was observed along sedimentary bedding and fractures throughout the outcrop. However, I observed that the highest values of RFS occur in fracture intersections and terminations. These are site of karst concentration. I finally suggest that there is a relationship between stress concentration and location of karst dissolution in carbonate rocks.

Simplified PRBMs of spatial compliant multi-beam modules for planar motion

PRBMs (pseudo-rigid-body models) have been becoming important engineering technologies/methods in the field of compliant mechanisms to simplify the design and analysis through the use of the knowledge body of rigid-body mechanisms coupling with springs. This article addresses the PRBMs of spatial multi-beam modules for planar motion, which are composed of three or more symmetrical wire/slender beams parallel to each other where the planar twisting DOF (degree of freedom) is assumed to be very small for specific applications/loading conditions. Simplified PRBMs are firstly proposed through replacing each beam in spatial multi-beam module with a rigid-body link plus two identical spherical joints at its two ends. The characteristics factor, bending stiffness and twisting stiffness for the spherical joint are determined. Load-displacement equations are then derived for a class of spatial multi-beam modules and general spatial multi-beam modules using the virtual work principle and kinematic relationships. Finally, nonlinear FEA (finite element analysis) is employed with comparisons with the PRBMs. The present PRBMs have shown the ability to predict the primary nonlinear constraint characteristics such as load-stiffening effect, cross-axis coupling in the two primary translational directions and buckling load.

Vibration based health monitoring of civil infrastructure using energy harvesting techniques

The use of structural health monitoring of civil structures is ever expanding and by assessing the dynamical condition of structures, informed maintenance management can be conducted at both individual and network levels. With the continued growth of information age technology, the potential arises for smart monitoring systems to be integrated with civil infrastructure to provide efficient information on the condition of a structure. The focus of this thesis is the integration of smart technology with civil infrastructure for the purposes of structural health monitoring. The technology considered in this regard are devices based on energy harvesting materials. While there has been considerable focus on the development and optimisation of such devices using steady state loading conditions, their applications for civil infrastructure are less known. Although research is still in initial stages, studies into the uses associated with such applications are very promising. Through the use of the dynamical response of structures to a variety of loading conditions, the energy harvesting outputs from such devices is established and the potential power output determined. Through a power variance output approach, damage detection of deteriorating structures using the energy harvesting devices is investigated. Further applications of the integration of energy harvesting devices with civil infrastructure investigated by this research includes the use of the power output as a indicator for control. Four approaches are undertaken to determine the potential applications arising from integrating smart technology with civil infrastructure, namely • Theoretical analysis to determine the applications of energy harvesting devices for vibration based health monitoring of civil infrastructure. • Laboratory experimentation to verify the performance of different energy harvesting configurations for civil infrastructure applications. • Scaled model testing as a method to experimentally validate the integration of the energy harvesting devices with civil infrastructure. • Full scale deployment of energy harvesting device with a bridge structure. These four approaches validate the application of energy harvesting technology with civil infrastructure from a theoretical, experimental and practical perspective.