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

user: 摘 要 通过HFACS模型，总结出非煤地下矿山人因事故的五个一级指标和二十三个二级指标，运用层次分析法对测评指标进行两两比较，描述了各个一级指标和二级指标的所占权重，采用专家打分取平均值的方法构造判断矩阵，用Matlab软件计算矩阵的最大特征值，一致性指标，得出层次分析排序的结果有满意的一致性，权系数的分配是合理的。因判断矩阵依赖于判断专家的知识结构、判断水平、个人喜好等主观因素影响，通过改进的模糊一致性判断矩阵构造方法进行修正，来验证专家的打分方法是否符合实际情况，经分析得到专家给出的优先矩阵有满意的模糊一致性，使可信度和实用性大大增强。 关键词 :HFACS模型，层次分析，判断矩阵，模糊一致性 Copyright © 2017 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. 引言 众所周知，我国是矿业大国，资源丰富，非煤矿山行业生产总值占全国GDP总值的1%左右，是国民经济发展的重要基础。但是，矿难时有发生，严重威胁着矿产资源的开发，各类矿难事故平均每年死亡约3000人，是矿山行业事故高发的国家之一，百万吨矿石死亡率超过美国、南非等矿业发达国家的30倍。各类伤亡事故造成的直、间接经济损失为GDP的1%~2.5%，照此计算，我国非煤矿山平均每年的事故损失约为11~26亿元 [ 1 ] 。据统计，2005年至2014年期间，全国非煤矿山累计发生事故17，919起，死亡21，025人，平均年发生事故1791起，死亡2102人 [ 2 ] ，是我国工业生产中伤亡事故最严重的行业之一，事故总量大、重特大事故比重高、死亡人数多、每年死亡人数仅次于煤矿事故和道路交通事故，在国内各行业中排第三 [ 3 ] 。2013年，地下矿山共发生事故378起、死亡480人，分别占事故总起数和死亡人数的57.4%，56.3% [ 4 ] ，可见地下矿山相对露天矿山危险性更大，伤亡更加严重，因此本文专门针对地下矿山进行研究。另外，通过权威部门近几十年来的统计调查，非煤地下矿山各类事故中，由于人的不安全行为而诱发的事故占到80%以上，其他原因主要指设备原因和自然灾害，这个数字已经成了不争的事实，因此我们当前的主要问题就是把人管好，约束和控制好人的行为，这个问题解决了，安全才有出路，那么如何解决人的问题，首先要深层次挖掘产生人因事故的本质原因，建立分析人因事故模型，进而针对原因展开工作。
system: [SYSTEM]

user: 2. HFACS模型指标体系 人因分析与分类系统(The Human Factors Analysis and Classification System， HFACS)是安全分析时广泛应用的人因分析工具之一，以航空为中心、在航海、铁道、煤矿、医疗等领域均被采用。使用该系统可以得出事故原因中的人为因素，并能从表层行为追溯到深层组织原因，无论在确定事故原因还是制定防范措施都具有重大意义 [ 5 ] 。但HFACS是以航空安全事故为背景设计的人因分析模型(表1)，针对非煤地下矿山而言，不能完全套用，本文首先根据我国非煤地下矿山的实际情况和特点对HFACS框架的详 Table 1 HFACS model index system [ 6 ] [ 7 ] [ 8 ] [9 一级指标 二级指标 具体描述 不安全行为 技能型失误 注意力不集中、记忆错误等 决策失误 判断失误、方案制定错误、采用错误手段等 认知失误 视觉错误、听觉错误、对设备、环境、人员认知错误、规范、区域认知错误等 习惯性违规 习惯性违反规章制度、未安全操作、指令执行等 特殊违规 情绪突变，饮酒、疾病而导致的违规 不利精神状态 精神疲倦、焦虑、丧失情境意识、疲劳作业等 不利生理状态 出现幻觉、方向感偏差、感冒等 不安全监督 不当的劳动组织 不合理的任务安排、人员搭配、工艺设计、生产计划等 未能纠正已知问题 违章操作、设备隐患、危险行为等没有及时纠正 监管违规 监督未按制度执行、未按相关法律监督、监督制度不合理等 安全投入不到位 技术、设备、人员培训、科技攻关等资源的配置和管理不足。 不安全行为的前提 生理与心理缺限 有限的夜间视力、体能不足、反应迟钝、不够迅速的理解和反应能力等 信息传递缺陷 人员沟通与协调不当、信息有误等 个体准备不达标 休息不足、知识缺乏、缺少安全警觉性等 设备缺陷 人机设计不当，防护缺陷，功能分配不当等 不利的环境因素 噪声、粉尘、振动、色彩、照明，气压等 监管不足 无有效监督、缺少监督政策、无监督制度、无监督人员、无监督培训等 组织影响 安全文化缺失 企业文化，安全意识、安全素养、责任心等缺失 安全教育培训不足 员工知识、技能、理念、行为等方面教育不足 危险源辨识不彻底 没有做到横向到边、纵向到底的危险识别 规章制度不落实 规章制度缺乏针对性且不落实 社会环境 经济发展水平 经济发展历史地位决定性作用 法制建设水平 法制建设水平决定着安全意识 表1. HFACS模型指标体系 [ 6 ] [ 7 ] [ 8 ] [ 9 ] 细内容进行全面的改良和补充，以提高HFACS分析的适用性，从而有利于HFACS框架在非煤地下矿山中的应用。
system: [SYSTEM]

user: 3. 模型分析 3.1. 评价因子集的建立 设地下非煤矿山安全行为等级的一级指标评价因子集为： 这里 为评价方面个数，分别为 —不安全行为、 —不安全行为的前提、 —不安全监督、 —组织影响、 —社会环境。 二级指标因子集记为 ，式子中， 为第 个评价方面评价因素的个数。记 为 应分配的权重值， 满足 ，且， 。 为对应 分配的权重值， 满足 ，且 。 层次分析法(Analytic Hierarchy Process， AHP)是萨蒂(Saaty， 1971)提出的一套决策方法，主要应用于不确定情况下及具有多个评估准则的决策问题。它根据各个测评指标的相对重要性来确定权重，通过对测评指标进行两两比较，同时运用萨蒂给出的标度法，使复杂的、无序的定性问题能够进行量化处理。表2为层次分析法对重要程度的划分情况。 表2反应了两个测评指标相对重要程度的得分，例如一个指标 相对于另一个指标 略微重要，则测评指标 相对测评指标 的比较得分为3，反过来测评指标 相对测评指标 的比较得分为1/3。 3.2. 判断矩阵的构造 采用专家打分取平均值的方法构造判断矩阵 ，用Matlab软件计算矩阵 的最大特征值 ，一致性指标为 ，查找相应的随机一致性指标 。对 ，saaty给出了 的值，如表3所示。 ，随机一致性比率 层次分析排序的结果有满意的一致性，权系数的分配是合理的。并且，最大特征值所对应的特征向量为 ， 做归一化处理可得 ，从计算结果可以看出，非煤地下矿山人因事故中，不安全行为和不安全监督占主要原因。作为作业人员，须加强对不安全行为的控制；作为监管部门，需要加大对不安全行为的监管力度。 同理，可以运用层次分析法得出二级综合指标，见表4。 从表4中，我们可以很清晰的看出，不安全行为中的决策失误占有绝对主导地位，不安全监督中的安全投入不到位和监督违规也同样占有主导地位。 从以上方法中可以看出，关键是构造判断矩阵，其一致性一般会受到专家的知识结构、判断水平、个人喜好差异等众多主观因素的影响，再加上人因事故本身的不确定性和模糊性，判断矩阵的一致性有时并不容易达到。下面我们采用改进的模糊一致性判断矩阵来验证专家的打分方法是否符合实际情况，其主要步骤是构造模糊互补矩阵，进一步得到模糊一致性判断矩阵，再通过两个检验的一致性标准进行评估。 决策者通过元素的两两比较重要性，通过0.1~0.9标度构造优先关系矩阵 ，该矩阵满足以下特点： 1) 2) 3) Table 2 Scale scal 相对重要程度 标度值 说明 同样重要 1 两比较方案的贡献程度同等重要 略微重要 3 经验与判断略微偏好某一方 明显重要 5 经验与判断强烈偏好某一方 非常重要 7 非常强烈偏好某一方 极端重要 9 有足够证据偏好某一方 相邻两程度的中间 2、4、6、8 需要折衷时使用 表2. 比例标度表 Table 3 The value of R 1 2 3 4 5 6 7 8 9 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 表3. 的值 Table 4 The weight of each index in the index syste 非 煤 地 下 矿 山 人 因 事 故 一级指标 二级指标 二级综合指标 不安全行为B1 (0.4759) 技能型失误C1 (0.185) 0.0880415 决策失误C2 (0.583) 0.2774497 认知失误C3 (0.075) 0.0356925 习惯性违规C4 (0.087) 0.0414033 特殊违规C5 (0.025) 0.0118975 不利精神状态C6 (0.032) 0.0152288 不利生理状态C7 (0.013) 0.0061867 不安全监督B2 (0.2636) 不当的劳动组织C8 (0.215) 0.056674 未能纠正已知问题C9 (0.055) 0.014498 监督违规C10 (0.358) 0.0943688 安全投入不到位C11 (0.372) 0.0980592 不安全行为的前提B3 (0.0538) 生理与心理缺陷C12 (0.131) 0.0070478 信息传递缺陷C13 (0.132) 0.0071016 个体准备不达标C14 (0.125) 0.006725 设备缺陷C15 (0.106) 0.0057028 不利的环境因素C16 (0.178) 0.0095764 监督不足C17 (0.328) 0.0176464 组织影响B4 (0.0981) 安全文化缺失C18 (0.089) 0.0087309 安全教育培训不足C19 (0.518) 0.0508158 危险源辨识不彻底C20 (0.075) 0.0073575 规章制度不落实C21 (0.318) 0.0311958 社会环境B5 (0.1087) 经济发展水平C22 (0.756) 0.0821772 法制建设水平C23 (0.244) 0.0265228 表4. 指标体系中各指标的权重 计算重要性排序指数 ，求模糊一致性判断矩阵元素为 ，这里的n为判断矩阵的阶数，得到模糊一致性判断矩阵为 ，为了验证专家给出的优先关系矩阵是否合理，我们给出两个一致性检验标准。 令 ，这里 ，通过大量实践，我们认为取 且 较为合理。 利用“和法”求得排序向量 ， ， ，略大于0.1，说明专家给出的优先矩阵有满意的模糊一致性，比较符合实际情况，可以采纳。
system: [SYSTEM]

user: 4. 结论 本文从非煤地下矿山的实际出发，依据HFACS模型，总结出非煤地下矿山人因事故的五个一级指标和二十三个二级指标，通过层次分析法对测评指标进行两两比较，同时运用萨蒂给出的1~9标度法，使复杂的、无序的定性问题能够进行量化处理，描述了各个一级指标和二级指标的所占权重，从中可以更加具体的反正各个指标的影响力度，比较客观的反应了非煤地下矿山人因事故的各种影响指标所占权重。而判断矩阵依赖于判断专家的知识结构、判断水平、个人喜好等主观因素影响，因此本文又通过改进的模糊一致性判断矩阵构造方法进行了修正，来验证专家的打分方法是否符合实际情况，其主要步骤是构造模糊互补矩阵，进一步得到模糊一致性判断矩阵，再通过两个检验的一致性标准进行评估，得出更为合理的评价指标，可信度和实用性大大增强。
system: [SYSTEM]

user: 摘 要 本文研究了特定形状的轴对称旋转体以不同初始条件在滑翔过程中，受马格努斯效应作用的运动特性。本文从理论上推导得到旋转体运动方程并进行数值求解，进一步预判特定初始条件下旋转体将出现特殊运动现象——“斤斗”运动，提取并分析了影响“斤斗”运动出现的参数和它们之间的相互关系。本文还从实验上进行了待定参数的测定与拟合，再现了理论预期的“斤斗”运动。 关键词 :马格努斯效应，轴对称旋转体，“斤斗”运动 Copyright © 2017 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. 引言 高速旋转的物体在空气中飞行时，其运动轨迹会发生弯曲。1852年德国物理学家海因里希·马格努斯(Heinrich Magnus)首次描述了这种空气动力学效应，因此该现象被称作马格努斯效应 [ 1 ] 。马格努斯效应在体育运动中极为常见，比如足球的“香蕉球”、排球的侧旋球、乒乓球的弧圈球，等等；因而对这种现象的研究能够促进相关体育项目的发展，提高运动员的竞技水平 [ 2 ] [ 3 ] [ 4 ] [ 5 ] 。时至今日，马格努斯效应的产生机制已经得到了较为完善的解释并在球类运动中得到广泛应用。另一方面，利用马格努斯效应滞空的飞行器作为一种长期存在的设想，尽管已有模型和实验进行了探究 [ 6 ] [ 7 ] [ 8 ] ，但至今并未实用化，因此对其滞空原理进行研究就成为了一个很有趣并具备实用价值的课题。这一问题也出现在2014年第28届国际青年物理学家锦标赛(28 th International Young Physicists’ Tournament)上。现有的相关文献大多只是对马格努斯飞行器进行模型设想和定性描述，或者仅对飞行器运动方程进行推导，鲜有文献将飞行器理论运动轨迹与实验轨迹进行对比分析；尚未有文献对飞行器在运动过程中出现的“斤斗”现象进行讨论。 本文的研究对象是由两个一次性纸杯底部对底部粘接而成的特定形状轴对称旋转体，如图1(a)所示。将橡皮筋缠绕在该旋转体表面并将其发射出去，则该旋转体能够借助马格努斯效应以特定轨迹在空中滑 图1. (a) 本文研究所用“马格努斯滑翔机”，系两个一次性纸杯底部粘合而成的轴对称旋转体；(b) 滑翔机受力分析图(图中所标物理量均为矢量) 翔一定距离，故称之为马格努斯滑翔机(下简称滑翔机)。本文将对其滑翔机制和运动形态进行较深入的理论探讨和实验验证，主要内容包括：受力分析得出滑翔机的运动方程；数值计算预测其运动轨迹以及对初始条件的依赖关系；实验测定待定参数，与理论预期进行对比；再现“斤斗”运动并分析该运动出现的条件。
system: [SYSTEM]

user: 2. 理论分析 滑翔机在空中飞行期间受到了3个力和1个力矩的作用，分别为马格努斯力F、重力G、空气阻力f以及空气阻力矩M r ，如图1(b)所示。由于马格努斯力与空气阻力的方向随滑翔机运动而改变，因此本文将滑翔机的运动在正交系中分解，分解为随质心的二维平动(可进一步分解为水平方向与竖直方向的直线运动)和绕滑翔机纵向对称轴的转动。 2.1. 平动空气阻力 考虑到空气是一种低黏流体，且通过自由落体运动实验测算，可以认为滑翔机在飞行过程中所受平动空气阻力f与飞行速度平方 成正比，即有 (1) 其中 为常数， 为空气密度， 为过滑翔机纵向对称轴的截面面积(即两个等腰梯形面积之和)。 2.2. 转动角速率 首先，已知牛顿粘性公式 ， (2) 其中， 为滑翔机受到的转动阻力， 为空气的粘滞系数， 为滑翔机运动速度切向分量的大小， 为半径为r处的切向速度梯度， 为滑翔机表面积。 引入一个近似模型来估算滑翔机角速度大小的变化情况。假设所用纸杯的平均半径为 ，由式(2)可求得滑翔机所受到的转动阻力矩大小为 (3) 方向与角速度方向相反。进一步，通过计算可得滑翔机的转动惯量为 ， (4) 其中， 为纸杯的面密度， 为纸杯高度。综上，由 ，可得滑翔机角速率随时间变化率为 . (5) 下面对 的数量级进行估算，以进一步化简公式，并得到可供实验参考的参数。用差分近似微分，可得 ， (6) 代入式(5)并做积分可得 ，(7) 其中， 为初始角速度， 为纸杯高度。 经测量，实验所用纸杯质量为6.6 g，高度为7.2 cm，纸杯底面直径为5.0 cm，开口直径为6.8 cm， 取为底面直径与开口直径的均值；空气粘滞系数取 。将以上数值换算为国际单位制并代入到式(7)进行计算化简，可得滑翔机转动角速率随时间变化关系为 ， (8) 其中常数量纲为[T −1 ]。实验中测得滑翔机滞空时间不超过10 s，由式(8)可知，滑翔机在整个滑翔过程中角速率变化量与初始值之比小于2%。综上可得一个重要结论，即针对本文所研究的系统，可以认为在滞空期间，滑翔机转动角速度大小保持初始值不变，据此本文在后续计算中忽略滑翔机飞行过程中转动角速度大小的变化。 2.3. 马格努斯力 在滑翔机垂直于对称纵轴方向上取任一截面，如图2所示。当滑翔机旋转时，由于粘性效应，其表面附流层气流将随之运动，流速大小可表示为 ，其中A是与滑翔机表面性质、内摩擦、粘滞性有关的常数，取值范围[0，1]。 用v表示滑翔机质心平动速度，u表示滑翔机表面转动切向线速度(即表面附流层气流切向流速)，由图3中所示的几何关系可得，滑翔机横截面某直径与表面圆周交点a、b处的合速度分别为 图2. 滑翔机垂直于对称纵轴的横截面图 图3. 组成滑翔机的单个纸杯沿对称轴的纵剖面图 (9) 设a， b两点的压强差为 ，由伯努利方程可得 (10) 其中， ， 分别为a，b点离地面的高度， 为所取截面半径。由图3得 ( 为纸杯的最大半径)，则 ；同时注意到 ，将力的微分 对面积 进行积分可得 ， (11) 即单位线元所受马格努斯力。由于a，b两点对称，因此积分范围为[0，π]。将 分别在 轴和 轴投影，可得 (12) 其中， ， 分别为滑翔机飞行速度的水平分量和竖直分量。由于组成滑翔机的两个纸杯对称，则可只对滑翔机一侧(即一个纸杯)进行受力分析，其所受马格努斯力的x轴分量和y轴分量分别为 (13) 其中， ， 为纸杯高度。 2.4. 滑翔机运动方程 综上，可到滑翔机的运动方程 (14) 其中， ，为平动空气阻力的水平分量； ，为平动空气阻力的竖直分量。由于该方程组没有解析解，下面将采用数学软件Matlab®编程求数值解，并将数值模拟结果与同参数同初始条件下的实验结果做比较。 相比于现有文献，本文理论分析在以下几个方面有所创新。首先，本文讨论了滑翔机在空气中运动时角速度的变化并得出了滑翔机在滞空期间其角速度的变化可以忽略的结论；其次，在受力分析中考虑了滑翔机自身重力的影响使得最终拟合结果更加精确；最后，待定参数A根据实验结果拟合确定(详见下文)，比现有文献 [ 2 ] [ 6 ] 直接假定A的取值更为合理。
system: [SYSTEM]

user: 3. 数值模拟、参数测定与实验对比 实验采用高速摄像机拍摄滑翔机的运动过程(摄影机已经通过砖墙测试以保证所用焦段没有镜头畸变)。将滑翔机两端开口用保鲜膜绷紧封住，在滑翔机正对摄影机的圆形侧面上用记号笔标出一条直径，以此作为滑翔机旋转角测量参考线。实验中，利用投影仪将直角坐标投影到墙上(如图4所示)，测得每一个投影坐标网格的边长为4.875 cm。以橡皮筋完全脱离滑翔机瞬间为初始时刻，该时刻滑翔机所在位置为坐标原点，侧面参考线与水平直线夹角为初始旋转角。实验中先后用480 fps、240 fps、120 fps等不同帧数的摄影模式进行拍摄，发现在480 fps模式下，参考线最为清晰，故采用此模式。将拍摄所得滑翔机空中运动视频导入图像处理软件Photoshop ® 进行逐帧处理，选取相邻的易于读取旋转角的画面进行读数，通过像素点与坐标网格长度的关系可计算得到滑翔机当前时刻的坐标，通过参考线可读取滑翔机当前时刻的旋转角；将坐标和旋转角对两帧画面之间的时间间隔进行差分即可求得相应时刻的平动速率和转动角速率。 实验还对式(1)中的待定参数c进行了测定，并对式(10)中的待定参数A进行拟合。c可通过滑翔机自由落体运动测得。实验中将滑翔机从一定高度无初速度释放并标记滑翔机下落高度为1 m的位置，对其整个下落过程进行拍摄，通过视频得到滑翔机下落1 m所用时间，代入运动方程 (15) 可解得c = 2.37。参数A则可采用拟合法来获得。将实验测得的初始值代入(14)式，A从0到1以0.01为间隔取值，在Matlab ® 中利用四阶龙格库塔法编程计算得出滑翔机的理论运动轨迹，与相同初始条件的实验运动轨迹进行比较，通过两条曲线的拟合程度来确定 的最终取值。结果显示，当A = 0.35时，两条曲线拟合度最好，如图5(a)所示。为了验证所得A值是否正确，图5(b)和(c)分别显示了相同参数值下水平位移与竖直位移随时间变化的实验和数值结果对比，从图中可以看出，实验与数值计算结果拟合得很好。 利用同一个滑翔机进行三组实验，通过改变橡皮筋缠绕圈数、松紧程度以及发射角度来改变初始条件。将实验结果与数值模拟结果(待定参数取值c = 2.37，A = 0.35)进行对比，如图6所示。可以看出，实验与数值结果符合得很好，有效地证明了理论分析和待定参数测定的正确性。 图4. 实验中的滑翔机发射装置图以及用于测量滑翔机运动坐标和速度的投影坐标网格 图5. 用拟合法确定参数A。实验中初始坐标为(0.45， 0.2)，单位：m，初始水平速度与竖直速度分别为5.2和2.5，单位：m/s。初始角速度为480/7，单位：rad/s 图6. 不同初始条件下实验结果和数值模拟结果对比。初始条件参数如图注所示
system: [SYSTEM]

user: 4. 滑翔机“斤斗”运动分析 数值模拟结果显示，在一定条件下，滑翔机会进行“斤斗”运动，即滑翔机的运动轨迹会出现回环。由式(13)可知，在忽略重力的情况下，马格努斯力的方向与速度方向垂直。数值计算可以验证重力项对马格努斯力的影响极小，数值小于速度项的5%，由此可认为马格努斯力的方向与速度方向垂直。马格努斯力使得速度方向朝该力的方向转动，速度方向的转动进而引起马格努斯力的同向转动，可类比洛伦兹力对点电荷的作用。如图1(b)所示，如果沿-z方向观察，可以发现马格努斯力与速度将同时逆时针旋转并保持相互垂直状态，因而一定条件下速度的水平分量会出现负值，即出现“斤斗”现象。改变实验的初始条件上可得到理论上预期的轨迹，如图7所示。 “斤斗”运动的出现与A的取值、初始时刻速度竖直分量与水平分量的比值以及角速率三个因素有关。A的取值与滑翔机的表面性质以及粘滞性有关，对于同一滑翔机以及同样的实验环境，可以认为是不变量。下面从改变初始时刻竖直速度与水平速度的比值以及初始角速率两个方面进行分析。通过控制变量法分别改变速度比值和角速率，利用Matlab ® 数值计算不同初始条件下水平速度的变化，如图8所 图7. 滑翔机出现“斤斗”运动的实验和理论对比。实验初始条件参数如图注所示 图8. “斤斗”运动出现条件研究。(a) 竖直速度与角速度保持不变，水平速度由小到大变化；(b) 水平速度与竖直速度保持不变，角速度由小到大变化。(v1表示水平速度，v2表示竖直速度) 示。当水平速度出现小于0的情况时，滑翔机即出现“斤斗”运动。 在图8(a)中，水平速度由6 m/s变化到42 m/s，水平速度与竖直速度的比值逐渐增大，可以看出滑翔机运动状态由“有斤斗区”向“无斤斗区”逐渐过渡，且角速度越大过渡周期越长。在8(b)中，角速度由55 π/s变化到100 π/s，可以看出，随着角速度的增大，滑翔机运动状态由“无斤斗区”向“有斤斗区”逐渐过渡。 由上述分析可得以下结论：1) 角速度大小、初始竖直速度与水平速度之比均与是否出现“斤斗”运动相关，在一个量保持不变的情况下，总可以通过改变另一个量来达到出现“斤斗”的效果；2) 在角速度增大的情况下，“斤斗”运动出现所需的竖直速度与水平速度的比值阈值减小；3) 在角速度不变的情况下，竖直速度与水平速度的比值越大越容易出现“斤斗”运动。
system: [SYSTEM]

user: 5. 总结 本文对马格努斯滑翔机进行了动力学分析，得到含待定参数c和A的运动方程，通过实验测定了待定参数c，以数值模拟与实验对比拟合的方法确定了参数A，并验证了理论分析的正确性；此外，本文还对滑翔机出现“斤斗”运动的条件进行了较为深入的分析。本文较为完整地研究了某一特定形状轴对称旋转体在空气中运动的马格努斯力效应问题，从理论推导、数值模拟和实验角度进行了较为完整的分析，得到了互相一致的结论。
system: [SYSTEM]

user: 摘 要 改善中国服务出口的地区市场结构对促进中国服务贸易的发展具有重要作用，而引力模型在该领域的研究中具有重要的意义和前景。本文从引力模型在服务贸易领域的适用性、修正和扩展引力模型的解释变量、改进引力模型的估计方法、特定因素对服务贸易的影响等方面对国内外的相关文献进行了系统的梳理和评述，以把握其目前动态，从而为后期的相关研究提供参考和借鉴。 关键词 :引力模型，服务贸易，文献综述 Copyright © 2017 by author 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. 引言 由于参与国际贸易的服务业具有高技术含量、高附加值、环保节能等特点，服务贸易已越来越成为各国特别是大国之间竞争的焦点。服务贸易作为我国经济发展的重点，虽已取得了快速的发展，但仍存在一些突出的问题，一是服务贸易长期处于逆差状态，特别是2008年金融危机之后，逆差还呈现快速增长的态势；二是服务贸易占我国对外贸易总额的比重偏低，长期以来只有10%，远低于20%的世界平均水平，另外，服务贸易出口占世界市场的份额也较低，只有4%左右。这些都表明中国服务贸易的发展还相对滞后。 针对中国服务贸易发展相对滞后的原因，一些学者从贸易管理体制、服务业自身发展、国外服务贸易经验借鉴等方面做出定性解释并给出应对思路，典型的有刘雪琴和梅新育(2004) [ 1 ] 、马凌远(2011) [ 2 ] 、李翀(2015) [ 3 ] 、赵景峰和陈策(2006) [ 4 ] 、朱惊萍(2010) [ 5 ] 、徐晓旭(2009) [ 6 ] 等。更多学者尝试从服务贸易竞争力出发，对中国服务贸易发展滞后进行定量解释，相关实证研究主要基于两种思路展开：第一种是依据各类竞争力指数对中国服务出口竞争力进行测度并从中判断出提升中国服务竞争力的发展方向和重点，代表性的研究有陈宪和殷凤(2008) [ 7 ] 、桑百川等(2014) [ 8 ] 、牟岚(2014) [ 9 ] 等。第二种是依据各种理论模型对影响中国金融服务贸易竞争力的因素及程度进行识别和测度，并为促进中国金融服务出口寻找突破口，代表性的研究有陈虹和章国荣(2010) [ 10 ] 、庄惠明等(2009) [ 11 ] 、宋加强和王强(2014) [ 12 ] 等。 从现有研究看，针对中国服务贸易发展滞后问题，很少有学者从改善中国服务贸易地区结构的视角展开分析 1 ，也就是说，中国的服务贸易对哪些国家有潜力？对哪些国家没有潜力？哪些行业有出口潜力？哪些行业没有出口潜力？什么因素决定了服务贸易的出口潜力？这些都是迫切需要研究的问题。而自20世纪60年代以来，引力模型在国际贸易的经验研究中获得了相当成功的运用，引力模型已被开发成为测算贸易潜力、鉴别贸易集团效果、分析贸易模式以及估计贸易壁垒边界成本的有力工具(盛斌和廖明中，2004) [ 13 ] ，但运用引力模型对双边服务贸易流量进行解释和分析的相关研究仍然相当有限，因此本文拟对该领域的国内外研究文献进行梳理，从而为后期的相关研究提供经验借鉴和参考。
system: [SYSTEM]

user: 4. 简单评论及展望 国内外学者关于服务贸易引力模型的研究虽然取得了很多成果，但在如下方面还存在不足：1) 在服务贸易的研究对象上，由于服务贸易统计数据的缺陷，国内外的大多研究主要以BOP统计数据为基础，缺乏对非常重要的商业存在模式服务贸易的研究，而根据WTO的统计，商业存在服务贸易要占到整个服务贸易的一半左右，忽略了它对整个服务贸易的研究就不够完整。2) 国内很多文献研究了某些特定因素对服务贸易的影响，而忽略了多因素的综合分析，比如重要的相对距离、贸易管制等因素，因此我们无法判断这些因素对中国服务贸易的影响，另外也容易产生遗漏变量偏差。3) 国内已有不少学者运用引力模型对中国服务贸易总体的出口潜力进行研究，但细分到具体行业具体模式的研究还比较少见。另外，在服务贸易壁垒的常见度量方法中，数量工具能比频度工具和价格工具提供更丰富的信息量，然而国内学者运用引力模型对中国不同贸易模式不同行业服务贸易面临的非关税壁垒进行数量研究的还非常少见。4) 在引力模型的估计方法上，国内学者虽然引进了一些国外的先进估计方法，但大多研究仍局限于运用截面数据对引力模型进行回归和分析，而忽视了“时间”对双边服务贸易流量的重要影响；另外，国内研究一般只基于单一方法(主要是OLS估计)对引力模型进行评估分析，而不注重多种模型评估方法的比较运用，当然也就无从判断不同的模型评估方法对现实情况的适应程度；还有，对服务贸易引力模型估计非常重要的内生性问题，国内大多研究都没有处理，从而导致估计结果的偏差。 针对服务贸易引力模型研究中的不足，在后期的相关研究中，这些方面是值得关注和加强的，比如关于商业存在模式服务贸易的出口潜力、影响因素、面临的贸易壁垒等问题，再比如运用引力模型对不同贸易模式不同行业服务贸易面临的非关税壁垒进行数量研究，还有通过各种方法对引力模型估计中的内生性问题进行处理等。
system: [SYSTEM]

user: 摘 要 网络舆情的传播主体在网络舆情的发生、发展和演化过程中具有重要作用，舆情的应对也离不开网络传播主体。网络意见领袖作为活跃主体，与网络舆情的演化有着密切的联系。方法层面有必要开展基于网络舆情传播主体的分类研究，掌握活跃主体的形成及在舆情传播中的作用；技术层面有必要基于信息分析方法研究网络意见领袖的发现，证明社会网络分析技术能有效用于舆情主体的挖掘过程及其量化分析。实证展示舆情主体的挖掘过程及其信息分析方法的应用。 关键词 :网络舆情传播，社会网络分析，信息技术，舆情分析 Copyright © 2017 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. 网络舆情传播主体及其分类 网络舆情传播主体是指所有能够连接网络，并可以自主发表意见、表达情感、参与主题讨论以及转发有价值观点的自然人与实体机构。网络舆情传播主体既可以是独立的自然人个体，如一位普通网民、专家学者、政府官员、专职媒体人等；也可以是一个组织机构，如一个新闻网站、官方组织、网络公关公司、军队医院等。 传播主体可分为领舆、白舆和蓝舆三大类。其中，“领舆”是指某个特定主题的权威代表，也就是舆情领袖，其对民众舆情有着特殊的、深刻的影响，民众乐于接受领舆的观点和意见，对大众有较强的影响力，可以引导民意；“白舆”是指舆情民众中的白领，其有良好的教育背景，比较稳定的立场和观点，他们关注政治和社会的热点问题，表达的立场与观点能对普通大众产生一定的影响；“蓝舆”是指舆情民众中的蓝领、网民中的大多数，他们不太关心与已无关的事件，相对缺乏主见，容易受周围人的影响 [ 1 ] 。根据网民表现形式，可将网民细分为潜水型、转发型、附和型、争论型以及领袖型等 [ 2 ] ，不同类型的网民在网络舆情传播过程中扮演了不同角色，承担着各自的作用。其中，潜水型网民表现为沉默式浏览网络内容，一般情况下不显露自己的情绪，不发布自已的观点与意见，他们是网络舆情的围观者，能增加网络的点击量和人气，但不是舆情的传播者，沉默是其最突出的特征；转发型网民以信息浏览为主，发现自身兴趣或认知程度高度一致的内容时，应用网络平台功能对其进行转贴、转发和共享等操作，该类网民的特征是对信息的选择性转发和不直接表态；附和型网民往往没有自己的主张，习惯用最简单语言或图形发表自己的观点和情绪，一般无建设性内容的言论；争论型网民思维敏捷，不满足作为旁观者或附和者，对特定事物有自己见解，并善于表达自己的情绪、态度和观点。潜水型网民不是网络舆情的传播者，但因为人数众多，容易聚合气场，围观形成力量势能；转发型网民承担网络舆情传播的作用，但对舆情本身的归转作用影响有限；附和型网民对舆情传播有一定的推动功能；争论型网民是网络舆情的主要发动者和推动者。
system: [SYSTEM]

user: 3. 基于信息分析方法的网络意见领袖的发现 网络如此浩瀚，意见领袖也多如牛毛，如何寻找和发现网络意见领袖，对应对网络舆情有现实指导作用。发现网络意见领袖的方法有以下4种：① 观察法，在网络新闻、博客、微博、论坛和微信等舆情传播内容载体中积极参与，经简单统计排名靠前的用户，该方法人工作业量大，精确度不高，一般在数据量不大情况下使用。② 统计学方法，根据舆情传播特点，对用户设置一定的参数，如在线时间、粉丝数量、关注度、发贴数、转贴数、跟贴数等，并对参数加权处理后综合排名，得分高者为意见领袖。③聚类分析法，可应用聚类方法对多维属性聚类可以发现意见领袖，但该方法对数据量较大的样本算法效率受到限制 [ 8 ] 。④ 社会网络分析法，下面将给予具体阐述。 3.1. 社会网络分析方法及其应用 实践表明，社会网络分析技术是一种有效的发现网络意见领袖手段 [ 9 ] 。社会网络是指社会行动者及其关系的集合。社会网络是由某些特定群体(人、企业和组织)间的社会关系构成相对稳定的关系网，如果用“点”表示社会行动者，“线”表示行动者之间的关系，社会网络就是由多个点和线组成的集合。社会网络分析技术就是研究“点”和“线”的关系，被广泛应用于社会学的量化研究中。日常研究中，人们关注行动者自身属性的研究，而社会网络分析技术认为社会不是由群体组成的，而是由网络关系组织的。社会网络中行动者之间关系的重要性大于行动者本身属性的重要性。社会网络分析技术就是把“关系”看成是分析单位，利用量化语言对社会结构中“行动者”和“关系”两个核心的结构进行描述。社会网络分析技术包括点和线连成的社群图以及由“社会行动者”与各个行动者之间“关系”形成的矩阵，通过对社群图和矩阵进行分析发现其中的规律。社会网络分析技术用点的度数、距离、直径和密度等参数来量化描述社群图中的社会关系。从“关系”的角度出发 [ 10 ] ，用点度中心度、中间中心性、接近中心性，以及群聚系数等指标来定量研究权力。 网络意见领袖在网络中的权力可以用下述指标进行量化 [ 11 ] ：① 意见领袖在舆情信息传播中的核心作用。一方面由此为中心向外界传播大量的舆情信息，意见领袖的点出度高；另一方面表现出多人对他的关注，即点入度高，点度中心度是其量化指标。② 意见领袖在网络中的控制作用强、权威高。意见领袖在不同的主体、群体和朋友圈之间传递舆情，在舆情传播过程中起到“桥”的作用，中间中心度是其量化指标。③ 意见领袖连接较多的群体。结构洞 [ 12 ] 加快网络信息传播的功能，一个节点控制的结构洞越多说明这个节点连接的子群越多。所以，可以通过社会网络分析技术 [ 13 ] 的点入度、点出度、距离、密度、点度中心度、中间中心性和结构洞等指标来研究和发现网络意见领袖。 3.2. 舆情主体的挖掘过程及其技术分析 以XT产妇死亡事件中天涯论坛的数据来阐述舆情主体的挖掘过程。该事件在“天涯杂谈”板块主题发贴和回贴最多，跟贴最多的贴子为2014年8月15日“XLYC”发起的“XT产妇死亡事件-小护士自己的见解-所有人都理智”主题，截止到2014年9月14日，该贴点击20438人次，发贴与回复721人次，本贴中共有161位发贴人，建立关系407次，参与人中仅有1次互动关系为262人次，参与人中有2次互动关系为70人次。大于2次的有48人次，占总系统的18.43%。将161名参与者进行编号，通过发贴与回复建立参与人之间的非对称关系矩阵，应用UCINET6.0工具软件进行计算，分析其社会网络总体密度值为0.0332，表示人群之间总体互动不高，这个社会网络关系分布符合幂律分布和无标度网络特征，也就是说在这个社会网络中只有少数舆情主体是网络图中的活跃分子，在网络中发挥关键作用。有5人的入度与出度均在前10名内，说明这5人在网络中活跃程度越高，其观点受到广泛关注，是舆情传播中的重要人物。 按时间顺序分时段来观察参与者行为来进行动态分析，互动频率大致可将事件分为：发生至18日、19~21日、22~25日以及26日以后4个时间段，形成4个分时社会网络关系图。4个时间段形成的社会网络密度均小于0.1，社会网络间的互动程度较低，但后期的密度值呈增强态势，表明随时间推移，参与主体逐渐减少，但是主体间的关联紧密程度却得到加强。各时段社会网络的平均直径在2~3之间，表示各主体仅需2~3个主体既可建立联系，说明舆情在整个网络中容易传播。网络的群聚系数总体偏低，提示网络中各主体的联系并不紧密。 通过对社会网络中度、密度、直径、中心性和聚类系数等指标的量化分析与可视化展示，科学阐述了网络舆情传播特征和演化特点，直观把握了舆情传播主体间的关系，动态揭示了舆情主体演变过程，同时还准确量化了舆情传播过程中的关键人物，进而从漫漫网海中发掘网络舆情演化过程中最为核心和活跃的舆情主体，挖掘意见领袖，对科学合理地引导和应对网络舆情有积极意义。 利用社会网络分析方法对网络舆情中的新闻媒体关系网络进行研究，对两个舆论场之间关系进行研究，能挖掘两个舆论场网络新闻媒体关系、演变规律和核心传播主体。运用Ucinet6.0软件中核心/边缘模型的CORE算法计算，判断核心成员在核心/边缘模型中的位置，得出各网络媒体的核心度。最终发现对舆情传播起关键作用的媒体有以下特点：① 以传统的主流媒体和大型商业网站为主，传统的主流媒体网络版媒体有中国青年网、新华网、人民网、央视网和凤凰网等，大型商业网站有搜狐、新浪、凤凰网、网易、腾讯等。② 传统的主流媒体以发布原创和新闻为主，而大型商业网站网络新闻则以转载其它网站新闻为主。③ 新的网络新闻媒体作用不可忽视，如新京报以企业化的经营架构，因其新闻速度快、原创新闻多、视角独特、剖析深入、言语犀利等特点，在舆情传播中表现非常活跃，对舆情的发展有很强的推动作用。
system: [SYSTEM]

user: 4. 小结 网络舆情的传播主体在网络舆情的发生、发展和演化过程中具有重要作用，舆情的应对也离不开网络传播主体。网络意见领袖和网络舆情的爆发、归转和演化有着密切的联系。意见领袖既能参与舆情事件，快速形成网络舆论压力，引发网络舆情爆发和动荡；同时也能设置议题，影响和引导舆论发展方向。实证表明卫生领域有官员、专家学者和大量医务工作者活跃在网络中，成为在卫生行业很有影响力的网络意见领袖。在舆情应对中，要善于挖掘意见领袖，正面发挥舆情领袖的舆情引导作用，大力培养网络评论员。而从不断发展变化的网络世界发现意见领袖，网络分析技术是一种有效挖掘和发现意见领袖的方法。
system: [SYSTEM]

user: 摘 要 用ADAMS软件建立旋转弹丸的仿真模型，基于翻转力矩矢量垂直于阻力面的特性给出了惯性坐标系中翻转力矩矢量分量的表达式，通过FORTRAN子程序实现了翻转力矩的作用。然后用仿真模型对各种转速和各种翻转力矩幅值情况下的弹丸陀螺稳定性进行了仿真实验，研究结果对弹丸的稳定性设计和研究具有参考价值。 关键词 :旋转弹丸，陀螺稳定性，ADAMS，仿真 Copyright © 2017 by author 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. 引言 弹丸的稳定飞行具有极其重要意义，要研究这个问题需要建立6自由度的弹丸运动方程组 [ 1 ] - [ 6 ] ，并进行数值求解，过程复杂而不直观。ADAMS是著名的机械系统动力学仿真软件，可以方便快捷地建立复杂机械系统的动力学模型，直观清晰地展示机械系统的动力学行为，在各行各业获得日益广泛的应用。本文希望将这一先进的动力学仿真工具用于弹丸陀螺稳定性的仿真研究，为弹丸飞行力学研究提供一种新的途径。旋转弹丸的模型建立比较简单，但弹丸上的翻转力矩的施加是一个难点，因为翻转力矩垂直于弹轴和弹丸质心速度矢量的平面，即阻力面，由于弹丸进动，阻力面不断旋转，因此，翻转力矩的方向是不断变化的，既不固定在弹丸上，又不固定在空间某一点。本文基于矢量叉乘公式得出了在惯性坐标系中翻转力矩矢量的分量表达式，通过FORTRAN语言编程，实现了翻转力矩的施加。最后应用ADAMS2010建立了旋转弹丸的仿真模型，进行了多工况仿真实验。研究结果对弹丸的稳定性设计和基于ADAMS的弹道仿真具有一定的参考价值。
system: [SYSTEM]

user: 2. 翻转力矩的施加 翻转力矩幅值的表达式如下 [ 1 ] [ 2 ] [ 7 ] ： (1) 其中ρ为空气密度，v为弹丸飞行速度，S为弹丸最大横断面积，l为弹丸特征长度， 为翻转力矩系数导数，δ为攻角。在攻角较小时(保持弹丸陀螺稳定性的攻角即为小攻角)， 只与v有关。在小攻角时(1)式可以简化为： (2) 显然上式中的k v 在小攻角时只与v有关。 在图1所示的弹丸模型中，翻转力矩方向的单位矢量表达式为： (3) 其中是 是翻转力矩方向的单位矢量， 是弹轴上的单位矢量， 是弹丸质心速度方向的单位矢量。坐标系oxyz是惯性坐标系， 、 构成的平面即为阻力面， 、 、 在惯性系中的分量分别是： 、 、 ，则有 (4) 图1. 弹丸模型 D 1 、D 2 是弹轴上两点，其在惯性系中的坐标分别是 、 ，则 (5) 其中 是两点间的距离。 在仅考虑弹丸的陀螺稳定性时，可以认为速度矢量的方向不变，即 是常矢量。在图1所示的弹丸上，质心速度矢量与惯性系 轴平行，则 ，这样(4)式简化为： (6) 根据(2)、(6)式可以在弹丸上施加有三个分量的翻转力矩矢量，由于用到弹轴上两点D 1 、D 2 的位置坐标，且弹轴的方位是不断变化的，需要调用ADAMS的位置坐标提取函数。通过FORTRAN子程序结合ADAMS的自定义函数，应用(2)、(6)式实现翻转力矩的施加。
system: [SYSTEM]

user: 3. 翻转力矩幅值为定值时的仿真 为了研究方便，先考虑翻转力矩幅值为定值时的情况，且在本文的仿真实验中，弹轴与弹丸质心速度矢量的夹角(即攻角)的初始值均设定为10˚。弹丸转速为10转/秒，翻转力矩幅值为0.3牛顿米，弹尖上一点(D2)在垂直弹丸质心速度矢量的平面内的轨迹如图2所示，攻角变化如图3所示。 图2表明弹丸做进动和章动合成运动，图3表明弹丸攻角不断减小，体现了弹丸运动的陀螺稳定性。对于一个给定的转速，通过大量仿真计算可以确定常数翻转力矩的稳定阈值，超过该阈值，弹丸就会失去旋转稳定性。为了揭示弹丸攻角的变化规律，将弹丸运动过程的仿真时间延长至250秒，如图4所示，显然攻角最终达到稳定值，弹丸转化为规则进动。 弹丸转速为定值10转/秒时，翻转力矩幅值分别取0.1牛顿米~0.5牛顿米的5个定值进行5次仿真，得到的弹丸进动周期和章动周期随翻转力矩幅值的变化曲线如图5所示。翻转力矩幅值为定值0.3牛顿米时，弹丸进动周期和章动周期随弹丸转速的变化曲线如图6所示。
system: [SYSTEM]

user: 4. 翻转力矩幅值随时间变化时的仿真 实际上翻转力矩随攻角和弹丸飞行速度变化，在研究弹丸陀螺稳定性时，可以认为在一个时间区间 图2. 弹尖点在平面内的轨迹 图3. 攻角变化曲线 内，弹丸的飞行速度变化较小，因此(2)式中的k v 可视为常数，这样翻转力矩只是攻角的线性函数： (7) 为了便于与上面的仿真结果做对比，将攻角为10˚时的翻转力矩幅值设定为0.1牛顿米，则翻转力矩幅值的表达式为： (8) 图4. 攻角变化曲线 图5. 进动与章动周期随力矩变化 弹丸转速为10转/秒(r/s)，弹尖点运动轨迹如图7 (50 s仿真)，攻角变化如图8。
system: [SYSTEM]

user: 5. 结论 1) 只要弹丸保持旋转稳定，翻转力矩幅值不变与变的情况下，弹丸最终都能达到规则进动状态(图4、图8)。 2) 弹丸转速一定时，弹丸进动周期随翻转力矩幅值增大而减小，章动周期基本不变(图5)。 图6. 进动与章动周期随转速变化 图7. 弹尖点在平面内的轨迹 3) 翻转力矩幅值一定时，弹丸进动周期随转速增大而增大，章动周期随转速增大而少许减小(图6)。 4) 弹丸的实际自转角可以高达100转/秒左右，本文仿真的弹丸转速最高为50转/秒，二者有较大差距，这样做的原因一方面是出于应用ADAMS进行弹丸陀螺稳定性的规律性和可行性探索，另一方面则是出于仿真效率的考虑，弹丸每转一周，ADAMS积分器至少需要积分50次才能保证仿真精度。 5) ADAMS可以实现弹丸的陀螺稳定性仿真，可以进行弹丸转速、翻转力矩和稳定性之间的关系研究。可以将上述弹丸陀螺稳定性的研究方法推广应用于弹丸追随稳定性研究，从而将ADAMS软件应用 图8. 攻角变化曲线 于弹丸飞行的全弹道仿真研究，为外弹道学研究开辟一条新途径。
system: [SYSTEM]

user: 摘 要 目的：观察不同年龄段纹眼线患者眼表状况。方法：随机选取眼科门诊纹眼线女性就诊患者60例(120只眼)，年龄24~65岁。行眼表疾病指数(ocular surface disease index， OSDI)问卷调查，裂隙灯检查，格瓦拉–眼表综合分析仪Keratograph 5M检查，根据OSDI干眼分级标准将患者分为干眼组(28例56眼)和正常组(32例64眼)，分别对两组的眼表和睑板腺状况进行评估与比较。结果：OSDI值：干眼组(27.02 ± 14.35)分，正常组OSDI值(7.29 ± 3.36)分；干眼组明显高于正常组(P < 0.001)；NIBUT：干眼组5.93 ± 3.07正常组：13.23 ± 4.89；泪河高度：干眼组(0.18 ± 0.04) mm，正常组(0.23 ± 0.03) mm，这两项指标在干眼组均明显低于正常组(P < 0.001)；角膜荧光素染色评分：干眼组(3.51 ± 1.67)分，正常组(1.23 ± 2.32)分；睑板腺缺失评分：干眼组(3.35 ± 1.25)分，正常组(0.65 ± 0.59)分，这两项指标在干眼组明显高于正常组(P < 0.001)；纳入的60例纹眼线患者的角膜荧光素染色评分与睑板腺的缺失评分呈正相关(rs = 0.838， P = 0.000)。结论：纹眼线患者干眼发病率为46.67%，睑板腺功能障碍是纹眼线患者发生干眼的重要病因之一，Keratograph 5M非侵入式眼表综合分析仪可无创性客观评价患者眼表状况。 关键词 :纹眼线，眼表综合分析仪，眼表，干眼，睑板腺 Copyright © 2017 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. 引言 纹眼线是用文刺器械将颜料用针刺入睫毛根部皮肤内，成外观连续、着色较深的线条，起到扩大、改变眼型，增加睫毛浓密感的作用，给人一种眼轮廓清晰有神的美感，其实质是一种创伤性皮肤着色术。这项曾经风靡一时的美容技术已开展二十余年，其间因纹眼线所造成的眼部损害现象屡见不鲜，如结膜炎、角膜炎、睑缘炎、麦粒肿、霰粒肿及干眼等眼病 [ 1 ] 。本研究通过眼表疾病指数(OSDI)问卷调查，裂隙灯检查，格瓦拉–眼表综合分析仪(Keratograph 5M)检查，评估纹眼线患者眼表状况。探讨Keratograph 5M非侵入性眼表综合分析仪对纹眼线患者的临床诊疗价值，为纹眼线后患者临床诊断干眼及评估睑板腺的形态改变提供参考依据。
system: [SYSTEM]

user: 2. 对象和方法 2.1. 对象 收集2015年8月至2016年8月在西安市第四医院眼科门诊就诊纹眼线的女性患者60例(120眼)，年龄24~65岁(40.73 ± 8.27岁)；纹眼线时间1~30年(13.87 ± 5.23年)。对所有研究对象进行OSDI问卷调查、非侵入性眼表综合分析仪Keratograph 5M检查及裂隙灯检查。根据OSDI评分将研究对象分为正常组和干眼组。所有检查均由同一名有经验的检查者实施。纳入标准：女性纹眼线患者。排除标准：病史，3个月内眼部无活动性炎症，局部未使用任何滴眼液；无角膜接触镜配戴史，无眼部激光或其他眼部手术操作史，无眼外伤病史，角膜缘干细胞异常或其他原因引起的眼表疾病如热烧伤或化学伤等。无其他影响泪液分泌的全身性疾病。 2.2. 方法 2.2.1. OSDI问卷调查 研究对象由同一名医师完成OSDI问卷调查 [ 2 ] ，完成时间为4~6 min。OSDI评分 = 所有得分总和 × 25/测评题目总数，总分为0~100分 [ 3 ] 。OSDI值为0~12分为正常；>12~22分为轻度干眼；>22~32分为中度干眼；>32~100分为重度干眼，根据OSDI评分将患者分为正常组和干眼组。 2.2.2. 非侵入性眼表综合分析仪Keratograph 5M检查 包括NIBUT、非侵入性泪河高度、睑板腺拍摄 [ 4 ] 。各项检查均先测右眼，再测左眼，连续测量三次，取平均值纳入统计。分别翻开研究对象上、下眼睑，用Keratograph 5M摄取睑板腺图像，观察睑板腺形态及缺失情况。并进行评分 [ 5 ] ：0分：无睑板腺的丢失；1分：丢失面积 < 睑板腺总面积的1/3；2分：丢失面积 ≥ 1/3~2/3；3分：丢失面积>2/3。采用上下眼睑的评分总和来计算每眼0~6分的评分。 2.2.3. 裂隙灯显微镜检查 包括睑板腺开口状况、睑板腺脂质分泌物性状、角膜荧光素染色评分 [ 6 ] [ 7 ] 。在裂隙灯显微镜下，嘱患者向上看，使用棉签拨开下眼睑，评价下睑中1/3区域，观察每条睑板腺的开口状况并进行分级：0级为正常，无开口阻塞，睑缘覆盖一层薄而平滑的油脂；1级为1~2个睑板腺开口闭塞或分泌物阻塞；2级为3~4个睑板腺开口阻塞，睑缘油脂较厚；3级为几乎一半睑板腺开口阻塞或狭窄；4级为超过一半睑板腺开口阻塞或狭窄且伴大量黏性分泌物。在裂隙灯显微镜下，嘱患者向上看，使用大拇指中等力度按压下眼睑中1/3区域，半定量评价睑板腺脂质分泌物性状并进行评分，共0~24分。0分：正常，睑缘覆盖清亮、透明睑脂；1分：睑脂污浊；2分：睑脂污浊伴碎屑；3分：睑脂稠厚呈牙膏状。使用荧光素染色观察角结膜染色情况，每眼角膜荧光素染色评分为角膜四个象限的总和，共0~12分。0分：无着染；1分：散在的轻微点状着染；2分：介于1~3分之间；3分：融合或丝状的严重着染。 2.3. 统计学处理 应用SPSS 19.0统计学软件对数据进行统计学分析，数据用均数 ± 标准差表示。使用Mann-Whitney U检验对各指标进行比较；对组间数据比较进行hapiro-Wilk检验，根据检测结果决定使用参数或者非参数检验；使用Spearman相关性分析评估两组中各种因素的相关性，两组间的数据评估使用卡方检验。所有的置信区间被设定为95%；P < 0.05为差异有统计学意义。
system: [SYSTEM]

user: 3. 结果 3.1. 干眼检出情况 本研究共纳入60例近纹眼线女性患者，年龄24~65 (40.73 ± 8.27岁)。OSDI检测出正常人群32例64眼(53.33%)，轻度干眼者8例16眼(13.33%)，中度干眼者15例30眼(25.0%)，重度干眼者5例10眼(8.33%)。根据OSDI干眼标准共检出干眼患者28例(46.67%)。将28例56眼干眼患者纳入干眼组，32例64眼正常人纳入正常组。 3.2. 干眼组与正常组的眼表参数比较 干眼组泪河高度为(0.18 ± 0.04) mm，正常组为(0.23 ± 0.03) mm，干眼组明显低于正常组(P < 0.001)。干眼组和正常组的OSD1值分别为(27.02 ± 14.35)分、(7.29 ± 3.36)分，角膜荧光素染色评分分别为(3.51 ± 1.67)分、(1.23 ± 2.32)分，这两项指标在干眼组均明显高于正常组(P < 0.001)，NIBUT小于正常组(P < 0.001) (见表1)。 3.3. 干眼组与正常组的睑板腺功能指标的比较 相对于正常组，干眼组睑板腺开口状况评分、睑板腺脂质分泌物性状评分和睑板腺缺失评分均明显增加(均为P < 0.0001，见表2)。 3.4. 干眼组与正常组的眼表参数的相关性分析 Spearman相关分析显示干眼组中角膜荧光素染色评分与睑板腺的缺失评分呈正相关(rs = 0.838， P = 0.000)。