Boundary condition identification of a bridge pier using experimental data and FEM modal updating

Over the past twenty years, new technologies have required an increasing use of mathematical models in order to understand better the structural behavior: finite element method is the one mostly used. However, the reliability of this method applied to different situations has to be tried each time. Since it is not possible to completely model the reality, different hypothesis must be done: these are the main problems of FE modeling. The following work deals with this problem and tries to figure out a way to identify some of the unknown main parameters of a structure. This main research focuses on a particular path of study and development, but the same concepts can be applied to other objects of research. The main purpose of this work is the identification of unknown boundary conditions of a bridge pier using the data acquired experimentally with field tests and a FEM modal updating process. This work doesn’t want to be new, neither innovative. A lot of work has been done during the past years on this main problem and many solutions have been shown and published. This thesis just want to rework some of the main aspects of the structural optimization process, using a real structure as fitting model.

Finite element analysis of aortic root dilation: a new procedure to reproduce pathology based on experimental data

Sinotubular junction dilation is one of the most frequent pathologies associated with aortic root incompetence. Hence, we create a finite element model considering the whole root geometry; then, starting from healthy valve models and referring to measures of pathological valves reported in the literature, we reproduce the pathology of the aortic root by imposing appropriate boundary conditions. After evaluating the virtual pathological process, we are able to correlate dimensions of non-functional valves with dimensions of competent valves. Such a relation could be helpful in recreating a competent aortic root and, in particular, it could provide useful information in advance in aortic valve sparing surgery.

Finite element study on fracture patterns in human scaphoid wrist bone due to free fall

Scaphoid is one of the 8 carpal bones found adjacent to the thumb supported proximally by Radius bone. During the free fall, on outstretched hand, the impact load gets transferred to the scaphoid at its free anterior end. Unique arrangement of other carpal bones in the palm is also one of the reasons for the load to get transferred to scaphoid. About half of the total load acting upon carpal bone gets transferred to scaphoid at its distal pole. There are about 10 to 12 clinically observed fracture pattern in the scaphoid due to free fall. The aim of the study is to determine the orientation of the load, magnitude of the load and the corresponding fracture pattern. This study includes both static and dynamic finite element models validated by experiments. The scaphoid model has been prepared from CT scans of a 27 year old person. The 2D slices of the CT scans have been converted to 3D model by using MIMICS software. There are four cases of loading studied which are considered to occur clinically more frequently. In case (i) the load is applied at the posterior end at distal pole whereas in case (ii), (iii) and (iv), the load is applied at anterior end at different directions. The model is given a fixed boundary condition at the region which is supported by Radius bone during the impact. Same loading and boundary conditions have been used in both static and dynamic explicit finite element analysis. The site of fracture initiation and path of fracture propagation have been identified by using max principal stress / gradient and max principal strain / gradient criterion respectively in static and dynamic explicit finite element analysis. Static and dynamic impact experiments were performed on the polyurethane foam specimens to validate the finite element results. Experimental results such as load at fracture, site of fracture initiation and path of fracture propagation have been compared with the results of finite element analysis. Four different types of fracture patterns observed in clinical studies have been identified in this study.

Finite element analysis for prediction of bone strength

Article preview View full access options BoneKEy Reports | Review Print Email Share/bookmark Finite element analysis for prediction of bone strength Philippe K Zysset, Enrico Dall'Ara, Peter Varga & Dieter H Pahr Affiliations Corresponding author BoneKEy Reports (2013) 2, Article number: 386 (2013) doi:10.1038/bonekey.2013.120 Received 03 January 2013 Accepted 25 June 2013 Published online 07 August 2013 Article tools Citation Reprints Rights & permissions Abstract Abstract• References• Author information Finite element (FE) analysis has been applied for the past 40 years to simulate the mechanical behavior of bone. Although several validation studies have been performed on specific anatomical sites and load cases, this study aims to review the predictability of human bone strength at the three major osteoporotic fracture sites quantified in recently completed in vitro studies at our former institute. Specifically, the performance of FE analysis based on clinical computer tomography (QCT) is compared with the ones of the current densitometric standards, bone mineral content, bone mineral density (BMD) and areal BMD (aBMD). Clinical fractures were produced in monotonic axial compression of the distal radii, vertebral sections and in side loading of the proximal femora. QCT-based FE models of the three bones were developed to simulate as closely as possible the boundary conditions of each experiment. For all sites, the FE methodology exhibited the lowest errors and the highest correlations in predicting the experimental bone strength. Likely due to the improved CT image resolution, the quality of the FE prediction in the peripheral skeleton using high-resolution peripheral CT was superior to that in the axial skeleton with whole-body QCT. Because of its projective and scalar nature, the performance of aBMD in predicting bone strength depended on loading mode and was significantly inferior to FE in axial compression of radial or vertebral sections but not significantly inferior to FE in side loading of the femur. Considering the cumulated evidence from the published validation studies, it is concluded that FE models provide the most reliable surrogates of bone strength at any of the three fracture sites.

Diagenesis of a light, tight-oil chert reservoir at the Ediacaran/Cambrian boundary, Sultanate of Oman

The Al Shomou Silicilyte Member (Athel Formation) in the South Oman Salt Basin shares many of the characteristics of a light, tight-oil (LTO) reservoir: it is a prolifi c source rock mature for light oil, it produces light oil from a very tight matrix and reservoir, and hydraulic fracking technology is required to produce the oil. What is intriguing about the Al Shomou Silicilyte, and different from other LTO reservoirs, is its position related to the Precambrian/Cambrian Boundary (PCB) and the fact that it is a ‘laminated chert‘ rather than a shale. In an integrated diagenetic study we applied microstructural analyses (SEM, BSE) combined with state-of-the-art stable isotope and trace element analysis of the silicilyte matrix and fractures. Fluid inclusion microthermometry was applied to record the salinity and minimum trapping temperatures. The microstructural investigations reveal a fi ne lamination of the silicilyte matrix with a mean lamina thickness of ca. 20 μm consisting of predominantly organic matter-rich and fi nely crystalline quartz-rich layers, respectively. Authigenic, micron-sized idiomorphic quartz crystals are the main matrix components of the silicilyte. Other diagenetic phases are pyrite, apatite, dolomite, magnesite and barite cements. Porosity values based on neutron density logs and core plug data indicate porosity in the silicilyte ranges from less than 2% to almost to 40%. The majority of the pore space in the silicilyte is related to (primary) inter-crystalline pores, with locally important oversized secondary pores. Pore casts of the silica matrix show that pores are extremely irregular in three dimensions, and are generally interconnected by a complex web or meshwork of fi ne elongate pore throats. Mercury injection capillary data are in line with the microstructural observations suggesting two populations of pore throats, with an effective average modal diameter of 0.4 μm. The acquired geochemical data support the interpretation that the primary source of the silica is the ambient seawater rather than hydrothermal or biogenic. A maximum temperature of ca. 45°C for the formation of microcrystalline quartz in the silicilyte is good evidence that the lithifi cation and crystallization of quartz occurred in the fi rst 5 Ma after deposition. Several phases of brittle fracturing and mineralization occurred in response to salt tectonics during burial. The sequences of fracture-fi lling mineral phases (dolomite - layered chalcedony – quartz – apatite - magnesite I+II - barite – halite) indicates a complex fl uid evolution after silicilyte lithifi cation. Primary, all-liquid fl uid inclusions in the fracturefi lling quartz are good evidence of growth beginning at low temperatures, i.e. ≤ 50ºC. Continuous precipitation during increasing temperature and burial is documented by primary two-phase fl uid inclusions in quartz cements that show brines at 50°C and fi rst hydrocarbons at ca. 70°C. The absolute timing of each mineral phase can be constrained based on U-Pb geochronometry, and basin modelling. Secondary fl uid inclusions in quartz, magnesite and barite indicate reactivation of the fracture system after peak burial temperature during the major cooling event, i.e. uplift, between 450 and 310 Ma. A number of fi rst-order trends in porosity and reservoir-quality distribution are observed which are strongly related to the diagenetic and fl uid history of the reservoir: the early in-situ generation of hydrocarbons and overpressure development arrests diagenesis and preserves matrix porosity. Chemical compaction by pressure dissolution in the fl ank areas could be a valid hypothesis to explain the porosity variations in the silicilitye slabs resulting in lower porosity and poorer connectivity on the fl anks of the reservoir. Most of the hydrocarbon storage and production comes from intervals characterized by Amthor et al. 114488 preserved micropores, not hydrocarbon storage in a fracture system. The absence of oil expulsion results in present-day high oil saturations. The main diagenetic modifi cations of the silicilyte occurred and were completed relatively early in its history, i.e. before 300 Ma. An instrumental factor for preserving matrix porosity is the diffi culty for a given slab to evacuate all the fl uids (water and hydrocarbons), or in other words, the very good sealing capacity of the salt embedding the slab.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p ≤ 0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Effect of boundary conditions on yield properties of human femoral trabecular bone

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7–72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient R2≥0.93. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.

(Table 1) Element analyses of DSDP Hole 15-149

An iridium anomaly has been found in coincidence with the known microtektite level in cores from Deep Sea Drilling Project site 149 in the Caribbean Sea. The iridium was probably not in the microtektites but deposited simultaneously with them; this could occur if the iridium was deposited from a dust cloud resulting from a bolide impact, as suggested for the anomaly associated with the Cretaceous-Tertiary boundary. Other workers have deduced that the microtektites are part of the North American strewn tektite field, which is dated at about 34 million years before present, and that the microtektite horizon in deep-sea cores is synchronous with the extinction of five radiolarian species. Mass extinctions also occur in terrestrial mammals within 4 million years of this time. The iridium anomaly and the tektites and microtektites are supportive of a major bolide impact about 34 million years ago.

Major and rare earth element composition of ODP Site 191-1179 sediments, northwest Pacific

Depth profiles of major element and rare earth element (REE) abundances in sediment samples (mainly siliceous ooze and clay) recovered from Holes 1179B and 1179C at Site 1179, Ocean Drilling Program Leg 191 (41.4°N, 159.6°E) were determined. The oxidation states of Mn and Ce were determined by X-ray absorption near-edge structure. Some geochemical indicators were tested, including the MnO/TiO2 ratios, a bivariate diagram of La/Ce vs. Al2O3/(Al2O3+Fe2O3), and other discrimination diagrams. The oxidation state of Mn is reduced Mn(II) in the depth profile below 0.60 meters below seafloor (mbsf), which is consistent with relatively low and high abundances of Mn in the sediments and pore waters, respectively. It is possible that the diagenetic effect on the oxidation state and abundance of Mn makes it difficult for the MnO/TiO2 ratio to reflect the depositional environment. The normalized ratio of La and Ce does not change very much with depth, suggesting that the diagenetic effect does not affect the REE signature in the sediments. On the diagram of La/Ce vs. Al2O3/(Al2O3+Fe2O3), the sediments studied here plot at the boundary of the pelagic and continental margin fields. This suggests that continental material has contributed to the sediment to some degree, even though Site 1179 is in a pelagic region of the northwestern Pacific Ocean, >1600 km from Japan.

Mineralogy and major element oxide composition of sediments from ODP Hole 172-1063D

Dansgaard-Oeschger (D-O) cycles in sediment at Site 1063 are characterized by distinct fluctuations in physical properties. Stadials are marked by low bulk density and interstadials by high bulk density. Compressional (P-)wave velocity is in phase with bulk density over some but not all depth intervals. Four of the D-O cycles straddling the oxygen isotope Stage 4/5 boundary have been studied in detail to understand the origin of the physical properties changes. Sediment on the Bermuda Rise is comprised of three main components: calcite, aluminosilicate minerals, and biogenic silica. Calcite concentrations vary from 1% to 43% of bulk sediment and are highest during interstadials. Aluminosilicate concentrations vary from 52% to 92% of bulk sediment and are highest during stadials. The major element ratios Al2O3/TiO2 and K2O/Al2O3 show increases across bulk density cycles, suggesting a change in the composition of aluminosilicates. This interpretation is supported by mineralogical analyses, which show a subtle change in clay composition. Biogenic silica concentrations vary from 0% to 23% of bulk sediment and are also highest during stadials. However, the abundance of silica varies significantly from one D-O cycle to another. Silt and fine sand abundance also increase during the first of the four stadials. This coarsening of sediment coincides with the increase in biogenic silica. The low grain density and high porosity associated with biogenic silica result in intervals of low bulk-sediment density. The abundance of biogenic silica closely matches P-wave velocity, suggesting that silica imparts a greater rigidity to the sediment.

Mineralogy and stable oxygen isotope ratios of the Cretaceous-Tertiary boundary

Results of detailed mineralogical, chemical, and oxygen isotope analyses of the clay minerals and zeolites from two Cretaceous-Tertiary (K/T) boundary regions, Stevns Klint, Denmark, and Deep Sea Drilling Project (DSDP) Hole 465A in the north central Pacific Ocean, are presented. In the central part of the Stevns Klint K/T boundary layer, the only clay mineral detected by x-ray diffraction is a pure smectite with > 95 percent expandable layers. No detrital clay minerals or quartz were observed in the clay size fraction in these beds, whereas the clay minerals above and below the boundary layer are illite and mixed-layer smectite-illite of detrital origin as well as quartz. The mineralogical purity of the clay fraction, the presence of smectite only at the boundary, and the d18O value of the smectite (27.2 ± 0.2 per mil) suggest that it formed in situ by alteration of glass. Formation from impact rather than from volcanic glass is supported by its major element chemistry. The high content of iridium and other siderophile elements is not due to the cessation of calcium carbonate deposition and resulting slow sedimentation rates. At DSDP Hole 465A, the principal clay mineral in the boundary zone (80 to 143 centimeters) is a mixed-layer smectite-illite with >=90 percent expandable layers, accompanied by some detrital quartz and small amounts of a euhedral authigenic zeolite (clinoptilolite). The mixed-layer smectite-illite from the interval 118 to 120 centimeters in the zone of high iridium abundance has a very low rare earth element content; the negative cerium anomaly indicates formation in the marine environment. This conclusion is corroborated by the d18O value of this clay mineral (27.1 ± 0.2 per mil). Thus, this mixed-layer smectite-illite formed possibly from the same glass as the K/T boundary smectite at Stevns Klint, Denmark.

Stable isotope and trace element geochemistry of carbonate sediments at DSDP Holes 87-577 and 6-47.2

Detailed analyses of well-preserved carbonate samples from across the Cretaceous/Tertiary boundary in Hole 577 have revealed a significant decline in the d13C values of calcareous nannoplankton from the Maestrichtian to the Danian Age accompanied by a substantial reduction in carbonate accumulation rates. Benthic foraminifers, however, do not exhibit a shift in carbon composition similar to that recorded by the calcareous nannoplankton, but actually increase slightly over the same time interval. These results are similar to the earlier findings at two North Pacific Deep Sea Drilling Project locations, Sites 47.2 and 465, and are considered to represent a dramatic decrease in oceanic phytoplankton production associated with the catastrophic Cretaceous/Tertiary boundary extinctions. In addition, the change in carbon composition of calcareous nannoplankton across the Cretaceous/Tertiary boundary at Hole 577 is accompanied by only minor changes in the oxygen isotope trends of both calcareous nannoplankton and benthic foraminifers, suggesting that temperature variations in the North Pacific from the late Maestrichtian to the early Danian Age were insignificant.

Major element concentrations in spinel, olivine and pyroxene from peridotites of ODP Hole 210-1277 of the Newfoundland margin

Serpentinized spinel peridotites of the Newfoundland margin drilled during ODP Leg 210 at Site 1277 have preserved, relic mineral compositions similar to the most depleted abyssal peridotites worldwide and different from those of the conjugate Iberian margin. The samples are derived from mass flows containing clasts of peridotite and gabbro and from in-situ basement, and are mostly mylonitic cpx-poor spinel harzburgites with Cr-rich spinels (Cr#0.35-0.66). Melting of the Newfoundland mantle occurred in the spinel peridotite field and probably exceeded the cpx-out phase boundary for some samples. Using proposed spinel peridotite melting models and experimentally derived phase diagrams, the Newfoundland harzburgites can be modeled as a residue after extraction of 14 to 20-25% melting. Basalts that are interleaved with mass flow deposits on top of the peridotite basement resemble normal to transitional mid-ocean ridge basalt. This, together with the unusually high Cr# of some spinel harzburgites suggest that the formation of basalts and partial melting of the underlying peridotite are not cogenetic. Among mantle samples some of the Newfoundland harzburgites approach mineral compositions of the Bay of island ophiolite and ophiolites from Japan that represent peridotites formed in an arc-setting. Thus, the peridotites drilled at Site 1277 may represent inherited (Caledonian or older) subarc mantle that was exhumed close to the ocean floor during the rifting evolution of the Atlantic. Compared to the spinel harzburgites from Newfoundland, the peridotites from the conjugate Iberian margin are, on average, less depleted and provide evidence for local equilibration in the plagioclase stability field. This can either be explained by an inherited, primary, Ca-richer composition of the Iberia peridotite, or, alternatively, by local melt impregnation and stagnation during continental rifting, and thus refertilizing previously depleted (arc-related) peridotite.

Stable organic carbon isotope record and total organic carbon content from the Cenomanian-Turonian boundary

Ocean Drilling Program (ODP) Leg 207 recovered expanded sections of organic-carbon-rich laminated shales on Demerara Rise (western tropical Atlantic). High-resolution organic carbon isotope and total organic carbon (TOC) records are presented, which span the Cenomanian-Turonian boundary interval (CTBI), including the Oceanic Anoxic Event (OAE) 2, from four sites oriented along a NW striking depth transect. These records represent the first high-resolution carbon isotope records across OAE 2 from the South American margin of the tropical Atlantic. Due to the scarcity of age significant fossils, the main purpose of this study was to develop a detailed carbon isotope stratigraphy in order to correlate the CTBI across the depth transect and to tie this to biostratigraphically well-defined sections in the Western Interior Basin (Pueblo, USA), boreal shelf seas (Eastbourne, England), and western Tethys (Oued Mellegue, Tunisia). All four sections studied document a 6 per mil increase of d13Corg values at the base of the CTBI, which is followed by an interval of elevated d13Corg values and a subsequent decrease. Our results supply an important stratigraphic base for subsequent paleoceanographic studies on Late Cenomanian to Early Turonian sediments from Demerara Rise and elsewhere.