一类非线性伪抛物型方程的伪谱解法

本文考虑了一类非线性伪抛物型方程的Fourier伪谱方法,建立了该方程的Fourier伪谱方法的半离散格式和全离散格式．并利用Sobolev空间的正交映射理论,给出了这两种格式的误差估计．最后针对全离散格式给出了数值算例,数值结果表明Fourier伪谱格式能正确加解密,且计算误差较小,效率较高,具有较好的稳定性,可用于提高热流密码体制的加解密效率．     

变系数情形下Criss-Cross三角形线性元的渐近展式与超收敛

本文讨论了二阶椭圆方程变系数情形下Criss-Cross三角形线性元的超收敛性质,得到了有限元的渐进展式、外推及高精度组合公式等结果.     

功能梯度材料板件三维分析的半解析梯度有限元法

将半解析有限元与梯度有限元相结合,形成一种半解析梯度有限元来求解功能梯度材料板件问题。该方法兼有有限元法的适应性强、程序统一,半解析有限元法的节省单元与计算工作量,梯度有限元法的适应构件内部材料性能任意梯度分布等特点,并实现用一维数值计算给出构件三维分析结果。算例分析表明了方法的精度、功能与上述特点,充分揭示了功能梯度材料板件力学响应的三维形态。半解析梯度有限元法可推广应用到其他功能梯度材料面结构的各类分析中。     

基于高阶位移模式的压电层合板动力有限元模型

建立了含压电片层合板的有限元动力学模型。以位于压电层上下表面处的电场强度和层间电压为未知量,给出了三次函数的电势分布模式,采用Reddy的高阶剪切理论描述板的位移场,假设板厚度方向的正应力为零给出了减缩的本构方程,采用有限元方法,基于Hamilton原理导出结构的动力学方程,然后用静态缩聚的方法压缩掉电场自由度和次要的位移自由度。最后用四边形矩形单元求解了一对称铺层和非对称铺层悬臂板的固有频率,并与ANSYS结果对比验证了本文模型的精确性。     

时-空混和元与多辛

指出所谓多辛相当于时间-空间混和坐标内的分别离散的规则网格;而时-空混和元则不用对时间、空间分别离散,而是组成混和的时空混和有限元网格,更加灵活些。     

求解带自由面渗流问题移动网格法的改进及其工程应用

采用有限元固定网格和移动网格相结合的方法,求解带自由面的坝体稳定渗流问题。本方法是针对一般坝体的结构,在有限元网格生成之后,又自动生成了网格移动所需的信息,实现了网格在迭代计算中随自由面的变化自动均匀移动;讨论了网格奇异的不同类型,提出了避免单元形状怪异的移动网格方法;在欧拉空间中分区记录坝体不同的渗透性能,在网格移动到新位置后,自动识别网格所在的物性空间,从而确定单元在新位置上的渗流系数。本文将该方法应用到吉林丰满混凝土重力坝防渗加固的实际工程当中,在考虑大坝原有防渗措施的前提下,利用数值分析,研究了防渗芯墙对大坝渗流性能的影响。计算结果表明,加芯墙后坝体自由面位置明显下移,湿水部分大为减少,坝基扬压力没有明显提升,下游出水点位置没有明显下降。计算过程表明,本文方法及措施具有很好的收敛性与稳定性,且收敛性与自由面初始位置的选择无关。     

高精度非线性格式WCNS的分析研究与其应用

首先将Fourier方法推广于高维方向研究了五阶精度WCNS的特性,并与其他高阶格式进行比较。分析结果表明WCNS的高精度特性普遍接近甚至好于迎风偏置五阶显式格式EUW-5与Pade′标准格式。然后开展了WCNS的应用研究,采用高效率的WCNS-E-5数值模拟了含强激波的高维复杂流场。算例包括二维高超声速边界层对自由流扰动波的吸收问题以及三维球头绕流问题。计算结果反映出WCNS-E-5对激波等间断的良好捕捉能力,图像清晰光滑,数据准确可靠。     

磁流变弹性体若干物理量的数值分析

应用有限元方法,考虑了颗粒的磁化饱和过程与非线性磁化过程,计算得到了磁流变弹性体中的磁场分布,进而研究了在不同磁场大小、不同颗粒体积比浓度下磁流变弹性体在成链方向的相对磁导率,计算结果和实验结果取得了一致。利用Maxwell应力张量,计算了磁场引起的磁流变弹性体的附加剪切模量。分析了颗粒体积比浓度、外加磁场对磁流变弹性体的磁致剪切模量的影响。研究了颗粒为旋转椭球形状时,颗粒的放置方式与其长短轴之比对磁流变弹性体的磁致剪切模量的影响。计算结果表明,磁流变弹性体的相对磁导率随颗粒体积比浓度的增大而增大,随磁场强度的增大而减小,颗粒的形状和放置方式对磁流变效应有很大的影响。     

桁架结构稳定分析的几何非线性欧拉稳定理论

目前,桁架结构稳定分析有两个理论:经典的特征值理论和几何非线性临界点理论。前者发展较早,在许多结构力学书中都有叙述。但是经过工程应用检验,它过高地估计了结构抗稳定的能力。后者是20年前提出的,提出者认为适用于所有的扁桁架。经过几年的研究,作者连续提出三个稳定理论:线性与非线性欧拉理论;非线性临界点-欧拉理论和两个计算方法,进行桁架结构的线性和非线性截面优化设计。经过理论研究和对前人文献中数值例题的计算和比较,发现特征值理论不符合实际,临界点理论只适用于大扁度的桁架,而不适用于一切扁桁架。研究结果提供了若干有用的结论,给出了各种稳定理论的适用范围,指出了国际文献中若干例题的正确和错误,补充了桁架结构的稳定分析理论,使之更为完整和实用了。     

离散元与壳体有限元结合的多尺度方法及其应用

在深入研究复杂结构和非均质材料冲击响应和破坏机理的过程中,往往遇到多尺度计算问题。本文尝试建立三维离散元与壳体有限元结合的多尺度方法用于处理圆柱壳问题,该方法采用三维离散元对感兴趣的局域进行局部模拟,利用平板壳体有限元进行整体模拟,采用一种特殊的过渡层使离散元区和有限元区能很好的衔接。我们将这一方法应用于激光辐照下充压柱壳的热/力耦合冲击破坏响应,得到的模拟结果与文献报道有较好的吻合。