skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Pseudo-bistability of viscoelastic shells
Viscoelastic shells subjected to a pressure loading exhibit rich and complex time-dependent responses. Here we focus on the phenomenon of pseudo-bistability, i.e. a viscoelastic shell can stay inverted when pressure is removed, and snap to its natural shape after a delay time. We model and explain the mechanism of pseudo-bistability with a viscoelastic shell model. It combines the small strain, moderate rotation shell theory with the standard linear solid as the viscoelastic constitutive law, and is applicable to shells with arbitrary axisymmetric shapes. As a case study, we investigate the pseudo-bistable behaviour of viscoelastic ellipsoidal shells. Using the proposed model, we successfully predict buckling of a viscoelastic ellipsoidal shell into its inverted configuration when subjected to an instantaneous pressure, creeping when the pressure is held, staying inverted after the pressure is removed, and eventually snapping back after a delay time. The stability transition of the shell from a monostable, temporarily bistable and eventually back to the monostable state is captured by examining the evolution of the instantaneous pressure–volume change relation at different time of the holding and releasing process. A systematic parametric study is conducted to investigate the effect of geometry, viscoelastic properties and loading history on the pseudo-bistable behaviour. This article is part of the theme issue 'Probing and dynamics of shock sensitive shells'.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We performed dynamic pressure buckling experiments on defect-seeded spherical shells made of a common silicone elastomer. Unlike in quasi-static experiments, shells buckled at ostensibly subcritical pressures, i.e. below the experimentally determined critical load at which buckling occurs elastically, often following a significant delay period from the time of load application. While emphasizing the close connections to elastic shell buckling, we rely on viscoelasticity to explain our observations. In particular, we demonstrate that the lower critical load may be determined from the material properties, which is rationalized by a simple analogy to elastic spherical shell buckling. We then introduce a model centred on empirical quantities to show that viscoelastic creep deformation lowers the critical load in the same predictable, quantifiable way that a growing defect would in an elastic shell. This allows us to capture how both the deflection at instability and the time delay depend on the applied pressure, material properties and defect geometry. These quantities are straightforward to measure in experiments. Thus, our work not only provides intuition for viscoelastic behaviour from an elastic shell buckling perspective but also offers an accessible pathway to introduce tunable, time-controlled actuation to existing mechanical actuators, e.g. pneumatic grippers. 
    more » « less

    The chemo-mechanical loading of rocks causes the dissolution and precipitation of multiple phases in the rock. This dissolution and precipitation of load-bearing mineral phases lead to the stress redistribution in neighboring phases, which in turn results in deformational changes of the sample composite. The aim of this study is to investigate the link between microstructural evolution and creep behavior of shale rocks subjected to chemo-mechanical loading through modeling time-dependent deformation induced by the dissolution-precipitation process. The model couples the microstructural evolution of the shale rocks with the stress/strain fields inside the material as a function of time. The modeling effort is supplemented with an experimental study where shale rocks were exposed to CO2-rich brine under high temperature and pressure conditions. 3D snapshots of the sample microstructure were generated using segmented micro-CT images of the shale sample. The time-evolving microstructures were then integrated with the Finite element-based mechanical model to simulate the creep induced by dissolution and precipitation processes independent of the intrinsic viscoelasticity/viscoplasticity of the mineral phases. After computation of the time-dependent viscoelastic properties of the shale composite, the combined microstructure model and finite element model were utilized to predict the time-dependent stress and strain fields in different zones of reacted shale.


    Determination of viscous behavior of shale rocks is key in wide range of applications such as stability of reservoirs, stability of geo-structures subjected to environmental forcing, underground storage of hazardous materials and hydraulic fracturing. Short-term creep strains in hydraulic fracturing can change stress fields and in turn can impact the hydraulic fracturing procedures(H. Sone & Zoback, 2010; Hiroki Sone & Zoback, 2013). While long-term creep strains can hamper the reservoir performance due to the reduction in permeability of the reservoir by closing of fractures and fissures(Du, Hu, Meegoda, & Zhang, 2018; Rybacki, Meier, & Dresen, 2016; Sharma, Prakash, & Abedi, 2019; Hiroki Sone & Zoback, 2014). Owing to these significance of creep strain, it is important to understand the viscoelastic/viscoplastic behavior of shales.

    more » « less
  3. We investigate aspects of the spherical squirmer model employing both large-scale numerical simulations and asymptotic methods when the squirmer is placed in weakly elastic fluids. The fluids are modelled by differential equations, including the upper-convected Maxwell (UCM)/Oldroyd-B, finite-extensibility nonlinear elastic model with Peterlin approximation (FENE-P) and Giesekus models. The squirmer model we examine is characterized by two dimensionless parameters related to the fluid velocity at the surface of the micro-swimmer: the slip parameter $\xi $ and the swirl parameter $\zeta $ . We show that, for swimming in UCM/Oldroyd-B fluids, the elastic stress becomes singular at a critical Weissenberg number, Wi , that depends only on $\xi$ . This singularity for the UCM/Oldroyd-B models is independent of the domain exterior to the swimmer, or any other forces considered in the momentum balance for the fluid – we believe that this is the first time such a singularity has been explicitly demonstrated. Moreover, we show that the behaviour of the solution at the poles is purely extensional in character and is the primary reason for the singularity in the Oldroyd-B model. When the Giesekus or the FENE-P models are utilized, the singularity is removed. We also investigate the mechanism behind the speed and rotation rate enhancement associated with the addition of swirl in the swimmer's gait. We demonstrate that, for all models, the speed is enhanced by swirl, but the mechanism of enhancement depends intrinsically on the rheological model employed. Special attention is paid to the propulsive role of the pressure and elucidated upon. We also study the region of convergence of the perturbation solutions in terms of Wi . When techniques that accelerate the convergence of series are applied, transformed solutions are derived that are in very good agreement with the results obtained by simulations. Finally, both the analytical and numerical results clearly indicate that the low- Wi region is more important than one would expect and demonstrate all the major phenomena observed when swimming with swirl in a viscoelastic fluid. 
    more » « less
  4. The bistable fluttering response of heavy inverted flags with different aspect ratios ( $AR$ ) is investigated to determine how the vortical structures affect the intermittent vibration response of the flag. A heavy inverted flag in a uniform flow may exhibit several response modes; amongst them are three major modes that occur over an extended velocity range: stationary, large-scale periodic oscillation and one-sided deflected modes. Significant hysteretic bistability is observed at the transition between these modes for all $AR$ , which is notably different from the conventional flag vibration with a fixed leading edge and free trailing edge where no hysteresis is observed at the lower $AR$ limit ( $AR<1$ ). The difference is associated with the distinct roles of vortices around the flag. Experiments with flags made of spring steel are conducted in a wind tunnel, where the flow speed is steadily increased and later decreased to obtain different oscillatory modes of the heavy inverted flags. The experimental results are used to validate the numerical model of the same problem. It is found that different critical velocities exist for increasing and decreasing flow velocities, and there is a sustained hysteresis for all $AR$ controlled by the initiation threshold and growth of the leading-edge and side-edge vortices. The effect of the vortices in the bistable oscillation regime is quantified by formulating a modal force partitioning approach. It is shown that $AR$ can significantly alter the static and dynamic vortex interaction with the flexible plate, thereby changing the flag's hysteresis behaviour and bistable response. 
    more » « less
  5. Abstract

    Bistable composite laminates have exhibited enormous potential in morphing and energy harvesting followed by a wide range of application in aerospace, power generation and automobile industries. This study presents the fatigue analysis of bistable laminates in terms of stiffness degradation and loss of bistability. Moisture saturation of the specimens are ensured by keeping them in a controlled laboratory environment for an extended period of time. Mass of the specimens have been measured to quantify the moisture saturation. Fatigue tests are performed at 1 Hz frequency, and R = −1 stress ratio which is the ratio of minimum stress to maximum stress. Specimens are tested for 3 million cycles in displacement control. Load-displacement plot from the test is divided into 5 stiffness regions. A piecewise study of each region has demonstrated good agreement with existing analytical model. Stiffness degradation in 4 regions corresponding to 2 stable configurations follows general trend for composites up to the second stage of damage in three stage damage progression model while the remaining region corresponding to unstable configuration is not considered in this analysis. Test results have been reproduced with minor discrepancy at the specified environmental and loading condition, ply configuration, and size of the laminate. Test protocols, results, and damage analysis presented in this study can be utilized to evaluate the fatigue performance of multistable CFRP structures subjected to higher amplitudes and frequencies.

    more » « less