A novel method for studying the buckling of nanotubes considering geometrical imperfections

Buckling of nanotubes has been studied using many methods such as molecular dynamics (MD), molecular mechanics, and continuum-based shell theories. In MD, motion of the individual atoms is tracked under applied temperature and pressure, ensuring a reliable estimate of the material response. The response thus simulated varies for individual nanotubes and is only as accurate as the force field used to model the atomic interactions. On the other hand, there exists a rich literature on the understanding of continuum mechanics-based shell theories. Based on the observations on the behavior of nanotubes, there have been a number of shell theory-based approaches to study the buckling of nanotubes. Although some of these methods yield a reasonable estimate of the buckling stress, investigation and comparison of buckled mode shapes obtained from continuum analysis and MD are sparse. Previous studies show that the direct application of shell theories to study nanotube buckling often leads to erroneous results. The present study reveals that a major source of this error can be attributed to the departure of the shape of the nanotube from a perfect cylindrical shell. Analogous to the shell buckling in the macro-scale, in this work, the nanotube is modeled as a thin-shell with initial imperfection. Then, a nonlinear buckling analysis is carried out using the Riks method. It is observed that this proposed approach yields significantly improved estimate of the buckling stress and mode shapes. It is also shown that the present method can account for the variation of buckling stress as a function of the temperature considered. Hence, this can prove to be a robust method for a continuum analysis of nanosystems taking in the effect of variation of temperature as well.

Phenomenology and control of buckling dynamics in multicomponent colloidal droplets

Self-assembly of nano sized particles during natural drying causes agglomeration and shell formation at the surface of micron sized droplets. The shell undergoes sol-gel transition leading to buckling at the weakest point on the surface and produces different types of structures. Manipulation of the buckling rate with inclusion of surfactant (sodium dodecyl sulphate, SDS) and salt (anilinium hydrochloride, AHC) to the nano-sized particle dispersion (nanosilica) is reported here in an acoustically levitated single droplet. Buckling in levitated droplets is a cumulative, complicated function of acoustic streaming, chemistry, agglomeration rate, porosity, radius of curvature, and elastic energy of shell. We put forward our hypothesis on how buckling occurs and can be suppressed during natural drying of the droplets. Global precipitation of aggregates due to slow drying of surfactant-added droplets (no added salts) enhances the rigidity of the shell formed and hence reduces the buckling probability of the shell. On the contrary, adsorption of SDS aggregates on salt ions facilitates the buckling phenomenon with an addition of minute concentration of the aniline salt to the dispersion. Variation in the concentration of the added particles (SDS/AHC) also leads to starkly different morphologies and transient behaviour of buckling (buckling modes like paraboloid, ellipsoid, and buckling rates). Tuning of the buckling rate causes a transition in the final morphology from ring and bowl shapes to cocoon type of structure. (C) 2015 AIP Publishing LLC.

Universal buckling kinetics in drying nanoparticle-laden droplets on a hydrophobic substrate

We provide a comprehensive physical description of the vaporization, self-assembly, agglomeration, and buckling kinetics of sessile nanofluid droplets pinned on a hydrophobic substrate. We have deciphered five distinct regimes of the droplet life cycle. Regimes I-III consists of evaporation-induced preferential agglomeration that leads to the formation of a unique dome-shaped inhomogeneous shell with a stratified varying-density liquid core. Regime IV involves capillary-pressure-initiated shell buckling and stress-induced shell rupture. Regime V marks rupture-induced cavity inception and growth. We demonstrate through scaling arguments that the growth of the cavity (which controls the final morphology or structure) can be described by a universal function.

Modulation of Buckling Dynamics in Nanoparticle Laden Droplets Using External Heating

Dynamics of contact free (levitated) drying of nanofluid droplets is ubiquitous in many application domains ranging from spray drying to pharmaceutics. Controlling the final morphology (macro to micro scales) of the dried out sample poses some serious challenges. Evaporation of solvent and agglomeration of particles leads to porous shell formation in acoustically levitated nanosilica droplets. The capillary pressure due to evaporation across the menisci at the nanoscale pores causes buckling of the shell which leads to ring and bowl shaped final structures. Acoustics plays a crucial role in flattening of droplets which is a prerequisite for initiation of buckling in the shell: Introduction of mixed nanocolloids (sodium dodecyl sulfate + nanosilica) reduces evaporation rate, disrupts formation of porous shell, and enhances mechanical strength of the shell, all of which restricts the process of buckling. Although buckling is completely arrested in such surfactant added droplets, controlled external heating using laser enhances evaporation through the pores in the shell due to thermally induced structural changes and rearrangement of SDS aggregates which reinitializes buckling in such droplets, Furthermore, inclusion of anilinium hydrochloride into the nanoparticle laden droplets produces ions which adsorb and modify the morphology of sodium dodecyl sulfate crystals and reinitializes buckling in the shell (irrespective of external heating conditions). The kinetics of buckling is determined by the combined effect of morphology of the colloidal particles, particle/aggregate diffusion rate within the droplet, and the rate of evaporation of water. The buckling dynamics leads to cavity formation which grows subsequently to yield final structures with drastically different morphological features. The cavity growth is controlled by evaporation through the nanoscate pores and exhibits a universal trend irrespective of heating rate and nanoparticle type.

Morphological transitions and buckling characteristics in a nanoparticle-laden sessile droplet resting on a heated hydrophobic substrate

In this work, we have established the evaporation-liquid flow coupling mechanism by which sessile nanofluid droplets on a hydrophobic substrate evaporate and agglomerate to form unique morphological features under controlled external heating. It is well understood that evaporation coupled with internal liquid flow controls particle transport in a spatiotemporal sense. Flow characteristics inside the heated droplet are investigated and found to be driven by the buoyancy effects. Velocity magnitudes are observed to increase by an order at higher temperatures with similar looking flow profiles. The recirculating flow induced particle transport coupled with collision of particles and shear interaction between them leads to the formation of dome shaped viscoelastic shells of different dimensions depending on the surface temperature. These shells undergo sol-gel transition and subsequently undergo buckling instability leading to the formation of daughter cavities. With an increase in the surface temperature, droplets exhibit buckling from multiple sites over a larger sector in the top half of the droplet. Irrespective of the initial nanoparticle concentration and substrate temperature, growth of a daughter cavity (subsequent to buckling) inside the droplet is found to be controlled by the solvent evaporation rate from the droplet periphery and is shown to exhibit a universal trend.

Response And Stability Of A Stochastic Beam-Column Using Stochastic Fem

A beam-column resting on continuous Winkler foundation and discrete elastic supports is considered. The beam-column is of variable cross-section and the variation of sectional properties along the axis of the beam-column is deterministic. Young's modulus, mass per unit length and distributed axial loadings of the beam-column have a stochastic distribution. The foundation stiffness coefficient of the Winkler model, the stiffnesses of discrete elastic supports, stiffnesses of end springs and the end thrust, are all considered as random parameters. The material property fluctuations and distributed axial loadings are considered to constitute independent, one-dimension uni-variate homogeneous real stochastic fields in space. The foundation stiffness coefficient, stiffnesses of the discrete elastic supports, stiffnesses of end springs and the end thrust are considered to constitute independent random variables. Static response, free vibration and stability behaviour of the beam-column are studied. Hamilton's principle is used to formulate the problem using stochastic FEM. Sensitivity vectors of the response and stability parameters are evaluated. Using these statistics of free vibration frequencies, mode shapes, buckling parameters, etc., are evaluated. A numerical example is given.

Energy-momentum conserving algorithm for nonlinear transient analysis within the framework of hybrid elements

This work deals with the formulation and implementation of an energy-momentum conserving algorithm for conducting the nonlinear transient analysis of structures, within the framework of stress-based hybrid elements. Hybrid elements, which are based on a two-field variational formulation, are much less susceptible to locking than conventional displacement-based elements within the static framework. We show that this advantage carries over to the transient case, so that not only are the solutions obtained more accurate, but they are obtained in fewer iterations. We demonstrate the efficacy of the algorithm on a wide range of problems such as ones involving dynamic buckling, complicated three-dimensional motions, et cetera.

Instability of laminated composite thin-walled open-section beams

Instability of thin-walled open-section laminated composite beams is studied using the finite element method. A two-noded, 8 df per node thin-walled open-section laminated composite beam finite element has been used. The displacements of the element reference axis are expressed in terms of one-dimensional first order Hermite interpolation polynomials, and line member assumptions are invoked in formulation of the elastic stiffness matrix and geometric stiffness matrix. The nonlinear expressions for the strains occurring in thin-walled open-section beams, when subjected to axial, flexural and torsional loads, are incorporated in a general instability analysis. Several problems for which continuum solutions (exact/approximate) are possible have been solved in order to evaluate the performance of finite element. Next its applicability is demonstrated by predicting the buckling loads for the following problems of laminated composites: (i) two layer (45°/−45°) composite Z section cantilever beam and (ii) three layer (0°/45°/0°) composite Z section cantilever beam.

Morphological transitions in polymer monolayers under compression

We present a systematic investigation of morphological transitions in poly vinylacetate Langmuir monolayers. On compression, the polymer monolayer is converted to a continuous membrane with a thickness of similar to 2-3 nm. Above a certain surface concentration the monolayer, on water, undergoes a morphological transition-buckling, leading to formation of striped patterns of period of lambda(b)similar to 160 nm, as determined from in situ grazing incidence small angle x-ray scattering measurements. The obtained value is much smaller than what has been typically observed for Langmuir monolayers on water or thin films on soft substrates. Using existing theories for buckling of fluidlike films on fluid substrates, we obtain very low values of bending rigidity and Young's modulus of the polymer monolayer compared to that observed earlier for lipid or polymeric monolayers. Since buckling in these monolayers occurs only above a certain surface concentration, we have looked at the possibility that the buckling in these films occurs due to changes in their mechanical properties under compression. Using the model of Huang and Suo of buckling of solidlike films on viscoelastic substrates, we find values of the mechanical properties, which are much closer to the bulk values but still significantly lower. Although the reduction could be along the lines of what has been observed earlier for ultrathin polymer film or surface layers of polymers, the possibility of micromechanical effects also determining the buckling in such polymer monolayers cannot be ruled out. We have provided possible explanation of the buckling of the poly vinylacetate monolayers in terms of the change in isothermal compression modulus with surface concentration.

Macromechanical behaviour of composite laminates

An improved higher order transverse shear deformation theory is employed to arrive at modified constitutive relations which can be used in the flexural, buckling and vibration analysis of laminated plates and shells. The strain energy for such systems is then expressed in terms of the displacements and the rotations for ready reference and use. Numerical values of vibration frequencies are obtained using this formulation employing Ritz's method of analysis. The results are compared with those available in the literature to validate the analysis presented.

Simultaneous material selection and geometry design of statically determinate trusses using continuous optimization

In this work, we explore simultaneous geometry design and material selection for statically determinate trusses by posing it as a continuous optimization problem. The underlying principles of our approach are structural optimization and Ashby’s procedure for material selection from a database. For simplicity and ease of initial implementation, only static loads are considered in this work with the intent of maximum stiffness, minimum weight/cost, and safety against failure. Safety of tensile and compression members in the truss is treated differently to prevent yield and buckling failures, respectively. Geometry variables such as lengths and orientations of members are taken to be the design variables in an assumed layout. Areas of cross-section of the members are determined to satisfy the failure constraints in each member. Along the lines of Ashby’s material indices, a new design index is derived for trusses. The design index helps in choosing the most suitable material for any geometry of the truss. Using the design index, both the design space and the material database are searched simultaneously using gradient-based optimization algorithms. The important feature of our approach is that the formulated optimization problem is continuous, although the material selection from a database is an inherently discrete problem. A few illustrative examples are included. It is observed that the method is capable of determining the optimal topology in addition to optimal geometry when the assumed layout contains more links than are necessary for optimality.

Stability Analysis of Composite Skew Plates Using a Dynamic Method

Stability analysis is carried out considering free lateral vibrations of simply supported composite skew plates that are subjected to both direct and shear in-plane forces. An oblique stress component representation is used, consistent with the skew-geometry of the plate. A double series, expressed in Chebyshev polynomials, is used here as the assumed deflection surface and Ritz method of solution is employed. Numerical results for different combinations of side ratios, skew angle, and in-plane loadings that act individually or in combination are obtained. In this method, the in-plane load parameter is varied until the fundamental frequency goes to zero. The value of the in-plane load then corresponds to a critical buckling load. Plots of frequency parameter versus in-plane loading are given for a few typical cases. Details of crossings and quasi degeneracies of these curves are presented.

Ultimate Strength of R.C. Wall Panels in Two-Way In-Plane Action

Test results of 24 reinforced concrete wall panels in two-way action (i.e., supported on all the four sides) and subjected to in-plane vertical load are presented. The load is applied at an eccentricity to represent possible accidental eccentricity that occurs in practice due to constructional imperfections. Influences of aspect ratio, thinness ratio, slendemess ratio, vertical steel, and horizontal steel on the ultimate load are studied. Two equations are proposed to predict the ultimate load carried by the panels. The first equation is empirical and is arrived at from trial and error fitting with test data. The second equation is semi-empirical and is developed from a modification of the buckling strength of thin rectangular plates. Both the equations are formulated so as to give a safe prediction of a large portion of ultimate strength test results. Also, ultimate load cracking load and lateral deflections of identical panels in two-way action (all four sides supported) and oneway action (top and bottom sides only supported) are compared.

Study of stress Relief Patterns In Diamond-Like Carbon-Films

An optical microscopy study of stress relief patterns in diamondlike carbon films is presented. Interesting stress relief patterns are observed which include the well known sinusoidal type, branching pattern and string of beads pattern. The last one is shown to relieve stresses under marginal conditions. Two new stress relief patterns are noted in the present study. One of them is of a sinusoidal shape with two extra branches at every peak position. The distribution of different stress relief forms from the outer edge of the films towards the interior is markedly dependent on film thickness. Our new patterns support the approach in which the stress relief forms have been analysed earlier using the theory of plate buckling.

Local variability and base sequence effects in DNA crystal structures.

The importance and usefulness of local doublet parameters in understanding sequence dependent effects has been described for A- and B-DNA oligonucleotide crystal structures. Each of the two sets of local parameters described by us in the NUPARM algorithm, namely the local doublet parameters, calculated with reference to the mean z-axis, and the local helical parameters, calculated with reference to the local helix axis, is sufficient to describe the oligonucleotide structures, with the local helical parameters giving a slightly magnified picture of the variations in the structures. The values of local doublet parameters calculated by NUPARM algorithm are similar to those calculated by NEWHELIX90 program, only if the oligonucleotide fragment is not too distorted. The mean values obtained using all the available data for B-DNA crystals are not significantly different from those obtained when a limited data set is used, consisting only of structures with a data resolution of better than 2.4 A and without any bound drug molecule. Thus the variation observed in the oligonucleotide crystals appears to be independent of the quality of their crystallinity. No strong correlation is seen between any pair of local doublet parameters but the local helical parameters are interrelated by geometric relationships. An interesting feature that emerges from this analysis is that the local rise along the z-axis is highly correlated with the difference in the buckle values of the two basepairs in the doublet, as suggested earlier for the dodecamer structures (Bansal and Bhattacharyya, in Structure & Methods: DNA & RNA, Vol. 3 (Eds., R.H. Sarma and M.H. Sarma), pp. 139-153 (1990)). In fact the local rise values become almost constant for both A- and B-forms, if a correction is applied for the buckling of the basepairs. In B-DNA the AA, AT, TA and GA basepair sequences generally have a smaller local rise (3.25 A) compared to the other sequences (3.4 A) and this seems to be an intrinsic feature of basepair stacking interaction and not related to any other local doublet parameter. The roll angles in B-DNA oligonucleotides have small values (less than +/- 8 degrees), while mean local twist varies from 24 degrees to 45 degrees. The CA/TG doublet sequences show two types of preferred geometries, one with positive roll, small positive slide and reduced twist and another with negative roll, large positive slide and increased twist.(ABSTRACT TRUNCATED AT 400 WORDS)