Numerical modelling of braided fibres for reinforced concrete

Fire has been always a major concern for designers of steel and concrete structures. Designing fire-resistant structural elements is not an easy task due to several limitations such as the lack of fire-resistant construction materials. Concrete reinforcement cover and external insulation are the most commonly adopted systems to protect concrete and steel from overheating, while spalling of concrete is minimised by using HPFRC instead of standard concrete. Although these methodologies work very well for low rise concrete structures, this is not the case for high-rise and inaccessible buildings where fire loading is much longer. Fire can permanently damage structures that cost a lot of money. This is unsafe and can lead to loss of life. In this research, the author proposes a new type of main reinforcement for concrete structures which can provide better fire-resistance than steel or FRP re-bars. This consists of continuous braided fibre rope, generally made from fire-resistant materials such as carbon or glass fibre. These fibres have excellent tensile strengths, sometimes in excess of ten times greater than steel. In addition to fire-resistance, these ropes can produce lighter and corrosive resistant structures. Avoiding the use of expensive resin binders, fibres are easily bound together using braiding techniques, ensuring that tensile stress is evenly distributed throughout the reinforcement. In order to consider braided ropes as a form of reinforcement it is first necessary to establish the mechanical performance at room temperature and investigate the pull-out resistance for both unribbed and ribbed ropes. Ribbing of ropes was achieved by braiding the rope over a series of glass beads. Adhesion between the rope and concrete was drastically improved due to ribbing, and further improved by pre-stressing ropes and reducing the slacked fibres. Two types of material have been considered for the ropes: carbon and aramid. An implicit finite element approach is proposed to model braided fibres using Total Lagrangian formulation, based on the theory of small strains and large rotations. Modelling tows and strands as elastic transversely isotropic materials was a good assumption when stiff and brittle fibres such as carbon and glass fibres are considered. The rope-to-concrete and strand-to-strand bond interaction/adhesion was numerically simulated using newly proposed hierarchical higher order interface elements. Elastic and linear damage cohesive models were used effectively to simulate non-penetrative 'free' sliding interaction between strands, and the adhesion between ropes and concrete respectively. Numerical simulation showed similar de-bonding features when compared with experimental pull-out results of braided ribbed rope reinforced concrete.

Modelling and studies of the spectral response of some optoelectronic components

Solar radiation takes in today's world, an increasing importance. Different devices are used to carry out spectral and integrated measurements of solar radiation. Thus the sensors can be divided into the fallow types: Calorimetric, Thermomechanical, Thermoelectric and Photoelectric. The first three categories are based on components converting the radiation to temperature (or heat) and then into electrical quantity. On the other hand, the photoelectric sensors are based on semiconductor or optoelectronic elements that when irradiated change their impedance or generate a measurable electric signal. The response function of the sensor element depends not only on the intensity of the radiation but also on its wavelengths. The radiation sensors most widely used fit in the first categories, but thanks to the reduction in manufacturing costs and to the increased integration of electronic systems, the use of the photoelectric-type sensors became more interesting. In this work we present a study of the behavior of different optoelectronic sensor elements. It is intended to verify the static response of the elements to the incident radiation. We study the optoelectronic elements using mathematical models that best fit their response as a function of wavelength. As an input to the model, the solar radiation values are generated with a radiative transfer model. We present a modeling of the spectral response sensors of other types in order to compare the behavior of optoelectronic elements with other sensors currently in use.

Numerical modelling of rock masses interacting with dams: the case of Ridracoli

A three-dimensional Direct Finite Element procedure is here presented which takes into account most of the factors affecting the interaction problem of the dam-water-foundation system, whilst keeping the computational cost at a reasonable level by introducing some simplified hypotheses. A truncated domain is defined, and the dynamic behaviour of the system is treated as a wave-scattering problem where the presence of the dam perturbs an original free-field system. The rock foundation truncated boundaries are enclosed by a set of free-field one-dimensional and two-dimensional systems which transmit the effective forces to the main model and apply adsorbing viscous boundaries to ensure radiation damping. The water domain is treated as an added mass moving with the dam. A strategy is proposed to keep the viscous dampers at the boundaries unloaded during the initial phases of analysis, when the static loads are initialised, and thus avoid spurious displacements. A focus is given to the nonlinear behaviour of the rock foundation, with concentrated plasticity along the natural discontinuities of the rock mass, immersed in an otherwise linear elastic medium with Rayleigh damping. The entire procedure is implemented in the commercial software Abaqus®, whose base code is enriched with specific user subroutines when needed. All the extra coding is attached to the Thesis and tested against analytical results and simple examples. Possible rock wedge instabilities induced by intense ground motion, which are not easily investigated within a comprehensive model of the dam-water-foundation system, are treated separately with a simplified decoupled dynamic approach derived from the classical Newmark method, integrated with FE calculation of dam thrust on the wedges during the earthquake. Both the described approaches are applied to the case study of the Ridracoli arch-gravity dam (Italy) in order to investigate its seismic response to the Maximum Credible Earthquake (MCE) in a full reservoir condition.

Optimizing fused deposition modelling: the impact of geometrical slicing parameters on mechanical properties

This comprehensive study explores the intricate world of 3D printing, with a focus on Fused Deposition Modelling (FDM). It sheds light on the critical factors that influence the quality and mechanical properties of 3D printed objects. Using an optical microscope with 40X magnification, the shapes of the printed beads is correlated to specific slicing parameters, resulting in a 2D parametric model. This mathematical model, derived from real samples, serves as a tool to predict general mechanical behaviour, bridging the gap between theory and practice in FDM printing. The study begins by emphasising the importance of geometric parameters such as layer height, line width and filament tolerance on the final printed bead geometry and the resulting theoretical effect on mechanical properties. The introduction of VPratio parameter (ratio between the area of the voids and the area occupied by printed material) allows the quantification of the variation of geometric slicing parameters on the improvement or reduction of mechanical properties. The study also addresses the effect of overhang and the role of filament diameter tolerances. The research continues with the introduction of 3D FEM (Finite Element Analysis) models based on the RVE (Representative Volume Element) to verify the results obtained from the 2D model and to analyse other aspects that affect mechanical properties and not directly observable with the 2D model. The study also proposes a model for the examination of 3D printed infill structures, introducing also an innovative methodology called “double RVE” which speeds up the calculation of mechanical properties and is also more computationally efficient. Finally, the limitations of the RVE model are shown and a so-called Hybrid RVE-based model is created to overcome the limitations and inaccuracy of the conventional RVE model and homogenization procedure on some printed geometries.

Effect Of Platform Connection And Abutment Material On Stress Distribution In Single Anterior Implant-supported Restorations: A Nonlinear 3-dimensional Finite Element Analysis.

Although various abutment connections and materials have recently been introduced, insufficient data exist regarding the effect of stress distribution on their mechanical performance. The purpose of this study was to investigate the effect of different abutment materials and platform connections on stress distribution in single anterior implant-supported restorations with the finite element method. Nine experimental groups were modeled from the combination of 3 platform connections (external hexagon, internal hexagon, and Morse tapered) and 3 abutment materials (titanium, zirconia, and hybrid) as follows: external hexagon-titanium, external hexagon-zirconia, external hexagon-hybrid, internal hexagon-titanium, internal hexagon-zirconia, internal hexagon-hybrid, Morse tapered-titanium, Morse tapered-zirconia, and Morse tapered-hybrid. Finite element models consisted of a 4×13-mm implant, anatomic abutment, and lithium disilicate central incisor crown cemented over the abutment. The 49 N occlusal loading was applied in 6 steps to simulate the incisal guidance. Equivalent von Mises stress (σvM) was used for both the qualitative and quantitative evaluation of the implant and abutment in all the groups and the maximum (σmax) and minimum (σmin) principal stresses for the numerical comparison of the zirconia parts. The highest abutment σvM occurred in the Morse-tapered groups and the lowest in the external hexagon-hybrid, internal hexagon-titanium, and internal hexagon-hybrid groups. The σmax and σmin values were lower in the hybrid groups than in the zirconia groups. The stress distribution concentrated in the abutment-implant interface in all the groups, regardless of the platform connection or abutment material. The platform connection influenced the stress on abutments more than the abutment material. The stress values for implants were similar among different platform connections, but greater stress concentrations were observed in internal connections.

The Basic Reproduction Number Obtained From Jacobian And Next Generation Matrices - A Case Study Of Dengue Transmission Modelling.

The basic reproduction number is a key parameter in mathematical modelling of transmissible diseases. From the stability analysis of the disease free equilibrium, by applying Routh-Hurwitz criteria, a threshold is obtained, which is called the basic reproduction number. However, the application of spectral radius theory on the next generation matrix provides a different expression for the basic reproduction number, that is, the square root of the previously found formula. If the spectral radius of the next generation matrix is defined as the geometric mean of partial reproduction numbers, however the product of these partial numbers is the basic reproduction number, then both methods provide the same expression. In order to show this statement, dengue transmission modelling incorporating or not the transovarian transmission is considered as a case study. Also tuberculosis transmission and sexually transmitted infection modellings are taken as further examples.

Photoelastic And Finite Element Analyses Of Occlusal Loads In Mandibular Body.

This study proposed to evaluate the mandibular biomechanics in the posterior dentition based on experimental and computational analyses. The analyses were performed on a model of human mandible, which was modeled by epoxy resin for photoelastic analysis and by computer-aided design for finite element analysis. To standardize the evaluation, specific areas were determined at the lateral surface of mandibular body. The photoelastic analysis was configured through a vertical load on the first upper molar and fixed support at the ramus of mandible. The same configuration was used in the computer simulation. Force magnitudes of 50, 100, 150, and 200 N were applied to evaluate the bone stress. The stress results presented similar distribution in both analyses, with the more intense stress being at retromolar area and oblique line and alveolar process at molar level. This study presented the similarity of results in the experimental and computational analyses and, thus, showed the high importance of morphology biomechanical characterization at posterior dentition.

Finite-element Analysis Of 3 Situations Of Trauma In The Human Edentulous Mandible.

Maxillofacial trauma resulting from falls in elderly patients is a major social and health care concern. Most of these traumatic events involve mandibular fractures. The aim of this study was to analyze stress distributions from traumatic loads applied on the symphyseal, parasymphyseal, and mandibular body regions in the elderly edentulous mandible using finite-element analysis (FEA). Computerized tomographic analysis of an edentulous macerated human mandible of a patient approximately 65 years old was performed. The bone structure was converted into a 3-dimensional stereolithographic model, which was used to construct the computer-aided design (CAD) geometry for FEA. The mechanical properties of cortical and cancellous bone were characterized as isotropic and elastic structures, respectively, in the CAD model. The condyles were constrained to prevent free movement in the x-, y-, and z-axes during simulation. This enabled the simulation to include the presence of masticatory muscles during trauma. Three different simulations were performed. Loads of 700 N were applied perpendicular to the surface of the cortical bone in the symphyseal, parasymphyseal, and mandibular body regions. The simulation results were evaluated according to equivalent von Mises stress distributions. Traumatic load at the symphyseal region generated low stress levels in the mental region and high stress levels in the mandibular neck. Traumatic load at the parasymphyseal region concentrated the resulting stress close to the mental foramen. Traumatic load in the mandibular body generated extensive stress in the mandibular body, angle, and ramus. FEA enabled precise mapping of the stress distribution in a human elderly edentulous mandible (neck and mandibular angle) in response to 3 different traumatic load conditions. This knowledge can help guide emergency responders as they evaluate patients after a traumatic event.

Keratoconus prediction using a finite element model of the cornea with local biomechanical properties

PURPOSE: The ability to predict and understand which biomechanical properties of the cornea are responsible for the stability or progression of keratoconus may be an important clinical and surgical tool for the eye-care professional. We have developed a finite element model of the cornea, that tries to predicts keratoconus-like behavior and its evolution based on material properties of the corneal tissue. METHODS: Corneal material properties were modeled using bibliographic data and corneal topography was based on literature values from a schematic eye model. Commercial software was used to simulate mechanical and surface properties when the cornea was subject to different local parameters, such as elasticity. RESULTS: The simulation has shown that, depending on the corneal initial surface shape, changes in local material properties and also different intraocular pressures values induce a localized protuberance and increase in curvature when compared to the remaining portion of the cornea. CONCLUSIONS: This technique provides a quantitative and accurate approach to the problem of understanding the biomechanical nature of keratoconus. The implemented model has shown that changes in local material properties of the cornea and intraocular pressure are intrinsically related to keratoconus pathology and its shape/curvature.

Stress distribution in the cervical region of an upper central incisor in a 3D finite element model

The aim of this study was to evaluate the stress distribution in the cervical region of a sound upper central incisor in two clinical situations, standard and maximum masticatory forces, by means of a 3D model with the highest possible level of fidelity to the anatomic dimensions. Two models with 331,887 linear tetrahedral elements that represent a sound upper central incisor with periodontal ligament, cortical and trabecular bones were loaded at 45º in relation to the tooth's long axis. All structures were considered to be homogeneous and isotropic, with the exception of the enamel (anisotropic). A standard masticatory force (100 N) was simulated on one of the models, while on the other one a maximum masticatory force was simulated (235.9 N). The software used were: PATRAN for pre- and post-processing and Nastran for processing. In the cementoenamel junction area, tensile forces reached 14.7 MPa in the 100 N model, and 40.2 MPa in the 235.9 N model, exceeding the enamel's tensile strength (16.7 MPa). The fact that the stress concentration in the amelodentinal junction exceeded the enamel's tensile strength under simulated conditions of maximum masticatory force suggests the possibility of the occurrence of non-carious cervical lesions such as abfractions.

Rare earth element and yttrium geochemistry applied to the genetic study of cryolite ore at the Pitinga Mine (Amazon, Brazil)

This work aims at the geochemical study of Pitinga cryolite mineralization through REE and Y analyses in disseminated and massive cryolite ore deposits, as well as in fluorite occurrences. REE signatures in fluorite and cryolite are similar to those in the Madeira albite granite. The highest ΣREE values are found in magmatic cryolite (677 to 1345 ppm); ΣREE is lower in massive cryolite. Average values for the different cryolite types are 10.3 ppm, 6.66 ppm and 8.38 ppm (for nucleated, caramel and white types, respectively). Disseminated fluorite displays higher ΣREE values (1708 and 1526ppm) than fluorite in late veins(34.81ppm). Yttrium concentration is higher in disseminated fluorite and in magmatic cryolite. The evolution of several parameters (REEtotal, LREE/HREE, Y) was followed throughout successive stages of evolution in albite granites and associated mineralization. At the end of the process, late cryolite was formed with low REEtotal content. REE data indicate that the MCD was formed by, and the disseminated ore enriched by (additional formation of hydrothermal disseminated cryolite), hydrothermal fluids, residual from albite granite. The presence of tetrads is poorly defined, although nucleated, caramel and white cryolite types show evidence for tetrad effect.

Unconstrained Finite Element for Geometrical Nonlinear Dynamics of Shells

This paper presents a positional FEM formulation to deal with geometrical nonlinear dynamics of shells. The main objective is to develop a new FEM methodology based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors not displacements and rotations. These characteristics are the novelty of the present work and avoid the use of large rotation approximations. A nondimensional auxiliary coordinate system is created, and the change of configuration function is written following two independent mappings from which the strain energy function is derived. This methodology is called positional and, as far as the authors' knowledge goes, is a new procedure to approximated geometrical nonlinear structures. In this paper a proof for the linear and angular momentum conservation property of the Newmark beta algorithm is provided for total Lagrangian description. The proposed shell element is locking free for elastic stress-strain relations due to the presence of linear strain variation along the shell thickness. The curved, high-order element together with an implicit procedure to solve nonlinear equations guarantees precision in calculations. The momentum conserving, the locking free behavior, and the frame invariance of the adopted mapping are numerically confirmed by examples. Copyright (C) 2009 H. B. Coda and R. R. Paccola.

Three-dimensional finite element thermal analysis of dental tissues irradiated with Er,Cr:YSGG laser

In the present study, a finite element model of a half-sectioned molar tooth was developed in order to understand the thermal behavior of dental hard tissues (both enamel and dentin) under laser irradiation. The model was validated by comparing it with an in vitro experiment where a sound molar tooth was irradiated by an Er,Cr:YSGG pulsed laser. The numerical tooth model was conceived to simulate the in vitro experiment, reproducing the dimensions and physical conditions of the typical molar sound tooth, considering laser energy absorption and calculating the heat transfer through the dental tissues in three dimensions. The numerical assay considered the same three laser energy densities at the same wavelength (2.79 mu m) used in the experiment. A thermographic camera was used to perform the in vitro experiment, in which an Er, Cr: YSGG laser (2.79 mu m) was used to irradiate tooth samples and the infrared images obtained were stored and analyzed. The temperature increments in both the finite element model and the in vitro experiment were compared. The distribution of temperature inside the tooth versus time plotted for two critical points showed a relatively good agreement between the results of the experiment and model. The three dimensional model allows one to understand how the heat propagates through the dentin and enamel and to relate the amount of energy applied, width of the laser pulses, and temperature inside the tooth. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.2953526]

Alpha element abundances and gradients in the Milky Way bulge from FLAMES-GIRAFFE spectra of 650 K giants

Aims. We present the analysis of the [alpha/Fe] abundance ratios for a large number of stars at several locations in the Milky Way bulge with the aim of constraining its formation scenario. Methods. We obtained FLAMES-GIRAFFE spectra (R = 22 500) at the ESO Very Large Telescope for 650 bulge red giant branch (RGB) stars and performed spectral synthesis to measure Mg, Ca, Ti, and Si abundances. This sample is composed of 474 giant stars observed in 3 fields along the minor axis of the Galactic bulge and at latitudes b = -4 degrees, b = -6 degrees, b = -12 degrees. Another 176 stars belong to a field containing the globular cluster NGC 6553, located at b = -3 degrees and 5 degrees away from the other three fields along the major axis. Stellar parameters and metallicities for these stars were presented in Zoccali et al. (2008, A&A, 486, 177). We have also re-derived stellar parameters and abundances for the sample of thick and thin disk red giants analyzed in Alves-Brito et al. (2010, A&A, 513, A35). Therefore using a homogeneous abundance database for the bulge, thick and thin disk, we have performed a differential analysis minimizing systematic errors, to compare the formation scenarios of these Galactic components. Results. Our results confirm, with large number statistics, the chemical similarity between the Galactic bulge and thick disk, which are both enhanced in alpha elements when compared to the thin disk. In the same context, we analyze [alpha/Fe] vs. [Fe/H] trends across different bulge regions. The most metal rich stars, showing low [alpha/Fe] ratios at b = -4 degrees disappear at higher Galactic latitudes in agreement with the observed metallicity gradient in the bulge. Metal-poor stars ([Fe/H] < -0.2) show a remarkable homogeneity at different bulge locations. Conclusions. We have obtained further constrains for the formation scenario of the Galactic bulge. A metal-poor component chemically indistinguishable from the thick disk hints for a fast and early formation for both the bulge and the thick disk. Such a component shows no variation, neither in abundances nor kinematics, among different bulge regions. A metal-rich component showing low [alpha/Fe] similar to those of the thin disk disappears at larger latitudes. This allows us to trace a component formed through fast early mergers (classical bulge) and a disk/bar component formed on a more extended timescale.

Neutron Activation Analysis: A Primary (Ratio) Method to Determine SI-Traceable Values of Element Content in Complex Samples

The metrological principles of neutron activation analysis are discussed. It has been demonstrated that this method can provide elemental amount of substance with values fully traceable to the SI. The method has been used by several laboratories worldwide in a number of CCQM key comparisons - interlaboratory comparison tests at the highest metrological level - supplying results equivalent to values from other methods for elemental or isotopic analysis in complex samples without the need to perform chemical destruction and dissolution of these samples. The CCOM accepted therefore in April 2007 the claim that neutron activation analysis should have the similar status as the methods originally listed by the CCOM as `primary methods of measurement`. Analytical characteristics and scope of application are given.