A hybrid finite element formulation for flexible multibody dynamics

This work deals with the transient analysis of flexible multibody systems within a hybrid finite element framework. Hybrid finite elements are based on a two-field variational formulation in which the displacements and stresses are interpolated separately yielding very good coarse mesh accuracy. Most of the literature on flexible multibody systems uses beam-theory-based formulations. In contrast, the use of hybrid finite elements uses continuum-based elements, thus avoiding the problems associated with rotational degrees of freedom. In particular, any given three-dimensional constitutive relations can be directly used within the framework of this formulation. Since the coarse mesh accuracy as compared to a conventional displacement-based formulation is very high, the scheme is cost effective as well. A general formulation is developed for the constrained motion of a given point on a line manifold, using a total Lagrangian method. The multipoint constraint equations are implemented using Lagrange multipliers. Various kinds of joints such as cylindrical, prismatic, and screw joints are implemented within this general framework. Hinge joints such as spherical, universal, and revolute joints are obtained simply by using shared nodes between the bodies. In addition to joints, the formulation and implementation details for a DC motor actuator and for prescribed relative rotation are also presented. Several example problems illustrate the efficacy of the developed formulation.

Asymptotic behavior and interpretation of virtual states: The effects of confinement and of basis sets

A real-space high order finite difference method is used to analyze the effect of spherical domain size on the Hartree-Fock (and density functional theory) virtual eigenstates. We show the domain size dependence of both positive and negative virtual eigenvalues of the Hartree-Fock equations for small molecules. We demonstrate that positive states behave like a particle in spherical well and show how they approach zero. For the negative eigenstates, we show that large domains are needed to get the correct eigenvalues. We compare our results to those of Gaussian basis sets and draw some conclusions for real-space, basis-sets, and plane-waves calculations. (C) 2016 AIP Publishing LLC.

A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.

Cold gas in cluster cores: global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (t(cool)/t(ff)). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.

Dynamic Characteristics of Spherically Converging Detonation Waves

The spherically converging detonation wave was numerically investigated by solving the one-dimensional multi-component Euler equations in spherical coordinates with a dispersion-controlled dissipative scheme. Finite rate and detailed chemical reaction models were used and numerical solutions were obtained for both a spherical by converging detonation in a stoichiometric hydrogen-oxygen mixture and a spherically focusing shock in air. The results showed that the post-shock pressure approximately arises to the same amplitude in vicinity of the focal point for the two cases, but the post-shock temperature level mainly depends on chemical reactions and molecular dissociations of a gas mixture. While the chemical reaction heat plays an important role in the early stage of detonation wave propagation, gas dissociations dramatically affect the post-shock flow states near the focal point. The maximum pressure and temperature, non-dimensionalized by their initial value, are approximately scaled to the propagation radius over the initial detonation diameter. The post-shock pressure is proportional to the initial pressure of the detonable mixture, and the post-shock temperature is also increased with the initial pressure, but in a much lower rate than that of the post-shock pressure.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Microstructure and mechanical properties of slowly cooled Cu47.5Zr47.5Al5

Cu47.5Zr47.5Al5 was prepared by arc melting and solidified in situ by suction casting into 2-5-mm-diameter rods under various cooling rates (200-2000 K/s). The microstructure was investigated along the length of the rods by electron microscopy, differential scanning calorimetry and mechanical properties were investigated under compression. The microstructure of differently prepared specimens consists of macroscopic spherical shape chemically inhomogeneous regions together with a low volume fraction of randomly distributed CuZr B2 phase embedded in a 2-7 nm size clustered "glassy-martensite" matrix. The as-cast specimens show high yield strength (1721 MPa), pronounced work-hardening behavior up to 2116 MPa and large fracture strain up to 12.1-15.1%. The fracture strain decreases with increasing casting diameter. The presence of chemical inhomogenities and nanoscale "glassy-martensite" features are beneficial for improving the inherent ductility of the metallic glass.

Synthesis of carbon 'onions' in water.

The fabrication of carbon nanomaterials usually calls for expensive vacuum systems to generate plasmas and yields are disappointingly low. Here we describe a simple method for producing high-quality spherical carbon nano-'onions' in large quantities without the use of vacuum equipment. The nanoparticles, which have C60 cores surrounded by onion-like nested particles, are generated by an arc discharge between two graphite electrodes submerged in water. This technique is economical and environmentally benign, and produces uncontaminated nanoparticles which may be useful in many applications.

Planar Thermocapillary Migration of Two Bubbles in Microgravity Environment

A theoretical investigation is performed on the thermocapillary motion of two bubbles in arbitrary configuration in microgravity environment under the assumption that the surface tension is high enough to keep the bubbles spherical. The two bubbles are dr

The effect of nucleation time on the growth of a microvoid in a viscoelastic material

In this paper, discussions are focused on the growth of a nucleated void in a viscoelastic material. The in situ tensile tests of specimens made of high-density polyethylene, filled with spherical glass beads (HDPE/GB) are carried out under SEM. The experimental result indicates that the microvoid nucleation is induced by the partially interfacial debonding of particles. By means of the Laplace transform and the Eshelby's equivalent inclusion method, a new analytical expression of the void strain at different nucleation times is derived. It can be seen that the strain of the nucleated void depends not only on the remote strain history, but also on the nucleation time. This expression is also illustrated by numerical examples, and is found to be of great usefulness in the study of damage evolution in viscoelastic materials.

Nonlinear Analysis of Oscillatory Indentation in Elastic and Viscoelastic Solids

Determining the mechanical properties at micro- and nanometer length scales using nanoindentation or atomic force microscopy is important to many areas of science and engineering. Here we establish equations for obtaining storage and loss modulus from oscillatory indentations by performing a nonlinear analysis of conical and spherical indentation in elastic and viscoelastic solids. We show that, when the conical indenter is driven by a sinusoidal force, the square of displacement is a sinusoidal function of time, not the displacement itself, which is commonly assumed. Similar conclusions hold for spherical indentations. Well-known difficulties associated with measuring contact area and correcting thermal drift may be circumvented using the newly derived equations. These results may help improve methods of using oscillatory indentation for determining elastic and viscoelastic properties of solids.

扩散处理热浸镀铝钢高温抗氧化行为的研究

Embedded cell model of three-dimensional two-phase composites

An embedded cell model is presented to obtain the effective elastic moduli and the elastic-plastic stress-strain relations of three-dimensional two-phase particulate composites. Each cell consists of an ellipsoidal inclusion surrounded by a finite ellipsoidal matrix that embedded in an infinite matrix. When both matrix and particle are elastic, the effective elastic moduli are derived which is an exact analytic formula without any simplified approximation that can be expressed in an explicit form. Further, the elastic-plastic stress-strain relations are obtained for spherical cells and oblate spheroid cells, in which the matrix is elastic and the particle is elastic-plastic. In addition, the macroscopic elastic-plastic constitutive relation of particle reinforced composites (PRC) is investigated by a systematic approach [1] in which the matrix is elastic-plastic and the particle is elastic.

The effect of straining mode on void growth

Large strain finite element method is employed to investigate the effect of straining mode on void growth. Axisymmetric cell model embedded with spherical void is controlled by constant triaxiality: loading,while plane-stress model containing a circular void is loaded by constant ratio of straining. Elastic-plastic material is used for the matrix in both cases. It is concluded that, besides the known effect of triaxiality, the straining mode which intensifies the plastic concentration around the void is also a void growth stimulator. Experimental results are cited to justify the computation results.