对两端支承圆柱体运动方程的简化

研究了两端支承圆柱体的稳定性和失稳振动.在方程中忽略了轴向运动和横向运动的高阶微量,从而进一步简化了两端支承圆柱体的运动微分方程.对六阶离散化后的方程进行了数值模拟,发现系统在u=14.28时发生颤振失稳,发生颤振失稳的平衡位置为q1=0.0946,与Y.Modarres-Sadeghi等人于2005年得到的结果基本吻合,说明本文的简化方法是合理的.     

非线性负阻尼极限环型角位移振荡运动定量分析

从某飞行器两次飞行记录的俯仰角θ、偏航角Ψ观测值出发,分别采用对Van der Pol方程描述的θ和Ψ的角位移振荡运动微分方程进行气动参数辨识,和对观测值θ-狋、Ψ-t曲线的外包络线进行参数拟合两种分析方法,对该飞行器的角位移振荡运动特性进行了定量分析。两种分析方法取得了较为一致的结果,并证明该飞行器飞行中出现的锥形振荡运动是典型的非线性负阻尼极限环型的振荡运动,获得了非线性负阻尼极限环型气动阻尼力矩的典型表达式。分析结果表明非线性负阻尼极限环类型的气动阻尼能够导致飞行器出现动不稳定。     

在速度平方控制项下可变方向悬臂梁的振动

本文讨论研究可改变方向的悬臂梁非线性系统参激振动并利用速度二次项对其振动控制。讨论中认为末端带集中质量块悬臂梁的第一阶模态具有决定性的影响，对系统位移响应采用单阶Galerkin截断方法处理将描述梁的偏微分方程转化为常微分方程。本文应用非线性梁方程中加入二次速度反馈项进行控制，对控制方程应用多尺度方法得到幅频关系并作稳定性分析，最后用数值方法分别讨论外谐振下线性阻尼项、梁非线性系项、二次速度控制项增益和激励幅对悬臂梁振幅在不同方向下的影响。     

简谐激励下弯曲薄板的非线性运动方程

推导了两端固定和两端简支的弯曲薄板在其支承简谐激励下的运动微分方程,并利用Galerkin方法对其进行了离散化处理。     

简谐激励下输流管道的非线性运动方程

本文推导了两端受扭转弹簧约束的简支输流管道在其支承简谐运动激励下的运动微分方程，并利用Ｇａｌｅｒｋｉｎ方法对其进行离散化处理。另外，针对于扭转弹簧刚度的两种极限情形①Ｋ１＝０，Ｋ２＝０（两端简支）和②Ｋ１→∞，Ｋ２→∞（两端固定）进行了讨论，还将运动方程进行了简化。     

基于三角剪切变形理论的对称层合板振动分析

根据三角剪切变形层合梁理论,推导出对称复合材料层合板的运动微分方程,利用逆复合二次径向基函数无网格配点法对运动微分方程进行离散,预测了对称复合材料层合板的自由振动特性。将不同材料参数、几何尺寸的层合板固有频率计算结果与相关文献中的结果进行对比,结果表明：用逆复合二次径向基函数离散的三角剪切变形理论在对称复合材料层合板自由振动分析方面具有一定的精度。     

第二类拉格朗日方程建模在自动武器上的应用

以高等动力学理论为基础,采用第二类拉格朗日方程建模方法,推导出典型二自由度机构运动微分方程,给出用该方法推导的某航空机关炮二自由度运动微分方程,并对模拟结果进行了分析。     

轴向力作用两端支承梁的固有特性和稳定性分析

推导了轴向力作用下的两端支承梁的运动微分方程,以两端支承梁的振型函数为假设振型函数,分析了不同扭转弹簧刚度下轴向力变化对梁固有频率的影响,导出了两端支承梁在临界轴向力作用下的解析表达式,讨论了扭转弹簧刚度对梁临界轴向力的影响特性.     

飞机起落架着陆撞击动力分析

根据现有油—气式缓冲器结构形式．针对模拟起落架着陆撞击过程的落震试验模型．考虑起落架航向刚度和阻尼的影响,依据撞击过程中各变量间的运动关系．建立了描述着陆过程中的运动微分方程组．并采用龙格库塔法对运动微分方程组进行求解。最后采用该方法对K8飞机主起落架(单腔)和某机起落架(双腔)落震试验过程进行了模拟．并与试验结果进行了比较。发现计算结果与试验结果吻合较好。     

两端固定输流管道在脉动内流作用下参数共振的最优控制

采用压电片作为控制执行元件,对两端固定输流管道在脉动流作用下的参数共振进行了最优控制.基于二次型性能指标最优控制原理,针对离散后的运动方程在零平衡点附近Jacobian矩阵的非时变部分设计了最优控制器.最后,通过数值仿真的方法对控制方案的有效性进行了验证.结果表明,该方案对输流管道在第一、第二振型次谐波共振区域及组合共振区域的参数共振的控制均有很好的效果.     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

飞行器自激振荡的流动物理与动力学机制

全局亚迭代耦合求解流体动力学方程和刚体动力学方程,研究轴对称飞行器单自由度俯仰失稳运动的非定常特征.高超音速锥-柱-裙飞行器的单自由度俯仰失稳运动发展为极限环形式,周期性运动伴随着波系结构的非定常变化;自激振荡源于静稳定构型的振动回复机制和尾裙激波流动迟滞效应所构成的阻尼机制.基于第二拉格朗日方程和虚功原理,导出能够描述迟滞现象的参数化非线性动力学模型.多尺度近似分析获得参数化运动特征:自激振动过程是拟简谐运动;静稳定构型才能出现自激振荡周期运动;平衡点阻尼是决定运动稳定特性的分叉参数;振幅特性与阻尼非线性相关,频率特性与刚度非线性相关.基于数值结果的参数辨识,非线性模型的理论分析、重构都与数值结果高度一致,从而有效地佐证了自激振荡建模研究的合理性.     

三维分块非结构化网格上的可压缩流计算方法

提出了分块非结构化网格生成方法,通过子块划分、子块网格划分、整体网格拓扑连接生成可局部加密的分块非结构化网格,既具备了非结构化网格对复杂区域的处理能力,又避免了分块结构化网格块间信息传递的复杂计算过程.通过将密度修正值引入压力修正方程,发展了非结构化网格上可压缩流动的SIMPLE算法.对孤立翼型及透平叶栅的流动进行了数值计算,压力系数、速度、折转角的数值结果与已有的实验值吻合较好,验证了数值方法的精度.      

飞行器失稳平面振荡运动的物理机制

基于全局亚迭代耦合求解非定常流体动力学方程和刚体动力学方程(CFD/RBD),研究动不稳定飞行器在自由俯仰与自由沉浮二自由度下自激发平面失稳运动的非定常特征。数值研究表明:超声速锥-柱-裙飞行器的平面失稳运动发展为极限环形式,并伴随着波系结构非定常变化;平面运动保持了自由俯仰基本运动特征,但同步自由沉浮使得极限环周期运动的振幅更小、频率更快;平面自由运动中飞行器绕靠近头部的"不动点"转动。基于第二拉格朗日方程和虚功原理,导出能够描述迟滞现象的参数化非线性动力学模型。多尺度近似分析(MTS)获得参数化运动特征:自激振动过程是拟简谐运动;平衡点阻尼是决定运动稳定特性的分叉参数;振幅特性与阻尼非线性相关,频率特性与刚度非线性相关;模型分析证实了平面自由运动的"不动点"现象并自洽地解释了沉浮自由度存在使得极限环振幅变小的动力学机制。非线性模型的理论分析、重构都与数值结果高度一致,从而有效地佐证了自激振荡建模研究的合理性。     

三维数值模拟再生冷却喷管的换热

为了解液体火箭发动机喷管再生冷却的换热特点，采用数值模拟的方法，对内喷管燃气、壁面和冷却液建立不同的三维控制方程，进行流动和传热的耦合计算。在计算中，假定喷管流动为冻结流动，考虑燃气向壁面的对流换热和辐射换热；采用二阶迎风格式离散控制方程，采用DO模型离散求解辐射换热方程，水蒸气的吸收系数根据Leckner公式计算。计算模型采用缩比热试车发动机，数值计算结果与实验结果吻合较好，较准确地模拟出了喷管的壁面热流密度，得到了喷管燃气和冷却液的流场和温度场，对高压再生冷却喷管的设计具有指导意义。     

高温下加热/冷却管内参与性介质的辐射与对流耦合换热研究

本文通过数值模拟方法研究了高温下加热/冷却管内层流充分发展流动的参与性介质热起始段的辐射与对流耦合换热。采用离散坐标法求解圆柱坐标系下的非线性积分微分辐射传递方程,采用有限差分法求解辐射与对流的耦合换热能量方程。通过计算,分析了不同壁温及入口温度下,介质导热系数、壁面发射率、介质散射反照率及散射相函数等因素对局部努谢尔数、壁面热流密度及有效加热/冷却长度的影响,并同时采用努谢尔数和壁面热流密度评价了耦合换热中辐射对对流换热的影响。      

A second-order numerical method for elliptic equations with singular sources using local flter

The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic equations with singular sources is introduced by employing a local kernel flter.In this method,the discontinuous equation is convoluted with the kernel function to obtain a more regular one.Then the original equation is replaced by this fltered equation around the singular points,to obtain discrete numerical form.The unchanged equations at the other points are discretized by using a central difference scheme.1D and 2D examples are carried out to validate the correctness and accuracy of the present method.The results show that a second-order of accuracy can be obtained in the fltering framework with an appropriate integration rule.Furthermore,the present method does not need any jump condition,and also has extremely simple form that can be easily extended to high dimensional cases and complex geometry.     

匀速及加速行驶时带后视镜轿车气动噪声的研究

应用计算流体力学软件FLUENT对带后视镜简化直背式轿车外部流场及车外某接收点处的气动噪声进行了数值计算和分析.脉动流场计算采用大涡模拟湍流模型,动量方程的离散格式采用有限中心差分格式.分别对两种匀速运动情况和一种加速运动进行了计算,由此研究不同车速、匀速和加速运动对接收点处脉动流场和气动噪声的影响.对于匀速运动,车速越大,接收点处流速越大,压力越小,脉动压力的脉动幅度越大,在每个频率上所对应的声压级基本上越大,总声压级也越大;在同一行驶速度下,加速运动比匀速运动在接收点处的流速要大,压力要小,在每个频率上所对应的声压级基本上要大,总声压级要大,声压级随频率的变化情况与匀速运动一样.      

边界元法分析回转体的自由扭振问题

 边界元法用于分析非线性问题和动力问题时,在最终得到的积分方程中同时包含边界积分和区域积分。通常的做法是对区域积分项采用常值有限元离散,这需要用很密的网格才能达到较高的精度。