舰艇抗冲击试验中爆距定位理论分析与计算

在舰艇抗冲击试验中,爆距直接影响到对目标的毁伤效果,爆距的误差范围也直接影响到试验数据的处理,因此,依据爆距的要求,需要对爆源布放沉底后进行准确的定位.准确计算出爆源布放至海底的位置尤为重要.本文对试验中爆源布放后布放索及浮漂在水中的受力情况进行了分析,建立了爆源布放索-浮漂系统受力的数学模型,采用Runge-Kutta方法和Newton迭代法计算出浮漂的漂移量,并给出了具体的算例.     

水下爆炸对沉箱重力式码头毁伤效应

通过开展不同爆距下水下爆炸对沉箱重力式码头模型毁伤效应试验,对水下荷载进行了采集分析,对毁伤因素、毁伤模式和毁伤机理开展了研究,初步讨论了爆距的影响。结果表明:试验中未形成完整的气泡脉动过程,荷载超压主要出现在冲击波传播阶段,爆炸冲击波、水底反射波和侧壁反射波是主要的毁伤因素;水下爆炸对沉箱重力式码头造成的毁伤程度大、模式多、机理复杂,主要毁伤部位为迎爆面外墙、迎爆侧管沟、封仓板和面板;爆距越近码头毁伤越严重;当爆距过近时,爆炸能量被迎爆面结构变形大量吸收,迎爆面毁伤程度大幅增大,码头其他部位毁伤程度增幅放缓。     

基于爆能等效原理大型模爆器燃气爆炸冲击加载的数值模拟

运用非线性显式动力有限元程序LS-DYNA,基于多物质Euler算法,对TNT炸药和乙炔-空气混合气体两种爆炸源在自由大气场中爆炸产生的冲击波荷载特征参数进行数值模拟,比较两种爆源产生的冲击波压力传播规律。基于爆能等效原理,按超压相等的原则给出了气体爆炸名义比例距离计算公式。结果表明,基于Euler算法可以较好地描述乙炔-空气混合气体爆炸空气冲击波传播规律,爆炸压力随着距爆源距离的增大而迅速衰减,且两种爆源产生的冲击波超压峰值误差随着冲击波传播距离的增大而逐渐减小。采用名义比例距离公式修正后,气体爆炸与炸药爆炸冲击波计算误差可以得到有效控制。当爆炸冲击波超压小于0.5MPa时,可以采用乙炔-空气混合气体代替化学炸药进行模爆器内爆炸实验加载。     

三峡工程基岩爆破振动特性的试验研究

以三峡工程坝基岩体开挖爆破为背景 ,在弱风化花岗岩底板内进行了 6次现场爆破试验 ,测量出距爆源不等距离处的 36组地震波形。通过频域与统计分析发现 ,岩石质点振动主频率与药量、距离成反比关系 ;地震波作用时间与爆破药量成正比关系 ,而与距离成反比关系。应用神经网络理论建立的基于爆破地震效应先验知识的网络模型 (PKFN)能很好地描述爆破地震波的衰减规律 ,其计算平均相对误差仅为 3 .5 %。用地震层析成象方法 (CT)较准确地测定出了岩体爆破松裂区边界 ,并结合PKFN模型得到了三峡工程坝基岩体的临界质点振动速度范围 13 .816 .6cm/s。     

水下爆炸对重力坝的毁伤效应及最优爆距

为研究不同爆距水下爆炸对重力坝的毁伤效应,并探讨是否存在“最优爆距”,基于离心模型试验建立了炸药-库水-空气-重力坝结构的全耦合数值模型,并设计了60组数值计算工况。不同工况水深均为600 mm,炸药量为2.2 g,重力坝模型几何比尺为1/80,包含5组爆深（50～250 mm）,每组爆深对应12组爆距,爆距范围为10～200 mm,相应比例爆距范围为0.077～1.54 m/kg1/3。对比分析了不同爆距水下爆炸对重力坝的毁伤程度,并定量比较了重力坝平均损伤、单元删除率、应力、应变等参数。结果表明,对于重力坝整体结构破坏,如重力坝整体弯曲导致的拉伸破坏,水下爆炸对重力坝的毁伤效应存在“最优爆距”,即随着爆距增加重力坝毁伤程度先增加后降低;与之类似,随着爆距的增加,重力坝上游坝面损伤区域的平均损伤、重力坝单元删除率、坝踵最大拉应力平均值和坝踵最大拉应变平均值先增加后降低且在40 mm爆距附近达到最大值。保持水深、炸药量和重力坝几何模型相同,5组不同爆深近水面水下爆炸对重力坝毁伤效应的“最优爆距”均在40 mm附近,表明近水面水下爆炸时爆深对“最优爆距”不存在显著影响。     

圆筒中炸药装药的爆速估算方法探讨

在圆筒试验中利用光光源药柱激发氩气发光所用时间与圆筒装药中爆轰波传到摄像位的时间之差来估算圆筒中炸药的爆速 ,并计算了两发TNT炸药圆筒中的爆速。结果表明 ,该方法准确程度较高 ,而且测试简便 ,能够应用于实际。     

水下爆炸冲击波载荷作用下冰层破碎特性及其影响因素

水下爆炸破冰是复杂的爆破工程,为了研究冰层在水下爆炸冲击波载荷作用下的破碎特性及规律,利用几何动力分析软件LS-DYNA对水下爆炸破冰的过程进行数值模拟,并将计算结果与实验结果进行对比,误差在8%以内,验证了数值模型的有效性。根据本文中的建模方法及建立的模型,计算不同的实验工况:实验场地环境不变,调整爆距分析不同爆距下冰层破碎特性;调整药量、爆距和冰厚,通过正交设计方法设计9组实验方案,应用灰色系统理论对3种因素进行分析,建立了各个因素与破冰半径之间的灰色关联度系数及灰色关联度。分析结果表明:药量为100 g,冰厚为29 cm,水深为2.9 m,爆距范围为0.3~1.5 m破冰的半径范围为0~1.1 m,最佳爆距范围在0.3~0.45 m之间;根据以上9种工况的分析可知,药量(100、200、300 g)、爆距(0.3、0.6、0.9 m)和冰厚(24、28、32 cm)对破冰半径的影响的主次关系依次为爆距、药量、冰厚。     

水下爆炸作用下基于声学的爆源子结构输入方法

无限域吸收边界和爆炸波输入是水下爆炸载荷数值模拟的两个关键问题.本文借鉴基于内部子结构的地震波动输入方法和爆源子结构多尺度分析方法,考虑水下爆炸载荷与地震载荷同属于波动问题,提出一种水下爆炸作用下爆源子结构的爆炸波输入方法,该方法首先通过对爆源区域进行水下爆炸波自由场的波场分解,并利用地震波的内部子结构输入方法,将该自由波场运动转换为等效爆炸波输入载荷,从而实现了水下爆炸问题中冲击波的输入;本文采用圆形爆源子结构区域,并采用AUTODYN软件中一维自由场模型计算水下爆炸作用下该区域的自由波场压力时程.进一步,基于连分式近似方法提出了一种模拟无限域水体辐射效应的高精度时域人工边界条件,该吸收边界可设置于离结构和爆源子结构较近处.本文提出的分析方法通过圆形爆源子结构进行水下爆炸载荷的转化,并采用高精度吸收边界大大减少计算区域,既保证了计算精度,又降低了单元数量,具有较高的计算效率和较强的实用性.通过数值算例分析,验证了本文模型与方法的准确性,模拟了水下爆炸作用下的冲击波和气泡脉动阶段测点的压力时程曲线,并研究了圆形结构对水下近场爆炸波散射效应的影响规律.     

超高压撞击流技术下制备破壁灵芝孢子的机理研究

在本研究的前期工作中,超高压超临界撞击流方法制备破壁灵芝孢子取得成功.为了进一步提高破壁率,本文研究了灵芝孢子在超高压撞击流下的破壁机理.从灵芝孢子粉破壁的过程出发,结合灵芝孢子粉的结构,分析了孢子在原料罐中及刚出喷嘴时孢壁的受力分析,超高压超临界流体气爆时孢壁的受力分析,以及孢子撞击靶板时的受力分析.分析结果表明:孢子在原料罐中承受超高压,产生开裂的起裂源.经过卸压气爆和超高速撞击,使得灵芝孢子破壁.研究结论可为工艺参数的优化提供重要依据.     

水下爆炸时爆源与压力探头间距离的测定

(一)前言水下爆炸的各主要参量如冲击波压力、冲量、能流密度等都是距爆源距离的函数,即这些参数的值都是指某一特定距离下的值.如果用传感器去测定某点的爆压等参量,则必须同时确切地测出该点与爆源的距离,否则所测得的 ...     

THE OPTIMAL SOLUTION FOR THE TWO-DIMENSIONAL CONVECTION-DIFFUSION MODEL IN BATCHING PRODUCT PIPELINE

将一种简化的自适应网格应用于顺序输送二维混油模型的计算中,网格的疏密程度由速度和浓度梯度控制且用动网格取代传统的固定网格。混油模型中采用混掺长度理论计算速度场和扩散系数,并考虑了油品黏度变化对混油的影响。提高传统数值解法的精度,减少计算量,使之适应于长输管线顺序输送的混油数值计算。通过在计算机上的数值模拟实验,证明了该方法的有效性,分析了雷诺数、输送顺序对混油的影响并得到混油随输送距离的变化规律。     

NONLINEAR LARGE DEFORMATION BENDING PROPERTIES OF THE TRACHEA1)

研究了人工气管和生物气管拉伸和弯曲的力学性能及仿真方法。根据气管拉伸实验数据,拟合建立了人工气管的线弹性和生物气管的非线性幂函数本构模型,其与实验的最大误差为18.8%。根据非线性弯曲理论,建立了气管大变形的弯曲变形方程,并得到了其解析解,计算的弯曲载荷与实验的最大误差为8.7%。基于显式有限元法提出了气管弯曲变形的有限元仿真方法,犬气管的弯曲角度的有限元计算结果与实验的误差为3%。     

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

STUDY ON THE PARAMETER VALUE DOMAIN OF MACRO-MICRO CONSTITUTIVE MODEL

宏微观耦合本构模型的参数识别往往通过反分析方法进行,为了使参数识别结果具有高的置信度,需要确定合适的参数取值范围. 基于动态再结晶过程的微观机理以及相应本构方程的数学特征,提出一个确定参数取值范围的方法. 首先详细给出考虑动态再结晶的黏塑性本构模型,并根据模型构造物理机理,提出通过6 步确定该模型参数取值范围的方法;其次,对300M 低合金钢进行不同温度、应变速率下的热变形试验,测试宏观的流动应力-应变数据及微观的组织数据;然后应用提出的方法,依据试验数据,确定参数取值范围;最后,基于确定参数取值范围中获得的知识,对模型进行局部修改,使模型模拟结果更接近实验结果.     

INVESTIGATION OF THE RELATION OF IMPORTANCE ANALYSIS INDICES FOR MODEL WITH CORRELATED INPUTS

为了根据不同需求找到输入变量相关情况下最优重要性测度分析方法,非常有必要对现有重要性测度方法之间的关系进行探究比较. 分别以不含交叉项和包含交叉项的二次多项式输出模型为例,解析推导了相关正态输入变量情况下基于协方差分解的重要性测度指标,包括总贡献、结构贡献和相关贡献. 进一步以所推导的解析结论为基础,理论推导了传统基于方差的重要性测度指标与基于协方差分解的重要性测度指标之间的关系.并从特殊二次多项式模型的研究结果对一般模型做出推断,然后从高维模型分解的角度验证所推断的结论,并详细阐述了不同重要性测度指标的优缺点. 最后结合具体算例深入地分析了各指标之间的关系,为输入变量相关情况下结构性能对输入特性的重要性分析与工程设计提供指导.     

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性-塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.      

饱和砂土中泥浆渗透的变形-渗流-扩散耦合计算模型

传统的泥浆渗透计算中没有考虑土体变形和浆液流速的影响.根据泥浆颗粒的质量守恒定律推导了耦合流速的浓度扩散方程,并通过在浓度方程中引入沉积系数进一步计算得到沉积颗粒的质量;同时,以沉积量作为耦合项对毕奥固结方程中的水量连续方程进行了修正,在此基础上建立了变形-渗流-扩散耦合的控制方程及其变分原理. 采用有限单元法求解基本方程,运用了时间增量法与直接迭代法,并利用一维试验验证计算方法的可靠性,并与赫齐格的经典模型的计算结果进行了比较,结果表明,本文建立的模型的计算结果可以较好地预测各组试验中颗粒的沉积规律,且吻合程度优于仅考虑颗粒对流和扩散的传统计算方法. 最后,将泥浆在槽壁中的渗透简化为二维问题并进行了计算,计算结果与工程认识相符合,泥浆的沉积填充效应随深度的增加而增大,施工时需要严格控制浅层作业段的机械垂直度;成槽机的下斗抓挖时机可以根据地层填充的致密程度进行计算,对现场施工具有一定的指导意义.      

豚鼠内耳前庭-半规管生物力学模型研究

人体维持平衡主要依靠前庭、视觉及本体感觉3个系统的相互协调来完成,其中前庭系统最为重要.因其结构细微且位置深在,传统方法难以满足现代前庭医学定位、定性、定量的研究需要,而生物数值模拟研究方法在现代耳科学研究中优势显著.建立精准的生物数值模型是其重点之一.本研究通过连续组织切片技术获取豚鼠内耳膜性结构的二维解剖数据,进一步建立同时包括前庭和3个半规管的宏观三维生物数值模型,其空间结构特征及尺寸与解剖学观察一致;数值模拟临床医学冷热实验,量化描述了不同环境温度激励下半规管内壶腹嵴顶位移、速度和压强等参数随时间变化的生物力学特征,其与临床观察到的眼震结果一致.总之,通过连续组织切片技术获取内耳膜性结构二维解剖数据并据此建立内耳三维生物数值模型的研究方法可行,所建立的生物模型可满足前庭-半规管平衡功能定位、定性、定量的研究需要.      

飞机机翼损伤的气动模型及故障诊断研究简

为了确保机翼损伤后飞机的飞行安全，提出了一种在线的故障诊断方法. 首先，根据输入输出特性，采用遗忘因子递推最小二乘法对飞机故障后的气动导数进行辨识，建立了机翼损伤故障的数学模型；然后，结合多模型方法和中心差分卡尔曼滤波器（central difference Kalman filter，CDKF）各自的优点，实现对机翼损伤的故障诊断，并采用强跟踪滤波器在线更新CDKF 的采样点，以增强CDKF 的自适应能力. 最后，通过仿真结果验证了本文所提方法的有效性.     

新型涡轮流量计的理论模型简

基于动量力矩定理和流体边界层理论,建立了一新型涡轮流量计的理论模型. 利用该理论模型可以分析叶轮几何参数对流量计计量性能的影响. 以一DN50 涡轮流量计为例,利用本文所建立的理论模型对该流量计的仪表系数进行了计算,并在一体积管标准装置对该流量计的仪表系数进行标定. 结果表明,计算值和实验值较为吻合,误差在±;3.5% 以内,验证了模型的有效性,从而为仪表的结构优化设计提供了理论依据.