基于无网格Galerkin法的穿甲侵彻数值模拟

将无网格Galerkin法与LS-DYNA软件相结合,建立了穿甲侵彻的分析模型,对三维球形弹体侵彻双层无限大钢靶板进行了数值模拟,解决了该过程中出现的负体积现象,并将结果与空穴膨胀理论和有限元法进行了对比,所得结果证实了无网格Galerkin法在模拟穿甲侵彻时的可行性与优越性.     

三维无网格Galerkin结构分析方法及其应用研究

针对无网格Galerkin法在三维复杂几何形状的结构分析中存在的刚度矩阵的稀疏存储实现难、六面体背景网格适应性差等问题,本文采用逐节点对法组装刚度矩阵,利用CSR格式存储刚度矩阵,提出了一种基于四面体背景积分的改进的三维无网格Galerkin法。通过采用罚函数法施加对称约束和给定位移值约束,并推导出了施加这两种位移约束的统一格式。利用所提算法完成了三维悬臂梁的计算,所得结果与其理论解相吻合;完成了轴流式风机轮毂的结构分析,得到的位移与应力分布结果与其有限元解相吻合。这表明本文所提方法能满足工程应用中的计算要求,并适用于具有复杂形状的几何模型分析。     

正交各向异性材料热弹性问题的三维无网格法计算模型及应用

建立了正交各向异性材料热弹性问题的三维无网格伽辽金(Element Free Galerkin, EFG)法计算模型。利用该计算模型对三维复合材料汽轮机叶轮和轴承座进行了热弹性分析,对比了材料方向角及热导率因子、热膨胀系数因子和拉压弹性模量因子不同组合情况下轴承座的最大热变形总位移和当量应力值,讨论了材料方向角及上述正交各向异性因子对热变形和当量应力的影响规律,并与各向同性材料进行了对比。结果表明:三维EFG模型的热变形总位移和当量应力相对误差范数分别比有限元法小0.1215%和0.1359%;材料方向角同时影响热变形的大小和方向,但对当量应力方向影响不大;正交各向异性材料因子主要影响热变形和当量应力的大小。在考虑热-机械载荷作用下的三维复合材料零件结构设计中,当以刚度或强度为主要需求时,材料方向角、热导率因子、热膨胀系数因子、拉压弹性模量因子分别在(45°~60°,8:1:4~10:1:5,(1/6):(1/5):1~(1/5):(1/4):1,(7/5):1:(9/5)~(3/2):1:2)或(0°~10°,(1/10):1:(1/5)~(1/8):1:(1/4),(1/5):1:(1/6)~(1/4):1:(1/5),1:(1/5):(1/10)~1:(1/4):(1/8))范围内取值能有效降低轴承座等结构的热变形和当量应力。     

耦合GPU与PCG的EFG法并行计算及应用研究

针对迭代法求解无网格Galerkin法中线性方程组收敛速度慢的问题,提出了一种耦合GPU和预处理共轭梯度法的无网格Galerkin法并行算法,在对其总体刚度矩阵、总体惩罚刚度矩阵进行并行联合组装的同时即可得到对角预处理共轭矩阵,有效地节省了GPU的存储空间和计算时间;通过采用四面体积分背景网格,提高了所提算法对三维复杂几何形状问题的适应性。通过2个三维算例验证了所提算法的可行性,且预处理共轭梯度法与共轭梯度法相比,其迭代次数最大可减少1686倍,最大的迭代时间可节省1003倍;同时探讨了加速比与线程数和节点个数之间的关系,当线程数为64时其加速比可达到最大,且预处理共轭梯度法的加速比与共轭梯度法相比可增大4.5倍,预处理共轭梯度法的加速比最大达到了88.5倍。     

积分背景网格对自适应EFG法拓扑优化的影响研究

采用应力能量范数作为误差指标,探讨了EFG法中积分背景网格对计算精度的影响,得到了合理划分背景网格的建议;建立了以节点密度为设计变量、以最小化柔度为优化日标的拓扑优化模型。采用以节点密度值为加点判据的自适应规则加点方案,开展了连续体结构的拓扑优化研究,该加点方案能有效地减少设计变量的个数,探讨了背景网格对拓扑优化结果的影响。算例结果表明,采用合适的背景网格不仅能进一步减少设计变量的个数,而且能够改善拓扑优化结果的光滑性,使计算效率和精度得到提高。     

无网格Galerkin法GPU加速并行计算及其应用

针对无网格Galerkin法计算耗时的问题,采用逐节点对法来组装刚度矩阵、共轭梯度法求解基于CSR格式存储的稀疏线性方程组,提出了一种利用罚函数法施加本质边界条件的EFG法GPU加速并行算法,给出了刚度矩阵和惩罚刚度矩阵的统一格式,以及GPU加速并行算法的流程图。编写了基于CUDA构架平台的GPU程序,且在NVIDIA GeForce GTX 660显卡上通过数值算例对所提算法进行了性能测试与分析比较,探讨了影响加速比的因素。算例结果验证了所提算法的可行性,并在满足计算精度的前提下,其加速比最大可达17倍;同时线性方程组的求解对加速比起决定性影响。     

基于虚荷载变量的形状优化和灵敏度分析

基于选择施加在结构“控制点”上的虚荷载作为优化设计变量,针对一种新的承受约束的形状优化数值方法进行了研究。借助于节点位移与虚荷载之间的线性关系,提出了一种新的计算灵敏度系数的解析方法。利用节点移动速度域概念构造了优化新形状产生的计算公式,以结构中节点的最大应力最小化作为优化目标,通过控制网格结点的最大位移量,较好地解决了单元网格在形状优化中的扭曲问题。对三个不同的实例成功地完成了形状优化。     

基于均匀化理论的管板有效弹性常数的研究

根据多尺度均匀化理论建立了管板有效弹性常数研究的有限元计算模型，在弹性平面状态下，计算结果与ASME规范中采用的有效弹性常数进行了比较，其相对误差在0．2％，然后对管板中胀接管子的加强作用进行了研究，得到了管板有效弹性常数与胀接管子在孔间带效率为0．4时的变化曲线。计算结果表明，用均匀化方法计算管板的有效弹性常数是可行的。     

基于物质点法的射弹侵彻双层隔板充液结构的数值模拟

在充液结构内设置双层隔板空气夹层可以有效减弱水锤效应带来的危害。为了研究双层隔板间距及位置对水锤效应的影响机理，通过已发表的试验数据对数值模拟方法进行验证后，基于物质点法开展了射弹侵彻双层隔板充液结构的数值模拟，分析了空穴形成过程、射弹剩余速度、固定点液体压力峰值、出入射壁及双隔板壁变形等情况。结果表明：随着双层隔板间距增大，充液结构的变形程度呈先减轻后加重的趋势；双层隔板位置越靠近入射壁，对压力冲击波传导的阻碍作用越强，充液结构的抗侵彻能力越好。