矩形薄膜和充气管的屈曲及后屈曲行为分析

 对矩形薄膜及充气管的屈曲及后屈曲过程进行了分析。首先基于矩形薄膜褶皱区域的应力平衡关系,得到了矩形剪切薄膜褶皱的波长及幅度的表达式,并考虑了微小的初始张力的作用。基于稳定性理论建立了薄膜褶皱的数值分析方法,通过实验的方法对分析结果进行了验证。分析表明矩形剪切薄膜的屈曲过程是分枝点平衡问题。对于充气管的分析是基于张力场理论建立了其扭转临界扭矩的表达式,并与数值分析的结果进行了对比分析,分析表明充气管的屈曲过程是极值点平衡问题。通过对比分析,两种结构形式的分析结果与数值分析的结果基本符合。建立的分析方法能够有效地进行薄膜结构的屈曲及后屈曲过程的分析。     

大型双曲冷却塔表面脉动风压随机特性——风压极值探讨

针对风压信号呈现的非高斯特性和传统极限估计方法的局限,首次提出基于保证率和相关性的极值估计方法——全概率迭代法进行冷却塔表面脉动风荷载极值分析,并和传统的峰值因子法及改进的Sadek-Simiu法计算结果进行对比验证。结果表明：全概率迭代法避开了对随机过程的高斯分布假定,相比传统的极值估计方法其结果更加真实可靠;表达风压极值中脉动分量的峰值因子数值沿着环向和子午向变化显著,如取为同一数值则偏于危险或过于保守;采用全概率迭代法得到的表面风压系数极值分布曲线与规范取值相比,迎风面和负压峰值区域极值偏小,背风区域极值偏大,且最小值对应角度相差约10°。     

Newton型迭代法在求解烧蚀边界条件控制方程中的应用

对于采用烧蚀手段进行防热的高超声速飞行器而言,利用有限速率烧蚀条件下的防热层热响应与流场耦合计算,则可以更加准确地预测包括气体烧蚀产物在内的化学非平衡流场特性。基于稳态能量平衡假设(SSEB),本文对有限速率烧蚀边界数值处理方法进行了初步探索。采用Newton迭代法对表面质量和能量平衡方程进行迭代时,为了避免Jacobi矩阵求逆困难,发展了一种基于对角化Jacobi矩阵的Newton型迭代法。以碳基石墨防热材料有限速率烧蚀边界为例,通过与相关文献的壁面温度、热流以及流场不同组分质量分数等计算结果的比较,表明了这种Newton型迭代法在数值求解有限速率烧蚀边界方面的有效性。     

陶瓷基FGM材料线形变厚度圆板的热后屈曲

建立沿半径线形变厚度陶瓷基功能梯度材料(functionally graded material,FGM)圆板在热环境中后屈曲控制微分方程,采用打靶法研究陶瓷二氧化锆基变厚度圆板在温度场中的后屈曲行为,给出了均匀和非均匀升温场中变厚度圆板的热后屈曲平衡路径特性曲线,讨论陶瓷梯度指标和厚度变化系数对后屈曲行为的影响。数值结果表明:非均匀升温下FGM圆板的热临界载荷小于均匀升温下的临界载荷;径向厚度的变化并不会影响热临界载荷的值,但会影响到屈曲以后的平衡路径;中心挠度随变厚度系数的增加而增加。     

模态密集柔性空间结构二阶内平衡降阶

 给出一种针对二阶线性系统方程直接进行降阶的二阶系统内平衡降阶方法。大型柔性空间结构动力学方程采用二阶线性微分方程描述,采用传统的一阶内平衡降阶方法降阶后的状态方程是一阶形式,破坏了原系统的二阶结构和物理意义。采用新方法降阶后的系统可以保持原系统二阶结构,同时可以进一步保持原系统的对称和正定特性。柔性空间结构系统级降阶的柔性模态方程通常为对角形式,针对这种特殊形式,系统可控和可观Gramian矩阵存在闭合解析解,给出了闭合解的具体表达形式。数值仿真结果表明,二阶内平衡降阶方法可以达到一阶内平衡方法一样的降阶精度,Gramian矩阵的闭合解析解可以大幅度提高Lyapunov方程求解速度。     

压气机过渡段造型及三维数值模拟

提出一种基于S形曲线压气机过渡段造型方法.该方法将过渡段造型归结为S形内壁曲线拐点相对位置,面积分布率极值及其极值点相对位置3个几何控制因素.并采用此方法构造了一系列压气机过渡段,并针对这些过渡段进行三维数值模拟.结果表明:面积分布率极值是影响过渡段性能最重要的因素;可以通过调整面积分布率极值来控制过渡段最大面积处相对马赫数,减小外壁气流附面层厚度及支板形成的低压尾迹区;同时,配合变化较陡的内壁造型和合理的面积分布率曲线极值点相对位置,可以改善外壁形状,抑制附面层变厚.对于所研究的过渡段,内壁拐点相对位置为0.18,面积分布率极值点相对位置为0.20,相对马赫数为0.65时,总压损失最小.      

关于病态线性方程组的迭代解法

本文基于简单迭代法的基础上，提供解病态线性方程组的一种迭代方法，并给出一些数值例子及应用于高雷诺数下计算对流扩散方程。     

风扇/增压级动平衡工装最佳支承跨距计算方法

某型航空发动机风扇增压级平衡工装采用悬臂支承结构,其支承跨距是该平衡工装的重要设计参数。针对该平衡工装的结构特点建立了最佳支承跨距的理论分析模型,给出了最佳支承跨距计算的解析表达式,并通过数值计算对最佳支承跨距进行了验证。结果表明:平衡工装最佳支承跨距的解析解与数值解一致,建立的最佳支承跨距计算方法可以为风扇增压级平衡工装的支承跨距的参数设计提供理论指导,提高平衡工装的平衡精度和平衡效果。     

铝锂合金加筋壁板剪切屈曲性能

为研究剪切载荷下2A97铝锂合金加筋壁板的屈曲与后屈曲行为,设计了加筋壁板和夹具,完成了壁板的剪切试验;得到了加筋壁板的失稳载荷、破坏载荷以及破坏模式;采用受剪板屈曲与张力场理论计算了加筋壁板的剪切屈曲失稳载荷;建立有限元数值计算模型对加筋壁板屈曲行为进行计算分析,并将数值结果与试验结果对比。结果表明：加筋壁板的屈曲模式为筋条间蒙皮的局部屈曲;加筋壁板的破坏模式为沿加载对角线方向蒙皮的凸起,破坏原因为蒙皮的塑性变形、撕裂以及筋条的扭转变形;利用张力场理论可以得到较准确的屈曲失稳载荷,与试验误差为6.56%;数值模拟得到的屈曲与破坏模式与试验吻合,失稳载荷和极限载荷与试验结果误差分别为1.22%和11.52%。     

全变差正则化方法对一类热传导反问题的应用

热传导终端反问题是典型的不适定问题,应用基于全变差的正则化方法来获得它的稳定数值解,将近似解的求解范围由连续函数空间扩大到了有界变差函数空间,并将固定点迭代法用于求解非线性的Euler方程,并进行了理论分析和数值实验。结果表明:该方法与算法具有高效稳定全局收敛等优点。     

含分层损伤复合材料层合板前后屈曲行为研究

采用基于Mindlin一阶剪切理论的四节点板单元,分析了简支和固支两种边界条件下含圆形分层复合材料层合板前后屈曲行为。结果表明,尽管浅部分层的存在使层合板的前屈曲临界载荷大幅度下降,但层合板仍能继续承受很大的后屈曲载荷,这种情况下的前屈曲是由于分层部分的局部行为造成的。同时,在相同的铺设和分层厚度及边界条件下,分层半径越大,后屈曲临界载荷越低,固支下的后屈曲临界载荷比简支下的临界载荷高。     

复合材料层合加筋板结构的后屈曲强度及破坏研究

 对于复合材料加筋层合板结构提出了一个逐步破坏的分析模型，它包含结构总体变形的几何非线性分析和结构内部材料局部损伤破坏的分析及其两者的相互作用。基于这一模型采用非线性有限元分析方法，研究了在压剪载荷作用下层合板结构的屈曲、后屈曲路径和破坏时极限载荷。     

复合材料层合板的二次失稳和破坏研究

 本文采用考虑横向剪切效应的三角形拟协调罚函数单元,推导了带有初始挠度的复合材料层合板在后屈曲阶段的增量平衡方程和相应的切线刚度矩阵,以及非线性特征方程,并以数值计算研究了单元的精度和收敛性,二次分支点的存在性和边界条件、加载方式及铺层次序对二次失稳和破坏载荷的影响。     

转动壳体外压屈曲的有限元分析

由于转动离心力的影响,壳体的屈曲临界压力随转速的提高而升高。本文推导了转动壳体线性屈曲的有限元公式,并用DMAP语言编制相应的程序修改MSC/NASTRAN的解题步骤考虑离心力完成了分析。还利用MSC/NASTRAN的几何非线性分析功能完成了有初始几何缺陷转动壳体的有限元分析,计算后屈曲的载荷—变形曲线,得出不同转速下的缺陷敏感度曲线。     

均匀受压三铰圆拱的屈曲

 本文基于经典线性理论,同时考虑了屈曲的径向位移和周向位移分量,由Trefftz稳定准则导得屈曲平衡微分方程及中间铰点的联接条件,并进一步导得求解非对称,不等截面均匀受压三铰圆拱屈曲载荷的超越方程,从而扩大了解题范围。将对称等截面情况的计算结果与Dinnik的结果相比较表明:忽略屈曲的周向位移分量,导致偏高的计算结果。     

层板分层后屈曲分析的一维模型

 将条形分层的复合材料层板看成两个独立的子层,未发生分层区域的层间基体由拉-压杆和剪切板模拟,从而建立了分析带有分层缺陷层板稳定性问题的一维力学模型;利用Rayleigh-Ritz能量法推导了模型的微分方程;通过计算,给出了两个子层刚度比在各种数值下、对应于不同分层长度,层板的局部或整体失稳载荷。     

用圆弧样条计算地形跟随技术中的期望航迹

从地形外形的采样数据设计飞机贴地飞行的期望航迹,然后由相应的最优控制器使飞机按此航迹飞行,这一地形跟随技术中的最优控制算法,日益得到重视和发展,并应用在综合地形跟随/回避技术中。 本文提出了采用圆弧样条计算期望航迹的方法,在几十种地形情况下进行了对比研究,证明是可行的。明显地减少了期望航迹计算的工作量,增加了系统实时计算的可能性。     

用边界元法计算变截面变轴力柱的屈曲

 1.变截面变轴力柱屈曲计算的边界元法 变截面变轴力柱的挠曲方程的一般形式为式中 d²/dx²[EI(x)dw/dx²]+p(x)d²x/dx²=0,E——弹性模量;W——挠度;I（x）——截面处的横截面惯性矩;P（x）——x截面处的轴向压力;x——轴向坐标。 （1）式是非线性常微分方程,本文采用微分方程分解法,将（1）式分解成两个二阶变系数常微分方程。     

推力矢量飞机平衡区的计算

通过对推力矢量控制下飞机动力学特性的分析，定义了飞机机动的平衡区并进行计算，确定了推力矢量在低速机动中对飞机平衡特性的决定性影响，并与常规飞机的平衡特性进行了比较分析。将以往用来分析飞机尾旋运动的分歧与突变理论（ＢＡＴＣＭ）推广到推力矢量飞机过失速平衡特性的计算中，确定了在不同的飞行状态及推力矢量系统配置下飞机过失速机动的平衡区。其结果有助于理解推力矢量系统的效用及其设计参数对飞机过失速机动能力的     

受压固支正交异性矩形分层的稳定性分析

在足够大的压力作用下,层合板中的分层（delamination）会产生屈曲和后屈曲。本文分别用力法和位移法求解了受压固支正交异性矩形分层的屈曲和后屈曲问题,把一组非线性偏微分方程构成的边值问题化为求解一组非线性代数方程组,最后给出了有关的数值结果。