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

user: 2. 范冲的史学成就及其原因 范冲(1067~1141)字元长，一字益谦，是范祖禹长子。范冲继承了家族的史学渊源，在史学上成就显著，成为了范氏家族第三代史学典范。作为家中长子，范冲受到父辈的学术影响较多，史学素养很高，并机缘巧合地得以传承父辈的史学事业，取得了令人瞩目的成就。 2.1. 范冲的史学成就 作为“三范修史”之一的范冲，因重修《神宗实录》而闻名于世，其史学成就毋庸置疑。 范冲在重修《神宗实录》时，最值得一提的是其在明显的政治立场影响下，依旧保持了严谨的态度。在重修实录的过程中，考虑到《神宗实录》几经删改，议论不一，恐他日无所质证，还编成《考异》一书，“明示去取，旧文以墨书，删去者以黄书，新修者以朱书”，即为“朱墨书” [ 4 ]。范冲的修史态度是严肃认真的，虽然删修的《神宗实录》偏私旧党，和其父范祖禹“大意止是尽书王安石过失” [ 5 ] 同一思路，但是他保留了墨本和朱本的旧文，使三次编修的《神宗实录》都得以保存，实属严谨之举。这样既明确了修史者的意图和历史责任，又可以使读者鉴识，他日可有所对比判断。 范冲参与修撰的《仙源庆系属籍总要》是范冲重要的史学成就，但是前人较少提及。宋朝的统治者十分重视宗室谱牒，专门设置宗正寺“掌奉宗庙、诸陵荐享，司宗室之籍” [ 6 ]。范冲作为宗正少卿，“掌叙宗派属籍，以别昭穆而定其亲疏”，由于南渡之时，宗正寺旧掌之书皆散失。绍兴四年(1134年)十月十日，范冲就编纂谱牒的重叠问题上书言：“渡江之后，散失簿籍，无凭照据重叠。欲乞自今后如遇训名，本寺撰讫，从本寺行下大宗司勘会，如有重叠，别行改撰” [ 6 ]。这一请求得到了肯定。由此可见，范冲治史严谨，尽力避免所编写的宗室谱牒出现重叠。绍兴五年(1135年)，范冲等编类修纂《仙源庆系属籍总要》，合仙源类谱、宗藩庆系录、属籍三者为一，竭尽全力做到“无愧于昔” [ 5 ]。 绍兴五年四月二十六日，范冲以《仙源庆系属籍总要》的名称和具体框架上书，将之前的三种谱牒合而为一，这是编纂体例上的创新，详细记载三祖之下宗室、宗妇、宗女姓名及其生卒、官爵、赐谥等各项内容 [ 6 ]。 在修纂的过程中，遇到困难及时请求解决，五月七日，范冲等言：“今来将已降指挥宗正司取索报到宗室等家状，与本寺丞孙纬收到从来宗支等文字一处参照，编类修纂……其使不足，申乞接续支降” [ 6 ]。范冲与宗正寺的同僚将所得到的资料认真核对，确保无遗漏，无错误。工作也严格按照规定进行，各种开支和人员配备也一应安排妥当。范冲参与编纂的《仙源庆系属籍总要》较为详尽地记载了宗室的系统，为后世修史提供了资料来源。 除此之外，范冲也奉命为司马光的《记闻》编类，《记闻》是司马光为编写《资治通鉴后记》所作的长编式的资料汇编，正如李焘所说的“今世所传《记闻》及《日记》并《朔记》，皆《(资治通鉴)后记》之具也”，在书写上对“所见所闻所传闻之异必兼存以求是” [ 7 ]。编类《记闻》无疑也为范冲重修《神宗实录》提供了史料来源，在这个过程中他也学习了司马光首创的编写历史的方法——长编法和考异法，并在具体的修撰中加以实践。 除了奉命编写史书之外，范冲独自撰写的史书还有《范太史遗事》1卷、《范祖禹家传》8卷、《宰相拜罢录》24卷以及《要语》1卷，可惜的是都已经亡佚，在其他史书中零星可见记载。 范冲在延续范氏家族史学的过程中发挥了承前启后的作用。其子范仲熊继承其史学传统，参与修史，著有《北记》，此书虽佚，但是其内容在《三朝北盟会编》中有所保留。范氏后学中有任馆阁职务者，侄子范仲艺秉承家风，历任著作郎、国史院编修官、实录院检讨官、实录院同修撰等官职，成为华阳范氏第四代史学典范。范仲黼历任著作佐郎、著作郎、秘书等职 [ 8 ]，不仅在理学上造诣颇高，也继承了家族的治史传统。 范冲能够重修实录，参与修纂皇族谱牒，说明他的史学才能得到了皇帝和时人的一致认可。奉命编修官方史书和编类司马光《记闻》的经历为范冲私修史书提供了史料来源以及方法指导，家族所传承的治史方法也在这个过程中得到了实践和进一步的丰富完善。 2.2. 范冲选择史学道路的原因 成都范氏家族是史学世家，范冲作为家中长子，受到父辈的影响更深，具有继承父辈事业的责任和使命。 父辈的交游圈对范冲也有着深刻的影响，因为曾祖父范镇和父亲范祖禹与司马光的交游，范冲的学术选择受到司马光的影响很深。范镇“平生与司马光相得甚欢，议论如出一口，且约生则互为传，死则作铭” [ 4 ]。从范镇开始，范氏家族就与司马光处在师友交游圈内了。司马光很早便已经认识范祖禹，曾说“自祖禹年未二十为举人时，臣已识之”，他十分欣赏范祖禹，认为其“智识明敏，而性行温良”，且“好学能文”，世所罕见，自己“远所不及”，在神宗面前举荐范祖禹，“使之供职秘省” [ 9 ]，祖禹得以入局与其同修《资治通鉴》，在学问道德上受到司马光指导，两人之间有频繁的书信往来，范祖禹曾感念司马光恩情：“教诲成就，义兼父师” [ 10 ]。在生活中司马光将范祖禹视为家人，“事无大小必与公议，至于家事，公休(司马光之子司马康)亦不自专，问于公而后行” [ 3 ]。 范祖禹与司马光在政见上也是一致的，他们都反对王安石新法。在元祐年间第一次修纂《神宗实录》时，范祖禹就“尽书王安石之过，以明神宗之圣” [ 4 ]，因此被王安石女婿蔡卞记恨，最终被贬死于岭表。在政见上，范祖禹和司马光又与理学家不谋而合，程颐曾经抨击王安石“心术不正，为害最大”，“坏了天下人心术” [ 5 ]。 理学家与司马光都严守礼法，并且有着同样的政治立场，理学圈子对司马光也敬重有加。出于对理学的敬仰，以及司马光对理学的偏向，范冲也与理学有着主动的接触。范冲是程颐的学生，曾经与程颐晚年高足张绎论道：“先生(程颐)昔在洛中，晚坐，张思叔诵‘逝者如斯夫’，范元长曰：‘此即是道体无穷’思叔曰：‘如是说，便不好。’先生曰：‘道须涵泳，方有自得’” [ 11 ]。可见范冲在理学上受到了程颐的指导，对道的看法也得到了程颐的肯定。 范冲肩负着司马光、父亲和理学家共同的愿景和期望——打倒王安石。高宗将重修重任交于范冲：“两朝大典，皆为奸臣所坏，故以属卿” [ 4 ]。心系旧党的范冲得到了机会，通过重修《神宗实录》给予王安石一个盖棺定论，完全否定王安石变法，认为“惟是直书安石之罪，则神宗成功圣德，焕然明白”。面对高宗问及有人重提行王安石法度，范冲的回答毫不留情：“今之背君父之恩，投拜而为盗贼者，皆合于安石之意，此所谓害天下人心术” [ 5 ]。事实上，范冲通过重修实录很好地完成了批判王安石这一任务。 范冲作为家中长子，在学术道路上更多地受到父辈的影响，加之自身较高的史学素养，机缘巧合地得以继承父辈的史学事业，并取得了令人瞩目的成就。与此同时，因为在礼法上一致的认识，以及史学上的师承关系，范冲的父辈与司马光建立起深厚的交情。再加上司马光倾向于理学，范冲热心钻研理学，范冲对司马光的钦佩之意更甚。长期的耳濡目染以及家族与司马光的接触，使得范冲受到了司马光史学的深刻影响。
system: [SYSTEM]

user: 3. 范温的文学成就及其原因 范温(约1070~1125)在文学上颇有建树，其代表作《潜溪诗眼》是宋代著名的诗学评论著作，能够体现范温的文学素养，也对后世产生了重要影响。范温能够在文学上取得如此突出的成就，主要是由于其从幼年开始就处在苏轼门人的交游圈中，以及自身的文学才能。 3.1. 范温的文学成就 范温继承了家族的文学渊源，并且深受苏轼、黄庭坚等文人影响，文学成就主要体现在诗学批评著作《潜溪诗眼》上，这本书由一条一条的论诗条目组成，对人和诗不作详细阐释，论诗尊杜重法，代表了江西诗派初期的诗学观点。范温也把自己的诗论见解、审美标准置于诗本身的记述之中，语言浅显易懂，通过《潜溪诗眼》可以更为清晰地知道当时的诗论宗旨。书中主要包括韵论和诗眼论，钱钟书先生最先注意到范温在韵论上的贡献，“吾国首拈‘韵’以通论书画诗文者，北宋范温其人也。因书画之‘韵’推及诗文之‘韵’，洋洋千数百言，匪特为‘神韵说’之挈纲要领，抑且为由画‘韵’而及诗‘韵’之转捩进阶” [ 12 ]。对范温有着颇高的评价。范温对韵论进行了发展，他总结了前人对韵的各种定义，并在分析的基础上进行了否定，提出了自己对韵的理解：“有余意之谓韵” [ 13 ]，同时将其包含范围扩大到诗、乐、书、画等多种领域。范温的韵学思想对后世也有着重要的影响。 《潜溪诗眼》一书已佚，但是郭绍虞先生整理出了辑佚本，可供学界研究。这本书深受黄庭坚的影响，“诗眼”二字就与黄庭坚诗学“字字有来处”的字句锤炼相契合。黄庭坚对这个门生也是甚为满意的，称赞“范侯年少百夫雄，言行一一无可柬” [ 14 ]。认为他的言行合乎君子之道。范温对于老师黄庭坚的诗论主张也十分推崇，《潜溪诗眼》“所述大抵多本师说，条绪井然” [ 15 ]。 范温在文学上的成就以及家族的文学传统对其后辈同样有着深刻的影响。在史学上颇有建树的范仲艺，其文名更盛。陈亮曾称赞他“文章议论为时宗工” [ 16 ]，而当时诗界著名人物杨万里也时常与之“握手论文”，对他的诗“口诵心藏而不能忘也” [ 17 ]，折服于他在诗学上的才能。他也与“其文擅于当世” [ 4 ] 的陈傅良交流颇多，两人经常一起谈论诗词文章。范仲艺赴任时，陈傅良以所赋之诗相送，他们之间不乏书信诗词往来。 [ 18 ] 同辈的范仲圭被时人王夷简赞赏为“吾蜀大儒，抱负醇正，学问鸿博，根本六经”，在文学上深有造诣，“诗于文章最为难事”，范仲圭“尤工于诗”，能够继承先辈遗风，“视其家三大老(范镇、范百禄、范祖禹)为无愧” [ 19 ]。 范温重视文学也影响了南宋诗坛。范温作为黄庭坚的传人，承袭了黄庭坚的风格，其成书于北宋时期的诗学评论著作《潜溪诗眼》对江西诗派产生了重要影响。范温是吕本中表叔，在吕本中的《紫薇诗话》中有所提及，吕本中在诗学上的成就在某些方面受范温的影响。范温颇为看重句法，认为“句法之学，自是一家工夫” [ 13 ]，诗歌字词的不同安排，能让诗歌展现出不同的韵味。范温在论句法时尤为推崇黄庭坚以及杜甫，以杜甫为祖，以黄庭坚为宗是江西诗派诗学批评的特点。吕本中在论句法时，也与范温一样多次提及黄庭坚和杜甫，认为“前人文章各自一种句法”，这与范温所主张的“句法之学，自是一家工夫”异曲同工 [ 13 ]。 由此可见，范温在诗学上深受黄庭坚影响，其诗学评论著作《潜溪诗眼》具有重要的文学价值，对南宋诗坛产生了深远影响。 3.2. 范温选择文学道路的原因 范温作为范祖禹的儿子，虽然没有选择史学道路，但是却继承了文学的血脉。其父范祖禹与苏轼关系匪浅，与苏门中人过从甚密。因为父亲范祖禹的关系，范温从小就处在与苏轼门人的师友交游圈中，受到苏轼门人的悉心指导。范祖禹曾“招李方叔(廌)教其子温辈” [ 20 ]，李廌是“苏门六君子”之一，其文深为苏轼赞赏。范温在其学术还未定型时，面对以文学称誉乡里的李廌，钦佩之情就会转化为行动，在文学上多下功夫，李廌对其的指导很可能为他在学术上选择文学奠定了基础。 范温在其诗学批评著作《潜溪诗眼》中多次提及苏轼，直接说到在苏轼的讲授下，“晓句法不当重叠” [ 13 ]。文中两次称赞“老坡作文，工于命意，必超然独立于众人之上” [ 13 ]。范温的学术受到了苏轼命意观的影响，在书中时常分析诗词的立意。苏轼在开封的时候，“每谒范纯夫，子孙环绕，投纸笔求作字” [ 2 ]，苏轼来到范祖禹的家中拜访，面对文学大家，范温和其他范氏子弟一样，对苏轼有着仰慕之情，自然不会错过讨教问询的好机会。苏轼曾在写给范温岳父秦观的信中提到“梦得之有子为不死也” [ 2 ]。对范温的欣赏溢于言表。 李石曾经说“苏范世好，兄弟雍睦” [ 21 ]，两家是同乡关系，但是更多的是文学上的惺惺相惜和长期交往所建立的深厚感情。苏轼对于范温十分关怀，其父范祖禹在绍圣五年(1098年)去世后，苏轼曾给范冲写了十三封书信。这些信并非写给范冲个人，而是写给兄弟二人。信中频繁出现“昆仲”以及“呈此纸令弟”，劝兄弟二人节哀，自立自爱，“惟昆仲金石乃心，困而不折，庶几先公之风没而不亡也” [ 2 ]。元符三年(1100年)，苏轼在信中苏轼哀叹“哀哉少游，痛哉少游，遂丧此杰耶!”范温的岳父秦观也去世了，两年之内父亲与岳父相继去世，范温受到的打击是巨大的。尽管苏轼也身处恶劣的环境，依然不忘写信劝慰范温。苏范两家之间经过长期的交游已经建立了深厚的感情，在遇到挫折时这种感情就显得更加难能可贵。 对于范温的文学修养提升最大的是黄庭坚，黄庭坚与眉山苏氏的联系密切，范温自然也会有很多接触到苏氏家族的机会。范温师承黄庭坚，跟随黄庭坚学诗，风格与黄庭坚几乎一致。吕本中曾谈到：“表叔范元实既从山谷学诗，要字字有来处” [ 22 ]。在婚姻方面，范温是秦观的女婿，并以此为荣，其学术也受到秦观的影响。黄庭坚和秦观都是苏门中的佼佼者，在文学上成就卓著，与张耒、晁补之合称为“苏门四学士”。在与苏轼门人的接触中，范温受到了苏氏学术的深刻影响。 除了受到苏门名士的文学指导，范温与同辈的苏学传人也不乏交流，苏轼的幼子苏过、秦观之子秦湛都与范温交往甚密。范温与秦湛交好，黄庭坚曾以诗“秦范波澜阔，笑陆海潘江” [ 14 ] 说明了两人关系密切。这句诗化用了南朝钟嵘《诗品》中的“陆才如海，潘才如江”，说明范温与秦湛常在一起探讨文学。两人之间的交游也是促成秦范两家联姻的一个重要原因，王文诰曾说“梦得即范祖禹，小范乃其次子，元长之弟也。少游至是始道其婚姻之故，其到雷不久审矣。其子湛，极有文名，自是秦范相合” [ 23 ]。两家实现婚姻上的结合之后，范温同秦湛接触更多，秦观去世后，范温与秦湛护送其丧北归，“秦少游之子湛，自古藤护丧北归，其婿范温候于零陵，同至长沙，适与山谷相遇” [ 24 ]。在护送途中，除了在失去亲人上的互相安慰，必定也有文学上的切磋。 范温在文学上能够有所成就，很可能是由于少年时期李廌为其奠定的文学基础，以及处在苏轼、黄庭坚、秦观等人所组成的以擅长文学著称的学术交游圈，当然，其自身的文学才能也不容忽视。
system: [SYSTEM]

user: 4. 余论 范冲与范温学术倾向的选择除了受到家学与交游圈的影响，与二人的性格也有一定的关系。编修史书需要史学家有严肃认真的态度，能够据实而书。范冲能够“专史笔”，主要是他“讲明是正”的功劳 [ 10 ]。他在重修《神宗实录》的过程中，并没有因为偏私旧党，就将原来的版本全部销毁，而是保留了墨本和朱本的旧文，可以看出范冲严谨的性格。胡昭曦先生认为李焘在修撰《续资治通鉴长编》时能够见到多次修撰的《神宗实录》，范冲存史功不可没 [ 25 ]。 范冲求实谨慎的治史态度在其编类司马光《记闻》时更为清晰。考虑到司马光“平生记录文字甚多”，战乱之中，“所存无几”。但是司马光颇具盛名，其“书已散落于世”，上至公卿，下至平民，皆有议论。在编类的过程中，范冲保留《记闻》的原貌，“据所录，疑者传疑，可正者正之；阙者从阙，可补者补之；事虽叠书，而文有不同者，两存之” [ 5 ]。为后世研究提供了更多的材料，是严谨的治史态度。 范冲曾写下“七不言，七不可”的座右戒以自警：“一、不言朝廷利害，边报差除……七、见人富贵，不可叹羡诋毁。凡此数事，有犯之者，足以见用意之不肖，于存心修身大有所害，因书以自警” [ 26 ]。范冲写下此座右戒，对自己在朝廷为官、日常生活、行为习惯、与人相处等方面严格要求，修身养性。座右戒体现了范冲的谨重，在言行举止上一丝不苟，与史学的严谨要求相契合。他“名德老成，极天下之选”，在给皇帝讲论《左氏春秋》时，“敷衍经旨，因以规讽”，得到皇帝赞赏“德行文学，为时正人” [ 4 ]。正直敢言是作为一名史学工作者必不可少的条件，能够详细记录，还原史实，并且给出客观的评价。史学也需要史官纵贯古今，范冲是一个博学之人。吕颐浩曾称赞范冲“世德硕学，徊翔列卿” [ 10 ]，钦佩之意溢于言表，高宗也曾提到范冲是“敦厚博硕之士”，对这个“名臣之后”评价颇高 [ 10 ]。 与史学严谨的性格不同，文学更为生动，喜欢文学之人大多活泼幽默、爱憎分明。苏轼就是一个典型的例子，苏轼为人放荡不羁，严守礼法的范祖禹就不能接受，“东坡好戏谑，言语或稍过，纯夫必戒之。东坡每与人戏，必祝曰：勿令范十三知” [ 3 ]。范温虽不似苏轼一般放荡不羁，但也具有文人幽默的个性： 温尝预贵人家会，贵人有侍儿，善歌秦少游长短句，坐闲略不顾，温亦谨，不敢吐一语。及酒酣欢洽，侍儿者始问：“此郎何人耶？”温逮起，叉手而对曰：“某乃‘山抹微云’女婿也。”闻者多绝倒 [ 27 ]。 可见范温谦恭儒雅，在侍儿歌唱秦少游诗词时，认真聆听，不发一言，在别人问及时，才表明自己的身份。但在介绍自己时，也不失幽默风趣，考虑到在座诸人对秦观的诗词更为熟悉，不直接点明身份，而是叉手而立，骄傲地用“山抹微云”来形容岳父秦观，达到了众人皆倾倒的目的。 范温也是一个不拘小节，敢说敢做之人。他对于时人的不当言行敢于提出批评：“我不解今时士大夫，不使人明目张胆直道而行，率要作匿情诡行，似王莽日事沽吊。是谁倡此！岂世美事耶？” [ 27 ] 可以看出范温是一个坦诚、爱憎分明之人。范温敢说敢做的性格对于《潜溪诗眼》的写作也具有重要的作用，敢于批评，善于评价。 虽然范冲与范温同为范祖禹后代，但是二人在学术上还是存在着较大的差异。范冲在史学上颇有建树，范温在文学上成就斐然。本文对二人文史异辙的情况及其原因的探析，还比较粗浅，比较单薄，希望今后有机会进一步发掘资料，做更深入的讨论。
system: [SYSTEM]

user: 摘要 极限思想是许多科学领域的重要思想之一，在函数求极限中大家经常会运用洛必达法则或者泰勒定理，为了解决某些数列极限问题，本文介绍了求数列极限的一种方法——Stolz定理，并对Stolz定理进行一些讨论。
system: [SYSTEM]

user: 摘要 采用大涡模拟(LES)的方法求解不可压缩流体连续方程和N-S方程，对格栅紊流场的均匀性，以及IEC阵风与紊流共同作用下的流场特性进行了数值模拟分析。通过分析格栅紊流流场特性参数，得出本文格栅流场在格栅下游4 m处的均匀性达到要求。IEC阵风与格栅紊流场结合较好，可用于后续桥梁断面特性的研究。
system: [SYSTEM]

user: 1. 引言 风工程研究中，常采用主动模拟和被动模拟两类方法来模拟自然环境中风的紊流特性。主动模拟分为振动格栅 [ 1 ]、振动尖塔 [ 2 ] 和多风扇风洞 [ 3 ] 等方法，被动模拟分为被动格栅 [ 4 ]、尖劈 [ 5 ] 和粗糙元 [ 6 ] 等方法。在风洞实验中，一般采用在风洞入口处布置被动格栅的方式来模拟均匀紊流场。张明月 [ 7 ] 利用风洞试验的手段研究了测点与格栅间距、格栅板宽度以及单元格栅边长等对风参数的影响规律。文水兵 [ 8 ] 通过改变格栅板宽度以及间距的方式，调试出了A、B、C三类地貌的紊流强度。 随着计算机技术的日渐成熟，计算流体动力学(Computational Fluent Dynamics)技术在风工程中的运用越来越多。且较风洞试验而言，数值模拟具有较强的灵活性、流场可观测和成本低等特点，成为了一种验证风洞试验的可靠手段。鲍泽辰等 [ 9 ] 将粗糙元、扰流杆以及格栅三者结合，采用LES的方法，模拟了C类地貌的风场特性。杨荣菲等 [ 10 ] 通过LES的方法，探讨了格栅稠度、来流风速以及格栅表面粗糙度，对格栅紊流场特性参数的影响。Torrano等 [ 11 ] 使用LES方法对被动格栅紊流场进行了数值模拟，并得出结论，y + 在30~300之间也可以得到较好的流场特性参数。验证了运用LES紊流模型模拟格栅紊流场的可行性。 理想化的阵风模型可分为TSI阵风模型和IEC阵风模型 [ 12 ] 等，但有研究表明将TSI阵风模型运用于桥梁抗风研究存在不足 [ 13 ]。IEC阵风模型是国际电工委员会制定的IEC61400-1标准中提出的，其主要特点表现为风速在短时间内减小到最小值，接着又急剧增大达到最大风速值，随后又在极短时间内迅速下降再次变为最小值，最后恢复为阵风发生前的速度。赵元元 [ 14 ] 利用流体力学前处理软件Gambit建立起1.0 MW风机的叶片的三维模型，并通过UDF用户自定义函数接口将极端运行阵风函数程序导入Fluent中作为流场入口条件，验证了CFD方法模拟极端运行阵风的正确性。Storey等 [ 15 ] 基于国际风轮机设计标准IEC61400-1建立了一种极端相关阵风模型，使用数值模拟中的LES方法模拟了一组风力机在极端瞬态风条件下的运行情况。基于此，本文将运用CFD (Computational Fluid Dynamic)中LES的方法，通过Fluent的UDF接口将IEC阵风模型作为流场入口风速条件，将格栅紊流场与IEC阵风模型结合，研究其可行性以及流场特性。
system: [SYSTEM]

user: 2. 紊流特性参数 紊流场的特性参数主要包括紊流强度、紊流积分尺度以及脉动风速功率谱等。 2.1. 紊流强度 紊流强度是紊流场中最简单但使用最多的参数，紊流强度是指风速根据空间和时间而变化的程度，通常定义为脉动风速的标准差和脉动风速的平均值之比。 I u = U r m s U ¯ (1) 式中： U r m s 为脉动风速的标准差， U ¯ 为脉动风速的平均值。 2.2. 紊流积分尺度 紊流积分尺度是度量紊流场中某一方向上漩涡尺寸平均值的一个指标。格栅紊流场中有包括顺流向、横向以及垂直三个方向的漩涡，每一个漩涡又有x、y、z三个方向的尺度，因此紊流积分尺度一共有九个。紊流积分尺度可由顺流向相邻两点的互相关函数定义，以顺流向的积分尺度为例，其表达式可写为： L u z = 1 U r m s 2 ∫ 0 ∞ R u 1 u 2 ( z ) d z (2) 式中： U r m s 2 是脉动风速 u 1 和 u 2 的方差； R u 1 u 2 ( z ) 是顺流向相邻两点 u 1 = ( x 1 ，   y 1 ，   z 1 ，   t ) 和 u 2 = ( x 1 ，   y 1 ，   z 1 + z ，   t ) 脉动风速的互相关函数，其中t是时间。但在处理数据过程中，常常引用泰勒尔紊流冻结假说 [ 16 ] 假设，即假设漩涡在紊流场中是以平均速度向下游迁移。此时，脉动风速 u 1 = ( z 1 + z ，   t ) 则可以改写为 u 1 = ( z 1 ，   t + τ ) ， z = U ¯ × τ 。公式也可改写为： L u z = 1 U r m s 2 ∫ 0 ∞ R u ( t ) d τ (3) 式中： R u ( t ) 是脉动风速 u 1 = ( z 1 ，   t + τ ) 的自相关函数。在数据处理过程中，由于自相关系数比较小，会导致泰勒尔假设会有比较大的误差。因此，通常取(3)式的积分上限为 R u ( t ) = 0.05 U r m s 2 。 2.3. 脉动风功率谱 脉动风功率谱表述的是，紊流的运动能量在不同频率上的分布情况。经过多年来国内外学者对理论与试验的研究，目前有多种类型的脉动风速谱，对于顺流向的脉动风功率谱，常用Von Karman谱 [ 17 ] 来描述脉动风速的特性。 n S u ( z ，   n ) u * 2 = 4 β f ( 1 + 70.8 f 5 ) 5 / 6 (4) f = n L u x U (5) 式中： S u ( z ，   n ) 是顺流向脉动风速功率谱密度函数；z是地面或水面以上的高度，单位(m)；n是顺流向风的脉动频率； u * 是摩阻速度（剪切速度），单位(m/s)； β 为摩阻速度修正系数； L u x 为紊流积分尺度，单位(m)；U为平均风速，单位(m/s)。
system: [SYSTEM]

user: 3. 计算模型 3.1. 格栅几何尺寸 便于将数值模拟结果与TJ-2风洞试验数据做对比，本文选用的格栅尺寸与风洞试验中的格栅尺寸相同，水平和竖向单元格栅板条宽度a均为15 cm，格栅板条围成正方形的中心孔，水平和竖向板条之间的净距b均为45 cm。格栅局部示意图如图1所示，水平和竖向均有四个格栅孔，格栅整体示意图如图2所示。 图1. 格栅局部示意图 图2. 格栅整体示意图 3.2. 模型建立及网格划分 数值模拟格栅紊流场计算域采用三维长方体计算域，宽度为2.4 m，高度为2.4 m，长度为10.5 m，格栅放置于距流场入口4.5 m处。坐标原点位于流场入口的左下方，顺流方向为z正方向，竖直方向为y轴，向上为正，水平方向为x轴，由笛卡尔右手法则确定。计算域示意图如图3所示。 图3. 计算域示意图 本文利用Fluent前处理软件Gambit对整个流场域进行网格划分。在格栅附近由于流动比较复杂，为保证结果的精确性，因此需要划分较精细的网格。贴体网格首层高度为0.75 cm，以1.02的增率沿z轴增加。在格栅附近采用非结构化网格对流场域进行划分，其余流场域采用结构化网格划分。格栅附近网格划分示意图如图4所示，整个计算模型网格数量约为330万。 图4. 格栅附近网格划分 本文采用大涡模拟(Large Eddy Simulation， LES)湍流模型模拟格栅紊流场。LES湍流模型是通过滤波的方式，把格栅流场中的旋涡分解为大尺度和小尺度旋涡。对于大尺度的旋涡，采用直接计算的方法。而对于小尺度的旋涡，采用压格子尺度(Sub-Grid Scale， SGS)模型进行模拟。在FLUENT求解器中，基于有限体积法(Finite-Volume Method， FVM)求解经过滤波处理后的不可压缩流体连续方程和Navier-Stokes方程。亚格子尺度模型采用Smagorinsky-Lilly亚格子尺度模型，Smagorinsky常数C s 取0.1 [ 11 ]。运用SIMPLE算法解决压力与速度耦合问题，压强方程采用二阶迎风差分格式(Second Order Upwind)，动量方程采用有界中心差分格式(bounded central differencing)。迭代时间步长取为0.001 s，共计算20，000步，即计算总时长为20 s。计算域的详细边界条件如表1所示。 Table 1 Boundary condition 位置 边界条件 入口 Velocity-inlet， U = 10 m/s 出口 Pressure-outlet 格栅 No slip wall 前后面 Symmetry 上下面 Symmetry 表1. 边界条件 监测面位于流场下游，距格栅4 m处。监测点间距为0.4 m，均匀布置于监测面内，监测点位置如图5所示。每个监测点记录z方向的风速，采样频率为1000 Hz，采样点数为20，000。
system: [SYSTEM]

user: 4. 紊流结果分析 本文从每个监测点获得的风速时程曲线进行分析，处理得到格栅紊流场的特性参数。由于流场在刚开始计算时不稳定，导致风速时程曲线的数据没有太大使用价值，因此本文在处理数据时，省去了前两秒的风速时程数据。典型风速时程曲线如图6所示。 图5. 监测点布置图(单位：cm) 图6. 风速时程曲线 图7为格栅紊流场的流线图，可以看出在格栅前为均匀流场。当经过格栅时，由于格栅条是离散布置的，因此，均匀流场被格栅打乱，在格栅后方形成了紊乱的流场。但经过一段距离的发展和扩散之后，流场又恢复了均匀性。 图7. 流场流线图 4.1. 紊流强度 格栅紊流场是由均匀流通过格栅装置后形成的，由于格栅条是离散布置的，因此均匀流刚通过格栅时，紊流特性是及其不均匀的。但旋涡经过一段距离的扩散，紊流特性的均匀性逐渐改善，形成各项同性的均匀紊流。本文通过对监测点的风速时程曲线分析得到的紊流强度，来验证格栅紊流场的均匀性。 表2给出了监测面上个点的坐标及紊流强度值，紊流强度最小值为13.4%，最大值为14.9%，平均值为14.3%，标准差为0.004。由此可知，该格栅流场在距格栅4 m的下游处，流场已经达到了较好的均匀性。本文数值模拟的紊流强度，与风洞试验 [ 8 ] 中的紊流强度I u = 15%相比，相对误差为4.67%。与试验数据较吻合，验证本文采用的数值模拟手段是可行的，该湍流场断面可用于开展节段模型的研究。 Table 2 Turbulence intensitie 坐标(m) I u (%) 坐标(m) I u (%) 坐标(m) I u (%) 坐标(m) I u (%) 坐标(m) I u (%) (0.4， 0.4) 14.7 (0.4， 0.8) 14.7 (0.4， 1.2) 14.9 (0.4， 1.6) 14.7 (0.4， 2.0) 14.9 (0.8， 0.4) 13.6 (0.8， 0.8) 13.8 (0.8， 1.2) 13.8 (0.8， 1.6) 14.6 (0.8， 2.0) 14.6 (1.2， 0.4) 14.4 (1.2， 0.8) 13.7 (1.2， 1.2) 13.4 (1.2， 1.6) 14.0 (1.2， 2.0) 14.6 (1.6， 0.4) 14.2 (1.6， 0.8) 14.3 (1.6， 1.2) 13.6 (1.6， 1.6) 14.2 (1.6， 2.0) 14.5 (2.0， 0.4) 14.4 (2.0， 0.8) 14.9 (2.0， 1.2) 14.0 (2.0， 1.6) 13.8 (2.0， 2.0) 14.2 表2. 紊流强度 4.2. 紊流积分尺度 本文同样对监测面上个监测点的风速时程信号进行分析，得到了各监测点的积分长度尺度，详细数据如表3所示。其中最大的紊流积分尺度为0.371 m，最小为0.231 m，平均值为0.279 m，标准差为0.0345，说明该格栅紊流场在此断面处已有足够的均匀性。与文献 [ 4 ] 中的结果相比较，相对误差约为16.7%，本文数值模拟结果偏大。其主要原因在于，本文数值模拟的计算条件较差，无法使用庞大数量的网格进行计算。因此导致在数值模拟时，不能捕捉到小尺度的旋涡。 Table 3 Turbulence integral scal 坐标(m) L x (m) 坐标(m) L x (m) 坐标(m) L x (m) 坐标(m) L x (m) 坐标(m) L x (m) (0.4， 0.4) 0.259 (0.4， 0.8) 0.371 (0.4， 1.2) 0.274 (0.4， 1.6) 0.351 (0.4， 2.0) 0.293 (0.8， 0.4) 0.236 (0.8， 0.8) 0.276 (0.8， 1.2) 0.277 (0.8， 1.6) 0.246 (0.8， 2.0) 0.252 (1.2， 0.4) 0.248 (1.2， 0.8) 0.258 (1.2， 1.2) 0.231 (1.2， 1.6) 0.261 (1.2， 2.0) 0.312 (1.6， 0.4) 0.299 (1.6， 0.8) 0.308 (1.6， 1.2) 0.237 (1.6， 1.6) 0.293 (1.6， 2.0) 0.305 (2.0， 0.4) 0.324 (2.0， 0.8) 0.258 (2.0， 1.2) 0.287 (2.0， 1.6) 0.262 (2.0， 2.0) 0.270 表3. 紊流积分尺度 4.3. 脉动风速谱 利用监测点测得的风速时程信号，对其进行傅里叶变换，编写将风速时程信号转化为脉动风速谱的程序，利用MATLAB软件进行处理，即可得到脉动风速谱。本文以顺风向，高度为0.4 m的监测点为例，得到脉动风速谱如图8所示。并Von Karman谱作比较，可以看出数值模拟得到的脉动风速谱，在低频段对比Von Karman谱偏小，但数值模拟的风场能量结构与Von Karman谱基本能够吻合。 图8. 脉动风速谱
system: [SYSTEM]

user: 5. IEC阵风与紊流场结合 本文采用Länger-Möller [ 12 ] 根据国际通用标准IEC61400-1建立的极端运行阵风(EOG)模型，用作模拟流场入口处的风速： U ( t ) = { U − 0.37 U g sin ( 3 π ( t − 2 ) T g ) ( 1 − cos ( 2π ( t − 2 ) T g ) )                         2 ≤ t ≤ 12.5 U                                                                                                                                                 t 取 其 他 值 (6) 式中：U g 为与阵风风速相关的参数；T g 为特征周期时间取10.5 s；U为均匀来流速度。 初始风速为10 m/s，持续时间为15 s时，IEC阵风模型风速时程曲线如图9所示。 图9. IEC阵风模型风速时程曲线 图10为阵风与格栅共同作用下，距格栅4 m的下游流场中心点处，顺流向的风速时程曲线。由于在t = 12.5 s后，IEC阵风风速持续为10 m/s，所以本文仅计算15 s的时长，迭代时间步长为0.001 s。由于刚开始计算时，流场尚未充分发展，因此0~1 s内的风速时程曲线振幅很大。在最大阵风风速时刻附近，紊流场内风速时程曲线振幅也较其他时刻偏大，但整体数值模拟的风速曲线与理论值吻合较好。 图11展示了格栅流场下游4 m处监测面上不同时刻顺流向的速度云图。t = 2.0 s和t = 12.5 s时速度横截面上平均速大致为10 m/s，t = 4.5 s和t = 10.0 s时速度横截面上平均速大致为7 m/s，t = 7.5 s时速度横截面上平均速大致为17 m/s。各时刻风速与IEC理论风速相对应，且监测面上的风速分布均匀，可说明IEC阵风与格栅紊流场结合较好，可在该断面上开展桥梁段面特性的研究。 图10. 阵风与紊流共同作用下的风速曲线 图11. 格栅下游4 m处不同时刻顺流向速度云图 图12是通过CFD后处理软件CFD-POST中的Q-criterion获得的涡量图。当空气经过格栅后，在格栅附近形成较小的旋涡结构，通过一段距离的发展和扩散，漩涡结构逐渐变大。可得出，积分长度尺度是沿顺流向逐渐增大。不同时刻，整个流场区域的漩涡结构分布一致，即在格栅附近为小旋涡结构，而后逐渐变大。IEC阵风的风速变化并不会破坏紊流场的漩涡分布，紊流场与IEC阵风能够很好地结合起来，可用于后续桥梁断面特性的研究。 图12. 不同时刻的涡量图 图13为流场中心线上不同时刻的静压图，任意时刻流场中静压仅在格栅附近有较大波动。在格栅前流场中心线上的静压为正，在靠近格栅处急剧上升。在格栅后流场中心线上的静压为负，并快速上升至零。流场中静压与风速密切相关，从IEC风速时程曲线上可以看出，t = 2.0 s和t = 12.5 s、t = 4.5 s和t = 10.0 s时的风速几乎是相同的。因此，图13(a)中t = 2.0 s和t = 12.5 s、t = 4.5 s和t = 10.0 s静压曲线几乎是重合的。图13(b)展示的是t = 2.0 s、t = 12.5 s和t = 7.5 s时流场中心线上的静压曲线，t = 7.5 s时对应的风速最大。因此，格栅前方的正压和格栅后方的负压，分别达到了整个模拟过程的最大值。相反，t = 4.5 s时流场中风速最小，流场中心线上的静压也最小。 图13. 流场中心线静压曲线图
system: [SYSTEM]

user: 6. 结论 本文采用LES湍流模型对格栅紊流场和IEC阵风与紊流共同作用下的流场进行了数值模拟，分析流场得到了以下结论： 1) 通过与风洞试验得到的结果对比，发现数值模拟得到的格栅紊流场特性满足要求。数值模拟的方法能够任意布置监测点，并且能直接观察到流场中不同时刻旋涡的分布情况。相对于风洞试验拥有更便于观测、更省时等优点。 2) 数值模拟得到的紊流场紊流强度与风洞试验吻合较好，紊流积分尺度稍微偏大，可能是计算条件不足、网格不够精细导致的，紊流脉动风速谱与Von Karman谱也基本能够吻合。可以利用数值模拟的方法，开展格栅紊流场特性的研究。 3) IEC阵风与紊流共同作用下的脉动风速时程曲线与IEC阵风理论风速曲线吻合较好，不同时刻流场中漩涡结构分布规律一致，不同时刻流场中心线上静压分布规律相同，说明IEC阵风与格栅紊流场结合较好。可用于后续桥梁断面特性的研究。
system: [SYSTEM]

user: 摘要 三轴立式加工中心是机械加工的重要装备，为提高其精度，本文利用齐次坐标变换建立了三轴立式加工中心的几何误差数学模型，开发了基于i5OS数控系统的误差补偿软件，为了简化误差检测时间和效率，选取对综合误差有主要贡献的直线度和定位误差进行补偿实验，结果表明，基于i5OS数控系统的综合误差实时性反馈补偿使数控机床的加工精度提高了38%以上，验证了模型的准确性和补偿软件的有效性。
system: [SYSTEM]

user: 1. 引言 精度是机床的重要指标，为了提高精度，“误差补偿”是目前常用的一种有效手段，所谓“误差补偿”是事先分析数控机床误差产生的机理，用数学的手段建立其误差数学模型，以此来预测机床工作时各空间位置的误差，最后利用软件处理的方式实现补偿，达到提高机床精度的目的 [ 1 ]。 数控机床综合误差数学模型是将数控机床空间误差用数学手段进行表达，从而能够预测其在数控机床整个空间位置上的误差。目前，主要手段是基于多体系统理论建立数控机床空间误差的数学模型。多体系统理论方法最大的优势在于其全面性以及通用性，将数控机床运动过程中的耦合因素进行综合考虑，适用于各种运动形式的数控机床 [ 2 ] [ 3 ] [ 4 ] [ 5 ] [ 6 ]。空间误差数学模型表达式的变量就是数控机床工作时的坐标值 [ 7 ]。 本文首先将利用齐次坐标法对三轴立式加工中心进行误差建模方法展开研究，采用激光涉仪实现误差的检测，最后基于国产数控系统i5OS开发出误差补偿系统并对模型进行有效性验证。
system: [SYSTEM]

user: 2. 数控机床几何误差建模 一般情况下，机床每一个移动轴均存在3项位移误差以及3项角度误差 [ 3 ]，并且三个移动轴两两之间还存在1项垂直度误差。三轴立式加工中心一共存在21项几何误差元素 [ 8 ]。本文采用沈阳机床i5 M1.4立式加工中心为研究对象，结构模型如图1所示，其运动形式为YXTZ型，即Y轴以及Z轴与床身固连，X轴与Y轴相连，工件台于X轴相连，而主轴以及刀具与Z轴相连。 按照多体系统理论描述为：① 刀具链分支：机床床身F→Z轴→主轴S→刀具T；② 工件链分支：机床床身F→Y轴→X轴→工件W。 以床身坐标系为参考坐标系，其刀具分支的运动链为：床身F→Z轴→主轴S→刀具T，则理想情况下，从床身坐标系到刀具坐标系齐次坐标变换关系如下： T F T = ∏ i = Z T T i − 1 i = T F Z T Z S T S T (1) 同理，工件分支的运动链为：机床床身M→Y轴→X轴→工件W，则理想情况下，从床身坐标系到工件坐标系齐次坐标变换关系如下： T F W = ∏ i = Y W T i − 1 i = T F Y T Y X T X W (2) 设刀尖点 t 0 在刀具坐标系中的位置坐标为 ( t x ， t y ， t z ) ，则理想情况下，刀尖点在工件坐标系中的位置坐标点 t w 为： [ t w T 1 ] = T W T [ t 0 T 1 ] = [ T F W ] − 1 T F T [ t 0 T 1 ] (3) 则进一步可以得到： [ t w T 1 ] = T X W − 1 T Y X − 1 T F Y − 1 T F Z T Z S T S T [ t 0 T 1 ] (4) 同理，在实际情况下，只需要在理想情况下齐次坐标变换过程中将各相邻体间误差矩阵引入表达式，即可得到实际情况下刀尖点在工件坐标系中的坐标点 t ′ w ，为： [ t ′ w 1 ] = T W T Δ T W T [ t 0 T 1 ] = ( ∏ i = Y W T i − 1 i Δ T i − 1 i ) − 1 ( ∏ i = Z T T i − 1 i Δ T i − 1 i ) [ t 0 T 1 ] (5) 因此， ε n = ( ε x ， ε y ， ε z ) 机床空间综合误差可以表达为： [ ε n T 0 ] = ( T i = Y T T i − 1 i ) − 1 ( T i = Z S T i − 1 i ) [ t 0 T 1 ] − ( T i = Y T T i − 1 i Δ T i − 1 i ) − 1 ( T i = Z S T i − 1 i Δ T i − 1 i ) [ t 0 T 1 ] (6) 由于单项误差元素都是微小量，忽略其高阶项，得到空间综合误差在X轴，Y轴，Z轴方向上的分量 [ 3 ]： { ε x = − δ x x − δ y x + δ z x − y ( θ x z + θ y z − ε x y ) − z ( θ x y + θ y y + ε x z ) ε y = − δ x y − δ y y + δ z y + x θ x z + z ( θ x x + θ y x − ε y z ) ε z = − δ x z − δ y z + δ z z − x θ x y + y ( θ x x + θ y x ) (7) 由上述数控机床的综合数学模型建模过程可知，刀具切削点的空间定位精度同时受到21项误差元素的影响。在机床误差的实际补偿中，如果对这21项误差元素都进行测量与补偿，往往会导致投入的时间与经济成本比较大，补偿效率不高。实际上，不同类型的机床由于其结构类型、运动行程大小、以及运动轴定位方式的不同，所面临的关键误差是不同的。对于立式加工中心而言，由于其运动轴行程短，角偏误差的影响往往比较小，可以忽略不计，定位误差以及直线度误差对误差有主要贡献。 图1. 沈阳i5 M1.4立式加工中心。(a) 机床结构示意图；(b) 机床误差检测实验现场
system: [SYSTEM]

user: 3. 数控机床误差测量 本文实验采用的误差测量仪器为Renishaw XL 80型激光干涉仪，对沈阳i5 M1.4立式加工中心X轴、Y轴进行误差测量。每个移动轴需要对线性定位误差以及直线度误差进行测量，需要用线性定位误差测量镜组以及直线度误差测量镜组。 实验中测量长度应为各移动轴的最长有效行程。但是受到被测数控机床结构的约束，各移动轴的测量行程分别为：X轴：0~500 mm、Y轴：−300~0 mm、Z轴：−300~0 mm。各移动轴从起始测点开始等间距选取51个测量点，且沿正、反向测量两次，再取其平均值。对于测量误差数据利用一次线性插值法进行单项误差建模。实验测量结果以及建模结果如图2所示。 图2. 实验测量结果以及建模结果。(a) X轴直线度和定位误差；(b) Y轴直线度和定位误差
system: [SYSTEM]

user: 4. 基于i5OS数控系统的机床综合误差实时补偿 沈阳机床(上海)智能系统研发设计有限公司i5OS数控系统通过Finch高速通信组件将机床CNC参数信息提供给上层应用软件访问，在上层软件中通过Finch高速通信组件访问CNC各轴位置反馈信息，再此时各轴位置信息写入误差模型中，处理计算得到此时误差补偿值，继而将其写入Finch高速通信组件，最后写入CNC运动控制模块，完成机床误差的补偿。Finch是一个主要用于i5OS本地以及局域网访问的高速通信组件，它位于应用软件中间层，主要依赖于动态链接库libfinch.so，具有很高的开放性。该组件不限制读取的参数数据，但是会限制参数的写入，根据CNC的运行状态和授权权限来决定能够写入的参数(图3)。 图3. Finch高速通信组件通信示意图 沈机(上海)智能系统研发设计有限公司i5OS数控系统通过Finch高速通信组件将机床CNC参数信息提供给上层应用软件访问，因此，在上层软件中通过Finch高速通信组件访问CNC各轴位置反馈信息，再此时各轴位置信息写入误差模型中，处理计算得到此时误差补偿值，继而将其写入Finch高速通信组件，最后写入CNC运动控制模块，完成机床误差的补偿。补偿原理如图4所示，这种补偿方法是一种基于数控系统的动态补偿，将数控机床空间误差与数控系统结合起来，通过外部补偿软件修正数控系统中的各轴位置反馈，实现对数控机床的空间误差补偿。 图4. 基于i5OS数控系统的位置反馈修正补偿方法
system: [SYSTEM]

user: 5. 补偿实验验证 为验证上述误差补偿方法的有效性，对沈阳i5 M1.4立式加工中心X轴以及Y轴联动时圆弧插补运动进行了一组补偿实验，实验仪器选择Renishaw QC20_W型球杆仪。 实验时球杆仪球心坐标为(170， −150， −292.578)，与主轴相连，测量半径为100 mm，X轴，Y轴进给速度设为300 mm/min，先逆时针方向再顺时针方向旋转360度，完成正向与反向圆弧插补运动的误差检测。 按照JIS B 6190-4:2008版本的标准分析补偿前后球杆仪所测的误差数据，两次测量所得的运动误差轨迹及其圆度偏差分析分别如图5与图6所示。补偿前，XY平面内的圆弧轨迹最大圆度误差为17.9 μm，X轴、Y轴的垂直度误差为0.0934 μm/mm；补偿后XY平面内圆弧轨迹的最大圆度误差为10.3 μm，X，Y轴的垂直度误差为0.0576 μm/mm。最大圆度误差补偿率为42.45%，垂直度误差补偿率为38.33%。 图5. 补偿前圆弧轨迹精度 图6. 补偿后圆弧轨迹精度 表1显示了沈阳i5 M1.4立式加工中心在XY平面内圆弧插补运动在补偿前后的整体效果。 Table 1 Overall effect of circular interpolation before and after motion compensatio 误差项 补偿前 补偿后 提高百分比 最大圆度误差 17.9 μm 10.3 μm 42.45% 垂直度误差 0.0934 μm/mm 0.0576 μm/mm 38.33% 表1. 圆弧插补运动补偿前后的整体效果
system: [SYSTEM]

user: 6. 结论 通过齐次坐标转换的数学方法，以沈阳机床立式加工中心为研究对象，建立了21项误差组成的机床误差综合模型，利用机床厂商提供的接口文件，开发出基于i5OS数控系统的误差综合实时补偿系统，为了简化误差检测过程、节省时间和提高效率，选取对综合误差有主要贡献的直线度和定位误差进行补偿实验，结果表明，基于i5OS数控系统的综合误差实时性反馈补偿使数控机床的加工精度提高了38%以上，验证了模型的准确性和补偿软件的有效性。
system: [SYSTEM]

user: 摘要 目的：探讨早产儿脑损伤发生的相关危险因素，为降低早产儿脑损伤的发生率及改善远期预后提供临床参考依据。方法：选取2019年8月~2020年8月在我院新生儿重症监护室治疗的55例脑损伤早产儿为观察组；另外选取同时期住院的70例非脑损伤的早产儿为对照组。比较两组母亲因素及早产儿情况，并采用多因素Logistic回归分析早产儿脑损伤发生的危险因素。结果：脑损伤组与无脑损伤组早产儿组孕母因素相比，除产前应用激素和母孕期感染外，分娩方式、羊水异常、胎膜早破、孕母贫血、胎盘异常、脐带异常、妊娠期糖尿病、妊娠期高血压等方面差异均无统计学意义(P < 0.05)。在早产儿组基本情况上，除生后低血糖、胎儿宫内窘迫、呼吸暂停外，新生儿脓毒症、新生儿贫血、动脉导管未闭、机械通气、生后1分钟Apgar评分、新生儿呼吸窘迫综合征和高碳酸血症差别均有统计学意义(P < 0.05)。对早产儿脑损伤的多因素Logistic回归分析，显示新生儿呼吸窘迫综合征、机械通气、母孕期感染、新生儿脓毒症和生后1分钟Apgar评分为脑损伤的独立危险因素，差别均有统计学意义(P < 0.05)。结论：早产儿脑损伤是多种因素相互作用的复杂结果，发生的相关危险因素包括孕母因素及新生儿因素。
system: [SYSTEM]

user: 1. 引言 近年来，随着新生儿重症监护室和的产科的救治水平不断提高，早产儿的存活率显著提高，同时相应并发症影响远期预后的情况也随之出现，早产儿脑损伤所引起的认知障碍、语言障碍、行为心理问题及脑瘫等严重后遗症越发受到重视。胎龄越小，脑损伤的发生率越高，程度越重。早产儿脑损伤主要包括脑室周围白质软化(periventricular leucomalacia， PVL)和脑室内出血(intraventricular haemorrhage， IVH)。本研究主要探讨引起该病的高危因素，以期有效地预防和治疗早产儿脑损伤，减轻社会及家庭的负担。
system: [SYSTEM]

user: 2. 资料与方法 2.1. 一般资料 选取2019年8月~2020年8月于我院NICU治疗的早产儿脑损伤患儿55例为观察组。胎龄小于32周，符合中华医学会儿科学分会新生儿学组制定的相关诊断标准(缺血性脑损伤局灶性缺血性脑损伤急性期MRI表现为半卵圆中心、侧脑室旁簇状或线状的短T1短T2信号，DWI示受累脑组织明显高信号；弥散性缺血性损伤仅DWI表现为侧脑室旁白质大片状高信号改变，神经元轴突损伤表现为大脑皮质、丘脑、小脑、基底核、脑干部位DWI高信号；出血急性期MRI表现为T1WI无异常信号或稍低信号，T2WI 呈高信号，亚急性期早期即出血(第4~7天)，T1WI呈高信号、T2WI稍低信号，慢性期(出血14 d以后) T1WI 和T2WI均为高信号)，同时结合颅脑彩超、头颅磁共振检查及临床症状后确诊，并随机抽取同期住院治疗的70例非早产儿脑损伤患儿为对照组。实验组与对照组胎龄分别为(29.32 ± 1.41)周和(30.57 ± 1.23)周，体重分别为(1.75 ± 0.33) kg和(1.81 ± 0.37) kg。排除标准：① 明确诊断存在遗传代谢病；② 先天性心脏病和其他先天性疾病引起的脑损伤；③ 生后24 h内家属放弃治疗。且所有实验均已通过医院伦理审查委员会审查。 2.2. 研究方法 统计两组母亲因素：高龄产妇、产前应用激素、保胎史、妊娠期糖尿病、妊娠期高血压疾病、妊娠期贫血、羊水异常、脐带异常、胎盘异常及胎膜早破；新生儿情况：细菌性脓毒症、贫血、胎儿宫内窘迫、Apgar评分、呼吸暂停、高碳酸血症、动脉导管未闭(PDA)、新生儿呼吸窘迫综合征(RDS)、低血糖、机械通气。 2.3. 统计学方法 采用SPSS 23.0统计软件进行数据处理，单因素分析中计量资料采用t检验，计数资料采用X 2 检验，早产儿脑损伤发生的危险因素采用多因素Logistic回归分析，P < 0.05为差异具有统计学意义。
system: [SYSTEM]

user: 3. 结果 3.1. 两组的孕母因素比较 脑损伤组与无脑损伤组早产儿组孕母因素相比，除产前应用激素和母孕期感染外，分娩方式、羊水异常、胎膜早破、孕母贫血、胎盘异常、脐带异常、妊娠期糖尿病、妊娠期高血压等方面差异均无统计学意义(P < 0.05)。见表1。 3.2. 两组新生儿情况比较 脑损伤组与无脑损伤组早产儿组基本情况相比，除生后低血糖、胎儿宫内窘迫、呼吸暂停外，新生儿脓毒症、新生儿贫血、动脉导管未闭、机械通气、生后1分钟Apgar评分、新生儿呼吸窘迫综合征和高碳酸血症差别均有统计学意义(P < 0.05)。见表2。 3.3. 早产儿脑损伤的多因素Logistic回归分析 对产前应用激素、母孕期感染、新生儿脓毒症、新生儿贫血、动脉导管未闭、机械通气、生后1分钟Apgar评分、新生儿呼吸窘迫综合征和高碳酸血症等高危因素进行多因素Logistic回归分析，显示新生儿呼吸窘迫综合征、机械通气、母孕期感染、新生儿脓毒症和生后1分钟Apgar评分为脑损伤的独立危险因素，差别均有统计学意义(P < 0.05)。见表3。 Table 1 Comparison of maternal factors between two group 组别 例数 产前激素 孕期感染 分娩方式 妊娠期高血压 羊水异常 是 否 是 否 顺产 剖宫 是 否 是 否 观察组 55 18 37 26 29 23 22 15 40 12 43 对照组 70 47 23 14 56 38 32 20 50 19 51 X 2 7.261 8.314 0.497 0.096 2.374 P值 0.025 0.022 0.524 0.894 0.186 组别 例数 胎膜早破 孕母贫血 胎盘异常 妊娠期糖尿病 脐带异常 是 否 是 否 是 否 是 否 是 否 观察组 55 18 37 12 43 11 44 10 45 19 36 对照组 70 23 47 12 58 17 53 15 45 24 46 X 2 0.003 0.264 0.296 0.286 0.003 P值 0.983 0.618 0.596 0.604 0.976 表1. 两组的孕母因素比较 Table 2 Comparison of newborn condition between two group 组别 例数 脓毒症 贫血 PDA 胎儿宫内窘迫 机械通气 是 否 是 否 是 否 是 否 是 否 观察组 55 24 31 36 19 12 43 3 52 27 28 对照组 70 16 54 35 35 6 64 5 65 18 52 X 2 7.156 4.731 9.315 1.349 8.486 P值 0.026 0.038 0.019 0.236 0.021 组别 例数 Apgar评分 RDS 呼吸暂停 高碳酸血症 低血糖 <7 ≥7 是 否 是 否 是 否 是 否 观察组 55 20 35 37 28 16 39 22 33 18 37 对照组 70 16 54 24 46 17 53 18 52 25 45 X 2 4.296 11.694 1.211 4.412 1.243 P值 0.041 0.011 0.281 0.035 0.513 表2. 两组新生儿情况比较 Table 3 Multiple Logistic regression analysis of risk factors of brain injury in preterm infant 危险因素 参考估计值 OR 95%CI P RDS 3.152 9.247 4.219~25.634 P < 0.05 机械通气 2.263 8.741 2.361~19.351 P < 0.05 孕期感染 1.156 4.394 1.263~11.237 P < 0.05 脓毒症 −1.647 0.193 0.059~0.511 P < 0.05 Apgar评分 −1.2634 0.513 0.276~0.821 P < 0.05 表3. 早产儿脑损伤影响因素的多因素Logistic回归分析