Nonlinear quasi-static finite element simulations predict in vitro strength of human proximal femora assessed in a dynamic sideways fall setup

Osteoporotic proximal femur fractures are caused by low energy trauma, typically when falling on the hip from standing height. Finite element simulations, widely used to predict the fracture load of femora in fall, usually include neither mass-related inertial effects, nor the viscous part of bone's material behavior. The aim of this study was to elucidate if quasi-static non-linear homogenized finite element analyses can predict in vitro mechanical properties of proximal femora assessed in dynamic drop tower experiments. The case-specific numerical models of thirteen femora predicted the strength (R2=0.84, SEE=540 N, 16.2%), stiffness (R2=0.82, SEE=233 N/mm, 18.0%) and fracture energy (R2=0.72, SEE=3.85 J, 39.6%); and provided fair qualitative matches with the fracture patterns. The influence of material anisotropy was negligible for all predictions. These results suggest that quasi-static homogenized finite element analysis may be used to predict mechanical properties of proximal femora in the dynamic sideways fall situation.

Anelastic strain recovery and elastic properties of ODP Leg 123 oceanic basalt samples

A knowledge of rock stress is fundamental for improving our understanding of oceanic crustal mechanisms and lithospheric dynamic processes. However, direct measurements of stress in the deep oceans, and in particular stress magnitudes, have proved to be technically difficult. Anelastic strain recovery measurements were conducted on 15 basalt core samples from Sites 765 and 766 during Leg 123. Three sets of experiments were performed: anelastic strain recovery monitoring, dynamic elastic property measurements, and thermal azimuthal anisotropy observations. In addition, a range of other tests and observations were recorded to characterize each of the samples. One common feature of the experimental results and observations is that apparently no consistent orientation trend exists, either between the different measurements on each core sample or between the same sets of measurements on the various core samples. However, some evidence of correspondence between velocity anisotropy and anelastic strain recovery exists, but this is not consistent for all the core samples investigated. Thermal azimuthal anisotropy observations, although showing no conclusive correlations with the other results, were of significant interest in that they clearly exhibited anisotropic behavior. The apparent reproducibility of this behavior may point toward the possibility of rocks that retain a "memory" of their stress history, which could be exploited to derive stress orientations from archived core. Anelastic strain recovery is a relatively new technique. Because use of the method has extended to a wider range of rock types, the literature has begun to include examples of rocks that contracted with time. Strong circumstantial evidence exists to suggest that core-sample contractions result from the slow diffusion of pore fluids from a preexisting microcrack structure that permits the rock to deflate at a greater rate than the expansion caused by anelastic strain recovery. Both expansions and contractions of the Leg 123 cores were observed. The basalt cores have clearly been intersected by an abundance of preexisting fractures, some of which pass right through the samples, but many are intercepted or terminate within the rock matrix. Thus, the behavior of the core samples will be influenced not only by the properties of the rock matrix between the fractures, but also by how these macro- and micro-scale fractures mutually interact. The strain-recovery curves recorded during Leg 123 for each of the 15 basalt core samples may reflect the result of two competing time dependent processes: anelastic strain recovery and pore pressure recovery. Were these the only two processes to influence the gauge responses, then one might expect that given the additional information required, established theoretical models might be used to determine consistent stress orientations and reliable stress magnitudes. However, superimposed upon these competing processes is their respective interaction with the preexisting fractures that intersect each core. Evidence from our experiments and observations suggests that these fractures have a dominating influence on the characteristics of the recovery curves and that their effects are complex.

Dynamic response of hydro power plants for providing secondary regulation reserves considering elastic water column effects

In this paper, the dynamic response of a hydro power plant for providing secondary regulation reserve is studied in detail. Special emphasis is given to the elastic water column effects both in the penstock and the tailrace tunnel. For this purpose, a nonlinear model based on the analogy between mass and momentum conservation equations of a water conduit and those of wave propagation in transmission lines is used. The influence of the plant configuration and design parameters on the fulfilment of the Spanish Electrical System Operator requirements is analysed

Nonlinear response to ultrasound of encapsulated microbubbles

The acoustic backscatter of encapsulated gas-filled microbubbles immersed in a weak compressible liquid and irradiated by ultrasound fields of moderate to high pressure amplitudes is investigated theoretically. The problem is formulated by considering, for the viscoelastic shell of finite thickness, an isotropic hyperelastic neo-Hookean model for the elastic contribution in addition to a Newtonian viscous component. First and second harmonic scattering cross-sections have been evaluated and the quantitative influence of the driving pressure amplitude on the harmonic resonance frequencies for different initial equilibrium bubble sizes and for different encapsulating physical properties has been determined. Conditions for optimal second harmonic imaging have been also investigated and some regions in the parameters space where the second harmonic intensity is dominant over the fundamental have been identified. Results have been obtained for albumin, lipid and polymer encapsulating shells, respectively.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.

A mathematical framework for new fault detection schemes in nonlinear stochastic continuous-time dynamical systems

n this work, a mathematical unifying framework for designing new fault detection schemes in nonlinear stochastic continuous-time dynamical systems is developed. These schemes are based on a stochastic process, called the residual, which reflects the system behavior and whose changes are to be detected. A quickest detection scheme for the residual is proposed, which is based on the computed likelihood ratios for time-varying statistical changes in the Ornstein–Uhlenbeck process. Several expressions are provided, depending on a priori knowledge of the fault, which can be employed in a proposed CUSUM-type approximated scheme. This general setting gathers different existing fault detection schemes within a unifying framework, and allows for the definition of new ones. A comparative simulation example illustrates the behavior of the proposed schemes.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.

Coherent excitation of nonlinear oscillators

Coherently driven, dissipative nonlinear oscillators,(driving kept permanently in phase with the oscillators response) are proposed as systems with interesting dynamics. Results for simple, preliminary examples, which do not show chaotic behavior, are briefly discussed.

Nonlinear Dynamics in experimental devices with compressed/expanded surfactant monolayers

A theory is provided for a common experimental set up that is used to measure surface properties in surfactant monolayers. The set up consists of a surfactant monolayer (over a shallow liquid layer) that is compressed/expanded in a periodic fashion by moving in counter-phase two parallel, slightly immersed solid barriers, which vary the free surface area and thus the surfactant concentration. The simplest theory ignores the fluid dynamics in the bulk fluid, assuming spatially uniform surfactant concentration, which requires quite small forcing frequencies and provides reversible dynamics in the compression/expansion cycles. In this paper, we present a long-wave theory for not so slow oscillations that assumes local equilibrium but takes the fluid dynamics into account. This simple theory uncovers the physical mechanisms involved in the surfactant behavior and allows for extracting more information from each experimental run. The conclusion is that the fluid dynamics cannot be ignored, and that some irreversible dynamics could well have a fluid dynamic origin

A nonlinear analysis of laying a floating pipeline on the seabed

In this article, a model for the determination of displacements, strains, and stresses of a submarine pipeline during its construction is presented. Typically, polyethylene outfall pipelines are the ones treated by this model. The process is carried out from an initial floating situation to the final laying position on the seabed. The following control variables are considered in the laying process: the axial load in the pipe, the flooded inner length, and the distance of the control barge from the coast. External loads such as self-weight, dead loads, and forces due to currents and small waves are also taken into account.This paper describes both the conceptual framework for the proposed model and its practical application in a real engineering situation. The authors also consider how the model might be used as a tool to study how sensitive the behavior of the pipeline is to small changes in the values of the control variables. A detailed description of the actions is considered, especially the ones related to the marine environment such as buoyancy, current, and sea waves. The structural behavior of the pipeline is simulated in the framework of a geometrically nonlinear dynamic analysis. The pipeline is assumed to be a two-dimensional Navier_Bernoulli beam. In the nonlinear analysis an updated Lagrangian formulation is used, and special care is taken regarding the numerical aspects of sea bed contact, follower forces due to external water pressures, and dynamic actions. The paper concludes by describing the implementation of the proposed techniques, using the ANSYS computer program with a number of subroutines developed by the authors. This implementation permits simulation of the two-dimensional structural pipe behavior of the whole construction process. A sensitivity analysis of the bending moments, axial forces, and stresses for different values of the control variables is carried out. Using the techniques described, the engineer may optimize the construction steps in the pipe laying process

Mechanical behavior and failure micromechanisms of hybrid 3D woven composites in tension

The deformation and failure micromechanisms of a hybrid 3D woven composite were studied in tension. Plain and open-hole composite coupons were tested in tension until failure in the fill and warp directions, as well as fiber tows extracted from the dry fabric and impregnated with the matrix. The macroscopic evolution of damage in the composite coupons was assessed by means of periodic unloading–reloading (to obtain the elastic modulus and the residual strain), whereas the microscopic mechanism were established by means of X-ray computed microtomography. To this end, specimens were periodically removed from the mechanical testing machine and infiltrated with ZnI-containing liquid to assess the main damage modes as a function of the applied strain. The experimental observations and the predictions of an isostrain model were used to understand the key factors controlling the elastic modulus, strength and notch sensitivity of hybrid 3D woven composites in tension. It was found that the full contribution of the glass fibers to the composite strength was not employed, due to the premature fracture of the carbon fibers, but their presence increased the fracture strain and the energy dissipated during fracture. Thus, hybridization of the 3D woven composite led to a notch-insensitive behavior as demonstrated by open-hole tests

Mechanical behavior and deformation micromechanisms of polypropylene nonwoven fabrics as a function of temperature and strain rate

The mechanical behavior and the deformation and failure micromechanisms of a thermally-bonded polypropylene nonwoven fabric were studied as a function of temperature and strain rate. Mechanical tests were carried out from 248 K (below the glass transition temperature) up to 383 K at strain rates in the range ≈10−3 s−1 to 10−1 s−1. In addition, individual fibers extracted from the nonwoven fabric were tested under the same conditions. Micromechanisms of deformation and failure at the fiber level were ascertained by means of mechanical tests within the scanning electron microscope while the strain distribution at the macroscopic level upon loading was determined by means of digital image correlation. It was found that the nonwoven behavior was mainly controlled by the properties of the fibers and of the interfiber bonds. Fiber properties determined the nonlinear behavior before the peak load while the interfiber bonds controlled the localization of damage after the peak load. The influence of these properties on the strength, ductility and energy absorbed during deformation is discussed from the experimental observations.

Railway Rapid Transit timetables with variable and elastic demand.

This paper focuses on the design of railway timetables considering a variable elastic demand profile along a whole design day. Timetabling is the third stage in the classical hierarchical railway planning process. Most of previous works on this topic consider a uniform demand behavior for short planning intervals. In this paper, we propose a MINLP model for designing non-periodic timetables on a railway corridor where demand is dependent on waiting times. In the elastic demand case, long waiting times lead to a loss of passengers, who may select an alternative transportation mode. The mode choice is modeled using two alternative methods. The first one is based on a sigmoid function and can be used in case of absence of information for competitor modes. In the second one, the mode choice probability is obtained using a Logit model that explicitly considers the existence of a main alternative mode. With the purpose of obtaining optimal departure times, in both cases, a minimization of the loss of passengers is used as objective function. Finally, as illustration, the timetabling MINLP model with both mode choice methods is applied to a real case and computational results are shown.

Numerical methods in Nonlinear Analysis of shell structures

The design of shell and spatial structures represents an important challenge even with the use of the modern computer technology.If we concentrate in the concrete shell structures many problems must be faced,such as the conceptual and structural disposition, optimal shape design, analysis, construction methods, details etc. and all these problems are interconnected among them. As an example the shape optimization requires the use of several disciplines like structural analysis, sensitivity analysis, optimization strategies and geometrical design concepts. Similar comments can be applied to other space structures such as steel trusses with single or double shape and tension structures. In relation to the analysis the Finite Element Method appears to be the most extended and versatile technique used in the practice. In the application of this method several issues arise. First the derivation of the pertinent shell theory or alternatively the degenerated 3-D solid approach should be chosen. According to the previous election the suitable FE model has to be adopted i.e. the displacement,stress or mixed formulated element. The good behavior of the shell structures under dead loads that are carried out towards the supports by mainly compressive stresses is impaired by the high imperfection sensitivity usually exhibited by these structures. This last effect is important particularly if large deformation and material nonlinearities of the shell may interact unfavorably, as can be the case for thin reinforced shells. In this respect the study of the stability of the shell represents a compulsory step in the analysis. Therefore there are currently very active fields of research such as the different descriptions of consistent nonlinear shell models given by Simo, Fox and Rifai, Mantzenmiller and Buchter and Ramm among others, the consistent formulation of efficient tangent stiffness as the one presented by Ortiz and Schweizerhof and Wringgers, with application to concrete shells exhibiting creep behavior given by Scordelis and coworkers; and finally the development of numerical techniques needed to trace the nonlinear response of the structure. The objective of this paper is concentrated in the last research aspect i.e. in the presentation of a state-of-the-art on the existing solution techniques for nonlinear analysis of structures. In this presentation the following excellent reviews on this subject will be mainly used.

Parametric analysis on bridge reinforced concrete elevated pile-cap foundations

Current design practices recommend to comply with the capacity protection principle, which pays special attention to ensuring an elastic response of the foundations under ground motion events. However, in cases such as elevated reinforced concrete (RC) pile-cap foundation typologies, this design criterion may lead to conservative designs, with excessively high construction costs. Reinforced concrete elevated pile-cap foundations is a system formed by a group of partially embedded piles connected through an aboveground stayed cap and embedded in soil. In the cases when they are subjected to ground motions, the piles suffer large bending moments that make it difficult to maintain their behavior within the elastic range of deformations. Aiming to make an in-depth analysis of the nonlinear behavior of elevated pile-cap foundations, a cyclic loading test was performed on a concrete 2x3 pile configuration specimen of elevated pile-cap foundation. Two results of this test, the failure mechanism and the ductile behavior, were used for the calibration of a numerical model built in OpenSees framework, by using a pushover analysis. The calibration of the numerical model enabled an in-depth study of the seismic nonlinear response of this kind of foundations. A parametric analysis was carried for this purpose, aiming to study how sensitive RC elevated pile-cap foundations are, when subjected to variations in the diameter of piles, reinforcement ratios, external loads, soil density or multilayer configurations. This analysis provided a set of ductility factors that can be used as a reference for design practices and which correspond to each of the cases analyzed.