Out-of-plane loading of anisotropic bimaterials and semi-infinite materials

Summary Analytical closed-form solutions are proposed in a rather compact form for the stress and displacement fields induced by out-of-plane loading of a semi-infinite anisotropic material with inclined strata. The solutions are then extended to include the case of a bimaterial with a planar interface. Several boundary conditions are considered for the interface which may be between two anisotropic half-planes with different elastic properties, or two different orientations of the strata in the same material.     

Green's functions for anisotropic piezoelectric bimaterials with a slipping interface

An explicit full-field expression of the Green's functions for anisotropic piezoelectric bimaterials with a slipping interface is derived. When the electro-elastic singularity reduces to a pure dislocation in displacement and electric potential, interaction energy between the dislocation and the bimaterials is obtained explicitly while the generalized force on the dislocation is given in a real form which is also valid for degenerate materials. The investigation demonstrates that the boundary conditions at the slipping interface between two piezoelectric materials will exert a prominent influence on the mobility of the dislocation. Project supported by the National Natural Science Foundation of China (No. 59635140).     

Generalized stroh formalism for anisotropic elasticity for general boundary conditions

The Stroh formalism is most elegant when the boundary conditions are simple, namely, they are prescribed in terms of traction or displacement. For mixed boundary conditions such as there for a slippery boundary, the concise matrix expressions of the Stroh formalism are destroyed. We present a generalized Stroh formalism which is applicable to a class of general boundary conditions. The general boundary conditions include the simple and slippery boundary conditions as special cases. For Green's functions for the half space, the general solution is applicable to the case when the surface of the half-space is a fixed, a free, a slippery, or other more general boundary. For the Griffith crack in the infinite space, the crack can be a slit-like crack with free surfaces, a rigid line inclusion (which is sometimes called an anticrack), or a rigid line with slippery surface or with other general surface conditions. It is worth mention that the modifications required on the Stroh formalism are minor. The generalized formalism and the final solutions look very similar to those of unmodified version. Yet the results are applicable to a rather wide range of boundary conditions.     

Radiation boundary conditions for finite element solutions of generalized wave equations

On the basis of the dispersion relation of the generalized linear wave equation we derive a radiation boundary condition (RBC) that explicitly incorporates the physical parameters of the governing equation into the form of the boundary condition. Using finite element techniques we investigate the properties of the generalized RBC by examining forced and unforced solutions to the telegraph and Klein-Gordon equations in one dimension. The results show that within the limits of the physical parameters of the problem the generalized RBC is an improvement over the Sommerfeld RBC when the governing equation contains additional terms that influence the propagation. These gains are achieved without introducing any computational overhead. A two-dimensional example suggests that the 1D findings can generalize to higher dimensions.     

The boundary layer solutions of the interface problem considering the strain gradient

According to the classical elastic theory, there is always a discontinuity of rotation angle on the interface between different materials. This illogic result can be overcome by the strain gradient plasticity theory. In the light of this theory, there is a group of boundary layer solutions near the interface, which have made important adjustment of the classical results. This project is supported by National Natural Science Foundation of China (19891180).     

Combined interface boundary condition method for coupled thermal simulations

A new procedure for modeling the conjugate heat‐transfer process between fluid and structure subdomains is presented. The procedure relies on higher‐order combined interface boundary conditions (CIBC) for improved accuracy and stability. Traditionally, continuity of temperature and heat flux along interfaces is satisfied through algebraic jump conditions in a staggered fashion. More specifically, Dirichlet temperature conditions are usually imposed on the fluid side and Neumann heat‐flux conditions are imposed on the solid side for the stability of conventional sequential staggered procedure. In this type of treatment, the interface introduces additional stability constraints to the coupled thermal simulations. By utilizing the CIBC technique on the Dirichlet boundary conditions, a staggered procedure can be constructed with the same order of accuracy and stability as those of standalone computations. Using the Godunov–Ryabenkii normal‐mode analysis, a range of values of the coupling parameter is found that yields a stable and accurate interface discretization. The effectiveness of the method is investigated by presenting and discussing performance evaluation data using a 1D finite‐difference formulation for each subdomain. Copyright © 2007 John Wiley & Sons, Ltd.     

Contact zone assessment for a fast growing interface crack in an anisotropic bimaterial

Summary A plane strain problem for a crack with a frictionless contact zone at the leading crack tip expanding stationary along the interface of two anisotropic half-spaces with a subsonic speed under the action of various loadings is considered. The cases of finite and infinite-length interface cracks under the action of a moving concentrated loading at its faces are considered. A finite-length crack for a uniform mixed-mode loading at infinity is considered as well. The associated combined Dirichlet-Riemann boundary value problems are formulated and solved exactly for all above-mentioned cases. The expressions for stresses and the derivatives of the displacement jumps at the interface are presented in a closed analytical form for an arbitrary contact zone length. Transcendental equations are obtained for the determination of the real contact zone length, and the associated closed form asymptotic formulas are found for small values of this parameter. It is found that independently of the types of the crack and loading, an increase of the crack tip speed leads to an increase of the real contact zone length and the correspondent stress intensity factor. The latter increase significantly for an interface crack tip speed approaching the Ragleigh wave speed.     

Analytical solutions for thermal forcing vortices in boundary layer and its applications

Using the Boussinesq approximation, the vortex in the boundary layer is assumed to be axisymmetrical and thermal-wind balanced system forced by diabatic heating and friction, and is solved as an initial-value problem of linearized vortex equation set in cylindrical coordinates. The impacts of thermal forcing on the flow field structure of vortex are analyzed. It is found that thermal forcing has significant impacts on the flow field structure, and the material representative forms of these impacts are closely related to the radial distribution of heating. The discussion for the analytical solutions for the vortex in the boundary layer can explain some main structures of the vortex over the Tibetan Plateau.     

Comparison of outflow boundary conditions for subsonic aeroacoustic simulations

Aeroacoustics simulations require much more precise boundary conditions than classical aerodynamics. Two classes of non‐reflecting boundary conditions for aeroacoustics are compared in the present work: the characteristic analysis‐based methods and the Tam and Dong approach. In the characteristic methods, waves are identified and manipulated at the boundaries, whereas the Tam and Dong approach use modified linearized Euler equations in a buffer zone near outlets to mimic a non‐reflecting boundary. The principles of both approaches are recalled, and recent characteristic methods incorporating the treatment of transverse terms are discussed. Three characteristic techniques—the original Navier–Stokes characteristic boundary conditions (NSCBC) of Poinsot and Lele and two versions of the modified method of Yoo and Im—are compared with the Tam and Dong method for four typical aeroacoustics problems: vortex convection on a uniform flow, vortex convection on a shear flow, acoustic propagation from a monopole, and acoustic propagation from a dipole. Results demonstrate that the Tam and Dong method generally provides the best results and is a serious alternative solution to characteristic methods even though its implementation might require more care than the usual NSCBC approaches. Copyright © 2011 John Wiley & Sons, Ltd.     

Rayleigh waves with impedance boundary conditions in anisotropic solids

The paper is concerned with the propagation of Rayleigh waves in an elastic half-space with impedance boundary conditions. The half-space is assumed to be orthotropic and monoclinic with the symmetry plane _x3=0$x_3=0$. The main aim of the paper is to derive explicit secular equations of the wave. For the orthotropic case, the secular equation is obtained by employing the traditional approach. It is an irrational equation. From this equation, a new version of the secular equation for isotropic materials is derived. For the monoclinic case, the method of polarization vector is used for deriving the secular equation and it is an algebraic equation of eighth-order. When the impedance parameters vanish, this equation coincides with the secular equation of Rayleigh waves with traction-free boundary conditions.     

经典梁热过屈曲问题的解析解

基于非线性经典梁理论,建立了控制轴向和横向变形的基本方程,将两个非线性方程化简为一个关于横向挠度的四阶非线性积分-微分方程。对于本文所考虑的三类边界条件,该方程与相应的边界条件构成了微分特征值问题;直接求解该问题,得到热过屈曲构形的解析解,该解是外加热载荷的函数。为考察热载荷以及边界条件的影响,根据得到的解析解给出了一些数值算例,讨论了梁过屈曲行为的性质。本文得到的解析解可用于验证或改进各类近似理论和数值方法。     

Exact solutions for variable-thickness inhomogeneous elastic plates under various boundary conditions

In this paper, an exact solution to the governing equations of the bending of a variable-thickness inhomogeneous rectangular plate is presented. The procedure is applicable to variable-thickness inhomogeneous rectangular plates with two opposite edges simply supported. The remaining ones subjected to a combination of clamped, simply supported, and free boundary conditions and between these two edges the plate may have varying thickness. The procedure is valuable in view of the fact that tables of deflections and stresses cannot be presented for variable-thickness inhomogeneous orthotropic plates as for uniform-thickness homogeneous isotropic plates even for commonly encountered loads because the results depend on the inhomogeneity coefficient and the orthotropic material properties instead of a single flexural rigidity. Numerical results, useful for the validation or otherwise of approximate solutions, are tabulated. The influences of the degree of the inhomogeneity, aspect ratio, thickness parameter and degree of non-uniformity on the deflections and stresses are investigated. This paper is partially supported by the Deanship of Scientific Research at King AbdulAziz University (Grant no. 172/427).     

On the boundary conditions at the macroscopic level

We study the problem of the boundary conditions specified at the boundary of a porous domain in order to solve the macroscopic transfer equations obtained by means of the volume-averaging method. The analysis is limited to the case of conductive transport but the method can be extended to other cases. A numerical study enables us to illustrate the theoretical results in the case of a model porous medium. Roman Letters _sf interfacial area of the s-f interface contained within the macroscopic system m² - A _sf interfacial area of the s-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - d _s , d _f microscopic characteristic length m - g vector that maps to _s, m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l REV characteristic length, m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for the f-phase at the f-s interface - n _e outwardly directed unit normal vector at the dividing surface - T ^* macroscopic temperature field obtained by solving the macroscopic equation (3), K - V averaging volume, m³ - V _s , V _f volume of the considered phase within the averaging volume, m³ - volume of the macroscopic system, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² Greek Letters _s , _f volume fraction - ratio of thermal conductivities - _s , _f thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding to T ^* , K - spatial deviation temperature K - error on the temperature due to the macroscopic boundary conditions, K - spatial average - ^s , ^f intrinsic phase average     

Analytical solutions for non-Newtonian fluid flows in pipe-like domains

This paper deals with some unsteady unidirectional transient flows of an Oldroyd-B fluid in unbounded domains which geometrically are axisymmetric pipe-like. An expansion theorem of Steklov is used to obtain exact solutions for flows satisfying no-slip boundary conditions. The well known solutions for a Navier-Stokes fluid, as well as those corresponding to a Maxwell fluid and a second grade one, appear as limiting cases of our solutions. The steady state solutions are also obtained for t→∞.     

Numerical boundary conditions for unsteady transonic flow calculations

In calculations of transonic flows it is necessary to limit the domain of computation to a size that is manageable by computers. At the boundary of the computational domain, boundary conditions are required to ensure a unique solution. Since wave solutions exist in the unsteady transonic flow field, incorrect boundary conditions may result in spurious reflections from the computational boundary. This may introduce errors into the solution. To prevent the spurious reflections, absorbing boundary conditions are often used on the computational boundary. In this paper we describe a method to derive absorbing boudary conditions for transonic calculations. We demonstrate both theoretically and numerically that the use of the absorbing boundary conditions will reduce the spurious reflections in the calculation.     

Analytical and numerical solutions of the Local Inertial Equations

Neglecting the convective terms in the Saint-Venant Equations (SVE) in flood hydrodynamic modelling can be done without a loss in accuracy of the simulation results. In this case the Local Inertial Equations (LInE) are obtained. Herein we present two analytical solutions for the Local Inertial Equations. The first is the classical instantaneous Dam-Break Problem and the second a steady state solution over a bump. These solutions are compared with two numerical schemes, namely the first order Roe scheme and the second order MacCormack scheme. Comparison between analytical and numerical results shows that the numerical schemes and the analytical solution converge to a unique solution. Furthermore, by neglecting the convective terms the original numerical schemes remain stable without the need for adding entropy correction, artificial viscosity or special initial conditions, as in the case of the full SVE.     

Analytical solutions for elastic binary nanotubes of arbitrary chirality

Analytical solutions for the elastic properties of a variety of binary nanotubes with arbitrary chirality are obtained through the study of systematic molecular mechanics. This molecular mechanics model is first extended to chiral binary nanotubes by introducing an additional outof-plane inversion term into the so-called stick-spiral model,which results from the polar bonds and the buckling of binary graphitic crystals. The closed-form expressions for the longitudinal and circumferential Young's modulus and Poisson's ratio of chiral binary nanotubes are derived as functions of the tube diameter. The obtained inversion force constants are negative for all types of binary nanotubes, and the predicted tube stiffness is lower than that by the former stick-spiral model without consideration of the inversion term, reflecting the softening effect of the buckling on the elastic properties of binary nanotubes. The obtained properties are shown to be comparable to available density functional theory calculated results and to be chirality and size sensitive. The developed model and explicit solutions provide a systematic understanding of the mechanical performance of binary nanotubes consisting of Ⅲ–Ⅴ and Ⅱ–Ⅵ group elements.     

Some refinements concerning the boundary conditions at the macroscopic level

The study of boundary effects initiated in a previous paper is continued. New assumptions regarding the geometrical structure of the boundary surface are introduced. Under these assumptions, it is shown that macroscopic Neumann conditions do not generally affect the determination of the macroscopic field in the case of the transport process considered — heat conduction. For this type of boundary condition, the boundary effect is generally confined within a thin layer near the boundary. When heat sources are taken into account within the porous domain, the result is different. In this case, making use of a Neumann boundary condition, expressed in terms of macroscopic variables, amounts to introducing an extra flux. Under normal circumstances, however, this additional flux is negligible.Roman Letters A cross-sectional area of a unit cell - A _e cross-sectional area of a unit cell at the boundary surface - A _sf interfacial area of the s-f interface contained within the averaging volume - surface area per unit volume (A _sf/ ) - A _sf interfacial area of the s-f interface contained within the macroscopic system - g closure vector - h closure vector - k heat transfer coefficient at the s-f interface - K_eff effective thermal conductivity tensor - _x unit cell length - n unit vector - n_e outwardly directed unit normal vector at the boundary - n_sf outwardly directed unit normal vector for thes-phase at f-s interface - q heat flux density - T ^* macroscopic temperature defined by the macroscopic problem - s closure variable - V volume of the macroscopic system - V boundary surface of the macroscopic domain - V ₁ macroscopic sub-surface of the boundary surface - x local coordinate Greek Letters _s,_f volume fraction - _s, gl_f microscopic thermal conductivities - true microscopic temperature - ^* microscopic temperature corresponding toT ^* - microscopic error temperature - vector defined by Equation (34) - < > spatial average     

Building generalized open boundary conditions for fluid dynamics problems

This paper deals with the design of an efficient open boundary condition (OBC) for fluid dynamics problems. Such problematics arise, for instance, when one solves a local model on a fine grid that is nested in a coarser one of greater extent. Usually, the local solution U^loc is computed from the coarse solution U^ext, thanks to an OBC formulated as , where B_h and B_H are discretizations of the same differential operator (B_h being defined on the fine grid and B_H on the coarse grid). In this paper, we show that such an OBC cannot lead to the exact solution, and we propose a generalized formulation , where g is a correction term. When B_h and B_H are discretizations of a transparent operator, g can be computed analytically, at least for simple equations. Otherwise, we propose to approximate g by a Richardson extrapolation procedure. Numerical test cases on a 1D Laplace equation and on a 1D shallow water system illustrate the improved efficiency of such a generalized OBC compared with usual ones. Copyright © 2012 John Wiley & Sons, Ltd.     

Analytical solutions of dynamic mode I crack under the conditions of displacement boundary

By the approaches of the theory of complex variable functions, the problems of dynamic mode I crack under the condition of displacement boundary are investigated. For this kind of dynamic crack extension problems with arbitrary index of self-similarity, the universal representations of analytical solutions are facilely deduced by the methods of self-similar functions. Analytical solutions of the stresses, displacements and stress intensity factors are readily acquired using the methods of self-similar functions. The problems studied can be very easily translated into Riemann–Hilbert problems and their closed solutions are gained rather straightforward in terms of this technique. According to corresponding material properties, the mutative rule of stress intensity factor was illustrated very well. Using those solutions and superposition theorem, the solutions of arbitrarily complex problems can be attained.