高速列车系统动力学耦合仿真方法技术方案

本发明专利技术涉及一种高速列车系统动力学耦合仿真方法，在基于接口的多设备参数、多领域间的相互作用关系及控制参数的环境中实现协同仿真基于仿真步长级的多向耦合控制，采用如下的手段，协调仿真步长不一致的各子系统实现耦合控制仿真。本发明专利技术全面考虑列车动力学、弓网动力学、空气动力学、车-线耦合动力学，车辆-空气耦合动力学等相关领域的影响，实现高速列车-线路-弓网-气动-牵引供电等领域间基于积分步长级的耦合计算，充分考虑高速列车与轨道线路、接触网、气流等周边环境的相互影响。

【技术实现步骤摘要】

本专利技术涉及仿真控制，尤其是多变量子系统实现耦合控制仿真

技术介绍
在高速运行环境下，高速列车的运行安全性是需要关注的主要目标之一，而这又取决于高速列车动力学决定的高速列车性能本身。高速列车运行环境复杂，动力学性能受到轨道线路、接触网、周边气流等的复杂耦合作用影响。然而在现有的高速铁路研究体系中，研究者们主要开展单子系统或局部耦合的仿真研究，单子系统的研究如高速列车动力学、高速弓网动力学、高速列车空气动力学等，对于高速列车耦合系统的研究多是基于两个系统的，比较典型的有高速列车和轨道线路(含桥梁)之间的车-线(桥)耦合动力学，高 速列车与周边流场间的流固耦合动力学，这些研究都只部分考虑了高速列车运行性能相关的影响因素，不能全面反映高速列车的运行品质。因此，为全面反映高速列车的运行品质，必须把高速列车与高速铁路其他相关子系统耦合起来进行作为整体进行高速列车耦合系统动力学仿真研究。而目前没有成熟的商业软件可以进行高速列车系统动力学的仿真计算研究，因此必须构建一种考虑多子系统相互作用的、且适合于当前高速列车性能仿真研究的耦合系统动力学仿真计算方法。
技术实现思路
鉴于现有技术的不足，本专利技术的目的是研究提供一种，使之克服现有技术的以上缺点。本专利技术的目的是通过如下的手段实现的。，在基于接口的多设备参数、多领域间的相互作用关系及控制参数的环境中实现协同仿真基于仿真步长级的多向耦合控制，采用如下的手段，协调仿真步长不一致的各子系统实现耦合控制仿真(I)根据子系统的仿真步长，将子系统划分为三层，第一层包括列车子系统、线路子系统和弓网子系统，该层的仿真步长记为S=5X10_5s，耦合仿真控制步长记为SI1=5X 10_5s，第二层为空气动力学子系统，该层的耦合仿真控制步长记为SI2=1 X 10_3s，第三层为牵引供电子系统，该层的耦合仿真控制步长记为SI3=1 X KT1S,并设置耦合仿真结束条件；(2)为每个层次设置一个累积仿真步长，第一层累积仿真步长记为SA1，第二层累积仿真步长记为SA2，第三层累积仿真步长记为SA3，并初始化SA1=0，SA2=0, SA3=0 ；(3)初始化各子系统间的耦合关系接口数据；(4)加载各子系统仿真模型，初始化各子系统；(5)各层累积仿真步长值增加第一层仿真步长值，并作如下判别I)如果第一层累积仿真步长值大于等于第一层耦合仿真控制步长值，即SAD=SIl,则更新第一层子系统间的耦合关系接口数据，激活第一层子系统进行单步长仿真计算，否则仿真退出并告警(告警是因为第一层是最小计算单元，如果该层无法计算，说明程序存在异常)；2)如果第二层累积仿真步长值大于等于第二层耦合仿真控制步长值，即SA2>=SI2，则更新第一和第二层子系统间的耦合关系接口数据，激活第二层子系统进行单步长仿真计算，否则跳到(7)，判别程序是否结束；3)如果第三层累积仿真步长值大于等于第三层耦合仿真控制步长值，即SA3>=SI3，则更新第一、第二和第三层子系统间的耦合关系接口数据，激活第三层子系统进行单步长仿真计算，否则跳到(7)，判别程序是否结束；(6)子系统被激活进行单步长仿真，各子系统被激活进行单步长仿真计算时，其计算流程相同，以列车子系统为例，其流程如下 I)激活列车子系统仿真计算；2)获取系统输入数据，系统输入数据是指该子系统上一次仿真计算结果中需要作为本次仿真的输入条件的数据，在第一步仿真计算时，系统输入数据为该子系统的初始化数据；3)获取耦合输入数据，系统耦合数据是指有耦合关系的相关子系统作用于该子系统的数据，作为该子系统的边界输入条件，由有耦合关系的相关子系统在上一次仿真计算时生成，在第一步仿真计算时，耦合输入数据为相关子系统的耦合关系接口数据的初始化数据；4)列车子系统单步长仿真计算，列车子系统利用给定的仿真步长、输入条件和边界条件，进行一次仿真计算，SAl的值增加列车子系统的仿真步长值；5)输出系统单步仿真结果数据，该结果数据包含两部分内容系统中间数据和系统分析数据，系统中间数据是指该子系统本次仿真计算结果中需要作为下一次仿真的输入条件的数据，系统分析数据是指该子系统本次仿真计算结果中用于后处理应用分析的数据；6)生成耦合输出数据，耦合输出数据是作用于有耦合关系的相关子系统的数据，经耦合作用更新相关子系统间的耦合关系接口数据后，作为这些子系统的边界输入条件；7)列车子系统计算等待，本次仿真计算完成后，列车子系统等待下一次仿真激活；(7)判断是否满足仿真结束条件，如果满足条件，则仿真结束退出，否则，返回到第(5)步,继续仿真。采用本专利技术方法，实现了高速列车各子系统间基于仿真步长级的多向耦合仿真，可以满足各子系统间关系复杂的耦合仿真应用需求。同时，本专利技术方法考虑了耦合计算稳定性及精确性保障措施，提高了耦合边界数据的精度，并最终提高整个协同仿真系统耦合计算的精度，将高速列车与周边环境作为一个整体进行研究，更准确地反映出高速列车运行时的动力学性能，提高了仿真研究的准确性。。附图说明图I耦合子系统在耦合器作用下的仿真流程2多层次的耦合仿真步长控制系统描述图图3分布式高速列车系统动力学耦合计算流程及仿真同步控制算法框图。具体实施例方式下面，结合附图，对本专利技术的内容作进一步的详述。在采用图3所示分布式高速列车系统动力学耦合计算流程及仿真同步控制中，耦合器(同图3中的中间件服务器)是联接各子系统实施耦合计算的关键。各子系统在耦合器的驱动下进行仿真的时空同步推进，通过求解高速列车系统动力学耦合关系模型，更新各子系统的仿真计算边界，驱动相关子系统间实施耦合作用，最终实现高速列车-线路-弓网-牵引供电-列车空气动力学子系统间基于积分步长级的耦合计算，达到对高速列车系统动力学进行计算机仿真的目标。耦合器的基本功能包括I)耦合计算流程控制功能可以控制高速列车各耦合计算子系统的启动、等待和 停止;2)耦合计算组织功能可以按用户所提的全耦合或部分耦合计算需求自适应组织相应的计算模型进行耦合计算；3)耦合计算协调功能由于高速列车各耦合子系统计算步长不一，可以按一定策略(如多层次的仿真步长同步协调策略及基于虚拟时间轴的仿真时长同步协调策略等)对各计算子系统的仿真推进过程进行协调，实现耦合计算功能；4)耦合边界数据处理功能在设定的耦合点进行耦合计算时，可以对高速列车各子系统的边界数据进行处理，为各子系统提供本步计算所需的输入数据，同时采用边界缓和过渡算法、分类数据插值算法等，保障耦合计算的稳定性和精确性；5)耦合数据通信功能可以和各耦合子系统进行实时通信，及时接收将各耦合子系统输出的数据，并将处理好的耦合数据及时发送至各耦合子系统以更新其输入。根据高速列车耦合系统动力学计算模型，在耦合器协调下的高速列车系统动力学仿真的一般流程为I)各子系统在一个积分步长内独立求解；2)向耦合器发送与其他子系统相关的计算结果；3)求解各子系统间耦合关系模型及耦合控制算法，并进行耦合数据处理；4)耦合器向各系统发送耦合数据，更新各子系统的计算边界；5)判断耦合系统计算终止条件；6)如果满足条件，系统退出；否则，返回第I)步。高速列车各耦合子系统的仿真计算流程是一致的，主要分为三个阶段计算等待及数据准备、单步积分计算、计算结果本文档来自技高网...

【技术保护点】
高速列车系统动力学耦合仿真方法，在基于接口的多设备参数、多领域间的相互作用关系及控制参数的环境中实现协同仿真基于仿真步长级的多向耦合控制，采用如下的手段，协调仿真步长不一致的各子系统实现耦合控制仿真：(1)根据子系统的仿真步长，将子系统划分为三层，第一层包括列车子系统、线路子系统和弓网子系统，该层的仿真步长记为S=5×10？5s，耦合仿真控制步长记为SI1=5×10？5s，第二层为空气动力学子系统，该层的耦合仿真控制步长记为SI2=1×10？3s，第三层为牵引供电子系统，该层的耦合仿真控制步长记为SI3=1×10？1s，并设置耦合仿真结束条件；(2)为每个层次设置一个累积仿真步长，第一层累积仿真步长记为SA1，第二层累积仿真步长记为SA2，第三层累积仿真步长记为SA3，并初始化SA1=0，SA2=0，SA3=0；(3)初始化各子系统间的耦合关系接口数据；(4)加载各子系统仿真模型，初始化各子系统；(5)各层累积仿真步长值增加第一层仿真步长值，并作如下判别：1)如果第一层累积仿真步长值大于等于第一层耦合仿真控制步长值，即SA1>=SI1，则更新第一层子系统间的耦合关系接口数据，激活第一层子系统进行单步长仿真计算，否则仿真退出并告警；2)如果第二层累积仿真步长值大于等于第二层耦合仿真控制步长值，即SA2>=SI2，则更新第一和第二层子系统间的耦合关系接口数据，激活第二层子系统进行单步长仿真计算，否则跳到(7)，判别程序是否结束；3)如果第三层累积仿真步长值大于等于第三层耦合仿真控制步长值，即SA3>=SI3，则更新第一、第二和第三层子系统间的耦合关系接口数据，激活第三层子系统进行单步长仿真计算，否则跳到(7)，判别程序是否结束；(6)子系统被激活进行单步长仿真，各子系统被激活进行单步长仿真计算时，其计算流程相同，以列车子系统为例，其流程如下：1)激活列车子系统仿真计算；2)获取系统输入数据，系统输入数据是指该子系统上一次仿真计算结果中需要作为本次仿真的输入条件的数据，在第一步仿真计算时，系统输入数据为该子系统的初始化数据；3)获取耦合输入数据，系统耦合数据是指有耦合关系的相关子系统作用于该子系统的数据，作为该子系统的边界输入条件，由有耦合关系的相关子系统在上一次仿真计算时生成，在第一步仿真计算时，耦合输入数据为相关子系统的耦合关系接口数据的初始化数据；4)列车子系统单步长仿真计算，列车子系统利用给定的仿真步长、输入条件和边界条件，进行一次仿真计算，SA1的值增加列车子系统的仿真步长值；5)输出系统单步仿真结果数据，该结果数据包含两部分内容：系统中间数据和系统分析数据，系统中间数据是指该子系统本次仿真计算结果中需要作为下一次仿真的输入条件的数据，系统分析数据是指该子系统本次仿真计算结果中用于后处理应用分析的数据；6)生成耦合输出数据，耦合输出数据是作用于有耦合关系的相关子系统的数据，经耦合器作用更新相关子系统间的耦合关系接口数据后，作为这些子系统的边界输入条件；7)列车子系统计算等待，本次仿真计算完成后，列车子系统等待下一次仿真激活；(7)判断是否满足仿真结束条件，如果满足条件，则仿真结束退出，否则，返回到第(5)步，继续仿真。...

【技术特征摘要】

【专利技术属性】
技术研发人员：邹益胜，丁国富，张卫华，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：