弹塑性摩擦接触问题形状设计灵敏度分析

提出了一种计算二维有限变形弹塑性摩擦接触问题形状设计灵敏度的算法. 采用主动集策略和mortar方法处理接触边线上的约束条件. 在mortar接触边线的切线和法线方向上采用相同的名义罚函数,提出基于名义罚函数的移动摩擦锥算法来正则化接触约束条件,发展了一种新的二维多体有限变形摩擦接触算法. 在此基础上, 通过将离散形式的摩擦接触问题控制方程对形状设计变量微分,得到了该路径相关问题的直接微分法解析设计灵敏度计算格式, 其节点位移灵敏度方程在每个增量步不用迭代、直接求解. 与国际上现有的二维多体有限变形摩擦接触问题的解析设计灵敏度算法相比,本算法不需分解为法向和切向推导,表达式较简洁,便于编程实现. 数值算例验证了算法的精度和有效性.      

移动简谐荷载作用下桥梁响应的高效计算

在计算移动荷载过桥问题中广泛使用的Newmark方法必须在每一时间步内限制荷载的大小和作用位置都不能改变。精细积分法虽然允许荷载的大小在每一时间步长内发生变化,但是仍假定其作用位置是不变的,未能采取措施以描述荷载沿着桥面的连续移动性。本文提出三种精细积分格式,在每一时间步内不但允许移动荷载的大小按简谐规律连续变化,而且模拟了简谐荷载在空间域的连续移动。通过与Newmark方法和简单问题的解析解进行数值比较,表明用本文提出的方法可以用较粗的结构单元和较大的时间步长而获得很高的计算精度。在精度相同的前提下,计算效率比Newmark方法可提高1~2个数量级。     

双重点最小二乘配点无网格法

配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。     

投影器的标定方法

随着近景摄影测量的发展,投影器的使用越来越频繁和广泛。为了更好地使用和充分地利用投影器,需要求解投影器的内方位元素,即对投影器进行标定。提出了一种方便有效的投影器标定的方法。只需要一架投影器、一部数码相机、一个平面格网和一台计算机,设备的连接和使用相对方便灵活。基于近景摄影测量的理论,结合二维直接线性变换(2D-DLT)和共线方程,运用图像处理的方法,利用光束法平差对投影器进行标定。实验及其结果证明了该标定方法的正确性和有效性。     

分段滤波用于消除脉搏血氧检测中多种运动干扰的研究

脉搏血氧饱和度检测仪可以无损、实时地检测动脉血液中的含氧百分比,从而快速对人体的呼吸状况以及心肺功能作出判断,但其检测精度极易受到测试部位运动的影响.人体测试部位的运动大致可以分为突发性运动干扰以及周期性的运动干扰.突发性干扰表现在局部脉搏波的波形突变上,而周期性的干扰体现在多个局部脉搏波的波形的周期性变化上.文章有针对性地提出了一种将微分阈值分段与迭带中值滤波相结合的方法来剔除突发扰动的方法,并且在上诉方法上再加上一级多点移动平滑以消除周期性运动干扰.实验证明该方法在消除多种运动伪差上有着较好的效果,且该方法计算简单,易于实现,有着一定的实际应用价值.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

Modified characteristic finite difference fractional step method for moving boundary value problem of percolation coupled system

For the coupled system with moving boundary values of multilayer dynamicsof fluids in porous media,a characteristic finite difference fractional step scheme appli-cable to the parallel arithmetic is put forward.Some techniques,such as the change ofregions,the calculus of variations,the piecewise threefold quadratic interpolation,themultiplicative commutation rule of difference operators,the decomposition of high orderdifference operators,and the prior estimates,are adopted.The optimal order estimatesin the l2norm are derived to determine the error in the approximate solution.This nu-merical method has been successfully used to simulate the flow of migration-accumulationof the multilayer percolation coupled system.Some numerical results are well illustratedin this paper.     

Development of a multi‐scale ocean model by using particle Laplacian method for anisotropic mass transfer

The direct injection of CO₂ into the deep ocean is one of the feasible ways for the mitigation of the global warming, although there is a concern about its environmental impact near the injection point. To minimize its biological impact, it is necessary to make CO₂ disperse as quickly as possible, and it is said that injection with a pipe towed by a moving ship is effective for this purpose. Because the injection ship moves over a spatial scale of O(10²km), a mesoscale model is necessary to analyse the dispersion of CO₂. At the same time, since it is important to investigate high CO₂ concentration near the injection point, a small‐scale model is also required. Therefore, in this study, a numerical model was developed to analyse CO₂ dispersion in the deep ocean by using a fixed mesoscale and a moving small‐scale grid systems, the latter of which is nested and moves in the former along the trajectory of the moving ship. To overcome the artificial diffusion of mass concentration at the interface of the two different grid systems and to keep its spatial accuracy almost the same as that in the small‐scale, a particle Laplacian method was adopted and newly modified for anisotropic diffusion in the ocean. Copyright © 2010 John Wiley & Sons, Ltd.     

A study of gridless method with dynamic clouds of points for solving unsteady CFD problems in aerodynamics

When solving unsteady computational fluid dynamics problems in aerodynamics with a gridless method, a cloud of points is usually required to be regenerated due to its accommodation to moving boundaries. In order to handle this problem conveniently, a fast dynamic cloud method based on Delaunay graph mapping strategy is proposed in this paper. A dynamic cloud method makes use of algebraic mapping principles and therefore points can be accurately redistributed in the flow field without any iteration. In this way, the structure of the gridless clouds is not necessarily changed so that the clouds regeneration can be avoided successfully. The spatial derivatives of the mathematical modeling of the flow are directly determined by using weighted least‐squares method in each cloud of points, and then numerical fluxes can be obtained. A dual time‐stepping method is further implemented to advance the two‐dimensional Euler equations in arbitrary Lagarangian–Eulerian formulation in time. Finally, unsteady transonic flows over two different oscillating airfoils are simulated with the above method and results obtained are in good agreement with the experimental data. Copyright © 2009 John Wiley & Sons, Ltd.