薄板的多变量小波有限元分析

基于区间B样条小波和广义变分原理,提出了多变量小波有限元法,构造了一种新的薄板多变量小波有限单元.由广义变分原理推导结构的多变量有限元列式,区间B样条小波尺度函数作为插值函数构造的多变量小波有限元法中,广义应力和应变被作为独立变量进行插值,避免了传统方法中应力应变求解的微分运算,减小了计算误差.区间B样条小波良好的数值...     

接触碰撞动力学的多变量选取方法

在柔性多体的接触碰撞动力学问题中,多变量方法基于附加约束的接触模型,将柔性体的变形用不同变量来描述：接触局部区域的变形用有限元节点坐标描述,非接触局部区域的变形用模态坐标描述,兼顾了计算精度和效率. 将该方法推广到三维空间碰撞问题,对两柔性杆纵向碰撞过程进行动力学仿真,数值结果与实验结果吻合良好,验证了该方法的有效性. 针对柔性体各自区域的变量如何选取的问题,研究了节点取法、模态阶数以及材料参数对计算结果精度的影响,寻找到合理的多变量选取方法,保证精度的同时使自由度得到最大程度的缩减.     

两类基于局部标架梁单元的闭锁缓解方法

对于大转动、大变形柔性体的刚柔耦合动力学问题,基于李群SE(3)局部标架(local frame formulation, LFF)的建模方法能够规避刚体运动带来的几何非线性问题,离散数值模型中广义质量矩阵与切线刚度矩阵满足刚体变换的不变性,可明显地提高柔性多体系统动力学问题的计算效率. 有限元方法中,闭锁问题是导致单元收敛性能低下的主要原因, 例如梁单元的剪切以及泊松闭锁.多变量变分原理是缓解梁、板/壳单元闭锁的有效手段. 该方法不仅离散位移场,同时离散应力场或应变场, 可提高应力与应变的计算精度. 本文基于上述局部标架,研究几类梁单元的闭锁处理方法, 包括几何精确梁(geometrically exact beam formulation, GEBF)与绝对节点坐标(absolute nodal coordinate formulation, ANCF)梁单元. 其中, 采用Hu-Washizu三场变分原理缓解几何精确梁单元中的剪切闭锁,采用应变分解法缓解基于局部标架的ANCF全参数梁单元中的泊松闭锁. 数值算例表明,局部标架的梁单元在描述高转速或大变形柔性多体系统时,可消除刚体运动带来的几何非线性, 极大地减少系统质量矩阵和刚度矩阵的更新次数.缓解闭锁后的几类局部标架梁单元收敛性均得到了明显提升.      

加筋路径驱动的板壳自适应等几何屈曲分析

加筋薄壁结构常被用于航空航天结构的轻量化设计.随着结构尺寸和几何特征的增加,需要更加精细的网格来满足分析精度的要求.传统的等几何方法采用NURBS张量积形式的拓扑结构,使得在分析过程中难以实现局部细化,而全局细化则会增加不必要的自由度.为了提升加筋板壳结构的数值分析精度和效率,提出一种基于RPHT (rational polynomial splines over hierarchical T-meshes)样条的加筋板壳自适应等几何屈曲分析方法.样条网格可以沿着加筋路径进行自适应的局部细化,有效提升低自由度下加筋板壳结构等几何屈曲分析的精度.首先,蒙皮和筋条分别采用RPHT样条曲面和NURBS样条曲线进行建模,几何建模与数值仿真采用统一的几何语言,实现建模与分析的一体化.其次,采用几何投影算法和样条插值算法实现筋条与蒙皮之间的高效高精度强耦合,并建立基于加筋路径驱动自适应网格细化方法.最后,曲线加筋板和网格加筋壳两个算例验证本方法的高效性和鲁棒性,通过与基于NURBS的等几何分析进行对比,本方法能够明显降低分析模型的总自由度.     

作大范围回转运动柔性梁斜碰撞动力学研究

为正确估计由于碰撞引起的多柔体系统动力学特性的变化,针对作大范围回转运动的柔性梁与一固定斜面发生斜碰撞的情况,在考虑刚柔耦合效应的多柔体系统动力学建模理论的基础上,利用假设模态法建立起重力场作用下的柔性梁一致线性化动力法向碰撞过程中系统的动力行为。基于Ｈｅｒｔｚ接触理论和非线性阻尼项建立法向碰撞接触模型,基于线性切向接触刚度建立柔性梁切向碰撞接触模型,提出的数值算法保证了计算结果的合理性,给出的仿     

非均匀矩形网格的局部网格细化LBM算法研究

传统的格子玻尔兹曼方法 (LBM),特别是基于均匀正方形网格的经典单松弛计算模型(SLBM),其算法鲁棒性和数值稳定性较差,限制了LBM的发展和应用.而网格细化策略可以有效缓解这一窘境,但是传统LBM中网格细化必然会导致计算效率骤降,计算设备要求攀高.为了解决这一问题,文章基于非均匀矩形网格结构,结合插值LBM算法的思路,在保证物面处和流动变化剧烈区域的局部网格细化以及计算精度的前提下,提出了25点拉格朗日插值LBM算法.以经典顶盖驱动方腔内流为算例,开展了包括不同网格分辨率和插值格式的对比分析研究.验证算例既包括了定常流动的数值模拟,也涉及了非定常周期性流动的求解.计算结果表明,相较于其他插值格式,拉格朗日插值格式表现优异;文章局部网格细化工作可以确保物面处及流动变化剧烈区域流动细节的捕捉;数值模拟算法可以为数值仿真提供可信的计算结果;同时大幅降低了总网格数量.因此很大程度上提升了计算效率;数值模拟方法鲁棒性较好,适用于包括定常和非定常流动的数值模拟.     

基于适合分析T样条的高阶数值流形方法

数值流形方法是一种非常灵活的数值计算方法,连续体的有限单元方法和块体系统的非连续变形分析方法只是这一数值方法的特例.数值流形方法中高阶位移函数的构造可通过提高权函数的阶次来实现,这种方法往往需要沿单元边界配置适当的边内节点,这些结点的出现增加了前处理的复杂性,特别是对于大型复杂的空间问题.另一方面,在数值流形方法中可通过缩小单元尺寸(h加密)来提高求解精度.当模拟裂纹扩展时,这种细化策略可用来克服裂纹尖端的奇异性.一个传统的解决方案是细化整个网格,但这会导致计算效率的显著降低.将适合分析的T样条(analysis-suitable T-spline,AST)引入数值流形方法中来建立高阶数值流形方法的分析格式,有效的避免了该问题的出现.AST样条基函数具有线性无关,单位分解,局部加密等许多重要性质,使得其非常适合用于工程设计及分析.在引入AST样条后,可通过改变数学覆盖的构造形式建立不同阶次的数值流形方法分析格式;AST样条自身的局部加密性质也使得数值流形方法中的数学网格局部加密更容易实现.算例结果表明:随着AST样条基函数阶次的提高,数值流形方法的计算结果有了明显的改善;基于AST样条基函数的数值流形方法在保持计算精度的前提下降低了自由度的数量.     

柔性梁刚柔耦合碰撞动力学及变形传播研究

运用柔性多体系统刚柔耦合动力学理论,研究了作大范围回转运动柔性梁的碰撞动力学问题.考虑柔性梁的横向变形,以及横向变形引起的纵向缩短项即非线性耦合变形项.采用基于Hertz接触理论及非线性阻尼理论的非线性弹簧阻尼模型来求解碰撞过程中产生的碰撞力,运用第二类拉格朗日方程建立了系统的刚柔耦合碰撞动力学方程.编制仿真软件进行动力学仿真计算,得到了碰撞力和系统动力学响应,对比分析了不同动力学模型对系统动力学响应的影响.同时研究了碰撞导致的柔性梁横向变形传播的波动特性.     

含面内转动自由度的广义协调曲面矩形扁壳元

扁壳单元中引入结点转角自由度可以在不增加结点的情况下,增加位移场的阶次,提高计算精度,从而显著地提高单元性能。同时在单元中引入泡状位移场,能有效地扩大了单元位移场的解空间,所构造的单元具有计算精度高、对计算网格畸变不敏感的优良特性。本文利用广义协调薄板单元RGC-12的位移函数作为扁壳元的法向位移,利广义协调矩形膜元的位移函数作为扁壳面的切向位移,通过附加面内转动自由度构造了一个具有24个自由度的4结点广义协调曲面矩形扁壳元GRC-S24。在此基础上再增加一个广义泡状位移,又构造了一个具有更高计算精度的曲面矩形扁壳元GRC-S24M。并通过实例分析对这两个单元的收敛性和精度进行了验证。     

柔性多体碰撞问题的多变量方法

针对柔性多体系统碰撞问题的特点, 提出了柔性多体系统碰撞问题的多变量方法. 在未碰撞阶段和分离后的阶段, 多变量方法考虑了物体大范围运动与变形的耦合.在碰撞阶段, 考虑非碰撞区域的变形, 对非碰撞区域用浮动参考系方法求解, 对发生碰撞的局部区域用非线性有限元方法求解. 非碰撞区域采用浮动基位形坐标和柔性体相对浮动基的变形模态坐标描述, 大大减缩了变量的维数,提高了计算效率. 碰撞局部区域采用非线性有限元节点坐标描述, 可以得到碰撞局部区域高精度的应力、应变响应, 而且可以反映碰撞局部区域的大变形、塑性等非线性响应. 该方法既能体现碰撞对系统大范围运动和长期动力学仿真的影响, 又能精确反映出碰撞发生的微小时间段内物体碰撞局部区域的应力、应变响应, 表现出碰撞过程中物体接触区域的演化历程.介绍了两杆正碰撞与铅垂面内双摆撞击自由圆盘的算例, 设计了两柔性杆正碰撞实验,多变量方法的数值仿真结果与实验测量结果很好地吻合.      

基于最佳质量网格的薄板问题的非协调流形方法

对于大部分非协调板单元,使用规则网格能得到很好的效果。但是,当网格不规则时,非协调元的数值特性将变得很差,甚至收敛性得不到保证。为解决网格依赖性问题,许多专家学者提出了改造单元,如拟协调元法和广义协调元法,这些方法能解决收敛性问题,但是数值实践证明没有一种单元能在所有情况下都具有良好的数值特性。考虑到流形方法采用两套完全独立的覆盖系统,可以用规则的数学网格来作为数学覆盖进行插值,取得最佳的插值效果,单元收敛性便能得到保证。再结合适用于流形方法的变分提法,建立起流形方法处理非规则物理边界非协调板单元的一般格式。以ACM 薄板单元为例,与ANSYS、拟协调元法和广义协调元法进行了对比,证明本文方法在处理具有曲线边界的薄板弯曲问题时具有收敛快和精度高等优势。     

基于李群局部标架的多柔体系统动力学建模与计算

多柔体系统动力学主要研究由多个具有运动学约束、存在大范围相对运动的柔性部件构成的动力学系统的建模、计算和控制.多柔体系统不仅具有柔体大变形导致的几何非线性,更具有大范围刚体运动引起的几何非线性,其非线性程度远高于计算结构力学所研究的几何非线性问题.本文基于李群局部标架(local frame of Lie group, LFLG),讨论如何发展一套新的多柔体系统动力学建模和计算方法体系, 具体内容包括:基于局部标架的梁、板壳单元,适用于长时间历程计算的多柔体系统碰撞动力学积分算法,结合区域分解技术的大规模多柔体系统动力学并行求解器, 以及若干验证性算例.上述基于李群局部标架的方法体系可在计算中消除刚体运动带来的几何非线性问题,使柔体系统的广义惯性力、广义弹性力及其雅可比矩阵满足刚体运动的不变性,使多柔体系统动力学与大变形结构力学相互统一,有望推动新一代多柔体系统动力学建模和计算软件的发展.      

Development of quadrilateral spline thin plate elements using the B-net method

The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the Bnet method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian coordinates. In this paper, a thin plate spline element is developed based on the spline element L8 and the refined technique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.     

Bending of a nanoplate with strain-dependent surface stress containing two collinear through cracks

Meccanica - The bending fracture problem of two thickness-through collinear cracks of equal length in a flexible nanoplate with surface stress is analyzed. Using the Kirchhoff thin plate theory...     

Efficient coupling of direct forcing immersed boundary-lattice Boltzmann method and finite element method to simulate fluid-structure interactions

Fluid-structure interactions (FSI) of rigid and flexible bodies are simulated in this article. For the fluid flow, multidirect forcing immersed boundary method (IBM) is adopted to capture the moving boundary, and lattice Boltzmann method (LBM) is used to evolve the flow field. Compared with our previous no-penetration IBM, less iterations are required in this work. In addition, larger velocity in lattice units can be used and the nonphysical force oscillations are suppressed due to the C³ 6-point kernel. Multi-relaxation-time collision operator and local grid refinement are also adopted in LBM to enhance the numerical stability and efficiency. To improve the efficiency of the FSI coupling algorithm, the mesh of the deformable structure can be coarser than the Lagrangian mesh using Newton-Cotes formulas to integrate the traction on the structure surface. A variety of benchmarks, including flow around a circular cylinder with Reynold number ranging from 20 to 200, forced oscillation of a circular cylinder, vortex-induced vibration (VIV) of an elastically mounted circular cylinder in two degrees of freedom, and VIV of an elastic cantilever beam attached to a circular cylinder, are carried out to evaluate the accuracy and stability of the present coupling algorithm. For the benchmark FSI problem considered in this article, a reduction of 54% of the calculation time is achieved using coarser structure mesh. As an application of the FSI coupling algorithm, the mechanism of an elastic beam in the wake of a circular cylinder is discussed.     

Nonlinear dynamics of three-dimensional belt drives using the finite-element method

In this paper, new nonlinear dynamic formulations for belt drives based on the three-dimensional absolute nodal coordinate formulation are developed. Two large deformation three-dimensional finite elements are used to develop two different belt-drive models that have different numbers of degrees of freedom and different modes of deformation. Both three-dimensional finite elements are based on a nonlinear elasticity theory that accounts for geometric nonlinearities due to large deformation and rotations. The first element is a thin-plate element that is based on the Kirchhoff plate assumptions and captures both membrane and bending stiffness effects. The other three-dimensional element used in this investigation is a cable element obtained from a more general three-dimensional beam element by eliminating degrees of freedom which are not significant in some cable and belt applications. Both finite elements used in this investigation allow for systematic inclusion or exclusion of the bending stiffness, thereby enabling systematic examination of the effect of bending on the nonlinear dynamics of belt drives. The finite-element formulations developed in this paper are implemented in a general purpose three-dimensional flexible multibody algorithm that allows for developing more detailed models of mechanical systems that include belt drives subject to general loading conditions, nonlinear algebraic constraints, and arbitrary large displacements. The use of the formulations developed in this investigation is demonstrated using two-roller belt-drive system. The results obtained using the two finite-element formulations are compared and the convergence of the two finite-element solutions is examined.     

The use of dynamic grid adaptation algorithms for the modelling of flow around a circular cylinder in sub‐critical flow regime

In the present study a dynamic grid adaptation (DGA) algorithm is utilized for predicting flow around a circular cylinder in sub‐critical flow regime at a Reynolds number of 1.4×10⁵. The reason for adopting a DGA algorithm is the unsteadiness of the flow field which makes a conventional mesh inefficient. The concept being adopted is to concentrate mesh refinement in regions with high gradients and high turbulent viscosity, while in the region further downstream where the flow is fully developed a coarser mesh will develop and turbulence is modelled with the large eddy simulation (LES) turbulence model. The aim of the study is to present an appropriate variable for mesh refinement, which accomplishes a high rate of mesh refinement in the region with high gradients. The new variable is a product of the local mesh cell size and the rate of strain and includes two additional variables to allow control over the refinement behaviour. The results are compared with experimental data at the corresponding Reynolds number and also with numerical results obtained with conventional mesh. It is demonstrated that DGA algorithms can give results of a very high quality for a mesh that is significantly smaller than for a conventional mesh. Copyright © 2003 John Wiley & Sons, Ltd.     

广义有限差分法求解Kirchhoff和Winkler薄板弯曲问题

本文将广义有限差分法用于数值计算Kirchhoff板和Winkler板的弯曲问题.广义有限差分法是基于最小二乘原理的一种区域型无网格方法. 相比于传统的网格类数值解法,广义有限差分法无需网格生成且无需数值积分.通过数值实验结果表明,广义有限差分法可以有效地求解两类薄板在不同横向荷载作用下的弯曲问题.