基于界面层模型的双周期纳米夹杂复合材料反平面问题研究

利用复变函数方法并结合双准周期Riemann边值问题理论，获得了含双周期分布非均匀相（夹杂/界面层）的复合材料在远场均匀反平面应力下弹性场的全场解答．该解答可用于对纳米夹杂复合材料的应力进行分析，结合平均场理论也用于预测纳米夹杂复合材料的有效性能．计算结果表明：当夹杂尺度在纳米量级时，应力和有效反平面剪切模量具有明显的尺度依赖性，并且随着夹杂尺寸的增加，趋近于不考虑界面效应时的结果；界面层厚度和性能对应力和有效反平面剪切模量明显变化时所对应的夹杂尺度范围和趋近于无界面效应结果的快慢有显著影响；当界面厚度足够薄时，界面层模型可用于模拟零厚度界面情况．     

A Cosserat shell model for interphases in elastic media

This study is concerned with the modeling of interphases in elastic media in general, and in composite materials in particular. The aim is to replace a boundary value problem consisting of a three-phase configuration, say that of fiber-interphase-matrix, by a simpler problem which involves the fiber and matrix only, plus certain matching conditions which simulate the interphase. The simplest of such known representations replaces a thin interphase by a “perfect contact interface” (a single surface) across which the displacements and tractions are assumed to be continuous. Another classical model replaces a thin and soft interphase by a “spring-type interface”, across which the tractions are continuous, but the displacement field undergoes a discontinuity. In the present paper, a Cosserat shell model of the interphase is derived which successfully models the original interphase in a unified manner, for the full range of its material parameters relative to those of the neighboring media. The model is derived in the setting of three-dimensional linear elasticity with small deformations and displacements. Comparisons with an existing exact solution of a coated fiber in an infinite matrix show that it performs extremely well even for moderately thick interphases.     

Contact analysis in the presence of an ellipsoidal inhomogeneity within a half space

Many materials contain inhomogeneities or inclusions that may greatly affect their mechanical properties. Such inhomogeneities are for example encountered in the case of composite materials or materials containing precipitates. This paper presents an analysis of contact pressure and subsurface stress field for contact problems in the presence of anisotropic elastic inhomogeneities of ellipsoidal shape. Accounting for any orientation and material properties of the inhomogeneities are the major novelties of this work. The semi-analytical method proposed to solve the contact problem is based on Eshelby’s formalism and uses 2D and 3D Fast Fourier Transforms to speed up the computation. The time and memory necessary are greatly reduced in comparison with the classical finite element method. The model can be seen as an enrichment technique where the enrichment fields from the heterogeneous solution are superimposed to the homogeneous problem. The definition of complex geometries made by combination of inclusions can easily be achieved. A parametric analysis on the effect of elastic properties and geometrical features of the inhomogeneity (size, depth and orientation) is proposed. The model allows to obtain the contact pressure distribution – disturbed by the presence of inhomogeneities – as well as subsurface and matrix/inhomogeneity interface stresses. It is shown that the presence of an inclusion below the contact surface affects significantly the contact pressure and subsurfaces stress distributions when located at a depth lower than 0.7 times the contact radius. The anisotropy directions and material data are also key elements that strongly affect the elastic contact solution. In the case of normal contact between a spherical indenter and an elastic half space containing a single inhomogeneity whose center is located straight below the contact center, the normal stress at the inhomogeneity/matrix interface is mostly compressive. Finally when the axes of the ellipsoidal inclusion do not coincide with the contact problem axes, the pressure distribution is not symmetrical.     

An efficient approximate numerical method for modeling contact of materials with distributed inhomogeneities

This research explores the influence of distributed non-interpenetrating inhomogeneities on the contact of inhomogeneous materials via a new efficient numerical model based on Eshelby’s Equivalent Inclusion Method. The half-space contact of a sphere with an inhomogeneous material is considered, and the solutions take into account interactions between all inhomogeneities. The efficiency and solution accuracy of the proposed method are demonstrated through comparative studies with those of an existing numerical method and the finite element method. The influence of spatial inhomogeneity orientations on the contact elastic field is investigated and parametric studies are conducted for the effect of arbitrarily distributed inhomogeneities on the stress field of the materials. The significance of the influences of inhomogeneity distribution parameters on the inverse volumetric stress integral is quantified and the corresponding data are fitted into selected several formulas as a step towards understanding the rolling contact fatigue life of the materials.     

Full Range Scattering Estimates and Their Application to Cloaking

We establish very precise estimates for the time harmonic scattering effects of an inhomogeneity. Our estimates are valid at all frequencies, and are independent of the contents of the inhomogeneity. The involved constants are independent of the frequency. We use these estimates to assess the effectivity of approximate electromagnetic cloaks constructed by so called “mapping techniques”.     

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.     

Time-harmonic dynamical stress field in a system comprising a pre-stressed orthotropic layer and pre-stressed orthotropic half-plane

The time-harmonic dynamical stress field in a system comprising a pre-stressed orthotropic layer and orthotropic half-plane is studied within the scope of the piecewise homogeneous body model utilizing the three-dimensional linearized theory of elastic waves in an initially stressed body. The main focus is on the influence of the mechanical properties of the constituent materials and the initial stresses present on the “resonance” values of the normal stress acting on the interface plane and on the “resonance” values of the frequency of the external point-located force. The numerical results are presented and discussed. In particular, it is shown that the values of the normal stress decrease with a decrease in the modulus of elasticity of the materials along the thickness of the covering layer.     

Void growth and cavitation in nonlinear viscoelastic solids

This paper discusses the growth of a pre-existing void in a nonlinear viscoelastic material subjected to remote hydrostatic tensions with different loading rates. The constitutive relation of this viscoelastic material is the one recently proposed by the present authors, which may be considered as a generalization of the non-Gaussian statistical theory in rubber elasticity. As the first order approx-imation, the above constitutive relation can be reduced to the “neo-Hookean” type viscoelastic one.Investigations of the influences of the material viscosity and the loading rate on the void growth, or on the cavitation are carried out. It is found that: (1) for generalized “inverse Langevin approximation” nonlinear viscoelastic materials, the cavitation limit does not exist, but there is a certain (remote) stress level at which the void will grow rapidly; (2) for generalized “Gaussian statistics” (neo-Hookean type) viscoelastic materials, the cavitation limit exists, and is an increasing function of the loading rate.The present discussions may be of importance in understanding the material failure process under high triaxial stress.     

Floating Bodies in Neutral Equilibrium

In his paper preceding in this issue, Finn proved that if the contact angle γ of a convex body B{\mathcal{B}} with a given liquid is π/2, and if B{\mathcal{B}} can be made to float in “neutral equilibrium” in the liquid in any orientation, then B{\mathcal{B}} is a metric ball. The present work extends that result, with an independent proof, to any contact angle in the range 0 < γ < π. Our result is equivalent to the general geometric theorem that if for every orientation of a plane, it can be translated to meet a given strictly convex body B{\mathcal{B}} in a fixed angle γ within the above range, then B{\mathcal{B}} is a metric ball.     

Behavior of Piezoelectric Beams under Axially Non-uniform Distributed Loading

The paper presents a generic solution methodology for a quasi-static homogeneous monoclinic piezoelectric beam under axially distributed electric and mechanical surface loads and body forces expressed as polynomials of degree K≥ 0 of the axis variable. (In the absence of any electrical loading, this problem is known as the Almansi–Michell problem). The stress and the electrical displacement components are presented as a set of polynomials of degree ≤K+2 of the axis variable (“solution hypothesis”) containing 4K unknown tip loading constants and 3K stress functions of two variables. The cases K=0,1 stand for uniform or linear distributed loads in the axis direction. The analysis is initiated by the Kth level and continues down to lower levels. The main result of this work generalizes the “elastic” solution given recently by O. Rand and the first author (2005). Examples of solutions for axially uniform distributed loads (K=0), and equilibrium in which the stress and the electrical displacement do not depend on the axis variable, are presented. The applications to constant body loads and a hydrostatic pressure are considered.      

Thermocapillary flow near a “cold corner”

The velocity field in a neighborhood of the point of contact between the free and solid boundaries is studied numerically for the problem of noncrucible zone melting in a two-dimensional model formulation. A distinct Prandtl boundary layer on the solid boundary and a Marangoni boundary layer on the free boundary and high gradients of the longitudinal velocity along the free boundary in the immediate vicinity of the “cold corner” are observed. It is found for the first time that with distance from the solid boundary, the velocity curve has a maximum, which is not typical of the ordinary flow near the solid boundary. Lavrent'ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 3, pp. 141–148, May–June, 1998.     

Crushability maps for structural polymeric foams in uniaxial loading under rigid confinement

A family of epoxy-based polymeric foams with various initial porosity levels was subjected to quasi-static uniaxial loading in rigid confinement (uniaxial strain) to investigate their crushability characteristics. Two issues were investigated. The first issue was the uniformity of deformation in a specimen as a function of porosity level by creating a grid of equally spaced thin stripes on the surface and by monitoring their pattern during the experiment. It was found that the higher the porosity of foam, the more non-uniform the deformation in the specimen. However, the localized non-uniform deformation did not affect the global stress-strain response, especially at large deformations. The second issue was the development of a new analysis tool, called “crushability map”. The purpose of the tool is to depict the evolution of porosity, bulk density and energy absorption as functions of applied strain, stress, and porosity. These maps can assist in characterizing the residual crushability or energy absorption capability of foams as a function of residual porosity. The maps can be used as a design tool for selection of suitable foams for a given application in conjunction with various design criteria.     

Convection of magnetic fluids in connected channels heated from below

The convective flows of a binary mixture in connected channels heated from below are studied experimentally. In contrast to homogeneous fluids, in magnetic colloids “hard” convection excitation, specific transient flows, and oscillatory convection regimes can be observed. The temperature fields and concentration inhomogeneities are measured.     

Nonlinearly Elastic Shell Models: A Formal Asymptotic Approach II. The Flexural Model

This paper is the sequel of Part I, in which the limiting displacement field of a thin shell when its thickness approaches zero is identified as the solution of a two‐dimensional nonlinear membrane shell model. When the geometry of the middle surface of the shell and the boundary conditions allow non‐zero “inextensional displacements”, the previous membrane limit model is not relevant. In this case, we show how to “update” the assumptions on the applied forces acting on the shell so that a limiting model can be derived by an asymptotic analysis. Furthermore, we identify this limit as the two‐dimensional nonlinear flexural shell model. (Accepted January 13, 1997)     

A new approach to predicting interfacial crack growth

I.Intr0ductionItisknownthatthenear-tipstressfieldforaninterfacecrackbetweendissimilarelasticsolidsisalinearcombinationoftwotypesofsingularities,namelyacoupledoscillatoryfieldscaledbyacomplexK(K=K1 iK,whichiscalledinterfacialstressintensityfactor)andanonos…     

Overall Solution-Difference Bounds on the Effects of Material Inhomogeneities

In this work variational bounds are developed on the difference between the solutions of two boundary-value problems involving linearly elastic material microstructure. Both problems have the same external geometry, boundary conditions and loading, however, one possesses no microstructure, i.e., it is homogeneous, while the other is “perturbed” with inhomogeneities. The bounds developed are solely in terms of the solution to the easier homogeneous material problem, the microstructure of the inhomogeneous problem and the given loading data. There are no unknown constants in the bounds. Furthermore, no calculation of the harder, typically intractable, inhomogeneous material problem is necessary. The bounds developed are obtained under no assumptions on the character of the microstructure of the inhomogeneous material problem, other than it be pointwise positive-definite, as well as under no assumptions on the external loading and geometry. This revised version was published online in August 2006 with corrections to the Cover Date.     

Change in the stress intensity factors for plane cracks in a transversely isotropic medium due to thermal loads

The feasibility of controlling the stress intensity factors for plane cracks of arbitrary form (distributed in a plane perpendicular to isotropy axis in a transversely isotropic material) which are subjected to symmetric mechanical loading by heating the material is demonstrated by using the congruence theorem and making an analogy between isotropic and transversely isotropic materials based on the theory of thermoelasticity. It is shown that a thermal load which fully compensates for the mechanical load can be created within the range in which the material behaves linearly. If it turns out to be technically impossible to create the necessary thermal load, a simpler temperature regime can be chosen that will “cancel” the mechanical force field with a certain factor of safety. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev. Translated from Prikladnaya Mekhanika, Vol. 35, No. 9, pp. 29–37, September, 1999.     

“Ion wind”, a gas-dynamic flow in the corona discharge,and its interaction with the external flow

The “ion wind”, a gas-dynamic flow in the corona discharge that arises owing to transfer of the ion component momentum to the neutral particles of an initially stationary gas, and its interaction with the external flow perpendicular to it are studied. A physico-mathematical model of the flows considered is proposed and the corresponding equations are analyzed numerically. The boundary conditions used for the electric quantities approximately model the conditions in the negative corona discharge between a thin corona-forming electrode and a plane grid electrode transparent to the gas.     

The Exact Theory of Plate Bending

A new approach is introduced for the analysis and calculation of homogeneous, isotropic elastic plates of constant thickness under arbitrary bending loads. This theory can be called “exact” because it leads to exact values of the generalized 2D quantities. Moreover, contrary to classical plate theories, it is not limited to relatively thin plates. This revised version was published online in June 2006 with corrections to the Cover Date.     

On testing of charpy specimens using the one-point bend impact technique

This paper reports our methodology and results for the assessment of the dynamic fracture energy of notched Charpy A508 steel specimens. The fracture tests consist of one-point bend impact applied to the specimen in contact with an instrumented bar. Fracture is caused by the inertia of the unsupported specimen only. The fracture energy is determined from the incident, reflected and single wire fracture gage signals. High-speed photographic recordings show that for all the specimens investigated in the “lower shelf” temperature regime, fracture occurs relatively early and prior to “taking off” of the bar by rigid body motion. It also confirms that the fracture gage readings indeed coincide with the formation of a crack from the notch tip. The present methodology is relatively easy to implement, and it allows the investigation of the fracture properties of materials at loading rates (and velocities) that are substantially higher than those achieved in a conventional Charpy test. Moreover, this test is attractive for modeling purposes since its boundary conditions are simple and well defined.