基于任意四边形单元的约束阻尼板的模态分析

约束层阻尼技术目前广泛应用于薄壁结构的减振降噪中,关于约束阻尼板的有限元分析基本采用矩形或三角形单元,但用于模拟不规则形状结构时,会带来形状拟合上的困难或精度上的不足。基于离散Kirchhoff理论和Layer-wise层合板理论,利用Hamilton原理推导了约束阻尼板的任意四边形单元,并在此基础上考虑粘弹性材料本构的频率相关性,给出了约束阻尼结构复特征值问题的迭代求解算法。数值算例对不同形状的约束阻尼板结构进行了模态分析,通过与解析解、实验结果及有限元结果的对比,表明了本文单元的有效性和对不规则形状结构的适用性。     

An analytical solution for the stress field and stress intensity factor in an infinite plane containing an elliptical hole with two unequal aligned cracks

The existing analytical solutions are extended to obtain the stress fields and the stress intensity factors (SIFs) of two unequal aligned cracks emanating from an elliptical hole in an infinite isotropic plane. A conformal mapping is proposed and combined with the complex variable method. Due to some difficulties in the calculation of the stress function, the mapping function is approximated and simplified via the applications of the series expansion. To validate the obtained solution, several examples are analyzed with the proposed method, the finite element method, etc. In addition, the effects of the lengths of the cracks and the ratio of the semi-axes of the elliptical hole (a/b) on the SIFs are studied. The results show that the present analytical solution is applicable to the SIFs for small cracks.     

Direct numerical simulation of turbulent flows through concentric annulus with circumferential oscillation of inner wall

Periodic wall oscillations in the spanwise or circumferential direction can greatly reduce the friction drag in turbulent channel and pipe flows. In a concentric annulus, the constant rotation of the inner cylinder can intensify turbulence fluctuations and enhance skin friction due to centrifugal instabilities. In the present study, the effects of the periodic oscillation of the inner wall on turbulent flows through concentric annulus are investigated by the direct numerical simulation (DNS). The radius ratio of the inner to the outer cylinders is 0.1, and the Reynolds number is 2 225 based on the bulk mean velocity U_m and the half annulus gap H. The influence of oscillation period is considered. It is found that for short-period oscillations, the Stokes layer formed by the circumferential wall movement can effectively inhibit the near-wall coherent motions and lead to skin friction reduction, while for long-period oscillations, the centrifugal instability has enough time to develop and generate new vortices, resulting in the enhancement of turbulence intensity and skin friction.     

非结构四边形网格下的一类保对称有限体元格式

针对定常扩散问题,在非结构四边形网格下,通过选取特殊的控制体和有限体元空间,建立两种保对称有限体元格式,在拟一致网格剖分下,当扩散系数光滑时,证明有限体元解函数在L²和H¹范数下均具有饱和误差阶.数值实验验证理论结果的正确性,同时表明新格式对扭曲大变形四边形网格、间断系数问题具有较强的适应性.在正交网格下,第二种格式对流(flux)函数在单元中心点的值还具有超逼近性.     

具有内部特征约束的四边形网格生成方法

提出一种具有内部特征约束的四边形网格自动生成方法,可以处理内部具有孔洞、约束线、约束点以及密度线、密度点等特征约束,满足数值分析中对网格生成的特殊要求。以区域分解法为基础,提出了确定最佳切割线的方法以及在切割线上生成过渡均匀网格节点的方法;将约束线、约束点以及密度线、密度点作为面积为零的孔洞,提出了对特征约束的处理方法;针对特征约束分隔区域的情况,提出了自动确定各子区域边界的方法。基于本文提出的方法已开发出具有特征约束的四边形网格自动生成程序,并成功应用于洪水分析系统和楼面设计分析系统。     

基于面积坐标与B网方法的四边形样条单元

传统等参元方法中, S型等参元完备阶较低,对网格畸变敏感, L型等参元具有高阶完备性但需要使用内部节点. 另外,由于引入等参变换, 采用数值积分可能导致总刚度矩阵出现奇异性.利用三角形面积坐标与B网方法建立了一类平面四边形的样条单元函数,它们的特点是满足协调条件, 克服网格畸变敏感性.其中8节点和12节点单元分别为2次和3次样条函数,对直角坐标分别具有二阶和三阶完备性, 高于相同节点的S型等参元.通过算例测试了这些样条单元, 并与等参元和其它四边形单元比较,数值结果显示了它们的高精度和有效性.      

一个不闭锁和抗畸变的四边形厚板元

构造一个彻底消除剪切闭锁现象并且对网格畸变不敏感的四边形厚薄板通用单元RPAQ。在方法上有三个特点：第一,在厚板挠度和转角的试函数中,采用了合理匹配方案,从而在源头上彻底消除了剪切闭锁现象;第二,采用四边形面积坐标,以代替通常的等参坐标,从而使网格畸变时仍然保持高精度;第三,采用广义协调元做法,使协调条件的采用灵活多样,并保证单元的收敛性。进行了一系列数值例题测试,表明单元RPAQ能自动消除闭锁现象,在由薄板到厚板的不同情况下,在各种网格畸变的情况下,都能体现出良好的精度和数值稳定性。     

基于非结构四边形网格的WENO格式

基于非结构四边形网格发展求解双曲守恒律的三阶加权基本无振荡（WENO）格式.针对任意非结构四边形网格选取重构模板，并给出基于线性多项式的三阶线性重构.但对于一般的非结构四边形网格，会出现非常大的线性权和负权，使得非线性重构的WENO格式对光滑问题也不稳定.本文给出一个处理非常大的线性权的优化重构方法，对优化后得到的负线性权采用分裂方法进行处理.对于非线性权，提出一种考虑局部网格和物理量间断的新光滑度量因子.采用优化重构方法和新的非线性权，当前的三阶WENO格式在质量很差的网格上也具有很好的稳定性.理论的三阶精度在数值精度测试算例中得到验证，同时一范数和无穷范数的误差绝对值不依赖于网格质量；具有强间断的数值结果证明了当前格式的有效性.     

Optimal Hierarchies for Quadrilateral Surfaces

Multiresolution representation of quadrilateral surface approximation (MRQSA) is a useful representation for progressive graphics transmission in networks. Based on two requirements: (1) minimum mean square error and (2) fixed reduction ratio between levels, this paper first transforms the MRQSA problem into the problem of solving a sequence of near-Toeplitz tridiagonal linear systems. Employing the matrix perturbation technique, the MRQSA problem can be solved using about 24mn floating-point operations, i.e. linear time, if we are given a polygonal surface with (2m–1)×(2n–1) points. A numerical stability analysis is also given. To the best of our knowledge, this is the first time that such a linear algebra approach has been used for solving the MRQSA problem. Some experimental results are carried out to demonstrate the applicability of the proposed method.     

Gears,pregears and related domains

We study conformal mappings from the unit disc to one-toothed gear-shaped planar domains from the point of view of the Schwarzian derivative. Gear-shaped (or “gearlike”) domains fit into a more general category of domains we call “pregears” (images of gears under Möbius transformations), which aid in the study of the conformal mappings for gears and which we also describe in detail. Such domains being bounded by arcs of circles, the Schwarzian derivative of the Riemann mapping is known to be a rational function of a specific form. One accessory parameter of these mappings is naturally related to the conformal modulus of the gear (or pregear) and we prove several qualitative results relating it to the principal remaining accessory parameter. The corresponding region of univalence (parameters for which the rational function is the Schwarzian derivative of a conformal mapping) is determined precisely.