黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.      

船舶与海洋平台结构冰载荷的高性能扩展多面体离散元方法

船舶与海洋平台结构的冰载荷是寒区海洋工程结构物设计中的关键参数,而离散元方法是有效计算结构冰载荷的重要手段. 本文采用基于闵可夫斯基和原理的扩展多面体离散元方法模拟船舶与海洋平台结构的相互作用过程. 其中,构造扩展多面体的近似包络函数并建立了基于优化模型的快速接触搜索算法;考虑单元间粘结作用的刚度软化过程建标识码元间的粘结-破碎模型. 同时,发展了 CPU-GPU 协同异构环境下的高性能并行算法. 为分析海冰与海洋结构作用中的冰载荷,采用ISO标准验证了扩展多面体离散元分析结构冰载荷的准确性. 采用离散元方法计算了船舶结构的冰载荷,研究了船舶结构表明的线载荷分布特点,并采用船舶结构冰阻力经验公式验证了计算结果的合理性. 采用离散元方法计算了平整冰区与多桩腿平台结构的相互作用,分析各桩腿上的冰载荷特点. 针对碎冰区的海冰管理过程,采用离散元方法分析了船舶结构绕行过程中的船舶和海洋平台结构冰载荷. 本文方法可有效应用于海洋结构冰载荷分析,能为极地船舶与海洋平台结构的设计和安全运行提供科学的分析手段.      

含非完整界面的功能梯度压电材料Ⅲ型裂纹问题

本文研究了含非完整界面的功能梯度压电复合材料的Ⅲ型裂纹问题．此裂纹垂直于非完整界面,采用弹簧型力电耦合界面模型模拟非完整界面．界面两侧材料的性质,如弹性模量、压电常数和介电常数均假定呈指数函数形式且沿着裂纹方向变化．运用积分变换法将裂纹面条件转换为奇异积分方程,并使用Gauss-Chebyshev方法对其进行数值求解．根据算例结果讨论了一些退化问题并分析了裂纹尖端强度因子与材料的非均匀系数和非完整界面参数的关系．     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程,结合广义Betti-Rayleigh互易等式与时域下的边界积分方程,推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解,将经过主部分析的积分核函数分解为静态和动态部分,其中动态积分核不具有奇异性。在裂纹前沿附近单元,采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换,采用配置点法计算超奇异积分,获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序,得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律,数值结果稳定且收敛速度快。     

基于Voronoi切割算法的碎冰区构造及离散元分析

离散元方法广泛应用于海冰,特别是碎冰区的动力过程及其对海洋结构作用过程的数值模拟。为构造碎冰区中的冰块几何特性,基于二维Voronoi图方法对计算域进行随机切割以生成碎冰区中冰块的几何形态,并采用球体单元对每个碎冰块单元进行填充,从而确定碎冰区的初始分布场。在采用Voronoi图进行碎冰区构造时,可对冰块尺寸、几何形态和密集度等海冰参数进行设定。为确定冰块的不同几何规则度,综合采用排斥法和扰动法以定量地控制碎冰块几何形态从完全随机分布到规则分布的连续变换。为分析不同几何规则度下碎冰块的几何特性概率分布规律,对计算域内冰块的面积和边数等参数进行统计分析,从而可更合理地参数化控制初始冰场中碎冰块的几何特性。在此基础上,本文基于粘接-破碎的球体离散元方法对不同冰况下锥体结构的冰荷载进行了数值计算,讨论分析了碎冰区的海冰密集度、冰块面积和几何规则度对冰载荷的影响。     

颗粒流的动力学模型和实验研究进展

首先叙述了颗粒流研究的背景和基本概念,接着概述了颗粒流动力学研究的三种基本数学力学模型:颗粒动理论(kinetic theory)模型、摩擦塑性模型和离散元模拟模型.介绍了它们的基本原理,并分析了各个模型的优缺点和适用范围.随之,对颗粒流的实验研究情况作了简要的介绍,列举了研究中几种常用的试验方法和观测手段,以及观测到的一些典型现象.最后,我们简述了颗粒流研究的主要困难,并提出了某些有待解决的研究课题.      

界面裂纹问题的解析—变分解法及其在复合材料盖板胶接件中的应用

本文根据穆斯海里什维利求解各向同性平面问题与列赫尼兹基求解各向异性平面问题的广义复变函数理论与本征函数展开法,分析了复合材料盖板胶接件的分层问题,得到了精确满足所有基本方程、裂纹表面边界条件与层间连续条件的位移场与应力场本征展开式.进而利用分区广义变分原理满足裂纹表面以外的边界条件并由此确定应力强度因子.     

多工况多约束下离散变量桁架结构的拓扑优化设计

提出了一个多工况下受应力、位移约束的离散变量桁架结构的拓扑优化方法，给出了结构拓扑形式变更时的约束处理方法及杆件删除策略，使基结构设计空间的维数不断降低，达到最优拓扑，避免了奇异解的出现，在算例中指出了截面离散集和位移约束对最优拓扑的影响；算例给出了满意的拓扑优化解。     

一种求解三维热辐射输运方程的整体预处理迭代方法及并行计算

为三维灰体热辐射输运方程的隐式离散纵标方法发展一个整体预处理迭代方法并研制并行程序.该方法采用组装线性代数方程组策略,同时求出所有离散方向上的辐射强度.借助预处理的Krylov子空间迭代法,避免复杂网格上扫描算法可能遇到的死锁问题,能够提高健壮性和计算效率.空间离散上采用一阶迎风有限体积格式.数值实验测试变形六面体网格...