Simulação do comportamento não linear de materiais quase-frágeis via elementos finitos com fissura embebida

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Rietveld analyses and piezoelectric properties of niobium doped bismuth titanate systems

Bi 4Ti 3- xNbxO 12 (BITNb) samples, with × ranging from 0 to 0.40 were obtained using a polymeric precursor solution. Rietveld analyses confirmed that the powders crystallize in an orthorhombic structure free of secondary phases with space group Fmmm. Raman analysis evidenced a sharp increase in the bands intensity located at 129 cm -1 and 190 cm -1 due the lattice distortion in BIT02Nb and BIT04Nb compositions. UV-vis spectra indicated that addition of niobium causes a reduction of defects in the BIT lattice due the suppression of oxygen vacancies located at BO-6 octahedral. Size and morphology of particles as well as electrical behavior of BIT ceramics were affected by addition of donor dopant. Polarization reversal was investigated by applying dc voltage through a conductive tip during the area scanning and was investigated by piezoresponse force microscopy (PFM). PFM measurements revealed a decrease in piezoelectric response with increasing Nb concentration originating from a reduced polarizability along the a-axis. High spontaneous polarization is noted for the less doped sample due the reduction of strain energy and pin charged defects after niobium addition. Copyright © 2010 American Scientific Publishers.

A finite element approach for predicting the ultimate rotation capacity of RC beams

This paper presents a numerical approach to model the complex failure mechanisms that define the ultimate rotational capacity of reinforced concrete beams. The behavior in tension and compression is described by a constitutive damage model derived from a combination of two specific damage models [1]. The nonlinear behavior of the compressed region is treated by the compressive damage model based on the Drucker-Prager criterion written in terms of the effective stresses. The tensile damage model employs a failure criterion based on the strain energy associated with the positive part the effective stress tensor. This model is used to describe the behavior of very thin bands of strain localization, which are embedded in finite elements to represent multiple cracks that occur in the tensioned region [2]. The softening law establishes dissipation energy compatible with the fracture energy of the concrete. The reinforcing steel bars are modeled by truss elements with elastic-perfect plastic behavior. It is shown that the resulting approach is able to predict the different stages of the collapse mechanism of beams with distinct sizes and reinforcement ratios. The tensile damage model and the finite element embedded crack approach are able to describe the stiffness reduction due to concrete cracking in the tensile zone. The truss elements are able to reproduce the effects of steel yielding and, finally, the compressive damage model is able to describe the non-linear behavior of the compressive zone until the complete collapse of the beam due to crushing of concrete. The proposed approach is able to predict well the plastic rotation capacity of tested beams [3], including size-scale effects.

Um modelo constitutivo de dano composto para simular o comportamento de materiais quase-frágeis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Controle ativo de estruturas reticuladas utilizando atuadores de pilha de pzt

Pós-graduação em Engenharia Mecânica - FEIS

Platinum-Based Nanoalloys PtnTM55-n (TM=Co, Rh, Au): A Density Functional Theory Investigation

Structural and electronic properties of the PtnTM55-n (TM = Co, Rh, Au) nanoalloys are investigated using density functional theory within the generalized gradient approximation and employing the all-electron projected augmented wave method. For TM = Co and Rh, the excess energy, which measures the relative energy stability of the nanoalloys, is negative for all Pt compositions. We found that the excess energy has similar values for a wide range of Pt compositions, i.e., n = 20-42 and n = 28-42 for Co and Rh, respectively, with the core shell icosahedron-like configuration (n = 42) being slightly more stable for both Co and Rh systems because of the larger release of the strain energy due to the smaller atomic size of the Co and Rh atoms. For TM = Au, the excess energy is positive for all compositions, except for n = 13, which is energetically favorable due to the formation of the core-shell structure (Pt in the core and Au atoms at the surface). Thus, our calculations confirm that the formation of core-shell structures plays an important role to increase the stability of nanoalloys. The center of gravity of the occupied d-states changes almost linearly as a function of the Pt composition, and hence, based on the d-band model, the magnitude of the adsorption energy of an adsorbate can be tuned by changing the Pt composition. The magnetic moments of PtnCo55-n decrease almost linearly as a function of the Pt composition; however, the same does not hold for PtRh and PtAu. We found an enhancement of the magnetic moments of PtRh by a few times by increasing Pt composition, which we explain by the compression effects induced by the large size of the Pt atoms compared with the Rh atoms.

Experimental investigation into the effect of stress on dissolution and growth of very soluble brittle salts in aqueous solution

Das Studium der Auflösungs- und Wachstumsprozesse an Feststoff-Flüssigkeits-Grenzflächen unter nicht-hydrostatischen Beanspruchungen ist wesentlich für das Verständnis von Defor-mationsprozessen, die in der Erde ablaufen. Unter diesen genannten Prozessen gehört die Drucklösung zu den wichtigsten duktilen Deformationsprozessen, von der Diagenese bishin zur niedrig- bis mittelgradigen metamorphen Bedingungen. Bisher ist allerdings wenig darüber bekannt, welche mechanischen, physikalischen oder chemischen Potentialenergie-Gradienten die Drucklösung steuern. I.a. wird angenommen, daß die Drucklösung durch Un-terschiede kristallplastischer Verformungsenergien oder aber durch Unterschiede der Normal-beanspruchung an Korngrenzen gesteuert wird. Unterschiede der elastischen Verformungs-energien werden dabei allerdings als zu gering erachtet, um einen signifikanten Beitrag zu leisten. Aus diesem Grund werden sie als mögliche treibende Kräfte für die Drucklösung vernachlässigt. Andererseits haben neue experimentelle und theoretische Untersuchungen gezeigt, daß die elastische Verformung in der Tat einen starken Einfluß auf Lösungs- und Wachstumsmechanismen von Kristallen in einer Lösung haben kann. Da die in der Erdkruste vorherrschenden Deformationsmechanismen überwiegend im elastischen Verformungsbereich der Gesteine ablaufen, ist es sehr wichtig, das Verständnis für die Effekte, die die elastische Verformung verursacht, zu erweitern, und ihre Rolle während der Deformation durch Drucklösung zu definieren. Die vorliegende Arbeit beschäftigt sich mit Experimenten, bei denen der Effekt der mechanisch kompressiven Beanspruchung auf Lösungs- und Wachstumsprozesse von Einzelkristallen unterschiedlicher, sehr gut löslicher, elastisch/spröder Salze untersucht wurde. Diese Salze wurden als Analoga gesteinsbildender Minerale wie Quarz und Calcit ausgewählt. Der Einfluß von Stress auf die Ausbildung der Oberflächenmikrostrukturen in einer untersättigten Lösung wurde an Kaliumalaun untersucht.Lösungsrillen (20 – 40 µm breit, 10 – 40 µm tief und 20 – 80 µm Abstand) entwickelten sich in den Bereichen, in denen die Beanspruchung im Kristall am größten war. Sie verschwanden wieder, sobald der Kristall entlastet wurde. Diese Rillen entwickelten sich parallel zu niedrig indizierten kristallographischen Richtungen und sub-perpendikular zu den Trajektorien, die der maximalen, lokalen kompressiven Beanspruchung entsprachen. Die Größe der Lösungsrillen hing von der lokalen Oberflächenbeanspruchung, der Oberflächenenergie und dem Untersättigungsgrad der wässrigen Lösung ab. Die mikrostrukturelle Entwicklung der Kristalloberflächen stimmte gut mit den theoretischen Vorhersagen überein, die auf den Modellen von Heidug & Leroy (1994) und Leroy & Heidug (1994) basieren. Der Einfluß der Beanspruchung auf die Auflösungsrate wurde an Natriumchlorat-Einzelkristallen untersucht. Dabei wurde herausgefunden, daß sich gestresste Kristalle schneller lösen als Kristalle, auf die keine Beanspruchung einwirkt. Der experimentell beobachtete Anstieg der Auflösungsrate der gestressten Kristalle war ein bis zwei Größenordnungen höher als theoretisch erwartet. Die Auflösungsrate stieg linear mit dem Stress an, und der Anstieg war um so größer, je stärker die Lösung untersättigt war. Außerdem wurde der Effekt der Bean-spruchung auf das Kristallwachstum an Kaliumalaun- und Kaliumdihydrogenphosphat-Ein-zelkristallen untersucht. Die Wachstumsrate der Flächen {100} und {110} von Kalium-alaun war bei Beanspruchung stark reduziert. Für all diese Ergebnisse spielte die Oberflächenrauhigkeit der Kristalle eine Schlüsselrolle, indem sie eine nicht-homogene Stressverteilung auf der Kristalloberfläche verursachte. Die Resultate zeigen, daß die elastische Verformung eine signifikante Rolle während der Drucklösung spielen kann, und eine signifikante Deformation in der oberen Kruste verursachen kann, bei Beanspruchungen, die geringer sind, als gemeinhin angenommen wird. Somit folgt, daß die elastische Bean-spruchung berücksichtigt werden muß, wenn mikrophysikalische Deformationsmodelle entwickelt werden sollen.

Influence of adhesive thickness on adhesively bonded joints under fatigue loading

Research on adhesive joints is arousing increasing interest in aerospace industry. Incomplete knowledge of fatigue in adhesively bonded joints is a major obstacle to their application. The prediction of the disbonding growth is yet an open question. This thesis researches the influence of the adhesive thickness on fatigue disbond growth. Experimental testing on specimens with different thickness has been performed. Both a conventional approach based on the strain energy release rate and an approach based on cyclic strain energy are provided. The inadequacy of the former approach is discussed. Outcomes from tests support the idea of correlating the crack growth rate to the cyclic strain energy. In order to push further the study, a 2D finite element model for the prediction of disbond growth under quasi-static loading has been developed and implemented in Abaqus. Numerical simulations have been conducted with different values of the adhesive thickness. The results from tests and simulations are in accordance with each other. According to them, no dependence of disbonding on the adhesive thickness has been evidenced.

Numerical model of the human cornea based on stromal microstructure: identification of mechanical properties for two age groups

The optical quality of the human eye mainly depends on the refractive performance of the cornea. The shape of the cornea is a mechanical balance between intraocular pressure and tissue intrinsic stiffness. Several surgical procedures in ophthalmology alter the biomechanics of the cornea to provoke local or global curvature changes for vision correction. Legitimated by the large number of surgical interventions performed every day, the demand for a deeper understanding of corneal biomechanics is rising to improve the safety of procedures and medical devices. The aim of our work is to propose a numerical model of corneal biomechanics, based on the stromal microstructure. Our novel anisotropic constitutive material law features a probabilistic weighting approach to model collagen fiber distribution as observed on human cornea by Xray scattering analysis (Aghamohammadzadeh et. al., Structure, February 2004). Furthermore, collagen cross-linking was explicitly included in the strain energy function. Results showed that the proposed model is able to successfully reproduce both inflation and extensiometry experimental data (Elsheikh et. al., Curr Eye Res, 2007; Elsheikh et. al., Exp Eye Res, May 2008). In addition, the mechanical properties calculated for patients of different age groups (Group A: 65-79 years; Group B: 80-95 years) demonstrate an increased collagen cross-linking, and a decrease in collagen fiber elasticity from younger to older specimen. These findings correspond to what is known about maturing fibrous biological tissue. Since the presented model can handle different loading situations and includes the anisotropic distribution of collagen fibers, it has the potential to simulate clinical procedures involving nonsymmetrical tissue interventions. In the future, such mechanical model can be used to improve surgical planning and the design of next generation ophthalmic devices.

Personalized X-Ray Reconstruction of the Proximal Femur via Intensity-Based Non-rigid 2D-3D Registration

This paper presents a new approach for reconstructing a patient-specific shape model and internal relative intensity distribution of the proximal femur from a limited number (e.g., 2) of calibrated C-arm images or X-ray radiographs. Our approach uses independent shape and appearance models that are learned from a set of training data to encode the a priori information about the proximal femur. An intensity-based non-rigid 2D-3D registration algorithm is then proposed to deformably fit the learned models to the input images. The fitting is conducted iteratively by minimizing the dissimilarity between the input images and the associated digitally reconstructed radiographs of the learned models together with regularization terms encoding the strain energy of the forward deformation and the smoothness of the inverse deformation. Comprehensive experiments conducted on images of cadaveric femurs and on clinical datasets demonstrate the efficacy of the present approach.

Quantum-Mechanical Investigation of Large Water Clusters

The PM3 quantum-mechanical method has been used to study large water clusters ranging from 8 to 42 water molecules. These large clusters are built from smaller building blocks. The building blocks include cyclic tetramers, pentamers, octamers, and a pentagonal dodecahedron cage. The correlations between the strain energy resulting from bending of the hydrogen bonds formed by different cluster motifs and the number of waters involved in the cluster are discussed. The PM3 results are compared with TIP4P potential and ab initio results. The number of net hydrogen bonds per water increases with the cluster size. This places a limit on the size of clusters that would fit the Benson model of liquid water. Many of the 20-mer clusters fit the Benson model well. Calculations of the ion cluster (H20)4o(H30+)2 reveal that the m/e ratio obtainable by mass spectrometry experiments can uniquely indicate the conformation of the 20 water pentagonal dodecahedron cage present in the larger clusters.

Efficient and Accurate Characterization of the Bergman Cyclization for an Expanded Substructure of Esperamicin A1

Incorporation of enediynes into anticancer drugs remains an intriguing yet elusive strategy for the design of therapeutically active agents. Density functional theory was used to locate reactants, products, and transition states along the Bergman cyclization pathways connecting enediynes to reactive para-biradicals. Sum method correction to low-level calculations confirmed B3LYP/6-31G(d,p) as the method of choice in investigating enediynes. Herein described as MI:Sum, calculated reaction enthalpies differed from experiment by an average of 2.1 kcal·mol−1 (mean unsigned error). A combination of strain energy released across the reaction coordinate and the critical intramolecular distance between reacting diynes explains reactivity differences. Where experimental and calculated barrier heights are in disagreement, higher level multireference treatment of the enediynes confirms lower level estimates. Previous work concerning the chemically reactive fragment of esperamcin, MTC, is expanded to our model system MTC2.

A Metric to Evaluate and Synthesize Distributed Compliant Mechanisms

Compliant mechanisms with evenly distributed stresses have better load-bearing ability and larger range of motion than mechanisms with compliance and stresses lumped at flexural hinges. In this paper, we present a metric to quantify how uniformly the strain energy of deformation and thus the stresses are distributed throughout the mechanism topology. The resulting metric is used to optimize cross-sections of conceptual compliant topologies leading to designs with maximal stress distribution. This optimization framework is demonstrated for both single-port mechanisms and single-input single-output mechanisms. It is observed that the optimized designs have lower stresses than their nonoptimized counterparts, which implies an ability for single-port mechanisms to store larger strain energy, and single-input single-output mechanisms to perform larger output work before failure.

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.

Subgrain boundaries and related microstructural features in EDML (Antarctica) deep ice core

Subgrain boundaries revealed as shallow sublimation grooves on ice sample surfaces are a direct and easily observable feature of intracrystalline deformation and recrystallization. Statistical data obtained from the EPICA Dronning Maud Land (EDML) deep ice core drilled in East Antarctica cannot detect a depth region of increased subgrain-boundary formation. Grain-boundary morphologies show a strong influence of internal strain energy on the microstructure at all depths. The data do not support the classical view of a change of dominating recrystallization regimes with depth. Three major types of subgrain boundaries, reflecting high mechanical anisotropy, are specified in combination with crystal-orientation analysis.