正交各向异性含液饱和多孔介质中应力波的传播

基于Biot的孔隙介质理论，研究了正交各向异性含液饱和多孔介质中应力波的传播特性．本文引入动态渗透率，导出了整个实频域内应力波传播的复特征方程及其解析解，给出了各种应力波成分的波速和衰减的解析表达武，计算了频散曲线和衰减曲线，并讨论了各类应力波之间的耦合关系及介质的各向异性对应力波传播的影响．     

移动简谐荷载作用下桥梁响应的高效计算

在计算移动荷载过桥问题中广泛使用的Newmark方法必须在每一时间步内限制荷载的大小和作用位置都不能改变。精细积分法虽然允许荷载的大小在每一时间步长内发生变化,但是仍假定其作用位置是不变的,未能采取措施以描述荷载沿着桥面的连续移动性。本文提出三种精细积分格式,在每一时间步内不但允许移动荷载的大小按简谐规律连续变化,而且模拟了简谐荷载在空间域的连续移动。通过与Newmark方法和简单问题的解析解进行数值比较,表明用本文提出的方法可以用较粗的结构单元和较大的时间步长而获得很高的计算精度。在精度相同的前提下,计算效率比Newmark方法可提高1~2个数量级。     

饱和-非饱和土壤中污染物运移过程的数值模拟

本文提出了一个模拟饱和 非饱和土壤中溶和污染物运移过程的数值模型．模拟的控制污染物运移的物理 化学现象包括：对流，机械逸散，分子弥散，吸附，蜕变，不动水效应．发展了一个修正的特征线Ｇａｌｅｒｋｉｎ方法以离散污染物运移过程的控制方程并导出了一个用于有限元方程求解的显式算法．数值例题结果表明所提出模型和算法的功能     

热结构稳态响应的耦合灵敏度分析方法

研究结构稳态热变形和热应力的灵敏度分析方法，给出了直接法和伴随法两种算法。考虑了温度场的耦合作用，在直接法中需要计算温度场对设计变量的导数，在伴随法中需要计算热载荷对温度场的导数。对尺寸和形状两类设计变量的灵敏度分析算例，验证了本文方法的精度。伴随法在应用程序中的实现，为大型结构优化提供了高效率的灵敏度计算方法。     

在轴向应力波传播和反射过程中弹性有限长圆柱壳非对称动态屈曲

本文讨论弹性有限长圆柱壳端部受冲击载荷作用，在轴向应力波传播和反射过程中的非对称动态屈曲问题。通过建立和求解扰动方程得到了动态屈曲的分叉条件，临界载荷和屈曲模态。数值结果表明：当壳壁厚不很薄时，轴对称屈曲临界载荷比非轴对称临界载荷高；反之，轴对称临界载荷会比非对称临界载荷低；由于应力波的反射，临界载荷降低，因而更容易发生屈曲，屈曲模态也有其不同特点。     

非齐次动力方程Duhamel项的精细积分

提出了不需要矩阵求逆运算的求解Duhamel积分项的精细积分方法.通过将精细积分法的关键思想--加法定理和增量存储--直接应用于Duhamel积分响应矩阵的求解,可给出当非齐次项分别为多项式、正弦/余弦以及指数函数等基本形式时Duhamel积分在计算机上的精确解.特别的,该算法不依赖于系统矩阵(或相关矩阵)的形态.当系统矩阵奇异或接近奇异时,其优越性更为显著.算例验证了该算法的有效性.     

旋转运动中弹性梁耦合振动的辛方法

研究旋转梁结构的弹性耦合振动问题。通过引入对偶体系，建立了解决该类问题的辛方法。在辛体系中描述旋转梁纵向和横向耦合振动控制方程，即哈密顿正则方程。进一步求解得到结构的固有振动频率及相应的振动模态，发现固有振动频率随转动角速度先升后降以及模态之间的某种转化规律。     

粘弹塑性统一本构模型的参数评估

针对粘弹塑性统一本构模型参效众多的特点和传统评估参效方法的缺点，应用非线性最小二秉法与敏感系数分析方法相结合，提出了一个比较系统的对粘弹塑性统一本构模型的参数进行评估的方法。对铜合金在室温下不同加载速率的试验曲线进行参数适配，表明本方法的合理性与正确性。为粘弹塑性统一本构模型的工程应用奠定了基础。     

纳米-Al2O3陶瓷粉末的预热爆炸压实实验研究

对纳米-Al2O3陶瓷粉末在预热温度为0.5Tm的条件下,以不同的爆炸压实压力实施了烧结实验。通过X射线衍射分析了烧结体的晶型,并用高分辨率扫描电子显微镜进行了微观组织观察。实验结果表明,纳米-Al2O3陶瓷粉末在爆炸压力为13.1 GPa时可得到晶粒间结合致密且晶粒度在100 nm左右的-Al2O3烧结体;而在爆炸冲击压力为9.35 GPa以下时,烧结体的晶型没有发生转变,仍然是-Al2O3晶型。     

海冰动力学数值模拟中改进的PIC方法

为了准确地模拟海冰的动力过程,需要建立精确有效的数值方法。本文结合质点网格法(PIC)和光滑质点流体动力学方法(SPH)发展了一种改进的PIC方法。该方法在欧拉坐标下对海冰动量方程进行差分计算,在拉格朗日坐标下进行海冰质点位移、厚度和密集度计算,并采用Gauss函数进行欧拉网格点与拉格朗日质点间海冰参数的交互插值。采用改进的PIC方法对规则区域内的海冰堆积过程进行了数值试验,对渤海海冰的动力过程进行了72小时数值模拟。计算结果均表明改进的PIC方法具有计算量小,计算结果平稳精确的优点,可很好地适用于海冰动力作用过程的数值模拟。