No Arabic abstract
We investigate the buckling and post-buckling properties of a hyperelastic half-space coated by two hyperelastic layers when the composite structure is subjected to a uniaxial compression. In the case of a half-space coated with a {it single} layer, it is known that when the shear modulus $mu_f$ of the layer is larger than the shear modulus $mu_s$ of the half-space, a linear analysis predicts the existence of a critical stretch and wave number, whereas a weakly nonlinear analysis predicts the existence of a threshold value of the modulus ratio $mu_s/mu_fapprox 0.57$ below which the buckling is super-critical and above which the buckling is sub-critical. It is shown in this paper that when another layer is added, a larger variety of behaviour can be observed. For instance, buckling can occur at a preferred wavenumber super-critically even if both layers are softer than the half-space although the top layer would need to be harder than the bottom layer. When the shear modulus of the bottom layer lies in a certain interval, the super-critical to sub-critical transition can happen a number of times as the shear modulus of the top layer is increased gradually. Thus, an extra layer imparts more flexibility in producing wrinkling patterns with desired properties, and our weakly nonlinear analysis provides a road map on the parameter regimes where this can be achieved.
We investigate the buckling under compression of a slender beam with a distributed lateral elastic support, for which there is an associated cost. For a given cost, we study the optimal choice of support to protect against Euler buckling. We show that with only weak lateral support, the optimum distribution is a delta-function at the centre of the beam. When more support is allowed, we find numerically that the optimal distribution undergoes a series of bifurcations. We obtain analytical expressions for the buckling load around the first bifurcation point and corresponding expansions for the optimal position of support. Our theoretical predictions, including the critical exponent of the bifurcation, are confirmed by computer simulations.
This article investigates the large deflection and post-buckling of composite plates by employing the Carrera Unified Formulation (CUF). As a consequence, the geometrically nonlinear governing equations and the relevant incremental equations are derived in terms of fundamental nuclei, which are invariant of the theory approximation order. By using the Lagrange expansion functions across the laminate thickness and the classical finite element (FE) approximation, layer-wise (LW) refined plate models are implemented. The Newton-Raphson linearization scheme with the path-following method based on the arc-length constraint is employed to solve geometrically non-linear composite plate problems. In this study, different composite plates subjected to large deflections/rotations and post-buckling are analyzed, and the corresponding equilibrium curves are compared with the results in the available literature or the traditional FEM-based solutions. The effects of various parameters, such as stacking sequence, number of layers, loading conditions, and edge conditions are demonstrated. The accuracy and reliability of the proposed method for solving the composite plates geometrically nonlinear problems are verified.
As a stiff polymer tumbles in shear flow, it experiences compressive viscous forces that can cause it to buckle and undergo a sequence of morphological transitions with increasing flow strength. We use numerical simulations to uncover the effects of these transitions on the steady shear rheology of a dilute suspension of stiff polymers. Our results agree with classic scalings for Brownian rods in relatively weak flows but depart from them above the buckling threshold. Signatures of elastoviscous buckling include enhanced shear thinning and an increase in the magnitude of normal stress differences. We discuss our findings in the light of past work on rigid Brownian rods and non-Brownian elastic fibers and highlight the subtle role of thermal fluctuations in triggering instabilities.
Recent theoretical advances offer an exact, first-principle theory of jamming criticality in infinite dimension as well as universal scaling relations between critical exponents in all dimensions. For packings of frictionless spheres near the jamming transition, these advances predict that nontrivial power-law exponents characterize the critical distribution of (i) small inter-particle gaps and (ii) weak contact forces, both of which are crucial for mechanical stability. The scaling of the inter-particle gaps is known to be constant in all spatial dimensions $d$ -- including the physically relevant $d=2$ and 3, but the value of the weak force exponent remains the object of debate and confusion. Here, we resolve this ambiguity by numerical simulations. We construct isostatic jammed packings with extremely high accuracy, and introduce a simple criterion to separate the contribution of particles that give rise to localized buckling excitations, i.e., bucklers, from the others. This analysis reveals the remarkable dimensional robustness of mean-field marginality and its associated criticality.
Oriented block copolymers exhibit a buckling instability when submitted to a tensile test perpendicular to the lamellae direction. In this paper we study this behavior using a coarse grained molecular dynamics simulation approach. Coarse grained models of lamellar copolymers with alternate glassy rubbery layers are generated using the Radical Like Polymerization method, and their mechanical response is studied. For large enough systems, uniaxial tensile tests perpendicular to the direction of the lamellae reveal the occurrence of the buckling instability at low strain. The results that emerge from molecular simulation are compared to an elastic theory of the buckling instability introduced by Read and coworkers. At high strain rates, significant differences are observed between elastic theory and simulation results for the buckling strain and the buckling wavelength. We explain this difference by the strain rate dependence of the mechanical response. A simple model that takes into account the influence of the strain rate in the mechanical response is presented to rationalize the results at low and moderate strain rates. At very high strain rates, cavitation takes place in the rubbery phase of the sample and limits the validity of the approach.