插值矩阵法分析正交各向异性板切口应力奇异性

基于切口尖端附近区域位移场渐近展开,提出了分析正交各向异性复合材料板切口奇异性的新方法.将位移场的渐近展开式的典型项代入弹性板的基本方程,得到关于正交各向异性板切口奇异性指数的一组非线性常微分方程的特征值问题;再采用变量代换法,将非线性特征问题转化为线性特征问题,用插值矩阵法求解获得的正交各向异性板切口若干阶应力奇异性指数和相应特征函数.该法可由相应的特征角函数对板切口的平面应力和反平面奇异特征值加以区分,并将计算结果与现有结果对照,表明了该文方法的有效性.     

双材料界面裂纹平面问题的半权函数法

应用半权函数法求解双材料界面裂纹的平面问题．由平衡方程、应力应变关系、界面的连续条件以及裂纹面零应力条件推导出裂尖的位移和应力场,其特征值为lambda及其共轭．设置特征值为lambda的虚拟位移和应力场,即界面裂纹的半权函数A·D2由功的互等定理得到应力强度因子KⅠ和KⅡ以半权函数与绕裂尖围道上参考位移和应力积分关系的表达式．数值算例体现了半权函数法精度可靠、计算简便的特点．     

裂纹与弹性夹杂的相互影响*

本文利用无限域上单根弹性夹杂和单根裂纹产生的位移和应力,将裂纹与弹性夹杂的相互影响问题归为解一组柯西型奇异积分方程,然后用此对夹杂分枝裂纹解答的奇性性态作了理论分析,并求得了振荡奇性界面应力场,对于不相交的夹杂裂纹问题,具体计算了端点的应力强度因子及夹杂上的界面应力,结果令人满意。     

正交异性双材料反平面界面端应力场分析

研究了正交异性双材料反平面平板搭接界面端问题,采用复合材料断裂复变方法,构造了特殊应力函数,通过求解一类广义重调和方程组的边值问题,推导出平板搭接界面端的应力场、位移场及应力强度因子的表达式,结果显示:反平面搭接界面端只有一个奇异性,上下材料常数比Γ＞0时,应力场具有幂次奇异性,且随着Γ增长,奇异指数趋于-1/2,并利用有限元算例分析验证了理论结果的正确性.     

深埋公路隧洞围岩应力和位移分布的复变函数解

为了优化公路隧道的设计和确保施工安全,必须明确公路隧洞开挖时围岩的力学行为.利用复变函数方法,通过保角映射函数把隧洞外域变换为单位圆外域.利用Cauchy(柯西)积分和留数定理求出两个应力函数,从而得到围岩的应力与位移的平面应变问题的解析解.结合曲墙马蹄形断面,通过数学软件MATLAB编程计算,分别给出了应力和仅考虑开挖引起的位移沿隧洞边和坐标轴方向的分布.利用有限元软件ANSYS建立二维平面应变模型,对理论推导得到的应力和位移的分布进行验证,数值解结果与近似解析解结果吻合性很好.研究结果表明:最大的环向应力发生在隧洞拱脚处,最大水平位移发生在拱腰处,最大的沉降和隆起分别发生在拱顶和仰拱中心处.沿坐标轴的正应力在隧洞附近变化较大,不一定在洞边取得最大值,离洞边不到10 m的距离,便分别趋于所加外荷载.位移值在洞边最大,随着离洞边距离的增大,逐渐单调趋于0.     

辛求解体系下功能梯度材料平面梁的完整解析解

采用辛弹性力学解法,求取弹性模量沿轴向指数变化,而Poisson比保持不变的功能梯度材料平面梁的完整解析解．通过求解被Saint-Venant原理覆盖的一般本征解,建立起完整的解析分析过程,进而给出平面梁位移和应力的精确分布规律．传统的弹性力学分析方法常常忽略被Saint-Venant原理覆盖的解,但这些衰减的本征解对材料的局部效应起着较大的影响作用,可能导致材料或结构的突然失效．采用辛求解方法,充分利用本征向量之间的辛共轭正交关系,得到了功能梯度材料梁的完整解析解．两个数值算例分别将功能梯度材料平面梁的位移和应力分布与相应均匀材料情形的结果进行比较,研究了材料非均匀性对位移和应力解的影响．     

引用复变量伪应力函数来解幂硬化材料平面应力问题

本文引用复变量伪应力函数将幂硬化材料平面应力问题的协调方程化为双调和方程,从而使此类有强化材料的弹塑性平面应力问题能像线弹性力学平面问题那样采用复变函数法进行求解.本文推导出了幂硬化材料平面应力问题的应力、应变及位移分量的复变函数表达式,可推广应用于满足全量理论的一股弹塑性平面应力问题.作为算例,文中给出了含圆孔幂硬化材料无限大板单向受拉问题的解答,并和有关文献用摄动法获得的同一问题的渐近解进行了比较.     

含非局部边界条件的奇异特征值问题的正解

本文考虑奇异特征值问题其中μ0,p∈(1/2,1]和λ[v]=∫_0~1v(t)dA(t)是C[0,1]上由Riemann-Stieltjes积分定义的一个线性泛函;函数g∈C(0,1)在t=0和/或t=1处可能有奇性,f在u=0处有奇性.本文首先研究Green函数的性质和先验估计,以及利用Krein-Rutman定理建立了线性算子第一特征值,最后联合不动点定理证明了特征值问题正解的存在性,同时给出了参数μ的取值区间.     

单裂纹问题的虚边界无网格伽辽金法分析

使用含裂纹复变基本解,虚边界无网格伽辽金法被进一步推广应用于弹性材料的单裂纹问题求解.为了清晰地说明单裂纹问题的虚边界元法实现过程,单裂纹问题的虚边界元法示意图、复变坐标平面下含裂纹问题的复变位移和复变面力基本解示意图被展示.含裂纹复变基本解,因自动满足裂纹处边界条件,故使用虚边界无网格伽辽金法计算单裂纹问题,无需在裂纹处布置节点或单元.给出含裂纹复变基本解中的Φ'(x)的详细表达式、裂纹左右裂尖应力强度因子的虚边界无网格离散公式,方便了其他学者使用本方法计算裂纹问题.数值计算两端受拉长方形钢板中心含有裂纹的应力强度因子的算例,计算结果证明了本方法的精确性与稳定性.     

指数型功能梯度材料平面问题热应力通解

研究了功能梯度材料平面问题的热应力场,首先引入热弹性位移势函数,得到温度场的应力解;然后引入Airy应力函数,通过求解功能梯度材料平面问题的基本方程,得到不考虑温度时的应力,叠加后得到平面问题的热应力通解.     

Evaluation of the stress singularities of plane V-notches in bonded dissimilar materials

According to the linear theory of elasticity, there exists a combination of different orders of stress singularity at a V-notch tip of bonded dissimilar materials. The singularity reflects a strong stress concentration near the sharp V-notches. In this paper, a new way is proposed in order to determine the orders of singularity for two-dimensional V-notch problems. Firstly, on the basis of an asymptotic stress field in terms of radial coordinates at the V-notch tip, the governing equations of the elastic theory are transformed into an eigenvalue problem of ordinary differential equations (ODEs) with respect to the circumferential coordinate θ around the notch tip. Then the interpolating matrix method established by the first author is further developed to solve the general eigenvalue problem. Hence, the singularity orders of the V-notch problem are determined through solving the corresponding ODEs by means of the interpolating matrix method. Meanwhile, the associated eigenvectors of the displacement and stress fields near the V-notches are also obtained. These functions are essential in calculating the amplitude of the stress field described as generalized stress intensity factors of the V-notches. The present method is also available to deal with the plane V-notch problems in bonded orthotropic multi-material. Finally, numerical examples are presented to illustrate the accuracy and the effectiveness of the method.     

Numerical solution of high-order differential equations by using periodized Shannon wavelets

In this paper, periodized Shannon wavelets are applied as basis functions in solution of the high-order ordinary differential equations and eigenvalue problem. The first periodized Shannon wavelets are defined. The second the connection coefficients of periodized Shannon wavelets are related by a simple variable transformation to the Cattani connection coefficients. Finally, collocation method is used for solving the high-order ordinary differential equations and eigenvalue problem. Some equations are solved in order to find out advantage of such choice of the basis functions.     

A coupled finite element-differential quadrature element method and its accuracy for moving load problem

In this paper a mixed method, which combines the finite element method and the differential quadrature element method (DQEM), is presented for solving the time dependent problems. In this study, the finite element method is first used to discretize the spatial domain. The DQEM is then employed as a step-by-step DQM in time domain to solve the resulting initial value problem. The resulting algebraic equations can be solved by either direct or iterative methods. Two general formulations using the DQM are also presented for solving a system of linear second-order ordinary differential equations in time. The application of the formulation is then shown by solving a sample moving load problem. Numerical results show that the present mixed method is very efficient and reliable.     

On the index of the boundary value problem for a homogeneous Hamiltonian system of differential equations

The index of the homogeneous self-adjoint boundary value problem for the Hamiltonian systems of ordinary differential equations is introduced. It is assumed that the system has a nontrivial solution. The relationship between the index of an eigenvalue of the nonlinear eigenvalue problem and the index of the corresponding homogeneous problem is established. Properties of the index of the problem and those of the eigenvalue are examined.     

Properties of the solutions of the fourth-order Bessel-type differential equation

The structured Bessel-type functions of arbitrary even-order were introduced by Everitt and Markett in 1994; these functions satisfy linear ordinary differential equations of the same even-order. The differential equations have analytic coefficients and are defined on the whole complex plane with a regular singularity at the origin and an irregular singularity at the point of infinity. They are all natural extensions of the classical second-order Bessel differential equation. Further these differential equations have real-valued coefficients on the positive real half-line of the plane, and can be written in Lagrange symmetric (formally self-adjoint) form. In the fourth-order case, the Lagrange symmetric differential expression generates self-adjoint unbounded operators in certain Hilbert function spaces. These results are recorded in many of the papers here given as references. It is shown in the original paper of 1994 that in this fourth-order case one solution exists which can be represented in terms of the classical Bessel functions of order 0 and 1. The existence of this solution, further aided by computer programs in Maple, led to the existence of a linearly independent basis of solutions of the differential equation. In this paper a new proof of the existence of this solution base is given, on using the advanced theory of special functions in the complex plane. The methods lead to the development of analytical properties of these solutions, in particular the series expansions of all solutions at the regular singularity at the origin of the complex plane.     

Accurate stress analysis of sandwich panels by the differential quadrature method

Sandwich structures are widely used in many engineering fields. It is possible but not easy for an engineering theory to recover all stresses accurately. In this paper, a modeling strategy is proposed to simplify the formulation. A classical sandwich panel is firstly divided into three parts, equations of the top and bottom face sheets are used as the boundary conditions of the two-dimensional core and then only the core needs to be analyzed by the differential quadrature method (DQM). In this way, both displacement and stress can be accurately obtained. Detailed formulations are worked out. Three boundary conditions and three types of loading, including the concentrated load regarded as a challenging problem for point discrete methods such as the DQM, are considered to investigate the effect of boundary conditions and loading on the distributions of displacement and stress. For verification, results are compared with theoretical solutions or/and numerical data. Presented data may be a reference for other investigators to develop more accurate engineering beam theory or new numerical method.     

Complete elastic–plastic stress asymptotic solutions near general plane V-notch tips

An efficient method is developed to determine the multiple term eigen-solutions of the elastic–plastic stress fields at the plane V-notch tip in power-law hardening materials. By introducing the asymptotic expansions of stress and displacement fields around the V-notch tip into the fundamental equations of elastic–plastic theory, the governing ordinary differential equations (ODEs) with the stress and displacement eigen-functions are established. Then the interpolating matrix method is employed to solve the resulting nonlinear and linear ODEs. Consequently, the first four and even more terms of the stress exponents and the associated eigen-solutions are obtained. The present method has the advantages of greater versatility and high accuracy, which is capable of dealing with the V-notches with arbitrary opening angle under plane strain and plane stress. In the present analysis, both the elastic and the plastic deformations are considered, thus the complete elastic and plastic stress asymptotic solutions are evaluated. Numerical examples are shown to demonstrate the accuracy and effectiveness of the present method.     

On certain properties of a nonlinear eigenvalue problem for Hamiltonian systems of ordinary differential equations

Properties of the eigenvalues are examined in a nonlinear self-adjoint eigenvalue problem for linear Hamiltonian systems of ordinary differential equations. In particular, it is proved that, under certain assumptions, every eigenvalue is isolated and there exists an eigenvalue with any prescribed index.     

Symmetry and singularity properties of second-order ordinary differential equations of Lie's type III

We extend the work of Abraham-Shrauner [B. Abraham-Shrauner, Hidden symmetries and linearization of the modified Painlevé-Ince equation, J. Math. Phys. 34 (1993) 4809-4816] on the linearization of the modified Painlevé-Ince equation to a wider class of nonlinear second-order ordinary differential equations invariant under the symmetries of time translation and self-similarity. In the process we demonstrate a remarkable connection with the parameters obtained in the singularity analysis of this class of equations.