Nonlinear Finite-Element Modeling of Batter Piles under Lateral Load

In this paper, a finite-element model is developed in which the nonlinear soil behavior is represented by a hyperbolic relation for static load condition and modified hyperbolic relation, which includes both degradation and gap for a cyclic load condition. Although batter piles are subjected to lateral load, the soil resistance is also governed by axial load, which is incorporated by considering the P-Δ moment and geometric stiffness matrix. By adopting the developed numerical model, static and cyclic load analyses are performed adopting an incremental-iterative procedure where the pile is idealized as beam elements and the soil as elastoplastic spring elements. The proposed numerical model is validated with published laboratory and field pile test results under both static and cyclic load conditions. This paper highlights the importance of the degradation factor and its influence on the soil resistance-displacement (p-y) curve, number of cycles of loading, and cyclic load response.

Modeling wrap faced reinforced soil retaining walls subjected to dynamic excitation

This paper describes some of the physical and numerical model tests of reinforced soil retaining walls subjected to dynamic excitation through uni-axial shaking tests. Models of retaining walls are constructed in a perspex box with geotextile reinforcement using the wrap around technique with dry sand backfill and instrumented with displacement sensors, accelerometers and soil pressure sensors. Numerical modelling of these shaking table tests is carried using FLAC. Numerical model is validated by comparing physical model results. Responses of wrap faced walls with different number of reinforcement layers are discussed from both the physical and numerical model tests. Results showed that the displacements are decreasing with the increase in number of reinforcement layers while acceleration amplifications are not affected significantly.

The Behaviour of Two-Phase Flow of DNAPL and Water through a Fractured Rock under Confining Pressure

This study presents the characterization of DNAPL and water flow in a fracture under confining pressure. A comprehensive mathematical model and the conditions under which DNAPL will enter an initially water-saturated deforming rock fracture are discussed. A numerical model with which to predict the quantity of each phase in terms of their saturations in deforming rock joint is developed. The effect of varying confining stresses on the traverse time of DNAPL across a fractured aquitard is studied. The sensitivity analysis for physical and hydraulic properties like initial fracture apertures, fracture dips, equivalent fracture aperture and confining pressures are performed and discussed.

Effect of pouring temperature on cooling slope casting of semi-solid Al-Si-Mg alloy

Present trend of semi-solid processing is directed towards rheocasting route which allows manufacturing of near-net-shape cast components directly from the prepared semi-solid slurry. Generation of globular equi-axed grains during solidification of rheocast components, compared to the columnar dendritic structure of conventional casting routes, facilitates the manufacturing of components with improved mechanical properties and structural integrity. In the present investigation, a cooling slope has been designed and indigenously fabricated to produce semi solid slurry of Al-Si-Mg (A356) alloy and successively cast in a metallic mould. The scope of the present work discusses about development of a numerical model to simulate the liquid metal flow through cooling slope using Eulerian two-phase flow approach and to investigate the effect of pouring temperature on cooling slope semi-solid slurry generation process. The two phases considered in the present model are liquid metal and air. Solid fraction evolution of the solidifying melt is tracked at different locations of the cooling slope, following Schiel's equation. The continuity equation, momentum equation and energy equation are solved considering thin wall boundary condition approach. During solidification of the liquid metal, a modified temperature recovery scheme has been employed taking care of the latent heat release and change of fraction of liquid. The results obtained from simulations are compared with experimental findings and good agreement has been found.

Eulerian two-phase flow simulation and experimental validation of semisolid slurry generation process using cooling slope

Experimental and numerical studies of slurry generation using a cooling slope are presented in the paper. The slope having stainless steel body has been designed and constructed to produce semisolid A356 Al alloy slurry. The pouring temperature of molten metal, slope angle of the cooling slope and slope wall temperature were varied during the experiment. A multiphase numerical model, considering liquid metal and air, has been developed to simulate the liquid metal flow along the cooling channel using an Eulerian two-phase flow approach. Solid fraction evolution of the solidifying melt is tracked at different locations of the cooling channel following Schiel's equation. The continuity, momentum and energy equations are solved considering thin wall boundary condition approach. During solidification of the melt, based on the liquid fraction and latent heat of the alloy, temperature of the alloy is modified continuously by introducing a modified temperature recovery method. Numerical simulations has been carried out for semisolid slurry formation by varying the process parameters such as angle of the cooling slope, cooling slope wall temperature and melt superheat temperature, to understand the effect of process variables on cooling slope semisolid slurry generation process such as temperature distribution, velocity distribution and solid fraction of the solidifying melt. Experimental validation performed for some chosen cases reveals good agreement with the numerical simulations.

Friction Welding of Thixocast A356 Aluminium Alloy

In this paper, a numerical model for friction welding of thixo-cast materials is developed, which includes a coupling of thermal effect and plastic deformation using a finite element method (FEM). As the constitutive equations for flow behavior of materials for a thixo-cast material are expected to be different from those of conventionally cast material of the same alloy, the necessary material data are experimentally determined from isothermal hot compression tests of the A356 thixocast alloy. The Johnson-Cook model has been employed to represent the flow behavior of the thixocast A356 alloy. The purpose of this FEM analysis is to provide better understanding of the friction welding process of thixo-cast material, and to obtain optimized process parameters before an actual welding is carried out.

Modeling the dynamic response of wrap-faced reinforced soil retaining walls

This paper describes the development of a numerical model for simulating the shaking table tests on wrap-faced reinforced soil retaining walls. Some of the physical model tests carried out on reinforced soil retaining walls subjected to dynamic excitation through uniaxial shaking tests are briefly discussed. Models of retaining walls are constructed in a perspex box with geotextile reinforcement using the wraparound technique with dry sand backfill and instrumented with displacement sensors, accelerometers, and soil pressure sensors. Results showed that the displacements decrease with the increase in number of reinforcement layers, whereas acceleration amplifications were not affected significantly. Numerical modeling of these shaking table tests is carried out using the Fast Lagrangian Analysis of Continua program. The numerical model is validated by comparing the results with experiments on physical models. Responses of wrap-faced walls with varying numbers of reinforcement layers are compared. Sensitivity analysis performed on the numerical models showed that the friction and dilation angle of backfill material and stiffness properties of the geotextile-soil interface are the most affecting parameters for the model response.

Study of dendrite growth in a rotational flow field

A numerical model to study the growth of dendrites in a pure metal solidification process with an imposed rotational flow field is presented. The micro-scale features of the solidification are modeled by the well-known enthalpy technique. The effect of flow changing the position of the dendrite is captured by the Volume of Fluid (VOF) method. An imposed rigid-body rotational flow is found to gradually transform the dendrite into a globular microstructure. A parametric study is carried out for various angular velocities and the time for merger of dendrite arms is compared with the order estimate obtained from scaling.

Optical read-out scheme based on Grated Waveguide Cantilever cavity resonance for interrogation of cantilever sensor arrays

The design and analysis of an optical read-out scheme based on a grated waveguide (GWG) resonator for interrogating microcantilever sensor arrays is presented. The optical system consisting of a micro cantilever monolithically integrated in proximity to a grated waveguide (GWG), is realized in silicon optical bench platform. The mathematical analysis of the optical system is performed using a Fabry-Perot interferometer model with a lossy cavity formed between the cantilever and the GWG and an analytical expression is derived for the optical power transmission as a function of the cantilever deflection which corresponds to cavity width variation. The intensity transmission of the optical system for different cantilever deflections estimated using the analytical expression captures the essential features exhibited by a FDTD numerical model.

Velocity Fluctuations in Helical Propulsion: How Small Can a Propeller Be

Helical propulsion is at the heart of locomotion strategies utilized by various natural and artificial swimmers. We used experimental observations and a numerical model to study the various fluctuation mechanisms that determine the performance of an externally driven helical propeller as the size of the helix is reduced. From causality analysis, an overwhelming effect of orientational noise at low length scales is observed, which strongly affects the average velocity and direction of motion of a propeller. For length scales smaller than a few micrometers in aqueous media, the operational frequency for the propulsion system would have to increase as the inverse cube of the size, which can be the limiting factor for a helical propeller to achieve locomotion in the desired direction.

Structural health monitoring of aerospace structural components using wave propagation based diagnostics

The paper discusses a wave propagation based method for identifying the damages in an aircraft built up structural component such as delamination and skin-stiffener debonding. First, a spectral finite element mode l (SFEM) is developed for modeling wave propagation in general built-up structures by using the concept of assembling 2D spectral plate elements. The developed numerical model is validated using conventional 2-D FEM. Studies are performed to capture the mode coupling,that is, the flexural-axial coupling present in the wave responses. Lastly, the damages in these built up structures are then identified using the developed SFEM model and the measured responses using the concept Damage Force Indicator (DFI) technique.

Transmission of elastic waves through a frictional boundary

Incident energy gets transmitted, reflected and absorbed across an interface in jointed rock mass leading to energy dissipation and alteration of waves. Wave velocities get attenuated during their propagation across joints and this behavior is studied using bender/extender element tests. The velocity attenuation and modulus reduction observed in experimental tests are modeled with three dimensional distinct element code and results are validated. Normal propagation of an incident shear wave through a jointed rock mass cause slip of the rock blocks if shear stress of wave exceeds the shear strength of the joint. As the properties of joint determine the transmission of energy across an interface, a parametric study is then conducted with the validated numerical model by varying the parameters that may determine the energy transmission across a joint using modified Miller's method. Results of the parametric study are analyzed and presented in the paper. (C) 2014 Elsevier Ltd. All rights reserved.

Modeling of the non-isothermal liquid droplet impact on a heated solid substrate with heterogeneous wettability

A comprehensive numerical investigation on the impingement and spreading of a non-isothermal liquid droplet on a solid substrate with heterogeneous wettability is presented in this work. The time-dependent incompressible Navier-Stokes equations are used to describe the fluid flow in the liquid droplet, whereas the heat transfer in the moving droplet and in the solid substrate is described by the energy equation. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the time-dependent incompressible Navier-Stokes equation and the energy equation in the time-dependent moving domain. Moreover, the Marangoni convection is included in the variational form of the Navier-Stokes equations without calculating the partial derivatives of the temperature on the free surface. The heterogeneous wettability is incorporated into the numerical model by defining a space-dependent contact angle. An array of simulations for droplet impingement on a heated solid substrate with circular patterned heterogeneous wettability are presented. The numerical study includes the influence of wettability contrast, pattern diameter, Reynolds number and Weber number on the confinement of the spreading droplet within the inner region, which is more wettable than the outer region. Also, the influence of these parameters on the total heat transfer from the solid substrate to the liquid droplet is examined. We observe that the equilibrium position depends on the wettability contrast and the diameter of the inner surface. Consequently. the heat transfer is more when the wettability contrast is small and/or the diameter of inner region is large. The influence of the Weber number on the total heat transfer is more compared to the Reynolds number, and the total heat transfer increases when the Weber number increases. (C) 2015 Elsevier Ltd. All rights reserved.

Three-dimensional warping in strip-based composite four-bar mechanisms

This work intends to demonstrate the effect of geometrically non-linear cross-sectional analysis of certain composite beam-based four-bar mechanisms in predicting the three-dimensional warping of the cross-section. The only restriction in the present analysis is that the strains within each elastic body remain small (i.e., this work does not deal with materials exhibiting non-linear constitutive laws at the 3-D level). Here, all component bars of the mechanism are made of fiber-reinforced laminates. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. Each component of the mechanism is modeled as a beam based on geometrically non-linear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and non-linear 1-D analyses along the three beam reference curves. The splitting of the three-dimensional beam problem into two- and one-dimensional parts, called dimensional reduction, results in a tremendous savings of computational effort relative to the cost of three-dimensional finite element analysis, the only alternative for realistic beams. The analysis of beam-like structures made of laminated composite materials requires a much more complicated methodology. Hence, the analysis procedure based on Variational Asymptotic Method (VAM), a tool to carry out the dimensional reduction, is used here. The representative cross-sections of all component bars are analyzed using two different approaches: (1) Numerical Model and (2) Analytical Model. Four-bar mechanisms are analyzed using the above two approaches for Omega = 20 rad/s and Omega = pi rad/s and observed the same behavior in both cases. The noticeable snap-shots of the deformation shapes of the mechanism about 1000 frames are also reported using commercial software (I-DEAS + NASTRAN + ADAMS). The maximum out-of-plane warping of the cross-section is observed at the mid-span of bar-1, bar-2 and bar-3 are 1.5 mm, 250 mm and 1.0 mm, respectively, for t = 0:5 s. (C) 2015 Elsevier Ltd. All rights reserved.

Effect of applied pressure on densification of monolithic ZrCx ceramic by reactive hot pressing

The effect of applied pressure on reactive hot pressing (RHP) of zirconium (Zr):graphite (C) in molar ratios of 1:0.5, 1:0.67, 1:0.8, and 1:1 was studied at 1200 degrees C for 60 min. The relative density achievable increased with increasing pressure and ranged from 99% at 4 MPa for ZrC0.5 to 93% for stoichiometric ZrC at 100 MPa. The diminishing influence of pressure on the final density with increasing stoichiometry is attributed to two causes: the decreasing initial volume fraction of the plastically deforming Zr metal which leads to the earlier formation of a contiguous, stress shielding carbide skeleton and the larger molar volume shrinkage during reaction which leads to pore formation in the final stages. A numerical model of the creep densification of a dynamically evolving microstructure predicts densities that are consistent with observations and confirm that the availability of a soft metal is primarily responsible for the achievement of such elevated densification during RHP. The ability to densify nonstoichiometric compositions like ZrC0.5 at pressures as low as 4 MPa offers an alternate route to fabricating dense nonstoichiometric carbides.