间接边界元法模拟含裂纹介质中弹性波的传播

介绍了间接边界元法的基本理论,利用二维弹性动力边界积分方程解决了地震波在含裂纹介质中的散射问题。编写了地震波在二维非均匀介质中传播的数值模拟程序,模拟了爆炸点源在含裂纹介质中的传播情况,分析了裂纹的分布及填充物对波传播的影响。结果表明,边界元法是研究含裂纹介质中弹性波传播的有效工具。     

冲击波载荷作用下固支正交各向异性薄板挠度特性分析

根据能量原理,采用层合板理论和薄板大变形理论,对四边固支的矩形正交铺层的纤维增强复合材料薄板在爆炸冲击波作用下的弹性大变形进行理论分析,得到了层合板的最大挠度计算公式。有限元数值模拟计算表明,理论计算结果与有限元数值模拟的结果吻合较好,因此该方法可用来对复合材料板在弹性范围内的最大挠度进行预测,同时认为根据薄板的最大挠度,可以计算板的应力(应变)分布和复合材料板失效时的爆炸载荷。     

基于粘弹性人工边界的点声源海底地震波ANSYS分析方法

航行舰船在海底岩土层中引起的弹性波被称为舰船地震波,该波主要由舰船低频辐射噪声引起,可用于识别舰船目标.采用大型有限元软件ANSYS中的流固耦合计算模块,结合粘弹性人工边界条件,对水平分层海洋环境下低频点声源海底地震波进行了数值计算.数值计算结果与波数积分方法结果的比较表明:两种方法相应的声传播损失曲线和频率特性曲线变化规律基本一致,海底表面竖直位移、水平位移、竖直加速度和水平加速度的频率特性曲线吻合较好,说明大型有限元软件ANSYS中的流固耦合模块可用于低频点声源海底地震波问题的数值计算.     

凝灰岩介质中爆破地震波传播的数值模拟研究

在核电站附近实施爆破开挖的一个重要问题就是控制爆破地震波的危害影响。研究爆破地震波在岩石中传播的规律和特征可为确定爆破方法和参数提供依据。本文在实测研究的基础上建立了数值计算模型,编制了有限元计算程序,对深孔爆破在凝灰岩介质中引起的地震波传播过程进行了数值模拟。计算中仅考虑弹性区域,爆破波对岩石介质的作用经粉碎区、破裂区衰减后成为输入荷载。计算结果与实测结果基本相符。说明该模型及方法可用于计算凝在     

伪谱法模拟各向异性/非均匀介质中弹性波的传播

数值模拟一直是研究弹性波在复杂介质中传播问题的有效手段之一。各向异性问题和非均匀问题日益引起关注。伪谱法是近年来发展较快、建立在弹性波动方程基础上的一种数值模拟方法。它的主要特征是将波动方程中的空间微分变换成频域中的乘法运算 ,从而只需要较少的空间格点就可以得到较高的计算精度。介绍了伪谱法及其稳定性条件 ;编制了一个适合于小型工作站或者PC机的模拟地震波在二维各向异性均匀 /非均匀介质中传播的数值模拟程序 ,模拟了各向异性 /非均匀介质中弹性波传播理论 ,分析了裂纹以及各向异性界面对波传播的影响     

炸药爆炸过程中电磁辐射数值模拟研究

为了给出炸药爆炸过程中产生电磁辐射的理论分析方法和特性,针对炸药爆炸过程中产生的电磁辐射进行数值模拟研究.采用热平衡电离模型描述爆炸过程中的电离气体,结合爆炸问题的数值模拟方法,得到爆炸场中的电导率分布,进而得到磁扩散率.基于非理想磁流体动力学方程组,计算在特定磁扩散率下的导电气体运动引起的电磁辐射,实现对炸药爆炸产生...     

圆柱壳在水下爆炸载荷作用下的抗冲击性模拟

利用LS-DYNA动力有限元程序对张效慈等人做的深水爆炸水动压力场对潜体结构的动态响应试验进行数值模拟，其模拟结果与实际具有较好的一致性；同时模拟了不同厚度圆柱壳在不同质量炸药爆炸情况下所表现的抗冲击性，为提高潜艇的抗冲击性提供参考。     

爆炸激波管管口稀疏波对试验段的影响

数值模拟了爆炸激波管不同隔离段长度时管口稀疏波对试验段超压的影响。为了准确高效地模拟试验段入口的超压曲线,采用了一种将一维球对称程序和三维程序相结合的计算方法,并在一维计算中利用爆炸相似律,采用小当量爆炸来模拟实际超压波形。计算结果表明,隔离段长度L的变化不影响超压峰值;L小于等于20m时,稀疏波的影响使得试验段超压的作用时间、比冲量减小;L大于等于30m时,稀疏波对试验段超压无影响。     

点火载荷下翼柱型装药结构完整性数值分析

针对固体火箭发动机翼柱型装药在点火载荷下结构完整性问题,提出了一种考虑三维损伤黏弹性本构性的数值模拟方法.基于商业有限元软件的二次开发平台,编写用户子程序,展开有限元数值计算,获取翼柱型装药在固化降温、点火载荷下装药结构的应力、应变与位移场分布.结果表明,在温度载荷与点火载荷下,药柱内表面变形模式是不同的.然而,不管在...     

三维超声速隔离段湍流内流场旋涡结构的数值模拟

从基于雷诺平均的N S方程出发 ,采用有限体积方法离散控制方程 ,数值模拟了三维超声速隔离段湍流内流场。计算中采用了二阶OC TVD差分格式、LU隐式算法和Baldwin Lomax代数湍流模型。数值结果与实验做了对比 ,并结合隔离段中的激波串结构分析了其横截面上旋涡结构的发展过程及不同外形条件下旋涡的不同结构。计算结果表明采用本文发展的方法模拟隔离段湍流流场是可行的 ,截面为正方形与长方形的隔离段内的涡旋结构截然不同     

舰船尾部振动分析整船混合有限元模型的集成方法研究

对舰船尾部振动分析中整船混合有限元模型的集成方法进行了研究．通过对现有的几种集成方法进行分析、比较,提出了一个较为简便、易于操作的集成方法,并用于对某型舰尾部振动的实例分析中,得出了满意的结果．本文结论可推广到采用整船混合有限元离散模型对舰船总体振动和局部振动的分析计算中．     

Semi-analytical method for calculating aeroelastic effect of profiled rod flying at high velocity

The key technique of a kinetic energy rod(KER) warhead is to control the flight attitude of rods. The rods are usually designed to different shapes. A new conceptual KER named profiled rod which has large L/D ratio is described in this paper. The elastic dynamic equations of this profiled rod flying at high velocity after detonation are set up on the basis of Euler-Bernoulli beam, and the aeroelastic deformation of profiled rod is calculated by semi-analytical method for calculating the vibration characteristics of variable cross-section beam. In addition, the aeroelastic deformation of the undeformed profiled rod and the aeroelastic deformation of deformed profiled rod which is caused by the detonation of explosive are simulated by computational fluid dynamic and finite element method(CFD/FEM), respectively. A satisfactory agreement of these two methods is obtained by the comparison of two methods. The results show that the semi-analytical method for calculating the vibration characteristics of variable cross-section beam is applied to analyze the aeroelastic deformation of profiled rod flying at high velocity.     

电磁轨道发射器动态发射过程的数值模拟

针对电磁轨道发射器动态发射过程的数值模拟问题,基于矢量磁位A和时间积分标量电位V,采用节点单元法,并选择运动坐标系描述运动问题,推导出动态发射下的电磁轨道发射器三维涡流场有限元离散方程。结合温度场控制方程,建立电磁轨道发射器的三维电磁-温度耦合有限元模型。针对一个电磁轨道发射器动态发射问题,对模型进行数值模拟,得到了动态发射下轨道及电枢的温度场分布特点及发射器电感梯度随时间的变化规律。计算结果表明,模型计算出的脉冲电流峰值、出口电流大小、出口速度等参数均与试验结果吻合较好,验证了所开发的有限元程序代码的正确性。     

随机梁的传递函数分析方法

该方法将传递函数方法与传统的随机摄动分析理论相结合 ,将随机场作Karhunen Loeve正交展开 ,对静态随机梁结构进行了分析 ,并计算了其可靠性指标。这样一来克服了传统的有限元数值方法计算量大的困难 ,确立了随机梁结构分析的一种解析方法。     

自动武器系统精度综合分析方法的研究

本文为探讨自动武器结构系统与射击精度的内在关系,以动态优化设计的思想,提出分析自动武器系统精度的一种综合方法,将武器的传统测试与理论模态分析、实验模态分析以及有限元结构动力分析有步骤、分阶段地有机结合、交叉运用,通过修改有限元结构模型达到控制武器射击精度与优化武器结构的目的。     

关于有限变形非弹性有限单元法的研究进展

对有限变形非弹性有限单元法研究发展中的一些普遍关注的研究课题,进行了回顾和展望。这些课题包括:不同参考构形下的有限元列式、客观性应力率的合理定义、有限变形非弹性本构方程的精确描述、率本构方程的增量客观数值积分、塑性流动的不可压缩性与自锁现象、以及交界面接触与摩擦的处理。     

螺旋桨激励水下艇体振动的试验及数值研究

为考虑螺旋桨真实激励特性,在循环水槽中开展艇尾伴流场中螺旋桨诱导艇尾脉动压力及螺旋桨激励水下艇体振动响应的测量试验。试验结果表明:脉动压力幅值在叶频处最大,且随螺旋桨负载增加而增大,随与螺旋桨之间距离的增大而减小,四叶桨脉动压力在尾翼后的高伴流区幅值较大,五叶桨则在尾翼之间的低伴流区幅值较大；大部分测点的振动响应幅值随螺旋桨负载增加而增大,但也存在叶频处幅值较小和未随负载增加而增大的情况；五叶桨激励引起的侧向振动较四叶桨有所增强而轴向振动有所减弱；特定谱峰频率处振动响应幅值呈一阶弯曲振型,其频率范围与有限元计算结果较为一致。综合采用计算流体动力学、有限元和模态叠加法建立螺旋桨激励水下艇体振动响应的数值计算方法,通过计算与试验结果的对比发现,该方法的计算结果与试验结果吻合较好,相比于采用单位简谐激励的谐响应分析方法更加接近真实情况。     

激光熔覆搭接过程残余应力应变的有限元计算

激光熔覆过程中产生的残余应力与应变对熔覆层的裂纹开裂倾向有重要影响。利用有限元法对激光多道搭接过程熔覆层的残余应力和应变进行计算,计算过程中主要考虑瞬时温度变化引起的热应力、σ-ε曲线、有限元网格、力学边界条件4个方面。热应力主要在温度场计算结果的基础上求得;σ-ε曲线采用的是线性强化材料的弹塑性曲线;有限元网格主要采用分区划分的形式进行。通过计算,得到了搭接后的熔覆层上关键点的残余应力和应变分别为1 400 MPa和2.2×10^-2。结果表明：搭接熔覆层交界处的点产生焊趾裂纹的倾向最大,第1道熔覆顶点处产生横向裂纹的倾向次之。     

细长壳体水声辐射问题的有限元结合边界元解法

采用三角形面积元,提出了一种计算任意形状的细长壳体水声辐射问题的有限元结合边界元解法,推导了三角形面积坐标下的质量矩阵。在水动力方面,根据单层势的概念,提出了声场速度势的边界积分表示式,给出了三角形面积坐标下的附加质量矩阵和附加阻尼矩阵,从而得以建立流固耦合的有限元动力响应方程。本方法的进一步发展可用于计算具有加强肋骨的双层壳体结构的声辐射问题以及船体在随机激励下的远场声功率谱。