摩擦接触问题的比例边界等几何B可微方程组方法

摩擦接触问题是计算力学领域最具挑战性的问题之一,接触系统的泛函具有非线性、非光滑的特点,导致接触算法的收敛性与精确性难以保证.因此将比例边界等几何分析(scaled boundary isogeometric analysis,SBIGA)与B可微方程组(B dierential equation,BDE)相结合,提出了求解二维摩擦接触问题的比例边界等几何B可微方程组方法.在比例边界等几何坐标变换的基础上,通过虚功原理推导了关于边界控制点变量的接触平衡方程,表示成B可微方程组形式的接触条件可被严格满足,求解B可微方程组的算法的收敛性有理论保证.此比例边界等几何B可微方程组方法(SBIGA-BDE)只需在接触体边界进行等几何离散,使问题降低一维,能精确描述接触边界,并可通过节点插入算法进行真实接触区域的识别.此外,由于几何建模和数值分析使用相同的基函数,节约了划分网格的时间.以赫兹接触问题和悬臂梁摩擦接触问题为例,通过与解析解及数值计算软件ANSYS计算结果进行对比,验证了该方法求解二维摩擦接触问题的有效性及高精度等特点.      

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量,引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感,不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应,较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料,基于非线性塑性冲击波模型发展了耐撞性反向设计方法,以维持冲击物受载恒定为目标,运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型,并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明,反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的,所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

梁柱简化模型在自立塔塔线体系动态响应中的适用性研究

对输电塔进行合理简化可以提高塔线体系动力学仿真的效率。本文给出自立塔梁柱简化模型的计算方法,并提出利用梁柱简化模型计算方法建立自立塔塔线体系整体模型,同时采用桁梁混合模型建立精细化塔线体系整体模型,对两种模型塔线体系静力特性及振型和固有频率等动力特性进行对比分析。以脱冰工况为例,采用生死单元技术将施加在输电线节点上的集中质量单元杀死来模拟脱冰,实现对塔线体系动力学响应的有限元模拟,研究塔线体系简化模型在动态响应中的适用性。结果表明,两种模型弯曲变形误差小,低阶的振型相同,固有频率值误差小,动力特性基本相同;脱冰工况下,自立塔节点位移和塔材内力时程曲线一致,在提高计算效率的情况下,能有效保证计算精度。     

易拉罐在轴-侧-扭-内压联合作用下的屈曲地貌

柱壳结构广泛应用于各个领域, 但由于其对初始缺陷较为敏感, 容易发生灾难性的屈曲失稳. 本文利用非线性有限元分析程序ABAQUS研究了柱壳屈曲问题, 并应用到了易拉罐的屈曲分析. 首先采用数值模拟的方法验证了Virot等学者的易拉罐屈曲试验结果, 然后为了获得屈曲的一些普适结果, 进一步考察了柱壳的屈曲表现. 对柱壳结构在不同载荷组合、不同几何参数作用下进行了细致分析. 为了讨论的直观, 本文绘制了柱壳结构在受到侧压-轴压载荷作用下外力-屈曲载荷-位移三维屈曲地貌图(称为landscape). 结果表明: 在侧压-轴压-扭转载荷作用下, 试件力-位移曲线出现了"cliff"(断崖)现象; 扭转载荷的施加不利于试件整体稳定性, 并造成了试件对初始缺陷的敏感性; 对于受到轴压-扭转载荷作用的试件, 本文定义承载力为零的平面为"sea level"(海平面)来区分试件破坏模式; 通过对不同边界条件的试件进行分析, 发现试件两端固定可以有效地增加结构的承载能力, 提高稳定性. 对柱壳结构内部充气可以大幅度提升结构的承载能力和稳定性, 减小对初始缺陷的敏感度.      

超声弹性成像应用在角膜疾病诊断中的有限元分析

弹性成像在医学成像领域具有广阔的应用前景.在本文中,有限元法被用于模拟角膜组织的超声弹性成像,通过改变角膜病变组织与正常组织杨氏模量的数值大小及分布,模拟角膜病变的不同情形,分析计算各种不同情形时的应变、应力和位移分布,并对仿真结果进行分析.当角膜正常组织与病变组织的杨氏模量之比为1∶4时,病变区域中心的应变为-0.0...     

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

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

磁头/磁盘滑动接触下磁盘温度及热退磁临界条件的研究

采用二维轴对称有限元模型计算磁头/磁盘滑动接触下,铝质磁盘的稳态温度和热应力场以及热退磁临界条件.结果表明:磁盘温度在极短时间内升至摩擦稳态值,然后缓慢线性升高到最终稳态值;经过充分热传导和热交换后磁盘的温度梯度较小,此时磁层内的热应力集中分布于磁盘固定端边缘附近;磁盘的稳态温度和热应力均随速度增大而增大,且载荷越大其值增大越快;热应力小于1.2 GPa时所对应的速度和载荷为安全工况;温升大于373 K时所对应的工况将导致磁盘退磁.     

薄板小波有限元理论及其应用

利用样条小波尺度函数构造了常用的三角形和矩形薄板单元的位移函数,得到了利用小波函数表示的形函数。采用合理的局部坐标,对单元进行压缩,使单元在局部坐标区间上有其值,成功地推导出了分域的三角形和矩形薄板小波有限元列式。在此基础上,提出了弹性地基薄板的小波有限元求解方法。通过两个算例对薄板的挠度和弯矩进行了计算,数值结果表明,求解结果具有收敛快、精度高的特点。     

PVDF涂层建筑膜材单元结构力学性能的试验研究

使用PVDF涂层建筑膜材制成的结构,在承受外载荷情况下,力学行为非常复杂。在本文中,重点研究其单元结构在承受外载荷情况下的力学性能。首先,对PVDF涂层建筑膜材进行拉伸试验,测定其在常温和低温环境下膜的拉伸力学性能参数。然后,以制成的三角框架结构膜单元为对象,在常温条件下,从有限元分析和实际加载试验两方面计算得出结构单元在均布载荷下的位移、应变和应力分布。最后,对其进行强度验证,得出采用膜单元构建大型建筑是安全的结论。试验方法和结果可以应用于膜建筑结构的设计。     

转向剂形成人工应力遮挡的实验研究及有限元分析

针对水力压裂过程中采用的人工应力遮挡控制缝高技术,通过室内物理模拟实验,在转向剂浓度一定的条件下考虑携带液与地层压差、渗流量、时间等因素,采用电测法对人工岩芯表面应变进行了测试;利用非线性有限元方法对裂缝尖端应力进行计算,研究了转向剂形成人工应力遮挡的作用效果。结果表明采用转向剂可以建立3~4MPa的人工应力遮档,对采用转向剂的压裂设计和施工参数提出了建议。