平头锥型回转体高速入水结构强度数值分析

为探究回转体在高速入水过程中的结构强度,基于非线性有限元LS-DYNA软件中流固耦合任意拉格朗日-欧拉(arbitrary Lagrangian-Eulerian, ALE)方法,分析了不同壁厚的回转体以100 m/s的初速度入水过程中的冲击力特性和结构强度。结果表明:数值计算得到的入水冲击压强峰值和速度衰减曲线与相应的理论值吻合较好,从而验证了数值方法的有效性;入水冲击载荷峰值出现在结构入水瞬间,结构入水后冲击载荷急剧变小且微小震荡;回转体的结构形式对其在高速入水过程中的结构强度有重要影响,尤其回转体头部厚度影响回转体结构强度,当回转体头部厚度为8 mm、后体壁厚大于2.5 mm时,可以保证回转体强度要求。     

约束端面对管内甲烷爆炸特性的影响

为研究不同约束端面下甲烷的爆炸特性,利用自行搭建的实验平台完成了多种约束端面下不同浓度甲烷的爆炸实验。研究表明:约束端面的性质对甲烷的爆炸特性有显著影响,约束端面的承压强度越高,甲烷的爆炸超压越大。单层PVC薄膜作用下,薄膜破裂,不会引起火焰与超压的振荡;而纸膜破裂后,管道内外气流的高速泄放和回流则会引起超压振荡,使火焰前锋波动并发生扭曲变形;两者共同作用时,PVC薄膜会阻碍气流的泄放与回流,加速超压衰减,抑制火焰和超压的振荡。然而,随着纸膜层数增加,破膜时管道内外形成的巨大压差会使约束端面完全破裂,降低PVC薄膜的抑制作用。当破膜难度达到一定程度时,约束端面作用下的泄压峰值成为不同浓度甲烷爆炸的最大超压峰值,且泄爆压力并不随甲烷浓度的改变而改变,因此不同浓度甲烷的爆炸超压在较高的泄爆压力下相同;此时,相同约束端面下不同浓度甲烷的压力振荡曲线在压力衰减的前半个周期内完全重合,管道内外的压差成为主导超压振荡的重要因素,而不同浓度甲烷的燃烧速率对超压振荡的影响则可以忽略不计。     

近壁声空泡溃灭微射流冲击流固耦合模型及蚀坑反演分析

超声场下液体环境中近壁空泡溃灭会产生强烈的微射流,为探究微射流冲击壁面流固耦合效应,利用流体力学及冲击动力学,考虑了率相关的J-C材料本构模型,建立并分析了微射流冲击壁面流固耦合三维模型,并通过超声空化试验和基于球形压痕试验理论的反演分析进行了验证。结果表明:微射流冲击下材料表面出现微型凹坑,凹坑深度由微射流速度和微射流直径共同决定且随其增大而增大,凹坑直径主要与微射流直径正相关,而凹坑径深比则主要与微射流速度负相关;壁面压强基本呈对称分布且最大压强出现在微射流冲击边缘;超声空化试验验证了微射流冲击下材料表面出现的微型凹坑,反演分析方法表明,在16~18的径深比下,微射流冲击强度为420~500MPa,对应的微射流速度为310~370m/s。试验及反演分析结果与理论分析结果相符,验证了流固耦合模型及反演分析方法的合理性及准确性,为后续工程应用中空化强度、微射流速度等的控制提供了理论参考。     

浸没边界-简化热格子Boltzmann方法研究及其应用

针对流固耦合传热问题,本文提出了一种基于浸没边界-简化热格子玻尔兹曼方法(immersed boundary method-simplified thermal lattice Boltzmann method,IB-STLBM)的耦合模型.不同于传统的格子玻尔兹曼方法使用分布函数演化流场和温度场,简化热格子玻尔兹曼方法(simplified thermal lattice Boltzmann method,STLBM)的演化过程不需要依赖分布函数,只涉及平衡态分布函数和非平衡态分布函数,能够直接演化宏观量,极大减小了计算过程中所占用的虚拟内存,简化了边界条件的实现方式,同时具有较高的稳定性.传统的浸没边界法对流场的计算采用欧拉网格,对固体边界采用拉格朗日网格,认为固体边界是对流场产生某种体积力.在应用浸没边界法时,汲取介观的思想,把固体的介入看作是对流场的干扰,打破了固体附近流体介观微团颗粒原始的平衡状态,这种干扰可以看作是在耦合边界上产生的一个非平衡项,可用非平衡态分布函数来表示.基于此,在模型中浸没边界法与简化热格子玻尔兹曼方法更紧密联系在一起,更大程度发挥二者的优点,整个计算过程更加简单直观,符合物理特性.通过对热圆柱绕流和内含热颗粒的封闭方腔自然对流问题的模拟以及对其结果的分析,验证了该算法在求解流固耦合传热问题的有效性和可行性.      

功能梯度球形水凝胶的化学力学耦合分析

作为一种新型智能材料,水凝胶具有特殊的化学力学耦合性能.采用功能梯度形式可使得水凝胶具有更好的适应性和可调控性.本研究中假设交联密度沿径向按幂函数规律变化,并基于水凝胶的大变形多场耦合一般理论,采用Flory-Huggins自由能函数,建立了功能梯度球形水凝胶在球对称情形的控制方程,并开展了功能梯度球形水凝胶在给定内压和化学势情行的非均匀大变形溶胀行为的理论研究.计算结果表明,不同梯度指数的球形水凝胶的内压、内孔半径曲线和内压、内表面径向伸长率曲线均呈现出一段稳定区间和另一段不稳定区间,说明内压超出某临界值会发生失稳并导致水凝胶的最终破坏.内压的临界值随梯度指数的增大而增大.研究表明,功能梯度球形水凝胶的材料参数(梯度指数、亲疏水特性、交联密度和溶剂分子的体积)和环境化学势对水凝胶溶胀行为具有重要的影响.在给定内表面压力的情况下,功能梯度球形水凝胶内表面的径向位移随梯度指数的改变接近为线性变化,而随其他参数的影响都呈现出明显的非线性.本研究有助于实现水凝胶智能结构和器件在复杂条件下的精准调控.}      

基于时域解耦算法的多液舱浮式结构物运动模拟

对于带有多个晃荡液舱的浮式结构物, 浮体的运动、外场水动力以及各舱内的液体晃荡力会实时相互决定, 发生复杂的耦合作用. 为准确模拟多液舱浮式结构物的运动, 本文引入一种有效的时域解耦算法. 该方法以模态分解法为基础, 通过对浮式结构物所受外域水动力和各液舱内非线性晃荡力进行模态分解, 最终形成时域解耦运动方程, 无需迭代求解过程即可显式计算浮式结构物的瞬时加速度. 该方法可避免传统迭代求解方法在迭代次数、截断误差和收敛特性等方面的不足, 减少解耦过程的计算耗时. 本文进一步结合边界元数值方法, 分别对单液舱浮式结构物和多液舱浮式结构物的工况开展数值模拟研究. 通过与单液舱浮式结构物的实验结果对比, 验证了本文时域解耦算法的有效性. 本文详细分析了晃荡力对单液舱浮式结构物运动的影响, 发现存在一个共振影响区间: 当外场波浪频率在该区间之外时, 可以在时域计算结果中观察到稳定的浮体运动; 在比该区间更低频的波况下, 液舱晃荡力与外场波浪力相位相反甚至可以相互抵消, 此时晃荡液舱的存在可以减弱浮体运动; 在比该区间更高频的波况下, 液舱内晃荡力与外场波浪力可以具有相同相位, 此时晃荡液舱的存在会加剧浮体的运动. 本文进一步研究了四液舱浮式结构物在波浪中的纵荡、垂荡和纵摇运动情况, 发现非线性液舱晃荡可对纵荡和纵摇运动产生影响, 但对垂荡运动影响很小.      

参外联合激励下的簇发振荡及其机理

当周期激励频率远小于系统固有频率时,会存在快慢耦合效应,与单项激励不同,参外联合激励不仅会导致快子系统平衡曲线和分岔行为的复杂化,也会产生一些特殊的非线性现象,为此,本文以两耦合Hodgkin-Huxley细胞模型为例,引入周期参外联合激励,探讨在频域不同尺度耦合时该系统的簇发振荡的特点及其分岔机制．通过建立相应的快慢子系统,得到慢变参数变化下的快子系统的各种分岔模式以及相应的分岔行为,结合转换相图,揭示耦合系统随激励幅值变化时的动力学行为及其机理．研究表明,在激励幅值较小时,系统表现为概周期振荡,两频率分别近似于快子系统平衡曲线由Hopf分岔引起的两稳定极限环的振荡频率．概周期解随激励幅值的增加进入簇发振荡,导致这些簇发振荡的主要原因是在慢变参数变化的部分区间内,存在唯一稳定的平衡曲线,使得系统的轨迹逐渐趋向该平衡曲线,产生沉寂态,并随着慢变参数的变化,由分岔进入激发态．同时,快子系统中参与簇发振荡的稳定吸引子随激励幅值的变化也会不同,导致不同形式的簇发振荡．另外,与单项激励下的情形不同,联合激励时快子系统的部分稳定吸引子掩埋在其它稳定吸引子内,从而失去对簇发振荡的影响．     

振动陀螺仪单向时滞耦合系统的摄动分析

针对振动陀螺仪单向时滞耦合系统,研究了角速度、耦合强度、耦合时滞与系统输出信号之间的关系,着重关注耦合时滞对待检测角速度测量结果的影响．本文采用摄动法对振动陀螺仪单向时滞耦合系统进行理论分析,得到了系统输出信号、系统振动边界和待检测角速度的近似表达式,并将数值和近似计算结果进行了对比．理论分析结果与数值结果比较吻合,验证了运用摄动法分析振动陀螺仪单向时滞耦合系统的有效性,也表明了在运用摄动法分析系统时具有局限性．     

耦合谐振电路系统中的Fano共振现象

常见的共振响应谱表现为近似对称分布的洛伦兹线型,而Fano共振则表现出明显的非对称的响应谱.尝试在普通物理实验层次展示Fano共振现象,建立基础实验与前沿研究的联系,通过构建RLC耦合谐振电路系统,在实验上展示了Fano共振现象,作为其中的一种特殊情况,还展示了类似电磁感应透明的现象,并通过理论公式的计算解释了实验结果.     

π-π Stacking and Magnetic Coupling Mechanism on a Mono-nuclear Mn(Ⅱ) Complex

A new mono-nuclear Mn(Ⅱ) complex [Mn(MPT)2(NCS)(HOCH3)]ClO4(MPT = 2-methoxyl-1,10-phenanthroline) has been synthesized with 2-methoxyl-1,10-phenanthroline and thiocyanate anion as ligands,and its crystal structure was determined by X-ray crystallography.The crystal data:monoclinic system,space group P21/c with a = 12.8849(17),b = 15.684(2),c = 14.2703(19) ,β = 92.126(2)°,V = 2881.9(7) 3,Z = 4,C28H24ClMnN5O7S,Mr = 664.97,Dc = 1.533 g/cm3,F(000) = 1364 and μ = 0.679 mm-1.In the crystal,MPT functions as a bidentate ligand and Mn(Ⅱ) ion assumes a distorted octahedral geometry.In the crystal there are two types of π-π stacking interactions among the adjacent complexes,which involves 1,10-phenanthroline rings.The theoretical calculations reveal that the two types of π-π stacking resulted in a weak anti-ferromagnetic coupling with 2J =-14.68 cm-1 and a weak ferromagnetic coupling with 2J = 0.70 cm-1,respectively.The theoretical calculations also reveal that there are both spin polarization and spin delocalization in the two π-π stacking systems,and the magnetic coupling signs are unable to explain with McConnell I spin-polarization mechanism.