DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 2. 案例分析 以下记录导师陈桂敏教授临床验案三则如下，且均征询患者并且取得患者同意。 案例一：王某，男，26岁，首诊时间：2018年3月5号，症状：结婚两年未育，性生活和谐，平素有心烦，手心发热，口干，小便夜间频数，腰膝酸软，尿道口有灼痛感，阴囊有刺痛感，面色尚可，大便干，舌红苔少有裂纹，脉弦细数，两天前做精液常规检查，精子数：46 × 10^6/L，存活率：17.3%，PR：13.6%，液化时间 > 60 min。(正常范围：精子数：≥ 39 × 10^6/L，存活率：85%~90%，PR（前向运动）：≥32%液化时间：液化时间 ≤ 30 min，完全液化时间 ≤ 60 min)。 西医诊断：弱精症精液不液化 中医诊断：不育症 证型：肾阴亏虚，络脉郁阻 治法：补肾填精，散结通络 方药：知母15 g山萸肉15 g 生地黄15 g 山药30 g 茯苓15 g泽泻15 g丹皮15 g茺蔚子15 g 锁阳10 g补骨脂15 g 菟丝子15 g 女贞子15 g 皂角刺15 g 连翘15 g穿山甲6 g红藤15 g 十四付水煎服 二诊：2018年3月20号，现心烦减轻，夜间小便次数减少，尿道口灼热感，刺痛感症状消失，口干口苦减轻，房事后腰膝酸软减轻，舌红苔白有裂纹，脉弦细，故原方去皂角刺，穿山甲，连翘，加西洋参15 g，麦冬15 g，继续服用14天。 三诊：2018年4月5号，主诉：上诉诸症均有明显减轻，舌质淡红，苔薄白，脉细，故上方加枸杞子15 g，太子参20 g，继续巩固服用7天并嘱咐夫妻在排卵期同房。患者4月25号来电，妻子已怀孕一月余。 病案分析：该患者精子活力低下故辨病为弱精子症，本案例患者心烦，手心发热，口干，皆为阴虚内热之候，虚热上扰心神故心烦，阴虚津液不足，故口干，肾阴亏虚，相火妄行，迫津外泄，故夜间尿频，肾精不足，腰膝失于濡养，故而酸软无力，舌脉也是一派阴虚内热。故本方用六味地黄丸合五子衍宗丸加减，补肾填精，滋阴补虚，皂角刺，连翘，穿山甲，红藤，取其散结通络，清热解毒之功，可促使液化时间缩短；加锁阳补肾阳，益精血，使阴得阳助则泉源不竭，改善精子活力，提高精子质量。西洋参，麦冬，补气养阴，益胃生津，有益于后天之本，枸杞，太子参，补肝肾，益精血，益气健脾，培元固本，起到一个巩固疗效的作用。 案例二：李某某，男，32岁，首诊时间：2017年10月23号，症状：现备孕二胎，平素有射精无力，射精快，怕冷，腰膝酸软，尤其房事时加剧，白天神疲乏力，小便频数，大便不成形，舌质淡，苔白，脉沉弦。一天前精液常规检查，精子数：58 × 10^6/L存活率：13.3%，PR：13.5%，液化时间 < 60 min。(正常范围：精子数：≥39 × 10^6/L，存活率：85%~90%，PR (前向运动)：≥32%液化时间：液化时间 ≤ 30 min，完全液化时间 ≤ 60 min)。 西医诊断：弱精症 中医辨证：不育症 证型：肾阳亏虚 治法：补肾助阳，固精缩尿 方药：淫羊藿10 g 巴戟天10 g 山萸肉15 g 生地黄15 g 山药30 g茯苓15 g泽泻15 g丹皮15 g 肉桂12 g附子12 g覆盆子10 g 韭菜子10 g 蛇床子15 g 枸杞15 g锁阳10 g菟丝子15 g 五味子15 g 益智仁10 g 14付水煎服 二诊：2017年11月8号，现上诉症状均匀明显改善，但怕冷大便稀溏仍存在，舌淡苔白，脉沉细，故原方加淫羊藿15 g，益智仁15 g，续服14剂。 三诊：2017年12月10号，主诉：上诉症状悉除，舌质淡苔白，脉沉，故嘱其停药，服金匮肾气丸两周。2018年1月20号来电，其妻已孕40天。 病案分析：该患者精子活力差，故西医诊断为弱精症，患者阳气不足，温煦乏源，故怕冷，阳气激发，兴奋不足，推动无力，故神疲乏力，腰为肾之府，阳气不足，故而酸软无力，肾阳虚，膀胱气化功能减弱，故而小便频数，舌脉即一派虚寒之征。故用桂附地黄丸合五子衍宗丸加减，桂附地黄丸，取其微微少火以生肾气之意，去桂枝换肉桂引火归元，淫羊藿、巴戟天兴奋性神经，补肾助阳，增加其温肾阳之功。五子衍宗丸补肾益精，使肾气生化有源。覆盆子，菟丝子，五味子，韭菜子，补肝肾且有固精缩尿之功，标本同治，蛇床子入肾经，温补肾阳之功巨，益智仁，暖肾固精缩尿，合全方经旨，温脾止泄摄唾，止泻之功著。后以金匮肾气丸，缓图其本，巩固疗效，以收全功。 案例三：余某，男，32岁，首诊时间：2017年10月24号，症状：右阴囊外侧疼痛一周，伴有勃起不坚，阴囊潮热，小便热赤涩痛，舌质紫红，苔白腻，脉弦细滑。精液常规检查，精子数：23.4 × 10^6/L，存活率：21.7%，PR：16.2%；液化时间 > 60 min。(正常范围：精子数：≥39 × 10^6/L，存活率：85%~90%，PR (前向运动)：≥32%液化时间：液化时间 ≤ 30 min，完全液化时间 ≤ 60 min)。 西医诊断：弱精症精液不液化 中医诊断：不育症 证型：湿热阻络 治法：清热利湿，散结通络 方药：车前子10 g 黄柏10 g 生地10 g 泽泻12 g 木通12 g 甘草梢12 g山慈菇10 g蜈蚣2条 山药15 g 郁金15 g 肉苁蓉10 g 穿山甲6 g 红藤15 g 皂角刺15 g 蒲公英15 g 黄精30 g 共14付水煎服 二诊：2017年11月8号，现上述诸症减轻，仍有勃起不坚，阴囊部刺痛感，小便涩痛不畅，舌质紫红，苔白，脉弦细，故原方加路路通15 g，滑石15 g，续服14剂。 三诊：2017年11月23号，主诉：诸症均有明显改善，舌质淡紫，苔白腻，故原方减黄柏，生地，甘草梢，山慈菇，加白术10 g，砂仁10 g，嘱咐患者再服14剂，后来电症状好转，需继续治疗。 病案分析：该患者精子活力差，故西医诊断为弱精症；液化时间超过60 min，为精液不液化。患者勃起不坚，阴囊潮热，宗筋驰软，不通则痛，故阴囊部刺痛，小便热赤涩痛，舌质紫红，苔白腻，脉弦细滑，证属湿热下注，互结交争，湿热阻络。用龙胆泻肝汤加减，清热利湿，导湿热于小便出，山慈菇，红藤，皂角刺，蒲公英，取其清热解毒，散结消肿排脓之功，蜈蚣，穿山甲，通络之功著，黄精，脾肺肾均补，具有明显改善精液质量的作用，郁金活血行气，清热利湿，肉苁蓉，主五劳七伤，补中，养五脏，强阴，益精气，多子，补精血，润肠燥 [ 3 ] 。路路通，利水通经，祛风通络，滑石，清热利湿通淋，利六腑之涩结，使留恋之湿热之邪尽除之，白术，砂仁，益气健脾，化湿和胃，鼓舞后天运化，使气血生化有源。
system: [SYSTEM]

user: 3. 结论 精液异常中医辨证论治确有奇效，临床多由肾虚与湿热瘀导致。由于精子生长周期是70天左右，所以我们在中医辨证治疗过程中不可操之过急，妄施温燥补益之品，过用寒凉清热之剂，从老师的医案中我们可以看出，要以平润为度，随证加减。除此之外，以肾为本，贯穿整个治疗弱精病的过程，责之于先天，也是我们基础加减化裁之根本。弱精症的饮食生活调理也很重要，饮食有节，忌食肥甘厚味，辛热炙煿之品，起居有常，方能药到病除之。
system: [SYSTEM]

user: 摘 要 目的：考察道德推脱与持续性欺骗行为的关系。方法：采用分层随机抽样法，使用《道德推脱问卷》和《持续性欺骗行为反应倾向问卷》在广州地区抽取两所高等院校共380名学生进行匿名问卷调查。结果：道德推脱与持续性欺骗行为反应倾向呈显著正相关，并能显著预测持续性欺骗行为反应倾向。结论：道德推脱越高的个体越可能发生持续性欺骗反应行为。 关键词 :道德推脱，欺骗行为，关系研究 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 欺骗作为常见的不道德行为之一，它不仅与市场经济发展过程中的市场调节失效、政府监管不到位以及法律执行力度等他律问题相关，还与个人的自律问题存在密切关联。本文即从个人自律的角度，采用问卷调查对持续性欺骗行为的发生、发展进行考察。 自我概念维持理论(self-concept maintenance)认为人们在进行欺骗行为前多半时候会面临两种相冲突的动机：避免惩罚和追求个人私欲，并试图在其中保持平衡，人们的平衡之道就在于如何从小奸小恶中获取个人好处，同时又不改变自己的诚实自我概念(Mazar， Amir， & Dan， 2008)。个体一方面想获得私欲，另一方面又想避免惩罚或维持自身的诚实自我概念。因此，我们认为即获得私欲又避免了惩罚的个体更可能在后期继续做出欺骗行为。欺骗的施行对个体来讲受到惩罚时更多的表现为一种道德惩罚。道德惩罚可以分为外在和内在道德惩罚两种形式(左高山，刘艳，席云鹏，2007)，其中内在道德惩罚也就是个体良知上的不安，表现为自责、内疚、羞愧等，而诸如此类的道德情绪，可以阻断不道德行为的发生和发展(任俊，高肖肖，2011)。Farnese，Tramontano，Fida & Paciello (2011)研究发现，进行欺骗行为的个体必须控制好羞愧、尴尬类具有抑制个体做出与社会行为规范相反行为功能的道德情绪，而控制好该类道德情绪则会有利于激活道德推脱的认知机制促进道德推脱的发生，如此以来，又反过来强化了欺骗行为。道德推脱(moral disengagement)，也译作道德脱离，是Bandura等人(1991)在道德自我的社会认知理论上发展而来，是指个体的一些特定的认知倾向，这些认知倾向包括重新定义自己的行为使其伤害性显得更小、最大程度地减少自己在行为后果中的责任和降低对受伤目标痛苦的认同(Bandura，2002；杨继平，王超，高玲，2010)，主要用于解释当个体做出不道德行为时却不会产生诸如内疚、自责等不良情绪的社会现象。道德自我的社会认知理论指出，大多数人都建立了个人道德行为标准，这些标准起着自我调节作用，可以引导良好的行为和制止不良行为，而与这些标准相违背的行为将导致个体的内疚和自责。然而，自我调节机制这一功能只有在被激活时才起作用。Bandura (1999)指出道德推脱主要在道德自我调节的过程中发挥着重要作用，能影响个体道德自我调节功能进行选择性的激活或失效。在违反个人道德标准时，个体就会通过道德推脱使道德的自我调节功能失效，以此摆脱内疚和自责。人们常说撒一个谎往往需要更多的谎言来去圆这个谎，那么当个体已经跨出欺骗行为的第一步，个体是“继续欺骗”还是“改过自新”停止进行更多欺骗行为，其核心影响因素是什么？前期已有研究发现道德推脱对学术欺骗、诚实互惠的亲社会行为等(刘珊，石人炳，2017；杨继平，王兴超，陆丽君，张力维，2010；Gabbiadini， Riva， Andrighetto， Volpato， & Bushman， 2014)，但并未区分初次欺骗行为与持续性欺骗行为。 基于以上分析，我们假设道德推脱可能是促使个体进行持续性欺骗行为的重要影响因素，并通过问卷调查考察了道德推脱与持续性欺骗行为倾向的关系。
system: [SYSTEM]

user: 2. 研究方法 2.1. 对象 在广州地区采用分层随机抽样法，抽取两所高等院校(一所大专，一所本科)的大一、大二共380名大学生进行问卷调查。回收有效问卷339份，有效回收率89.2%。其中大专生149名，本科生190名；男生106人，女生232人(1人性别信息缺失)；城市和农村户口分别为127和210人(2人户口所在地信息缺失)。 2.2. 工具 2.2.1. 道德推脱问卷 采用我国学者王兴超(2011)年修订的适用于大学生的《道德推脱问卷》。问卷共32个条目，釆用Likert 5点计分法，得分越高表示道德推脱水平越高。该问卷由道德辩护、委婉标签、有利比较、责任转移、责任分散、结果扭曲、非人性化、以及责任归因这八个道德推脱机制组成。修订后的中文版道德推脱问卷信效度良好，其内部一致性系数为0.92，分半信度系数均在0.80以上。总体而言，修订的道德推脱问卷具有良好的信效度。本研究中的Cronbach α系数为0.909。 2.2.2. 持续性欺骗行为反应倾向量表(Response-to-Failure Cheating Scale) 采用Zemack-Rugar等人(2012)编制的持续性欺骗行为反应倾向量表(RTF-Cheating Scale)，量表由三个假设情境组成，每个情境有9个测量项目，内容包括个体在发现已经发生了欺骗行为的情况下可能出现的认知、情绪和行为三种反应倾向，共计27个项目，量表包括“继续欺骗倾向”(维度一)和“改过倾向”(维度二)两个维度。维度一包括9个项目，为对自我控制失败的适应倾向，表现为“破罐子破摔”，如“我认为我已经作弊了也就没必要诚实到底”、“我认为既然我已经在某些问题上作弊了，我可能还是继续作弊”等；维度二包括18个项目，为对自我控制失败的纠正，主要通过内疚、自责、责任感等情绪情感来体现，表现出“改过自新”的倾向，以继续维持或提高自我控制目标，如“我对于初始问题的作弊有强烈的罪恶感”、“我觉得对开始的问题抄了答案是很不负责任的”等。量表釆用Likert7点计分法，从“一点都不可能”到“非常可能”依次记为1到7分，持续性欺骗行为反应倾向的得分等于“继续欺骗倾向维度分”减去“改过倾向维度分”，得分越高表明继续欺骗的反应倾向越高。本研究的量表翻译采用回译方法进行，首先由三名心理学研究生分别对量表进行翻译，并进行讨论确定最终翻译文本，再请一位心理学博士生对文本进行回译，以尽量保证问卷忠于原版。本研究中该量表的一致性系数为0.887。 2.3. 数据采集及处理 由一名心理学研究生做主试，在任课老师的组织协调下，以班级为单位进行集体施测，统一宣读指导语后开始作答，并当场回收所有问卷。所有数据采用SPSS20.0进行统计分析。
system: [SYSTEM]

user: 3. 结果 3.1. 共同方法偏差检验 由于研究采用问卷法，且都是由大学生在同样的教室环境下进行填答，有可能存在共同方法偏差(common method biases)的问题。一方面在进行问卷施测时，采取程序控制，在班级环境下统一测试并当场回收，强调问卷的匿名性及保密性；另一方面，在数据分析时采用Harman单因素检验进行统计控制，结果显示，共有13个因子特征根值大于1，且第一个因子的解释变异量为20.53%，小于40%，因此可以不计共同方法偏差对本研究结果的影响。 3.2. 道德推脱与持续性欺骗行为反应倾向的相关分析 对道德推脱与持续性欺骗行为反应倾向的相关分析表明(见表1)：道德推脱8个维度及其总分均与持续性欺骗行为反应倾向呈显著正相关，说明道德推脱越高的个体越可能发生持续性欺骗反应行为，初步证明了本研究的假设。 3.3. 道德推脱对持续性欺骗行为反应倾向的回归分析 为了进一步考察道德推脱对持续性欺骗行为反应倾向的影响，在控制了其他变量的情况下，采用层次回归分析继续探讨道德推脱对持续性欺骗行为反应倾向的预测作用：第一步将人口统计学变量(学校类型、性别、年级和户口所在地)引入到回归方程，第二步引入道德推脱，结果见表2)，在控制年级和性别变量之后，道德推脱对持续性欺骗行为反应倾向有显著的正向预测作用，对其变异的解释率新增了4.4%，达到7.5%。这一结果进一步表面道德推脱水平越高的个体越可能做出持续性欺骗行为。 Table 1 Correlation analysis of moral disengagement and sustained decei 1 2 3 4 5 6 7 8 9 1.道德辩护 2.委婉标签 0.585 ** 3.有利比较 0.559 ** 0.481 ** 4.责任转移 0.519 ** 0.505 ** 0.674 ** 5.责任分散 0.495 ** 0.552 ** 0.621 ** 0.644 ** 6.结果扭曲 0.453 ** 0.527 ** 0.591 ** 0.537 ** 0594 ** 7.非人性化 0.573 ** 0.538 ** 0.462 ** 0.518 ** 0.523 ** 0.457 ** 8.责备归因 0.526 ** 0.430 ** 0.549 ** 0.511 ** 0.514 ** 0.482 ** 0.515 ** 9.道德推脱总分 0.778 ** 0.758 ** 0.792 ** 0.786 ** 0.801 ** 0.738 ** 0.759 ** 0.740 ** 10.持续性欺骗行为反应倾向 0.224 ** 0.211 ** 0.185 ** 0.244 ** 0.237 ** 0.155 ** 0.201 ** 0.227 ** 0.277 ** 表1. 道德推脱与持续性欺骗行为反应倾向的相关分析 注： * p < 0.05； ** p < 0.01； *** p < 0.001；下同。 Table 2 Regression analysis of moral disengagement and sustained decei 变量 持续欺骗行为反应倾向 第一步 第二步 第一步 年级 −0.134 ** −0.099 性别 0.202 *** 0.163 ** 第二步 道德推脱 0.218 *** F 12.683 *** 14.037 *** R 0.274 0.345 ΔR 0.075 *** 0.044 *** 表2. 道德推脱对持续性欺骗行为反应倾向的回归分析
system: [SYSTEM]

user: 4. 讨论 本研究的调查结果显示，道德推脱各个维度及其总分均与持续的欺骗行为反应倾向呈显著的正相关。进一步回归分析表明道德推脱可以正向预测持续性的欺骗行为反应倾向。这个结果提示我们，道德推脱水平越高的个体越有可能进行持续性的欺骗反应行为。 我们认为主要是因为个体通过道德推脱而对欺骗行为进行认知重建以忽视单个欺骗行为对自己或他人的有害影响，减少自身的负性自我情绪、降低个体原有的道德或目标行为水平，从而增加了个体“继续犯错”的可能性，表现出持续的欺骗行为。换言之，个体在进行了欺骗行为后如果是内疚、羞愧预期占主导，那么个体就会规避这种情绪，产生较高的道德自律水平，以减少欺骗行为的发生。而道德推脱则是一种让个体减少甚至摆脱失控后的内疚、羞愧等自我责备情绪的认知倾向，正因如此，当个体第一次发生某类型的欺骗行为后，通过道德推脱机制的激活和运用，阻碍了个体对该欺骗行为的反思，达到保护相对诚实的自我概念，而通过这种方式所带来的最终结果即为发生持续性的欺骗行为。 本研究的不足之处在于，通过问卷调查法只能反应道德推脱与持续性欺骗行为倾向的相关关系，并对其作出可能性推论，对其间的因果关系及认知神经发生机制尚不明确。因此，后期可结合心理行为实验法及脑科学技术(如ERP技术、FMRI技术等)，从认知加工的时间进程及脑生理基础角度对道德推脱影响持续性欺骗行为的发生、发展机制作进一步分析探讨。
system: [SYSTEM]

user: 1. 引言及主要结果 拓扑是研究几何图形或空间在连续改变形状后还能保持不变的一些性质的一个学科。有了具体的拓扑定义之后，在二十世纪三十年代便有人提出有限集上的拓扑计数问题，即是探究拓扑数量与元素个数之间的关系表达式。但若采用穷举法求取一个n元集合 X n 的拓扑数，其时间复杂度为 O ( 2 2 n ) 。故而，从(n − 1)元集合的拓扑结构性质来递推n元集合的拓扑结构显得事半功倍。 而在近现代对拓扑计数问题的研究成果中，国内有相关的研究： 郭志勇在1991和1992年发表了 [ 1 ] [ 2 ] ，定义了多种特殊的拓扑并进行相关讨论； 何昌和熊明在1995年发表了 [ 3 ] ，给出了一个有限集合上的拓扑数估计公式； 张震于2001年发表的 [ 4 ] ，主要研究了n元有限集上互不同胚拓扑的计数问题； 2006年，洪彩霞在郭志勇的文章基础上对问题进行深入探讨，并发表了 [ 5 ] ； 岳崇山，张贺和景海斌等人在2007年发表的 [ 6 ] ，给出了一个拓扑存在性的必要条件：设τ为一n元集合上的拓扑，则满足 2 n − n + 2 ≤ | τ | < 2 n 的拓扑不存在； 赵婷婷等人在2013年和2014年发表的 [ 7 ] [ 8 ] ，主要讨论了如何计算给定有限集合上的拓扑数的算法； 洪楠坤于2016年发表的 [ 9 ] 中，主要研究了n元有限集上的拓扑数 T ( n ) 和 T 0 拓扑个数 T 0 ( n ) 的一些 性质，得到 T ( n ) = ∑ i = 1 n − 1 ( n i ) T ( i ) + B ( n ) + C ( n ) ，其中， C ( n ) 表示所定义的拓扑在n元有限集上的拓扑数， B ( n ) 表示n元集合上的σ-代数的个数。 基于上述相关文献资料的基础，得到本文研究成果： 针对n元有限集的拓扑数探究，本文通过研究两种变换：α变换和β变换，得到各种分类的情况： { ϕ ， X n − 1 ， X n } ∪ A 1 、 { ϕ ， X n − 1 ， X n } ∪ A 2 、 { ϕ ， X n − 1 ， X n } ∪ A 1 ∪ A 2 、 { ϕ ， X n − 1 ， X n ， { x n } } ∪ A 2 、 { ϕ ， X n − 1 ， X n ， { x n } } ∪ A 1 ∪ A 2 。 再考虑交并运算封闭的情况，排除掉不能构成拓扑的情况有：当 A 1 ≠ { ϕ } 时， { ϕ ， X n − 1 ， X n } ∪ A 2 ， { ϕ ， X n − 1 ， X n ， A 1 ∘ ， A 2 • } ， { ϕ ， X n − 1 ， X n ， { x n } } ∪ A 2 ， { ϕ ， X n ， X n − 1 ， A 1 ∘ ， A 2 • ， { A 1 ∘ ∪ x n } ， { x n } } 均不构成拓扑。 从而对各种情况进行递推计数，得到(表1)。其中，未知量 Q 1 ， Q 2 等代替无法直接求取出计数结果的拓扑情况(所有符号定义见表2)。 Table 1 Counting formula 变换α的总计数公式 2 n − 4 + 2 P n − 1 + 2 T 1 ， n − 2 U p + 1 n − 2 + W 1 ， n − 2 + Q 1 变换β的总计数公式 P n − 1 + 3 ∑ k = 1 n − 2 C n − 1 k P k + 3 ⋅ 2 n − 1 − 5 + 2 T 1 ， n − 2 U p + 1 n − 2 + W 1 ， n − 2 + Q 2 总的计数公式 2 n − 4 + 2 T 1 ， n − 2 U p + 1 n − 2 + W 1 ， n − 2 + Q 1 + P n − 1 + 3 ∑ k = 1 n − 2 C n − 1 k P k + 3 ⋅ 2 n − 1   − 5 + 2 T 1 ， n − 2 U p + 1 n − 2 + W 1 ， n − 2 + Q 2 表1. 计数公式
system: [SYSTEM]

user: 2. 预备知识 本节先给出需要用到的相关定义与符号概念： 定义2.1 设X是一个非空集合，τ为集合X上的一族子集，称τ为X上的一个拓扑，当且仅当： 1) X和空集 ϕ 在τ中； 2) τ中任意个元素的并集仍在τ中； 3) τ中有限个元素的交集仍在τ中。 于是称，集合X和它的拓扑τ共同构成了一个拓扑空间，记为(X，τ)；称τ中的元素为这个空间的开集 [ 10 ] 。 定义2.2 称一个拓扑空间(X，τ)满足 T 0 公理，若对于任何不同两个点 x ， y ∈ X ，存在τ中的开集包含 其中一个，而不包含另一个，此时称τ为X上的 T 0 拓扑。本文用 T 0 ( n ) 表示n元有限集合上的 T 0 拓扑个数，其中 n ≥ 2 。 定义2.3 设X是一个非空集合， Ω 是X的一族子集，我们称 Ω 是X上的一个σ代数，当且仅当： 1) ϕ ∈ Ω ； 2) 若 A ∈ Ω ，则A的补集 X \ A ∈ Ω ； 3) Ω 中任意个元素的并仍属于X。 定义2.4 集合S上的一个分割是指S的一族两两互不相交的非空子集，它们的并是S，且排除掉分割情况为 { ϕ ， S } ；n元有限集合 X n 的分割数是指 X n 上所有可能的分割的个数，记为 X n 。即， X n 表示n元 有限集 { x 1 ， x 2 ， ⋯ ， x n } 的分割成各个子集 C a 1 ， C a 2 ， ⋯ ， C a k ( k ≥ 2 )，使得 ∩ s = 1 k C a s = ϕ 且 ∪ s = 1 k C a s = X n − 1 ，同时 C a s ≠ ϕ 成立的情况种数。 利用贝尔数 B n (Bell Number)直接求得： X n ∗ = B n − 1 = ∑ k = 1 n − 1 ( n − 1 k ) B k = ∑ k = 2 n S ( n ， k ) ，其中， S ( n ， k ) 是第二类Stirling数，其值为把基数为n的集划分为正好k个非空集的方法的数目。 另外，相关符号概念整理如表2所示。
system: [SYSTEM]

user: 3. 变换描述 当n = 1时，拓扑数为1。当 n ≥ 2 时，记 n − 1 元有限集的拓扑个数为 P n − 1 ，设其形为 { ϕ ， X n − 1 } ∪ A 1 ，其中 X n − 1 = { x 1 ， x 2 ， ⋯ ， x n − 1 } ， A 1 ⊆ 2 X n − 1 ，这里 { ϕ ， X n − 1 } ∪ A 1 包含了 { ϕ ， X n − 1 } 的情况，但 X n − 1 ∉ A 1 。 记 X n = X n − 1 ∪ { x n } ，设 A 2 是将 A 1 中的所有元素分别加入 x n 之后所成的集合。 针对n元有限集的拓扑数探究，将n − 1元素拓扑情况递推到n元素拓扑情况，考虑并运算封闭，可进行两种变换，把n元有限集的拓扑情况分类成α变换和β变换的结果(表3)。 Table 2 Implications of correlation symbol 符号 表示含义 相关表达式 A 1 作为n − 1元集的幂集的子集， 形为 { ϕ ， X n − 1 } ∪ A 1 为n − 1元有限集的拓扑 A 1 ⊆ 2 X n − 1 ，其中 X n − 1 = { x 1 ， x 2 ， ⋯ ， x n − 1 } ； 这里 { ϕ ， X n − 1 } ∪ A 1 包含了 { ϕ ， X n − 1 } 的情况， 但 X n − 1 ∉ A 1 A 2 将 A 1 中的所有元素分别加入 x n 之后所成的集合 \ X n ∗ n元有限集 { x 1 ， x 2 ， ⋯ ， x n } 的分割成各个子集 C a 1 ， C a 2 ， ⋯ ， C a k ( k ≥ 2 )，使得 ∩ s = 1 k C a s = ϕ 且 ∪ s = 1 k C a s = X n − 1 ，同时 C a s ≠ ϕ 成立的情况种数 X n ∗ = ∑ k = 2 n S ( n ， k ) C j A 1 中的元素 \ B i A 2 中的元素 B i = C i ∪ { x n } A 1 ∘ A 1 中只有单个子集的情况 A 1 ∘ ⊆ X n − 1 ， 且 A 1 ∘ ≠ ϕ ， A 1 ∘ ≠ X n − 1 A 1 • \ A 1 • = X n − 1 \ A 1 ∘ A 2 • \ A 2 • = X n \ A 1 ∘ C j p ⊕ 将“ A 1 中的元素为：若为有p个元素的集合”，归为一类，记为 C j p ⊕ C j k + 1 ⊕ 比 C j k ⊕ 的元素个数多1 B i q ⊕ 将“ A 2 中的元素为：若为有q个元素的集合”，归为一类，记为 B i q ⊕ B i k + 1 ⊕ 比 B i k ⊕ 的元素个数多1 表2. 相关符号含义 Table 3 Transform classificatio 变换α α-1 { ϕ ， X n − 1 ， X n } ∪ A 1 α-2 { ϕ ， X n − 1 ， X n } ∪ A 2 α-3 { ϕ ， X n − 1 ， X n } ∪ A 1 ∪ A 2 α-4 { ϕ ， X n − 1 ， X n ， { x n } } ∪ A 2 α-5 { ϕ ， X n − 1 ， X n ， { x n } } ∪ A 1 ∪ A 2 变换β β-1 { ϕ ， X n } ∪ A 1 β-2 { ϕ ， X n } ∪ A 2 β-3 { ϕ ， X n ， { x n } } ∪ A 2 β-4 { ϕ ， X n } ∪ A 1 ∪ A 2 β-5 { ϕ ， X n ， { x n } } ∪ A 1 ∪ A 2 表3. 变换分类
system: [SYSTEM]

user: 6. 待研究问题 对于n元有限集上的拓扑数，还有无法直接求取出的计数结果： 1) 对于α变换中的(α-3-3-2)情况，若 C r 1 ， C r 2 ， ⋯ ， C r m ， C j p ⊕ 等各自所含元素个数有相同的情况，或 B z 1 ， B z 2 ， ⋯ ， B z q 等各自所含元素个数有相同的情况； 也即未知量 Q 1 的求取。 2) 对于β变换中的(β-4-3-2)情况，若 C r 1 ， C r 2 ， ⋯ ， C r m ， C j p ⊕ 等各自所含元素个数有相同的情况，或 B z 1 ， B z 2 ， ⋯ ， B z q 等各自所含元素个数有相同的情况； 也即未知量 Q 2 的求取。
system: [SYSTEM]

user: 致 谢 作者衷心感谢华南师范大学数学科学学院的赵浩老师，给予我细心的指导与帮助！
system: [SYSTEM]

user: 摘 要 万有引力定律不含时间，不少人认为是超距作用理论，引力可瞬间传播。狭义相对论指出，光速是物理作用速度上限。科学家一直在寻找测量引力速度方法。2012年11月《科学通报》报道中国科学家用重力固体潮观测数据，在时间域分析，得出引力光速传播的结论。本文从重力固体潮理论、计算、天体位置确定等方面，指出固体潮理论隐含引力速度与光速相同。由固体潮观测值与理论值符合，可得出引力相互作用传播速度与光速相同。本文根据武汉站超导重力仪长序列数据分离出月亮M 2 波和太阳S 2 波相对于理论值的相位差均接近于零，符合固体潮理论的结果，再次得到引力以光速传播结论。由日月位置须作光行差校正的天文常规操作，从逻辑上可得出：引力相互作用传播速度必定与光速相同。 关键词 :地球重力学，引力相互作用传播速度，重力固体潮，光速，超导重力仪 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 牛顿“万有引力定律”(1687年)是17世纪自然科学的伟大成果， 它把地面上物体运动和天体运动的规律统一了起来。牛顿之后三百多年，科学家始终没有测量出作为基本物理常数的引力相互作用传播速度(简称引力速度)。 在引力和电磁这两种长程相互作用中，人们对电磁相互作用规律的了解更深一些。在麦克斯韦电磁理论中，光速 C = 1 / ε μ 。即光速与电磁理论中的真空介电常数和真空磁导率常数之间有简单的关系。所以，光速(电磁波速)的确定，对电磁学理论的建立，有重要作用。引力速度测定，对引力相互作用理论的完善是十分必要的。 引力速度没有测定，加之牛顿“万有引力定律”中不显含时间，不少人认为引力可以瞬间传播。而狭义相对论指出，光速是所有物理作用速度上限 [ 1 ] 。另一些科学家则认为引力作用以光速传播。基本物理常数认识上存在着严重分歧。 对大自然基本常数之一的物体间引力作用传播速度测量无疑是一个亟待解决的重大科学问题。但引力作用极其微弱且无法屏蔽、传播速度又快，依靠现有技术在实验室中测量引力速度没有可能。 罗默(Ole Christensen Roemer)利用木卫食现象首次测出光速。引力速度是否也能利用太阳系中行星、卫星在引力相互作用下的运动效应测量？——如果在运动天体上有能够观测到另一天体对其引力变化的仪器，且有足够高的时间分辨率，引力速度就有测出的可能。 地球物理学为探索地球内部结构和探矿研制出的高精度重力仪，在地球平台上已经能够清晰记录到日月引力的变化，且具有足够高的时间分辨率。 中国科学院地质与地球物理研究所汤克云团队从中国地震局西藏狮泉河站和新疆乌什站PET型弹簧重力仪固体潮观测数据，比较固体潮理论值、引力瞬间传播理论值与观测值相同相位的时间滞后，用时间域资料得出引力以光速传播，相对误差5%的结论 [ 2 ] [ 3 ] 。 如果汤克云这项工作的结论正确，就是中国科学家在国际上首次测得引力速度，无疑是物理学的一个基础性重要结果。 对汤克云团队的研究结论，有支持的、也有反对和怀疑的或不能判断真伪的。为对这项工作的科学性达成共识，中科院组织了两次评议研讨会，有9位院士、11位研究员参会。会上支持和反对的院士都有。比较一致的看法是，问题最终在天体位置(视位置与真位置)的确定。研讨会未能达成共识，反对意见较多，但也不能证明汤的这项工作就是错的。 在维基百科的Gravity词条下，关于Speed of gravity一段，目前只介绍了两项关于测量引力速度的工作：“2012年12月，中国的一个研究团队宣布，他们已经测量了地球潮汐的相位滞后时间，这似乎证明了引力的速度等于光速；……在2017年10月，LIGO和室女座探测器在伽玛射线卫星和光学望远镜两秒钟内接收到来自同一方向的引力波信号。这证实了引力波的速度与光速相同” [ 4 ] 。 被推选为2003年度世界十大科技进展候选项目，由理论物理学家Kopeikin和射电天文学家Fomalont利用木星掩食射电类星体机会，用射电望远镜阵观测类星体位置变化，宣称首次测得引力以光速传播的结果(误差20%) [ 5 ] ，却没有列入维基百科该词条。这说明，中国科学家用测量地球潮汐相位滞后时间测量引力速度的科学思路，在国际上并未被否定。 本文从重力固体潮理论、计算、天体位置如何确定等方面加以述说，如何从台站固体潮观测值与理论值的符合，得出引力相互作用的传播速度与光速相同。 汤克云用乌什和狮泉河站弹簧重力仪一天的观测数据，在时间域中分析得出引力以光速传播的结论。本文根据中国科学院测量与地球物理研究所武汉站观测精度更高的超导重力仪3年连续观测数据，用调和分析方法分离出月亮M 2 波和太阳S 2 波，根据它们各自的观测值相对于理论值的相位差均接近于零，在频率域中分析再次得到引力以光速传播的结论。 本文还从日月位置须作光行差校正这一天文常规操作，从逻辑上指出：引力速度只能与光速相同。 下面先简单介绍日月引力引起地球重力固体潮的计算表达。
system: [SYSTEM]

user: 2. 重力固体潮理论值的计算表达 重力固体潮由日月引力变化造成，测站的重力固体潮理论值由牛顿万有引力定律计算。 万有引力定律表示为如下公式 F = − G M m r m 2 (1) 以月球和地球为例，m为月球质量，M为地球质量，r m 为月心与地心之间的距离，F为月球与地球的相互引力，G为引力常数。如图1所示。 图1. 地球与月球坐标位置示意图。τ，δ m ，Ψ，r m ，m，d m 分别为月球时角、赤纬、月球极角、月心到地心距离、月球的质量、月心到A点的距离；当图1中月球改为太阳时，上述参数用τ s ，δ s ，Ψ s ，r s ，m s ，d s 表示 A ( r ， θ ， L ) 点的月球引力位 V ( A ) 为 [ 6 ] [ 7 ] ： V ( A ) = − G m d m = − G m ∑ 0 ∞ r n r m n + 1 P n ( cos ψ ) (2) 勒让德多项式 P n ( cos ψ ) 是A点的纬度 φ ( φ = π 2 − θ ) 和天体在观测时刻的赤纬 δ m 、时角 τ 和极角 ψ 的函数。 A点的引潮力位 T ( A ) 等于天体对地球上单位质点的引力位 V ( A ) 加上公转时单位质点的惯性离心力位 Q ( A ) ： T ( A ) = V ( A ) + Q ( A ) = − G m r m ∑ 2 ∞ ( r r m ) n P n ( cos ψ ) (3) 重力潮汐 Δ g 是引潮力位在垂直方向的导数： Δ g = − ∂ T ∂ r (4) 在观测点A记录到的重力固体潮是观测时刻日月的质量、时角、赤纬、极角、到地心距离、到A点距离的函数。 计算地球某测站A于某时间段的重力固体潮理论值，只须在计算程序中输入测站的经纬度、高程和所计算时间段的起始和终止时刻。 图2是新疆乌什站PET重力仪2009年7月21~22日的重力观测值经低通滤波滤去地脉动噪声后的 图2. 新疆乌什站PET型重力仪记录的2009年7月21~22日的重力固体潮观测值 重力固体潮。重力潮幅约280 μgal (微伽)，只有地球表面重力值的三百五十万分之一。 理论值(乘以潮汐因子值)与观测值之间的差别很小，两条曲线几乎重叠在一起，经多次局部放大才看出差别。 用重力固体潮观测确定引力速度，测站测到的重力固体潮与日月引力之间是什么关系？
system: [SYSTEM]

user: 3. 重力固体潮与日月引力的关系 按定义，测站A的重力固体潮是日月对A点的引潮力在垂直方向的分量。日月对A点的引潮力则是日月对A点引力与对地心引力之矢量差。 重力固体潮与日月引力之间是什么关系？下面以太阳为例，先看太阳对A点的引潮力和太阳引力的关系。 设想地球绕太阳轨道的偏心率为零，太阳对地球质心的引力就是常量。太阳到测站A的距离d m 因为地球自转而变化，A受到的引力出现周期性变化。地球上单位质量受到太阳引力约0.6 gal (伽)，太阳对A点的引潮力幅度只有约±50 μgal，约为太阳引力万分之一。 可见，太阳对A点的引潮力就是扣除了太阳引力占99.99%的不变(偏心率不为零时为极缓慢变化)部分后，因地球自转而产生的变化部分，此部分引力的垂直分量，就是重力固体潮，日波、半日波是主要成分。 高精度重力仪可以记录到日月引力作用下，一天中双峰双谷、有清晰相位特征的重力潮，记录曲线讯噪比可超过50分贝(图2)。 正是日月引力的这个微小的变化部分，让探寻引力相互作用的传播速度成为可能。 与光(电磁波)不同，引力是一种相互作用，不能脱离相互作用着的两端的源。对于太阳的引力作用，t时刻测站A的重力潮值，由测站A与地心的位置及日月质心的位置决定。下面以太阳为例，太阳的位置该如何确定？
system: [SYSTEM]

user: 4. 太阳视位置需作周年光行差校正到真位置 其实，万有引力定律并非不含时间，时间项出现在r中。距离r由相互作用的两质量体的时空坐标计算，时间是必然要参与到运算中的。下面以太阳位置的确定来说明。 由光学仪器在t时刻观测得到太阳在天球上的位置，称为太阳视位置 X S L ( t ) 。太阳光到达地球需要的时间为 Δ t S ≈ 8 m .311 分钟，约500秒。太阳光发出时刻的真实位置称为太阳真位置。 在天文观测中，天体视位置要作光行差校正得到天体的真位置，是因为光速有限，而地球以约30 km/s的速度绕日公转。周年光行差角约 2 0 ″ • 5 角秒。 假设，一年中太阳黄经变化是匀速的，在 Δ t S ≈ 8 m .311 分钟中，太阳黄经 λ S ( t ) 变化 = ( 360 × 3600 × 8.311 ) / ( 365.256 × 1440 ) ≈ 2 0 ″ • 5 角秒。 可见，太阳的光行差改正 2 0 ″ • 5 角秒是由于太阳视位置与真位置之间的时间差 Δ t S ≈ 8 m .311 分钟引起的。计算A测站t时刻重力固体潮理论值，太阳位置要用经光行差校正后的真位置 X S R ( t − Δ t s ) 的时空值。 太阳视位置 X S L ( t ) 与太阳真位置 X S R ( t − Δ t s ) 的黄经相差 2 0 ″ • 5 角秒，时间相差为 Δ t S ≈ 8 m .311 时分。 月亮光到地球的时间 Δ t M 约1 s .28时秒，在重力仪目前观测精度下，难以用重力固体潮资料，把月亮的真位置和视位置分辨出来。因此可以忽畧1 s .28时秒的时差。月亮的真位置与视位置可以认为重合。
system: [SYSTEM]

user: 6. 光行差校正操作，逻辑上隐含：引力速度只能与光速相同 固体潮理论在确定日月位置时需作光行差改正的操作，这一操作从逻辑上再次认定引力速度与光速相同。 太阳作为一个物质实体，同时发出光和引力。太阳的“光源位置”和“引力源位置”自然是一致、重合的。 因为光速有限，地球又以30 km/s速度运动，我们通过光学观测获得太阳的“光源视位置”需要作光行差改正，校正到“光源真位置”。我们认为校正后太阳的“光源位置”和“引力源位置”仍然是一致、重合的。 我们可能忽视了，太阳“视位置”需要作光行差校正的原因是：光速有限。若光速无限，光行差校正角就等于零。对引力也一样，只有引力速度有限，才有“引行差”校正的必要。 我们可以设想，太阳发射“光子”，也发射“引力子”。若要光行差改正后太阳的“光源位置”与“引力源位置”仍重合，只有引力速度与光速相同才可能。 若光速有限引力速度却无限，光行差改正就会出现太阳的“光源位置”移动了 2 0 ″ • 5 角秒，“引力源位置”却没有移动，光行差改正后太阳的“引力源位置”与“光源位置”就不再重合。 所以，固体潮理论中，天体位置作光行差校正操作，认为校正后的“光源真位置”就是“引力源真位置”，实际上已默认引力速度与光速相同。 总之，若光速有限引力速度却无限，光行差改正就会出现太阳的“光源”位置与“引力源”位置不重合的悖论。 狭义相对论指出，光速是物理作用速度上限。若有某种相互作用速度超过光速，就会出现违反“因果律”等各种逻辑上不自洽的境况。读者不难发现，在我们构建的引力瞬间传播固体潮理论模型中，就存在这类窘境。 结论是：引力相互作用的传播速度只能与光速相同。
system: [SYSTEM]

user: 7. 结论 科学界一直存在对引力速度的不同看法。 文献 [ 11 ] 依据1805年Laplace对月球轨道的稳定性分析，断定引力传播速度至少大于7 × 10 6 c [ 12 ] 。而Flandern于1998年对脉冲双星PSR1534+12的分析给出引力速度的下限值为2 × 10 10 c [ 13 ] 。而狭义相对论指出，光速是所有物理作用速度的上限。 引力相互作用传播速度成为物理学亟待确定的基本常数。 精密重力仪能够观测到日、月引力导致的重力固体潮幅度及相位变化，引力相互作用的传播速度就有可能测量出来。这一点首先为从事重力固体潮观测且关心引力物理的中国地球物理学家意识到。汤克云团队用重力固体潮观测数据，在时间域中分析，得出引力以光速传播的结论。本文在频率域的分析，再次确认引力速度与光速相同。 现任国际地潮中心指导委员会成员的许厚泽院士来信指出：“固体潮理论是在牛顿力学框架下，基于引力相互作用的牛顿万有引力定律建立的。现在所有的精密固体潮观测结果和理论计算结果是完全符合的，也是自恰的。但由于地球对引潮力的响应并非完全弹性，地球介质的滞弹性将导致此相位差有一微小延迟。” 本文的分析指出：固体潮理论隐含了引力相互作用传播速度与光速相同。正如许厚泽院士所指出的：“现在所有的精密固体潮观测结果和理论计算结果是完全符合的”，得出引力速度与光速相同的结论是有大量观测数据为依据的。 2017年10月，LIGO和室女座探测器在伽玛射线卫星和光学望远镜两秒钟内接收到来自同一方向的引力波信号。这证实了引力波的速度与光速相同。 引力波和引力，有联系也有差别。引力波是波，可脱离“源”传播；引力是相互作用，不能脱离“源”存在。在我们周围，接触更多的是引力的作用。当前，天体力学、人造卫星、宇宙飞船都在应用万有引力定律，但引力瞬间传播的牛顿万有引力定律存在“超距作用”和“不满足洛伦兹协变”两大缺陷。我国科学家在国际上首次测得引力传播速度，是物理学的一个基础性重要结果，据此，物理学家就能正确理解牛顿引力公式，克服上述逻辑缺陷，构建更完善的引力理论。
system: [SYSTEM]

user: 摘 要 评估司机能否在合理时间内安全疏散是动车组列车司机室设计中的一项重要内容。本文基于对现有轨道交通领域和公共建筑领域相关安全疏散标准及文献，研究提出了动车组列车安全疏散的评估依据和评估内容；结合通用的人员安全疏散时间判据图，提出了动车组列车所需安全疏散时间构成项及相关合理的取值；最后基于仿真工具Java，以一款司机室的安全疏散仿真分析为例对相关取值的应用进行了介绍。 关键词 :安全疏散，动车组列车，列车司机，评判准则 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 对人员安全疏散的研究源于20世纪初，当时常采用观察描述、访问研究等定性分析的方法进行研究；从20世纪70年代开始，安全疏散问题的研究开始聚焦人的逃生行动能力研究；从20世纪80年代以后，借助实验手段，安全疏散研究开始从动力学及计算机模拟方向进行研究，并建立了多种关于安全疏散行为规律的数学模型；从20世纪90年代初至今，研究开始聚焦到随机行为规律的定量研究，研究人员逐渐提出了各种研究安全疏散仿真的实际模型和手段 [ 1 ] [ 2 ]。 公共场所安全疏散系统的研究和应用相当广泛，在银行、信息、交通、商业等领域均需要进行大量仿真模拟 [ 3 ]。有很多专用的工具可用于对安全疏散系统进行仿真分析，如 GPSS-PC，SimScript，MathLab等 [ 4 ]。可是这些仿真工具普及程度较低，兼容性较差，很多仿真必须借助C语言进行编程，总体难度较大，维护也比较困难 [ 5 ]。Java较上述GPSS-PC，SimScript，MathLab等仿真工具而言，具有系统兼容性强，操作方便，具备跨平台特性的优点 [ 5 ] [ 6 ]。基于Java的这些优越性，本文选择它作为仿真工具来进行动车组列车安全疏散仿真分析，并以3X动车组 [ 7 ] 司机室安全疏散仿真分析为例进行了介绍。
system: [SYSTEM]

user: 2. 列车安全疏散研究和评判准则 2.1. 仿真分析的依据 用于列车安全疏散仿真分析的评判依据包括：《GB50490-2009城市轨道交通技术规范》 [ 8 ] ；《GB50157-2013地铁设计规范》 [ 9 ] ；《SFPE: Handbook of Fire Protection Engineering》 [ 10 ] ；《建筑防火性能化设计》 [ 11 ] ；《BSDD240 Fire Safety Engineering in Buildings》 [ 12 ] ；《TB/T 3091-2008 铁路机车司机职业健康检查规范》 [ 13 ]。 2.2. 仿真分析的内容 司机室是司机操作列车作业的空间，为了保证司机在紧急情况或发生事故时能够安全疏散，《EN 45545-4: Railway applications—Fire protection of railway vehicles—Part 4: Fire safety requirements of railway rolling stock design》 [ 14 ] 规定机车司机室两侧每侧均应设置1处尺寸不小于500 mm高、400 mm宽的紧急逃生窗。因此，在紧急情况下，司机可以有2条安全疏散线路：1) 从司机室两侧窗之一疏散；2) 司机从最近的侧门疏散。 2.3. 安全疏散涉及的标准内容分析 2.3.1. 动车组列车可用安全疏散时间的确定 目前，国内外均无关于动车组列车车内人员疏散的时间规定，因此缺少用于评估可用安全疏散时间(Available Safety Egress Time， ASET) [ 15 ] 的依据。但是，我国《GB50490-2009城市轨道交通技术规范》 [ 8 ] 中7.3.2条规定内容和《GB50157-2013地铁设计规范》 [ 9 ] 中19.1.19条规定内容均提出，在紧急情况发生的情况下，应保证将一列进站列车的预测最大载客量以及站台上的候车乘客在6 min (360 s)内全部撤离至安全区域。鉴于动车组列车运营于地面上，而且列车外面一般是开阔的露天环境，因此将乘客离开动车组列车即视为达到安全区域。基于上述分析，本文以这2个标准文件作为依据，将可用安全疏散时间(ASET)确定为360 s (6 min)。 2.3.2. 人员安全疏散时间判据图 图1是常用的人员安全疏散时间判据图 [ 11 ] [ 15 ] [ 16 ] [ 17 ]。从图中可知，所需安全疏散时间(Required Safety Egress Time， RSET)是指从灾情发生时刻起到人员安全疏散到安全区域的时间，包括灾情探测警报时间t alarm 、预动作时间t pre 和人员疏散移动时间t move ，其中预动作时间又包括觉察时间t reg 和反应时间t resp 两部分，其关系可以用下式表达 [ 11 ] [ 17 ] [ 18 ] [ 19 ] ： t RSET = t alarm + t pre + t move               = t alarm + ( t reg + t resp ) + t move (1) 图1. 人员安全疏散时间判据图 1) 探测警报时间(t alarm ) 探测警报时间包括监测出灾害的时间t det 和发出警报通知人员的时间。探测灾害有两种方法：自动探测系统和人工探测方法。对于由自动探测装置探测灾害的时间主要由装置的性能、灾源的强度以及灾源距离装置的距离决定，探测系统越先进，探测灵敏度越高，探测时间也越短。人工探测方法主要是内部人员通过视觉和嗅觉等感官通道来感知灾害的发生。如火灾的探测，一般靠近灾源近的人员易感觉到，但大多数情况下并不能有效察觉出灾害的发生。如自动检测装置探测火灾，考虑火灾发展20 s之后火源附近温度和烟气浓度已很高，探测装置即可探测到火灾，因此铁路领域通常将探测警报时间取为20 s [ 16 ] ；但文献 [ 17 ] 也提出从安全角度考虑探测警报时间应为60 s (1 min)。 2) 预动作时间(t reg + t resp ) 预动作时间是指在接到报警信息后，采取具体疏散行动之前所需要的时间，包括信息确认时间、行为反应时间 [ 11 ] [ 17 ] [ 18 ] [ 19 ]。预动作时间的长短与选用的报警设备有直接关系，另外还与人员的年龄、行为特征、对建筑物的熟悉程度、人员的灵敏性以及人员的集群特征相关 [ 20 ]。对于客运列车，一般可采用现场广播的方式进行报警，加之不同的功能区域都设有一定数量的乘务人员，所以人员的响应时间通常较短。表1是英国《BSDD240 Fire Safety Engineering in Buildings》 [ 12 ] 中根据统计数据和经验推荐的各种用途建筑内采用不同火灾广播系统时的人员预动作时间。由于车辆上采用的是现场广播系统，而且车辆上通常会配备乘务人员，因此人员的预动作时间可参考表1中选取为60 s (1 min)。 3) 司机疏散运动时间(t move ) 疏散运动时间是指从人员开始做出疏散行动到疏散至安全区域的时间 [ 17 ] [ 18 ] [ 19 ] ，即发生灾情后车内人员从车厢经由通道、出口疏散到车厢外。所有人员安全撤出列车车厢的时间取决通道的顺畅、人员行为、车门宽度等多种因素 [ 16 ] [ 17 ]。人在灾情中的行为十分复杂，不同年龄、性别、对灾情的认知程度对人的逃生行为有很大差别，而这个差别主要体现在疏散速度上 [ 21 ]。《SFPE: Handbook of Fire Protection Engineering》 [ 10 ] 和《建筑防火性能化设计》 [ 11 ] 给出了人员疏散的步行速度，详见表2。 我国的动车组列车均采用单司机驾驶，《TB/T 3091-2008铁路机车司机职业健康检查规范》 [ 13 ] 对我国铁路机车司机的身体健康状态均有明确规定，所有执行驾驶任务的司机均属于身心健康的人员。由于长期从事驾驶任务，因此列车司机对列车的安全疏散线路应比较熟悉和清楚。考虑到列车内的空间较一般建筑的空间狭窄，因此本次分析采用水平疏散中的平均步行速度1 m/s作为司机的疏散平均速度。此外，由于司机从司机室逃生窗口内翻窗出去比较困难，所耗时较长，因此根据经验将逃离司机室侧窗过程的时间设置为60 s。 通常，若人员疏散所需时间(t RSET )小于可用安全疏散时间(t ASET )则视为人员能够安全疏散 [ 21 ]。综合考虑到上述各项因素，将探测报警时间(取60 s)、司机的预动作时间(取60 s)、司机以1 m/s的疏散平均速度，从各逃生出口离开仿真对象的时间总和小于6 min (360 s)视为设计方案满足疏散要求。也即仿真获得的人员疏散所需时间(t RSET )满足下式： t RSET < 36 0 s (2) 也即：t alarm + (t reg + t resp ) + t move < 360 s；其中，t alarm = 60 s，t reg + t resp = 60 s；因此t move < 240 s即为合格。
system: [SYSTEM]

user: 3. 列车安全疏散仿真分析案例 3.1. 仿真模型 用于仿真分析的司机室车体外观如图2所示。司机室车厢内部可分为3个部分：司机室、会议室和VIP室。司机室是司机的作业环境，两侧各设置1个满足安全疏散尺寸的逃生窗口；会议室包含8个座 Table 1 The preparation time of different fire broadcasting syste 建筑物用途 建筑物特性 预动作(响应)时间/min 报警系统类型 W 1 W 2 W 3 办公楼、商业、厂房和学校 建筑物的人员处于清醒状态 熟悉建筑物及报警系统和疏散措施 <1 3 >4 商店、展览馆、博物馆、休闲中心等 建筑物的人员处于清醒状态 不熟悉建筑物及报警系统和疏散措施 <2 3 >6 住宅或寄宿学校 建筑物的人员处于睡眠状态 熟悉建筑物及报警系统和疏散措施 <2 4 >5 旅馆或公寓 建筑物的人员处于睡眠状态 不熟悉建筑物及报警系统和疏散措施 <2 4 >6 医院、疗养院及其他社会公共福利设施 有相对数量的人员需要帮助 <3 5 >6 表1. 不同火灾广播系统时的人员预动作时间 注：W 1 为现场广播，来自消防控制室；W 2 为事先录制好的声音广播系统；W 3 为采用警铃、警笛或其他类型报警装置的报警系统。 Table 2 The walk speed of different safety egress condition 人员特点 平均步行速度(m/s) 水平疏散 楼梯向上疏散 楼梯向下疏散 对建筑的位置、通行线路不熟悉的人员 1.00 0.45 0.60 对建筑的位置、通行线路熟悉且身心健康的人员 1.30 0.58 0.78 不能自力行动的人员 0.80 0.33 0.50 表2. 人员疏散步行速度 图2. 用于仿真分析的司机室车内布置 位以及2张折叠桌，并与司机室相邻；VIP车厢包含16个座位，与会议室间隔一组侧门(左、右各1个侧门)。在不考虑旅客影响的情况下，本文分别对司机从单侧侧窗(安全疏散出口1)和从就近的侧门(安全疏散出口2)疏散的2种情况进行仿真评估分析。 3.2. 司机从司机室两侧窗之一疏散 图3是利用Java对司机从司机室两侧窗之一(由于两侧对称分布，因此分析一侧即可)的左侧逃生窗口疏散的仿真分析(图3中座椅与侧窗间的蓝色线条为安全疏散路径)，通过仿真可以得出，司机从座椅到逃生窗处(平均速度取值1 m/s)的移动时间(t move )约为2 s。 因此，司机从司机室的左侧逃生窗口疏散并离开列车的移动时间为：t move = 2 s + 60 s (司机翻爬侧窗的时间) = 62 s。t move < 240 s，也即表明该司机室侧窗设计符合安全疏散要求。 3.3. 司机从侧门疏散 距离司机座椅的最近的侧门位于会议室和VIP室之间，司机从侧门疏散必须经过会议室。图4是利用Java对司机从最近的侧门疏散的仿真分析(图4中的蓝色线条为安全疏散路径)，通过仿真可以得出，司机从侧门疏散的移动时间(t move )约为11.8 s。 因此，司机从就近的侧门疏散并离开列车的移动时间为：t move = 11.8 s。t move < 240 s，也即表明该司机室就近的侧门设计符合安全疏散要求。 图3. 司机从司机室左侧安全疏散窗疏散时的路径和时间 图4. 司机从就近的侧门疏散时的路径和时间