流体饱和多孔介质动力问题的显式时域解法

为了简化分析,Zienkiewicz等基于Biot理论,在忽略流体相对于土骨架运动的加速度条件下,建立了以土骨架位移u和孔隙流体压力p为基本变量的u-p格式饱和两相介质动力方程。针对该u-p方程,在空间上,采用伽辽金法有限元离散,并结合对角化形式的质量矩阵和流体压缩矩阵,忽略相邻结点间的惯性和流体压缩量间的耦合作用。在时域内,基于杜修力等提出的显式算法和Euler预估-校正法,建立了一种具有二阶精度的全显式时域积分法。采用一维饱和土模型,对比提出算法的数值解与Simon方法的解析解,发现两者吻合良好,验证了本文方法的正确性。并分析了饱和土二维动力问题,以及渗透系数和排水条件对饱和土动力响应的影响。     

一种新的高阶弹簧-阻尼-质量边界——-无限域圆柱对称波动问题

提出一种描述力-位移时间卷积关系的高阶弹簧-阻尼-质量模型,并将其作为人工边界条件直接应用于弹性动力学无限域圆柱对称运动问题的时域数值求解. 该人工边界条件不存在旁轴近似、多次透射等位移型外推人工边界条件普遍存在的高、低频失稳问题;与黏性、黏弹性边界等应力型人工边界条件相比,它具有高阶精度,且是严格高、低频双渐近的,也可以退化到黏性、黏弹性边界;该边界可以像黏性、黏弹性边界一样利用商用有限元软件中内置的并联弹簧-阻尼器、质量单元和时间积分求解器在商用软件中方便地实现,便于研究人员和工程师应用. 分析的几个简单数值算例也验证了该边界条件的上述优点.      

水下爆炸作用下基于声学的爆源子结构输入方法

无限域吸收边界和爆炸波输入是水下爆炸载荷数值模拟的两个关键问题.本文借鉴基于内部子结构的地震波动输入方法和爆源子结构多尺度分析方法,考虑水下爆炸载荷与地震载荷同属于波动问题,提出一种水下爆炸作用下爆源子结构的爆炸波输入方法,该方法首先通过对爆源区域进行水下爆炸波自由场的波场分解,并利用地震波的内部子结构输入方法,将该自由波场运动转换为等效爆炸波输入载荷,从而实现了水下爆炸问题中冲击波的输入;本文采用圆形爆源子结构区域,并采用AUTODYN软件中一维自由场模型计算水下爆炸作用下该区域的自由波场压力时程.进一步,基于连分式近似方法提出了一种模拟无限域水体辐射效应的高精度时域人工边界条件,该吸收边界可设置于离结构和爆源子结构较近处.本文提出的分析方法通过圆形爆源子结构进行水下爆炸载荷的转化,并采用高精度吸收边界大大减少计算区域,既保证了计算精度,又降低了单元数量,具有较高的计算效率和较强的实用性.通过数值算例分析,验证了本文模型与方法的准确性,模拟了水下爆炸作用下的冲击波和气泡脉动阶段测点的压力时程曲线,并研究了圆形结构对水下近场爆炸波散射效应的影响规律.     

钢筋混凝土界面的近场动力学方法

近场动力学方法已被广泛用于钢筋混凝土的开裂破坏研究,传统近场动力学方法的控制方程与参数是基于同种均质材料的能量方程确定,在处理不同种材料之间的相互作用时,无法合理反映其界面的力学行为.针对这一问题,通过分析钢筋混凝土界面的黏结-滑移机理,提出了近场动力学界面区材料点的相互作用模型,发展了考虑钢筋混凝土界面黏结的键基近场动力学方法.基于键基近场动力学与连续介质力学的能量密度等效方法,提出了界面微弹性参数的确定方法;根据钢筋肋间混凝土的应力分布规律,获得界面材料点域半径与受限楔体半径的等效关系;利用界面黏结-滑移曲线峰值应力对应的滑移变形,给出了界面临界拉伸常数确定方法.通过与2组钢筋混凝土构件的拉拔试验对比,验证了发展的界面近场动力学方法,并开展了不同条件下钢筋混凝土构件的数值试验.结果表明,发展的近场动力方法能够合理反映钢筋直径、锚固长度、混凝土强度以及肋间距对钢筋混凝土界面黏结行为的影响,体现了所提方法的合理性与优越性.