Superelastic shape memory alloy cables: Part II – Subcomponent isothermal responses

This paper constitutes the second part of our experimental study of the thermo-mechanical behavior of superelastic NiTi shape memory alloy cables. Part I introduced the fundamental, room temperature, tensile responses of two cable designs (7 × 7 right regular lay, and 1 × 27 alternating lay). In Part II, each cable behavior is studied further by breaking down the response into the contributions of its hierarchical subcomponents. Selected wire strands were extracted from the two cable constructions, and their quasi-static tension responses were measured using the same experimental setup of Part I. Consistent with the shallow wire helix angles in the 7 × 7 construction, the force–elongation responses of the core wire, 1 × 7 core strand and full 7 × 7 cable were similar on a normalized basis, with only a slight decrease in transformation force plateaus and slight increase in plateau strains in this specimen sequence. By contrast, each successive 1 × 27 component (1 × 6 core strand, 1 × 15 strand, and full cable) included an additional outer layer of wires with a larger number of wires, greater helix radius, and deeper helix angle, so the normalized axial load responses became significantly more compliant. Each specimen in the sequence also exhibited progressively larger strains at failure, reaching 40% strain in the full 1 × 27 cable.Stress-induced phase transformations involved localized strain/temperature and front propagation in all of the tested 7 × 7 components but none of the 1 × 27 components aside from the 1 × 27 core wire. Stereo digital image correlation measurements revealed finer features within a global transformation front of the 1 × 7 core strand than the 7 × 7 cable, consisting of an staggered pattern of individual wire fronts that moved in lock-step during elongation. Although the 1 × 27 multi-layer strands exhibited temperature/strain localizations in a distributed pattern during transformations, the localizations did not propagate and their cause was traced back to contact indentations (stress concentrations) arising from the cable’s fabrication. The normalized axial torque responses of the multi-layer 1 × 27 components during transformation were distinctly non-monotonic and complex, due to the alternating handedness of the layers. Force and torque contributions of individual wire layers were deduced by subtracting 1 × 27 component responses, which helped to clarify the transformation kinetics within each layer and explain the unusual force and torque undulations seen in the 1 × 27 cable response of Part I.     

Measurement of dynamic surface tension for liquid metal by capillary jet method

Most liquid metals are easily oxidized, causing their surface tensions to vary over time. In this study, a novel method for measuring the dynamic surface tension of a liquid metal during oxidization is proposed. This method is based on the fact that the descent trajectory of a capillary jet ejected horizontally from a small orifice at slow speed depends on the surface tension as well as inertial and gravitational forces. We derive a theoretical model to predict the jet trajectory and determine the dynamic surface tension of the liquid metal by matching the predicted theoretical trajectory to the experimentally measured one. Actual measurements for Wood's alloy at various oxygen concentrations demonstrate that the surface tension decreases to an equilibrium value on a time scale of 10–11 ms at oxygen densities of less than 1.5%, whereas it increases on a time scale of 10–17 ms at greater oxygen densities. This method has the ability to measure dynamic surface tension on a time scale of 1.5–55 ms at an accuracy of ±1.5 ms, making it suitable for measuring the dynamic surface tension of liquid metal in an oxidation process.     

Buckling behavior of Kelvin open-cell foams under [0 0 1], [0 1 1] and [1 1 1] compressive loads

This paper describes buckling modes and stresses of elastic Kelvin open-cell foams subjected to [0 0 1], [0 1 1] and [1 1 1] uniaxial compressions. Cubic unit cells and cell aggregates in model foams are analyzed using a homogenization theory of the updated Lagrangian type. The analysis is performed on the assumption that the struts in foams have a non-uniform distribution of cross-sectional areas as observed experimentally. The relative density is changed to range from 0.005 to 0.05. It is thus found that long wavelength buckling and macroscopic instability primarily occur under [0 0 1] and [0 1 1] compressions, with only short wavelength buckling under [1 1 1] compression. The primary buckling stresses under the three compressions are fairly close to one another and almost satisfy the Gibson–Ashby relation established to fit experiments. By also performing the analysis based on the uniformity of strut cross-sectional areas, it is shown that the non-uniformity of cross-sectional areas is an important factor for the buckling behavior of open-cell foams.     

Analytical modeling of synthetic fiber ropes subjected to axial loads. Part I: A new continuum model for multilayered fibrous structures

Synthetic fiber ropes are characterized by a very complex architecture and a hierarchical structure. Considering the fiber rope architecture, to pass from fiber to rope structure behavior, two scale transition models are necessary, used in sequence: one is devoted to an assembly of a large number of twisted components (multilayered), whereas the second is suitable for a structure with a central straight core and six helical wires (1 + 6). The part I of this paper first describes the development of a model for the static behavior of a fibrous structure with a large number of twisted components. Tests were then performed on two different structures subjected to axial loads. Using the model presented here the axial stiffness of the structures has been predicted and good agreement with measured values is obtained. A companion paper (Ghoreishi, S.R. et al., in press. Analytical modeling of synthetic fiber ropes, part II: A linear elastic model for 1 + 6 fibrous structures, International Journal of Solids and Structures, doi:10.1016/j.ijsolstr.2006.08.032) presents the second model to predict the mechanical behavior of a 1 + 6 fibrous structure.     

Superelastic shape memory alloy cables: Part I – Isothermal tension experiments

Cables (or wire ropes) made from NiTi shape memory alloy (SMA) wires are relatively new and unexplored structural elements that combine many of the advantages of conventional cables with the adaptive properties of SMAs (shape memory and superelasticity) and have a broad range of potential applications. In this two part series, an extensive set of uniaxial tension experiments was performed on two Nitinol cable constructions, a 7 × 7 right regular lay and a 1 × 27 alternating lay, to characterize their superelastic behavior in room temperature air. Details of the evolution of strain and temperature fields were captured by simultaneous stereo digital image correlation and infrared imaging, respectively. Here in Part I, the nearly isothermal, superelastic responses of the two cable designs are presented and compared. Overall, the 7 × 7 construction has a mechanical response similar to that of straight wires with propagating transformation fronts and distinct stress plateaus during stress-induced transformations. The 1 × 27 construction, however, exhibits a more compliant and stable mechanical response, trading a decreased force for additional elongation, and does not exhibit transformation fronts due to the deeper helix angles of the layers. In Part II that follows, selected subcomponents are dissected from the two cable’s hierarchical constructions to experimentally break down the cable’s responses.     

Large-eddy simulations of chevron jet flows with noise predictions

Hybrid large-eddy type simulations for chevron nozzle jet flows are performed at Mach 0.9 and Re = 1.03 × 10⁶. Without using any subgrid scale model (SGS), the numerical approach applied in the present study is essentially implicit large-eddy simulation (ILES). However, a Reynolds-averaged Navier–Stokes (RANS) solution is patched into the near wall region. This makes the overall solution strategy hybrid RANS–ILES. The disparate turbulence length scales, implied by these different modeling approaches, are matched using a Hamilton–Jacobi equation. The complex geometry features of the chevron nozzles are fully meshed. With numerical fidelity in mind, high quality, hexahedral multi-block meshes of 12.5 × 10⁶ cells are used. Despite the modest meshes, the novel RANS–ILES approach shows encouraging performance. Computed mean and second-order fluctuating quantities of the turbulent near field compare favorably with measurements. The radiated far-field sound is predicted using the Ffowcs Williams and Hawkings (FW–H) surface integral method. Encouraging agreement of the predicted far-field sound directivity and spectra with measurements is obtained.     

The high strain mechanical response of the wet Kelvin model for open-cell foams

Surface Evolver software was used to create the three-dimensional geometry of a Kelvin open-cell foam, to simulate that of polyurethane flexible foams. Finite Element Analysis (FEA) with 3D elements was used to model large compressive deformation in the [0 0 1] and [1 1 1] directions, using cyclic boundary conditions when necessary, treating the polyurethane as an elastic or elastic–plastic material. The predicted foam Young’s moduli in the [0 0 1] direction are double those of foams with uniform Plateau border cross-section edges, for the same foam density and material properties. For compression in the [1 1 1] direction, the normalized Young’s modulus increases from 0.9 to 1.1 with foam relative density, and the predicted stress–strain relationship can have a plateau, even for a linearly-elastic polymer. As the foam density increases, the predicted effects of material plasticity become larger. For foam of relative density 0.028, edge-to-edge contact is predicted to occur at a 66% strain for [1 1 1] direction compression. The foam is predicted to contract laterally when the [1 1 1] direction compressive strain exceeds 25%.     

Fatigue of polycrystalline copper with different grain sizes and texture

Fatigue experiments of polycrystalline copper with different grain sizes and texture were conducted under tension–compression, torsion, and non-proportional loading. The grain sizes ranged from 10 μm to 2 mm. The stress–strain response was found to be a strong function of the grain size and texture. A plasticity-based critical plane multiaxial fatigue criterion was used to predict the fatigue lives of the polycrystalline copper. It was found that the criterion was able to correlate all the experimental results with one single set of material constants. This indicates that the fatigue failure of the material under consideration is dominated by the fatigue resistance of the grains with an insignificant influence of the grain boundaries on the fatigue of the polycrystalline material. It was found that the fatigue model with the material constants obtained from fatigue experiments can be applied to predict failure under monotonic torsion.     

Modeling and measuring visco-elastic properties: From collagen molecules to collagen fibrils

Collagen is the main structural protein in vertebrate biology, determining the mechanical behavior of connective tissues such as tendon, bone and skin. Although extensive efforts in the study of the origin of collagen exceptional mechanical properties, a deep knowledge of the relationship between molecular structure and mechanical properties remains elusive, hindered by the complex hierarchical structure of collagen-based tissues. Understanding the viscoelastic behavior of collagenous tissues requires knowledge of the properties at each structural level. Whole tissues have been studied extensively, but less is known about the mechanical behavior at the submicron, fibrillar and molecular level. Hence, we investigate the viscoelastic properties at the molecular level by using an atomistic modeling approach, performing in silico creep tests of a collagen-like peptide. The results are compared with creep and relaxation tests at the level of isolated collagen fibrils performed previously using a micro-electro-mechanical systems platform. Individual collagen molecules present a non-linear viscoelastic behavior, with a Young's modulus increasing from 6 to 16 GPa (for strains up to 20%), a viscosity of 3.84±0.38 Pa s, and a relaxation time in the range of 0.24–0.64 ns. At the fibrils level, stress–strain–time data indicate that isolated fibrils exhibit viscoelastic behavior that could be fitted using the Maxwell–Weichert model. The fibrils showed an elastic modulus of 123±46 MPa. The time-dependent behavior was well fit using the two-time-constant Maxwell–Weichert model with a fast time response of 7±2 s and a slow time response of 102±5 s.     

Nanoparticle debonding strength: A comprehensive study on interfacial effects

The most appealing feature of nanofilled polymers is the perspective of obtaining surprisingly high mechanical properties at low nanofiller volume fractions. The knowledge of nanostructure–property relationships is however essential for the design of these materials.In the present work, a model for the critical hydrostatic tension related to nanoparticle debonding is presented. The model accounts for some important issues inherently related to the nanoscale with particular reference to surface elastic stresses on the nanoparticle periphery and the emergence of a zone of altered chemistry surrounding the nanoparticle. The analytical solution suggests that the range of nanoparticle radii where interfacial effects do affect the solution is limited to the nanometer scale. In more details, considering the interphase and surface elastic properties used in the analysis, it has been found that for stiff particles with radius between 10 nm and 100 nm (silica, alumina and other metal oxide nanoparticles) the prominent role is played by the interphase elastic properties. Surface elastic constants were found to have, instead, only a negligible effect.     

Temperature dependent bulge test for elastomers

The bulge test is a particularly convenient testing method for characterizing elastomers under biaxial loading. In addition, it is convenient to utilize this test for validating material models in simulation due to the heterogeneous strain field induced during inflation. During the bulge test the strain field for elastomers covers uniaxial tension at the border to pure shear and equibiaxial tension at the pole. Elastomeric materials exhibit a hyperelastic material behavior, with a dependency on temperature and loading rate. The temperature effect on the mechanical behavior during biaxial loading is considered in the present study. A bulge test setup combined with a temperature chamber is developed in order to characterize this effect, and an exemplary temperature dependent characterization of a poly(norbornene) elastomer is performed with this setup. The equibiaxial stress–strain curves measured at 60 °C, 20 °C and −20 °C are presented.     

Resolving the stratification discrepancy of turbulent natural convection in differentially heated air-filled cavities. Part III: A full convection–conduction–surface radiation coupling

The present study concerns an air-filled differentially heated cavity of 1 m × 0.32 m × 1 m (width × depth × height) subject to a temperature difference of 15 K and is motivated by the need to understand the persistent discrepancy observed between numerical and experimental results on thermal stratification in the cavity core. An improved experiment with enhanced metrology was set up and experimental data have been obtained along with the characteristics of the surfaces and materials used. Experimental temperature distributions on the passive walls have been introduced in numerical simulations in order to provide a faithful prediction of experimental data. By means of DNS using spectral methods, heat conduction in the insulating material is first coupled with natural convection in the cavity. As heat conduction influences only the temperature distribution on the top and bottom surfaces and in the near wall regions, surface radiation is added to the coupling of natural convection with heat conduction. The temperature distribution in the cavity is strongly affected by the polycarbonate front and rear walls of the cavity, which are almost black surfaces for low temperature radiation, and also other low emissivity walls. The thermal stratification is considerably weakened by surface radiation. Good agreement between numerical simulations and experiments is observed on both time-averaged fields and turbulent statistics. Treating the full conduction–convection–radiation coupling allowed to confirm that experimental wall temperatures resulted from the coupled phenomena and this is another way to predict correctly the experimental results in the cavity.     

A review of new wire arrays with open and closed magnetic configurations at the 1.6 MA Zebra generator for radiative properties and opacity effects

The studies emphasize investigation of plasma formation, implosion, and radiation features as a function of two load configurations: compact multi-planar and cylindrical wire arrays. Experiments with different Z-pinch loads were performed on 1.6 MA, 100 ns, Zebra generator at University of Nevada, Reno. The multi-planar wire arrays (PWAs) were studied in open and closed configurations with Al, Cu, brass, Mo and W wires. In the open magnetic configurations (single, double, triple PWAs) magnetic fields are present inside the arrays from the beginning of discharge, while in closed configurations (prism-like PWA) the global magnetic field is excluded inside before plasma flow occurs. The new prism-like PWA allows high flexibility in control of implosion dynamics and precursor formation. The spectral modeling, magneto-hydrodynamic (MHD) and wire ablation dynamic model (WADM) codes were used to describe the plasma evolution and plasma parameters. Experimentally observed electron temperature and density in multiple bright spots reached 1.4 keV and 5 × 10²¹ cm^?3, respectively. Two types of bright spots were observed. With peak currents up to 1.3 MA opacity effects became more pronounced and led to a limiting of the X-ray yields from compact cylindrical arrays. Despite different magnetic energy to plasma coupling mechanisms early in the implosion a comparison of compact double PWA and cylindrical WA results indicates that during the stagnation stage the same plasma heating mechanism may occur. The double PWA was found to be the best radiator tested at University scale 1 MA generator. It is characterized by a combination of larger yield and power, mm-scale size, and provides the possibility of radiation pulse shaping. Further, the newer configuration, the double PWA with skewed wires, was tested and showed the possibility of a more effective X-ray generation.     

Mechanics of the rate-dependent elastic–plastic deformation of glassy polymers from low to high strain rates

A combined experimental and analytical investigation has been performed to understand the mechanical behavior of two amorphous polymers—polycarbonate and poly(methyl methacrylate)—at strain rates ranging from 10⁻⁴ to 10⁴ s⁻¹. This range in strain rates was achieved in uniaxial tension and compression tests using a dynamic mechanical analyzer (DMA), a servo-hydraulic testing machine, and an aluminum split-Hopkinson pressure bar. DMA tension tests were used to characterize the viscoelastic behavior of these materials, with focus on the rate-dependent shift of material transition temperatures. Uniaxial compression tests on the servo-hydraulic machine (10⁻⁴ to 1 s⁻¹) and the split-Hopkinson pressure bar (10³ to 10⁴ s⁻¹) were used to characterize the rate-dependent yield and post-yield behavior. Both materials were observed to exhibit increased rate sensitivity of yield under the same strain rate/temperature conditions as the β-transition of the viscoelastic behavior. A physically based constitutive model for large strain deformation of thermoplastics was then extended to encompass high-rate conditions. The model accounts for the contributions of different molecular motions which become operational and important in different frequency regimes. The new features enable the model to not only capture the transition in the yield behavior, but also accurately predict the post-yield, large strain behavior over a wide range of temperatures and strain rates.     

CFD simulation of a gas–solid fluidized bed with two vertical jets

A computational fluid dynamics (CFD) model is used to investigate the hydrodynamics of a gas–solid fluidized bed with two vertical jets. Sand particles with a density of 2660 kg/m³ and a diameter of 5.0 × 10^?4 m are employed as the solid phase. Numerical computation is carried out in a 0.57 m × 1.00 m two-dimensional bed using a commercial CFD code, CFX 4.4, together with user-defined Fortran subroutines. The applicability of the CFD model is validated by predicting the bed pressure drop in a bubbling fluidized bed, and the jet detachment time and equivalent bubble diameter in a fluidized bed with a single jet. Subsequently, the model is used to explore the hydrodynamics of two vertical jets in a fluidized bed. The computational results reveal three flow patterns, isolated, merged and transitional jets, depending on the nozzle separation distance and jet gas velocity and influencing significantly the solid circulation pattern. The jet penetration depth is found to increase with increasing jet gas velocity, and can be predicted reasonably well by the correlations of Hong et al. (2003) for isolated jets and of Yang and Keairns (1979) for interacting jets.     

Experimental study on the compressive behaviors of CFRP tube-encased concrete columns

Nine square concrete columns including 6 CFRP/ECCs and 3 concrete columns are prepared, which have cross-section of 200 mm × 200 mm and height of 600 mm. The CFRP tubes with fibers oriented at hoop direction were manufactured to have 3 or 5 layers of CFRP with 10 mm, 20 mm, or 40 mm rounding corner radii at vertical edges. A 100 mm overlap in the direction of fibers was provided to ensure proper bond. Uniaxial compression tests were conducted to investigate the compressive behavior. It is evident that the CFRP tube confinement can improve the behavior of concrete core, in terms of axial compressive strength or axial deformability. Test results show that the stress-strain behavior of CFRP/ECCs vary with different confinement parameters, such as the number of confinement layers and the rounding corner radius.     

Hybrid RANS/LES computations of plane impinging jets with DES and PANS models

The qualities of a DES (Detached Eddy Simulation) and a PANS (Partially-Averaged Navier–Stokes) hybrid RANS/LES model, both based on the k–ω RANS turbulence model of Wilcox (2008, “Formulation of the k–ω turbulence model revisited” AIAA J., 46: 2823–2838), are analysed for simulation of plane impinging jets at a high nozzle-plate distance (H/B = 10, Re = 13,500; H is nozzle-plate distance, B is slot width; Reynolds number based on slot width and maximum velocity at nozzle exit) and a low nozzle-plate distance (H/B = 4, Re = 20,000). The mean velocity field, fluctuating velocity components, Reynolds stresses and skin friction at the impingement plate are compared with experimental data and LES (Large Eddy Simulation) results. The k–ω DES model is a double substitution type, following Davidson and Peng (2003, “Hybrid LES–RANS modelling: a one-equation SGS model combined with a k–ω model for predicting recirculating flows” Int. J. Numer. Meth. Fluids, 43: 1003–1018). This means that the turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the eddy viscosity formula. The k–ω PANS model is derived following Girimaji (2006, “Partially-Averaged Navier–Stokes model for turbulence: a Reynolds-Averaged Navier–Stokes to Direct Numerical Simulation bridging method” J. Appl. Mech., 73: 413–421). The turbulent length scale in the PANS model is constructed from the total turbulent kinetic energy and the sub-filter dissipation rate. Both hybrid models change between RANS (Reynolds-Averaged Navier–Stokes) and LES based on the cube root of the cell volume. The hybrid techniques, in contrast to RANS, are able to reproduce the turbulent flow dynamics in the shear layers of the impacting jet. The change from RANS to LES is much slower however for the PANS model than for the DES model on fine enough grids. This delays the break-up process of the vortices generated in the shear layers with as a consequence that the DES model produces better results than the PANS model.     

Computational homogenization of periodic beam-like structures

This paper is concerned with the computation of the effective elastic properties of periodic beam-like structures. The homogenization theory is used and leads to an equivalent anisotropic Navier–Bernoulli–Saint-Venant beam. The overall behavior is obtained from the solution of basic cell problems posed on the three-dimensional period of the structure and solved using three-dimensional finite element implementation.This procedure is first applied to two corrugated zigzag and sinus beams subjected to in-plane loading. Next, the axial elastic properties of a stranded ‘6 + 1’ wire-cable are computed. The effective properties values obtained appear to be very close to analytical reference results showing the efficiency of the approach.     

Numerical simulation and analysis of the effect of rain and surface property on wind-turbine airfoil performance

The performance of wind turbines is significantly affected by the atmospheric condition of their operating environment. Because rain is a common phenomenon in many parts of the world, understanding its effect on the performance of wind turbines provides valuable information in determining the site for a new wind farm. We developed a multiphase computational fluid dynamics (CFD) model to estimate the effect of rain by simulating the actual physical process of rain droplets forming a water layer over the blades by coupling the Lagrangian Discrete Phase Model (DPM) and the Eulerian Volume of Fluid (VOF) models. We applied our model to a wind-turbine blade airfoil and studied the effect of rain for different rainfall rates in addition to the effect of surface tension and surface property of the airfoil. We observed that, at low rainfall rates, the performance of the airfoil is highly sensitive to the rainfall rate. However, if the rainfall rate is high enough to immerse most of the airfoil surface under water, a further increase in the rainfall rate does not have a substantial effect on the performance of the airfoil.     

Energy loss at bends in the pneumatic conveying of fly ash

An accurate estimation of the total pressure drop of a pipeline is important to the reliable design of a pneumatic conveying system. The present paper presents results from an investigation into the modelling of the pressure drop at a bend in the pneumatic conveying of fly ash. Seven existing bend models were used (in conjunction with solids friction models for horizontal and vertical straight pipes, and initial acceleration losses) to predict the total pipeline pressure drop in conveying fly ash (median particle diameter: 30 μm; particle density: 2300 kg/m³; loose-poured bulk density: 700 kg/m³) in three test rigs (pipelines with dimensions of 69 mm inner diameter (I.D.) × 168 m length; 105 mm I.D. × 168 m length; 69 mm I.D. × 554 m length). A comparison of the pneumatic conveying characteristics (PCC) predicted using the seven bend models and experimental results shows that the predicted total pipeline PCC and trends depend on the choice of bend model. While some models predict trends that agree with the experimental results, other models predicted greater bend pressure drops for the dense phase of fly ash than for the dilute phase. Models of Pan, R. (1992). Improving scale-up procedures for the design of pneumatic conveying systems. Doctoral dissertation, University of Wollongong, Australia, Pan, R., & Wypych, P.W. (1998). Dilute and dense phase pneumatic conveying of fly ash. In Proceedings of the sixth International Conference on Bulk Materials Storage and Transportation (pp. 183–189), Wollongong, NSW, Australia and Chambers, A.J., & Marcus, R.D. (1986). Pneumatic conveying calculations. In Proceedings of the second International Conference on Bulk Materials Storage and Transportation (pp. 49–52), Wollongong, Australia reliably predicted the bend losses for systems conveying fly ash over a large range of air flows.