任意平面区域的高质量有限元网格自动剖分

实现了任意平面区域的高质量三角形有限元网格自动剖分。初始化时利用广度优先搜索查找孔及凹槽中需删除的三角形 ,不需定义有向边界 ,可十分方便地定出实际剖分区域 ;给出了一种方便、快捷的点定位法 ,可大大提高程序运行速度 ;利用Delaunay网格优化算法实现了对三角形形状的控制 ,可使三角形最小角达到 35° ;分割坏三角形时 ,按最小角由小到大的顺序处理 ,使三角形网格在满足质量要求的同时网格数量较少     

非结构网格质量对梯度重构及无粘流动模拟精度的影响

综合利用理论分析和数值测试手段,研究了非结构格心型有限体积离散中梯度重构算法的计算精度,分别给出了非结构算法中常用的基于Green-Gauss公式(GG方法)和基于Least squares方法(LSQ方法)的梯度重构方法达到至少一阶精度的条件。其中,GG方法在面积分低阶项不能互相抵消的情况下,要求面心插值精度达到至少二阶;而LSQ方法对于任意网格均能实现梯度重构一阶精度。在各向同性网格上的梯度重构精度数值测试验证了数学推导结论;进一步通过制造解方法量化无粘流动数值离散误差,结合网格收敛性测试研究了网格质量(网格点随机扰动、网格弯曲度和网格倾斜度等因素)以及网格类型(三角形和四边形)对无粘流动模拟精度的影响,验证了理论分析结论。     

一种自适应轴对称FEM-SPH耦合算法及其在高速冲击模拟中的应用

为了准确、高效地模拟高速冲击问题,提出了一种自适应轴对称有限元(FEM)-光滑粒子流体动力 学(SPH)耦合算法。该算法在初始时刻全部采用FEM 计算,在动态变形过程中自动将畸变单元转化为粒 子,采用SPH 计算。该算法采用一种新的耦合算法实现单元与粒子间的高精度耦合,并应用最小内角转化准 则和单元分组转化方式实现单元向粒子的自动转化。计算了几个典型的高速冲击问题:首先,通过计算应力 波传播测试了新的单元-粒子耦合算法的精度;然后,通过计算泰勒杆问题验证了自适应耦合算法及相应程序 的正确性;最后,计算了弹体侵彻铝板和混凝土板。结果表明:自适应耦合算法计算精度好且效率高,适合模 拟高速冲击问题。     

高速冲击过程数值分析的再生核质点法

利用新型的无网格方法-再生核质点方法对高速冲击过程进行数值模拟，给出了控制 方程，分析中引入Bordner-Partom本构模型来实现材料高速冲击条件下的大应变和高应变 率的特性，提出了新的接触面处理方法，介绍了速度配点法来实现接触面及其它本质边界的 速度条件. 数值模拟表明，再生核质点法能够准确模拟高速冲击过程，具有较高的计算精度.     

冲击载荷作用下结构动态响应分析的罚函数有限元法

基于罚函数法，根据接触－冲击问题的边界条件，运用双线性参数法推导出用有限元模拟三维接触－冲击问题的具体算法，并对杆和梁在冲击载荷下的动态响应进行了仿真分析，取得较好计算结果。     

基于单元相交的混合网格精确守恒插值方法

基于网格切割思想,发展了二维/三维混合网格条件下的单元相交算法,可精确计算任意两个多边形/多面体的交集。在此基础上,实现了基于单元相交(CIB/DC)的精确守恒插值算法。二维和三维验证算例表明,该方法能够保证插值过程中计算域内物理量的严格守恒,且具有比常规二阶插值更高的精度。     

用逐点插入法自动生成全四边形的自适应有限元网格

本文给出一种有限元网格全自动划分的方法。方法包括三部分：１）用单元尺寸场控制生成网格的疏密分布；２）用基于Ｄｅｌａｕｎａｙ三角化过程的逐点插入法生成纯三角形网格；３）将三角形网格转化为全四边形网格。相应的网格生成器具有良好的用户界面。只须最小限度的边界描述即可自由、快速地生成单元形状良好、疏密分布任意的网格。     

用逐点插入法自动生成全四边形的自适应有限元网络

本文给出一种有限元网格全自动划分的方法。方法包括三部分；１．用单元尺寸控制生成网格的疏密分布１；２．用基于Ｄｅｌａｕｎａｙ三角化过程的逐点插入法生成纯三角形网格；３．将三角形网格转化为全四边形网格。相应的网格生成器具成良好的用户界面。只须最小限度的边界描述即可自由，快速地生成单元形状良好，疏密分布任意的网格。     

高导无氧铜杆的冲击拉伸断裂系列实验以及基于样本体积单元的分析

采用Hopkinson装置和一种基于一级气体炮的高速冲击拉伸断裂装置,研究了无刻槽高导无氧铜 (OFHC)杆在一系列冲击拉伸速度下的断裂。当冲击拉伸速度大于40m/s时,断裂位置总在冲击拉伸端附 近,此速度被确定为OFHC的实验临界冲击拉伸速度。一种受单轴冲击拉伸荷载的、中心含椭球空穴的样本 体积单元被用于数值模拟所含空穴的增长与失稳的过程。OFHC的J-C与Z-A 本构关系用于描述基体材料 的动态响应。讨论了空穴失稳条件并提出以空穴形状演化为判据,比较了空穴失稳时的样本体积单元平均径 向应变与无刻槽杆的冲击断裂应变。也用这种样本体积单元模型分析了OFHC的实验临界冲击拉伸速度。     

任意多连通域网格自动生成算法及其应用

给出了一种任意平面域内三角形网格自动生成的算法。本算法是根据下面两个事实而设计出来的：１．任意平面多连通域总可以剖分成若干个三角形单元。２．平面域内的活节点数目随着单元剖分总会趋近于零。文中还介绍了几种网格性能改善的方法。整个程序结构紧凑，需要人工输入的数据量少，生成的网格质量较高，而且算法稳定性能好。     

隧道爆破近区振动的预测方法

将比例距离的概念应用于隧道爆破振动分区中,以速度衰减曲线斜率的大小作为分区计算的依据。在此基础上,提出采用BP小波神经网络的方法预测爆破近区振速,以棋盘山隧道实测数据验证模型的可行性。结合泉厦高速公路大坪山隧道工程实例,对临近既有隧道形式的隧道爆破地震波传播规律进行分析,并对近区振速进行预测。分析表明:（1）在无实测值时,隧道爆破分区可按比例距离大致划分为:比例距离 < 5.0为爆破近区;5.0≤比例距离≤9.0为爆破中区;比例距离 > 9.0为爆破远区。（2）BP小波神经网络爆破近区预测模型不仅适用于新建分离式隧道,也适用于临近既有隧道的新建小净距隧道。研究成果对复杂环境下的隧道钻爆施工具有一定的指导意义。     

惰性介质对钝感炸药爆轰的约束效应

首先用改进冲击波极曲线理论分析惰性介质对钝感炸药爆轰的约束作用类型。改进冲击波极曲线基于爆轰ZND模型建立在前导冲击波上,并且未反应炸药采用JWL状态方程,惰性介质采用p（ρ,T）形式状态方程。理论考察声速小于炸药CJ爆速且压缩性不同的6种典型惰性介质约束情况。然后用带三项式Lee-Tarver化学反应率的二维Lagrange流体力学方法数值模拟考察约束相互作用。数值考察约束介质的影响因素有:压缩性、厚度、典型双层介质组合约束。从数值结果看出,由介质压缩性的不同给出的约束作用方式共7种:其中6种出现在介质声速小于炸药CJ爆速条件下,可运用冲击波极曲线理论;另外一种出现在介质声速大于炸药CJ爆速条件下,不能使用冲击波极曲线理论。同时,介质厚度、双层介质组合方式也能够影响爆轰前导冲击波阵面形状以及爆轰化学反应流动状态。     

高流态地质聚合物混凝土的高应变率动态压缩变形特性

以矿渣、粉煤灰为原材料,以NaOH、Na₂CO₃为碱激发剂,制备了强度等级为C30的高流态地质聚合物混凝土（highly-fluidized geopolymer concrete,HFGC）,运用波形整形技术改进了?100mm SHPB实验装置,通过参数控制保证应力均匀和恒应变率加载,对HFGC开展了动态压缩实验,分析了HFGC在冲击压缩荷载下的变形特性。HFGC属于应变率敏感材料和脆性材料,高应变率作用下,HFGC的典型应力应变曲线包括压实挤密阶段、弹性阶段和软化、屈服阶段。在10~100s-1的应变率范围内,HFGC的峰值应变ε_c随应变率的变化表现出明显的冲击韧化效应,ε_c随应变率的升高先增大后减小,满足二次函数关系ε_c＝-1.2×10-6+1.6×10-4+0.001 7,变形特性变化的临界应变率为66.7s-1。HFGC的动态弹性模量均低于其在准静态下的弹性模量。     

弹体高过载硬回收测量技术的实验探讨

为提供引信服役的力学参数,达到引战配合目的,弹体过载(减速度)的测量十分必要。采用硬回 收法测量了高速侵彻混凝土的先进钻地弹缩比弹过载时间历程。实验弹最大直径48mm,计及过载测量装 置,弹体质量约1.8kg,撞击速度1142m/s。过载测试系统模块化和集成化设计,包括过载测量模块、电路模 块、信号存储模块、电源模块和外部计算机模块,传感器为压阻式。量程在(-8~21)104g,可测量过载幅值 约20104g。     

用CSTBD试样确定砂岩的动态起裂和扩展韧度

利用大直径霍普金森压杆径向冲击中心直裂纹巴西圆盘(CSTBD)砂岩试样,完成Ⅰ型动态断裂实验。利用实验-数值方法确定了不同动态加载率下砂岩的动态起裂韧度;结合实验-数值法以及普适函数确定了不同裂纹扩展速度下砂岩的动态扩展韧度。为验证普适函数法和实验-数值法的有效性,将实验所得结果与其他学者的研究成果进行了对比分析,得到了相同的规律。所确定的岩石动态起裂韧度和动态扩展韧度分别随动态加载率的提高和裂纹扩展速度的提高而增加。     

碱类型对地聚合物混凝土应变率效应影响的对比研究

采用经波形整形技术改进后的?100 mm SHPB实验装置开展了强度等级均为C30的NS激发矿渣粉煤灰基地聚合物混凝土(NaOH and sodium silicate-activated slag and fly ash based geopolymer concrete,NSSFGC)和NN激发矿渣粉煤灰基地聚合物混凝土(NaOH and Na₂CO₃-activated slag and fly ash based geopolymer concrete,NNSFGC)的动态压缩实验,并对比分析了在冲击荷载作用下的应变率效应。结果表明:NSSFGC和NNSFGC的峰值应力均随应变率的增加而越大。这表明,GC为应变率敏感材料;在动态压缩状态下,NSSFGC和NNSFGC的应变率敏感阈值分别为51.82、28.89 s-1;GC的应变率敏感性强于普通混凝土;NN型激发剂更有利于发挥GC的整体强度特性。由此可知,NNSFGC的应变率敏感性明显强于NSSFGC。     

RDX/TNT、HMX/TNT炸药Dn(κ)关系研究

研究了不同配比HMX/TNT炸药的Dn()关系,并与不同配比RDX/TNT 炸药的Dn()关系进 行了对比。结果表明,配比相同时HMX/TNT炸药Dn()关系曲线总是在RDX/TNT炸药Dn()关系曲线 上方;RDX/TNT和HMX/TNT炸药的Dn()关系曲线随着炸药中TNT含量增加而下移。引入变量Rt 分 别将两种混合炸药的Dn()关系推广为Dn(,Rt)关系,给定变量Rt 就可以通过Dn(,Rt)关系得出相应未标 定的Dn()关系,Dn(,Rt)关系不仅可以减少实验标定次数,并且能满足DSD方法计算精度要求。     

前驱体相对物质的量的变化对氢氧气相爆燃制备纳米SiO2颗粒的影响

以SiCl4为前驱体,以H2 和O2 的混合气体为爆源,通过改变前驱体的相对物质的量,考察前驱体 的相对物质的量对产物的影响。并通过XRF、TEM 和XRD等测量手段对纳米SiO2粉体进行测试分析:通过 XRF对产物进行成分分析,测得所得产物主要由SiO2组成,且纯度达到99.9%以上;由XRD结果得出,生成 的SiO2产物属于典型的不定形结构;由TEM 图片可以得出,生成的产物粒径大小大部分在50~100nm 之 间,且呈圆球状。随着前驱体相对物质的量的增加,所得到的产物粒子逐渐出现粘结,分散性下降,前驱体相 对物质的量应该控制在1.25以内。     

复合射孔爆燃气体压裂裂缝起裂扩展研究

结合线弹性断裂力学的裂缝尖端应力强度因子判据,建立了复合射孔爆燃气体压裂裂缝的起裂扩展模型,通过建立与多个变量相关的缝内气体压力分布函数,利用迭代法实现了模型的数值求解,获得了缝内气体压力分布随时间的动态变化规律,并分析了不同特征参数对裂缝起裂扩展与止裂过程的影响。实例计算结果表明:(1)随着裂缝扩展的进行,爆燃气体流动尖端与裂缝尖端经历了由重合到不重合再到重合的过程;(2)地应力越大,裂缝起裂扩展越困难,爆燃气体有效致裂作用时间越短,最终得到的裂缝扩展长度也越小;(3)初始裂缝越长,裂缝更容易起裂扩展,爆燃气体能量利用率越高,裂缝扩展更长;(4)岩石断裂韧性的改变对裂缝起裂、止裂和裂缝扩展长度没有明显的影响;(5)升压速率越小,爆燃气体有效致裂作用时间越长,最终裂缝扩展也更长,但对裂缝起裂压力与止裂压力几乎没有影响。     

超高性能水泥基复合材料抗多次侵彻性能研究

采用普通混凝土成形工艺和热水养护制度,制备出了抗压强度120~250 MPa的超高性能水泥基 复合材料。在获得优异静态力学性能的基础上,进行了超高性能水泥基复合材料抗多次侵彻研究,分析了侵 彻过程中的靶体破坏形态。根据获得的侵彻深度数据,在J.T.Gomez等提出的多次侵彻模型基础上,修正了 与靶体抗压强度有关的多次侵彻深度计算经验方程。