一种无条件稳定的时域推进法与有限元混合分析瞬态热传导问题

本文将一种新的无条件稳定的时域推进法与有限元结合,用于分析瞬态热传导问题。根据微-积分型热传导控制方程,对时间变量在小区闻内插值积分,变初-边值问题为一系列离散时刻的边界值问题,再应用有限元法求解之。由于推进中的每一时刻的解,都严格满足原给定的初始条件,就消除了累积误差的影响。此外,这种时域推进法对时间变量的数值积分也优于数值微分。因而,与现行的基于微分型控制方程,采用时-空有限元模型同时离散,再逐步求解的直接积分算法相比,计算精度可大大提高。不难预料,当求解较长时间后的瞬态值时,本算法的优越性会更加明显。     

固体火箭复合结构喷管传热的隐式计算方法

采用适用于多层复合结构的当量热容瞬态热传导方程,通过坐标变换解决多个移动边界问题,用ADI方法推导了二维隐式求解格式,在与烧蚀相耦合的情况下获得了复合结构喷管的温度分布。用隐式方法计算有效地减少了计算机时,对工程设计大有裨益。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

叶轮机转子内流动的三维跨声Navier-Stokes方程组数值计算

应用有限体积法及Baldwin-Lomax代数紊流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛。采用隐式算法对对流通量进行求解,而扩散通量采用显式算法。采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算扩散通量。关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解。应用该算法模拟了转子NASARotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等等。在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验略有差别。      

叶轮机转子内流动的三维跨声速Navier-Stokes方程组数值计算及试验验证

应用有限体积法及Baldwin-Lomax代数湍流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛.采用隐式算法对对流通量进行求解,而耗散通量采用显式算法.采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算耗散通量.关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解.应用该算法模拟了转子Rotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等.在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验有较大的差别.     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

基于HDCS-E8T7格式的气动噪声数值模拟方法研究进展

发展了气动噪声高精度数值模拟方法,空间导数离散采用了HDCS-E8T7格式及精度与之匹配的边界格式,时间离散方法采用了高精度多步龙格库塔方法和高精度隐式双时间步方法,发展了高精度对接边界算法,将方法推广至多块对接网格以满足解决复杂几何构型问题的需要,采用隐式大涡的概念处理可能出现的湍流问题。在此基础上,研究了几何守恒律对计算结果的影响,展示了复杂网格中高精度计算满足几何守恒律的重要性,完成了等熵涡、双圆柱散射、串列柱翼构型和喷嘴射流等典型噪声问题的求解,所得计算结果展示了所发展的模拟方法具有良好的预测精度和解决复杂构型气动噪声问题的潜力。     

篦齿容积效应对盘腔空气系统瞬态特性的影响

针对航空发动机中空气系统的瞬态响应问题,以多级封严盘腔为研究对象,开展了考虑篦齿容积效应下的空气系统瞬态计算分析研究.在研究中,采用一维瞬态流体网络计算方法,重点分析了封严篦齿和旋转盘腔的容积效应,并建立了相应的一维瞬态数学模型,通过迭代求解瞬态网络守恒方程组,得到了不同边界工况下盘腔内流体参数随时间响应规律.研究结果...     

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响      

固体发动机药柱热线粘弹性边界元分析

采用遗传积分形式的本构关系,推导了热线粘弹性材料时域边界元公式。由于把与时间有关的剪切模量G(t)表示成Pnony级数的形式,因而可用递推公式较方便地计算粘弹性记忆应力。本文提出的方法比以前的对应原理方法有较广的应用范围。为了验证提出的方法,编制了二维粘弹性边界元程序。两个算例的数值结果表明:本方法具有精度高,计算机解题快,占内存少等优点。     

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分。同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式。对静止介质中的二维圆柱声散射,数值解在研究的频率与观察点处与理论解一致。采用数值方法计算二维圆柱绕流诱发的半空间声场基本解,结果显示半空间边界强化了声源的声辐射能力,小于马赫数0.2的介质运动对声传播的影响可以忽略。     

SPACE-TIME CONSERVATION SCHEME IN NONORTHOGONAL CURVILINEAR COORDINATES

ＳＰＡＣＥ┐ＴＩＭＥＣＯＮＳＥＲＶＡＴＩＯＮＳＣＨＥＭＥＩＮＮＯＮＯＲＴＨＯＧＯＮＡＬＣＵＲＶＩＬＩＮＥＡＲＣＯＯＲＤＩＮＡＴＥＳＬｉｕＨａｉｔａｏ（刘海涛）,ＸｕＪｉａｎｚｈｏｎｇ（徐建中）,ＨｕａｎｇＸｉｕｑｉａｎ（黄修乾）（Ｉｎｓｔｉｔｕｔｅｏ...     

基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型.该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在两者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解.这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率.该模型先得到了Levine-Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响.      

Boundary discontinuous Fourier analysis of thick beams with clamped and simply supported edges via CUF

This paper presents an analytical solution for static analysis of thick rectangular beams with different boundary conditions.Carrera's Unified Formulation (CUF) is used in order to consider shear deformation theories of arbitrary order.The novelty of the present work is that a boundary discontinuous Fourier approach is used to consider clamped boundary conditions in the analytical solution,unlike Navier-type solutions which are restricted to simply supported beams.Governing equations are obtained by employing the principle of virtual work.The numerical accuracy of results is ascertained by studying the convergence of the solution and comparing the results to those of a 3D finite element solution.Beams subjected to bending due to a uniform pressure load and subjected to torsion due to opposite linear forces are considered.Overall,accurate results close to those of 3D finite element solutions are obtained,which can be used to validate finite element results or other approximate methods.     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

复合材料层合结构频率约束的可靠性优化设计(英文)

The reliability based optimization （RBO） issue of composite laminates trader fundamental frequency constraint is studied. Considering the tmcertainties of material properties, the frequency constraint reliability of the structure is evaluated by the combination of response surface method （RSM） and finite element method. An optimization algorithm is developed based on the mechanism of laminate frequency characteristics, to optimize the laminate in terms of the ply amount and orientation angles. Numerical examples of composite laminates and cylindrical shell illustrate the advantages of the present optimization algorithm on the efficiency and applicability respects. The optimal solutions of RBO are obviously different from the deterministic optimization results, and the necessity of considering material property uncertainties in the composite structural frequency constraint optimization is revealed.     

Hoff型夹层板振动分析的边界积分方程法

 本文利用静力问题的基本解建立了Hoff型夹层板振动分析的边界积分方程。在数值实现过程中,边界上和区域内的未知函数均釆用线性插值以提高计算精度。数值结果表明本文方法具有未知量少和计算精度高等优点。对于低阶频率,本文结果与精确解的误差不大于1％。     

航天器贮箱变质量流固耦合系统的动力学响应

研究了航天器贮箱变质量流固耦合系统的动力学特性.根据虚拟质量法(VMM),结合边界元法和有限元法构建了系统的动力学模型,建模过程中着重考虑了质量变化对系统的动态影响.通过Newmark直接积分法计算出贮箱变质量系统的振动响应.结果表明:由于系统质量减少,引起了系统振动频率的增大,并产生一个附加负阻尼.系统的振动频率的范围可以通过系统质量的范围确定.变质量引起的附加负阻尼的大小与系统的质量变化率成正比.对于系统的横向振动,质量减少引起的附加负阻尼在整个过程对系统振动的影响比较稳定,对于系统的纵向振动,质量减少引起的附加负阻尼对系统振动的影响随时间的增加而增大.