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

提高用有限无法进行桨叶动力分析、特别是内力计算的精度是本文工作的主要目的。本文建立了变剖面旋转梁协调单元族,给出了位移函数为三次、五次、七次幂多项式时变剖面旋转梁协调元的几何刚度矩阵显表达式。文中还提出了直接用节点位移计算内力的“动刚度法”。作为例子,用所给出的有限元和方法计算了桨叶弯曲振动固有频率、振型及模态弯矩。计算结果表明,工作可以达到预期的目的。     

短粗梁的固有振动计算

用有限单元法计算,选取总线位移、弯曲转角和剪切转角为节点自由度,舍弃常剪力假设,总线位移取三次函数,剪切转角为线性函数,满足几何边界条件。由简支梁、悬臂梁和自由梁的计算结果表明,精度是比较高的。     

分布传递函数方法的梁杆几何非线性分析

本文提出了一种基于动坐标的、适用于梁大变形分析的传递函数方法，在本方法中，一个复杂系统被分为若干一维子系统，各子系统根据梁的有限变形理论进行推导，推导中把轴力看作子系统的参数，引入状态向量后，可以将子系统的控制方程改写为含参状态空间的形式，然后利用传递函数方法求解，系统的合成借鉴了有限元的处理方法，引入了动坐标关系，使方法可以用于对大变形的分析。文后给出了算例表明本文方法有很高的精度和很好的收敛性。     

弹性转动约束复合材料层板的中等大挠度

本文首先用虚位移原理推导出以位移形式表达的Reddy型高阶剪切变形理论复合材料层板的非线性控制方程及相应的边界条件。选定的5个位移函数均满足弹性转动约束边界条件，用Galerkin方法把无量纲化之后的控制方程转化为一组非线性代数方程组，用线性化的方法和可调节参数的修正迭代法求解这组方程。最后求出了不同复合材料的挠度和弯矩值。     

塑性动态断裂过程中梁的塑性铰转动因子的确定方法

运用理论分析方法和实测试验方法，探讨了梁在塑性动态断裂过程中，其塑性铰转动因子的变化规律，提出了断裂过程中塑性铰几何转动因子的概念及其确定方法，并以３点弯曲梁理论分析解为例，计算了梁在塑性动态断裂过程中，其塑性铰转动因子的变化规律。该计算结果与３点弯曲梁动态断理解实验实测相一致。理论分析进一步明确了梁的塑性铰转动因子具有几何与物理两种不同的本质含意，并可直接用于修正梁的塑性动态断裂的理论分析解。     

转动旋翼模型动柔度的试验研究

介绍测定转动旋翼模型动柔度的试验方案、试验程序，综述以作者自行设计的旋翼模型进行试验所得到的结果及其与理论计算结果的比较。结果表明，试验方案与方法是可行的，理论计算公式是正确的。     

二维几何非线性问题的线有限元解法

提出了在线有限元方法中考虑几何非线性二维问题的理论公式和计算步骤，并以完全的拉格朗日表述为基础，首次导出了用线有限元方法进行结构分析中大变形问题载荷增量法的基本公式，对计算过程中的重要方面：包括对非线性微分方程线性化的方法，分步迭代的步骤等进行了详细讨论，以这些公式为基础，编制了相应的计算程序，并给出了算例，结果是令人满意的。     

用连续质量有限元素法求解动力系统的固有特性

本文简略地回顾了假定位移函数的有限元素法求解动力系统固有特性存在的问题,提出了连续质量有限元素法,并将其推广到适用于链式结构的传递矩阵法。具体推导了均匀的和线性变截面的一维元素(包括杆、轴和梁元素)的动刚度矩阵和传递矩阵。最后用其对各种梁的固有特性作了一些对比算例。     

纤维增强塑料拉挤型材弯曲强度研究

研究了纤维增强塑料 ( FRP)拉挤型材的弯曲强度。试验结果表明 ,FRP拉挤型材在大跨距 (跨高比 l/h>1 0 )下的弯曲破坏模式通常是梁的纵向开裂 ;FRP拉挤型材由于剪切变形较大 ,强度计算时不能按普通各向同性梁理论进行计算 ,应考虑剪切变形的影响。为此 ,本文应用复合材料薄壁梁理论给出了 FRP拉挤型材在大跨距下弯曲强度计算的半经验公式     

大机动无人机飞行轨道计算中两个问题的分析

本文讨论无人机作大机动飞行时,轨道计算中必须考虑的两个问题。 一个问题是按照大机动无人机的油箱结构,考虑燃料消耗和前后左右的波动所引起的无人机重心位置的变化。这个变化引起的作用在无人机上的力矩的变化对无人机进入和退出大机动时的飞行参数有明显的影响。另一个问题是da/dt的计算,文中推导了准确计算的公式,按此公式考虑无人机转动受下洗气流的影响,可以准确地计算出阻尼力矩的影响。 在轨道计算中计及这两个问题后,不仅可以计算出无人机在大机动飞行中高度和速度的波动,而且可以计算出滚转角的波动,这和实际飞行情况是比较接近的。     

基于流体速度势的非线性晃动问题的ALE有限元分析

基于流体速度势，本导出了非线性晃动问题的ALE有限元格式。通过引入ALE运动学描述，使网格的运动独立于流体的运动，通过采用与贮箱固连的非惯性参考系建立控制方程，使网格只须相对贮箱更新，因此运动边界和自由液面的跟踪算法不仅非常简单，而且不会引起网格畸变．本分别采用有限元法和有限差分法进行空间和时间离散，中的算例分别模拟了贮箱受横向激励和俯仰激励时流体的非线性晃动问题，所得结果证实了本方法的有效性与可靠性。中的方法还可进一步推广用于流体—结构耦合问题．     

变截面梁的弯矩和位移计算法

变截面梁的控制方程为变线数微分方程 ,一般只能用各种近似方法求解 ,而且比较复杂。对于梁式构件 ,本文不求解变系数方程 ,而是用力法 ,在力法方程中 ,主、副系数和自由项的莫尔积分使用积分表 ,可以方便而且精确地求出变截面梁的弯矩和位移     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

多视图场景确定下的空间点位置视觉标定方法

针对空间特征点三维位置标定问题,提出一种多视图场景确定下的空间点位置视觉标定方法。该方法基于计算机视觉原理,通过在场景内放置棋盘模板确定图像场景位姿参数,而后根据多视图图像下的空间点图像坐标和场景位姿参数求解空间点三维位置坐标。试验结果表明,该方法下的空间点标记位置与真实位置间的距离误差小于1 mm,适用于高精度要求下的特征点位置标定。     

设计灵敏度分析的迭代模态法

在结构动态有限元数学模型修正、结构动力学修改以及结构动态设计中，灵敏度分析是十分重要的一环，往往也是主要计算工作量所在。如何在保证精确度前提下减少灵敏度分析工作量和计算时间，具有重要意义．本文提出的迭代模态法，可以大大提高特征向量对设计参数导数（模态灵敏度）的计算效率。本文还从理论上证明经典模态法和修改模态法只是迭代模态法的特殊情况。文中还给出了ＦＯＲＫ，３Ｄ－ＦＲＡＭＥ两种结构为实例的计算机仿真，并与已有两种模态方法进行比较。理论分析与计算机仿真表明，本方法不仅能充分保证特征向量灵敏度分析的精度，而且大大提高了计算速度．     

三维位移场测量的贴片云纹干涉法

三维位移场测量的贴片云纹干涉法是用多束准直相干光同时入射到贴有全息胶片的试件表面,并记录由其相互干涉形成的三组频率不同、方向不同的全息光栅。若对试件加载前后两次曝光,可得到含有试件三维位移场(u,v,w)的三组云纹图。利用光栅衍射方向的差异,容易在三个彼此分开的方向上观察到三组云纹,由此可定量分析试件的三维位移,无需进行复杂的光学滤波处理。文中给出了严格的理论推导和实验验证。     

翼型测力天平的研制

为研究翼型结冰状态下的载荷变化规律,对翼型进行测力实验.研制一台三分量翼型测力天平.天平设计的突出难点为设计载荷极不匹配,阻力测量载荷很小.设计时,针对翼型载荷的特点,采用片梁组合方式测量阻力,四柱梁测量法向力及俯仰力矩,并采用有限元方法对结构进行了验证计算.主要介绍天平的设计、校准及应用情况.     

用于非定常粘性流动计算的无网格算法(英文)

把无网格算法发展用于求解涉及动边界的非定常粘性流动问题。在处理粘性流动的布点问题时,通过在远离物面的区域采用无粘流动计算时所用的各向同性点云,而在物面附近引入各向异性点云,从而在保证物面法线方向布点较密以准确模拟边界层的同时,有效控制了布点总量,减少了计算时间。对于动边界问题,本文在上述方法获得的初始布点的基础上,采用基于扰动衰减规律的点云移动技术,快速得到物面运动到其他位置时流场求解所需点云。在所获得的点云上,采用一种带权系数的二次极小曲面逼近方法来离散Navier-Stokes方程的空间导数,权系数的引入使得流场求解时对各向同性点云和各向异性点云无需分开考虑,而是采用统一的求解方法,从而简化了编程求解。用发展的无网格算法,结合Navier-Stokes方程求解双时间推进方法,并耦合Spalart-Allmaras湍流模型,先成功地模拟出NLR7301翼型后缘摆动诱发的非定常流,并通过与实验比较,验证了本方法.接着进行了NACA0012失速模拟,给出了动态失速过程,展示出用本文方法处理复杂动边界问题的效果。