幂律流体两平行圆板间径向扩散层流进口段流动阻力的分析

在文献￣［２，３］的基础上，将其方法推广到幂律流体流动，首先导出幂律流体两平行圆板间径向扩散层流边界层的动量积分方程和能量积分方程。通过对动量积分方程的求解，对幂律流体的进口段长度以及进口段效应压力损失和流量修正系数进行了分析计算，并且讨论了幂律流体的流动指数ｎ与进口处修正雷诺数Ｒｅ＿１对进口段长度和压力损失系数的影响。特别是当ｎ＝１时，本文的解与文献￣［３］中的解完全一致，因而间接验证了本文结果的可靠性。     

裂纹面无磨擦接触的弹塑性界面裂尖应力场

本文通过对弹塑性幂硬化双材料界面裂纹尖端应力场的高阶渐近分析，获得了裂纹面无磨擦接触的裂尖一阶和二阶应力场解答，位移场在界面处呈现交叉匹配是本文解答的一个重要特点，最后结果表明，当界面上下材料的硬化指数之差大于１时（即ｎ１－ｎ２＞１时），二阶应力场角分布为一常数解；而当０＜ｎ１－ｎ２≤１时，二阶应力场角分布函数则随θ变化而变化。     

裂纹面无摩擦接触的弹塑性界面裂尖应力场

本文通过对弹塑性幂硬化双材料界面裂纹尖端应力场的高阶渐近分析，获得了裂纹面无摩擦接触的裂尖一阶和二阶应力场解答，位移场在界面处呈现交叉匹配是本文解答的一个重要特点．最后结果表明，当界面上下材料的硬化指数之差大于１时，（即ｎ１－ｎ２＞１时），二阶应力场角分布为一常数解，而当０＜ｎ１－ｎ２≤１时，二阶应力场角分布函数则随θ变化而变化。     

幂软化损伤材料Ⅲ型动态裂纹尖端的渐近场

利用幂软化损伤模型，对Ⅲ型裂纹进行了详细的研究。给出了本构方程及渐近方程。位移、应变和应力用对数级数展开，揭示了场的渐近特性。其结果表明，应力和应变分别具有如下的奇异性：σ－（ｌｎＲ／ｒ）＾１／ｎ＋１，ε－（ＬｎＲ／ｒ）＾ｎ／ｎ＋１。     

梯形断面明渠丁坝绕流水力特性三维数值模拟

采用两相流混合模型,并选取RNG k-ε湍流模型封闭两相流时均方程,对梯型断面明渠非淹没式丁坝绕流水力特性进行了三维数值模拟。采用有限体积法离散计算区域,求解速度与压力耦合方程组时使用半隐式SIMPLE(Semi-Implicit Method for Pressure-Linked Equations)算法,模拟自由水面时采用了VOF(Volume of Fluid)法。对丁坝后不同的回流长度进行了分析比较,并将模拟结果与实测资料进行了对比验证,结果表明两者吻合较好,相对误差小于8%,说明该模型能够很好地模拟明渠丁坝绕流的水力特性分布规律。     

延时神经网络混沌现象的研究

本文研究一个具有部分延时性能的三个神经元的不对称神经网络系统的动力学行为，讨论了当控制参数变化时网络呈现的各种现象，计算了李雅普人指数研究了阵发混沌现象，最后讨论了全延时神经网络的吸引子。     

应变损伤材料I型动态扩展裂纹尖端场的渐近分析

研究了应变损伤材料Ｉ型动态扩展的裂纹尖端场。假定材料服从Ｊ２流动理论，且损伤规律以幂律应变软化的规律给出。对于塑性区引进了应力函数φ，ψ０借助于动力学方程的分析，给出了渐近方程及数值解。结果表明，对于可压缩材料Ｉ型平面应变尖端场是完全由塑性区组成，没有弹性卸载区。在裂纹尖端附近，应力和应变分别具有如下的奇异性：σ￣（ｌｎＲ／ｒ）＾－ｎ／ｎ＋１，ε￣（ｌｎＲ／ｒ）＾１／ｎ＋１。     

弯曲井眼中受压管柱的屈曲分析

本文通过微元体的变形和受力分析，导出了弯曲井眼中受压管柱的屈曲方程──一个合参数。的四阶非线性常微分方程。ε是一个综合考虑轴压（Ｆ０）、管柱截面抗弯刚度（ＥＪ）、井眼轴线轨迹曲车半径（Ｒ）以及共眼有效半径（ｒ）的无因次参数。本文利用线性化方法及小参数摄动法求得了屈曲方程的一次分叉点ε１及二次分叉点ε２。当ε＞ε１时，管柱处于稳定状态，对应于屈曲方程的零解；当ε２＜ε≤ε１时，管柱处于正弦屈曲状态，对应于屈曲方程的周期解；当ε≤ε２时，管柱处于螺旋屈曲状态，对应于屈曲方程的螺线解。本文给出的弯曲井眼中受压管柱的两个临界屈曲载荷（与两个分叉声、对应），发生螺旋屈曲时的螺距等理论结果与他人发表的实验结果吻合较好。     

相同雷诺数下串列双圆柱绕流的仿真分析

为研究均匀水流场中串列排布的柱群之间的干涉影响,本文以三维串列双圆柱为例,通过计算流体力学(CFD)软件FLUENT15.0中双方程k-ε模型,分析模拟了双圆柱所受平均阻力、平均升力、后柱周向压力、斯特劳哈尔数等水动力特性。结果表明:在雷诺数为Re=2×10~4的串列双圆柱绕流中,两圆柱中心间距L与圆柱直径D的比值为L/D=4时,后柱受前柱绕流尾流影响大,明显高于单圆柱绕流的平均阻力系数,后柱的周向压力值也随前柱尾流的摆动呈现显著的不对称性;当L/D=8时,前柱绕流尾流对后柱影响逐渐减弱;当L/D=12时,两圆柱之间的相互干扰几乎可以忽略,可以看作是相互独立的单圆柱绕流。最后,计算的斯特劳哈尔数与单圆柱绕流对应的斯特劳哈尔数相近且仿真数值在计算数值范围之内,验证了整个仿真分析的准确性,也进一步说明了双圆柱绕流的柱群的干涉影响。双圆柱间距越大,前、后柱之间的干涉影响越弱。     

轻载径向滑动轴承中Taylor涡动的产生和影响研究

本文用原始变量法直接求解了三维的Ｎ－Ｓ方程，计算分析了高速旋转有限长圆柱轴承中油膜层流失稳产生涡动的临界Ｔａｙｌｏｒ数及流场、压力场和摩擦阻力。轴承端部泄油量等的变化。结果表明，在有限长同心圆柱轴承中，随着轴旋转速度的提高，轴承磨擦阻力线性增大，油膜层流失稳出现的涡动增加轴承摩擦阻力并减少轴承端部泄油量，油膜层流失稳后，轴承长度方向均匀地排列着一些流体涡，涡动的强度从轴承中间截面向轴承端部逐渐减弱     

光滑拉伸试件中不同初始形状孔洞长大的有限元模拟

本文对具有不同硬化指数（ｎ＝０．０５，ｎ＝０．１，ｎ＝０．２）的幂硬化材料的光滑拉伸试样的拉伸变形过程进行了有限元模拟，通过有限元计算和经Ｂｒｉｄｇｍａｎ修正分别得到了试样变形过程中心部应力三维度随应变的变化情况；在此基础上运用控制体胞宏观应力三维度的方法，对含不同初始形状孔洞的体胞模型进行了有限元分析，计算结果表明：（１）孔洞初始形状和材料硬化特性对试样的拉伸破坏过程有重要影响；（２）Ｂｒｉｄｇ     

统一二阶矩模型用于模拟旋流湍流两相流动

用统一二阶矩模型（ＵＳＭ）模拟了旋流数为０４７和１５的气粒两相流动，并和实验结果以及ｋ ε ｋｐ模型的模拟结果进行了对比．研究结果表明，提高旋流数减小了轴向速度反流区，增大了切向速度似固核区．ＵＳＭ和ｋ ε ｋｐ模型预报旋流数为０４７时的两相速度场差别不大，并都和实验结果接近，但前者预报的旋流数为１５的两相速度场比后者有改进，在两种情况下，前者都能揭示出后者无法预报的两相湍流各向异性规律．     

准周期激励Duffing系统中奇异非混沌吸引子的判别方法

奇异非混沌动力学是非线性动力学领域中的新课题。本文以准周期激励Duffing振子为例,对其产生的奇异非混沌吸引子(strange nonchaotic attractors, SNAs)进行分析。通过三维庞加莱截面和定量方法如傅里叶变换、李雅普诺夫指数、李雅普诺夫维数、关联维数和盒维数检测SNAs是否存在。研究结果表明,傅里叶变换无法判断混沌与奇异非混沌行为。而李雅普诺夫指数、李雅普诺夫维数可以作为检测系统混沌与非混沌指标。关联维数和盒维数显著表明系统奇异与非奇异性,从而阐明适用于准周期驱动Duffing振子中存在SNAs的判别方法,并为其他类似系统检测SNAs提供指导。     

超高雷诺数下矩形截面气动特性研究

为研究工程常见的超高雷诺数下矩形柱流场特性和气动特性,以常见的起重机单箱梁截面(梁高h=1m,雷诺数Re=6.85×105)为例,采用雷诺平均RANS Realizable k-ε方法,对中等湍流强度来流下矩形柱绕流进行数值模拟。首先进行网格无关性验证,然后选取宽高比为0.6、0.8、1.0、1.5、2.0共5组不同截面的矩形柱结构进行流场分析与研究。通过对比不同宽高比柱体截面周围的气动特性参数,以及尾流区的尾涡、回流参数等特征量表明:阻力系数、平均背压系数、升力系数均方根、斯托罗哈数等均随着宽高比增大而降低,尾涡宽度和距离随着宽高比增大而减小,回流长度和最大回流速度随宽高比增大而增大。同时,通过与相关雷诺数(Re=104～105)的风洞实验和仿真结果进行对比,发现其在数值上和参考值相差很小,说明高雷诺数对绕流参数的影响很小,可为相关工程设计提供参考。     

圆柱绕流流场结构的大涡模拟研究

为进一步揭示绕流现象的形成机理,本文分别对处于层流稳态区、尾流过渡区、剪切层转换区Re分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

泥沙运动的混沌吸引子

本文应用河流泥沙运动的实测资料以及重构相空间的方法，研究了河流泥沙运动的混沌特征。结果发现在河流泥沙运动中混沌吸引子是存在的。就本文给出的实例而言，它们的关联维数和最大的李雅普诺夫指数分别为６．６和０．０１３左右，上述结果对河流泥沙运动模式的评估，泥沙观测方案的设计以及河流泥沙可预测性问题的研究是有重要参考价值的。     

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

非定常不可压N-S方程的最小二乘算子分裂有限元法数值求解

采用最小二乘算子分裂有限元法求解非定常不可压N-S(Navier-Stokes)方程,即在每个时间层上采用算子分裂法将N-S方程分裂成扩散项和对流项,这样既能考虑对流占优特点又能顾及方程的扩散性质。扩散项是一个抛物型方程,时间离散采用向后差分格式,空间离散采用标准Galerkin有限元法。对流项的时间项采用后向差分格式,非线性部分用牛顿法进行线性化处理,再用最小二乘有限元法进行空间离散,得到对称正定的代数方程组系数矩阵。采用Re=1000的方腔流对该算法的有效性进行检验,表明其具有较高的精度,能够很好地捕捉流场中的涡结构。同时,对圆柱层流绕流进行了数值研究,通过流线图、压力场、阻力系数、升力系数及斯特劳哈数等结果的分析与对比,表明本文算法对于模拟圆柱层流绕流是准确和可靠的。     

亚临界雷诺数下圆柱绕流的大涡模拟

本文应用Smagorinsky涡粘性模式和二阶精度的有限体积法对圆柱绕流的流场进行大涡模拟．求解了非正交曲线坐标系下的N-S方程,对雷诺数为100和20000的工况进行了计算．计算结果与实验及动力涡粘性模式的结果进行了比较,表明计算对于层流及高亚临界雷诺数的湍流流动是合理的     

血液入口段两相流动的层流NS方程数值模拟

本文采用球形刚性颗粒悬浮于牛顿流体的血液模型，基于有相间滑移的两相流动层流Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，采用相间滑移算法（Ｉｎｔｅｒ－ＰｈａｓｅＳｌｉｐＡｌｇｏｒｉｔｈｍ）对圆管内定常轴对称血液入口两相流动问题进行了计算，获得了与实验结果吻合的计算结果，较好地模拟血液流动中红细胞的径向迁移现象，结果表明，采用二相流动模型研究血液流动是一种有前途的方法。