建立了预估谐振载荷作用下结构振动疲劳寿命的损伤力学-有限元方法。首先根据损伤热力学原理,构建了损伤演化方程,建立光滑试件在恒幅应变交变载荷作用下寿命预估方法;进一步由损伤力学守恒积分原理,得到恒幅重复载荷作用下应力与寿命的关系式;然后根据标准件疲劳试验结果,拟合得到损伤演化方程中的材质参数;最后利用APDL语言编程对ANSYS软件进行了二次开发,借助ANSYS软件对谐振载荷作用下结构振动疲劳裂纹萌生寿命进行预估。作为算例,本文利用该方法预估了LC9CgS1铝合金梁谐振载荷作用下疲劳裂纹萌生寿命。

本文运用有限体积方法，对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响，选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离，Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题，流动特征在很大程度上取决于间距的大小。且间距存在一临界值，间距比从小于临界值变化到大于临界

对于具有局部非线性的多自由度动力系统，提出一种有效方法，该方法将线性自由度转换到模态空间中，并对其进行缩减，而将非一自由度仍保留在物理空间中，避免了在数值分析中求非线性因素时的坐标转换。

用各向异性模型定义NURBS曲线的变形能,基于所提出的变形能模型,用有限元法对NURBS曲线表达的设计边界进行模态计算,然后,将设计边界用模态向量的线性组合参数化表示.将这种基于边界特征向量的几何形状表示方法应用于优化参数定义,提出了一种适用于结构形状优化的自适应几何精化方法,它有效地将p型有限元分析、优化方法和设计边界的形状表达集成在一起.

研究在平面声波斜入射情况下，无穷大双周期加筋的微穿孔薄板结构的振动响应及吸声性能。首先在马大猷和Takahashi微穿孔板声学理论基础上，建立了微穿孔加筋薄板结构振动的半解析模型；然后应用傅里叶变换及空间波数分析方法，将周期加筋微穿孔薄板的振动位移及辐射声压表示为频域内波数的分量迭加形式；最后通过对波数分量进行求解，并利用傅里叶逆变换得到双周期加筋的微穿孔薄板的振动响应及吸声系数表达式。计算结果表明，薄板的弯曲振动在水中对吸声的影响较大，空气中仅对轻质穿孔板的低频吸声效果有一定影响；同时微穿孔率对周期加筋薄板吸声系数的影响明显，通过改变穿孔率和加筋周期等可有效地提高水中微穿孔薄板结构的吸声性能。

由双控制体模型的三个控制方程入手，将其无量纲化。并分解为零阶和一阶无量纲方程。由零阶方程求得定常流场。继而求解一阶方程得到密封的动特性参数，其中一阶方程以简洁的矩阵形式表达，操作更为简便。本文双控制体模型是较为成熟且得到认可的分析方法，双控制体改进模型具有更为清晰的模型及更为简洁的算法。最后与传统模型计算结果比较．验证了本文算法的准确性。随后分析了几种密封参数对其动力特性的影响．揭示某些规律并给出一些建议。

提出一种改进的栅格法网格剖分算法，该算法综合了栅格法和Delaunay算法的优点，生成了四边形和三角形混合单元，保证了整体性能与局部性态最优，满足汽车覆盖件冲压成形模拟软件对有限元网格质量的需求。实践证明，本文提出的改进的栅格法网格剖分算法可行有效，生成网格的质量高。程序容易实现。该算法已经成功集成到商业化冲压成形模拟软件KMAS中。

基于数值方法实现快速有效地分析评价大跨桥梁颤振稳定性。针对强迫振动法识别颤振导数试验中的大跨桥梁二维节段模型，利用RNGk-ε湍流模型并采用有限体积法求解桥梁模型绕流二维不可压缩流体Navier-Stokes方程，通过计算桥梁断面模型在周期运动中少数离散时刻的气动力，利用最小二乘法计算颤振导数，采用SCANLAN方法求解颤振临界风速，最终进行颤振稳定性分析。通过该方法计算出了丹麦Great Belt East桥和我国虎门大桥的颤振临界风速。计算结果与已有试验结果十分接近，进而验证了该数值方法的有效性和可靠性。

基于Reddy提出的板高阶剪切变形简化理论，研究了含界面脱粘损伤压电复合材料层合板非线性动力稳定性问题。首先，建立了分层模型。推导了考虑几何非线性、阻尼效应、纵向惯性力和力一电耦合效应的Mathieu方程，并且给出了该方程解的解析表达式。其次，通过典型算例讨论了界面脱粘损伤以及反馈控制力对该层合板动力不稳定区域、纵向、横向共振频率和最大“牵引”深度的影响。由典型算例讨论可知：随着层合板界面脱粘损伤的扩大,其动力稳定性能逐渐减弱，其中在损伤较小时，反馈控制力对智能结构几乎没有影响；而在损伤比较大的情况下，反馈控制力将能有效地减少动力不稳定区域重合面积。

用张量运算推导了弹性应变梯度轴对称问题的基本公式。建立了应变梯度轴对称不协调元的弱连续条件，进一步建立了满足弱连续条件的应变梯度轴对称18-DOF三角形单元（BCIZ＋ART9）。其中BCIZ满足线性应变C^0连续，用于计算应变ε；ART9满足常曲率C^1弱连续，用于计算应变梯度刁。数值结果表明该单元通过C^0-1分片检验并能体现材料的尺度效应。