Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

Benefit of Cup Medialization in Total Hip Arthroplasty is Associated With Femoral Anatomy.

BACKGROUND: Medialization of the cup with a respective increase in femoral offset has been proposed in THA to increase abductor moment arms. Insofar as there are potential disadvantages to cup medialization, it is important to ascertain whether the purported biomechanical benefits of cup medialization are large enough to warrant the downsides; to date, studies regarding this question have disagreed. QUESTIONS/PURPOSES: The purpose of this study was to quantify the effect of cup medialization with a compensatory increase in femoral offset compared with anatomic reconstruction for patients undergoing THA. We tested the hypothesis that there is a (linear) correlation between preoperative anatomic parameters and muscle moment arm increase caused by cup medialization. METHODS: Fifteen patients undergoing THA were selected, covering a typical range of preoperative femoral offsets. For each patient, a finite element model was built based on a preoperative CT scan. The model included the pelvis, femur, gluteus minimus, medius, and maximus. Two reconstructions were compared: (1) anatomic position of the acetabular center of rotation, and (2) cup medialization compensated by an increase in the femoral offset. Passive abduction-adduction and flexion-extension were simulated in the range of normal gait. Muscle moment arms were evaluated and correlated to preoperative femoral offset, acetabular offset, height of the greater trochanter (relative to femoral center of rotation), and femoral antetorsion angle. RESULTS: The increase of muscle moment arms caused by cup medialization varied among patients. Muscle moment arms increase by 10% to 85% of the amount of cup medialization for abduction-adduction and from -35% (decrease) to 50% for flexion-extension. The change in moment arm was inversely correlated (R(2) = 0.588, p = 0.001) to femoral antetorsion (anteversion), such that patients with less femoral antetorsion gained more in terms of hip muscle moments. No linear correlation was observed between changes in moment arm and other preoperative parameters in this series. CONCLUSIONS: The benefit of cup medialization is variable and depends on the individual anatomy. CLINICAL RELEVANCE: Cup medialization with compensatory increase of the femoral offset may be particularly effective in patients with less femoral antetorsion. However, cup medialization must be balanced against its tradeoffs, including the additional loss of medial acetabular bone stock, and eventual proprioceptive implications of the nonanatomic center of rotation and perhaps joint reaction forces. Clinical studies should better determine the relevance of small changes of moment arms on function and joint reaction forces.

Impedance boundary conditions for pseudo-spectral time-domain methods in room acoustics

The pseudo-spectral time-domain (PSTD) method is an alternative time-marching method to classicalleapfrog finite difference schemes in the simulation of wave-like propagating phenomena. It is basedon the fundamentals of the Fourier transform to compute the spatial derivatives of hyperbolic differential equations. Therefore, it results in an isotropic operator that can be implemented in an efficient way for room acoustics simulations. However, one of the first issues to be solved consists on modeling wallabsorption. Unfortunately, there are no references in the technical literature concerning to that problem. In this paper, assuming real and constant locally reacting impedances, several proposals to overcome this problem are presented, validated and compared to analytical solutions in different scenarios.

Impedance boundary conditions for pseudo-spectral time-domain methods in room acoustics

The Pseudo-Spectral Time Domain (PSTD) method is an alternative time-marching method to classical leapfrog finite difference schemes inthe simulation of wave-like propagating phenomena. It is based on the fundamentals of the Fourier transform to compute the spatial derivativesof hyperbolic differential equations. Therefore, it results in an isotropic operator that can be implemented in an efficient way for room acousticssimulations. However, one of the first issues to be solved consists on modeling wall absorption. Unfortunately, there are no references in thetechnical literature concerning to that problem. In this paper, assuming real and constant locally reacting impedances, several proposals toovercome this problem are presented, validated and compared to analytical solutions in different scenarios.

Solving higher-dimensional continuous time stochastic control problems by value function regression

The paper develops a method to solve higher-dimensional stochasticcontrol problems in continuous time. A finite difference typeapproximation scheme is used on a coarse grid of low discrepancypoints, while the value function at intermediate points is obtainedby regression. The stability properties of the method are discussed,and applications are given to test problems of up to 10 dimensions.Accurate solutions to these problems can be obtained on a personalcomputer.

Effect of tibial tray inclination on the femoro-tibial contact pattern of a mobile total knee arthroplasty

Introduction: The posterior inclination of the tibial component is an important factor that can affect the success of total knee arthroplasty. It can reduce the posterior impingement and thus increase the range of flexion, but it may also induce instability in flexion, anterior impingement between the polyethylene of postero-stabilizing knee prosthesis, and anterior conflict with the cortical bone and the stem. Although the problem is identified, there is still a debate on the ideal inclination angle and the surgical technique to avoid an excessive posterior inclination. The aim of this study was to predict the effect of a posterior inclination of the tibial component on the contact pattern on the tibial insert, using a numerical musculoskeletal model of the knee joint. Methods: A 3D finite element model of the knee joint was developed to simulate an active and loaded squat movement after total knee arthroplasty. Flexion was actively controlled by the quadriceps muscle and muscle activations were estimated from EMG data and were synchronized by a feedback algorithm. Two inclinations of the tibial tray were considered: a posterior inclination of 0° or 10°. During the entire range of flexion, the following quantities were calculated: the tibiofemoral and patello-femoral contact force, and the contact pattern on polyethylene insert. The antero-posterior displacement of the contact pattern was also measured. Abaqus 6.7 was used for all analyses. Results: The tibio-femoral and patello-femoral contact forces increased during flexion and reached respectively 4 and 7 BW (bodyweight) at 90° of flexion. They were slightly affected by the inclination of the tibial tray. Without posterior inclination, the contact pattern on the tibial insert remained centered. The contact pressure was lower than 5 MPa below 60° of flexion, but exceeded 20 MPa at 90° of flexion. The posterior inclination displaced the contact point posteriorly by 2 to 4 mm. Conclusion: The inclination of the tibial tray displaced the contactpattern towards the posterior border of the tibial insert. However, even for 10° of inclination, the contact center remained far from the posterior border (12 mm). There was no instability predicted for this movement.

Early stage disc degeneration does not have an appreciable affect on stiffness and load transfer following vertebroplasty and kyphoplasty.

Vertebroplasty and kyphoplasty have been reported to alter the mechanical behavior of the treated and adjacent-level segments, and have been suggested to increase the risk for adjacent-level fractures. The intervertebral disc (IVD) plays an important role in the mechanical behavior of vertebral motion segments. Comparisons between normal and degenerative IVD motion segments following cement augmentation have yet to be reported. A microstructural finite element model of a degenerative IVD motion segment was constructed from micro-CT images. Microdamage within the vertebral body trabecular structure was used to simulate a slightly (I = 83.5% of intact stiffness), moderately (II = 57.8% of intact stiffness), and severely (III = 16.0% of intact stiffness) damaged motion segment. Six variable geometry single-segment cement repair strategies (models A-F) were studied at each damage level (I-III). IVD and bone stresses, and motion segment stiffness, were compared with the intact and baseline damage models (untreated), as well as, previous findings using normal IVD models with the same repair strategies. Overall, small differences were observed in motion segment stiffness and average stresses between the degenerative and normal disc repair models. We did however observe a reduction in endplate bulge and a redistribution in the microstructural tissue level stresses across both endplates and in the treated segment following early stage IVD degeneration. The cement augmentation strategy placing bone cement along the periphery of the vertebra (model E) proved to be the most advantageous in treating the degenerative IVD models by showing larger reductions in the average bone stresses (vertebral and endplate) as compared to the normal IVD models. Furthermore, only this repair strategy, and the complete cement fill strategy (model F), were able to restore the slightly damaged (I) motion segment stiffness above pre-damaged (intact) levels. Early stage IVD degeneration does not have an appreciable effect in motion segment stiffness and average stresses in the treated and adjacent-level segments following vertebroplasty and kyphoplasty. Placing bone cement in the periphery of the damaged vertebra in a degenerative IVD motion segment, minimizes load transfer, and may reduce the likelihood of adjacent-level fractures.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.

Local character of magnetic coupling in ionic solids

Magnetic interactions in ionic solids are studied using parameter-free methods designed to provide accurate energy differences associated with quantum states defining the Heisenberg constant J. For a series of ionic solids including KNiF3, K2NiF4, KCuF3, K2CuF4, and high- Tc parent compound La2CuO4, the J experimental value is quantitatively reproduced. This result has fundamental implications because J values have been calculated from a finite cluster model whereas experiments refer to infinite solids. The present study permits us to firmly establish that in these wide-gap insulators, J is determined from strongly local electronic interactions involving two magnetic centers only thus providing an ab initio support to commonly used model Hamiltonians.

An Experimental Validation of a Statistical-Based Damage Detection Approach, January 2011

In this work, a previously-developed, statistical-based, damage-detection approach was validated for its ability to autonomously detect damage in bridges. The damage-detection approach uses statistical differences in the actual and predicted behavior of the bridge caused under a subset of ambient trucks. The predicted behavior is derived from a statistics-based model trained with field data from the undamaged bridge (not a finite element model). The differences between actual and predicted responses, called residuals, are then used to construct control charts, which compare undamaged and damaged structure data. Validation of the damage-detection approach was achieved by using sacrificial specimens that were mounted to the bridge and exposed to ambient traffic loads and which simulated actual damage-sensitive locations. Different damage types and levels were introduced to the sacrificial specimens to study the sensitivity and applicability. The damage-detection algorithm was able to identify damage, but it also had a high false-positive rate. An evaluation of the sub-components of the damage-detection methodology and methods was completed for the purpose of improving the approach. Several of the underlying assumptions within the algorithm were being violated, which was the source of the false-positives. Furthermore, the lack of an automatic evaluation process was thought to potentially be an impediment to widespread use. Recommendations for the improvement of the methodology were developed and preliminarily evaluated. These recommendations are believed to improve the efficacy of the damage-detection approach.

Strengthening of Existing Single Span Steel Beam and Concrete Deck Bridges, HR-238, Part II, 1985

The unifying objective of Phases I and II of this study was to determine the feasibility of the post-tensioning strengthening method and to implement the technique on two composite bridges in Iowa. Following completion of these two phases, Phase III was undertaken and is documented in this report. The basic objectives of Phase III were further monitoring bridge behavior (both during and after post-tensioning) and developing a practical design methodology for designing the strengthening system under investigation. Specific objectives were: to develop strain and force transducers to facilitate the collection of field data; to investigate further the existence and effects of the end restraint on the post-tensioning process; to determine the amount of post-tensioning force loss that occurred during the time between the initial testing and the retesting of the existing bridges; to determine the significance of any temporary temperature-induced post-tensioning force change; and to develop a simplified design methodology that would incorporate various variables such as span length, angle-of-skew, beam spacing, and concrete strength. Experimental field results obtained during Phases II and III were compared to the theoretical results and to each other. Conclusions from this research are as follows: (1) Strengthening single-span composite bridges by post-tensioning is a viable, economical strengthening technique. (2) Behavior of both bridges was similar to the behavior observed from the bridges during field tests conducted under Phase II. (3) The strain transducers were very accurate at measuring mid-span strain. (4) The force transducers gave excellent results under laboratory conditions, but were found to be less effective when used in actual bridge tests. (5) Loss of post-tensioning force due to temperature effects in any particular steel beam post-tensioning tendon system were found to be small. (6) Loss of post-tensioning force over a two-year period was minimal. (7) Significant end restraint was measured in both bridges, caused primarily by reinforcing steel being continuous from the deck into the abutments. This end restraint reduced the effectiveness of the post-tensioning but also reduced midspan strains due to truck loadings. (8) The SAP IV finite element model is capable of accurately modeling the behavior of a post-tensioned bridge, if guardrails and end restraints are included in the model. (9) Post-tensioning distribution should be separated into distributions for the axial force and moment components of an eccentric post-tensioning force. (10) Skews of 45 deg or less have a minor influence on post-tensioning distribution. (11) For typical Iowa three-beam and four-beam composite bridges, simple regression-derived formulas for force and moment fractions can be used to estimate post-tensioning distribution at midspan. At other locations, a simple linear interpolation gives approximately correct results. (12) A simple analytical model can accurately estimate the flexural strength of an isolated post-tensioned composite beam.

Design Manual for Strengthening Single-Span Composite Bridges by Post-Tensioning Part 3, HR-238, 1985

The authors have post-tensioned and monitored two Iowa bridges and have field tested the post-tensioning of a composite bridge in Florida. In order to provide the practical post-tensioning distribution factors given in this manual, the authors developed a finite element model of a composite bridge and checked the model against a one-half scale laboratory bridge and two actual composite bridges, one of which had a 45 deg skew. Following a brief discussion of this background research, this manual explains the use of elastic, composite beam and bridge section properties, the distribution fractions for symmetrically post-tensioned exterior beams, and a method for computing the strength of a post-tensioned beam. Also included is a design example for a typical, 51.25-ft (15.62-m) span, four-beam composite bridge. Moments for Iowa Department of Transportation rating trucks, H 20 and HS 20 trucks, have been tabulated for design convenience and are included in the appendix.

Field Monitoring of Curved Girder Bridges with Integral Abutments, TPF-5(169), 2014

Nationally, there are questions regarding the design, fabrication, and erection of horizontally curved steel girder bridges due to unpredicted girder displacements, fit-up, and locked-in stresses. One reason for the concerns is that up to one-quarter of steel girder bridges are being designed with horizontal curvature. There is also an urgent need to reduce bridge maintenance costs by eliminating or reducing deck joints, which can be achieved by expanding the use of integral abutments to include curved girder bridges. However, the behavior of horizontally curved bridges with integral abutments during thermal loading is not well known nor understood. The purpose of this study was to investigate the behavior of horizontal curved bridges with integral abutment (IAB) and semi-integral abutment bridges (SIAB) with a specific interest in the response to changing temperatures. The long-term objective of this effort is to establish guidelines for the use of integral abutments with curved girder bridges. The primary objective of this work was to monitor and evaluate the behavior of six in-service, horizontally curved, steel-girder bridges with integral and semi-integral abutments. In addition, the influence of bridge curvature, skew and pier bearing (expansion and fixed) were also part of the study. Two monitoring systems were designed and applied to a set of four horizontally curved bridges and two straight bridges at the northeast corner of Des Moines, Iowa—one system for measuring strains and movement under long term thermal changes and one system for measuring the behavior under short term, controlled live loading. A finite element model was developed and validated against the measured strains. The model was then used to investigate the sensitivity of design calculations to curvature, skew and pier joint conditions. The general conclusions were as follows: (1) There were no measurable differences in the behavior of the horizontally curved bridges and straight bridges studied in this work under thermal effects. For preliminary member sizing of curved bridges, thermal stresses and movements in a straight bridge of the same length are a reasonable first approximation. (2) Thermal strains in integral abutment and semi-integral abutment bridges were not noticeably different. The choice between IAB and SIAB should be based on life – cycle costs (e.g., construction and maintenance). (3) An expansion bearing pier reduces the thermal stresses in the girders of the straight bridge but does not appear to reduce the stresses in the girders of the curved bridge. (4) An analysis of the bridges predicted a substantial total stress (sum of the vertical bending stress, the lateral bending stress, and the axial stress) up to 3 ksi due to temperature effects. (5) For the one curved integral abutment bridge studied at length, the stresses in the girders significantly vary with changes in skew and curvature. With a 10⁰ skew and 0.06 radians arc span length to radius ratio, the curved and skew integral abutment bridges can be designed as a straight bridge if an error in estimation of the stresses of 10% is acceptable.

Nonlinear Pile Behavior in Integral Abutment Bridges, HR-227, Part 1, 1982

The highway departments of all fifty states were contacted to find the extent of application of integral abutment bridges, to survey the different guidelines used for analysis and design of integral abutment bridges, and to assess the performance of such bridges through the years. The variation in design assumptions and length limitations among the various states in their approach to the use of integral abutments is discussed. The problems associated with lateral displacements at the abutment, and the solutions developed by the different states for most of the ill effects of abutment movements are summarized in the report. An algorithm based on a state-of-the-art nonlinear finite element procedure was developed and used to study piling stresses and pile-soil interaction in integral abutment bridges. The finite element idealization consists of beam-column elements with geometric and material nonlinearities for the pile and nonlinear springs for the soil. An idealized soil model (modified Ramberg-Osgood model) was introduced in this investigation to obtain the tangent stiffness of the nonlinear spring elements. Several numerical examples are presented in order to establish the reliability of the finite element model and the computer software developed. Three problems with analytical solutions were first solved and compared with theoretical solutions. A 40 ft H pile (HP 10 X 42) in six typical Iowa soils was then analyzed by first applying a horizontal displacement (to simulate bridge motion) and no rotation at the top and then applying a vertical load V incrementally until failure occurred. Based on the numerical results, the failure mechanisms were generalized to be of two types: (a) lateral type failure and (b) vertical type failure. It appears that most piles in Iowa soils (sand, soft clay and stiff clay) failed when the applied vertical load reached the ultimate soil frictional resistance (vertical type failure). In very stiff clays, however, the lateral type failure occurs before vertical type failure because the soil is sufficiently stiff to force a plastic hinge to form in the pile as the specified lateral displacement is applied. Preliminary results from this investigation showed that the vertical load-carrying capacity of H piles is not significantly affected by lateral displacements of 2 inches in soft clay, stiff clay, loose sand, medium sand and dense sand. However, in very stiff clay (average blow count of 50 from standard penetration tests), it was found that the vertical load carrying capacity of the H pile is reduced by about 50 percent for 2 inches of lateral displacement and by about 20 percent for lateral displacement of 1 inch. On the basis of the preliminary results of this investigation, the 265-feet length limitation in Iowa for integral abutment concrete bridges appears to be very conservative.

Investigation of Two Bridge Alternatives for Low Volume Roads Volume 2 of 2 Concept 2: Beam In Slab Bridge, HR-382, 1997

Recent reports have indicated that 23.5% of the nation's highway bridges are structurally deficient and 17.7% are functionally obsolete. A significant number of these bridges are on the Iowa secondary road system where over 86% of the rural bridge management responsibilities are assigned to the counties. Some of the bridges can be strengthened or otherwise rehabilitated, but many more are in need of immediate replacement. In a recent investigation (HR-365 "Evaluation of Bridge Replacement Alternatives for the County Bridge System") several types of replacement bridges that are currently being used on low volume roads were identified. It was also determined that a large number of counties (69%) have the ability and are interested in utilizing their own forces to design and construct short span bridges. After reviewing the results from HR-365, the research team developed one "new" bridge replacement concept and a modification of a replacement system currently being used. Both of these bridge replacement alternatives were investigated in this study, the results of which are presented in two volumes. This volume (Volume 2) presents the results of Concept 2 - Modification of the Beam-in-Slab Bridge. Concept 1 - Steel Beam Precast Units is presented in Volume 1. Concept 2 involves various laboratory tests of the Beam-in-Slab bridge (BISB) currently being used by Benton County and several other Iowa counties. In this investigation, the behavior and strength of the BISB were determined; a new method of obtaining composite action between the steel beams and concrete was also tested. Since the Concept 2 bridge is primarily intended for use on low-volume roads, the system can be constructed with new or used beams. In the experimental part of the investigation, there were three types of laboratory tests: push-out tests, service and ultimate load tests of models of the BISB, and composite beam tests utilizing the newly developed shear connection. In addition to the laboratory tests, there was a field test in which an existing BISB was service load tested. An equation was developed for predicting the strength of the shear connection investigated; in addition, a finite element model for analyzing the BISB was also developed. Push-out tests were completed to determine the strength of the recently developed shear connector. A total of 36 specimens were tested, with variables such as hole diameter, hole spacing, presence of reinforcement, etc. being investigated. In the model tests of the BISB, two and four beam specimens [L=9,140 mm (30 ft)] were service load tested for behavior and load distribution data. Upon completion of these tests, both specimens were loaded to failure. In the composite beam tests, four beams, one with standard shear studs and three using the shear connection developed, were tested. Upon completion of the service load tests, all four beams were loaded to failure. The strength and behavior of the beams with the new shear connection were found to be essentially the same as that of the specimen with standard shear studs.