简介：计算机病毒的存在,使得很多计算机无法正常运行,并造成巨大的损失.针对这种情况,提出一类新型计算机病毒的最优控制.控制首先对原有的计算机病毒模型进行改进,加入控制项,提出最优控制问题,并证明最优控制的存在性,然后利用庞德里亚金的极小值原理进行理论分析,最后进行数值模拟.数值模拟结果表明,在没有控制的条件下,原模型得出基本再生数大于1,说明存在地方病平衡点,并且最终会导致病毒爆发,而数值模拟揭示,运用有效的控制策略能够更好地抑制计算机病毒的传播.

简介：利用计算流体力学软件Ansys/FlotranCFD,首先对粘性不可压缩流体的固定圆管绕流进行了数值模拟,然后结合逐步积分法完成了同时考虑纵横两向弹性支撑圆管绕流旋涡脱落诱导振动的数值模拟,并通过快速傅立叶变换,得到了弹性支撑圆管和固定圆管的升力及弹性支承圆管横向位移响应的功率谱.通过计算结果分析,得出了一些有价值的结论,可供从事具有圆管绕流构件设备设计的工程技术人员参考.

简介：气体炮以其优良的性能,在兵器弹道环境模拟领域得到了较好的应用.首先分析了典型弹药的发射环境参数,并找出了其特征值.然后运用气体动力学的相似理论建立了气体炮的内弹道模型,并在计算机上进行了数值模拟.为后续气体炮的结构设计、相关设备的选型、气体炮参数的调整提供了理论依据.

简介：针对日益受到关注的液体晃动问题，提出了一种基于浅水波理论的研究方案．该方案采用浅水波理论而非势流理论导出系统控制方程，并通过哈密顿体系表达；利用中心有限差分法和Stormer-Yerlet算法进行空间和时间离散；模拟了不同初值条件下的液体晃动情况并对比分析了影响系统非线性响应的主要因素．结果表明，基于浅水波理论能有效解决液体晃动问题；与Euler格式对比，Stormer-Verlet算法精度较高；除共振外对于系统非线性响应的影响容器初始位移比初始速度更显著；非共振情况一定条件下，充液容器运动过程中液体晃动能起到阻尼作用．

简介：基于拉格朗日描述的柔性多体系统动力学理论,采用绝对节点坐标有限元方法描述液体大变形运动,开展铁路液罐车内液体晃动模拟研究.本方法能够模拟液体自由表面的连续性变化,并适用于研究具有复杂外形容器的内部液体晃动问题.基于流体力学牛顿体基础理论,推导液体粘性方程和满足体积不可压缩的条件方程;采用基于绝对节点坐标方法描述的实体单元进行液体网格划分;采用罚函数方法描述液体与罐体之间的接触关系,组建液体-罐体耦合多体系统动力学方程.仿真计算液罐车内液体的横向和纵向晃动行为,发现液体自由表面形状呈非线性变化,不同断面处的高度和形状不同.

简介：为了研究零质量射流的作用机理和流场结构，发展了一套面向二维零质量射流的非结构化动网格模拟方法：采用控制容积法，引入动网格控制方程，并与任意曲线坐标系下矩阵形式的时均可压缩N-S方程组联合求解，迭代过程中采用弹性类推法进行动态网格更新．基于此方法，对二维零质量射流进行数值模拟，对计算获得的流场涡线和流线分布进行了分析和讨论，并与其他学者类似算例进行了比较，表明该方法能够合理揭示零质量射流的流场结构和作用机理，可实现二维零质量射流的数值模拟．

简介：根据Hamilton原理和假设模态离散化，建立了柔性梁绕定点旋转过程的动力学模型，并通过假设模态法对方程进行了离散处理．在此基础上，基于中心差分法的计算原理，提出了柔性梁动力学方程的子循环计算公式，分别建立了其同步更新格式和子步更新格式，在子循环积分过程中，通过步长修正保证了计算的精度和稳定性．计算结果表明：所提出的算法能在保持合适的精度要求下，有效地提高响应的计算效率，并且通过积分步长的修正可以提高计算稳定性，有效地处理了方程的刚性问题．

简介：基于日本东京工艺大学风洞试验数据,针对低矮四坡房屋展开研究,验证了RNG模型对四坡低矮建筑表面风压研究的可靠性.利用此模型深度研究檐口外伸长度与出挑高度对四坡低矮屋面表面风压的影响规律.结果表明：檐口的长度及出挑高度对屋面风压影响明显,当檐口外伸长度不变,随檐口竖向高度增加,迎风面风吸力随之增大,当檐口竖向高度不变,高度为0.5m、1.0m时,随着檐口外伸长度的增加,迎风面负风压减小,檐口外伸长度为1.5m,竖向高度为1.0m为最有利于房屋表面抗风设计,结论可为台风多发地区低矮民居设计提供建议.

简介：利用加性掩盖和函数调制两种混沌加密方式对模拟信号进行加密，分别从幅值和频率两方面分析加性掩盖方式和函数调制方式，对比两种加密方式加密效果，了解两种加密方式的差异．计算结果表明：函数调制方式在幅值和频率的范围上都好于加性掩盖方式的幅值和频率范围，函数调制方式比加性掩盖方式更具安全性．

简介：将振型向量与向量空间余弦因子的概念相结合,提出了度量模型整体相关度的定量指标.即计算模型和试验模型的整体相关度指标VM,以此来修改计算模型局部参数的方法.算例表明,相对于传统的方法而言,该方法能更好的表达出计算模型和试验模型之间的相关关系,并且能有效的用于工程实际.

简介：在一类高维映射中实现了由Iooss等人提出的映射不变圈的算法．首先分析了不变圈的分岔条件，然后通过Fredholm择一方法分析了在计算不变圈过程中出现的一类方程解的存在性，再根据不变圈上映射到自身的不变性，通过分析振幅各阶项的系数，最终在一高维映射中实现了不变圈的计算。

简介：研究了重物对圆板的冲击问题.采用伽辽金原理及拉普拉斯变换推导出了物体对圆板的冲击力解析表达式.通过数值实例,讨论了圆板半径、板厚、缓冲垫刚度、重物下落高度、重物质量等因素对重物对圆板冲击力影响,并绘出了冲击力随时间的变化曲线.算例表明:用该法求冲击力问题,不但比传统的Hertz接触理论更接近真实情况,而且计算简便,便于工程设计人员应用.

简介：在Birkhoff框架下，采用离散变分方法研究了非Hamilton系统一Hojman—Urrutia方程的数值解法，并通过和传统的Runge—Kutta方法进行比较，说明了在Birkhoff框架下研究这类不具有简单辛结构的非Hamilton系统可以得到更可靠和精确的数值结果．

简介：为了保证空间站正常运行所需要的大量能源及发射入轨时的较小重量代价,我国空间站采用大面积柔性电池翼的设计方案.大柔性电池翼因其结构复杂且频率极低,在轨运行时又要承受多次大机动过程产生的载荷,各承力组件的强度校核和参数设计成为影响整个空间站任务成败的重要环节.本文从工程应用角度出发建立了空间站柔性电池翼的线性简化模型,计算并分析了柔性电池翼在轨大机动过程中各关键组件的载荷,计算结果可对大柔性电池翼的强度校核和减重设计起到一定指导作用.

简介：在状态空间下,将线性陀螺系统微振动问题导向哈密顿体系,可以得到一组加权共轭辛正交关系和模态展开定理.利用这种特点构造了陀螺系统模态摄动计算式与灵敏度计算式,从而解决了拉格朗日体系下陀螺系统模态摄动分析与灵敏度计算的困难,算例显示了文中计算方法的有效性.

简介：提出了一种快速计算变截面铁木辛柯梁横向振动特性的方法．基于铁木辛柯梁理论建立的变截面梁的横向振动方程，其梁的截面参数如有效剪切面积、密度、弯曲刚度、转动惯量等沿梁轴线连续或非连续变化；首先将变截面梁等效为多段均匀阶梯梁；然后基于相邻两段连接处的位移（位移、转角）和力（弯矩、剪力）连续条件，建立相邻两段模态函数间相互关系，并递推出首段段与末段模态函数相互关系，利用边界条件得到相应特征方程，使用Newton—Raphson方法计算其固有频率；最后针对梁常见边界条件，得到计算变截面铁木辛柯梁横向振动固有频率特征方程的具体形式．用该方法计算-变截面梁在常见边界条件下前三阶固有频率．将计算结果同有限元计算结果进行比较，验证所提方法的有效性．然后与欧拉-伯努利梁计算结果比较，验证了本文方法求解短粗梁固有频率具有更好适用性．

简介：在外弹道数据处理中，奇异点处理、特征点求取与随机误差削弱都是精度估计的关键环节．本文首先利用小波变换在处理奇异点、特征点、噪声消除方面的优势，对观测数据进行基于小波变换的分解、融合、重构处理，剔除奇异点，查找特征点，削弱随机误差．其次利用节点自由分布B样条描述导弹运动轨迹，使该弹道确定方法转化为关于求解导弹轨道样条表示参数和测量系统误差的多模融合的非线性优化问题，采用非线性最优化方法，进而得到待估参数的最优估计，完成弹道的最佳逼近．仿真结果表明，该技术应用在奇异点处理、特征点提取与随机误差削弱方面效果较好，多模融合算法能减少计算量，且能切实提高参数估计精度．

简介：大型钢筋混凝土结构做长期的结构健康监测,需要综合考虑混凝土的收缩徐变影响.对应变传感器采集的原始应变数据做适当的修正,可提高监测的准确度.本文以广州新电视塔结构健康监测中应变修正的计算为例,阐述了美国混凝土学会（ACI）推荐的收缩徐变分析理论在结构健康监测中的应用,并结合实验室的试件测试数据,修正徐变和收缩的具体计算公式,最终得到修正后的实测应变与理论应变吻合得很好.

简介：通过设置涡核模型的角度条件,使涡核模型在极限状态下仍保持收敛,进而改进了利用Biot-Savart定律计算直线涡元对空间任意一点诱导速度的模型;桨叶气动模型采用Weissinger-L升力面理论模拟;自由尾迹的求解采用PIPC松弛迭代算法结合具有二阶精度的CB2D时间步进算法.利用上述模型和算法对某型号旋翼的尾迹进行数值计算,结果显示：利用改进涡核模型计算的桨尖涡径向位移收缩更加明显,这与实际情况更加接近;利用改进涡核模型得到的自由尾迹结果与实验数据吻合的更好.上述结论可以证明,新建立的旋翼自由尾迹模型提高了原有模型的准确性.

简介：基于经典的Magnus级数方法提出了一个简单有效的四阶近似积分格式，用于求解一般非线性动力学系统．它是一种几何积分方法，能保持精确解的许多定性性质，并且该方法只包含二个或三个指数矩阵的乘积，避免了通常的Magnus级数方法涉及的复杂的交换子运算．数值算例显示该方法是有效的。