用混合标架法处理药柱尾部人工脱粘滑动边界条件

介绍了采用混合标架法处理药柱尾部人工脱粘滑动边界条件以及对总刚方程进行对称处理求解的方法，用此法对某药柱模型进行了轴向加速度载荷作用下的应力分析，并与有关文献作了比较的分析，证明计算结果更接近实际。     

复约束方法计算斜齿圆柱齿轮结构特征值问题

根据齿轮结构的循环对称性，以仅含有一个完整轮齿的扇形子结构作为齿轮有限元计算模型，用计算循环对称结构特征值问题的复约束方法，通过求解一个扇形子结构的特征值问题得到整个齿轮结构的特征值和特征向量。计算了４个不同参数齿轮的低阶固有频率和振型，并与实验结构数据比较，获得了满意的结果。这证明复约束法既极大地减少了计算所需内存和机时，又可建立包含轮齿和轮体在内的计算模型，保证了齿轮结构特征值问题计算的完整性和准确性。     

改进节点积分的无单元Galerkin法及其在流动问题中的应用

无单元Galerkin法需要在背景网格上积分,计算量大,且在求解对流占优问题时会出现非物理的数值伪振荡现象。为此,基于局部Taylor展开思想,采用节点处的局部Taylor展开计算积分,建立了局部Taylor展开积分无单元Galerkin法。该方法同时解决了标准的无单元Galerkin法计算量大和对流占优时会出现数值伪振荡的问题。一维定常对流扩散方程和二维Burgers方程的求解说明了该方法的有效性,且计算效率远高于无单元Galerkin法。     

循环对称结构特征值问题的复数域子空间迭代法

提出求解在循环对称结构固有特性问题的研究中出现的复Ｈｅｒｍｉｔｅ阵特征值问题的复数域子空间迭代法。根据Ｔｈｏｍａｓ等人提出的理论，生成循环对称结构任一子结构的复约束质量阵和复约束刚度阵。根据本文提出的办法，生成一个包含ｑ个相互质量正交归一的复向量组成的初始里兹迭代复向量组，并将原广义特征值问题的求解空间投影到由这ｑ个复向量张成的Ｋｒｙｌｏｖ空间上可很方便地在ｑ维Ｋｒｙｌｏｖ空间上进行广义复Ｈｅｒｍｉｔｅ矩阵特征值问题的求解。     

复合材料翼面结构气动弹性剪裁

 <正> 1.引言 本文采用了混合设计变量。对复材蒙皮,由于其连续性,我们用双三次多项式来描述各定向铺层的厚度,将多项式系数做为设计变量,对于其他部件如肋、腹板等,我们仍用其有限单元的截面积或厚度做为设计变量。整个结构采用有限元法做结构分析。固有振型采用高效率的谱变换Lanczos法计算,非定常气动力用亚音速偶极子格网法计算,用可行方向法求解优化问题。     

曲边单元间断有限元方法求解二维Euler方程

考虑到间断有限元方法对边界的敏感性,采用基于八节点曲边四边形单元的间断有限元方法求解了Euler方程的圆柱绕流问题。详细推导了八节点四边形单元的变换关系,给出了Jacobi矩阵行列式的具体表达式。对比直边单元和曲边单元的计算结果,采用曲边单元后,计算结果符合Euler方程的无粘假设,得出了八节点四边形单元间断有限元方法求解Euler方程是合适的结论。     

间断探测器在间断Galerkin方法中的应用

在二维非结构网格上,对高阶间断Galerkin方法求解定常二维欧拉方程进行数值研究。为了很好地抑制解在流场间断附近处产生的伪振荡,采用了间断探测器和限制器相结合的方法,并分别对Krivodonova间断探测器和基于激波物理特征的激波探测器进行了比较和研究。数值结果表明：如果流场中只存在激波时,两者探测效果相似,且后者更适合应用于求解激波问题。如果流场更加复杂,即存在多种间断时,前者依然可以比较准确地用来探测间断单元,能够很好抑制流场间断附近处产生的伪振荡,而后者失效。     

基于EBE策略求解流场问题的有限元方程解法研究

本文首先阐述了用有限元法求解流场问题的基于单元接单元（EBE）技术的并行双共轭梯度法（Bi-CGSTAB)。然后引入网格分维概念，并基于网格分维分析和比较了直接法等三种算法的存储量和操作数。最后在国家性能计算中心（成都）的曙光2000并行系统上，采用MPI消息传递实现了此种算法，计算了定常三维空穴流，取了较好的加速比。得到如下结论：EBE迭代法不必生成和存储总刚矩阵，在很大程度上节省了内存，并且它具有良好的并行性，所以EBE的迭代法能有效提高求解流场问题的规模和速度。     

含支承松动故障的弹用涡扇发动机整机振动建模与机匣响应特征

针对某型弹用涡扇发动机结构特点,建立了一种转子-支承-机匣整机模型,对转子与机匣采用有限元梁模型,支承采用集总质量模型,引入支承松动故障模型,利用数值积分方法求解耦合系统的响应.基于机匣加速度信号,研究了对称刚度以及不对称刚度模型下松动故障的冲击特征分析.结果表明:①松动故障所引发的机匣加速度时域波形具有上下不对称性冲击特征以及频谱中出现倍频特征;②仿真计算结果与实际弹用涡扇发动机试车数据时域波形特征以及频谱特征非常一致,而且验证了不对称刚度松动故障模型更适合弹用涡扇发动机的松动故障建模.      

缘板及叶冠间的接触状况及其对叶片位移和应力的影响

1.计算方法及其实验验证 本文用分析法见文献[1],这里以WS9高压二级涡轮转子叶片(简称WS9叶片)为例来说明方法的特点。该转子轮盘上均匀安装110片叶片,组成具有循环对称接触结构的叶片环,有接触表面的循环对称结构。据此,可将图1所示的一片叶片取作基本段分析,但同时必须在该叶片一般有限元刚度方程中引入缘板及叶冠的循环对称接触表面(叶冠的A1侧面与A2侧面,以及缘板的B1侧面与B2侧面)处的循环对称接触     

周期对称性在叶盘结构瞬态响应求解中的应用

研究了利用结构周期对称性计算叶盘结构瞬态响应的方法.首先用ANSYS得到扩充扇区的刚度和质量矩阵,通过周期对称边界上节点的对应关系得到扩充扇区到基本扇区的转换矩阵.利用叶盘结构中各扇区所受载荷的特点,得到任意时刻载荷的空间分布,并对其进行Fourier展开,以得到周期对称模型的载荷.采用Newmark法进行时域积分,利...     

结构分析程序中几个问题的处理方法

JY851程序是一个复合材料、金属混合结构应力分析及强度校核程序。从1986年投入使用以来,已推广到几个厂、所使用,计算了几个型号的机翼、垂尾、方向舵等复合材料结构。现就编写元素刚度矩阵中所遇到的几个问题的处理方法作一简单的介绍。 1.复合材料板弹性模量阵的计算 复合材料板弹性模量阵A、B、D的计算中,关键是求偏轴模量阵[ō],层压板中第     

轮盘-叶片组件的统一应力分析

本文用含接触面的循环对称计算模型进行轮盘-叶片组件的统一有限无应力分析。此法考虑了榫头与榫槽间的三维接触传力,利用循环对称性,不人为假定轮盘切开处的边界条件,可同时得到轮盘及叶片中的应力分布。还用另外五种简化计算模型来分析同一计算对象。将所得结果同含接触面的循环对称模型的计算结果作比较,评价了简化模型的适用性。     

转子分区循环对称接触应力分析

 通过光弹实验表明了端齿连接转子受力的分区循环对称性。针对这~特性建立了各循环对称分区的循环对称接触有限元应力分析方法,并利用叠加原理将各循环对称分区间过渡区化为几个小规模的循环对称问题。     

基于子结构法的大型结构数值敏度计算技术

针对大型结构直接精确敏度分析方法解方程规模大、求解时间长的问题,发展一套基于有限元子结构法的处理技术。大型结构中,数值敏度计算的伪载荷仅与涉及设计变量的结构单元有关,因此存在大量零元的可能性。根据此特点通过节点重排将与设计变量有关的节点位移排到总位移列阵序的后面,按重排后的顺序投放刚度矩阵,然后对其进行区域分块,并聚缩得到规模较小的矩阵。算例计算表明:利用此矩阵求解结构位移场的导数,不但保持了精确法的精度,而且在优化设计的敏度计算中,这些较小规模的数据又被多次使用,从而显著提高计算效率。     

旋翼桨叶动力分析及内力计算中有限元法的改进

提高用有限元法进行桨叶动力分析、特别是动内力计算的精度是本文工作的主要目的。本文建立了变剖面旋转梁协调单元族;提出利用单元动刚度矩阵直接由节点位移计算桨叶内力（简称为计算内力的“动刚度法”）。作为例子,用所导出的不同精度的有限单元和建议的内力计算方法计算了具有不连续结构特性的变剖面桨叶弯曲振动固有频率、振型和模态内力。计算结果表明,工作可以达到预期的目的。     

三维正交机织复合材料单胞特征单元及其应用

 基于三维正交机织复合材料,提出了一种特征单元均匀化分析方法。用“特征单元”来表示能够反映单胞的几何构造细节和材料构造细节的单元。首先用一般有限元方法对三维正交机织复合材料单胞进行分析得到总体刚度矩阵,利用该矩阵的特征向量建立了单胞的特征单元。通过比较特征单元的刚度矩阵元素和利用传统均匀化方法得到的刚度矩阵元素,说明了特征单元能够反应单胞的几何构造和材料构造细节。为了验证特征单元的有效性,分别用特征单元、传统均匀化方法和一般有限元方法计算了三维正交机织复合材料梁的固有频率,结果表明特征单元精度比传统均匀化方法高,而计算量比一般有限元方法大幅度降低。     

EBE技术在结构分析中的应用——EBE—CG方法及其并行实现

共轭梯度法(CG法)是求解大型稀疏有限元方程组的有效迭代解法。针对大规模结构分析问题,提出了不形成总刚度矩阵时的EBE-CG法及其并行处理方法。在国内规模最大的大型MIMD分布式系统上的初步数值试验结果表明,并行的EBE-CG法是很有效的。     

用试验数据修正刚度矩阵

讨论用试验数据修正刚度矩阵的问题。依据特征方程、刚度矩阵的对称性和半正定性,利用代数特征值反问题的理论和方法,研究了这个问题解的存在性和唯一性,提出了修正刚度矩阵的一个新方法。用这个方法修正的刚度矩阵不仅满足特征方程,而且是唯一的对称半正定矩阵。