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

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

内爆荷载作用下钢筋混凝土安全壳的非线性响应分析

安全壳内部爆炸和恐怖袭击可能导致放射性物质的泄露,会造成严重的灾难。运用非线性有限元 软件LS-DYNA,采用分离式建模方法和流构耦合算法,建立钢筋混凝土安全壳有限元模型,对钢筋混凝土安 全壳遭受内部爆炸进行了数值分析,研究了内爆炸冲击波作用下安全壳内的动力响应,并与Henrych公式和 TM5-1300结果进行了比较。结果表明,在比例距离为3.258m/kg1/3、爆炸源距基底48m 的内爆作用下,数 值模拟能保持较好的精度。     

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

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

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

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

双颗粒沉降的动力学分析

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

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

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

压电材料平面裂纹尖端场的杂交应力有限元分析

基于复势理论和杂交变分原理建立了一种适用于力电耦合分析的杂交应力有限元模 型. 给出了建立刚度矩阵的主要公式和推导过程，单元内的位移场和应力场采用满足平 衡方程的复变函数级数解，假设的复变函数级数解事先精确满足裂纹的无应力和电位移法向 分量为零的条件，单元外边界的位移场假设按抛物线变化, 单元的刚度矩阵采用Gauss积分的方法得出. 通过对力电耦合裂尖场的数值计算验证了程序 的正确性和单元的有效性，同时也用所得结果校验了理论解.     

计算材料科学中桥域多尺度方法的若干进展

材料科学中存在固有的多尺度特性,桥域多尺度方法是在宏观尺度(如连续介质力学)中引入不同的细微观尺度的计算区域,乃至纳米尺度的分子动力学、量子力学计算区域,将不同尺度的研究方法通过一定的数学模型耦合在一起。该方法既能节约计算成本,又能保证所研究问题的物理特性。本文对多尺度方法的基本概念、跨尺度桥域多尺度方法的发展、基本原理、耦合方法和离散方程进行了讨论,给出了几个应用算例,并在最后进行了总结,展望了今后的可能发展方向。     

平面弹性力学问题的离散元法

根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.     

弹塑性摩擦接触问题形状设计灵敏度分析

提出了一种计算二维有限变形弹塑性摩擦接触问题形状设计灵敏度的算法. 采用主动集策略和mortar方法处理接触边线上的约束条件. 在mortar接触边线的切线和法线方向上采用相同的名义罚函数,提出基于名义罚函数的移动摩擦锥算法来正则化接触约束条件,发展了一种新的二维多体有限变形摩擦接触算法. 在此基础上, 通过将离散形式的摩擦接触问题控制方程对形状设计变量微分,得到了该路径相关问题的直接微分法解析设计灵敏度计算格式, 其节点位移灵敏度方程在每个增量步不用迭代、直接求解. 与国际上现有的二维多体有限变形摩擦接触问题的解析设计灵敏度算法相比,本算法不需分解为法向和切向推导,表达式较简洁,便于编程实现. 数值算例验证了算法的精度和有效性.