微分求积法的特性及其改进

微分求积法已在科学和工程计算中得到了广泛应用。然而,有关时域微分求积法的数值稳定性、计算精度即阶数等基本特性,仍缺乏系统性的分析结论。依据微分求积法的基本原理,推导证明了微分求积法的权系数矩阵满足V-变换这一重要特性;利用微分求积法和隐式Runge-Kutta法的等值性,证明了时域微分求积法是A-稳定、s级s阶的数值方法。在此基础上,为进一步提高传统微分求积法的计算精度,利用待定系数法和Padé逼近,推导出了一类新的s级2s阶的微分求积法。数值计算对比结果验证了所提出的新微分求积法比传统的微分求积法具有更高的计算精度。     

残余力向量法在结构损伤识别中的应用研究进展

工程结构的损伤识别技术对于把握结构工作状态及评估结构的安全性与正常使用性能具有重要的意义。近年来基于残余力向量法的损伤识别技术受到了关注并取得了一定的研究成果。文章从基于残余力向量法的损伤识别技术、残余力向量法和灵敏度分析方法相结合、残余力向量法的改进、残余力向量法和人工神经网络技术的结合、残余力向量法和智能算法的融合等5个方面综述了目前国内外基于残余力向量法进行结构损伤识别研究的成果。并根据残余力向量法应用上存在的问题展望了应用残余力向量法进行结构损伤识别时在如何减小误差;如何克服测试信息不完备的影响;如何进行实际工程损伤识别的研究以及残余力向量法的改进以及残余力向量法和智能算法结合等方面的发展趋势。     

声子晶体中SH波缺陷模的相干叠加理论

利用波的相干叠加原理推导出一维掺杂声子晶体中SH波缺陷模的透射率公式和频率公式,即建立了缺陷模的相干叠加法。将相干叠加法与转移矩阵法和共振理论进行了比较研究,结果表明缺陷模的相干叠加法具备转移矩阵法和共振理论各自的优点,又克服了转移矩阵法和共振理论各自的不足。相干叠加法是研究一维掺杂声子晶体中SH波缺陷模的一种更有效的方法。     

试谈微机汉字库的建库方法及其基本程序

本文介绍了各种微机汉字库的造字建库方法和基本程序,所给出的6个IBMBASICA造字源程序,都是行之有效的。本文介绍的造字建库方法包括点阵汉字的点阵法、单模法和双模法、通用向量汉字的点阵法、转换法、单模法和双模法以及AutoCAD汉字的格栅法和移植法。最后,本文提供了绘图汉字库的开发过程。     

热结构稳态响应的耦合灵敏度分析方法

研究结构稳态热变形和热应力的灵敏度分析方法，给出了直接法和伴随法两种算法。考虑了温度场的耦合作用，在直接法中需要计算温度场对设计变量的导数，在伴随法中需要计算热载荷对温度场的导数。对尺寸和形状两类设计变量的灵敏度分析算例，验证了本文方法的精度。伴随法在应用程序中的实现，为大型结构优化提供了高效率的灵敏度计算方法。     

疲劳短裂纹扩展速率的电测技术

木文研究了间接直流电位法(IDC法)和交流电位法(AC法)用于监测疲劳短裂纹扩展速率的有关技术问题。文中提出了一种用电磁场边界元法对电位法进行标定的理论标定方法。研究了有关的试验标定方法。并对IDC法和AC法的适用范围、灵敏度和稳定性等做了详细的比较分析。     

竖向载荷作用下框架结构力学模型的简化与计算

针对结构力学与工程结构设计中力学模型和载荷作用方式的简化,分别采用一次性加载法、模拟施工三加载法和分层法,研究了竖向载荷作用下框架结构力学模型的简化与计算,分析了三种方法的计算假定和计算结果,给出了合理适用范围。一次性加载法更适用于抗连续性倒塌分析和建筑结构改造设计计算,模拟施工三加载法更适用于新建工程的设计计算,分层法作为一种简化的手算方法可作为理解结构力学原理的基本方法。     

竖向载荷作用下框架结构力学模型的简化与计算

针对结构力学与工程结构设计中力学模型和载荷作用方式的简化,分别采用一次性加载法、模拟施工三加载法和分层法,研究了竖向载荷作用下框架结构力学模型的简化与计算,分析了三种方法的计算假定和计算结果,给出了合理适用范围。一次性加载法更适用于抗连续性倒塌分析和建筑结构改造设计计算,模拟施工三加载法更适用于新建工程的设计计算,分层法作为一种简化的手算方法可作为理解结构力学原理的基本方法。     

一种考虑旋转位移分量的数字图像相关法及其精度分析

提出了一种新的考虑旋转位移分量的数字图像相关法[CD2]旋转相关匹配法,其特点是先对图像子区进行预旋转,再进行相关匹配、相关系数、整像素和亚像素位移的计算。在物体产生旋转位移分量时,该方法可有效提高变形前后图像子区间的相关系数和位移计算精度。通过斜拉伸、剪切和悬臂梁弯曲实验,对旋转相关匹配法和直接相关匹配法进行了对比研究,并用灰度梯度法计算旋转后新图像子区和变形后图像子区间的亚像素位移。结果表明,在旋转位移分量小于10°的情况下,相比于直接相关匹配法,旋转相关匹配法可以得到更高的相关系数;在用灰度梯度法计算亚像素位移条件下,位移精度与试件位移不产生旋转时的直接相关匹配相当。     

流固耦合不可压物质点法及其在晃动问题中的应用

作为一种混合拉格朗日欧拉法,物质点法在流固耦合问题中具有重要的应用前景。对于自由液面的流动问题,基于物质点法框架已建立了弱可压物质点法和完全不可压物质点法,但在处理流固耦合问题时遇到了困难。弱可压物质点法由于采用可压缩状态方程,导致求解时间步长过小,压力振荡严重,产生了非物理的飞溅现象;完全不可压物质点法基于投影算法和不可压条件,消除了弱可压物质点法的压力振荡,提高了时间步长,但难以处理移动边界问题。基于变分形式的投影算法提出了一种新型流固耦合不可压物质点法,得到了体积加权的压力泊松方程PPE(Pressure Poisson Equation),解决了完全不可压物质点法无法处理不规则边界和移动边界的问题。采用流固耦合不可压物质点法研究了运动刚体容器中的液体晃动问题,并与已有实验和数值结果进行对比,验证了算法的正确性和精度。     

On unsteady unidirectional flows of a second grade fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows produced by the sudden application of a constant pressure gradient or by the impulsive motion of one or two boundaries. Exact analytical solutions for these flows are obtained and the results are compared with those of a Newtonian fluid. It is found that the stress at the initial time on the stationary boundary for flows generated by the impulsive motion of a boundary is infinite for a Newtonian fluid and is finite for a second grade fluid. Furthermore, it is shown that initially the stress on the stationary boundary, for flows started from rest by sudden application of a constant pressure gradient is zero for a Newtonian fluid and is not zero for a fluid of second grade. The required time to attain the asymptotic value of a second grade fluid is longer than that for a Newtonian fluid. It should be mentioned that the expressions for the flow properties, such as velocity, obtained by the Laplace transform method are exactly the same as the ones obtained for the Couette and Poiseuille flows and those which are constructed by the Fourier method. The solution of the governing equation for flows such as the flow over a plane wall and the Couette flow is in a series form which is slowly convergent for small values of time. To overcome the difficulty in the calculation of the value of the velocity for small values of time, a practical method is given. The other property of unsteady flows of a second grade fluid is that the no-slip boundary condition is sufficient for unsteady flows, but it is not sufficient for steady flows so that an additional condition is needed. In order to discuss the properties of unsteady unidirectional flows of a second grade fluid, some illustrative examples are given.     

The singular behaviour in the vicinity of intersection between the body and free surface

The singular behaviour in the vicinity of intersection between the body and free surface is presented.It is shown that in the linear regime the singularity of velocity potential for transient problem is in d~2|nd.The singular behaviour for harmonic problem is the same as the result for the transient problem.In particular,the singularity for the harmonic problem with infinite frequency is in d~2 lnd for velocity potential(d is the distance between field point and intersection).     

扰动本构方程及粘弹流体拉伸流动不稳定性

提出“准衰退记忆”新概念,发展了非牛顿流体扰动本构理论,并研究了粘弹流体拉伸流动的不稳定性规律     

Orthogonal arrays: An introduction and their application in optimizing underrelaxation factors in ‘simple’-based algorithm

The primary aim of this paper is to demonstrate how the ‘design-of-experiments’ techniques which are successful in physical experiments could also be adapted to a numerical simulation code. As an example this technique is applied to a general finite difference code used for predicting three-dimensional turbulent recirculating flows. Here the equations for velocities and continuity are solved using the algorithm called SIMPLE, which stands for semi-implicit method for pressure-linked equations. Physical modelling of turbulence is taken care of by means of kinetic energy and turbulence dissipation rate equations. The objective is to optimize the underrelaxation factors of primary and secondary flow variables so that the number of iterations required for convergence is minimum. This is done by the orthogonal array technique (a particular type of design-of-experiment technique). The geometry considered for this purpose is that of a simple gas turbine can combustor and the study is restricted to the isothermal non-reacting condition. Tests are carried out on three different grid configurations. In each case the underrelaxation factor for velocities contributed most to speed up the rate of convergence. Also, for each grid configuration the underrelaxation factor settings for minimum iterations for convergence was found to be same. Hence it is proposed that when doing grid independence tests for any similar flow situation, all the underrelaxation factors could be optimized on coarse grids.     

无网格局部强弱法求解不规则域问题

无网格局部彼得洛夫-伽辽金(meshless local Petrov-Galerkin,MLPG)法是一种具有代表性的无网格方法,在计算力学领域得到广泛应用.然而,这种方法在边界上需执行积分运算,通常很难处理不规则求解域问题.为了克服MLPG法的这种局限性,提出了无网格局部强弱(meshless local strong-weak,MLSW)法.MLSW法采用MLPG法离散内部求解域,采用无网格介点(meshless intervention-point,MIP)法施加自然边界条件,并采用配点法施加本质边界条件,避免执行边界积分运算,可适用于求解各类复杂的不规则域问题.从理论上讲,这种结合式方法,既保持了MLPG法稳定而精确计算的优势,同时兼备配点型方法在处理复杂结构问题时简洁而灵活的优势,实现了弱式法和强式法的优势互补.此外,MLSW法采用移动最小二乘核(moving least squares core,MLSc)近似法来构造形函数,是对传统移动最小二乘(moving least squares,MLS)近似法的一种改进.MLSc使用核基函数代替通常的基函数,有利于数值求解的精确性和稳定性,而且其导数近似计算变得更为简单.数值算例结果初步表明:这种新方法实施简单,求解稳定、精确,表现出适合工程运用的潜力.     

跨音涡轮叶栅底压的无粘模型

本文提出一个无粘的底部压力模型。证明了在给定反压的情况下,超音出口的涡轮叶栅的叶片尾缘的底压以及叶型损失可由该模型确定。与无粘的流场计算相结合,可用该模型计算二维叶栅的底压。对模型的参数研究表明:具有正曲率的喉道下游叶背型面将有利于减少损失;而出口为音速时流动则最为不利。     

Generalized Airy Stress Functions

The classical Airy stress function in planar elastostatics cannot, in general, be a smooth function for multiply connected domains. Moreover, if a non-null body force field is active the classical Airy representation for the stress is not complete. Here, a generalized Airy representation for the stress is presented which preserves smoothness and which is complete. The generalized form identifies the explicit additional pieces that are needed for completeness in multiply connected domains and when a body force field is present.     

An integral representation of the surface-impedance tensor for incompressible elastic materials

A fundamental result in anisotropic elasticity and surface-wave theory is the integral representation for the surface-impedance tensor first derived by Barnett and Lothe in 1973. However, this representation is only valid for compressible materials but not valid for incompressible materials. In this paper the corresponding integral representation for the surface-impedance tensor valid for incompressible materials is derived and is used to establish the uniqueness of surface-wave speed and to obtain an expression for the tensor Green's function for the infinite space. Mathematics subject classifications (2000) 74B05, 74B15, 74B20, 74J15     

Anisotropic damage–plasticity model for concrete

A material model for concrete is proposed here within the framework of a thermodynamically consistent elasto-plasticity–damage theory. Two anisotropic damage tensors and two damage criteria are adopted to describe the distinctive degradation of the mechanical properties of concrete under tensile and compressive loadings. The total stress tensor is decomposed into tensile and compressive components in order to accommodate the need for the above mentioned damage tensors. The plasticity yield criterion presented in this work accounts for the spectral decomposition of the stress tensor and allows multiple hardening rules to be used. This plastic yield criterion is used simultaneously with the damage criteria to simulate the physical behavior of concrete. Non-associative flow rule for the plastic strains is used to account for the dilatancy of concrete as a frictional material. The thermodynamic Helmholtz free energy concept is used to consistently derive dissipation potentials for damage and plasticity and to allow evolution laws for different hardening parameters. The evolution of the two damage tensors is accounted for through the use of fracture-energy-based continuum damage mechanics. An expression is derived for the damage–elasto-plastic tangent operator. The theoretical framework of the model is described here while the implementation of this model will be discussed in a subsequent paper.     

Green's function for SH-waves in a cylindrically monoclinic material

Green's function for SH-waves in a cylindrically monoclinic material is considered for impulsive and time-harmonic sources. Closed form expressions for the Green's function are derived for a few limited values of anisotropic parameters. A very interesting time development of the wave front shape is illustrated and the wave front singularity is discussed for the transient SH-wave. Contours of the displacement amplitude for the time-harmonic wave are also shown.