VIBRATIONAL CHARACTERISTICS OF ZIGZAG, ARMCHAIR AND CHIRAL CANTILEVER SINGLE-WALLED CARBON NANOTUBES

Finite element analysis has been performed to study vibrational characteristics of cantilever single walled carbon nanotubes. Finite element models are generated by specifying the C-C bond rigidities, which are estimated by equating energies from molecular mechanics and continuum mechanics. Bending, torsion, and axial modes are identified based on effective mass for armchair, zigzag and chiral cantilever single walled carbon nanotubes, whose Young's modulus is evaluated from the bending frequency. Empirical relations are provided for frequencies of bending, torsion, and axial modes.

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.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

On the invariant representation of spin tensors with applications

The invariant representation of the spin tensor defined as the rotation rate of a principal triad for a symmetric and non-degenerate tensor is derived on the basis of the general solution of a linear tensorial equation. The result can be naturally specified to study the. spin of the stretch tensors and to investigate the relations between various rotation rate tensors encountered frequently in modern continuum mechanics. A remarkable formula which relates the generalized stress conjugate to the generalized strain in Hill's sense. to Cauchy stress, is obtained in invariant form through the work conjugate principle. Particularly, a detailed discussion on the time rate of logarithmic strain and its conjugate stress is made as the principal axes of strain arc not fixed during deformation.

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Atomistic Simulation On Size-Dependent Yield Strength And Defects Evolution Of Metal Nanowires

A simple derivation based on continuum mechanics is given, which shows the surface stress is critical for yield strength at ultra-small scales. Molecular dynamics (MD) simulations with modified embedded atom method (MEAM) are employed to investigate the mechanical behaviors of single-crystalline metal nanowires under tensile loading. The calculated yield strengths increasing with the decrease of the cross-sectional area of the nanowires are in accordance with the theoretical prediction. Reorientation induced by stacking faults is observed at the nanowire edge. In addition. the mechanism of yielding is discussed in details based on the snapshots of defects evolution. The nanowires in different crystallographic orientations behave differently in stretching deformation. This study on the plastic properties of metal nanowires will be helpful to further understanding of the mechanical properties of nanomaterials. (C) 2009 Elsevier B.V. All rights reserved.

Compressible flows at small Reynolds numbers

The problem of the slow viscous flow of a gas past a sphere is considered. The fluid cannot be treated incompressible in the limit when the Reynolds number Re, and the Mach number M, tend to zero in such a way that Re ~ o(M^2 ). In this case, the lowest order approximation to the steady Navier-Stokes equations of motion leads to a paradox discovered by Lagerstrom and Chester. This paradox is resolved within the framework of continuum mechanics using the classical slip condition and an iteration scheme that takes into account certain terms in the full Navier-Stokes equations that drop out in the approximation used by the above authors. It is found however that the drag predicted by the theory does not agree with R. A. Millikan's classic experiments on sphere drag.

The whole question of the applicability of the Navier-Stokes theory when the Knudsen number M/Re is not small is examined. A new slip condition is proposed. The idea that the Navier-Stokes equations coupled with this condition may adequately describe small Reynolds number flows when the Knudsen number is not too large is looked at in some detail. First, a general discussion of asymptotic solutions of the equations for all such flows is given. The theory is then applied to several concrete problems of fluid motion. The deductions from this theory appear to interpret and summarize the results of Millikan over a much wider range of Knudsen numbers (almost up to the free molecular or kinetic limit) than hitherto Believed possible by a purely continuum theory. Further experimental tests are suggested and certain interesting applications to the theory of dilute suspensions in gases are noted. Some of the questions raised in the main body of the work are explored further in the appendices.

Extracting material response from simple mechanical tests on hardening-softening-hardening viscoplastic solids

Compliant foams are usually characterized by a wide range of desirable mechanical properties. These properties include viscoelasticity at different temperatures, energy absorption, recoverability under cyclic loading, impact resistance, and thermal, electrical, acoustic and radiation-resistance. Some foams contain nano-sized features and are used in small-scale devices. This implies that the characteristic dimensions of foams span multiple length scales, rendering modeling their mechanical properties difficult. Continuum mechanics-based models capture some salient experimental features like the linear elastic regime, followed by non-linear plateau stress regime. However, they lack mesostructural physical details. This makes them incapable of accurately predicting local peaks in stress and strain distributions, which significantly affect the deformation paths. Atomistic methods are capable of capturing the physical origins of deformation at smaller scales, but suffer from impractical computational intensity. Capturing deformation at the so-called meso-scale, which is capable of describing the phenomenon at a continuum level, but with some physical insights, requires developing new theoretical approaches.

A fundamental question that motivates the modeling of foams is ‘how to extract the intrinsic material response from simple mechanical test data, such as stress vs. strain response?’ A 3D model was developed to simulate the mechanical response of foam-type materials. The novelty of this model includes unique features such as the hardening-softening-hardening material response, strain rate-dependence, and plastically compressible solids with plastic non-normality. Suggestive links from atomistic simulations of foams were borrowed to formulate a physically informed hardening material input function. Motivated by a model that qualitatively captured the response of foam-type vertically aligned carbon nanotube (VACNT) pillars under uniaxial compression [2011,“Analysis of Uniaxial Compression of Vertically Aligned Carbon Nanotubes,” J. Mech.Phys. Solids, 59, pp. 2227–2237, Erratum 60, 1753–1756 (2012)], the property space exploration was advanced to three types of simple mechanical tests: 1) uniaxial compression, 2) uniaxial tension, and 3) nanoindentation with a conical and a flat-punch tip. The simulations attempt to explain some of the salient features in experimental data, like
1) The initial linear elastic response.
2) One or more nonlinear instabilities, yielding, and hardening.

The model-inherent relationships between the material properties and the overall stress-strain behavior were validated against the available experimental data. The material properties include the gradient in stiffness along the height, plastic and elastic compressibility, and hardening. Each of these tests was evaluated in terms of their efficiency in extracting material properties. The uniaxial simulation results proved to be a combination of structural and material influences. Out of all deformation paths, flat-punch indentation proved to be superior since it is the most sensitive in capturing the material properties.

Material model of lung parenchyma based on living precision-cut lung slice testing.

We describe a novel constitutive model of lung parenchyma, which can be used for continuum mechanics based predictive simulations. To develop this model, we experimentally determined the nonlinear material behavior of rat lung parenchyma. This was achieved via uni-axial tension tests on living precision-cut rat lung slices. The resulting force-displacement curves were then used as inputs for an inverse analysis. The Levenberg-Marquardt algorithm was utilized to optimize the material parameters of combinations and recombinations of established strain-energy density functions (SEFs). Comparing the best-fits of the tested SEFs we found Wpar = 4.1 kPa(I1-3)2 + 20.7 kPa(I1 - 3)3 + 4.1 kPa(-2 ln J + J2 - 1) to be the optimal constitutive model. This SEF consists of three summands: the first can be interpreted as the contribution of the elastin fibers and the ground substance, the second as the contribution of the collagen fibers while the third controls the volumetric change. The presented approach will help to model the behavior of the pulmonary parenchyma and to quantify the strains and stresses during ventilation.

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113-2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent -1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.

Evolution of Localized Damage Zone in Heterogeneous Media

Evolution of localized damage zone is a key to catastrophic rupture in heterogeneous materials. In the present article, the evolutions of strain fields of rock specimens are investigated experimentally. The observed evolution of fluctuations and autocorrelations of strain fields under uniaxial compression demonstrates that the localization of deformation always appears ahead of catastrophic rupture. In particular, the localization evolves pronouncedly with increasing deformation in the rock experiments. By means of the definition of the zone with high strain rate and likely damage localization, it is found that the size of the localized zone decreases from the sample size at peak load to an eventual value. Actually, the deformation field beyond peak load is bound to suffer bifurcation, namely an elastic unloading part and a continuing but localized damage part will co-exist in series in a specimen. To describe this continuous bifurcation and localization process observed in experiments, a model on continuum mechanics is developed. The model can explain why the decreasing width of localized zone can lead stable deformation to unstable, but it still has not provided the complete equations governing the evolution of the localized zone.

基于连续介质力学的块体单元离散弹簧法研究

提出一种基于连续介质力学的块体单元离散弹簧法,并将其应用于地质体渐进破坏过程的数值模拟研究。该方法以连续介质力学的基本理论为基础,以八节点正六面体单元为离散对象,根据有限元的形函数理论及正六面体单元刚度矩阵的具体表达式,将正六面体单元离散成为12根棱弹簧,并推导各弹簧的弹簧力计算公式,给出各弹簧的弹簧力分量的物理意义、对应的弹簧刚度及相关联的位移。在弹性部分的计算结果与传统的有限元计算结果一致的基础上,在弹簧中引入Mohr-Coulomb准则及拉伸破坏准则,进行弹簧破坏的判断。如果块体内的弹簧全部断裂,便将其视为散体。在判断破坏时采用双重判断模式,块体单元的判断主要用于选取弹簧潜在的破坏面,离散弹簧的判断用于计算弹簧上的真实弹簧力。最后,通过相关算例验证引入破坏准则后所得计算结果的合理性。块体单元离散弹簧法的实质是通过12根离散弹簧将单元转化为结构,因此可以通过分析各弹簧的断裂情况研究块体的内部破坏特征,更可以通过各弹簧的先后断裂过程来研究地质体的渐进破坏过程。

CFD analysis of one-dimensional infiltration in vadose zone

This study presents a CFD analysis constructed around PHYSICA, an open framework for multi-physics computational continuum mechanics modelling, to investigate the water movement in unsaturated porous media. The modelling environment is based on a cell-centred finite-volume discretisation technique. A number of test cases are performed in order to validate the correct implementation of Richard's equation for compressible and incompressible fluids. The pressure head form of the equation is used together with the constitutive relationships between pressure, volumetric water content and hydraulic conductivity described by Haverkamp and Van Genuchten models. The flow problems presented are associated with infiltration into initially dry soils with homogeneous or layered geologic settings. Comparison of results with the problems selected from literature shows a good agreement and validates the approach and the implementation.

Simulation de la nage anguilliforme

Ce document traite premièrement des diverses tentatives de modélisation et de simulation de la nage anguilliforme puis élabore une nouvelle technique, basée sur la méthode de la frontière immergée généralisée et la théorie des poutres de Reissner-Simo. Cette dernière, comme les équations des fluides polaires, est dérivée de la mécanique des milieux continus puis les équations obtenues sont discrétisées afin de les amener à une résolution numérique. Pour la première fois, la théorie des schémas de Runge-Kutta additifs est combinée à celle des schémas de Runge-Kutta-Munthe-Kaas pour engendrer une méthode d’ordre de convergence formel arbitraire. De plus, les opérations d’interpolation et d’étalement sont traitées d’un nouveau point de vue qui suggère l’usage des splines interpolatoires nodales en lieu et place des fonctions d’étalement traditionnelles. Enfin, de nombreuses vérifications numériques sont faites avant de considérer les simulations de la nage.

On the regular-geometric-figure solution to the N-body problem

The regular-geometric-figure solution to the N-body problem is presented in a very simple way. The Newtonian formalism is used without resorting to a more involved rotating coordinate system. Those configurations occur for other kinds of interactions beyond the gravitational ones for some special values of the parameters of the forces. For the harmonic oscillator, in particular, it is shown that the N-body problem is reduced to N one-body problems.