轴对称一维六方准晶圆柱的精化理论

将轴对称圆柱的精化分析推广到一维六方准晶中轴对称圆柱的研究当中。利用准调和函数的Bessel算子函数表示以及一维六方准晶中的通解,在不做任何预先假设的情况下,给出了一维六方准晶中轴对称圆柱的精化理论。首先,根据准调和函数的Bessel算子函数表示,利用三个一维待定函数,表示出声子场和相位子场的位移场和应力分量。再根据非齐次边界条件,推导出柱面受径向外载时的精化方程。通过舍弃高阶项,推导出了在径向方向受到柱面载荷的近似解。     

The general solution of one-dimensional hexagonal quasicrystal

A general solution is given for three-dimensional dynamic problems in 1D hexagonal quasicrystalline material with Laue class 6/m_hmm. In the general solution, the phonon and phason displacements and stresses are controlled by two displacement functions Ψ and F, which satisfy a second-order partial differential equation, and a sixth-order partial differential equation, respectively. The general solution for static equilibrium problems is treated as a special case.     

Locally-exact homogenization theory for transversely isotropic unidirectional composites

The locally-exact homogenization theory for unidirectional composites with square periodicity and isotropic phases proposed by Drago and Pindera [18] is extended to architectures with hexagonal symmetry and transversely isotropic phases. The theory employs Fourier series representation for the displacement fields in the fiber and matrix phases in the cylindrical coordinate system that satisfies exactly the equilibrium equations and continuity conditions in the unit cell's interior. The inseparable exterior problem involves satisfaction of periodicity conditions for the hexagonal unit cell geometry demonstrated herein to be readily achievable using the previously introduced balanced variational principle for square geometries. This variational principle plays a key role in the employed unit cell solution, ensuring rapid convergence of the Fourier series coefficients with relatively few harmonic terms, yielding converged homogenized moduli and local stress fields with little computational effort. The solution's stability is illustrated using the dilute case which is shown to reduce to the Eshelby solution regardless of the employed number of harmonic terms. Comparison with published results and predictions of a finite-volume based homogenization in a wide fiber volume range and different fiber/matrix modulus contrast validates the approach's accuracy, and its utility is demonstrated through rapid local stress recovery in a multi-scale application. This extension completes the development of the theory for three important classes of unidirectional reinforcement arrays, thereby providing an efficient alternative to finite-element based homogenization techniques or approximate micromechanical schemes, as well as an efficient standard against which other methods may be compared.     

CRACK TIP PLASTICITY OF A THERMALLY LOADED PENNY-SHAPED CRACK IN AN INFINITE SPACE OF 1D QC

The present work is concerned with a penny-shaped Dugdale crack embedded in an infinite space of one-dimensional(1D) hexagonal quasicrystals and subjected to two identical axisymmetric temperature loadings on the upper and lower crack surfaces. Applying Dugdale hypothesis to thermo-elastic results, the extent of the plastic zone at the crack tip is determined.The normal stress outside the plastic zone and crack surface displacement are derived in terms of special functions. For a uniform loading case, the corresponding results are presented by simplifying the preceding results. Numerical calculations are carried out to show the influence of some parameters.     

An Approximate Analytical Solution for Counter-Current Spontaneous Imbibition

An approximate analytical solution is provided for one-dimensional, counter- current, spontaneous imbibition of a wetting phase (water) into a semi-infinite porous medium. The solution is based on the assumption that a similarity solution exists for the displacement process. This assumption, in turn, rests on the assumption that the set of relative permeability and capillary pressures curves are unique functions of saturation and do not depend on the nature of the displacement. It further rests on the assumption that the saturation at the imbibition face does not vary with time. It is demonstrated that the solution is in agreement with results obtained from experiments and also numerical analyses of these experiments. The experiments utilize cylindrical samples with the radial surface and one end-face sealed, and with counter-current imbibition occurring at the open end-face. The stage of the experiment that is modeled by the present solution is the period before the imbibition front contacts the sealed end-face. An important finding of the present analysis is that the pressure upstream of the advancing invasion front is a constant. A second, improved solution is also presented; this solution is an iterative, series solution of an integral-differential equation. It converges to a stable solution in very few terms.     

Solvability on boundary-value problems of elasticity of three-dimensional quasicrystals

Weak solution (or generalized solution) for the boundary-value problems of partial differential equations of elasticity of 3D (three-dimensional) quasicrystals is given, in which the matrix expression is used. In terms of Korn inequality and theory of function space,we prove the uniqueness of the weak solution.This gives an extension of existence theorem of solution for classical elasticity to that of quasicrystals,and develops the weak solution theory of elasticity of 2D quasicrystals given by the second author of the paper and his students.     

Characteristic equation solution strategy for deriving fundamental analytical solutions of 3D isotropic elasticity

A simple characteristic equation solution strategy for deriving the fundamental analytical solutions of 3D isotropic elasticity is proposed. By calculating the determinant of the differential operator matrix obtained from the governing equations of 3D elasticity, the characteristic equation which the characteristic general solution vectors must satisfy is established. Then, by substitution of the characteristic general solution vectors, which satisfy various reduced characteristic equations, into various reduced adjoint matrices of the differential operator matrix, the corresponding fundamental analytical solutions for isotropic 3D elasticity, including Boussinesq-Galerkin (B-G) solutions, modified Papkovich-Neuber solutions proposed by Min-zhong WANG (P-N-W), and quasi HU Hai-chang solutions, can be obtained. Furthermore, the independence characters of various fundamental solutions in polynomial form are also discussed in detail. These works provide a basis for constructing complete and independent analytical trial functions used in numerical methods.     

Transverse vibrations of arbitrary non-uniform beams

Free and steady state forced transverse vibrations of non-uniform beams are investigated with a proposed method, leading to a series solution. The obtained series is verified to be convergent and linearly independent in a convergence test and by the non-zero value of the corresponding Wronski determinant, respectively. The obtained solution is rigorous, which can be reduced to a classical solution for uniform beams. The proposed method can deal with arbitrary non-uniform Euler-Bernoulli beams in principle, but the methods in terms of special functions or elementary functions can only work in some special cases.     

Boundary conditions for axisymmetric circular cylinder in one-dimensional hexagonal quasicrystals

Through generalizing the method of a decay analysis technique determining the interior solution developed by Gregory and Wan, a set of necessary conditions on the end-data of axisymmetric circular cylinder in one-dimensional (1D) hexagonal quasicrystals (QCs) for the existence of a rapidly decaying solution is established. By accurate solutions for auxiliary regular state, and using the reciprocal theorem and general solution, these necessary conditions for the end-data to induce only a decaying elastostatic state (boundary layer solution) will be translated into appropriate boundary conditions for the circular cylinder with axisymmetric deformations in 1D hexagonal QCs. The results of the present paper enable us to establish a set of correct boundary conditions, most of which are obtained for the first time. Furthermore, corresponding necessary conditions for the case of axisymmetric circular cylinder with transversely isotropic are derived subsequently, and their isotropic elastic counterparts are also obtained. The accuracy of mixed boundary conditions of isotropic axisymmertic circular cylinder is proved at last.     

Elastic Anisotropic Material with Purely Longitudinal and Transverse Waves

The simplest form of the matrix of elasticity moduli of an anisotropic material conducting purely longitudinal and transverse waves with an arbitrary direction of the wave normal is obtained. A generic solution of equations in displacements is represented in terms of three functions satisfying independent wave equations. In the case of planar deformation, this solution yields a complex representation coinciding with the Kolosov–Muskhelishvili formulas for an isotropic material. The formulas in the present work also determine an anisotropic material with Young's modulus identical for all directions, as in an isotropic medium.     

Three-dimensional analytical solution for a rotating disc of functionally graded materials with transverse isotropy

Based on the basic equations for axisymmetric problems of transversely isotropic elastic materials, the displacement components are expressed in terms of polynomials of the radial coordinate with the five involved coefficients, named as displacement functions in this paper, being undetermined functions of the axial (thickness) coordinate. Five equations governing the displacement functions are then derived. It is shown that the displacement functions can be found through progressive integration by incorporating the boundary conditions. Thus a three-dimensional analytical solution is obtained for a transversely isotropic functionally graded disc rotating at a constant angular velocity.The solution can be degenerated into that for an isotropic functionally graded rotating disc. A prominent feature of this solution is that the material properties can be arbitrary functions of the axial coordinate. Thus, the solution for a homogeneous transversely isotropic rotating disc is just a special case that can be easily derived. An example is finally considered for a special functionally graded material, and numerical results shows that the material inhomogeneity has a remarkable effect on the elastic field.     

弹性力学广义平面应力的精化理论

从Papkovich-Neuber通解出发，导出了弹性力学广义平面应力问题的精化理论控制方程．若取零阶近似表示，方程即为通常弹性力学平面应力问题的平衡方程．对略去h^2及以上高阶项的八阶方程给出了通解．     

A two-step back–forth algorithm for the solution of the unsteady adjoint equations

A new algorithm for the solution of the unsteady adjoint equations is proposed in this article, aiming at overcoming the excessive computational cost and memory requirements of the conventional adjoint approach for the optimisation of unsteady problems in computational mechanics. The total cost is equal to four times the cost of the unsteady state solution, which is twice the cost of the conventional backward-in-time adjoint calculation but the memory requirements are very small, equivalent to those of a steady-state problem, while stability is acceptable. The proposed algorithm is validated in the case of the 1D unsteady Burgers equation with non-smooth source terms.     

Some exact expressions for the temporal evolution of long Rossby waves on a beta-plane

Some exact expressions are derived to describe the temporal evolution of forced Rossby waves in a two-dimensional beta-plane configuration where the background flow has constant zonal-mean velocity. The meridional length scale of the problem is assumed to be small relative to the zonal length scale and so the long-wave limit of zero aspect ratio is taken. In the case where the background flow velocity is zero, an exact solution is obtained in terms of generalized hypergeometric functions. A late-time asymptotic approximation is obtained and it shows that the solution oscillates with time and its amplitude goes to zero in the limit of infinite time. In the case of a non-zero background flow velocity, the solution is evaluated using two different procedures which give two equivalent expressions in terms of different generalized hypergeometric functions. The late-time asymptotic behaviour is investigated and it is found that the solution approaches a steady state in the limit of infinite time.We also derive a solution in the form of an asymptotic series expansion for the more general situation where a Rossby wave packet is generated by a zonally-localized boundary condition comprising a continuous spectrum of wavenumbers or Fourier modes. The exact solutions found here can be used as leading-order solutions in weakly-nonlinear analyses and other studies involving more realistic configurations for time-dependent Rossby waves or wave packets.     

横观各向同性三维热弹性力学通解及其势理论法

通过引入两个位移函数，对用位移表达的运动平衡方程作了简化．利用算子理论，严格地导出了横观各向同性非耦合热弹性动力学问题的通解．对于静力学问题，通解的形式可进一步简化成用4个准调和函数来表示．具体考察了横观各向同性体内平面裂纹上下表面有对称分布温度作用的问题，推广了势理论方法，导出了一个积分方程和一个微分-积分方程．针对币状裂纹表面受均布温度作用情形，给出了具体的解。     

Several three-dimensional solutions for transversely isotropic functionally graded material plate welded with circular inclusion

The problem of a transversely isotropic functionally graded material (FGM) plate welded with a circular inclusion is considered. The analysis starts with the generalized England-Spencer plate theory for transversely isotropic FGM plates, which expresses a three-dimensional (3D) general solution in terms of four analytic functions. Several analytical solutions are then obtained for an infinite FGM plate welded with a circular inclusion and subjected to the loads at infinity. Three different cases are considered, i.e., a rigid circular inclusion fixed in the space, a rigid circular inclusion rotating about the x-, y-, and z-axes, and an elastic circular inclusion with different material constants from the plate itself. The static responses of the plate and/or the inclusion are investigated through numerical examples.     

一类指数矩阵函数及其应用

研究了一阶常微分方程组特解的精细积分方法. 针对非齐次项为多项式、指数函数以及二者的乘积的情况,在Duhamel积分形式特解的基础上,引入了一类指数矩阵函数. 通过该类函数的线性组合即可表达出非齐次方程的特解. 建立了该类指数矩阵函数的一种高效递推算法,并在此基础上实现了特解的精细积分. 由于特解的积分过程能充分利用通解精细积分过程的中间量,因此两个精细积分过程能有机地结合起来,形成了一种高效、统一的广义精细积分法. 对上述递推算法做了进一步优化,并给出了通用的计算公式.算例结果证明了该方法的有效性.      

Boundary Conditions for Plate Bending in One-dimensional Hexagonal Quasicrystals

For plate bending in one-dimensional (1D) hexagonal quasicrystals (QCs), the reciprocal theorem and the general solution for QCs media are applied in a novel way to obtain the appropriate stress and mixed boundary conditions accurate to all order for plates of general edge geometry and loadings. Through generalizing the method developed by Gregory and Wan, a set of necessary conditions on the edge-data for the existence of a rapidly decaying solution is established. The prescribed data must satisfy these conditions in order that they should generate a decaying state. When a set of stress edge-data or mixed edge-data is imposed on the plate edge, these decaying state conditions for the case of axisymmetric deformation of 1D hexagonal QC plates are derived explicitly. They are then used for the correct formulation of boundary conditions for the plate theory solution (or the interior solution). Furthermore, in the absence of phonon–phason fields coupling effect, corresponding necessary conditions for the case of transversely isotropic elastic plates are derived subsequently, and their isotropic elastic counterparts are also obtained.      

A general solution to some plane problems of micropolar elasticity

We obtain a general solution to the field equations of plane micropolar elasticity for materials characterized by a hexagonal or equilateral triangular structure. These materials exhibit 3-fold symmetry in the plane and the elastic response is isotropic. Utilizing two displacement potential functions, the solution is obtained in terms of two analytic functions and a third function satisfying the modified homogeneous Helmholtz equation. Expressions for the two-dimensional components of displacement, stress, and couple stress, along with the resultant force on a contour, are presented. We observe that micropolar effects are most significant in material regions subjected to large deformation gradients. Specific results are presented for the classical crack problem, the half plane loaded uniformly on the surface, Flamant's problem, and the circular cylinder compressed by equal and opposite concentrated forces.     

二维复合材料结构模态分析的边界元法

将边界元方法用于分析二维复合材料结构的自由振动模态，利用特解处理体积力（惯性力）仅需静态基本解就可求解问题，对一各向同性悬臂梁，用该法得到的结果与用有限元或各向性边界元法得到的结果符合得很好，但该法可解各向异性问题，对层状复合材料简支梁，用该法得到了数值结果与用一维层状复合材料梁的理论解的比较表明，当结构的长厚比大于２０时，二者符合得很好，当结构的长厚比小于２０时，一维层状材料梁的理论将产生很大的