基于物理模型的实时海面模拟_邱捷
基于时间卷积网络的长江下荆江航道水位预测

第 6 期水 利 水 运 工 程 学 报No. 6 2023 年 12 月HYDRO-SCIENCE AND ENGINEERING Dec. 2023 DOI:10.12170/20220809002李港,李有为,舒章康,等. 基于时间卷积网络的长江下荆江航道水位预测[J]. 水利水运工程学报,2023(6):84-92. (LI Gang, LI Youwei, SHU Zhangkang, et al. Water level prediction of lower Jingjiang Waterway in Yangtze River based on temporal convolution network[J]. Hydro-Science and Engineering, 2023(6): 84-92. (in Chinese))基于时间卷积网络的长江下荆江航道水位预测李港1, 2,李有为1, 2,舒章康3,张宇3,王江1, 2,查伟3(1. 长江航道勘察设计院(武汉)有限公司,湖北武汉 430040; 2. 国家内河航道整治工程技术研究中心,湖北武汉 430040; 3. 南京水利科学研究院水灾害防御全国重点实验室,江苏南京 210029)摘要: 航道水位的精准预测对保障船舶通航安全具有重要意义。
以长江下荆江航道为研究区域,采用2019—2020年水文时间序列数据为训练集,2021年数据为验证集,构建基于时间卷积网络(TCN)的长江下荆江水位变化预测模型,并与基于长短时记忆神经网络(LSTM)和支持向量机(SVM)的水位预测模型进行计算精度的对比。
结果表明:不同站点TCN对应的最优输入时间窗口存在一定差异,监利站、调弦口站及石首站的最优输入时间窗口分别为前2 d、前2 d和前3 d;TCN在2021年下荆江各站点水位预测结果的纳什系数(E NS)和决定系数(R2)均高于0.995,均方根误差(E RMS)小于0.21 m,整体预测效果优于LSTM,两者预测精度均较高,均显著优于SVM;但随着预测时间尺度的增加,水位预测精度整体呈下降趋势;TCN模型各站点枯水期大部分时段的水位预测绝对误差小于0.2 m,这表明TCN在航道水位预测领域具有较好的应用潜力。
舰船与海面红外热像仿真

舰船与海面红外热像仿真
陈彬;谭建宇
【期刊名称】《节能技术》
【年(卷),期】2015(033)006
【摘要】模拟舰船与海面红外热像可以直观定性地分析舰船目标与海面背景的对比特征,判别舰船目标识别的特征部位,具有强烈的实际军事需求和应用背景.本文首先根据坐标变换及投影法并结合大气传输模型编写程序模拟了8 ~12 μm波段舰船的红外热像.可以看出烟囱出口是舰船探测最主要的识别部位.然后使用Cox-Munk粗糙海面模型和遮蔽函数并结合蒙特卡洛方法编写程序模拟不同风速条件下海面的红外热像,得到海面平均辐射亮度随风速增大而减小的结论.最后将舰船目标与海面背景的红外图像进行融合,分析不同探测距离条件下海面背景下舰船目标的红外热像特征,计算结果表明,随着探测距离的增大,舰船整体辐射亮度减小.【总页数】5页(P535-538,545)
【作者】陈彬;谭建宇
【作者单位】哈尔滨工业大学(威海)汽车工程学院,山东威海264209;哈尔滨工业大学(威海)汽车工程学院,山东威海264209
【正文语种】中文
【中图分类】U675.79
【相关文献】
1.基于O GRE的海面舰船目标红外仿真方法 [J], 杨敏;李敏;易亚星;杨壹斌;王亚楠;卢羿
2.海面舰船目标SAR雷达成像仿真研究 [J], 纪纲;时立攀
3.海面舰船目标检测优化算法仿真 [J], 王培元;周建军;王日胜;陈杰
4.针对海面舰船目标的RBM和DBS雷达成像仿真研究 [J], 纪纲;陈明
5.机载红外热像仪探测海面温度场的仿真计算 [J], 金晟;张晓晖;夏珉
因版权原因,仅展示原文概要,查看原文内容请购买。
基于联邦学习的船舶AIS轨迹预测算法

基于联邦学习的船舶AIS轨迹预测算法郑晨俊;曾艳;袁俊峰;张纪林;王鑫;韩猛【期刊名称】《计算机工程》【年(卷),期】2024(50)2【摘要】联邦学习是一种可以在弱通信环境下有效解决数据孤岛问题的分布式机器学习方法。
针对海上船舶轨迹实时预测问题,提出基于Fedves联邦学习框架与卷积神经网络-门控循环单元(CNN-GRU)模型的船舶轨迹预测算法(E-FVTP)。
根据Fedves联邦学习框架,通过规范客户端数据集规模以及客户端正则项,在保持原有客户端数据特征的前提下,减小数据非独立同分布特征对全局模型的影响,加快收敛速度。
面向海洋通信资源短缺场景,建立基于船舶自动识别系统(AIS)数据的CNN-GRU船舶轨迹预测模型,解决了船舶终端设备计算能力不足的问题。
在MarineCadastre开源和舟山海洋船舶航行AIS数据集上的实验结果表明,E-FVTP 在预测误差比集中式训练降低40%的情况下,收敛速度提升67%、通信代价降低76.32%,可实现复杂海洋环境中船舶轨迹的精确预测,保障海上交通安全。
【总页数】10页(P298-307)【作者】郑晨俊;曾艳;袁俊峰;张纪林;王鑫;韩猛【作者单位】杭州电子科技大学计算机学院;杭州电子科技大学复杂系统建模与仿真教育部重点实验室;数据安全治理浙江省工程研究中心;杭州电子科技大学圣光机联合学院【正文语种】中文【中图分类】TP181【相关文献】1.基于卡尔曼滤波算法船舶AIS轨迹估计研究2.基于多项式卡尔曼滤波的船舶轨迹预测算法3.基于联邦学习的船舶AIS轨迹谱聚类算法研究4.基于联邦随机森林的船舶AIS轨迹分类算法5.基于联邦深度强化学习的多无人机轨迹规划算法因版权原因,仅展示原文概要,查看原文内容请购买。
基于CFD

第8卷㊀第6期2023年11月气体物理PHYSICSOFGASESVol.8㊀No.6Nov.2023㊀㊀DOI:10.19527/j.cnki.2096 ̄1642.1088基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析李泳德ꎬ㊀郭㊀力ꎬ㊀季㊀辰(中国航天空气动力技术研究院ꎬ北京100074)CharacterizationofTransonicAerodynamicDampingofRocketsBasedonCFD/CSDCouplingLIYong ̄deꎬ㊀GUOLiꎬ㊀JIChen(ChinaAcademyofAerospaceAerodynamicsꎬBeijing100074ꎬChina)摘㊀要:随着新型大推力火箭的发展ꎬ弯曲模态频率的不断降低ꎬ以及流动分离和跨声速飞行时产生的激波震荡等因素ꎬ其在跨声速飞行过程中更容易出现非定常振动发散ꎮ文章以某带助推的运载火箭模型为研究对象ꎬ通过数值计算获取火箭强迫振动时的气动阻尼ꎬ并对影响火箭气动阻尼的因素进行了分析ꎮ包括结构节点位置㊁振动振幅大小㊁脉动压力等ꎮ研究表明:助推主要起到增大气动阻尼的作用ꎻ前节点主要影响收缩段的气动阻尼ꎻ振动振幅大小和脉动压力对气动阻尼的影响可忽略不计ꎮ关键词:气动阻尼ꎻ数值计算ꎻ跨声速ꎻ气动弹性ꎻ运载火箭㊀㊀㊀收稿日期:2023 ̄09 ̄25ꎻ修回日期:2023 ̄10 ̄23第一作者简介:李泳德(1995 ̄)㊀男ꎬ工学硕士ꎬ助理工程师ꎬ主要研究方向为气动弹性分析ꎮE ̄mail:562064169@qq.com通信作者简介:季辰(1982 ̄)㊀男ꎬ工学博士ꎬ研究员ꎬ主要研究方向为气动弹性力学ꎮE ̄mail:jichen167@hotmail.com中图分类号:V475.1㊀㊀文献标志码:AAbstract:Withthedevelopmentofnewhigh ̄thrustrocketsꎬthedecreasingfrequencyofthebendingmodesoftherocketꎬaswellasthefactorssuchasflowseparationandshockoscillationsgeneratedduringtransonicflightmakeitmorepronetonon ̄constantvibration.Inthispaperꎬalaunchvehiclemodelwithboostwastakenastheresearchobjectꎬandtheaerody ̄namicdampingoftherocketduringforcedvibrationwasobtainedthroughnumericalcalculation.Thefactorsaffectingtheaerodynamicdampingoftherocketwereanalyzedꎬincludingthepositionofstructuralnodesꎬthemagnitudeofvibrationam ̄plitudeꎬpulsatingpressureandsoon.Thestudyshowsthattheboostmainlyplaystheroleofincreasingaerodynamicdamp ̄ingandthefrontnodemainlyaffectstheaerodynamicdampingofthecontractionsection.Thevibrationamplitudesizeandthepulsatingpressurehaveanegligibleeffectontheaerodynamicdamping.Keywords:aerodynamicdampingꎻnumericalcalculationꎻtransonicꎻaeroelasticityꎻlaunchvehicle引㊀言通常情况下人们认为气动力对火箭的振动起到阻尼作用ꎬ即气动阻尼为正值ꎮ然而随着大推力火箭发展ꎬ火箭的长细比逐渐加大ꎬ导致弯曲刚度越来越小ꎬ同时为了满足有效载荷的外形要求ꎬ火箭头部整流罩尺寸不断加大ꎬ后续箱体的直径却保持不变ꎬ形成了典型的锤头体外形ꎮ国内外大量的火箭研制经验表明[1 ̄9]ꎬ对于此类锤头体外形火箭的气动设计ꎬ必须要进行动态气动载荷与动态气弹稳定性分析ꎬ否则设计的疏忽可能会导致火箭结构出现毁灭性的破坏进而导致发射失败ꎮ目前常用的衡量气弹稳定性的方法是通过风洞试验来获取气动阻尼系数ꎮ早在1963年ꎬ美国国家航空航天局Ames研究中心(NASAAmesRe ̄searchCenter)采用半刚性模型开展试验研究[10]ꎬ获取火箭头部的气动阻尼来评估其稳定性ꎬ但这只能用来模拟火箭弯曲振型前节点之前部分的结构动力学特性ꎮ直到兰利研究中心(NASALangleyResearchCenter)开发了全弹性模型气动阻尼试验气体物理2023年㊀第8卷技术ꎬ其可以模拟整体的结构动力学特性以及气动外形ꎬ并应用于多款运载火箭研制[11 ̄15]ꎮ国内ꎬ中国航天空气动力技术研究院对气动阻尼问题开展过较多的研究[16 ̄20]ꎬ从模型设计方法㊁模型制作工艺㊁试验机构设计和数据处理等诸多方面ꎬ逐步改进实现了从半刚性模型到全弹性模型的过渡ꎬ并在多个型号上得到验证ꎮ然而通过风洞试验研究气动弹性问题ꎬ技术难度大ꎬ试验成本高ꎬ同时几乎不可能开展全尺寸试验ꎮ因此通过数值计算的方法开展相关研究是另一种重要的手段ꎮ刘子强等[21]实现了通过数值计算确定气动阻尼系数的技术和方法ꎬ并与试验结果进行对比ꎬ证实了该方法的可靠性ꎮ冉景洪等[22]通过模态数据结合准定常理论的方法分析了减阻杆加后体这一弹性结构的气动阻尼ꎬ结果表明减阻杆造成的分离流会对后体的气动阻尼系数产生影响ꎮ朱剑等[23]针对新一代捆绑式运载火箭发展了非结构网格下的气动阻尼计算方法ꎬ并分析了攻角㊁Mach数等参数对气动阻尼的影响ꎮ本文在之前的计算方法[23]的基础上采用IDDES模型ꎬ考虑脉动压力的影响ꎬ通过强迫振动的方式ꎬ针对捆绑式运载火箭的某一特定模态进行数值计算仿真ꎬ研究前节点位置ꎬ振动振幅ꎬ脉动压力等参数对气动阻尼的影响规律ꎮ1㊀计算方法图1为本文所用的捆绑式运载火箭的计算模型ꎬ是典型的锤头体结构ꎮ在跨声速阶段ꎬ其头部会产生激波造成激波边界层干扰ꎬ而在锤头体外形的过渡段会出现气流分离ꎮ为探究各部分气动阻尼的变化ꎬ将整个箭体分为头部㊁过渡段㊁弹身3个部分ꎮ图1㊀表面网格及区域划分Fig.1㊀Surfacegridandregiondivision1.1㊀流场仿真模型本文分别用Reynolds平均法(Reynolds ̄averagedNavier ̄StokesꎬRANS)和改进的延迟分离涡模拟(improveddelayeddetached ̄eddysimulationꎬID ̄DES)[24 ̄25]进行计算ꎬ在RANS方程中ꎬ将变量分为平均值和波动值两部分ꎬ对于速度分量有ui=ui+uᶄi其中ꎬi=1ꎬ2ꎬ3ꎬui和uᶄi分别代表平均量和波动量ꎬ对于压强和其他标量也采用类似的形式ꎬ将这种形式代入连续性方程和动量方程中ꎬ并写成张量形式∂ρ∂t+∂∂xi(ρui)=0(1)∂∂t(ρui)+∂∂xj(ρuiuj)=∂p∂xi+∂∂xjμ∂ui∂xj+∂uj∂xi-23δij∂uk∂xkæèçöø÷éëêêùûúú+∂∂xj(-ρuᶄiuᶄj)(2)其中ꎬiꎬjꎬk可分别取1ꎬ2ꎬ3ꎻρ是密度ꎻt是时间ꎻ当i=j时δij取0ꎬ否则取1ꎮ式(1)㊁(2)是RANS方程ꎬ由方程可知RANS方法将湍流脉动对平均流动的作用模化为Reynolds应力项即-ρuᶄiuᶄjꎬ之后采用湍流模型进行封闭ꎬ本文采用的湍流模型为SSTk ̄ω模型ꎬ其输运方程为∂∂t(ρk)+∂∂xi(ρkui)=∂∂xjΓk∂k∂xjæèçöø÷+Gk-Yk∂∂t(ρω)+∂∂xi(ρωui)=∂∂xjΓω∂ω∂xjæèçöø÷+Gω-Yω其中ꎬk和ω分别代表湍流动能和湍流耗散率ꎬΓk和Γω分别代表k和ω的有效扩散系数ꎬGk和Gω分别代表k和ω的生成率ꎬYk和Yω分别代表k和ω的耗散率ꎮ因此RANS方法只能计算大尺度的平均流动ꎬ本文采用IDDES方法计算脉动压力对气动阻尼的影响ꎮIDDES方法是由分离涡模拟(detached ̄eddysimulationꎬDES)方法改进而来ꎬ其本质思想与DES方法相同ꎬ是想以网格尺度和模型中的特征尺度隐式划分RANS和大涡模拟(large ̄eddysimulationꎬLES)区域ꎬ使其既能处理RANS方法无法得到的脉动场ꎬ也能降低LES方法在模拟高Reynolds数流动时所需的计算资源ꎮ区别在于当边界层较厚或者分离区域较窄时ꎬDES方法会出现如模型应力损耗(modeledstressdepletionꎬMSD)ꎬ网格诱导分离(grid ̄inducedseparationꎬGIS)以及对数层不匹配(logarithmic ̄layermismatchꎬLLM)问题[24]ꎬ而IDDES模型通过改良计算区域划分ꎬ结合延迟分离涡模拟(delayeddetached ̄eddysimulationꎬDDES)和03第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析壁面模型大涡模拟(wall ̄modeledlarge ̄eddysimula ̄tionꎬWMLES)ꎬ定义新的长度尺度解决了这些问题ꎬ具体公式详见文献[25]ꎮ流场网格如图2㊁图3所示ꎬ边界层采用棱柱层结构ꎬ并调整第1层网格高度使得y+小于1ꎬ远场部分采用六面体结构网格ꎬ与边界层的过渡层采用非结构网格ꎮ整体网格单元数量为4.2ˑ106ꎮ图2㊀y方向截面网格示意图Fig.2㊀Schematicdiagramofcross ̄sectionalgridinthey ̄direction图3㊀x方向截面网格示意图Fig.3㊀Schematicdiagramofcross ̄sectionalgridinthex ̄direction物面边界条件为无滑移壁面条件ꎬ远场采用压力远场边界条件ꎬ湍流模型采用SSTk ̄ω模型ꎬ采用密度基求解ꎬ气体黏性采用Sutherland定律ꎬ空间离散采用2阶迎风格式ꎬ对流通量采用Roe格式ꎮ1.2 结构分析模型结构与流场耦合分析过程中ꎬ结构部分可以采用模态方法描述ꎮ结构模态可以通过有限元方法与结构模态试验方法获得ꎮ本文采用有限元分析结果获得的模态ꎬ图4所示为结构的前3阶模态ꎬ本文只分析计算结果中气动阻尼最小的第2阶模态ꎮ(a)f=1.200Hz(b)f=2.460Hz(c)f=2.957Hz图4㊀结构的前3阶模态Fig.4㊀Firstthreemodesofthestructure由于火箭结构外形简单ꎬ一般不考虑其扭转影响ꎬ因此可以将其简化为简单的梁模型ꎬ这样就可以给出其模态振动方程q㊆i+2biωiq˙i+ω2iqi=fi(3)式中ꎬqi为第i阶模态的广义位移ꎬbi为第i阶模态的结构阻尼系数ꎬωi为第i阶模态的固有频率ꎬ13气体物理2023年㊀第8卷fi为第i阶模态下质量归一化的广义气动力ꎮ若将fi按照Taylor展开并略去高阶项ꎬ可以将其转化为气动阻尼项与气动刚度项的形式ꎬ则式(3)可写为q㊆i+2(bi+Bi)ωiq˙i+(Ki+1)ω2iqi=0(4)式中ꎬBi为气动阻尼系数ꎬKi为气动刚度系数ꎬ研究表明[26]ꎬ气动刚度相对于结构刚度为小量可以忽略不计ꎬ而在计算中结构阻尼往往设置为0ꎬ因此气动阻尼可以直接反映其气弹稳定性ꎮ1.3㊀气动阻尼分析原理气动阻尼的分析可以采用强迫振动或者自由振动的方式进行ꎬ这两种方法获得的时域数据不同ꎬ提取气动阻尼的方式也不同ꎮ强迫振动方法初始演化过程较短ꎬ因此计算量较小ꎬ同时能够分析某一种振动形式的气动阻尼ꎬ明确该振动形式是收敛还是发散ꎮ分析过程中能够获得不同部位与部件的气动阻尼ꎮ但是对于多模态相互作用引起的发散(例如颤振)较难预测ꎮ自由振动方法需要一定的自由演化时间才能够对时域数据进行分析ꎬ不过自由振动方法能够获得最能够吸收能量的模态及其振动频率ꎮ对于本研究所关注的问题ꎬ气动载荷对结构振动的过程中气动阻尼的影响较大ꎬ而对气动刚度与气动质量影响较小ꎬ即结构的固有振动频率受到来流的影响较小ꎬ其稳定性问题主要由气动阻尼的正㊁负引起ꎬ所以采用强迫振动方法分析ꎮ强迫振动下结构做简谐模态振动qi(t)=Asin(ωit)式中ꎬA表示振动的振幅ꎬ将其代入计算气动力的公式中[21]并做正交积分可得Bi=ʏl0Bx(x)dx=-1MiAω2iTʏl0ʏt0+Tt0G(xꎬt)cos(ωit)dtdx(5)式中ꎬMi为第i阶模态的模态质量ꎬT为整数倍周期ꎬG为广义气动力ꎮ根据式(5)便可以得到局部或分区域的气动阻尼ꎮ1.4㊀耦合计算流程首先进行模态分析ꎬ以确定结构的模态频率与振型ꎬ用以设计强迫振动的频率和振幅ꎮ非定常流场计算前先进行定常流场计算ꎬ来加快非定常计算的演化速度并增强收敛性ꎬ结构节点位移通过径向基函数(RBF)插值方法[27]映射到气动网格节点上ꎬ来进行网格的变形ꎬ这里径向基函数选用WendlandC2ꎬ如下所示φ(x)=(1-x)4(4x+1)最后将计算出来的广义力提取出来ꎬ截取演化完毕的整数倍周期ꎬ进行气动阻尼计算ꎮ耦合计算流程图如图5所示ꎮ图5㊀耦合计算流程图Fig.5㊀Flowchartofcoupledcalculation2㊀结果分析与讨论2.1㊀流场分析结果计算的来流Mach数范围为0.7~1.2ꎮ其中中截面的压力分布如图6所示ꎮ可以看出在头部出现了膨胀波以及跨声速激波ꎬ在过渡段存在流动分离ꎬ随着Mach数的增大ꎬ头部低压区域逐渐扩张ꎬ并且能明显看到ꎬ在流动再附的位置产生了再附激波ꎮ(a)Ma=0.7023第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析(b)Ma=0.75(c)Ma=0.80(d)Ma=0.85(e)Ma=0.88(f)Ma=0.90(g)Ma=0.92(h)Ma=0.96(i)Ma=0.9833气体物理2023年㊀第8卷(j)Ma=1.00(k)Ma=1.05(l)Ma=1.10图6㊀不同Mach数下的中截面压力分布Fig.6㊀PressuredistributioninthemiddlesectionatdifferentMachnumbers2.2 气动阻尼分布通过上述流场分析ꎬ可以看出火箭不同部位流动结构并不相同ꎬ在头部与箭身上ꎬ流动主要为附着流动ꎬ而在过渡段会出现较为复杂的波系结构以及流动分离ꎮ针对不同的流动结构随流向站位x的变化ꎬ设该位置上广义力与广义位移的相位差为φ(x)ꎬ并且简谐振动没有引入其他模态的广义力ꎬ则广义力的表达式为G(xꎬt)=Fgen sin[ωt+φ(x)]+F0(6)其中ꎬFgen为广义力的振动幅度ꎬF0为广义力的常数偏移量ꎮ将式(6)代入到式(5)中得到B(x)=-FgenMAω2Tʏt0+Tt0sin[ωt+φ(x)]cos(ωt)dt其中ꎬ广义力的常数偏移量F0的积分为0ꎬ因此省略ꎮ通过将等式中的正弦函数部分进行和差化积得到B(x)=-FgenMAω2Tʏt0+Tt0sin(ωt)cos[φ(x)]cos(ωt)dt+[ʏt0+Tt0sin[φ(x)]cos(ωt)cos(ωt)dt](7)式(7)中第1部分在整个周期中的积分为0ꎬ只有第2部分保留ꎬ因此得到B(x)=-Fgensin[φ(x)]MAω2Tʏt0+Tt0cos2(ωt)dt(8)式(8)中积分部分恒为正值ꎬ决定整个气动阻尼的部分只有相位角φ(x)的正弦值sin[φ(x)]ꎬ为了能够更加直观地获得相位角与气动阻尼B之间的关系ꎬ须将符号转化为对应的正弦函数转角ꎬ根据正弦关系ꎬ此转角为πꎬ因此得到B(x)=-Fgen(x)sin[φ(x)+π]MAω2Tʏt0+Tt0cos2(ωt)dt(9)图7为气动阻尼变化曲线ꎬ可以看出随着Mach数的增大ꎬ整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ图7㊀有助推时气动阻尼变化曲线Fig.7㊀Aerodynamicdampingchangecurvewithboost根据式(9)ꎬ得到相位角与气动阻尼B之间的关系为:当φ(x)ɪ(-πꎬ0)时ꎬ相位角滞后ꎬ气动阻尼B为负值ꎻ当φ(x)ɪ(0ꎬπ)ꎬ相位角提前ꎬ43第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析气动阻尼B为正值ꎻ为当φ(x)=0时ꎬ无相位角差别ꎬ气动阻尼B为0ꎮ在过渡段上ꎬ复杂的波系结构以及流动分离ꎬ使得气动力与结构位移之间会出现较为明显的迟滞现象ꎬ从而导致相位角φ(x)ɪ(-πꎬ0)ꎬ由此在过渡段上产生了负的气动阻尼ꎮ计算过程中的广义力与广义位移随时间变化曲线如图8所示ꎬ可以看出所有工况计算结果都表现良好ꎬ需要注意的是在非定常计算初期ꎬ演化的不完全导致广义力存在一些突变异常的结果ꎬ计算气动阻尼时须剔除ꎬ选择后面演化完全的周期ꎮ本文计算了9个周期ꎬ剔除了第1个周期出现的错误结果ꎬ采用后8个周期进行气动阻尼分析ꎮ强迫运动振幅为芯级直径的0.5%ꎮ(a)Ma=0.70㊀㊀㊀(b)Ma=0.75(c)Ma=0.80㊀㊀㊀(d)Ma=0.85(e)Ma=0.88㊀㊀㊀(f)Ma=0.9053气体物理2023年㊀第8卷(g)Ma=0.92㊀㊀㊀(h)Ma=0.96(i)Ma=0.98㊀㊀㊀(j)Ma=1.00(k)Ma=1.05㊀㊀㊀(l)Ma=1.10图8㊀不同工况下的广义力与广义位移随时间变化曲线Fig.8㊀Timedependentcurvesofgeneralizedforceandgeneralizeddisplacementunderdifferentoperatingconditions2.3㊀气动阻尼影响因素2.3.1㊀有无助推对气动阻尼的影响捆绑式运载火箭相比于传统的运载火箭ꎬ最大的区别就是在尾部四周捆绑了助推器ꎬ使得其流场特性变得复杂ꎬ因此须分析其对气动阻尼的影响ꎮ图7㊁图9分别为有无助推时气动阻尼变化曲线ꎬ可以看出随着Mach数的增大整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过63第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ对比两个图可知ꎬ助推主要起增大气动阻尼的作用ꎮ还可以看出有无助推情况下头部的气动阻尼变化很小ꎬ意味着在箭体尾部施加控制很难影响到头部的气动阻尼ꎬ特别是在超声速流场中ꎮ图9㊀无助推时气动阻尼变化曲线Fig.9㊀Aerodynamicdampingchangecurvewithoutboost2.3.2㊀前节点位置影响为了考察前节点位置变化对气动阻尼的影响ꎬ在保持振动频率不变㊁头部最大振型位置与振幅不变的条件下移动前节点ꎬ变化后的振型如图10所示ꎮ(a)Frontnodeafterthetransitionregion(b)Frontnodeinthetransitionregion(c)Frontnodebeforethetransitionregion图10㊀前节点变化后的振型Fig.10㊀Vibrationmodeafterthechangeofformernode根据对计算结果的分析分别获得了不同前节点位置的整体气动阻尼对比与过渡段气动阻尼对比ꎬ如图11㊁图12所示ꎬ可以看出前节点位置的改变并没有影响整体气动阻尼随Mach数增大而增大的趋势ꎬ且前节点在过渡段上与过渡段前的整体气动阻尼相差不大ꎬ而前节点在过渡段后的整体气动阻尼要高于另两种情况ꎬ因此过渡段与头部放在同一侧有助于提高气动阻尼ꎮ过渡段的气动阻尼会随着前节点的变化发生剧烈改变ꎬ前节点在过渡段前后随Mach数增大的变化规律相反ꎬ节点前后的振动相位变化导致不同节点位置过渡段的振动相位不同ꎬ进而导致气动阻尼发生变化ꎮ图11㊀不同节点位置的整体气动阻尼Fig.11㊀Overallaerodynamicdampingatdifferentnodepositions图12㊀不同节点位置的过渡段气动阻尼Fig.12㊀Aerodynamicdampingofthetransitionregionatdifferentnodepositions2.3.3㊀强迫振动振幅大小对气动阻尼的影响为了考察强迫振动振幅大小对气动阻尼的影响ꎬ在保证流场结构不发生改变的前提下ꎬ振动振幅分别为原来的一半和两倍ꎬ根据工程经验ꎬ如果振幅超过芯级直径的5%ꎬ则须考虑流场结构改变所造成的影响ꎮ图13㊁图14分别为不同振幅下的整体与头部气动阻尼ꎮ73气体物理2023年㊀第8卷图13㊀不同振幅下整体气动阻尼Fig.13㊀Overallaerodynamicdampingatdifferentamplitudes图14㊀不同振幅下头部气动阻尼Fig.14㊀Aerodynamicdampingoftheheadregionatdifferentamplitudes可以发现改变振幅无论是对整体气动阻尼还是头部气动阻尼来说变化都很小ꎬ这意味着气动阻尼的大小主要取决于气动力与结构振动的相位差ꎬ不依赖于振动幅度的大小ꎮ2.3.4㊀脉动压力对气动阻尼的影响为了模拟出脉动压力的影响ꎬ采用IDDES方法对火箭气动阻尼进行计算ꎬ计算来流Mach数为0.92ꎬ计算过程中的广义力与广义位移如图15所示ꎬ相较于图8可以看出广义力随时间变化曲线并不光滑ꎬ脉动压力的存在导致广义力由多个频率叠加而成ꎮ由于第2阶模态的频率为2.46Hzꎬ而由分离流㊁激波振荡等引起的脉动压力频率往往远大于此频率ꎬ因此这里选择3.5Hz为分界ꎬ将高于3.5Hz的部分视为由抖振脉动压力引起的广义力ꎬ低于3.5Hz的部分视为强迫振动引起的广义力ꎬ通过低通滤波把高于3.5Hz的广义力滤掉ꎬ可以获得由强迫振动引起的广义力与广义位移变化曲线ꎬ如图16所示ꎬ通过此广义力计算的气动阻尼为2.08ɢꎮ同样地ꎬ进行高通滤波将低于3.5Hz的广义力滤掉ꎬ可以获得由抖振脉动压力引起的气动阻尼为(2.94ˑ10-3)ɢꎬ由此得到脉动压力引起的气动阻尼变化为0.14%ꎬ可以忽略不计ꎮ同时使用RANS方法计算的气动阻尼为2.07ɢꎬ与IDDES的计算结果相比误差约为(2.94ˑ10-3+2.08-2.07)/2.07ʈ0.48%ꎬ这说明针对气动阻尼的模拟ꎬ抖振引起的脉动压力对气动阻尼的计算结果影响很小ꎬ起主要作用的还是广义力的变化ꎬ该变化由强迫振动引起的结构边界变化所导致ꎮ图15㊀基于IDDES的广义力与广义位移变化曲线Fig.15㊀VariationcuresofgeneralizedforceandgeneralizeddisplacementbasedonIDDES图16㊀滤波后的广义力与广义位移变化曲线Fig.16㊀Variationcuresofgeneralizedforceandgeneralizeddisplacementvariationcurveafterfiltering3㊀结论本文通过数值计算方法研究了火箭的气动阻尼特性ꎮ根据流动特征分析与理论推导ꎬ发现火箭过渡段几何外形的收缩导致该区域出现复杂的分离与激波结构ꎬ从而造成了气动力相对于结构振动83第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析相位的滞后ꎬ导致了该区域为气动负阻尼ꎬ即气动不稳定性的主要来源ꎮ在此机理的基础上ꎬ分析了前节点位置㊁振动振幅㊁脉动压力等因素对气动阻尼的影响规律ꎮ可以得出以下结论:1)助推增加了正阻尼区域的面积ꎬ从而相对于没有助推的构型起到了增加气动阻尼的作用ꎮ2)前节点位置的改变对过渡段气动阻尼影响很大ꎬ节点前后的振动方向相反ꎬ导致节点在过渡段前后的气动阻尼变化规律也截然相反ꎬ将过渡段与头部区域放在节点的同一侧有助于增加气动阻尼ꎮ3)在不改变流场结构的前提下ꎬ改变振动的振幅ꎬ气动力也会产生相应幅度的变化ꎬ因此结构振幅对气动阻尼的影响可忽略不计ꎮ4)高频部分的广义力对气动阻尼的贡献很小ꎬ即结构振动引起的广义力变化对气动阻尼起主要作用ꎬ而脉动压力对计算气动阻尼影响不大ꎬ可忽略不计ꎮ参考文献(References)[1]㊀CoeCF.Steadyandfluctuatingpressuresattransonicspeedsontwospace ̄vehiclepayloadshape[R].NASATMX ̄503ꎬ1961.[2]ColeSRJrꎬHenningTLꎬRaineyAG.NASAspaceve ̄hicledesigncriteria[R].NASASP ̄8001(REV)ꎬ1964. [3]EricssonLEꎬRedingJP.Analysisofflowseparationeffectsonthedynamicsofalargespacebooster[J].Jour ̄nalofSpacecraftandRocketsꎬ1965ꎬ2(4):481 ̄490. [4]RedingJPꎬEricssonLE.Effectofaeroelasticconsidera ̄tionsonseasat ̄Apayloadshrouddesign[J].JournalofSpacecraftandRocketsꎬ1981ꎬ18(3):241 ̄247. [5]程镇煌.宇航飞行器跨音速气动弹性问题探讨[J].上海航天ꎬ1997(6):16 ̄21.ChengZH.Discussionontransonicpneumaticelasticityofspacecraft[J].AerospaceShanghaiꎬ1997(6):16 ̄21(inChinese).[6]倪嘉敏.我国运载火箭气动设计回顾[C].近代空气动力学研讨会论文集ꎬ2005.NiJM.ReviewofaerodynamicdesignofChinaᶄslaunchvehicle[C].Proceedingsofthemodernaerodynamicssymposiumꎬ2005(inChinese).[7]臧涛成ꎬ胡焕性.大长细比弹箭弹性效应研究综述[J].弹道学报ꎬ1999ꎬ11(3):89 ̄93ꎬ96.ZangTCꎬHuHX.Areviewofgreatslendernessratioprojectileelasticeffectresearch[J].JournalofBallisticsꎬ1999ꎬ11(3):89 ̄93ꎬ96(inChinese). [8]吴志刚ꎬ杨超.细长体弹箭的气动弹性问题与研究方法[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.WuZGꎬYangC.Aeroelasticityproblemsandresearchmethodsofslenderbodymissiles[C].9thNationalSym ̄posiumonAeroelasticityꎬ2005(inChinese). [9]张贺ꎬ黄晓鹏.弹性细长旋转弹箭运动稳定性问题的研究进展[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.ZhangHꎬHuangXP.Researchprogressonthestabilityofelasticslenderrotatingprojectiles[C].9thNationalSymposiumonAeroelasticityꎬ2005(inChinese). [10]ColeHAJr.Dynamicresponseofhammerheadlaunchvehiclestotransonicbuffeting[R].NASATND ̄1982ꎬ1963.[11]BartelsREꎬWiesemanCDꎬMineckRE.ComputationalaeroelasticanalysisoftheAreslaunchvehicleduringas ̄cent[R].AIAA2010 ̄4374ꎬ2010.[12]AzevedoJLF.Aeroelasticanalysisoflaunchvehiclesintransonicflight[J].JournalofSpacecraftandRocketsꎬ1989ꎬ26(1):14 ̄23.[13]SinclairAꎬFlowersG.Low ̄orderaeroelasticmodeloflaunch ̄vehicledynamics[R].AIAA2010 ̄7725ꎬ2010. [14]DotsonKW.Transientcouplingoflaunchvehiclebendingresponseswithaerodynamicflowstatevariations[J].JournalofSpacecraftandRocketsꎬ2001ꎬ38(1):97 ̄104.[15]ColeSRꎬHenningTL.Buffetresponseofahammerheadlaunchvehiclewind ̄tunnelmodel[J].JournalofSpacecraftandRocketsꎬ1992ꎬ29(3):379 ̄385.[16]崔尔杰.流固耦合力学研究与应用进展[C].钱学森科学贡献暨学术思想研讨会论文集ꎬ2001.CuiEJ.Researchandapplicationprogressoffluid ̄struc ̄tureinteractionmechanics[C].ProceedingsofSeminarofQianXuesenScientificContributionsandAcademicThoughtsꎬ2001(inChinese).[17]冯明溪ꎬ王志安.火箭跨音速动导数和抖振实验[J].宇航学报ꎬ1987(1):55 ̄63.FengMXꎬWangZA.Experimentsoftransonicderiva ̄tivesandbuffetingofrocket[J].JournalofAstronauticsꎬ1987(1):55 ̄63(inChinese).[18]白葵ꎬ冯明溪.弹性模型实验技术[J].流体力学实验与测量ꎬ1999ꎬ13(1):38 ̄42.BaiKꎬFengMX.Aeroelasticmodelandthebuffetex ̄perimentaltechnique[J].ExperimentsandMeasurementsinFluidMechanicsꎬ1999ꎬ13(1):38 ̄42(inChinese). [19]JiCꎬRanJHꎬLiFꎬetal.Theaerodynamicdamping93气体物理2023年㊀第8卷testofelasticlaunchvehiclemodelintransonicflow[C].Proceedingsofthe64thInternationalAstronauticalCon ̄gressꎬ2013.[20]季辰ꎬ吴彦森ꎬ何岗ꎬ等.运载火箭气动阻尼风洞试验研究[C].第十二届全国空气弹性学术交流会论文集ꎬ2011.JiCꎬWuYSꎬHeGꎬetal.Experimentalstudyonaero ̄dynamicdampingwindtunneloflaunchvehicle[C].Pro ̄ceedingsofthe12thNationalAeroelasticityConferenceꎬ2011(inChinese).[21]刘子强ꎬ白葵ꎬ毛国良ꎬ等.锤头体弹性振动跨音速气动阻尼系数的确定[J].宇航学报ꎬ2002ꎬ23(6):1 ̄7.LiuZQꎬBaiKꎬMaoGLꎬetal.Thedeterminationofaerodynamicdampingonhammerheadlaunchvehiclesattransonicspeeds[J].JournalofAstronauticsꎬ2002ꎬ23(6):1 ̄7(inChinese).[22]冉景洪ꎬ刘子强ꎬ胡静ꎬ等.减阻杆气动阻尼研究[J].力学学报ꎬ2014ꎬ46(4):636 ̄641.RanJHꎬLiuZQꎬHuJꎬetal.Researchofaero ̄dampingforbluntwithspike[J].ChineseJournalofThe ̄oreticalandAppliedMechanicsꎬ2014ꎬ46(4):636 ̄641(inChinese).[23]朱剑ꎬ冉景洪ꎬ吴彦森ꎬ等.捆绑式运载火箭的气动阻尼数值计算方法[C].第十三届全国空气弹性学术交流会论文集.哈尔滨:中国力学学会ꎬ中国空气动力学会ꎬ2013.ZhuJꎬRanJHꎬWuYSꎬetal.Numericalcalculationmethodforaerodynamicdampingofbundlelaunchvehicles[C].Proceedingsofthe13thNationalAeroelasticityCon ̄ference.Harbin:ChineseSocietyofTheoreticalandAp ̄pliedMechanicsꎬChinaAeromechanicsSocietyꎬ2013(inChinese).[24]SpalartPR.Detached ̄eddysimulation[J].AnnualReviewofFluidMechanicsꎬ2009ꎬ41:181 ̄202. [25]GritskevichMSꎬGarbarukAVꎬSchützeJꎬetal.Devel ̄opmentofDDESandIDDESformulationsforthek ̄ωshearstresstransportmodel[J].FlowꎬTurbulenceandCombustionꎬ2012ꎬ88(3):431 ̄449.[26]季辰ꎬ吴彦森ꎬ侯英昱ꎬ等.捆绑式运载火箭跨声速气动阻尼特性试验研究[J].实验流体力学ꎬ2020ꎬ34(6):24 ̄31.JiCꎬWuYSꎬHouYYꎬetal.Experimentalstudyofaerodynamicdampingcharacteristicsofalaunchvehiclewithboostersintransonicflow[J].JournalofExperimentsinFluidMechanicsꎬ2020ꎬ34(6):24 ̄31(inChinese). [27]AllenCꎬRendallTCS.UnifiedapproachtoCFD ̄CSDinterpolationandmeshmotionusingradialbasisfunctions[R].AIAA2007 ̄3804ꎬ2007.04。
虚拟海洋环境中的舰船运动的仿真

虚拟海洋环境中的舰船运动的仿真
孙建华;邱志强
【期刊名称】《江苏科技大学学报(自然科学版)》
【年(卷),期】2001(015)004
【摘要】采用源汇法和STF切片法,计算了在不同浪向下舰船横向运动的水动力模型系数和波浪扰动力(矩),并存入数据库.提出采用四阶龙格-库塔法求解舰船运动水动力模型的线性微分方程组,获得了舰船随时间变化的各自由度运动幅值.采用虚拟现实技术,创建了天空和海面波浪等三维海洋虚拟场景,实现了具有真实感的、随时间变化的三维舰船运动(横摇、艏摇、横荡)在该场景中的动态仿真.
【总页数】5页(P15-19)
【作者】孙建华;邱志强
【作者单位】中船重工集团第719研究所;中船重工集团第719研究所
【正文语种】中文
【中图分类】TP317.4;U66.32
【相关文献】
1.虚拟仿真技术在舰船研制评估中的应用研究 [J], 冯佰威;李俊华;汪玉;高兴华;蔡薇
2.海洋环境中船舶运动规律的可视化仿真方法研究 [J], 肖凡;张涛;万峻
3.虚拟仿真技术在水声工程国防特色专业中的应用—以复杂海洋环境中的声呐被动探测性能分析实验为例 [J], 殷敬伟; 生雪莉; 高明生; 周天; 郭龙祥; 韩笑; 朱广平; 陈文剑
4.虚拟仿真实验在教学中的尝试与应用——以脊柱运动康复生物力学虚拟仿真实验为例 [J], 王钟梅;宋雅伟
5.虚拟海洋环境中潜艇六自由度运动的视景仿真 [J], 赵新华;孙尧;卢志忠
因版权原因,仅展示原文概要,查看原文内容请购买。
211133086_相似预报原理之再认识:动力统计相似集合预报理论及其对登陆台风降水预报的应用研究进

doi:10.11676/qxxb2023.20220064气象学报相似预报原理之再认识:动力统计相似集合预报理论及其对登陆台风降水预报的应用研究进展*任福民1贾莉1吴彩铭1丁晨晨1张大林1,2贾作3马蕴琦1邱文玉4 REN Fumin1 JIA Li1 WU Caiming1 DING Chenchen1 ZHANG Dalin1,2 JIA Zuo3 MA Yunqi1 QIU Wenyu41.中国气象科学研究院灾害天气国家重点实验室,北京,1000812.美国马里兰大学大气和海洋科学系,马里兰州,207423.中船航海科技有限责任公司,北京,1000704.南京信息工程大学,南京,2100441.State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences,Beijing 100081,China2.Department of Atmospheric and Oceanic Science,University of Maryland,College Park,Maryland 20742,USA3.CSSC Marine Technology Co.,Ltd.,Beijing 100070,China4.Nanjing University of Information Science and Technology,Nanjing 210044,China2022-04-24收稿,2022-10-30改回.任福民,贾莉,吴彩铭,丁晨晨,张大林,贾作,马蕴琦,邱文玉. 2023.相似预报原理之再认识:动力统计相似集合预报理论及其对登陆台风降水预报的应用研究进展.气象学报,81(2):193-204Ren Fumin, Jia Li, Wu Caiming, Ding Chenchen, Zhang Dalin, Jia Zuo, Ma Yunqi, Qiu Wenyu. 2023. Advances in dynamic-statistical analog ensemble forecasting and its application to precipitation prediction of landfalling typhoons: A renewed understanding. Acta Meteorologica Sinica, 81(2):193-204Abstract Combining dynamical and statistical methods is one of the important ways to improve weather and climate prediction. A key issue is how to effectively combine numerical model results with historical bining the above two methods with an analog method is an important direction for future improvement of weather and climate prediction, although the analog method is still limited to similarity assumptions and lacks solid physical basis. Based on the initial condition problem of a perfect model, this work proposes the concept of initial condition perturbation of the perfect model and develops a Dynamical-Statistical-Analog Ensemble Forecast (DSAEF) theory. The DSAEF theory shows not only why the analogue-based forecast can be conducted, but also how it can be conducted.That is, the perfect model is used to produce forecasts and an ensemble prediction scheme is used to achieve the forecast accuracy.Based on the DSAEF theory, the DSAEF_LTP (Landfalling Typhoon Precipitation, LTP) model has been developed. This model includes the following four steps:(Ⅰ) forecast typhoon track, (Ⅱ) Construct generalized Initial Value (GIV), (Ⅲ) identify analogs from historical observations, and (Ⅳ) produce an ensemble forecast of typhoon precipitation.The GIV is constructed by physical variables that affect LTP. The DSAEF_ LTP model has the characteristic of sustainable development, which can be improved by introducing new variables or refining the existing parameters of the model. At present, the model versions 1.0 and 1.1 have been released.In version 1.0, GIVs include three physical variables, i.e., typhoon track, landfall season and typhoon intensity. In the version 1.1, two extra improved parameters of "similarity region" and "ensemble scheme" are added. The most recent version 1.1 has shown significant improvements in the performance for forecasting LTP. Three large-sample forecast experiments using version 1.1 show that, compared ECMWF, CMA-GFS, NCEP-GFS and SMS-WARMS (Shanghai regional model), version 1.1*资助课题:国家重点研究发展计划项目(2019YFC1510205)、海南省南海气象防灾减灾重点实验室开放基金项目(SCSF202202)、国家自然科学基金面上项目(42275037)。
长江口外海域波浪场数值模拟
长江口外海域波浪场数值模拟邱桔斐;马越;徐新华【摘要】采用MIKE21的SW模块,选用麦莎台风和2006年1月的一次冷空气,对不同极端天气系统下长江口外海域的波浪场进行模拟,计算结果与实测值吻合较好,分析说明不同天气系统下波浪场的分布具有显著不同的特征,同时验证了模型在该海域有较好的适用性.【期刊名称】《水运工程》【年(卷),期】2011(000)010【总页数】4页(P11-14)【关键词】长江口附近;波浪场;台风;冷空气【作者】邱桔斐;马越;徐新华【作者单位】上海东海海洋工程勘察设计研究院,上海200137;上海东海海洋工程勘察设计研究院,上海200137;上海东海海洋工程勘察设计研究院,上海200137【正文语种】中文【中图分类】TV139.2+2长江口与杭州湾位于副热带季风气候区,冬季盛行偏北风、夏季盛行偏南风,季节性变化十分明显。
风向的季节变化必然引起浪向的变化,附近海域盛行浪向与盛行风向颇为一致,冬季以偏北向浪为主,夏季以偏南向浪为主,春秋两季为浪向交替过渡季节。
长江口和杭州湾海域最大波高往往产生在台风期间,其次在寒潮大风天气时。
在浙江海门以北、上海附近沿海登陆和在长江口外海面过境的台风,都会在本海域引起较大的台风浪。
台风浪的大小主要与台风路径、台风强弱及其移动速度等条件有关,同时还受到地形和水深的影响。
在该海域,前人也采用过不同波浪模式对该海域的波浪场进行数值模拟,其中胡克林等[1]以0215号鹿沙台风和0216号森拉克台风为例,将SWAN模型应用到长江口附近海域,进行台风浪的数值模拟;隋洪波[2]应用中国海洋大学开发的“二维浅水海浪谱数值模式”,选择东向入射波进行了口门及航道内较具代表性的5个点的模拟,初步分析了该数值模式在长江口海域的应用;冯卫兵[3]利用定常缓坡方程计算了长江口海域波浪的折射绕射,并分析了涨急落急两种流态对波浪场的影响;孔令双等[4]采用MIKE21 SW 波浪谱模型对长江口海域在凤凰台风期间的波浪场进行了模拟计算。
基于ARM9的海底观测网节点电路硬件设计
基于ARM9的海底观测网节点电路硬件设计吕斌;杜立彬;贺海靖;曲君乐;王秀芬【摘要】The paper devises a control system for ARM9 based sea floor observatory network node circuit to satisfy the real-time acquisition requirement of submarine observation data. This system employs the modular design method of an embedded system to construct a circuit system including extended UART communications circuit, analog-to-digital conversion circuit, ethernet communications circuit, voltage level conversion circuit and USB interface circuit. Simulation and application results prove that this node circuit system has better stability and reliability, and easy to be monitored and maintained by a shore station.%为满足海底观测实时获取数据的功能需求,设计了一种基于ARM9的海底观测网节点电路控制系统。
采用嵌入式系统硬件模块化设计,构建了包括节点串口扩展通讯、ADC 数据采集、以太网通讯、电平转换及USB接口的电路系统。
现场试验证明,该节点电路运行稳定可靠,且易于岸站监控与维护。
【期刊名称】《山东科学》【年(卷),期】2012(025)001【总页数】6页(P59-64)【关键词】海底观测网;节点电路系统;ARM9【作者】吕斌;杜立彬;贺海靖;曲君乐;王秀芬【作者单位】山东省海洋环境监测技术重点实验室,山东省科学院海洋仪器仪表研究所,山东青岛266001;山东省海洋环境监测技术重点实验室,山东省科学院海洋仪器仪表研究所,山东青岛266001;山东省海洋环境监测技术重点实验室,山东省科学院海洋仪器仪表研究所,山东青岛266001;山东省海洋环境监测技术重点实验室,山东省科学院海洋仪器仪表研究所,山东青岛266001;山东省海洋环境监测技术重点实验室,山东省科学院海洋仪器仪表研究所,山东青岛266001【正文语种】中文【中图分类】TP274.2当前的海洋观测,尤其是深海环境的海洋观测通常都是基于船载设备或定点投放观测设备来实施的,这种观测是短期的、离散的、断续的观测[1]。
海上无人系统跨域协同运用与技术发展
海上无人系统跨域协同运用与技术发展邱志明, 孟祥尧, 马 焱, 王 亮, 肖玉杰(海军研究院, 北京, 100442)摘 要: 海上无人系统跨域协同是未来海上无人系统的发展趋势和重要的应用方式。
随着各种海上无人系统的快速发展和在世界局部战争冲突中的应用, 如何更好地使用海上无人系统跨域完成作战任务成为研究的重点。
文中以不同空间域的海上无人系统为研究对象, 梳理总结了当前海上无人系统以及国外海上无人系统跨域运用的发展现状。
重点针对海上无人系统的跨域协同运用基本原理和方法进行了分析, 提出了海上无人系统跨域协同运用的关键问题, 并在此基础上梳理了技术发展中需要重点关注的关键技术。
最后提出了未来发展的几点启示, 以期为海上无人系统的跨域运用和技术发展提供参考和借鉴。
关键词: 海上无人系统; 跨域协同; 技术发展中图分类号: TJ6; U674 文献标识码: R 文章编号: 2096-3920(2024)02-0184-10DOI: 10.11993/j.issn.2096-3920.2024-0053Cross-domain Collaborative Application and Technology Development ofMaritime Unmanned SystemsQIU Zhiming, MENG Xiangyao, MA Yan, WANG Liang, XIAO Yujie(Naval Research Institute, Beijing 100442, China)Abstract: Cross-domain collaboration of maritime unmanned systems is the future development trend and important way of application of maritime unmanned systems. With the rapid development of various maritime unmanned systems and their application in local wars and conflicts in the world, how to better use maritime unmanned systems to complete cross-domain combat missions has become the focus of research. With maritime unmanned systems in different spatial domains as the research objects, the current maritime unmanned systems and the development status of cross-domain application of maritime unmanned systems in other countries were summarized. The basic principles and methods of cross-domain collaborative application of maritime unmanned systems were analyzed, and the key issues of cross-domain collaborative application of maritime unmanned systems were put forward. On this basis, the key technologies requiring attention in the development of technologies were sorted out. Finally, several enlightenments for future development were put forward, so as to provide a reference for the cross-domain application and technology development of maritime unmanned systems.Keywords: maritime unmanned system; cross-domain collaboration; technology development0 引言在新一轮科技革命和世界海军强国竞争的双重驱动下, 无人系统在海战中的地位日益突出。
海浪可视化仿真及军事应用
海浪可视化仿真及军事应用
龚琳;顾大权
【期刊名称】《气象水文海洋仪器》
【年(卷),期】2008(000)002
【摘要】阐述了计算机图形学领域中海浪的可视化仿真方法.详细介绍了基于几何多边形,基于波分析,基于谱分析,基于物理模型,基于粒子系统及基于光影处理的海浪建模方法.归纳了提高海浪仿真实时性的主要措施及研究动向.分析了现有海浪仿真模型的优缺点,预测了其发展趋势.阐述了海浪的可视化仿真方法在军事领域的应用.【总页数】5页(P48-52)
【作者】龚琳;顾大权
【作者单位】解放军镇江船艇学院,镇江,212003;解放军理工大学,气象学院,南京,211101
【正文语种】中文
【中图分类】TP391.9
【相关文献】
1.基于浅水海浪数值模型的近岸海浪仿真 [J], 肖剑波;胡大斌;胡锦晖
2.基于 OpenGL 的实时三维海浪可视化仿真 [J], 沈祉怡;赵冰化;张晓蕾
3.基于海浪谱的风暴潮海浪数值仿真与功率谱估计 [J], 陈小健;吴庆;张道明
4.基于Gerstner函数和JONSWAP海浪谱的海浪仿真 [J], 杨志钢;孙鹏
5.基于OpenGL的三维浅水海浪可视化仿真 [J], 石贱弟;姜昱明
因版权原因,仅展示原文概要,查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2
海洋模型原理
2. 1 海洋运动特征 海洋的表面运动其实是由一系列的波组成的 , 当海风作用于海面时 , 海面会产生主方向与海风同 向的波族。这些波的主要运动方向和海风同向, 但 是由于海面本身的流体性质, 这些波之间会通过黏 性阻力等内在作用力相互作用 , 产生如漩涡等自然 现象。并且内在的作用力会使得每个波都会沿与 海风不完全平行的方向进行扩散。在过去, 为了渲 染海面 , 通常使用一系列波长 , 振幅和频率都不相 同的正弦波作用于一个 M @ N 的顶点序列 , 每个顶 点的高度根据通过该点的所有正弦波在时刻 t 的 偏移进行叠加。这种技术能够生成外观上比较真 实的海面 , 而且性能的损耗可以控制在合理的范围 内。但是这种技术并不是在真正模拟海面, 它无法
根据海面所受到的力的变换来反映出真实的交互。 与本文所介绍的技术相比, 后者能够实时的根据海 面所受力的变化改变海面的状态 , 比如当风向改变 时, 海面所 产生的波浪也会相 应的改变方向 和大 小。 为了模拟海面的波浪变化, 需要计算出波浪的 波长和角度分布。这里将波长作为向量, 沿这不同 的角度分布, 被分解成数个分量波长的向量和。直 接计算波浪的波长非常困难 , 但是可以通过计算波 浪的频率来间接计算波长。下面将介绍 4 种常用 的物理模型 , 并分析其中一种模型 , 推导出详细的 计算步骤。 2. 2 海洋运动物理模型 海洋物理模型是海洋学家为了研究海洋而利 用数学方法抽象出来的一系列数学模型。这些模 型都是根据观察经验来建立的, 因此每个模型都有 一定的适用范围。没有模型能在所有的情况下正 确描述海面的运动特征。由于本文只着重于介绍 如何渲染出视觉真实的海面效果 , 所以将介绍 4 种 用于模拟海洋表面波形运动的物理模型。这四种 模型对于海面所受到的力不是特别大时, 能正确模 拟海面与外力之间的交互行为, 但是当外力变的很 大时 ( 如飓风 ) , 这种模型将可能产生错误的结果。 2. 2. 1 Pierson- Moskowitz model 模型 这个是目前最为著名的海洋模型 , 它将海面波 浪频率的分布表示为海面能量对频率的函数。该 模型可以真实的模拟出海面通过受到各种力( 如海 风) 而产生的外形改变。 E PM ( f ) = fm Aù g 2 5 ù ( ) 4 ) ( 1) 4 5 ù exp (4 f ( 2 ù P) f
1
当前主流海洋模型概述
目前, 国内外在研究海浪模拟时选择的海浪模
¹
收到本文时间 : 2005 年 9 月 15 日 基金项目 : 信息产业部电子信息 产业发展基金项目 ; 四川省科技厅研究项目
第 34 卷 ( 2006) 第 1 期
计算机与数字工程
5
率分布数组, 有 JONSW AP 谱模型所求得的海面斜率 比有其他模型所求得范围更广。 JONSW AP 谱模型 在反映有一定风浪的海面的实际情况时具有较大的 优势。JONSW AP 谱模型的缺点: JONSW AP 谱模型所 需要输入的可变参数较多 , 计算量相对较大, 需要时 间较长。另外, 它在模拟海面起伏较小的波浪时效 果还不太好。 近来 有 不 少 人 研 究 海 浪 的 模 拟 方 法。如: Peachey 采用正弦函数和二次函数的线性组合来模 拟波浪的外形, Foumier 采用参数曲面来模拟波峰卷 曲的波浪, 童若锋等则对雨点波、 涟漪和紊乱的短峰 波进行分类造型 , 模拟波浪的大致形状。这类方法 总的来说比较简单直观, 计算量小, 但是产生的效果 不好。另一类方法一般是从著名的流体动力学方程 组- - Navier Stockes 方程组出发, 通过一系列近似 和简化, 用数值方法求解流场来实现。如 : Kass 通过 对浅水方程加以简化来模拟水波纹动画, Chen 采用 数值迭代方法求解二维的 Navier Stockes 方程 , 徐迎 庆等从描述明渠不稳定流的 Saint Venant 方程组出 发来模拟波浪 , Foster 等则直接用数值方法求解三维 流场从而更加真实全面地模拟了流体运动。这类方 法相对而言, 效果比较真实, 使用范围较广泛 , 但计 算量往往很大, 效率较低。 由于算法复杂性和数据庞大性 , 海洋模拟在实 时系统中最难以解决的问题就是效率问题。降低 模型复杂度可以明显的提高系统效率, 但是同时也 降低了视觉效果。如何在保证模型真实度的前提 下降低实时系统的硬件要求成为当前图形研究工 作者能量关于频率 和角度的分布关系。 结合 Pierson - Moskow it 方程 ( 1) , 得到以下方程: E (f , H ) = EPM (f ) ù D( f , H) = ù exp(2ù 2. 3. 2 振幅计算方程 通过以上方程, 可以计算出波浪的频率, 为了计 算出波浪的振幅, 需要通过傅立叶变换进行转换: F( k ) = E( f , H ) dk kù df ( 10) Aù g 2 ( 2 ù P) 4f 5
6 - (f - f m ) )] } ( 3) 2 ù R2 ù f 2 m 其中 C= 3. 3 及 R = 0. 07( f < f m ) 、 R = 0. 09( f < f m) 。 从上式可以看出 , Jonswap model 其实是 Pierson [ exp ( - Moskowitz 的改进版, 但是同时也需要更多的计 算时间。 2. 2. 3 Kruseman model 模型 类似 Jonswap 模型, 该模型也计算出比 Pierson - Moskowitz 略高一点的顶点频率。 但是它用另外 一种方式来表示能量与波频率之间的关系。 Aù f - 5 Ek (f ) = a ù f - 6 (f - b ù f m ) 1- b m 0 fm Ff b ùf m F f F f m 0 F f F b ùf m
4
计算机与数字工程
第 34 卷
基于物理模型的实时海面模拟
邱 捷 陈雷霆
成都 610054 ) ( 电子科技大学 计算机科学与工程学院 摘
¹
要 : 海洋模拟是计算机图形学中的一 个重要分支。在过去的数年中 , 由于海洋模型 本身的复杂 和计算机 硬件性能
的限制 , 无法实时的模拟出逼真的海面效果。但随着计算机硬 件性能 的飞速发 展和物 理工程 等领域 对海面 模型研 究的深 入 , 在最近几年出现了多种实时渲染海面的技术。该文将介 绍一种 利用海 洋物理 模型模 拟海面运 动现象 的技术。 这种技 术的实时不是特别复杂 , 而且能够拥有良好的 性能 , 非常适合在实时三维引擎中生成大面积外观真 实的海面。 关键词 : 海洋模型 实时 计算机图形 三维引擎 流体模拟 中图分类号 :TP391. 9
型有: Pierson- Moscowitz 谱模型、 JONSW AP 谱模型。 通过我们前期对这些模型在计算机三维图形系统中 模拟的研究有如下研究结论。Pierson- Moscowitz 谱 模型的优点 : 随着风速的增大, 谱曲线下面积亦增 大, 且谱宽度亦增大。当风速增大时, 海浪谱的高频 部分增多 , 海浪变得更混沌、 更无规则性。该波谱经 改进后被第十一届国际船模水池会议定为标准海浪 谱。Pierson- Moscowitz 谱模型的缺点 : 只包括海浪 有效波高参数, 不足以表征复杂的海浪情况; 只适用 于充分发展的海浪, 不包括海底对海浪传播方向和 振幅的影响。JONSW AP 谱模型的优点: JONSW AP 谱 模型是目前国际上通用的标准。在诸多反映海浪的 模型中, JONSW AP 谱模型所反映的海面起伏相对较 大, 在现实应用中亦更能够接近实际的情况。此模 型的可变参数有: 风速、 风向、 风程 , 比其他的谱模型 更为精细 , 可以较为真实地反映海面多变的情况, 并 且能够比较合理地模拟海面的粗糙程度。JONSW AP 谱模型可以结合快速傅立叶变换技术求得海面的斜
0 . 13 ù g 。 g 是重力加速度 9. v 80665/ m 2 , v 是海风的速度。 f m 是海面波浪的顶点 其中 a = 0. 008、 fm = 频率。 通过 f m = 0. 13 ù g , 可以用散布方程求出一 v 定风速下波浪的波长。 散布 方程: X( k ) = = 2 ù Pù f 2. 2. 2 Jonswap model 模型 Jonswap 模型和 Pierson - Moskowit z 模型在形 式上很接近, 同样需要先计算出波浪的顶点频率。 但是其计算结果 比 Pierson - Moskowit z 要略 高一 些, 并且转换成能量函数后结构更加复杂。 E J ( f ) = EPM (f ) ù exp{ ln( C )ù g ùk = gù( 2ù P K ) ( 2)
Real Time Ocean Simulation Based on Physical Ocean Model
Qiu Jie Chen Leiting
( School of Computer Science and Engineering , UEST, Chengdu 610054) Abstract: Ocean Simulation is an important part of computer graphic. In the past year, since the hardware limition, there is not a way to render ocean with realtime system for goodlooking. But with the hardware running more faster than before, rendering ocean with reality appearance and acceptable performance becomes possible. T his article will discuss a technology based on some physical ocean model. It is not very difficult and has a fast speed. Key words: Ocean Model , Realtime Render, Computer Graphics, 3D Engine, Ocean Simulation Class number: TP391. 9