基于矩阵误差的流动自适应求解

为了构造一个网格沿流场特征方向进行分布,以保证整个流场的误差尽量均匀化的网格系统,本文对可压缩流动的自适应求解进行了研究.研究内容包括基于误差矩阵的流场特征探测方法和网格自适应方法.网格自适应方法包含网格移动、网格加密、网格稀疏以及这3种方法的混合.基于误差矩阵的流场探测方法能找到流场特征的主方向, 通过网格自适应使得网格按主方向进行分布,可以减少该区域的误差而得到高精度的解.以边为单位能使得该流场的误差尽可能均匀化.本文以捕捉激波为例,对该方法的过程及应用进行了描述,结果验证了该方法的可靠性.     

三角网格模型重建中的误差分析与显示

根据实物模型的坐标测量，数据重建模拟的三角网格曲面表示在实践中有着广泛的应用。在许多对精度有较高要求的应用中，对重建网格模型的误差分析是必不可少的。本文提出并实现了一种对重建三角网格模型进行误差分析的方法，该方法首先对网格模型空间进行划分以提高算法的效率，通过计算测量点与网格模型间的最短距离，得到网格模型的误差分布，最后将计算得到的三角网格模型中每个顶点的误差度量转换为相应的颜色表示，并用彩色云图对分析结果进行直观的显示。     

三维网格动态细化技术在切削仿真中的应用

为解决切削加工有限元仿真模型求解速度与精度难以平衡的问题,提出了一种适用于三维有限元模型的网格动态细化算法。该算法的主要功能包括网格细化区域判断、网格细化以及新旧网格物理场传递。采用Python语言对Abaqus软件进行二次开发,将该算法应用于Ti2AlNb钛合金车削加工仿真之中,并最终通过实验验证了仿真模型的准确性。与采用局部网格细化的仿真模型计算结果对比,采用网格动态细化技术的仿真模型求解切削力误差增加了7.3%,求解最大应力误差增加0.2%,求解速度提升174.2%。实现了在保证仿真计算精度的基础上,有效提高运算速度。     

基于有限元网格变形的飞机外形曲面修改

根据变形后的有限元模型修改原始的飞机外形,从而自动得到变形以后的产品模型.考虑到有限元网格与原始曲面之间的误差及细小特征上的差别,本文采用两种方法来重新生成产品的曲面:(1)通过有限元网格节点的变形直接修改原始曲面,(2)考虑原始曲面与初始有限元模型的误差直接对变形后的有限元网格进行曲面的重构.实例表明重构以后的模型具有良好的外形,能有效地恢复有限元网格中忽略掉的特征,同时曲面之间的拓扑结构能够与原始模型保持一致.     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

圆柱壳板非线性屈曲的分叉解

采用商用有限元分析软件MSC/Nastran,对两直边铰支、其余两边自由的各向同性圆柱壳板中心受横向集中载荷下的非线性屈曲和后屈曲过程进行了分析,简单地采用了不对称网格技术就可获得该基准问题的分叉屈曲解,计算得到分叉屈曲解的极限载荷比基准解极限载荷小10%左右。对采用了多种不对称网格划分得到的模型进行了分析,以研究不同网格划分对最终结果的影响。结果表明,虽然不对称网格技术十分简单,但是也有一定的局限性,文中给出了能得到分叉屈曲解的网格划分方法,还分析了类似的复合材料圆柱壳板问题。研究结果表明,即使采用对称网格划分也能够得到该问题的分叉屈曲解,并阐述了其原因。     

一种网格并行自适应处理方法

给出了一种可压流数值模拟的网格并行自适应处理方法.首先,发展了一种简单的误差估算方法.然后,以初网格上获得的当地误差估算值为权,应用递归谱对剖分方法(Recursive spectral bisection,RSB)划分初网格,以解决负载平衡问题并保证各子域上总体误差近似相等.接着,以误差值为判据对网格进行了并行的自适应处理.最后,在整合的总体网格上应用Galerkin有限体积法分区并行求解Euler方程.为验证方法的可靠性,给出了数值模拟结果.     

壁面冷却的三股混合气流的数值计算

对双层壁扩压器与波瓣喷管组合的三股气流的引射混合进行了流场温度场的数值计算。计算域划分为各子域 ,通过求解拉普拉斯方程 ,分别生成三维贴体曲线坐标网格 ,边界网格加密且波瓣处正交 ,各子域组合成整个网格。为避免因波瓣造成的网格强烈的非正交而引起解的发散 ,本文采用了同位网格 SIMPLEC计算方法和Chen-Kim修正的 k-ε湍流模型。固体区域采用大粘性解决流固耦合 ,还结合了线性欠松弛和拟瞬态欠松弛。结果证明 :双层壁间有外界冷气流被引射进入 ,形成壁面的冷却气流 ,相对单层壁扩压器 ,双层壁扩压器的壁面温度明显降低。对于小扩压角 ( 0～ 5°) ,随着扩压角的增大 ,壁面平均温度降低 ;但当扩压角较大时 ( 5～ 1 0°) ,壁面平均温度则增加     

用于冲击/发散双层壁冷却数值模拟的源项法模型

为了避免冷却孔内的网格划分从而减小计算量,将源项法模型应用到冲击加发散双层壁冷却结构的数值模拟中,并分析了网格设置和负能量源离散方式的影响。结果表明：源项法模型能较准确地预测综合冷效的变化趋势;在本文计算范围内,当发散孔入口/出口处面上平均面网格数超过一定量时,数值模拟结果不再随入口/出口处面上网格密度而变化;发散孔热侧壁面附近第一层网格高度对计算结果影响较大,较合适的高度设置是小于或等于0.05倍发散孔直径;固体内负能量源对应源项应均匀地加载在多个体单元上。     

基于ANSYS的齿根应力和齿面变形的合理分析

根据齿轮啮合原理,准确求出复杂的齿廓曲线;采用MATALB—APDL混合建模方法,创建精确的齿轮有限元模型。在不同网格密度等级下,分别对齿轮的齿面变形、齿根应力分布进行分析。对照经验公式的计算结果,得出如下结论：在齿轮数值模拟中,由于有限元软件计算中舍入误差的存在,网格划分合理与否将直接影响仿真程度的高低;网格密度等级对齿面最大变形的仿真结果有较大影响,而对齿根应力影响甚微。     

基于Voronoi序列采样的加筋壁板优化设计

为提高加筋壁板结构轻质优化效率,提出一种基于Voronoi序列采样的加筋壁板优化设计方法.建立了该结构的参数化模型,分析了加筋壁板蒙皮与桁条腹板的网格规模对结构承载能力及失稳模式的影响规律,并依据桁条腹板高厚比进行了网格划分,进而平衡后屈曲分析精度与效率;然后,提出了基于探索策略和开发策略的序列近似优化方法,其中Vor...     

基于平方距离极小化方法用C-C细分曲面拟合三角网格

基于平方距离极小化方法(SDM),给出了用C-C细分曲面重构有特征的、任意拓扑三角网格模型的算法.首先识别特征并结合人工交互的方式在初始网格上进行四边域划分,然后直接在域面上对数据点进行近似参数化和分区,以域顶点为顶点构造初始控制网格.建立局部坐标系并优化每个数据点的参数值,基于SDM建立拟合方程.循环进行控制顶点的反求和待拟合数据点的参数校正,直至达到给定的误差要求.与传统的最小二乘法拟合(LSQ)相比,本文方法的逼近精度(二阶逼近)要高得多.     

一种解Dirichlet问题的虚拟区域法及其在广义Stokes问题中的应用(英文)

讨论了求解线性 Dirichlet问题的虚拟区域法及其在广义 Stokes问题中的应用 .这种方法通过 Lagrange乘子技术来处理 Dirichlet边界条件 ,因而非常适用于粘性流动问题中的无滑移边界条件。为提高解的精度 ,我们根据事后误差估算对网格进行自适应加密。Mini-element离散被应用于广义 Stokes问题的求解。最后 ,我们给出一些数值结果以验证这种针对 Dirichlet边界条件的偏微分程解法。     

二维应力集中问题的多尺度RKPM自适应分析

借助于“小波”理论,再生核质点方法RKPM(Reproducing Kernel Particle Method)以将形状函数及求得的结构响应分解为多个尺度。本文对线弹性二维应力集中问题进行了双尺度分解,并由各应力分量计算得到的高梯度点作为误差指示,实现了该方法的h型自适应分析。并且提出了一种新的方法——“四象限法”对高梯度区域进行加密,计算结果表明自适应后的解的精度更高,从而证明了这种自适应无网格方法的有效性。     

网格整体加筋贮箱圆筒壳结构优化设计

为完成网格整体加筋贮箱圆筒壳的结构优化设计,本文首先开展了网格整体加筋贮箱圆筒壳的有限元模型参数化研究;然后基于参数化模型,进行了网格加筋角度的优化设计和加筋尺寸与壁面厚度尺寸优化设计,讨论了网格尺寸和加筋角度与结构屈曲稳定性和质量之间的关系;并利用实验设计的方法构造响应面代理模型对网格尺寸和加筋角度进行了优化设计.通过以上参数化建模、代理模型构建和优化设计,得到了一套较为完整可行的网格整体加筋贮箱圆筒壳结构优化设计方法.最后还对研究过程中发现的问题进行了总结.     

AUFS格式在多组分反应流场模拟中的应用

对AUFS格式进行了讨论,并将其应用于高超声速多组分化学反应流场的数值模拟.超声速燃烧流场复杂性使得发展基于非结构网格的高精度逆风格式数值求解方法成为必然.本文中控制方程为轴对称多组分Eu-ler方程,空间分裂采用基于非结构网格的有限体积逆风格式AUFS方法,时间推进为考虑了化学反应特征时间效应的显式5步Runge-Kutta方法.化学反应模型采用经典的H2/Air的7组分8反应模型.对钝头体激波诱导亚爆轰及超爆轰流场进行了数值模拟,得到了与实验及文献相符合的计算结果.     

刚架模型驱动网格曲面保特征变形

将刚架结构受力变形原理应用于网格曲面的变形设计中,提出一种新的曲面造型方法。通过建立与网格拓扑关系一致的刚架模型,运用刚架模型驱动网格曲面变形:即对刚架节点施加载荷或添加几何约束反算出刚架节点载荷,在载荷作用下刚架变形驱动网格曲面变形。针对几何约束变形问题,建立以节点位移为变量,由载荷最小和节点位移最小组成的目标函数,组合使用两种不同的最小化模型,实现不同的变形需求。在变形过程中,采用增大网格曲面特征区域所关联单元体弹性模量的方法,达到整体变形中局部形状特征的保持,避免了传统方法中繁琐的特征信息记录公式。实验结果表明,该变形方法直观方便,可以产生理想的变形效果。     

基于网格划分的ADS-B地面站信号覆盖及选址分析

针对通用航空飞行活动的安全监视问题,合理的广播式自动相关监视（Automatic dependent surveillance-broadcast,ADS-B）地面站选址规划能够提高空域的监视范围,基于网格划分分析ADS-B地面站对真实高度的监视信号覆盖范围并进行选址分析。考虑地球曲率和地形遮蔽对信号传播距离的影响,采用数字地理高程数据建立地理模型,将空域网格化处理,用Xdraw算法计算网格覆盖性,确定ADS-B地面站对真实高度的最大理论覆盖范围,构建地面站最大覆盖模型进行空间规划布局。最后,对ADS-B地面站信号覆盖情况进行仿真,并以湖北省为例,计算出满足该区域覆盖需求的地面站布局。为实际工作中的地面站布局规划提供了依据。     

三维有限元网格的自动划分

本文作者从研究结构的数学模型入手,讨论了建立三维结构有限元计算模型的方法。成功地获得了一些具有应用代表性的三维板壳和块体结构有限元计算网格的自动划分原理和算法。     

矩形渡槽槽体结构绕流流场数值模拟研究

基于二维非定常不可压缩雷诺平均N-S方程,采用RNG(renormalization group)k-ε两方程湍流模型,对矩形渡槽在空槽和满槽情况下的二维绕流特性进行了分析.计算网格为四边形结构网格,采用有限体积法对微分方程进行离散,应用SIMPLE算法解决压强-速度耦合关系.引入空腔流动研究成果对空槽和满槽风载体形系数基本相同的原因进行了理论分析,揭示了空槽时槽体内复杂的旋涡流动结构.