《数值流形法》序

得知郑宏教授的《数值流形法》即将脱稿,激动之情难以言表.我完全同意郑教授将流形法的特点归结为“升阶”和“切割”.受郑教授的委托,我愿在此将建立流形法的历程做一简单的回忆.其实早在1985年———也就是在我完成不连续变形分析法(DDA)之前,“升阶”在我心中已完全成熟.“升阶”似乎能取得一致收敛,弥补有限元只能按能量收敛的缺憾.我没有急于发表是因为一直没能想到令我满意的解决“不连续”的办法.起初曾想过在DDA块体上引入高阶多项式,或在块体上再划网格,也就是用分片多项式来提高DDA块体的变形精度.但我对这些方案都不甚满意,因为实际计算时多项式的阶数会受到限制,引入网格又会招致有限元既有的毛病——解的精度受控于网格质量.     

一个基于非单调技术的超记忆梯度法

本文给出一个修正的非单调线搜索策略,并结合该策略提出一个求解无约束优化问题的超记忆梯度算法.该算法的主要特点是:在每一次迭代中,它所产生的搜索方向总是满足充分下降条件.这一特性不依赖于目标函数的凸性以及方法所采用的线搜索策略.在较弱的条件下,该方法具有全局收敛和局部R-线性收敛性.数值实验表明了该方法的有效性.     

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。     

阴离子-π复合物(I-·C6F6)的成键分析

最近Anst?ter, et al. 发表了对阴离子-π复合物(I^-·C₆F₆)的第一个定量谱学测量[J. Am. Chem. Soc. 141, 6132 (2019)]，认为成键作用中相关作用占41%，静电作用占23%，得出相关作用占主导的结论. 本研究表明，该文献的“静电作用”中混入了Pauli排斥作用，后者在数值上抵消了前者的作用. 在复合物I^-·C₆F₆中，发现静电作用是相关作用的两倍多，因而阴离子-π复合物中仍应是静电作用占主导.     

光纤布拉格光栅传感器应变传递双向耦合分析

粘贴于受弯基体表面的光纤布拉格光栅传感器测量应变与基体真实应变之间存在误差,因此研究光纤布拉格光栅传感器的变形机理、分析测量应变与真实应变之间的关系是目前的研究热点.首先研究光纤布拉格光栅传感器与基体之间的相互作用机理,然后,利用有限元解、实验值和理论解进行对比验证,并分析产生误差的原因.最后,通过参数分析研究弹性模量、厚度、粘结长度等参数对光纤布拉格光栅传感器测量效果的影响.结果表明有限元解、实验值和理论解具有相同的变化趋势,有限元解与理论解的误差在2%以内,测量值与理论解的误差在7%以内.平均应变传递率随着基体弹性模量的增大、粘结长度的增长而逐渐增大,随粘结层弹性模量的减小、粘结层厚度的增大而逐渐减小.该理论对应用于受弯基体应变测量的光纤布拉格光栅传感器的设计具有一定的指导意义.     

平面波输入下海水-海床-结构动力相互作用分析

海洋工程结构的地震反应分析是保证海洋工程结构地震安全的重要环节.由于其所处的复杂环境, 该问题涉及到流固耦合和土-结相互作用.本文基于海水、饱和海床、基岩流固耦合统一计算框架,采用Davidenkov模型和修正的Masing准则考虑饱和海床的非线性,在脉冲SV波垂直入射下, 进行了海域场地和海洋工程结构的动力响应分析. 首先,对比分析了线性自由场和非线性自由场输入情形的海域场地非线性反应,结果表明线性自由场输入时反应不合理,自由场分析和场地分析应该采用相一致的本构模型. 然后,对比分析了海床分别为线性和非线性情形时,海域场地以及海水-海床-结构体系的反应特征. 与线性海床情形相比,非线性对海床反应的影响主要由如下两方面因素控制: 一方面,非线性导致饱和海床模量减小, 饱和海床与基岩间的波阻抗比减小,由基岩到饱和海床间的反射系数和透射系数增加, 导致反应增大; 另一方面,非线性导致阻尼加大, 使海床反应减小. 对于本文算例而言,阻尼对非线性海床结果的影响占主导作用.      

考虑制动力影响的大型客车舒适性分析

利用有限元法和动力平衡原理,建立了具有层间接触的沥青路面和13个自由度的大型客车人体三维模型。以不平顺作为激励,随机模拟车流分布,在不同制动情况下,从时域和频域两方面分析人-车-路耦合振动下车辆和人体的动力响应。时域分析采用客观评价标准和主观烦恼率相结合的方法,频域分析考虑人体共振、人的心理和生理因素。结果表明:随着制动力增大,人体的竖向加速度幅值不变,俯仰、侧倾加速度幅值都增大,但是俯仰加速度增幅远大于侧倾加速度的情况,人的舒适性变差,烦恼率增高;在紧急制动时,主频段对人体共振、心理和生理产生的影响很小,但是次主频段与人体某些器官固有频段重合,对人体共振、心理和生理产生的影响是不容忽视的。     

水平井降压法和热激法水合物开采对海底边坡稳定性的影响

天然气水合物广泛赋存在深海沉积物孔隙中,被认为是具有巨大开发潜力的未来绿色能源之一,引起全球的关注.深海水合物开采将造成含水合物储层的强度劣化,可能产生孔压积聚,诱发海床失稳.本文基于边坡稳定极限平衡分析框架,引入考虑水合物开采热-流-化学耦合过程的数值分析模型,研究水合物开采对海底边坡稳定性的影响.采用TOUGH+HYDRATE热-流-化学耦合分析程序,模拟了采用水平井降压法和热激法开采深海水合物的过程,分析了水合物分解锋面扩展和瞬态孔压演变的规律,并通过SLOPE/W程序采用极限平衡分析方法计算水合物开采过程及停采后的海底边坡安全系数,分析开采井位置和开采方法对海底边坡稳定性的影响.研究表明,对于存在致密盖层的细砂储层陡坡,单水平井降压开采过程中,由于孔压降低,土体有效应力增加,边坡稳定性显著提高,当开采井布设在坡体中部时,边坡稳定性提高最为明显;停采后,由于水合物分解导致土体黏聚强度降低,且孔压逐渐回升到静水压状态,导致边坡稳定性下降,最危险滑弧通过水合物分解区.若采用双水平井热激法开采,开采过程与停采后的最危险滑弧始终通过水合物分解区,由于开采过程中温度升高,井周孔压显著上升,导致边坡安全系数明显下降,存在诱发滑坡的风险.