Hybrid stress analysis by digitized photoelastic data and numerical methods

A method is presented that determines photoelastic isochromatic values at the nodal points of a grid mesh which in turn is generated by a computer program that accepts digitized input. Values of σ₁ - σ₂ are computed from the digitized fringe orders. The Laplace equation is solved to separate the principal stresses at each nodal point. The method is extended to digitize isoclinics. Subsequently, σ _x - σ _y and τ _xy are calculated to be used as starting values for the solution of the pertaining partial differential equations to enhance convergence. For further accelerating the rate of convergence, superfluous boundary conditions are added from the digitized data; significant improvement is demonstrated. Estimated values of σ _x - σ _y from the digitized data are further used in conjunction with the solution of the Laplace equation to determine the state of stress without solving the boundary value problems. Paper was presented at the 1988 SEM Spring Conference on Experimental Mechanics held in Portland, OR on June 5–10.     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

Numerical simulation of the interaction of gravity waves with elastic ice floes

The self-consistent motion of a fluid and elastically oscillating plates partially covering the fluid is simulated numerically in the linear approximation. The problem is reduced to the simultaneous solution of the Laplace equation for the fluid and the equation of elastic plate oscillations for the ice. The numerical and analytical solutions, the latter obtained from an integral equation containing the Green’s function, are compared. To solve the problem numerically, the boundary element method for the Laplace equation and the finite element method for the equation describing the elastic plate are proposed. The coefficients of transmission and reflection of surface gravity waves from the floating plates are calculated. It is shown that the solution may be quasi-periodic with characteristics determined by the initial values of the wave and ice-floe parameters. The ice floes may exert a filtering effect on the surface wave spectrum, essentially reducing its most reflectable components. Sankt-Peterburg. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 123–131, May–June, 2000.     

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

Boundary linear integral method for compressible potential flows

A boundary linear integral method based on Green function theory has been developed to solve the full potential equation for subsonic and transonic flows. In this integral method, potential values in the flow region are determined by potential values represented by boundary integrals and a volume integral. The boundary potential values are obtained by implementing the boundary integrals along boundary segments where a linear potential relation is assumed. The volume integral is evaluated in a grid generated by finite element discretization. The volume integral is evaluated only outside the body. Therefore there is no extra boundary treatment required for evaluation of the volume integral. The source term is assumed to be constant in an element integral volume. The volume integral needs to be evaluated only once and can be stored in computer memory for further usage.     

A PIECEWISE–LINEARIZED METHOD FOR ORDINARY DIFFERENTIAL EQUATIONS: TWO–POINT BOUNDARY–VALUE PROBLEMS

Piecewise-linearized methods for the solution of two-point boundary value problems in ordinary differential equations are presented. These problems are approximated by piecewise linear ones which have analytical solutions and reduced to finding the slope of the solution at the left boundary so that the boundary conditions at the right end of the interval are satisfied. This results in a rather complex system of non-linear algebraic equations which may be reduced to a single non-linear equation whose unknown is the slope of the solution at the left boundary of the interval and whose solution may be obtained by means of the Newton–Raphson method. This is equivalent to solving the boundary value problem as an initial value one using the piecewise-linearized technique and a shooting method. It is shown that for problems characterized by a linear operator a technique based on the superposition principle and the piecewise-linearized method may be employed. For these problems the accuracy of piecewise-linearized methods is of second order. It is also shown that for linear problems the accuracy of the piecewise-linearized method is superior to that of fourth-order-accurate techniques. For the linear singular perturbation problems considered in this paper the accuracy of global piecewise linearizat ion is higher than that of finite difference and finite element methods. For non-linear problems the accuracy of piecewise-linearized methods is in most cases lower than that of fourth-order methods but comparable with that of second-order techniques owing to the linearization of the non-linear terms.     

A finite element analysis of mach reflection by using the Boussinesq equation

The numerical analysis of ‘Mach reflection’, which is the reflection of an obliquely incident solitary wave by a vertical wall, is presented. For the mathematical model of the analysis, the two-dimensional Boussinesq equation is used. In order to solve the equation in space, the finite element method based on the linear triangular element and the conventional Galerkin method is applied. The combination of explicit and semi-implicit schemes is employed for the time integration. Moreover, one of the treatments for the open boundary condition, in which the analytical solution of the linearized Boussinesq equation in the outside domain is linked to the discrete values of velocity and water elevation in the inside domain, is applied for the modeling of the Mach reflection problem. © 1998 John Wiley & Sons, Ltd.     

Axisymmetric non-Fourier temperature field in a hollow sphere

The non-Fourier axisymmetric (2+1)-dimensional temperature field within a hollow sphere is analytically investigated by the solution of the well-known Cattaneo–Vernotte hyperbolic heat conduction equation. The material is assumed to be homogeneous and isotropic with temperature-independent thermal properties. The method of solution is the standard separation of variables method. General linear time-independent boundary conditions are considered. Ultimately, the presented solution is applied to a (1+1)—as well as a (2+1)—dimensional problem, and their respective non-Fourier thermal behavior is studied. The present solution can be reduced to special cases of interest by choosing appropriate boundary conditions parameters. Dedicated to Prof. Gholamali Atefi, with appreciation and admiration on the occasion of his 65th birthday.     

薄板弯曲问题的虚边界元-最小二乘法

本文提出求解任意形状的薄板弯曲问题的虚边界元-最小二乘法。本法首先利用薄板弯曲平衡方程的格林函数和离开实际边界上分布的未知的横向荷载和法向弯矩函数建立满足实际边界条件的积分方程;然后采用最小二乘法和沿虚边界分段离散化的待定的分布横向荷载和法向弯矩函数得到求上述积分方程离散化数值解的线性代数方程组。导出了一系列的数值积分的公式,并求解了许多例题,数值结果说明本法完全避免了奇异积分及其复杂的处理方法和耗时的运算,而且在边界及其附近区域解的精度比普通边界元(以后简称边界元)法大大地提高了。     

An analytic separation method for photoelasticity

A method for separating principal stresses in photoelasticity is presented. This method is based upon the series solution of Laplace's equation and the determination of the unknown coefficients arising in this series by a least-squares numerical technique. By selecting an adequate number of terms in the series, the representation of the boundary values of the first stress invariant can be established as accurately as the initial photoelastic data. This form of representation of the first stress invariant at interior points in the region is moe accurate than the boundary values employed.     

The meshless analog equation method: I. Solution of elliptic partial differential equations

A new purely meshless method for solving elliptic partial differential equations (PDEs) is presented. The method is based on the principle of the analog equation of Katsikadelis, hence its name meshless analog equation method (MAEM), which converts the original equation into a simple solvable substitute one of the same order under a fictitious source. The fictitious source is represented by multiquadric radial basis functions (MQ-RBFs). The integration of the analog equation yields new RBFs, which are used to approximate the sought solution. Then inserting the approximate solution into the PDE and boundary conditions (BCs) and collocating at the mesh-free nodal points results in a system of linear equations, which permit the evaluation of the expansion coefficients of the RBFs series. The method exhibits key advantages over other RBF collocation methods as it is highly accurate and the coefficient matrix of the resulting linear equations is always invertible. The accuracy is achieved using optimal values for the shape parameters and the centers of the multiquadrics as well as of the integration constants of the analog equation, which are obtained by minimizing the functional that produces the PDE. Without restricting its generality, the method is illustrated by applying it to the general second order 2D and 3D elliptic PDEs. The studied examples demonstrate the efficiency and high accuracy of the developed method.     

Dual boundary element analysis of wave scattering from singularities

The dual boundary element method is used to obtain an efficient solution of the Helmholtz equation in the presence of geometric singularities. In particular, time-harmonic waves in a membrane which contains one or more fixed edge stringers (or cracks) are investigated. The hypersingular integral equation is used in the procedure to ensure a unique solution for the problem with a degenerate boundary. The method yields a solution for the entire membrane as well as the dynamic stress intensity factor. Numerical results are presented for a circular membrane containing a single edge stringer, two edge stringers and an internal stringer. Also, the first three critical wave numbers of the membrane with the homogeneous boundary condition are determined, and the dynamic stress intensity factors are found for problems with the nonhomogeneous boundary condition. Good agreement is found after comparing the results with exact solutions, and with results obtained using DtN-FEM and ABAQUS.     

Analytical solution for the cylindrical bending vibration of piezoelectric composite plates

An analytical solution for the cylindrical bending vibrations of linear piezoelectric laminated plates is obtained by extending the Stroh formalism to the generalized plane strain vibrations of piezoelectric materials. The laminated plate consists of homogeneous elastic or piezoelectric laminae of arbitrary thickness and width. Fourier basis functions for the mechanical displacements and electric potential that identically satisfy the equations of motion and the charge equation of electrostatics are used to solve boundary value problems via the superposition principle. The coefficients in the infinite series solution are determined from the boundary conditions at the edges and continuity conditions at the interfaces between laminae, which are satisfied in the sense of Fourier series. The formulation admits different boundary conditions at the edges of the laminated piezoelectric composite plate. Results for laminated elastic plates with either distributed or segmented piezoelectric actuators are presented for different sets of boundary conditions at the edges.     

The general solution for axial symmetrical bending of nonhomogeneous circular plates resting on an elastic foundation

In this paper,a new method,the exact analytic method,is presented on the basis of stepreduction method.By this method,the general solution for the bending of nonhomogenouscircular plates and circular plates with a circular hole at the center resting,on an elastfcfoundation is obtained under arbitrary axial symmetrical loads and boundary conditions.The uniform convergence of the solution is proved.This general solution can also be applieddirectly to the bending of circular plates without elastic foundation.Finally,it is onlynecessary to solve a set of binary linear algebraic equation.Numerical examples are givenat the end of this paper which indicate satisfactory results of stress resultants anddisplacements can be obtained by the present method.     

Extension of domain‐free discretization method to simulate compressible flows over fixed and moving bodies

This paper is the first endeavour to present the local domain‐free discretization (DFD) method for the solution of compressible Navier–Stokes/Euler equations in conservative form. The discretization strategy of DFD is that for any complex geometry, there is no need to introduce coordinate transformation and the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior dependent points are updated at each time step to impose the wall boundary condition by the approximate form of solution near the boundary. Some points inside the solution domain are constructed for the approximate form of solution, and the flow variables at constructed points are evaluated by the linear interpolation on triangles. The numerical schemes used in DFD are the finite element Galerkin method for spatial discretization and the dual‐time scheme for temporal discretization. Some numerical results of compressible flows over fixed and moving bodies are presented to validate the local DFD method. Copyright © 2006 John Wiley & Sons, Ltd.     

The complex variable boundary element method for solving sandwich plate problems

The objective of this paper is to develop a new complex variable boundary element method for sandwich plates of Reissner's type and Hoff's type. The general solution of Helmhotz equation in complex field is given. Based on the Vekua's complex integral representation of the analytic function, the new boundary integral equations are formulated. The density function in the integral equation is determined directly by boundary element method. Some standard examples are presented, and the results of numerical solutions are accurate everywhere in the plate. The approach presented is only applicable for bounded simply connected regions. The project is supported by the National Science Foundation of China.     

High-order methods for differential equations with large first-derivative terms

A method is developed to solve elliptic singular perturbation problems. Examples are presented in one and two dimensions for both linear and non-linear problems. In particular, examples are presented for fluid flow problems with boundary layers. In the one-dimensional case an approximating equation is developed using just three points. The method first presented is a fourth-order approximation but is extended to become a higher-order method. Results are included for the fourth-, sixth-, eighth- and tenth-order methods. The results are first compared with results found by Segal in an article about elliptic singular perturbation problems. The elliptic singular perturbation problems are compared with a method by Il'in and also with central and backward difference schemes from Segal's article. There was only one case where the results in Segal's paper were as accurate as the results presented in this paper. However, in this case the method used by Segal did not give accurate values for a second problem presented. The results are also compared with results given by Spalding and by Christie. The method of this paper was also tested on the solution of some non-linear diffusion equations with concentration-dependent diffusion coefficients. The results were superior to results presented by Lee and by Schultz. Finally, the method is extended to several two-dimensional problems. The method developed in this paper is accurate, easy to use and can be generalized to other problems.     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

Transient 1D transport equation simulated by a mixed Green element formulation

New discrete element equations or coefficients are derived for the transient 1D diffusion–advection or transport equation based on the Green element replication of the differential equation using linear elements. The Green element method (GEM), which solves the singular boundary integral theory (a Fredholm integral equation of the second kind) on a typical element, gives rise to a banded global coefficient matrix which is amenable to efficient matrix solvers. It is herein derived for the transient 1D transport equation with uniform and non-uniform ambient flow conditions and in which first-order decay of the containment is allowed to take place. Because the GEM implements the singular boundary integral theory within each element at a time, the integrations are carried out in exact fashion, thereby making the application of the boundary integral theory more utilitarian. This system of discrete equations, presented herein for the first time, using linear interpolating functions in the spatial dimensions shows promising stable characteristics for advection-dominant transport. Three numerical examples are used to demonstrate the capabilities of the method. The second-order-correct Crank–Nicolson scheme and the modified fully implicit scheme with a difference weighting value of two give superior solutions in all simulated examples. © 1997 John Wiley & Sons, Ltd.