Transformation field method for calculating the effective properties of isotropic graded composites having arbitrary shapes

A method of transformation field is developed to estimate the effective properties of graded composites whose inclusions have arbitrary shapes and gradient profiles by means of a periodic cell model. The boundary-value problem of graded composites having arbitrary inclusion shapes is solved by introducing the transformation field into the inclusion region. As an example, the effective dielectric response of isotropic graded composites having arbitrary shapes and gradient profiles is handled by the transformation field method (TFM). Moreover, TFM results are validated by the exact solutions of isotropic graded spherical inclusions having a power-law profile and good agreement is obtained in the dilute limit. Furthermore, it is found that the inclusion shapes and the parameters of the gradient profiles can have profound effect on the effective properties of composite systems at high concentration of inclusions.

Nine microsatellite DNA primers for Hippophae rhamnoides ssp sinensis (Elaeagnaceae)

Hippophae rhamnoides ssp. sinensis occurs mainly in the and regions of northwest China. The wood stands of this subspecies play an important role in maintaining the local ecosystems in these regions. In addition, the genetic characteristics are essential to understand the historical range changes of this subspecies and its morphological differentiation with other subspecies. In this study, we developed nine microsatellite loci for this subspecies for the first time. We used the combining biotin capture method to enrich AG/CT/AC/GT/CG/GTG/CCA microsatellites. Twenty-six microsatellites were isolated from the enriching library and nine of them were found to be polymorphic through screening 12 distantly distributed individuals. The number of alleles per locus ranged from three to twelve and expected heterozygosity from 0.2659 to 0.4767, respectively. We further performed cross-priming tests in another subspecies and two congeneric species. These firstly isolated loci will provide a useful tool to investigate the genetic structure of this subspecies and its morphological differentiation from the other subspecies.

Isolation and characterization of microsatellite DNA primers in Juniperus przewalskii Kom (Cupressaceae)

Juniperus przewalskii (Cupressaceae) is a dominant tree species endemic to the northeast Qinghai-Tibetan Plateau. This species plays an important role in maintaining the arid ecosystem in this region. However, natural distributions of this species have been declined. In order to develop effective conservation methods, it is important to know the distribution of the genetic diversity within and among populations. In this study, we developed nine new microsatellite loci for this species. We used the combining biotin capture method to enrich AG/CT/AC/GGT microsatellites. The polymorphisms of each locus were further assessed in 12 individuals from four geographically distant populations. The number of alleles per locus varied from three to six and expected heterozygosity ranged from 0.58 to 0.70. These loci together provide a useful tool to investigate the genetic diversity of this species. In addition, all markers have been crossly checked in the other four congeneric species.

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.

Relationship Between Contact Stiffness, Contact Depth, and Mechanical Properties for Indentation in Linear Viscoelastic Solids Using Axisymmetric Indenters

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. Although the same expres

Three-Dimensional Network Model And Its Parameter Calibration

A material model, whose framework is parallel spring-bundles oriented in 3-D space, is proposed. Based on a discussion of the discrete schemes and optimum discretization of the solid angles, a 3-D network cell consisted of one-dimensional components is developed with its geometrical and physical parameters calibrated. It is proved that the 3-D network model is able to exactly simulate materials with arbitrary Poisson ratio from 0 to 1/2, breaking through the limit that the previous models in the literature are only suitable for materials with Poisson ratio from 0 to 1/3. A simplified model is also proposed to realize high computation accuracy within low computation cost. Examples demonstrate that the 3-D network model has particular superiority in the simulation of short-fiber reinforced composites.

Elastic fields of interacting elliptic inhomogeneities

This paper presents a method for the calculation of two-dimensional elastic fields in a solid containing any number of inhomogeneities under arbitrary far field loadings. The method called 'pseudo-dislocations method', is illustrated for the solution of interacting elliptic inhomogeneities. It reduces the interacting inhomogeneities problem to a set of linear algebraic equations. Numerical results are presented for a variety of elliptic inhomogeneity arrangements, including the special cases of elliptic holes, cracks and circular inhomogeneities. All these complicated problems can be solved with high accuracy and efficiency.

Numerical Solutions for the Transient Flow in the Homogenous Closed Circle Reservoirs

There are many fault block fields in China. A fault block field consists of fault pools. The small fault pools can be viewed as the closed circle reservoirs in some case. In order to know the pressure change of the developed formation and provide the formation data for developing the fault block fields reasonably, the transient flow should be researched. In this paper, we use the automatic mesh generation technology and the finite element method to solve the transient flow problem for the well located in the closed circle reservoir, especially for the well located in an arbitrary position in the closed circle reservoir. The pressure diffusion process is visualized and the well-location factor concept is first proposed in this paper. The typical curves of pressure vs time for the well with different well-location factors are presented. By comparing numerical results with the analytical solutions of the well located in the center of the closed circle reservoir, the numerical method is verified.

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

Influence of van der Waals and Casimir Forces on Electrostatic Torsional Actuators

The influence of van der Waals (vdW) and Casimir forces on the stability of the electrostatic torsional nanoelectromechanical systems (NEMS) actuators is analyzed in the paper. With the consideration of vdW and Casimir effects, the dependence of the critical tilting angle and pull-in voltage on the sizes of structure is investigated. And the influence of vdW torque is compared with that of Casimir torque. The modified coefficients of vdW and Casimir torques on the pull-in voltage are, respectively, calculated. When the gap is sufficiently small, pull-in can still take place with arbitrary small angle perturbation because of the action of vdW and Casimir torques even if there is not electrostatic torque. And the critical pull-in gaps for two cases are, respectively, derived.

Regularity of strain distribution in short-fiber/whisker reinforced composites

Based on studies on the strain distribution in short-fiber/whisker reinforced metal matrix composites, a deformation characteristic parameter, lambda is defined as a ratio of root-mean-square strain of the reinforcers identically oriented to the macro-linear strain along the same direction. Quantitative relation between lambda and microstructure parameters of composites is obtained. By using lambda, the stiffness moduli of composites with arbitrary reinforcer orientation density function and under arbitrary loading condition are derived. The upper-bound and lower-bound of the present prediction are the same as those from the equal-strain theory and equal-stress theory, respectively. The present theory provides a physical explanation and theoretical base for the present commonly-used empirical formulae. Compared with the microscopic mechanical theories, the present theory is competent for stiffness modulus prediction of practical engineering composites in accuracy and simplicity.

Dynamic equations for curved submerged floating tunnel

In virtue of reference Cartesian coordinates, geometrical relations of spatial curved structure are presented in orthogonal curvilinear coordinates. Dynamic equations for helical girder are derived by Hamilton principle. These equations indicate that four generalized displacements are coupled with each other. When spatial structure degenerates into planar curvilinear structure, two generalized displacements in two perpendicular planes are coupled with each other. Dynamic equations for arbitrary curvilinear structure may be obtained by the method used in this paper.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

On the initial unloading slope in indentation of elastic-plastic solids by an indenter with an axisymmetric smooth profile

A simple relationship between the initial unloading slope, the contact area, and the elastic modulus is derived for indentation in elastic-plastic solids by an indenter with an arbitrary axisymmetric smooth profile. Although the same expression was known to hold for elastic solids, the new derivation shows that it is also true for elastic-plastic solids with or without work hardening and residual stress. These results should provide a sound basis for the use of the relationship for mechanical property determination using indentation techniques. (C) 1997 American Institute of Physics.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.