Model to describe the mode I fracture of steel fiber reinforced ultra-high performance concrete

Ultra-high performance fiber reinforced concrete (UHPFRC) has arisen from the implementation of a variety of concrete engineering and materials science concepts developed over the last century. This material offers superior strength, serviceability, and durability over its conventional counterparts. One of the most important differences for UHPFRC over other concrete materials is its ability to resist fracture through the use of randomly dispersed discontinuous fibers and improvements to the fiber-matrix bond. Of particular interest is the materials ability to achieve higher loads after first crack, as well as its high fracture toughness. In this research, a study of the fracture behavior of UHPFRC with steel fibers was conducted to look at the effect of several parameters related to the fracture behavior and to develop a fracture model based on a non-linear curve fit of the data. To determine this, a series of three-point bending tests were performed on various single edge notched prisms (SENPs). Compression tests were also performed for quality assurance. Testing was conducted on specimens of different cross-sections, span/depth (S/D) ratios, curing regimes, ages, and fiber contents. By comparing the results from prisms of different sizes this study examines the weakening mechanism due to the size effect. Furthermore, by employing the concept of fracture energy it was possible to obtain a comparison of the fracture toughness and ductility. The model was determined based on a fit to P-w fracture curves, which was cross referenced for comparability to the results. Once obtained the model was then compared to the models proposed by the AFGC in the 2003 and to the ACI 544 model for conventional fiber reinforced concretes.

Development of a continuum mechanics model of passive skeletal muscle

Skeletal muscle force evaluation is difficult to implement in a clinical setting. Muscle force is typically assessed through either manual muscle testing, isokinetic/isometric dynamometry, or electromyography (EMG). Manual muscle testing is a subjective evaluation of a patient’s ability to move voluntarily against gravity and to resist force applied by an examiner. Muscle testing using dynamometers adds accuracy by quantifying functional mechanical output of a limb. However, like manual muscle testing, dynamometry only provides estimates of the joint moment. EMG quantifies neuromuscular activation signals of individual muscles, and is used to infer muscle function. Despite the abundance of work performed to determine the degree to which EMG signals and muscle forces are related, the basic problem remains that EMG cannot provide a quantitative measurement of muscle force. Intramuscular pressure (IMP), the pressure applied by muscle fibers on interstitial fluid, has been considered as a correlate for muscle force. Numerous studies have shown that an approximately linear relationship exists between IMP and muscle force. A microsensor has recently been developed that is accurate, biocompatible, and appropriately sized for clinical use. While muscle force and pressure have been shown to be correlates, IMP has been shown to be non-uniform within the muscle. As it would not be practicable to experimentally evaluate how IMP is distributed, computational modeling may provide the means to fully evaluate IMP generation in muscles of various shapes and operating conditions. The work presented in this dissertation focuses on the development and validation of computational models of passive skeletal muscle and the evaluation of their performance for prediction of IMP. A transversly isotropic, hyperelastic, and nearly incompressible model will be evaluated along with a poroelastic model.

In vitro biomechanical testing of a micro external skeletal fixator

OBJECTIVE To biomechanically test the properties of three different Universal Micro External Fixator (UMEX™) configurations with regard to their use in very small animals (<5kg) and compare the UMEX system to the widely used IMEX External Skeletal Fixation (SK™) system in terms of stiffness, space needed for pin placement and weight. METHODS Three different UMEX configurations (type Ia, type Ib, and type II modified) and one SK configuration type Ia were used to stabilize Delrin plastic rods in a 1 cm fracture gap model. These constructs were tested in axial compression, craniocaudal bending, mediolateral bending, and torsion. Testing was conducted within the elastic range and mean stiffness in each mode was determined from the slope of the linear portion of the load-deformation curve. A Kruskal Wallis one-way analysis of variance on ranks test was utilized to assess differences between constructs (p <0.05). RESULTS The UMEX type II modified configuration was significantly stiffer than the other UMEX configurations and the SK type Ia, except in craniocaudal bending, where the SK type Ia configuration was stiffer than all UMEX constructs. The UMEX type Ia configuration was significantly the weakest of those frames. The UMEX constructs were lighter and smaller than the SK, thus facilitating closer pin placement. CONCLUSIONS Results supported previous reports concerning the superiority of more complex constructs regarding stiffness. The UMEX system appears to be a valid alternative for the treatment of long-bone fractures in very small animals.

Improved Localization of Implanted Subdural Electrode Contacts on MRI Using an Elastic Image Fusion Algorithm in Invasive EEG Recording

BACKGROUND: Accurate projection of implanted subdural electrode contacts in presurgical evaluation of pharmacoresistant epilepsy cases by invasive EEG is highly relevant. Linear fusion of CT and MRI images may display the contacts in the wrong position due to brain shift effects. OBJECTIVE: A retrospective study in five patients with pharmacoresistant epilepsy was performed to evaluate whether an elastic image fusion algorithm can provide a more accurate projection of the electrode contacts on the pre-implantation MRI as compared to linear fusion. METHODS: An automated elastic image fusion algorithm (AEF), a guided elastic image fusion algorithm (GEF), and a standard linear fusion algorithm (LF) were used on preoperative MRI and post-implantation CT scans. Vertical correction of virtual contact positions, total virtual contact shift, corrections of midline shift and brain shifts due to pneumencephalus were measured. RESULTS: Both AEF and GEF worked well with all 5 cases. An average midline shift of 1.7mm (SD 1.25) was corrected to 0.4mm (SD 0.8) after AEF and to 0.0mm (SD 0) after GEF. Median virtual distances between contacts and cortical surface were corrected by a significant amount, from 2.3mm after LF to 0.0mm after AEF and GEF (p<.001). Mean total relative corrections of 3.1 mm (SD 1.85) after AEF and 3.0mm (SD 1.77) after GEF were achieved. The tested version of GEF did not achieve a satisfying virtual correction of pneumencephalus. CONCLUSION: The technique provided a clear improvement in fusion of pre- and post-implantation scans, although the accuracy is difficult to evaluate.

Bone Volume Fraction and Fabric Anisotropy Are Better Determinants of Trabecular Bone Stiffness than Other Morphological Variables

As our population ages, more individuals suffer from osteoporosis. This disease leads to impaired trabecular architecture and increased fracture risk. It is essential to understand how morphological and mechanical properties of the cancellous bone are related. Morphologyelasticity relationships based on bone volume fraction (BV/TV) and fabric anisotropy explain up to 98% of the variation in elastic properties. Yet, other morphological variables such as individual trabeculae segmentation (ITS) and trabecular bone score (TBS) could improve the stiffness predictions. A total of 743 micro-computed tomography reconstructions of cubic trabecular bone samples extracted from femur, radius, vertebrae and iliac crest were analysed. Their morphology was assessed via 25 variables and their stiffness tensor (inline image) was computed from six independent load cases using micro finite element analyses. Variance inflation factors were calculated to evaluate collinearity between morphological variables and decide upon their inclusion in morphology-elasticity relationships. The statistically admissible morphological variables were included in a multi-linear regression modelling the dependent variable inline image. The contribution of each independent variable was evaluated (ANOVA). Our results show that BV/TV is the best determinant of inline image (inline image=0.889), especially in combination with fabric (inline image=0.968). Including the other independent predictors hardly affected the amount of variance explained by the model (inline image=0.975). Across all anatomical sites, BV/TV explained 87% of the variance of the bone elastic properties. Fabric further described 10% of the bone stiffness, but the improvement in variance explanation by adding other independent factors was marginal (<1%). These findings confirm that BV/TV and fabric are the best determinants of trabecular bone stiffness and show, against common belief, that other morphological variables do not bring any further contribution. These overall conclusions remain to be confirmed for specific bone diseases and post-elastic properties.

Development of a balanced experimental–computational approach to understanding the mechanics of proximal femur fractures

The majority of people who sustain hip fractures after a fall to the side would not have been identified using current screening techniques such as areal bone mineral density. Identifying them, however, is essential so that appropriate pharmacological or lifestyle interventions can be implemented. A protocol, demonstrated on a single specimen, is introduced, comprising the following components; in vitro biofidelic drop tower testing of a proximal femur; high-speed image analysis through digital image correlation; detailed accounting of the energy present during the drop tower test; organ level finite element simulations of the drop tower test; micro level finite element simulations of critical volumes of interest in the trabecular bone. Fracture in the femoral specimen initiated in the superior part of the neck. Measured fracture load was 3760 N, compared to 4871 N predicted based on the finite element analysis. Digital image correlation showed compressive surface strains as high as 7.1% prior to fracture. Voxel level results were consistent with high-speed video data and helped identify hidden local structural weaknesses. We found using a drop tower test protocol that a femoral neck fracture can be created with a fall velocity and energy representative of a sideways fall from standing. Additionally, we found that the nested explicit finite element method used allowed us to identify local structural weaknesses associated with femur fracture initiation.

End caps prevent nail migration in elastic stable intramedullary nailing in paediatric femoral fractures: a biomechanical study using synthetic and cadaveric bones.

End caps are intended to prevent nail migration (push-out) in elastic stable intramedullary nailing. The aim of this study was to investigate the force at failure with and without end caps, and whether different insertion angles of nails and end caps would alter that force at failure. Simulated oblique fractures of the diaphysis were created in 15 artificial paediatric femurs. Titanium Elastic Nails with end caps were inserted at angles of 45°, 55° and 65° in five specimens for each angle to create three study groups. Biomechanical testing was performed with axial compression until failure. An identical fracture was created in four small adult cadaveric femurs harvested from two donors (both female, aged 81 and 85 years, height 149 cm and 156 cm, respectively). All femurs were tested without and subsequently with end caps inserted at 45°. In the artificial femurs, maximum force was not significantly different between the three groups (p = 0.613). Push-out force was significantly higher in the cadaveric specimens with the use of end caps by an up to sixfold load increase (830 N, standard deviation (SD) 280 vs 150 N, SD 120, respectively; p = 0.007). These results indicate that the nail and end cap insertion angle can be varied within 20° without altering construct stability and that the risk of elastic stable intramedullary nailing push-out can be effectively reduced by the use of end caps.

A class of exact Navier-Stokes solutions for homogeneous flat-plate boundary layers and their linear stability

We introduce a new boundary layer formalism on the basis of which a class of exact solutions to the Navier–Stokes equations is derived. These solutions describe laminar boundary layer flows past a flat plate under the assumption of one homogeneous direction, such as the classical swept Hiemenz boundary layer (SHBL), the asymptotic suction boundary layer (ASBL) and the oblique impingement boundary layer. The linear stability of these new solutions is investigated, uncovering new results for the SHBL and the ASBL. Previously, each of these flows had been described with its own formalism and coordinate system, such that the solutions could not be transformed into each other. Using a new compound formalism, we are able to show that the ASBL is the physical limit of the SHBL with wall suction when the chordwise velocity component vanishes while the homogeneous sweep velocity is maintained. A corresponding non-dimensionalization is proposed, which allows conversion of the new Reynolds number definition to the classical ones. Linear stability analysis for the new class of solutions reveals a compound neutral surface which contains the classical neutral curves of the SHBL and the ASBL. It is shown that the linearly most unstable Görtler–Hämmerlin modes of the SHBL smoothly transform into Tollmien–Schlichting modes as the chordwise velocity vanishes. These results are useful for transition prediction of the attachment-line instability, especially concerning the use of suction to stabilize boundary layers of swept-wing aircraft.

Fluid-structure interaction of an elastic pipe immersed inside a coaxial rigid pipe: a model for the Tullio phenomenon

An axisymmetric, elastic pipe is filled with an incompressible fluid and is immersed in a second, coaxial rigid pipe which contains the same fluid. A pressure pulse in the outer fluid annulus deforms the elastic pipe which invokes a fluid motion in the fluid core. It is the aim of this study to investigate streaming phenomena in the core which may originate from such a fluid-structure interaction. This work presents a numerical solver for such a configuration. It was developed in the OpenFOAM software environment and is based on the Arbitrary Lagrangian Eulerian (ALE) approach for moving meshes. The solver features a monolithic integration of the one-dimensional, coupled system between the elastic structure and the outer fluid annulus into a dynamic boundary condition for the moving surface of the fluid core. Results indicate that our configuration may serve as a mechanical model of the Tullio Phenomenon (sound-induced vertigo).

Burrowing in Marine Muds by Crack Propagation: Kinematics and Forces

The polychaete Nereis virens burrows through muddy sediments by exerting dorsoventral forces against the walls of its tongue-depressor- shaped burrow to extend an oblate hemispheroidal crack. Stress is concentrated at the crack tip, which extends when the stress intensity factor (K-I) exceeds the critical stress intensity factor (K-Ic). Relevant forces were measured in gelatin, an analog for elastic muds, by photoelastic stress analysis, and were 0.015 +/- 0.001 N (mean +/- s.d.;N= 5). Measured elastic moduli (E) for gelatin and sediment were used in finite element models to convert the forces in gelatin to those required in muds to maintain the same body shapes observed in gelatin. The force increases directly with increasing sediment stiffness, and is 0.16 N for measured sediment stiffness of E=2.7x10(4) Pa. This measurement of forces exerted by burrowers is the first that explicitly considers the mechanical behavior of the sediment. Calculated stress intensity factors fall within the range of critical values for gelatin and exceed those for sediment, showing that crack propagation is a mechanically feasible mechanism of burrowing. The pharynx extends anteriorly as it everts, extending the crack tip only as far as the anterior of the worm, consistent with wedge-driven fracture and drawing obvious parallels between soft-bodied burrowers and more rigid, wedge-shaped burrowers(i.e. clams). Our results raise questions about the reputed high energetic cost of burrowing and emphasize the need for better understanding of sediment mechanics to quantify external energy expenditure during burrowing.

Physical properties and elastic constants of ODP Holes 137-504B and 140-504B upper crustal rocks

Seismic velocities have been measured at confining pressures of 100 MPa and 600 MPa for sheeted dike samples recovered during Ocean Drilling Program Legs 137 and 140. The compressional- and shear-wave velocities show an increase with depth at Hole 504B, which is in sharp contrast to the atmospheric pressure velocity measurements performed as part of the shipboard analyses. Rocks exposed to different types of alteration and fracture patterns show distinct changes in their physical properties. The seismic reflectors observed on the vertical seismic profile (VSP) experiment performed during Leg 111 may have been caused by low velocity zones resulting from alteration. The amount of fracturing and hydrothermal alteration in several zones also may have contributed to the acoustic impedance contrast necessary to produce the E5 reflector. Poisson's ratios calculated from laboratory velocity measurements show several low values at depths ranging from 1600 mbsf to 2000 mbsf, which tends to follow similar trends obtained from previous oceanic refraction experiments. A comparison of physical properties between samples recovered from Hole 504B and ophiolite studies in the Bay of Islands and Oman shows a good correlation with the Bay of Islands but significant differences from the measurements performed in the Oman complex.

Measurements of fracture toughness applying a four-point-bending technique of ice core B34 from Kohnen Station, Dronning Maud Land, Antarctica

The critical fracture toughness is a material parameter describing the resistance of a cracked body to further crack extension. It is an important parameter to simulate and predict the break-up behaviour of ice shelves from calving of single icebergs to the disintegration of entire ice shelves over a wide range of length scales. The fracture toughness values are calculated with equations that are derived from an elastic stress analysis. Additionally, an X-ray computer tomography (CT scanner) was used to identify the density as a function of depth. The critical fracture toughness of 91 Antarctic inland ice samples with densities between 840 to 870 kg/m**3 has been determined by applying a four-point-bending technique on single edge v-notched beam samples. The examined ice core was drilled 70m north of Kohnen Station, Dronnning Maud Land (75°00' S, 00°04' E, 2882 m).

Elastic (stress-strain) halo associated to ion-induced nano-tracks in lithium niobate: role of crystal anisotropy

The elastic strain/stress fields (halo) around a compressed amorphous nano-track (core) caused by a single high-energy ion impact on LiNbO3 are calculated. A method is developed to approximately account for the effects of crystal anisotropy of LiNbO3 (symmetry 3m) on the stress fields for tracks oriented along the crystal axes (X, Y or Z). It only considers the zero-order (axial) harmonic contribution to the displacement field in the perpendicular plane and uses effective Poisson moduli for each particular orientation. The anisotropy is relatively small; however, it accounts for some differential features obtained for irradiations along the crystallographic axes X, Y and Z. In particular, the irradiation-induced disorder (including halo) and the associated surface swelling appear to be higher for irradiations along the X- or Y-axis in comparison with those along the Z-axis. Other irradiation effects can be explained by the model, e.g. fracture patterns or the morphology of pores after chemical etching of tracks. Moreover, it offers interesting predictions on the effect of irradiation on lattice parameters

Linear global instability of non-orthogonal incompressible swept attachment-line boundary layer flow

Instability of the orthogonal swept attachment line boundary layer has received attention by local1, 2 and global3–5 analysis methods over several decades, owing to the significance of this model to transition to turbulence on the surface of swept wings. However, substantially less attention has been paid to the problem of laminar flow instability in the non-orthogonal swept attachment-line boundary layer; only a local analysis framework has been employed to-date.6 The present contribution addresses this issue from a linear global (BiGlobal) instability analysis point of view in the incompressible regime. Direct numerical simulations have also been performed in order to verify the analysis results and unravel the limits of validity of the Dorrepaal basic flow7 model analyzed. Cross-validated results document the effect of the angle _ on the critical conditions identified by Hall et al.1 and show linear destabilization of the flow with decreasing AoA, up to a limit at which the assumptions of the Dorrepaal model become questionable. Finally, a simple extension of the extended G¨ortler-H¨ammerlin ODE-based polynomial model proposed by Theofilis et al.4 is presented for the non-orthogonal flow. In this model, the symmetries of the three-dimensional disturbances are broken by the non-orthogonal flow conditions. Temporal and spatial one-dimensional linear eigenvalue codes were developed, obtaining consistent results with BiGlobal stability analysis and DNS. Beyond the computational advantages presented by the ODE-based model, it allows us to understand the functional dependence of the three-dimensional disturbances in the non-orthogonal case as well as their connections with the disturbances of the orthogonal stability problem.

The scenario of two-dimensional instabilities of the cylinder wake under EHD forcing: A linear stability analysis

We propose to study the stability properties of an air flow wake forced by a dielectric barrier discharge (DBD) actuator, which is a type of electrohydrodynamic (EHD) actuator. These actuators add momentum to the flow around a cylinder in regions close to the wall and, in our case, are symmetrically disposed near the boundary layer separation point. Since the forcing frequencies, typical of DBD, are much higher than the natural shedding frequency of the flow, we will be considering the forcing actuation as stationary. In the first part, the flow around a circular cylinder modified by EHD actuators will be experimentally studied by means of particle image velocimetry (PIV). In the second part, the EHD actuators have been numerically implemented as a boundary condition on the cylinder surface. Using this boundary condition, the computationally obtained base flow is then compared with the experimental one in order to relate the control parameters from both methodologies. After validating the obtained agreement, we study the Hopf bifurcation that appears once the flow starts the vortex shedding through experimental and computational approaches. For the base flow derived from experimentally obtained snapshots, we monitor the evolution of the velocity amplitude oscillations. As to the computationally obtained base flow, its stability is analyzed by solving a global eigenvalue problem obtained from the linearized Navier–Stokes equations. Finally, the critical parameters obtained from both approaches are compared.