基于点球弹簧修匀法的混合网格变形并行算法A parallel algorithm for hybrid mesh deformation based on Vertex-Ball Spring Smoothing

针对二维/三维混合网格,提出基于点球弹簧修匀法的并行网格变形算法。按特定模板将混合网格中的非三角形/四面体单元分解成三角形/四面体单元。针对每个内部节点及其相邻节点建立相应的子弹簧系统,并通过增加Ball-Vertex弹簧避免弹簧系统的塌陷问题。由于点球弹簧法在计算中逐点对网格内部节点进行计算,在计算过程中具有良好的弱耦合性质,因此有利于算法并行化。在并行化时仅需对网格进行虚拟分区操作,不必进行复杂的几何分区操作,同时避免了混合网格不同单元之间的兼容性问题。该方法适用于具有复杂外形的大规模混合网格的变形问题,能够显著提高网格变形的效率,同时具有良好的适应性。     

Investigating cross-wind stability of high-speed trains with large-scale parallel CFD

ABSTRACT

The side-wind loading on a simplified train model at scale 1:25 is investigated by parallel large eddy simulation (LES) with incompressible solvers from the OpenFOAM package and a novel dynamically adaptive, parallel LES-type lattice Boltzmann method (LBM) implemented in our own AMROC framework. It is found that the new LBM code provides more accurate time-averaged force predictions, while compute times are reduced.     

非匹配结点的有限单元等参插值方法及其应用研究

针对两零件的异构网格单元结点在接触界面不能相互匹配导致结点属性不能连续过渡和传递的问题,提出非匹配结点的有限单元等参插值方法,通过构建非匹配结点的形函数和修正原结点的形函数,将结点属性值的影响范围限制在可控的局部区域,从而实现两异构网格结点属性在接触界面的连续过渡和传递。通过两个异构的四边形单元网格的结点属性在接触界面的过渡实例和啮合齿轮的接触分析应用,验证了该方法的正确性和有效性。     

Nodeless variable finite element method for heat transfer analysis by means of flux-based formulation and mesh adaptation

Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method. The solution accuracy is further improved by implementing an adaptive meshing technique to generate finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems are: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to a moving heat source. The English text was polished by Yunming Chen.     

Radial electroelastic nonstationary vibration of a hollow piezoceramic cylinder subject to electric excitation

The paper studies the radial nonstationary vibration of a piezoceramic cylinder polarized throughout the thickness and subjected to a dynamic electric load. A numerical algorithm for solving an initial–boundary-value problem using mesh-based approximations and difference schemes is developed. The dynamic electroelastic state of the cylinder subjected to a constant potential difference applied instantaneously is analyzed Translated from Prikladnaya Mekhanika, Vol. 45, No. 2, pp. 30–35, February 2009.     

两种湍流模型时域颤振计算方法研究

采用时域计算分析方法进行了机翼跨音速颤振特性研究。在结构运动网格的基础上,采用格点格式有限体积方法进行空间离散和双时间全隐式方法进行时间推进求解雷诺平均N-S方程。针对流动粘性分别应用了SST湍流模型和SSG雷诺应力模型,通过对跨音速标模算例AGARD445.6机翼的计算结果与实验值的对比分析,其中应用SST湍流模型得到的颤振速度与实验值最为接近,特别是在跨音速段平均相对误差在3%以内;并且计算结果整体上反映了跨音速颤振"凹坑"物理特性,验证了方法的有效性。     

混凝土三维细观接触面模型数值模拟与CT试验验证

对混凝土三维随机骨料模型进行了改进,在骨料和砂浆之间加入了"面-面接触单元"模拟两者的接触特性,编制了相应的命令程序,称新模型为混凝土接触面模型.采用双折线损伤演化模型进行计算,对数值模拟结果进行了分析研究.从混凝土破坏过程图和荷载-位移曲线图两方面与CT试验结果进行了比较,表明试件破坏时裂纹的萌生、扩展过程与CT试验...     

基于双重自适应的六面体网格再生成方法

提出了一种以栅格法为基本方法,基于几何特征和物理场量双重自适应的六面体网格再生成方法。首先,依据旧网格的表面曲率和几何特征,采用基于栅格法的几何自适应网格再生成方法,生成密度受控的基础网格;然后,将旧网格的物理场量传递到基础网格中;最后,采用有限元误差估计方法对新网格单元的计算误差进行估计,对误差较大的单元进行加密,减...     

基于非局部塑性模型的应变局部化理论分析及数值模拟

通过求解一个第二类Fredholm方程,得到了基于非局部塑性软化模型的应变局部化问题理论解,结果表明,只有在当采用过非局部修正形式的非局部塑性软化模型才能得到应变局部化解,且得到的塑性应变分布和荷载响应依赖于所引入的特征长度及过非局部权参数。通过一维应变局部化有限元数值解,验证了非局部理论的引入能克服计算结果的网格敏感...     

轴对称高速碰撞问题的拉格朗日无结构网格有限体积法的并行计算

给出计算轴对称高速碰撞问题的拉格朗日无结构三角形网格有限体积法的并行格式，并给出以小巨型机ＡｌｉａｎｔＦＸ／４０为目标计算机的算例数值模拟结果和效率分析