延性单晶非均匀变形的三维有限元模拟

对延性单晶在拉伸载荷作用下的应变局域化和颈缩等非均匀变形过程进行了三维有限元数值模拟。将相关晶体塑性本构模型及一种新的数值积分方法补充到ABAQUS6.1商用有限元软件中。该方法的特点是，利用晶体塑性的动力学方程，获得一个关于晶体弹性变形梯度的演化方程，采用半隐式积分方案进行求解。本文推导出一种新的应力变本构矩阵。按此方式更新本构矩阵，计算速度和计算稳定性大大提高。加载方式，边界条件和变形程度等因素影响着滑移系的启动状况，这是平面模型所不能预测的。本文利用三维有限元方法模拟了不同取向下滑移系的启动状况，全面地考虑了FCC单晶材料12个可能滑移系在变形过程中的启动状况，合理地模拟了FCC面心立方单晶沿不同取向加载时晶轴旋转导致的应变局域化和颈缩等非均匀变形过程。     

多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程,是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．      

多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程，是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．     

椭圆截面曲线管道内二次流动的Galerkin解

推导了任意空间曲线直角坐标系内的张量流体力学方程．采用Galerkin方法和计算机符号运算技术求解了椭圆截面螺旋管道（电括弯管、扭管）内的二次流动．计算结果表明了Galerkin方法的适用和有效性，克服了振动法的小参数局限．对所得结果的分析揭示了椭圆截面螺旋管道内流动的特性，给出了k,r和Re较大情况的研究结论．     

水岩耦合变形破坏过程及机理研究进展

坝址水力劈裂、煤矿和隧道突水以及水力压裂工程都是岩体损伤演化积累诱发破坏 (灾变) 的过程,然而能否正确认识和评价这些过程取决于人们对水岩耦合变形破坏过程和机理的研究程度. 阐述了岩石渗流 --变形 --破坏特性及机理研究的进程, 侧重于物理实验成果. 综述了水岩耦合变形破坏过程分析模型和数值计算, 特别是有限单元法处理裂纹的思路, 并介绍了作者课题组近年来所取得的一些代表性成果. 给出了当前理论 (数值) 研究中所面临的几个具体问题和挑战, 并对今后研究工作进行了力所能及的展望.      

跨断层埋地管线-土接触非连续变形分析

跨断层埋地管线系统主要由管线及其周围土体两种不同介质组成,具有分析介质不连续性特性.为实现对管-土接触介质的不连续性分析,采用非连续变形分析与有限元分析相结合的方法,将跨逆断层埋地管线系统从实际工作状态中取出一部分作为分析对象.利用有限元分析方法对管线及其周围土体进行网格划分;而管-土之间的相互作用,利用非连续变形分析中的不连续介质接触处理方法实现模拟.通过模型的数值分析,研究了逆断层作用下管-土之间的非连续变形相互作用状态,验证了利用非连续变形分析与有限元方法相结合解决埋地管线-土接触的可行性和有效性,为管-土相互作用分析提供了新的研究思路和研究方法.     

花岗岩单轴压缩侧向变形及脆性破坏机制实验研究

为研究花岗岩侧向变形及脆性破坏机制,对花岗岩试件进行单轴压缩实验。利用动态应变采集系统、数字散斑相关方法(DSCM)和显微观测手段,记录并分析花岗岩试件在单轴压缩过程中的宏观侧向应变、局部侧向应变以及破裂面形貌,并与水泥砂浆试件的破坏过程对比,讨论了花岗岩脆性破坏机制。实验与分析结果表明:(1)花岗岩试件在加载初期发生侧向收缩变形,产生并发展于压密阶段,消失于线弹性阶段初期,这主要由于试件内部裂纹闭合造成的;此后,宏观侧向应变持续增长,当侧向应变与轴向应变之比接近0.5时试件破坏;(2)在峰值载荷前很长一段时间内,局部侧向应变在一定范围内波动,临近试件破坏时局部侧向应变最大值和最小值均出现较大幅度的波动,二者差值迅速增大,试件不均匀程度增大,最终导致试件破坏;(3)在峰值载荷前有无塑性屈服阶段是峰值载荷后脆性破坏程度的重要影响因素,而宏观裂纹的贯通程度是峰值载荷后应力降大小的决定因素。     

数字图像相关方法中匹配及过匹配形函数的误差分析

基于图像子区的数字图像相关方法需采用合适的形函数来近似目标图像子区的真实变形。由于实际测量时目标子区的局部变形往往是未知的,实际采用的不同阶次(零阶、一阶和二阶)的形函数不可避免地产生误匹配(欠匹配和过匹配)问题,从而引入位移测量的系统或随机误差。尽管由欠匹配形函数引起的系统误差已被充分认识,由过匹配形函数引起的位移误差仍缺少理论解释。本文首先推导出采用一阶和二阶形函数的数字图像相关方法的随机误差理论公式,随后采用一系列数值实验验证了理论公式的准确性。结果显示:过匹配形函数不会引入额外的系统误差,但会增加随机误差,且二阶形函数的随机误差是一阶形函数的二倍。考虑到由欠匹配一阶形函数引入的系统误差往往远大于过匹配二阶形函数的随机误差,因此在未能确知变形的情况下,推荐使用二阶形函数。     

双颗粒沉降的动力学分析

采用格子Boltzmann-虚拟区域方法对雷诺数范围为50≤Re≤200的双颗粒自由沉降进行了直接数值模拟。首先,研究发现双颗粒最终沉降的位置分别在通道的1/4和3/4位置附近,且与颗粒的初始间距、雷诺数以及通道宽度无关。其次,重点研究了颗粒沉降过程中所受的侧向力(与沉降方向垂直),首次揭示了侧向力的振动频率与雷诺数呈二次关系,且单颗粒的结果始终小于双颗粒的结果;研究还发现侧向力的振动频率与通道宽度近似呈幂律函数关系,且幂律指数与雷诺数有关,雷诺数越大,幂律指数的绝对值越小。最后还研究了雷诺数及通道宽度对侧向力振动振幅的影响。     

刃型位错和圆弧形刚性线夹杂的干涉效应

位错和夹杂的干涉效应对于理解材料的强化和韧化机理具有十分重要的意义。文中研究了晶体材料中刃型位错和多条共圆弧刚性线夹杂的干涉作用。利用Riemann—Schwarz反照原理和复势函数的奇性主部分析技术，得到了问题的一般解答；对于只含一条刚性线夹杂的情况，给出了复势函数的封闭形式解。由Peach-Koehler公式求出了作用在刃型位错上的位错力，并讨论了圆弧形刚性线夹杂对位错力的影响规律，发现弧形刚性线对刃型位错有很强的排斥作用。本文解答不但可作为格林函数获得任意分布位错的相应解答，而且可以用于研究刚性线夹杂和任意形状裂纹的干涉效应问题。