Numerical solution of the incompressible Navier–Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier–Stokes equations is presented. The scheme is constructed on a staggered‐mesh grid system. The convection terms are discretized with a fifth‐order‐accurate upwind compact difference approximation, the viscous terms are discretized with a sixth‐order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth‐order difference approximation on a cell‐centered mesh. Time advancement uses a three‐stage Runge–Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Deteriorated Geometrical Stiffness for Higher Order Finite Elements with Application to Panel Flutter

Higher order elements were first design for linear problems where, in certain situations, they present advantages over the lower order elements. A method to efficiently extend their use to geometrical nonlinear problems as panel flutter and postbuckling behavior is presented. The chaotic and limit-cycle oscillations of an isotropic plate are obtained based on direct integration of the discretized equation of motion. The plate is modeled using the von Karman theory and the geometrical nonlinearities are separated in a nonlinear term of the first kind which manifests especially in the prebuckling and buckling regimes, and a nonlinear term of the second kind which is responsible for the postbuckling behavior. A fifth order, fully compatible element has been used to model thin plates while the inplane loads where introduced through a membrane element. The aerodynamics was modeled using the first order 'piston theory. The method introduces the concept of a deteriorated form of the second geometric matrix which is equivalent to neglecting higher order terms in the strain energy of the plate. This allows for a drastic reduction in the computational effort with no observable loss of accuracy. Well established results in the literature are used to validate the method.     

A NEW METHOD TO DETERMINE THE R-CURVE OF CERAMICS-FRACTURE GRATING METHOD

The etching technique of the single-lined zero-thickness specimen grating is applied tothe surface of the SiC fiber toughening Si_3N_4 ceramic composite specimen.The position of the crackand the crack length during the process of crack extension when the load is applied and gradually in-creased can be determined by recording the output voltage value of the Wheatstone bridge in which theceramic specimen with the fracture grating on is located.The crack-growth-resistance(R-curve) of thismaterialis thus obtained.     

陶瓷阻力曲线测定的新方法——断裂栅法

将零厚度单向试件栅刻蚀工艺应用于晶须增韧氮化硅陶瓷基复合材料表面，通过记录试件所在桥路的输出电压值，可得到随载荷不断增大裂纹逐渐向前扩展过程中每一时刻的裂尖位置和裂纹长度，从而得到该种材料的阻力曲线。     

FORCED OSCILLATIONS OF BOUNDARY VALUE PROBLEMS OF HIGHER ORDER FUNCTIONAL PARTIAL DIFFERENTIAL EQUATIONS

ＦＯＲＣＥＤＯＳＣＩＬＬＡＴＩＯＮＳＯＦＢＯＵＮＤＡＲＹＶＡＬＵＥＰＲＯＢＬＥＭＳＯＦＨＩＧＨＥＲＯＲＤＥＲＦＵＮＣＴＩＯＮＡＬＰＡＲＴＩＡＬＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＳＪｉｎＭｉｎｇＺｈｏｎｇ（靳明忠），ＤｏｎｇＹｉｎｇ（董莹），Ｌｉ...     

复合随机振动分析的扩阶系统方法

提出了随机结构系统反应的子空间次序正交分解的思想．基于这一思想，文中导出了一类用于考虑随机激励的随机结构复合随机振动分析的扩阶系统方法，从而可以应用传统的确定性结构随机振动各种分析方法求解复合随机振动问题．作为特例，本文给出了使用模态分析法求解的过程．将文中算例与随机模拟结果相比较，证实了本文思想与方法的实用性．     

基于能量辨识的非线性静电-结构系统的自由度缩减建模方法

给出了一种基于系统能量函数辨识的静电致动微薄板系统自由度缩减建模方法.从Von Karman应变-位移关系式出发,推导出以广义模态坐标为变量的系统动能、应变能以及电容函数的函数表达式.为了将应变能以及电容函数写成广义模态坐标的多变量多项式形式,利用一系列经静态非线性结构有限元计算的结果,拟合得到多变量多项式的未知系数.由Lagrangian方程获得原系统的自由度缩减模型.利用该模型对器件的静/动态特性进行仿真,其计算费用很低.与有限元结果比较,验证了建模方法的正确性.     

TRANSFER MATRIX METHOD FOR ANALYZING VIBRATION AND DAMPING CHARACTERISTICS OF ROTATIONAL SHELL WITH PASSIVE CONSTRAINED LAYER DAMPING TREATMENT

The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.     

弹性静力问题的无网格弱-强形式结合法

基于改进的移动最小二乘(MLS)二阶导数近似,建立了一种求解弹性静力问题的无网格弱-强形式结合法(MLS-MWS)。该方法采用节点离散求解域,通过MLS构造形函数,将求解域划分为边界域和内部域,并分别使用控制方程的局部弱形式和强形式来建立离散系统方程。对强形式中涉及的近似函数二阶导数计算,提出了一种将其转化为求两次一阶导数的方法,与传统方法相比,该方法计算简单、精度高。MLS-MWS法结合了弱、强形式无网格法的优点,Neumann边界条件容易满足,并且只需在边界区域进行积分。文中应用该方法分析了两个弹性力学平面问题,分析结果表明本文方法具有良好的精度和收敛性。     

On the development of a triple‐preserving Maxwell's equations solver in non‐staggered grids

We present in this paper a finite difference solver for Maxwell's equations in non‐staggered grids. The scheme formulated in time domain theoretically preserves the properties of zero‐divergence, symplecticity, and dispersion relation. The mathematically inherent Hamiltonian can be also retained all the time. Moreover, both spatial and temporal terms are approximated to yield the equal fourth‐order spatial and temporal accuracies. Through the computational exercises, modified equation analysis and Fourier analysis, it can be clearly demonstrated that the proposed triple‐preserving solver is computationally accurate and efficient for use to predict the Maxwell's solutions. Copyright © 2009 John Wiley & Sons, Ltd.