Wrinkled and slack membranes: nonlinear 3D elasticity solutions via smooth DMS-FEM and experiment

The smooth DMS-FEM, recently proposed by the authors, is extended and applied to the geometrically nonlinear and ill-posed problem of a deformed and wrinkled/slack membrane. A key feature of this work is that three-dimensional nonlinear elasticity equations corresponding to linear momentum balance, without any dimensional reduction and the associated approximations, directly serve as the membrane governing equations. Domain discretization is performed with triangular prism elements and the higher order (C1 or more) interelement continuity of the shape functions ensures that the errors arising from possible jumps in the first derivatives of the conventional C0 shape functions do not propagate because the ill-conditioned tangent stiffness matrices are iteratively inverted. The present scheme employs no regularization and exhibits little sensitivity to h-refinement. Although the numerically computed deformed membrane profiles do show some sensitivity to initial imperfections (nonplanarity) in the membrane profile needed to initiate transverse deformations, the overall patterns of the wrinkles and the deformed shapes appear to be less so. Finally, the deformed profiles, computed through the DMS FEM-based weak formulation, are compared with those obtained through an experiment on an ultrathin Kapton membrane, wherein wrinkles form because of the applied boundary displacement conditions. Comparisons with a reported experiment on a rectangular membrane are also provided. These exercises lend credence to the feasibility of the DMS FEM-based numerical route to computing post-wrinkled membrane shapes. Copyright (c) 2012 John Wiley & Sons, Ltd.

Elasto-plastic analysis of jointed rocks using discrete continuum and equivalent continuum approaches

Results from elasto-plastic numerical simulations of jointed rocks using both the equivalent continuum and discrete continuum approaches are presented, and are compared with experimental measurements. Initially triaxial compression tests on different types of rocks with wide variation in the uniaxial compressive strength are simulated using both the approaches and the results are compared. The applicability and relative merits and limitations of both the approaches for the simulation of jointed rocks are discussed. It is observed that both the approaches are reasonably good in predicting the real response. However, the equivalent continuum approach has predicted somewhat higher stiffness values at low strains. Considering the modelling effort involved in case of discrete continuum approach, for problems with complex geometry, it is suggested that a proper equivalent continuum model can be used, without compromising much on the accuracy of the results. Then the numerical analysis of a tunnel in Japan is taken up using the continuum approach. The deformations predicted are compared well against the field measurements and the predictions from discontinuum analysis. (C) 2012 Elsevier Ltd. All rights reserved.

Geometrically non-linear dynamics of composite four-bar mechanisms

This work intends to demonstrate the importance of a geometrically nonlinear cross-sectional analysis of certain composite beam-based four-bar mechanisms in predicting system dynamic characteristics. All component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). 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. For the thin rectangular cross-sections considered here, the 2-D cross-sectional non-linearity is also overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the non-linear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the non-linear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict multi-body dynamic responses more quickly and accurately than would otherwise be possible. The analysis methodology can be viewed as a three-step procedure: First, the cross-sectional properties of each bar of the mechanism is determined analytically based on an asymptotic procedure, starting from Classical Laminated Shell Theory (CLST) and taking advantage of its thin strip geometry. Second, the dynamic response of the non-linear, flexible four-bar mechanism is simulated by treating each bar as a 1-D beam, discretized using finite elements, and employing energy-preserving and -decaying time integration schemes for unconditional stability. Finally, local 3-D deformations and stresses in the entire system are recovered, based on the 1-D responses predicted in the previous step. With the model, tools and procedure in place, we identify and investigate a few four-bar mechanism problems where the cross-sectional non-linearities are significant in predicting better and critical system dynamic characteristics. This is carried out by varying stacking sequences (i.e. the arrangement of ply orientations within a laminate) and material properties, and speculating on the dominating diagonal and coupling terms in the closed-form non-linear beam stiffness matrix. A numerical example is presented which illustrates the importance of 2-D cross-sectional non-linearities and the behavior of the system is also observed by using commercial software (I-DEAS + NASTRAN + ADAMS). (C) 2012 Elsevier Ltd. All rights reserved.

VAM applied to dimensional reduction of non-linear hyperelastic plates

This work aims at dimensional reduction of non-linear isotropic hyperelastic plates in an asymptotically accurate manner. The problem is both geometrically and materially non-linear. The geometric non-linearity is handled by allowing for finite deformations and generalized warping while the material non-linearity is incorporated through hyperelastic material model. The development, based on the Variational Asymptotic Method (VAM) with moderate strains and very small thickness to shortest wavelength of the deformation along the plate reference surface as small parameters, begins with three-dimensional (3-D) non-linear elasticity and mathematically splits the analysis into a one-dimensional (1-D) through-the-thickness analysis and a two-dimensional (2-D) plate analysis. Major contributions of this paper are derivation of closed-form analytical expressions for warping functions and stiffness coefficients and a set of recovery relations to express approximately the 3-D displacement, strain and stress fields. Consistent with the 2-D non-linear constitutive laws, 2-D plate theory and corresponding finite element program have been developed. Validation of present theory is carried out with a standard test case and the results match well. Distributions of 3-D results are provided for another test case. (c) 2012 Elsevier Ltd. All rights reserved.

Poincare invariant quantum field theories with twisted internal symmetries

Following up the work of 1] on deformed algebras, we present a class of Poincare invariant quantum field theories with particles having deformed internal symmetries. The twisted quantum fields discussed in this work satisfy commutation relations different from the usual bosonic/fermionic commutation relations. Such twisted fields by construction are nonlocal in nature. Despite this nonlocality we show that it is possible to construct interaction Hamiltonians which satisfy cluster decomposition principle and are Lorentz invariant. We further illustrate these ideas by considering global SU(N) symmetries. Specifically we show that twisted internal symmetries can provide a natural-framework for the discussion of the marginal deformations (beta-deformations) of the N = 4 SUSY theories.

Biomechanics of Substrate Boring by Fig Wasps: Role of Zinc in Insect Cuticle

There are many biomechanical challenges that a female insect must meet to successfully oviposit and ensure her evolutionary success. These begin with selection of a suitable substrate through which the ovipositor must penetrate without itself buckling or fracturing. The second phase corresponds to steering and manipulating the ovipositor to deliver eggs at desired locations. Finally, the insect must retract her ovipositor fast to avoid possible predation and repeat this process multiple times during her lifetime. From a materials perspective, insect oviposition is a fascinating problem and poses many questions. Specifically, are there diverse mechanisms that insects use to drill through hard substrates without itself buckling or fracturing? What are the structure-property relationships in the ovipositor material? These are some of the questions we address with a model system consisting of a parasitoid fig wasp - fig substrate system. To characterize the structure of ovipositors, we use scanning electron microscopy with a detector to quantify the presence of transition elements. Our results show that parasitoid ovipositors have teeth like structures on their tips and contain high amounts of zinc as compared to remote regions. Sensillae are present along the ovipositor to aid detection of chemical species and mechanical deformations. To quantify the material properties of parasitoid ovipositors, we use an atomic force microscope and show that tip regions have higher modulus as compared to remote regions. Finally, we use videography to show that ovipositors buckle during oviposition and estimate the forces needed to cause substrate boring based on Euler buckling analysis. Such methods may be useful for the design of functionally graded surgical tools.

Effects of different modes of hot cross-rolling in 7010 aluminum alloy: part I. evolution of microstructure and texture

The current study describes the evolution of microstructure and texture in an Al-Zn-Mg-Cu-Zr-based 7010 aluminum alloy during different modes of hot cross-rolling. Processing of materials involves three different types of cross-rolling. The development of texture in the one-step cross-rolled specimen can be described by a typical beta-fiber having the maximum intensity near Copper (Cu) component. However, for the multi-step cross-rolled specimens, the as-rolled texture is mainly characterized by a strong rotated-Brass (Bs) component and a very weak rotated-cube component. Subsequent heat treatment leads to sharpening of the major texture component (i.e., rotated-Bs). Furthermore, the main texture components in all the specimens appear to be significantly rotated in a complex manner away from their ideal positions because of non-symmetric deformations in the two rolling directions. Detailed microstructural study indicates that dynamic recovery is the dominant restoration mechanism operating during the hot rolling. During subsequent heat treatment, static recovery dominates, while a combination of particle-stimulated nucleation (PSN) and strain-induced grain boundary migration (SIBM) causes partial recrystallization of the grain structure. The aforementioned restoration mechanisms play an important role in the development of texture components. The textural development in the current study could be attributed to the combined effects of (a) cross-rolling and inter-pass annealing that reduce the intensity of Cu component after each successive pass, (b) recrystallization resistance of Bs-oriented grains, (c) stability of Bs texture under cross-rolling, and (d) Zener pinning by Al3Zr dispersoids.

Room-temperature anelasticity and viscoplasticity of Cu-Zr bulk metallic glasses evaluated using nanoindentation

Anelastic and viscoplastic characteristics of Cu50Zr50 and Cu65Zr35 binary bulk metallic glasses at room temperature were examined through nanoindentation creep experiments. Results show that both the deformations are relatively more pronounced in Cu50Zr50 than in Cu65Zr35, and their amount increases with the loading rate. The results are analyzed in terms of the influences of structural defects and loading rate on the room temperature indentation creep.

Dynamics of a flexible splitter plate in the wake of a circular cylinder

Rigid splitter plates in the wake of bluff bodies are known to suppress the primary vortex shedding. In the present work, we experimentally study the problem of a flexible splitter plate in the wake of a circular cylinder. In this case, the splitter plate is free to continuously deform along its length due to the fluid forces acting on it; the flexural rigidity (EI) of the plate being an important parameter. Direct visualizations of the splitter plate motions, for very low values of flexural rigidity (EI), indicate periodic traveling wave type deformations of the splitter plate with maximum tip amplitudes of the order of I cylinder diameter. As the Reynolds number based on cylinder diameter is varied, two regimes of periodic splitter plate motions are found that are referred to as mode I and mode II, with a regime of aperiodic motions between them. The frequency of plate motions in both periodic modes is found to be close to the plane cylinder Strouhal number of about 0.2, while the average frequencies in the non-periodic regime are substantially lower. The measured normalized phase speed of the traveling wave for both periodic modes is also close to the convection speed of vortices in the plane cylinder wake. As the flexural rigidity of the plate (EI) is increased, the response of the plate was found to shift to the right when plotted with flow speed or Re. To better capture the effect of varying EI, we define and use a non-dimensional bending stiffness, K*, similar to the ones used in the flag flutter problem, K*=EI/(0.5 rho(UL3)-L-2), where U is the free-stream velocity and L is the splitter plate length. Amplitude data for different EI cases when plotted against this parameter appear to collapse on to a single curve for a given splitter plate length. Measurements of the splitter plate motions for varying splitter plate lengths indicate that plates that are substantially larger than the formation length of the plane cylinder wake have similar responses, while shorter plates show significant differences.

High Speed Imaging of Generation and Collapse of Multielectron Bubbles in Liquid Helium

We present experimental results on the generation and collapse of multielectron bubbles in liquid helium. By applying voltage pulses to a tungsten tip above the surface of the liquid, millimetre sized deformations were formed. Using high speed photography, we have imaged the disintegration of these deformations into bubbles of sizes ranging from ten to few hundred microns. At temperatures less than 2 K, the bubbles split into smaller bubbles and then disappeared in a time scale of few milliseconds. Smaller bubbles were formed at temperatures around 3 K, but were visible for more than hundreds of milliseconds. Although we have not been able to measure their charge directly, some of these bubbles responded to electric fields, implying these were indeed multielectron bubbles. With the existing theoretical picture, it is not possible to understand the strong dependence of the lifetime of multielectron bubbles on temperature.

A reexamination of some puzzling results in linearized elasticity

In this paper, we analyse three commonly discussed `flaws' of linearized elasticity theory and attempt to resolve them. The first `flaw' concerns cylindrically orthotropic material models. Since the work of Lekhnitskii (1968), there has been a growing body of work that continues to this day, that shows that infinite stresses arise with the use of a cylindrically orthotropic material model even in the case of linearized elasticity. Besides infinite stresses, interpenetration of matter is also shown to occur. These infinite stresses and interpenetration occur when the ratio of the circumferential Young modulus to the radial Young modulus is less than one. If the ratio is greater than one, then the stresses at the center of a spinning disk are found to be zero (recall that for an isotropic material model, the stresses are maximum at the center). Thus, the stresses go abruptly from a maximum value to a value of zero as the ratio is increased to a value even slightly above one! One of the explanations provided for this extremely anomalous behaviour is the failure of linearized elasticity to satisfy material frame-indifference. However, if this is the true cause, then the anomalous behaviour should also occur with the use of an isotropic material model, where, no such anomalies are observed. We show that the real cause of the problem is elsewhere and also show how these anomalies can be resolved. We also discuss how the formulation of linearized elastodynamics in the case of small deformations superposed on a rigid motion can be given in a succinct manner. Finally, we show how the long-standing problem of devising three compatibility relations instead of six can be resolved.

Biomechanics of substrate boring by fig wasps

Female insects of diverse orders bore into substrates to deposit their eggs. Such insects must overcome several biomechanical challenges to successfully oviposit, which include the selection of suitable substrates through which the ovipositor can penetrate without itself fracturing. In many cases, the insect may also need to steer and manipulate the ovipositor within the substrate to deliver eggs at desired locations before rapidly retracting her ovipositor to avoid predation. In the case of female parasitoid ichneumonid wasps, this process is repeated multiple times during her lifetime, thus testing the ability of the ovipositioning apparatus to endure fracture and fatigue. What specific adaptations does the ovipositioning apparatus of a female ichneumonoid wasp possess to withstand these challenges? We addressed this question using a model system composed of parasitoid and pollinator fig wasps. First, we show that parasitoid ovipositor tips have teeth-like structures, preferentially enriched with zinc, unlike the smooth morphology of pollinator ovipositors. We describe sensillae present on the parasitoid ovipositor tip that are likely to aid in the detection of chemical species and mechanical deformations and sample microenvironments within the substrate. Second, using atomic force microscopy, we show that parasitoid tip regions have a higher modulus compared with regions proximal to the abdomen in parasitoid and pollinator ovipositors. Finally, we use videography to film wasps during substrate boring and analyse buckling of the ovipositor to estimate the forces required for substrate boring. Together, these results allow us to describe the biomechanical principles underlying substrate boring in parasitoid ichneumonid wasps. Such studies may be useful for the biomimetic design of surgical tools and in the use of novel mechanisms to bore through hard substrates.

STM verification of the reduction of the Young's modulus of CdS nanoparticles at smaller sizes

We demonstrate the first STM evaluation of the Young's modulus (E) of nanoparticles (NPs) of different sizes. The sample deformation induced by tip-sample interaction has been determined using current-distance (I-Z) spectroscopy. As a result of tip-sample interaction, and the induced surface deformations, the I-z curves deviates from pure exponential dependence. Normally, in order to analyze the deformation quantitatively, the tip radius must be known. We show, that this necessity is eliminated by measuring the deformation on a substrate with a known Young's modulus (Au(111)) and estimating the tip radius, and afterwards, using the same tip (with a known radius) to measure the (unknown) Young's modulus of another sample (nanoparticles of CdS). The Young's modulus values found for 3 NP's samples of average diameters of 3.7, 6 and 7.5 nm, were E similar to 73%, 78% and 88% of the bulk value, respectively. These results are in a good agreement with the theoretically predicted reduction of the Young's modulus due to the changes in hydrostatic stresses which resulted from surface tension in nanoparticles with different sizes. Our calculation using third order elastic constants gives a reduction of E which scales linearly with 1/r (r is the NP's radius). This demonstrates the applicability of scanning tunneling spectroscopy for local mechanical characterization of nanoobjects. The method does not include a direct measurement of the tip-sample force but is rather based on the study of the relative elastic response. (C) 2014 Elsevier B.V. All rights reserved.

Dynamics of intermittent force fluctuations in vesicular nanotubulation

Irregular force fluctuations are seen in most nanotubulation experiments. The dynamics behind their presence has, however, been neither commented upon nor modeled. A simple estimate of the mean energy dissipated in force drops turns out to be several times the thermal energy. This coupled with the rate dependent nature of the deformation reported in several experiments point to a dynamical origin of the serrations. We simplify the whole process of tether formation through a three-stage model of successive deformations of sphere to ellipsoid, neck-formation, and tubule birth and extension. Based on this, we envisage a rate-softening frictional force at the neck that must be overcome before a nanotube can be pulled out. Our minimal model includes elastic and visco-elastic deformation of the vesicle, and has built-in dependence on pull velocity, vesicle radius, and other material parameters, enabling us to capture various kinds of serrated force-extension curves for different parameter choices. Serrations are predicted in the nanotubulation region. Other features of force-extension plots reported in the literature such as a plateauing serrated region beyond a force drop, serrated flow region with a small positive slope, an increase in the elastic threshold with pull velocity, force-extension curves for vesicles with larger radius lying lower than those for smaller radius, are all also predicted by the model. A toy model is introduced to demonstrate that the role of the friction law is limited to inducing stick-slip oscillations in the force, and all other qualitative and quantitative features emerging from the model can only be attributed to other physical mechanisms included in the deformation dynamics of the vesicle. (C) 2014 AIP Publishing LLC.

Dielectric elastomers: Asymptotically-correct three-dimensional displacement field

Asymptotically-accurate dimensional reduction from three to two dimensions and recovery of 3-D displacement field of non-prestretched dielectric hyperelastic membranes are carried out using the Variational Asymptotic Method (VAM) with moderate strains and very small ratio of the membrane thickness to its shortest wavelength of the deformation along the plate reference surface chosen as the small parameters for asymptotic expansion. Present work incorporates large deformations (displacements and rotations), material nonlinearity (hyperelasticity), and electrical effects. It begins with 3-D nonlinear electroelastic energy and mathematically splits the analysis into a one-dimensional (1-D) through-the-thickness analysis and a 2-D nonlinear plate analysis. Major contribution of this paper is a comprehensive nonlinear through-the-thickness analysis which provides a 2-D energy asymptotically equivalent of the 3-D energy, a 2-D constitutive relation between the 2-D generalized strain and stress tensors for the plate analysis and a set of recovery relations to express the 3-D displacement field. Analytical expressions are derived for warping functions and stiffness coefficients. This is the first attempt to integrate an analytical work on asymptotically-accurate nonlinear electro-elastic constitutive relation for compressible dielectric hyperelastic model with a generalized finite element analysis of plates to provide 3-D displacement fields using VAM. A unified software package `VAMNLM' (Variational Asymptotic Method applied to Non-Linear Material models) was developed to carry out 1-D non-linear analysis (analytical), 2-D non-linear finite element analysis and 3-D recovery analysis. The applicability of the current theory is demonstrated through an actuation test case, for which distribution of 3-D displacements are provided. (C) 2014 Elsevier Ltd. All rights reserved.