粘弹性板的蠕变屈曲研究

研究了有初始小挠度单向受压粘弹性板在两种民政部下的蠕变屈曲问题。第一种情况是建立控制方程时忽略横向剪切变形；第二种情况是不忽略横向剪切变形。通过求解控制方程得到了挠度随时间扩的解析表达式，从而又得到了瞬时临界载荷和持久临界载荷，并讨论了横向剪切变形对临界载荷的影响。     

正交各向异性粘弹性薄板的蠕变屈曲研究

研究有初始挠度的单向受压粘弹性板的蠕变屈曲问题，在建立控制方程时利用了经典薄板理论。通过求解控制方程得到挠度随时间变化的解析表达式，又通过分析得出了具有不同长宽比的粘弹性薄板的瞬时临界载荷及持久临界载荷。     

柔性系绳平衡状态临界绳长研究

为了确定柔性系绳系统在平衡状态下的临界绳长,在考虑大气阻力、系绳伸展性、系绳质量等因素的情况下,基于系绳系统轨道运动建立了系绳平衡状态方程.为确定求解系绳平衡方程所需初值条件,分析了载荷受力情况并建立其平衡状态方程,进而分析给出了系绳长度达到临界值的条件以及临界绳长的计算方法.仿真计算了不同载荷质量、系绳直径和载荷弹道系数这三种情况下系绳的临界绳长,分别拟合得到这三种情况下临界绳长的计算方程,并分析了不同载荷质量、系绳直径和载荷弹道系数对系绳临界绳长的影响.      

含椭圆分层缝纫复合材料层板的局部屈曲研究

建立了分析含椭圆分层的缝纫增强复合材料层板压缩屈曲问题的三维模型.通过应用瑞利-里兹法研究了分层椭圆几何参数和缝线等效弹性系数对压缩屈曲载荷的影响.分析结果表明,缝纫对提高临界曲屈载荷的作用十分明显,随着缝线等效弹性系数的增大,临界曲屈载荷迅速提高;分层椭圆的长轴垂直于载荷方向时,临界曲屈载荷最小,并且随长轴的长度增大时,临界曲屈载荷减小并逐渐收敛.     

板料不均匀拉伸起皱临界条件的分析

用弹性稳定理论能量法求得了一对边固支另一对边简支矩形板在任意载荷作用下失稳的临界方程,得到了板料不均匀拉伸诱导压应力起皱的临界载荷系数,并以此为基础,证明了方板对角拉伸与矩形板不均匀拉伸失稳临界点的一致性,以及方板对角拉伸试验(YBT)作为不均匀拉伸起皱现象模拟试验的可行性。     

用区间牛顿法跟踪结构的后屈曲路径

结构后屈曲分析传统上都是采用弧长控制法来处理.考虑到区间牛顿法求解非线性方程组的全局收敛性质,本文建议将这种方法用于跟踪结构后屈曲平衡路径.文内建立了将区间牛顿法用于有限元控制方程的增量求解迭代格式.在低载荷水平下,为减少迭代次数采用了小区间半径.当时间步接近临界点时,方法可以自动搜索出屈曲前和屈曲后的2个平衡构形,然后在屈曲后构形基础上跟踪临界点后的路径,这样就可以避免在临界载荷处由于刚度矩阵奇异而使迭代不能继续的数值困难.典型的数值算例证实了新方法的有效性,同时也表明算法能毫不费力地处理所谓"位移回弹"(snap-back)问题.     

复合材料夹层结构的逐渐破坏计算模型

建立了可以预测面内破坏、分层破坏和芯层破坏的复合材料夹层结构逐渐破坏模型.利用能量变分原理和非线性几何方程,以反对称失稳为例导出了考虑剪切效应夹层板的非线性稳定性控制方程组,并利用广义傅氏级数法对此高阶偏微分非线性方程组进行数值求解,在此基础上利用蔡-胡及最大剪应力强度理论及刚度退化规律得到复合材料夹层板的临界载荷、后屈曲路径、破坏过程和强度极限.      

薄板的后屈曲损伤分析与疲劳寿命预估

针对薄板在面内压缩载荷作用下的后屈曲损伤问题进行了研究,并进一步考虑屈曲与疲劳损伤的耦合作用,预估了薄板的疲劳寿命.首先建立了薄板的有限元模型,通过线性屈曲分析得到屈曲临界载荷和屈曲模态,进而采用大变形理论,将线性屈曲的一阶屈曲模态作为初始位移扰动,进行薄板的非线性屈曲分析,得到屈曲临界载荷.其次,根据损伤力学理论与方法建立了薄板材料在单次加载过程中的损伤演化方程,并根据材料疲劳试验结果进行参数识别,获取损伤演化参数.根据非线性屈曲分析结果和损伤演化方程进行了后屈曲损伤分析.最后,考虑疲劳载荷的作用,基于损伤力学理论,采用有限元数值方法求解,考虑每次加载引起的损伤与后屈曲应力应变场分析的耦合作用,通过反复迭代计算,给出了结构疲劳寿命.本研究为工程结构的后屈曲损伤分析以及考虑后屈曲损伤的疲劳寿命分析提供了一种新方法和实现手段.      

曲梁弹性波的传播

用特征线法求出了脉冲载荷作用下的控制曲梁运动的Ｍｏｒｌｅｙ方程的数值解，并用Ｐｌａｓｓ的直梁结果验证了曲率趋于零的极限情形，应用性地研究了曲梁中部受冲击载荷作用的弹性波传播问题，计算结果表明：该数值计算方法简便易行、效率高、有良好的收敛性和稳定性。     

复合材料层板在弹性支持下的非线性稳定性

利用能量变分原理和非线性几何方程，建立具有弹性约束的复合材料层板，在面内载荷作用下的非线性稳定性控制方程组，通过运用广义傅里叶级数法对控制方程组进行求解，得到相应的截荷－中心挠度曲线。计算结果表明：非对称层合板在面内载荷作用下，一般产生的是向变形，但在一定的铺层顺序，边界条件、受载方式等情况下，也可能发生分叉型屈曲。     

Hydrothermal instability of thermocapillary convection in large-Prandtl-number liquid bridges under microgravity

Linear stability analysis was performed to study the mechanism of transition of thermocapillary convection in liquid bridges with liquid volume ratios ranging from 0.4 to 1.2, aspect ratio of 0.75 and Prandtl number of 100. 2-D governing equations were solved to obtain the steady axi-symmetric basic flow and temperature distributions. 3-D perturbation equations were discretized at the collocation grid points using the Chebyshev-collocation method. Eigenvalues and eigenfunctions were obtained by using the Q–R method.The predicted critical Marangoni numbers and critical frequencies were compared with data from space experiments. The disturbance of the temperature distribution on the free surface causes the onset of oscillatory convection. It is shown that the origin of instability is related to the hydrothermal origin for convections in large-Prandtl-number liquid bridges.     

基于润滑理论的二维积冰数值模拟

将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统Messinger模型的模拟结果以及冰风洞试验结果进行对比,低温明冰条件下采用当前方法计算的翼型表面冰形接近于Messinger模型的模拟冰形曲线;相对高温条件下的计算结果比传统Messinger模型更为精确.低速积冰环境下Messinger模型难以模拟的传输线表面积冰同样可以采用水膜-冰层模型进行有效 预测.     

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

软管锥套式空中加油系统建模与特性分析

针对自主空中加油研究中软管锥套运动模型过于简化,现有建模方法又存在成本高、计算量大、软管长度恒定等缺陷的问题,根据集中参数法原理,提出了一种长度可变的多级串联理想单摆系软管锥套运动模型.同时考虑软管收放、加油机牵连运动、重力、定常流、大气扰动、加油机尾流等内外部因素,推导了迭代形式的变长度软管锥套三维运动方程.由摆长约束导出了求解软管拉力的代数线性方程组,进而给出了模型稳定性证明和适用条件.通过数值仿真,测试了锥套阻力、平衡位置、软管收放等稳态特性,分析了加油机滚转运动、尾涡流场对软管锥套的动态影响以及软管甩鞭现象的产生机理.数值仿真结果验证了所建模型的有效性.      

火焰筒结构局部热屈曲分析方法

根据局部热屈曲的失效机理,考虑局部温度和压应力水平两个热屈曲主要因素,从以应力函数和挠度函数为变量的受热板控制方程出发,推导出薄壁结构热屈曲失稳边界的判据表达式.依据确定的热载荷和设计的压力载荷,确定火焰筒结构的局部热屈曲危险部位,并将子模型技术引入到热屈曲分析中,确定危险部位的临界温度和临界压应力,以此对火焰筒结构进行局部热屈曲判别.通过工程实例的初步验证,计算与试验结果吻合良好,对火焰筒结构设计具有重要的工程应用价值.     

空间运动体上梁的三维动力学建模和仿真

对附着在空间运动体上梁的动力学问题进行了研究,运用Lagrange方法建立了中心刚体作任意三维运动时,梁作横向振动和纵向振动的动力学方程,考虑了柔性梁的横向弯曲变形所带来的纵向位移,这种耦合变形使所得到的振动方程包含了动力刚化效应,这是与传统的梁振动模型的不同之处,最后通过仿真算例验证了该方法.