大型水下浮体强度分析计算方法研究

采用内插算法开发了流固耦合网格转换软件HUSTM2M并提出了大型水下浮体在极限工况下强度分析的计算方法.建立大型水下浮体的流体计算网格以及结构计算网格,利用流体计算软件FLUENT进行流场分析,用流固耦合网格转换软件HUSTM2M将流体网格上的压力值转换到结构网格上,采用软件ABAQUS对结构物进行强度有限元分析.最后采用该计算方法分析由武钢研制的大型水下浮体在南海海域极限工况下的强度.算例表明该计算方法较为实用,具有潜在的工程应用价值.     

结构拓扑优化中变量连接算法研究与软件实现Research and software implementation of variable connection algorithm for structural topology optimization

结构拓扑优化的变量连接,是通过对设计变量之间添加约束关系,从而得到特定的拓扑优化构型,使得优化结果能够满足工程上的特殊要求和工艺制造技术的限制。针对拓扑优化中的几类过滤形式及灵敏度分析,给出了考虑变量连接的计算公式;基于自主研发的SiPESC软件集成化平台,在SiPESC ．TOPO拓扑优化模块上进行二次开发,构建了拓扑优化的变量连接算法框架,其核心思想是基于面向对象设计方法和软件设计模式,实现算法与数据分离。详细阐述了变量连接的作用方式,以及软件框架通用接口设计方案,并通过数值算例验证了其在静力问题、动力问题和热传导问题上的可行性。     

结构拓扑优化中变量连接算法研究与软件实现

结构拓扑优化的变量连接,是通过对设计变量之间添加约束关系,从而得到特定的拓扑优化构型,使得优化结果能够满足工程上的特殊要求和工艺制造技术的限制。针对拓扑优化中的几类过滤形式及灵敏度分析,给出了考虑变量连接的计算公式;基于自主研发的SiPESC软件集成化平台,在SiPESC.TOPO拓扑优化模块上进行二次开发,构建了拓扑优化的变量连接算法框架,其核心思想是基于面向对象设计方法和软件设计模式,实现算法与数据分离。详细阐述了变量连接的作用方式,以及软件框架通用接口设计方案,并通过数值算例验证了其在静力问题、动力问题和热传导问题上的可行性。     

基于SiPESC平台的弹塑性分析的软件框架

基于开放式结构有限元分析软件系统SiPESC.FEMS,针对弹塑性非线性有限元分析理论与流程,采用工厂模式和构造器模式等经典软件设计模式,提出了一种基于面向对象设计方法的通用小变形率无关弹塑性分析软件框架。该框架利用算法与数据模型分离的思想,通过抽象出数据模型类、算法类和算法参数类,支持各功能模块的动态替换和组装,具备开放性和可扩展性。同时,基于SiPESC平台工程数据库的数据管理能力,实现了弹塑性分析中产生的大量历史数据以统一的接口与管理方式进行存取,为大规模问题的求解提供了高效、稳定的支撑。通过多类型实体单元、多类本构模型和多种规模的计算实例,验证了该计算框架在弹塑性有限元程序开发中的通用性、有效性和灵活性。     

一类面向高阶精度自适应流动计算的流场插值方法

网格自适应技术和高阶精度数值方法是提升计算流体力学复杂问题适应能力的有效技术途径. 将这两项技术结合需要解决一系列技术难题, 其中之一是高阶精度流场插值. 针对高阶精度自适应流动计算, 提出一类高精度流场插值方法, 实现将前一迭代步网格中流场数值解插值到当前迭代步网格中, 以延续前一迭代步中的计算状态. 为实现流场插值过程中物理量守恒, 该方法先计算新旧网格的重叠区域, 然后将物理量从重叠区域的旧网格中转移到新网格中. 为满足高阶精度要求, 先采用k-exact最小二乘方法对旧网格上的数值解进行重构, 获得描述物理量分布的高阶多项式, 随后采用高阶精度高斯数值积分实现物理量精确地转移到新网格单元上. 最后, 通过一个具有精确解的数值算例和一个高阶精度自适应流动计算算例验证了本文算法的有效性. 第一个算例结果表明当网格规模固定不变时, 插值精度阶数越高, 插值误差越小; 第二个算例显示本文方法可以有效缩短高精度自适应流动计算的迭代收敛时间.      

基于优化模型的结构声学耦合界面载荷传递算法

针对结构声学耦合系统的界面载荷传递问题,提出了一种基于约束最优化模型的局部参数插值算法。耦合界面处结构和流体单元表面坐标通过各自的形状函数进行插值,把界面载荷在互不匹配的网格节点间的传递问题转化为一个点到用自然坐标表示的有限边界曲面的最小距离问题,以便利用成熟稳定的优化算法对其进行高效求解。与已有方法相比,该算法在耦合界面单元为曲面的情况下仍能保持较高的计算精度。本文给出的数值算例验证了本算法的有效性和可靠性。     

薄膜结构流固耦合的CFD数值模拟研究

基于弱耦合分区求解策略,在CompaqVisualFortran6．5环境下搭建了薄膜结构三维流固耦合效应的CFD数值模拟平台。程序采用模块化编程思想,主要包含几何建模、流体分析、结构分析和数据交换四个模块。其中几何建模模块采用自行编制的膜结构找形分析程序,流体分析模块采用经过二次开发的计算流体力学软件FLUENT6．0,结构分析模块采用自行编制的膜结构动力分析程序MDLFX;在数据交换模块中,编制了基于薄板样条法的插值计算程序,以实现流固交界面上不同区域网格间的数据传递问题,编制了基于代数法和迭代法的动网格变形程序,以实现流固耦合运算中的动网格更新。基于该软件平台,对单向柔性屋盖和鞍形膜结构屋盖进行了流固耦合数值模拟,验证了方法的有效性。     

格子玻尔兹曼法多块网格的交界结构优化研究

基于格子玻尔兹曼方法LBM(Lattice Boltzmann Method)对多块网格方法(Multi-Block)的粗细网格交界结构进行了研究,提出了一种新的优化处理方案。解决了原有网格交界结构存在的三个问题,即两套插值运算造成的程序结构复杂的问题,存储前几个时间步的节点流场数据以备插值运算造成内存浪费的问题和基于时间插值结果进行空间插值计算造成插值误差积累的问题。用一次多点二维空间插值的方式,将原方法的空间和时间双插值,简并成一次空间插值。通过对经典的非定常的圆柱绕流算例和定常的标准顶盖方腔驱动流算例的仿真模拟,验证了交界面处质量、动量及应力的连续性以及网格交界面数据过渡的流畅度,最终验证了改进方法的正确性。数值模拟结果表明,改进后多块算法可实现局部网格细化,进一步推动LBM方法在实际工程问题中的应用。     

近场水下爆炸瞬态强非线性流固耦合无网格数值模拟研究

近场水下爆炸涉及多相流体的掺杂耦合以及结构的大变形、损伤和断裂等瞬态强非线性现象, 传统的网格算法在模拟近场水下爆炸时面临结构网格畸变、多相界面捕捉精度不足等难题, 鉴于此, 本文建立了完全无网格的近场水下爆炸冲击波和气泡全物理过程瞬态强非线性流固耦合动力学模型. 流体采用基于黎曼求解器的光滑粒子流体动力学(SPH)方法求解, 结构采用重构核粒子法(RKPM)求解, 并基于法向通量边界条件实现流固耦合. 为提高SPH对流场间断的求解精度, 引入黎曼问题思想并结合MUSCL重构算法, 为解决流场粒子体积变化剧烈导致的精度下降问题, 应用了自适应粒子分割与合并方法. 为模拟水下爆炸对结构造成的损伤断裂, 基于退化实体几何表述, 采用Lemaitre损伤算法, 建立了RKPM壳结构断裂损伤模型. 依据所建立的SPH-RKPM流固耦合模型, 对近场水下爆炸冲击波传播、气泡脉动与射流以及结构毁伤进行了模拟, 将得到的冲击波载荷、气泡演化以及结构响应与实验值和其他数值解对比, 验证了当前建立的SPH-RKPM流固耦合模型的有效性和精度, 并给出了水下爆炸载荷特性及其对结构的流固耦合毁伤机制与规律, 旨在为近场水下爆炸载荷预报提供理论和基础性技术支撑, 为毁伤威力评估和舰船防护结构设计提供参考.      

结合并行策略与距离减缩法的搜索算法及其在CAA中的应用

为提高流场与声场信息传递效率,建立了一种耦合MPI并行策略与改进距离减缩法的搜索算法。在完成点搜索后,采用一种适合于结构和非结构网格的形函数插值算法进行流场插值,实现了流场信息从流场网格到声场网格的快速传递。针对网格点搜索算法效率的验证,选用二维30P30N三段翼为研究对象,在不同的子进程数下进行对比分析。结果表明,相对于传统的距离减缩算法,本文提出的改进算法能有效提高CAA网格点的搜索效率,并且效率随着子进程数的增加而提高。在此基础上,将建立的流场与声场信息传递技术模块应用于CAA方法中,并对二维NACA0012翼型的后缘噪声问题进行计算分析,计算结果反映基于流场与声场信息快速传递算法的CAA方法能有效模拟宽频气动噪声问题。     

基于SiPESC通用试验设计和代理模型算法构架

基于面向服务的开放软件平台SiPESC,针对试验设计和代理模型共性特征,采用面向服务插件编程的设计方法和软件设计模式,研发了试验设计和代理模型通用算法构架。构架的核心思想是将算法、数据模型和任务管理相分离,形成独立的服务,从而实现算法通用性。整个构架中,试验设计和代理模型分别由五个基本服务构成。重点阐述基本插件所代理功能的抽象过程和通用接口的设计思想,给出算法构架的使用步骤及扩展方式。利用该构架试验设计开发了均匀试验设计、正交试验设计、析因试验设计及中心复合试验设计;代理模型实现了响应面模型、径向基函数模型及Kriging模型,并对部分算法进行数值验证。研究工作表明,算法构架适用于通用试验设计和代理模型算法,可方便进行动态扩展,具备良好的开放性和重用性。     

Computational algorithms for tracking dynamic fluid–structure interfaces in embedded boundary methods

A robust, accurate, and computationally efficient interface tracking algorithm is a key component of an embedded computational framework for the solution of fluid–structure interaction problems with complex and deformable geometries. To a large extent, the design of such an algorithm has focused on the case of a closed embedded interface and a Cartesian computational fluid dynamics grid. Here, two robust and efficient interface tracking computational algorithms capable of operating on structured as well as unstructured three‐dimensional computational fluid dynamics grids are presented. The first one is based on a projection approach, whereas the second one is based on a collision approach. The first algorithm is faster. However, it is restricted to closed interfaces and resolved enclosed volumes. The second algorithm is therefore slower. However, it can handle open shell surfaces and underresolved enclosed volumes. Both computational algorithms exploit the bounding box hierarchy technique and its parallel distributed implementation to efficiently store and retrieve the elements of the discretized embedded interface. They are illustrated, and their respective performances are assessed and contrasted, with the solution of three‐dimensional, nonlinear, dynamic fluid–structure interaction problems pertaining to aeroelastic and underwater implosion applications. Copyright © 2012 John Wiley & Sons, Ltd.     

特征值摄动法在利用动响应结构小损伤检测中的应用

目前在使用遗传算法或神经网络方法进行结构动力学损伤检测,需要基于少量的在线测量损伤结构数据和大量的数值仿真数据来实现,其中通过有限元方法来获得仿真数据的巨大计算量是动力学结构损伤检测方法发展中所面临的一个重要问题。本文在建模方面应用近年来提出的调整单元刚度模拟损伤的先进方法,以保证在损伤前后结构自由度数目不变;在此基础上应用特征值摄动法来减少损伤检测中计算量,并通过对复合材料层合板响应信号的小波分析验证了使用一阶矩阵摄动在有效降低计算量的同时,可以获得对损伤检测而言足够准确的响应信号。     

Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems

Embedded boundary methods for CFD (computational fluid dynamics) simplify a number of issues. These range from meshing the fluid domain, to designing and implementing Eulerian‐based algorithms for fluid–structure applications featuring large structural motions and/or deformations. Unfortunately, embedded boundary methods also complicate other issues such as the treatment of the wall boundary conditions in general, and fluid–structure transmission conditions in particular. This paper focuses on this aspect of the problem in the context of compressible flows, the finite volume method for the fluid, and the finite element method for the structure. First, it presents a numerical method for treating simultaneously the fluid pressure and velocity conditions on static and dynamic embedded interfaces. This method is based on the exact solution of local, one‐dimensional, fluid–structure Riemann problems. Next, it describes two consistent and conservative approaches for computing the flow‐induced loads on rigid and flexible embedded structures. The first approach reconstructs the interfaces within the CFD solver. The second one represents them as zero level sets, and works instead with surrogate fluid/structure interfaces. For example, the surrogate interfaces obtained simply by joining contiguous segments of the boundary surfaces of the fluid control volumes that are the closest to the zero level sets are explored in this work. All numerical algorithms presented in this paper are applicable with any embedding CFD mesh, whether it is structured or unstructured. Their performance is illustrated by their application to the solution of three‐dimensional fluid–structure interaction problems associated with the fields of aeronautics and underwater implosion. Copyright © 2011 John Wiley & Sons, Ltd.     

Efficient geometrical parametrisation techniques of interfaces for reduced-order modelling: application to fluid–structure interaction coupling problems

We present some recent advances and improvements in shape parametrisation techniques of interfaces for reduced-order modelling with special attention to fluid–structure interaction problems and the management of structural deformations, namely, to represent them into a low-dimensional space (by control points). This allows to reduce the computational effort, and to significantly simplify the (geometrical) deformation procedure, leading to more efficient and fast reduced-order modelling applications in this kind of problems. We propose an efficient methodology to select the geometrical control points for the radial basis functions based on a modal greedy algorithm to improve the computational efficiency in view of more complex fluid–structure applications in several fields. The examples provided deal with aeronautics and wind engineering.     

A generic framework for time-stepping partial differential equations (PDEs): general linear methods,object-oriented implementation and application to fluid problems

Time-stepping algorithms and their implementations are a critical component within the solution of time-dependent partial differential equations (PDEs). In this article, we present a generic framework – both in terms of algorithms and implementations – that allows an almost seamless switch between various explicit, implicit and implicit–explicit (IMEX) time-stepping methods. We put particular emphasis on how to incorporate time-dependent boundary conditions, an issue that goes beyond classical ODE theory but which plays an important role in the time-stepping of the PDEs arising in computational fluid dynamics. Our algorithm is based upon J.C. Butcher's unifying concept of general linear methods that we have extended to accommodate the family of IMEX schemes that are often used in engineering practice. In the article, we discuss design considerations and present an object-oriented implementation. Finally, we illustrate the use of the framework by applications to a model problem as well as to more complex fluid problems.     

基于开放式结构有限元系统 SiPESC.FEMS平板壳单元研发

基于开放式结构有限元系统SiPESC.FEMS的单元计算模块的设计模式,研发设计一种通用的平板壳单元计算框架。考虑板壳单元的组合关系和程序编制过程中的重用性及灵活性等特点,采用了软件设计中的构造器(Builder)模式实现不同的组合单元。本框架具有很好的通用性和可扩展性,为有限元程序研发提供了一个新的方式;同时,系统能够处理复杂荷载和边界条件,并可灵活实现不同类型单元的组合分析。本文利用此方法构造五种平板壳单元,通过数值算例分析对比讨论其性能,为选取合适的平板壳单元类型进行结构数值分析提供参考。     

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.