液饱和多孔介质中三维应力波的传播

为了深入研究液饱和多孔介质中应力波的传播，提出了三维两相细观计算模型．基于此模型。应用Galerkin余量法并计及流-固耦合界面的耦合效应，利用直接耦合的技术，开发了三维流-固混合显式动力有限元计算程序．在此基础上对冲击载荷作用下液饱和多孔介质中三维应力波的传播现象进行了数值模拟，并详细讨论了孔隙率，孔隙形状等因素对应力波传播主导波形的影响．     

受冲击弹性圆板动应力集中的数值计算

1 引言在纵向冲击载荷作用下,圆柱体及圆板内产生的应力波经传播及在其自由表面上的反射,将使其某些局部产生动应力集中,进而导致脆性破坏.Kolsky 曾通过冲击实验对此问题进行了较为系统的研究.笔者曾通过数值计算分析了Kolsky 的试验中圆柱内产生的各应力成分的变化情况及动应力集中现象.求得的动应力集中的发生位置同Kolsky的试验结果基本一致.本文采用同样的数值计算方法,以立体图形表示了受冲击弹性圆板中发生动应力集中的部位,并研究了不同的冲击载荷及其作用范围,圆板的尺寸及材料特性等因素对动应力集中的发生及其发生部位的影响.     

三维热弹性应力波传播问题的数值分析

为求解三维动态热弹性问题提出了一种基于特征线法的数值方法，并对矩形截面的热弹性棒在冲击载荷和热冲击作用下的三维应力波传播过程进行了数值模拟。     

离散元法研究的评述

介绍了离散元法的基本理论、计算方法及其应用的现状和最新进展.从离散元法的离散模型特点及便于甄别与其它数值计算方法的关系的角度给予离散元法一个比较宽松的定义.在此基础上阐明了离散元方法与刚体-弹簧模型(rigid body spring model, RBSM)方法,不连续变形分析(discontinuous deformation analysis, DDA)方法,分子动力学(moleculardynamics, MD)方法,三维离散元(discrete meso-element dynamicmethod, DM$^2$)方法及无网格方法(meshless method)等数值计算方法的关系, 并讨论了离散元法研究中亟待解决的问题和今后的发展方向.      

Analysis of High Accuracy Conservation-Element and Solution-Element Schemes

We propose a technique to construct higher accuracy conservation element and solution element (CE/SE) schemes.A second-order CE/SE scheme is established and its stability condition is examined on the basis of the von Neumann necessary condition. From the viewpoints of accuracy and efficiency, the applied range of the one-order CE/SE scheme and the high-order ones is discussed through some numerical examples.     

横观各向同性体应力波传播问题的特征分析

本文根据广义特征理论，对横观各向同性体轴对称弹性波传播问题进行了理论分析。提出了一个简化特征分析的方法，这种方法将成为应用特征线法求解种类各向异性应力波问题的一个突破口。     

脆性材料动态破坏过程的数值模拟

根据连续介质力学基本原理和离散元法基本思想，针对脆性材料动态破坏过程问题，推导出基于圆盘单元的正交各向异性离散元二维计算模型，并利用此模型计算了各向异性薄板在冲击载荷下的应力波传播问题。通过将计算结果与LS－DYNA程序相比较，表明了本文提出的离散元模型的精确性。另外，对钢弹侵彻下滑凝土圆板破坏过程这个典型的发生连续介质向非连续介质转化的动力学问题进行了数值模拟和动画显示，并同侵彻实验中出现的几种破坏形式进行了分析和对比，从而证明了本计算方法是计算和模拟脆性材料出现从连续介质向非连续介质转变的动态破坏问题的有力工具。     

弹粘塑性圆管中动应力集中的数值计算研究

本文采用沿次特征线积分方法,求解纵向冲击载荷作用下弹粘塑性圆管的动应力集中问题。数值计算结果表明:动应力集中的发生位置取决于圆管的壁厚与外环半径之比,圆管内各点最大拉应力中的极值依赖于冲击载荷的升载时间及圆管的材料特性。本文研究了避免发生动应力集中的条件,并同相对应的圆柱体的数值计算结果进行了比较。     

高速冲击问题的离散元法数值模拟

以能量等效为原则,构建了新的求解弹塑性轴对称问题的离散元模型。通过引入断裂准则,使模型不仅可以应用于连续介质问题,还可以应用于从连续介质向非连续介质的动态转化过程。最后,通过模拟钢板受冲击载荷产生层裂的过程,说明模型在高速冲击问题中的适用性及其有效性。     

离散元法在求解三维冲击动力学问题中的应用

提出了三维连结型离散模型,建立了可实现连结型模型(用于连续介质)-接触型模型(用于非连续介质)转化的三维离散元计算程序,用来模拟连续介质转变为非连续介质的力学过程．利用该计算程序对冲击载荷下混凝土块体内(连续体情况下)的应力波传播过程进行了数值模拟．将计算结果的数值与LS-DYNA程序计算的结果进行比较,验证了该计算程序的计算精度．在此基础上,模拟了混凝土块体的动态破坏(连续介质向非连续介质转化)过程．其计算结果可用动画显示,得到的破坏形式与由实验得到的破坏形式相近．两个算例说明该离散元模型及其计算程序是模拟计算伴随有连续介质向非连续介质转变的动态破坏问题的有力工具．