卵形杆弹对铝靶的斜侵彻

在Lagrange有限元基础上,扼要介绍了一种适用于三维斜侵彻数值模拟计算的滑移面处理技术。该方法放弃了传统滑移计算中从单元的设置,而代之以预设从节点,从而避免三维斜侵彻数值计算中主从单元的相交、滑移面的识别、修正与再定义的困难。在确保计算精度的同时,有效地提高了计算效率。卵形杆弹对铝靶侵彻的系列数值模拟计算表明,无论对于正撞击或斜撞击,本文中所获结果与实验结果都有良好的一致性,这说明本文中所述方法和所建程序的合理性和有效性,为侵彻贯穿过程的数值分析提供了一种实用和有效的手段。     

应用杆束系统研究水泥砂浆节理面的压剪动特性

基于分离式Hopkinson压杆系统提出一套Hopkinson束杆装置,研究水泥砂浆节理面在压剪复合加载下的动态界面滑移特性。应用三根小杆作为接收杆,分别采用单杆和三杆两种入射方式,对含有节理面的水泥砂浆试件进行冲击加载。在节理面发生滑动时,由三根接收杆测试得到水泥砂浆试件不同位置滑移状态,进而得出节理面的整体滑移速度。此研究初步揭示了节理面在发生滑移时的局部滑移状态以及整体滑移状态,为节理面的动摩擦特性研究提供了较可靠的实验技术。     

一种基于物理机制的后继屈服面演化模型

考虑滑移为晶体塑性的主要变形机制,推广单晶滑移理论至多晶体,建立滑移构元模型,研究了后继屈服面的移动、变形等演化规律。基于滑移构元的包氏效应和潜在硬化的影响,后继屈服面的尾部缩进,变得较“平坦”,而前端突出,形成“高曲率区域”。为了验证模型的有效性,本文分别预测了单轴拉伸、纯扭转和拉—扭组合比例/非比例加载下的后继屈服面演化过程,与已有实验结果符合一致。关键词:滑移构元模型,后继屈服面,包氏效应,潜在硬化,非比例加载      

输电塔螺栓搭接滑移过程的数值模拟

螺孔直径一般都大于螺栓外径,因此当横向载荷超过一定数值时,螺栓连接件之间将发生相对滑动,这一过程就是螺栓滑移.在输电塔的结构分析中一般不考虑复杂的螺栓滑移,以致在煤矿采空区或冻胀区等可能导致输电塔塔基非均匀沉降的地区中,传统结构分析软件未能很好地判断出输电塔的运行状态,因此对螺栓滑移过程的准确分析是十分重要的.本研究通过有限元方法,建立输电塔角钢螺栓搭接结构的三维有限元模型,对已有的单螺栓滑移试验进行数值模拟,首次通过数值模拟成功还原已有试验的结果,并进行了搭接结构的滑移试验和模拟验证.通过数值模拟及与已有试验结果的对比,分析了影响滑移曲线的因素,并初步给出各模拟参数的取值范围.     

基于滑移理论的超疏水表面减阻性能分析

采用滑移理论计算和数值模拟方法研究超疏水表面的减阻性能.层流流动状态的数值模拟和理论计算结果的一致性较好,并可证明滑移速度与通道内流体速度之比与无量纲压降比的数量级相当;湍流状态的数值仿真结果表明,在目前可实现的滑移速度范围内,超疏水表面对水下航行器的流体阻力影响不显著.     

含摩擦滑移铰平面多刚体系统动力学的数值算法

基于LuGre摩擦模型和线性互补问题(LCP)的数值算法,给出了具有双边约束含摩擦滑移铰平面多体系统动力学的数值算法.首先,根据滑移铰的特点,当间隙充分小时,将其视为双边约束,给出了滑移铰中滑道作用于滑块上的法向接触力的互补关系;LuGre摩擦模型能有效地描述机械系统中的黏滞与滑移运动,将该模型用于描述滑块与滑道间的摩擦力.其次,结合Baumgarte约束稳定化方法,应用第一类Lagrange方程,建立了该多体系统的动力学方程,给出了Lagrange乘子与滑移铰中作用于滑块上的法向接触力的关系式.然后,将滑块与滑道间多种接触状态的判断以及作用于滑块上的法向接触力的计算转换为线性互补问题的求解,并用常微分方程的数值算法求解该多体系统的动力学方程.最后,通过数值仿真算例揭示了滑移铰中滑块的黏滞与滑移现象,以及滑块在滑道内的多种接触状态;另外,在文中分别采用Coulomb干摩擦模型和LuGre摩擦模型,对算例中的某些工况进行了数值仿真,并且分别用本文方法得到的数值仿真结果与已有方法得到的数值仿真结果对比,表明了本文给出的方法的有效性.      

含摩擦滑移铰平面多刚体系统动力学的数值算法

基于LuGre摩擦模型和线性互补问题(LCP)的数值算法,给出了具有双边约束含摩擦滑移铰平面多体系统动力学的数值算法.首先,根据滑移铰的特点,当间隙充分小时,将其视为双边约束,给出了滑移铰中滑道作用于滑块上的法向接触力的互补关系;LuGre摩擦模型能有效地描述机械系统中的黏滞与滑移运动,将该模型用于描述滑块与滑道间的摩擦力.其次,结合Baumgarte约束稳定化方法,应用第一类Lagrange方程,建立了该多体系统的动力学方程,给出了Lagrange乘子与滑移铰中作用于滑块上的法向接触力的关系式.然后,将滑块与滑道间多种接触状态的判断以及作用于滑块上的法向接触力的计算转换为线性互补问题的求解,并用常微分方程的数值算法求解该多体系统的动力学方程.最后,通过数值仿真算例揭示了滑移铰中滑块的黏滞与滑移现象,以及滑块在滑道内的多种接触状态;另外,在文中分别采用Coulomb干摩擦模型和LuGre摩擦模型,对算例中的某些工况进行了数值仿真,并且分别用本文方法得到的数值仿真结果与已有方法得到的数值仿真结果对比,表明了本文给出的方法的有效性.     

表面荷载作用下半无限层状介质中的裂纹分析

分析了半无限层状介质中的正方形裂纹。层状材料的层面互相平行,外部荷载作用在边界面上,正方形裂纹平行于层面。基于Yue基本解的数值方法和线弹性断裂力学叠加原理,首先采用一种数值方法获得无裂纹半无限层状介质的应力场,然后将计算得到的应力按叠加原理施加在裂纹面上,并采用另一种数值方法计算此情形下裂纹面的间断位移,最后采用裂纹面的间断位移计算应力强度因子。结果显示:I型和II型应力强度因子的变化与裂纹所处的位置关系密切;层状介质中的裂纹张开和滑移受到不同介质存在的影响,进而影响到裂纹的应力强度因子。建议的数值方法可用于分析复杂荷载作用下层状介质中裂纹的断裂力学特性。     

可压缩材料平面应变滑移线场的数值解和边界条件

简要介绍了可压缩材料平面应变滑移线理论。导出了滑移线场应力和速度数值解。处理了应力边界条件，阐明了应力和速度间断规律。举例说明了边界条件处理和数值解步骤。     

侧壁边界条件对平面二维模型流态的影响

侧壁无滑移、部分滑移与完全滑移三种不同的侧壁边界模拟方法将影响边界约束水域的潮流流态模拟结果.通过平底矩形水池的流场数值模拟,分析了侧壁无滑移与完全滑移两种情况下水池内的流态,讨论了水池入水口附近产生旋涡以及在侧壁无滑移时产生次生旋涡的机理,特别是侧壁部分滑移时不同侧壁粗糙度下水池内产生旋涡及次生旋涡的现象.给出了不同...     

基于边界元法求解三维弹性摩擦接触问题

边界元法求解三维摩擦接触问题,其中一个关键点在于如何确定滑移方向。即当出现相对滑移时,滑移方向如何确定。当前常采用的方法是,粘结点利用切向面力得到滑移方向,滑移点利用切向相对位移得到滑移方向。不过该方法难以保证收敛性。针对这一问题,本文采用滑移方向预测技术得到滑移方向。即以后出现相对滑移时,滑移方向采用预测技术中得到的滑移方向。由于摩擦接触问题和历史加载相关,本文采用增量法求解。不同摩擦系数下的数值结果都证明了本文算法的有效性和收敛性及滑移方向预测技术的有效性。     

一般各向异性单侧接触界面上波的反射和折射

研究简谐弹性波在一般各向异性介质单侧接触界面上的反射和折射问题．利用Fouier分析方法将非线性Coulomb摩擦接触边界波动问题化为一组代数方程．给出了确定局部分离、滑移和粘着区的思路和方法及各区域的解；讨论了出现界面局部分离和滑移的条件．对特定材料组合情况进行了详细数值计算，给出了界面力、相对滑移速度、张开位移、高频谐波的反射折射系数等特征参量；考察了平面和反平面波动的耦合及整体滑移等．其中关于高频谐波的结果可对已有实验结果给出很好的定性解释．在大多数情况下，即使对摩擦系数无穷大的粘滞接触界面，分离区端部也总是存在一个很小的滑移区。     

铜单晶拉伸试样表面滑移带痕迹的晶体塑性分析

采用建立在晶体塑性理论基础上的晶体塑性有限变形计算方法,针对铜单晶试样单轴拉伸过程中晶体滑移在试样表面留下的滑移带痕迹进行了数值研究.作者利用三维有限元模拟不同取向铜单晶试样的拉伸变形,通过晶体塑性滑移面与试样表面交线的几何分析,得到了试样在不同晶向拉伸下不同滑移系启动造成的试样表面滑移痕迹,并对数值计算的试样表面滑移痕迹作了初步讨论.所得结果表明晶体塑性理论能够用于单晶试样拉伸试验的表面滑移带痕迹形成的分析.     

纳米尺度下烷烃润滑薄膜分层现象和速度滑移机制的分子模拟研究

采用非平衡分子动力学模拟方法,研究了纳米尺度下十六烷润滑薄膜的分层现象和速度滑移现象,重点考察了剪切速度对速度滑移现象的影响规律,并对其中的微观机理作出了解释.研究结果表明:在铁壁面的限制情况下润滑薄膜出现了分层现象,当润滑薄膜厚度超过50?时,其中间区域呈现出体相均质流体的特征.润滑薄膜层间滑移和界面滑移的临界剪切速度分别为5.5和7.5?/ps,随着剪切速度增加,界面滑移程度增强,而层间滑移程度减弱.润滑薄膜第1和2层十六烷分子层间桥接分子数目决定层间滑移程度,随着剪切速度增加,桥接分子数目也相应增多,层间滑移程度随之减弱.     

过渡区气体微尺度流动的格子Boltzmann模拟

采用格子Boltzmann方法(LBM)对过渡区内的微尺度气体流动进行了模拟研究. 在已有滑移区微流动LBM模型中引入Knudsen层速度修正,选取合适的修正函数表达式并依据动理论确定了可调参数的合理取值. 在边界条件的处理格式上,采用了适合过渡区模拟的高阶滑移边界的替代格式来捕捉过渡区微流动的滑移速度,避免了直接求解高阶速度导数项的数值困难. 通过对两类不同的微流动进行模拟的结果表明:与数值解吻合得较好,尤其是对Kn>0.5微流动滑移速度的预测,与已有LBM的模拟结果相比有明显的提高.      

基于非结构任意多边形网格的滑移面计算技术

基于任意多边形网格管理体系,针对流体多介质问题的数值模拟,发展了拉氏方法滑移面计算技术.文章给出了滑移线设置的数据结构,滑移线上主从点速度与位置的计算格式,及节点滑移后引起界面上点、相关网格邻域关系变化的算法.该滑移计算技术避免了传统算法中由于以模拟法(重叠或分离网格)代替直接法(拼接网格)而造成几何守恒律被破坏的缺陷.数值例子验证了该算法的可行性,体现了算法无缝连接的特点.     

边坡大距离整体滑移的模拟研究

边坡整体滑移模拟是迄今未被解决的一个课题。文章通过构制一个专门的滑面单元, 采用数值模拟的原理实现了对粘滑类边坡失稳的大距离滑移模拟; 采用分析模拟法实现了对脆性滑移类边坡失稳滑移的模拟, 并开发了相应的软件。     

纤维增强混凝土材料的界面剪应力分布研究

针对纤维与混凝土界面的破坏过程,提出了几种简化的粘结-滑移本构模型,以双线性局部粘结-滑移本构模型为基础,在受力平衡和变形协调的基本原理基础上,推导了纤维脱粘过程中界面剪应力的解析解.采用弹簧粘结单元,通过数值方法模拟了纤维与混凝土之间的粘结-滑移过程,给出了纤维与混凝土界面脱粘过程中界面剪应力的分布、变化情况.对解析解、有限元计算结果和试验结果之间的差异进行了对比分析,验证了简化模型的合理性和有效性.     

地下爆炸波作用下基底滑移隔震大跨结构的动力响应分析

对由地下隧道内发生意外爆炸引起的地下爆炸波作用下基底滑移隔震大跨结构考虑土一结构相互作用（SSI）的动力响应进行了全面的分析。在SSI分析中,分别采用有限元、无穷元和杆件体系建立了地下衬砌及周围有限土区、远场无限土区以及基底滑移隔震大跨结构的计算模型;在基底滑移隔震分析中,提出用连续摩擦力模型取代传统库仑摩擦力模型,以避免跟踪啮合滑移的过渡边界并可减小计算误差。最后,通过与Lysmer透射边界法和修正阻尼透射边界法的算例结果对比分析,验证了本文方法的有效性,其数值结果揭示了地下爆炸波作用下大跨结构的基底滑移隔震效果以及不同隔震参数对隔震效果的影响,从而为城市地下隧道建设相关的理论研究和工程实践提供了重要的依据。     

基于能量极值原理的单晶体塑性滑移的有限变形分析

对FCC单晶体的率无关弹塑性力学响应的本构关系进行了数值模拟.用一个基于能量极值原理的数值计算方法来处理复杂的多面塑性问题.这种算法可以有效地模拟单晶体多滑移系的启动.在这一理论框架下,增量的应力应变关系可以从所构造的能量函数中推导出来.对于滑移系激活情况的判定则可转化为求解活动约束的非线性数学规划问题.通过对时间的离散,此问题又可细化为逐步二次规划问题,并采用有效集法来搜索启动滑移系,进而求得弹塑性本构关系.数值结果表明该方法具有稳定、收敛、可行的特点.在数值计算的基础上研究了单轴拉伸下晶体的不同取向对单晶体硬化程度和滑移系激活情况的影响.