风能是清洁的可再生能源，总量丰富，是目前商业开发前景最好的可再生能源之一。海上风场稳定性高、空间广阔，十分适合大规模风电工程建设。随着水深的增加，为节约建设成本，近海风电逐步采用浮式风机替代固定式风机。浮式风机支撑平台大多采用单柱式、张力腿式和半潜式支撑平台基础。

报警管理及Web服务器上运行AAS软件，包括Web服务软件，该服务器通过防火墙与报警管理数据采集Buffer机和公司管理信息网相连，同时实现了安全隔离。该服务器上运行实时报警模块，为中心控制室的客户端提供实时报警及报表查询功能。Web服务器通过防火墙获取报警与事件数据，为公司管理信息网上的客户端提供实时报警及报表查询功能。

细婆说：“香就好。慢慢吃，莫着急呀。”说完，给罗爹爹也端过一碗，又说：“罗爹爹，今天没有去东湖打拳？”

浮式风机系统的动力响应呈现明显的非线性耦合特性；如何分析和预报其复杂动力响应，近年来已经成为业界学者的研究重点之一。Wayman[1]忽略非线性和瞬态效应，在频域下分析张力腿式浮式风机动力特性。Ma[2]等利用FAST软件研究了浮式风机运动性能。Matha[3]考虑非线性效应，在时域下模拟研究张力腿式风机的极限载荷和疲劳载荷。Nematbakhsh[4]认为张力腿张力主要受纵摇响应激励影响，且纵摇运动比纵荡、垂荡含有更多非线性成分。

目前较为广泛使用的分析软件，例如美国可再生能源中心开发的FAST，张力腿结构通常被简化为悬链线模型，这种简化虽然使计算简便，但是也导致针对张力腿结构的动力响应特性模拟不准确。截至目前为止，尚无广为认可的软件开发工具具备对张力腿式浮式风机系统进行复杂动力响应分析的能力。

本研究根据张力腿结构的结构和动力特性，利用自主开发的DARwind计算程序，采用谱方法对浮式风机系统的张力腿结构建立耦合动力模型，在时域下分析张力腿式浮式风机耦合动力响应。论文最后以UMAINE设计的张力腿式浮式风机为研究对象，将DARwind程序计算结果与试验数据对比，验证本文提出的理论模型和编制计算程序的可靠性。

1 张力腿耦合建模理论

张力腿式浮式风机是由风机、塔架、平台和张力腿等多个系统构成的复杂多刚体系统，各构件相互连接，采用多刚体模型能更好地反映各模块多自由度的运动特性。由于张力腿结构的特殊性，在建立浮式风机系统模型时，将张力腿锚泊系统同其它载荷分开建模，其余外载荷主要包括惯性和重力载荷、气动载荷及水动载荷。将锚泊系统单独建立数学模型，综合考虑张力腿与平台耦合处导缆孔的位移、速度、加速度，耦合建模求解，浮式风机可认为是锚泊力和其余外部环境载荷共同作用下的浮式系统。

近几年因为发展速度的不断加快，对于各种矿物资源的需求量也随之增加，这就导致出现了众多中小型矿山开采企业。这些企业的经济实力较弱，技术和设备落后，往往采用的是传统的技术和设备，导致各种安全事故发生。矿山开采企业重视PLC控制技术，可提高矿山开采电气设备的稳定性和安全性，降低开采成本，提高开采效率。

1.1 张力腿结构动力响应数学模型

  

图1 张力腿微段Fig.1 Tendon element

张力腿是轴向刚度很大的细长体，底端锚泊处位移为零，顶端导缆孔处随平台运动，相当于强迫位移。为了模拟张力腿动力响应，将张力腿基函数选取为一系列具有完备性的三角级数展开式的叠加，满足上述两端的边界条件，再用哈密顿原理求解非线性耦合常微分方程，在此基础上编程计算张力腿响应。通过三角级数展开也称为谱方法，计算精度较高，能较好地模拟张力腿张力的变化。此外，由于张力腿细长构件的特性，粘性力不可忽略，研究中采用Morison公式计算。

张力腿单根未伸长长度为L，均为垂直布置，z向即为构件中性轴伸长方向。如图1所示，水深为H时，张力腿l处微段上端点初始坐标为(X，Y，－H＋l)，受力状态下产生的位移为 (x，y，z)，则受力后坐标为 (X ＋ x，Y ＋ y， － H ＋ l＋ z)，张力方向为(xl，yl，1 ＋ zl)。

在l处取dl微段，受力后位移为(dx，dy，dz)，由于z向刚度很大，dz≪dx≈dy，dl近似沿z向，故

 

平台的立柱、浮筒满足L/D＞5的细长体假定，也需计算粘性力。将立柱、浮筒沿长度方向分割成多个截面，利用Morison公式(10)计算粘性力。

将应变式Taylor展开，得

 

重力为平台、塔架及风机的重力之和，惯性载荷则表现为广义惯性力。

 

从应变可以看出，张力腿形变伸长量不仅与轴向垂直位移有关，还同水平位移有关，类似于柔性体，动力响应呈现非线性，因而需要单独建立张力腿动力控制方程，研究张力腿动力响应。

张力腿建模时，利用谱方法，将受力状态下节点位移(x，y，z)看作强迫位移下一系列正弦级数的叠加，可以表达为

 

式中:x，y，z代表张力腿上沿轴向在l处选取各点的x，y，z三个方向上的位移；xf，yf，zf为张力腿顶端导缆孔处的强迫位移，根据平台参考点位移由方向余弦阵转换得出，由平台运动决定；l为张力腿轴向坐标，l＝0代表底端锚泊点，l＝L 代表顶端导缆孔系泊点； ai，bi，ci(i＝ 1，2，…，5) 为三角展开系数。

4）课文背景知识共享：教师引导学生分小组搜集和整理与课文内容相关的背景知识，并制作成PPT，在课堂上分组依次上台分享课文背景知识。由每组的小组长分配工作，组织所有组员们搜集信息，整理信息和上台展示成果。教师将根据每一位同学的表现给予评价、鼓励，并给予评分。通过这一任务活动，彻底改变教师直接灌输的传统方式，让学生主动参与教学、构建知识，共享自学信息。

水动力基于势流假定，计算包括入射势引起的傅汝德－克雷洛夫力、由绕射势引起的绕射力、由辐射势引起的辐射力以及静水回复力、细长构件的粘性力。利用WAMIT[8]计算出频域下的附加质量系数、阻尼系数、静水回复力系数和波浪激励力响应函数，转化为浮式风机系统受到的水动力，频域下计算公式为:

张力腿系统总势能由保守外力的势能及弹性势能组成。张力腿各分段的重力、浮力为保守外力，在对应z位移下产生广义的重力势能。考虑张力腿弹性势能时，张力腿为细长构件，轴向刚度很大，弯曲变形、扭转变形相较于拉压变形对势能影响很小，弯曲变形在中间段极小，仅在边界端内力有较大梯度，对整体外力并无影响，可以忽略不计。将应变表达式(4)代入系统的弹性势能，可以表示为:

 

式中:为非线性项，耦合编程时独立求解，对各广义坐标的偏导数在时间步上独立积分。

系统的动能可以表示为:

 

作用在张力腿上的外载荷主要包括水平方向的粘性力和惯性力，根据Morison公式，粘性力及惯性力可表示为:

 

式中:Uwj为波浪来流速度；uj为系统对应于j自由度的速度。

至此，张力腿系统的动能、势能及外载荷已表达完毕，可利用哈密顿(Hamilton)原理[5]组装方程，推导张力腿动力控制方程。哈密顿原理可表述为

 

式中:T是系统的动能；V是系统的总势能，包括应变能及任何保守外力的势能；Wde是作用在系统上的非保守力(包括阻尼力及任意外载荷)所做的功；δ是在指定时间区间内所取的变分。代入式(8)、(9)、(10)，离散简化得

 

 

以式(12)为例化简如下，由分部积分得

 

哈密顿原理中假定变分δai在积分限t1和t2处为0，所以将式(15)代入式(12)化简得

 

在RAO工况波浪激励下，平台进入稳态振荡后，激励运动幅值较小，相对速度也较小，使得与速度平方成正比的粘性力幅值很小，而RAO主要受到一阶波浪力和辐射阻尼的影响。从图7中可以看出:纵荡RAO数值计算结果同试验数据吻合度较好，DARwind数值计算在固有频率0.024 Hz处出现响应峰值，试验中由于无法模拟此低频部分，无法给出试验对比数据，但DARwind计算得出的固有频率0.024 Hz与UMAINE设计值完全一致。纵荡RAO计算中，张力腿粘性力对平台运动影响较小，主要是线性波浪力成分贡献较大，RAO吻合较好说明DARwind程序中一阶波浪力模拟计算较为准确。数值计算同试验结果的纵荡RAO趋势一致，但仍有一定误差，这是由于质量模型有偏差造成的，试验中传感器线缆布置对质量分布的干扰，使得RAO计算时的MCK控制方程中的M项并不准确，与数值计算会产生一定偏差。

 

将式(10)代入式(17)，整理得

 

式中:Uwx，Uwy为水质点速度在x，y方向分量；为非线性项，耦合编程对各广义坐标的偏导数在时间步上独立积分求解；

 

同理可得

根据式(18)、(19)、(20)建立了 ai，bi，ci控制方程，再由四阶 Runge-Kutta 数值迭代计算，可得出 ai，bi，ci以及应变ε，进而求得张力腿锚泊张力。

1.2 其余动力载荷模型

其余外部环境外载荷主要包括惯性载荷和重力载荷、空气动力载荷以及水动力载荷。详细的系统动力载荷模型参见文献[6]，这里仅作简单介绍。

1.2.1 惯性载荷和重力载荷

又zl≪1，二阶小量可忽略，故

1.2.2 空气动力载荷

系统所受气动力采用叶素动量理论(BEM)进行计算，BEM在一维叶轮动量理论和二维机翼理论基础上，又考虑了流场动量变化及叶片局部载荷，通过迭代计算相对风速、攻角以及轴向、切向诱导因子来计算各风速、叶片转速和偏角下的气动载荷[7]。叶片来流速度，诱导系数为:

 

 

  

图2 叶素分段受力图Fig.2 Local force on the blade element

式中:a是轴向诱导因子，a′是周向诱导因子，σ＝是旋翼自身的参数(B代表桨叶的数量，c(r)代表局部桨叶段弦长，r为该处与轮毂的距离)，ϕ是叶片旋转平面与来流速度的夹角，CL和CD是升力系数和阻力系数，ω是叶片旋转角速度。求解诱导系数是一个反复迭代的过程，直到收敛到一个合理值即可。叶素分段受力如图2所示，气动载荷由下式积分得到:

 

风机运行时，气动弹性变形、塔尖振动以及桨叶单元动力的非定常状态都需要被考虑进去，同时一些修正方法也被引入到BEM中来提高其准确性。对于无限数目桨叶风机，引入普朗特叶尖损失因子修正；对于高轴向诱导因子风机，引入Glauert修正；对于桨叶面偏航引入相应的偏航模型修正。

1.2.3 水动力载荷

在我国的部分地区，开展的民俗旅游，都出现了同质化现象，而造成这一现象的主要原因，则是开发者过度重视经济利益，而纷纷效仿别的地区的开发模式，并且将本地区原有的民俗文化活动摒弃，造成了本地区民俗文化的失真，而原本引进过来的民俗文化，又因与当地历史背景不符，从而产生了徒有其形的旅游项目，让慕名而来的旅者大失所望，影响了业界的口碑。

 

式中:fWj为对应于j自由度的波浪激励力；ϕI为入射波速度势；为对应于j自由度的绕射势；ρw为海水密度；ω为入射波频率；s0为浮式平台静水湿表面面积；μjk(∞)为对应于j自由度下无穷频的附加质量系数；Kjk(t－τ)为时延函数卷积项。

应变因为 z≪ x≈y≪ l，所以 zl≪ xl≈ yl≪l

2 时域耦合计算程序

张力腿浮式风机是典型的多刚体系统，本文建立起张力腿浮式风机的多刚体耦合计算模型，通过计算空气动力模块、水动力模块、张力腿等模块，得出对应的力及力矩，再利用凯恩方程[9]，计算整体系统运动方程求解系统动力响应的时域结果。详细的浮式风机系统动力响应计算参见文献[6]，这里仅简单介绍。

求解多刚体系统响应的凯恩动力学方程为

 

式中:为第l点处的主动内力(矩)元，为第l点处铰的相对平动速度对第j个广义速率uj的偏速度，为第k点处铰的相对转动角速度对第j个广义速率uj的偏角速度；为作用在各刚体上主动力的主矢，为刚体惯性力的主矢，为刚体惯性力对O点的主矩。

这里拥有众多古老的教堂和华丽的维多利亚时代的建筑，也拥有优美的自然风光。荷里路德宫作为爱丁堡的名片之一，营造了一个独特的人文语境，记述了苏格兰民族的悲怆与辉煌，为爱丁堡增添了无尽的人文魅力。

在多体系统运动学基础上，根据各构件连接关系递推计算参考点速度以及各构件角速度，计算每个构件的广义主动力和广义惯性力，并累加组成多刚体系统动力学方程。

计算张力腿锚泊张力时，采用谱方法通过两端边界条件建立正弦基函数，展开5项，计算张力腿的动能、势能，考虑了张力腿细长体的结构特性，由Morison公式计算了粘性力，代入哈密顿原理建立了张力腿运动的控制方程；最后采用四阶Runge-Kutta数值方法迭代求解张力腿中的张力载荷。

水利水电工程施工合同示范文本执行中的问题及对策……………………………………………………… 章 辉(4.32)

在时域下求解张力腿平台浮式风机系统耦合运动响应，计算流程如图3所示。在每一时刻，由上一时刻的整机系统运动方程，求解得出参考点的加速度，积分得到这一时刻的位移、速度，通过方向余弦阵转换为张力腿导缆孔处运动位移、速度，再由张力腿模块计算得到张力，作为外力传回整机系统运动方程计算这一时刻的加速度，反复迭代计算，直到获得下一时刻的系统动力响应参数。利用这个方法，计算出张力腿式浮式风机动力响应的时域结果。研究中根据图3所表达的流程，编制计算程序DARwind，完成张力腿浮式风机系统动力响应时域耦合计算。

  

图3 时域计算流程Fig.3 Time-domain simulation

  

图4 张力腿平台结构尺寸Fig.4 Structure parameters of the tension-leg platform

3 实验对比验证

使用美国可再生能源中心NREL和缅因州立大学UMAINE设计的张力腿式浮式风机作为研究对象，与实验数据对比，开展DARwind程序验证工作。试验中，采用傅汝德数相似准则按1∶50缩尺比制作模型，相应数据均已换算为实尺度结果。

3.1 浮式风机主要参数

研究对象采用 NREL提出的 5MW风机[10]和UMAINE设计的TLP平台[11]，张力腿平台上立柱及浮筒主要结构尺寸如图4所示。详细参数可参见文献[11]。

张力腿式浮式风机及支撑平台的主要参数见表1。

 

表1 张力腿式浮式风机主要参数Tab.1 Parameters of TLP wind turbine proposed by UMAINE

  

风机参数TLP平台参数额定功率/MW　5.0　整机总质量/kg　1 361 000风机质量/kg　350 000　排水量/kg　2 840 000桨叶长度/m　61.5　吃水/m　30.0轮毂半径/m　1.5　水深/m　200.0塔架质量/kg　249 718　张力腿直径/m　0.6塔架高度/m　77.6　张力腿刚度/MN　7 430.0

3.2 固有周期验证

为验证DARwind程序的可靠性，需要对平台运动固有周期及阻尼进行验证。文献[11-12]中给出了UMAINE张力腿平台固有周期的设计值与试验值，对比分析试验数据与DARwind计算结果。

以纵荡为例，选取参考文献[13]中的试验数据，在自由衰减工况下将数值计算结果试验数据进行对比。从图5可以看出，DARwind程序纵荡自由衰减在幅值和周期上与试验数据吻合较好，DARwind通过加入粘性阻尼综合考虑了立柱、浮筒、张力腿等细长体的粘性力。实验值阻尼相较于程序计算结果偏大，主要是由于DARwind中阻尼系数估量的不准确造成，由于TLP中细长体结构包含棱台等复杂构件，阻尼经验值在选取时稍显偏小，使得纵荡自由衰减同周期内幅值稍大，到达平衡状态历时偏长。

3.3 刚度验证

张力腿作为浮式平台的外载荷，轴向刚度很大，与平台运动耦合相关。平台导缆孔处的位移、速度、加速度影响张力腿动力控制方程的求解，而张力腿张力又影响平台的运动。为了验证张力腿对平台运动的影响，以刚度曲线的形式与试验进行对比。

  

图5 纵荡自由衰减对比Fig.5 Surge free decay comparison

  

图6 刚度对比Fig.6 Stiffness comparison

验证张力腿式浮式风机系统刚度时，强迫位移不一定满足垂荡和纵荡(set-down)之间的对应关系，所以取锚链由于强迫位移被拉伸的瞬态效应后的位移同张力关系，得到系统刚度曲线，与试验数据[11]对比如图6所示。图中DARwind程序同试验数据得到的刚度曲线都近似成直线，小位移下刚度呈现线性。张力腿轴向刚度很大，纵荡位移较小，从图中可以看出，刚度大致为85 kN/m，与试验曲线吻合较好，这说明张力腿上线性回复力模拟得较为准确，确保张力腿动态锚泊特性与试验相符。张力腿式浮式风机水平面内运动与刚度关联很大，小位移下DARwind计算结果同试验数据符合较好。

7.2.12根腐病①加强管理，保持树体生长健壮。多施有机肥，病树可施入赛众28微肥2～3公斤，使土壤透气性良好，根系在生长良好的条件下抗病力强。尽量不用山洪等污水灌溉。②加强早期落叶病、腐烂病的防治，避免引起树体衰弱，诱发根腐病。③对发病树及周围可被侵染的植株，用甲基硫菌灵800倍，硫酸铜300～400倍液灌根。连片发病时有必要挖深60厘米、宽20～30厘米的隔离沟，防止传播蔓延。

3.4 RAO工况验证

为了验证DARwind程序水动力性能计算的准确性，模拟试验[11]中的白噪声工况计算UMAINE TLP纵荡RAO曲线，数据对比结果如图7所示。

空间外差光谱仪捕获到观测物质的干涉图信息，利用干涉图与复原光谱间的傅里叶变换关系提取出目标光谱信号.传统的傅里叶变换法提取的光谱存在分辨率低、估计谱线浮动剧烈、混叠误差严重等缺陷[7-8].简小华等利用现代谱估计中的多重信号分类(Multiple Signal Classification，MUSIC)算法对传统傅里叶变换单色光谱信号进行提取，光谱复原效果很好，但在复色光谱复原方面还存在一系列问题，而且MUSIC算法对于信号空间维数p的取值十分敏感，虽然手工调试能获取最佳维数p值，但过程繁琐、耗时长，当需要处理数据量较大时，手工调试将无法完成工作 [9-15].

由于变分δai的任意性，所以仅当式(16)被积分项表达式之和为零时方程才能始终成立，即

  

图7 纵荡RAO对比Fig.7 Surge RAO comparison

4 结 语

针对张力腿动力特性，以导缆孔处位移、速度为边界条件，根据谱方法建立基函数自主编程，得到了张力腿时域耦合动力模型，利用DARwind程序计算动力响应。谱方法精度较高，在推导得出张力腿动力控制方程的基础上，DARwind程序能够较好模拟地张力腿动力特性，求解浮式风机的整机耦合响应。

以UMAINE TLP为研究对象开展动力响应数值计算，同试验数据进行了自由衰减、锚泊刚度、RAO等工况的对比，张力腿张力及浮式平台固有周期及运动特性都与实验吻合度较高，验证了程序的适用性，可用于张力腿式浮式风机的设计及响应分析。张力腿轴向刚度大，与平台纵摇等高频响应耦合振动，相较于将张力腿结构简化为准静态悬链线模型，动力耦合的张力腿控制模型将张力腿视为复杂运动下的动态结构，能更准确地描述张力腿同浮式风机间的耦合影响。

综上所述，通过阅读大量有关CT辐射防护研究资料，回顾分析笔者认为对CT辐射防护方法总结如下：（1）对CT扫描的剂量严格规定，确保参数设置的灵活性，综合患者的实际情况确保医师操作技术的准确数量，确保图像质量从而能够达到辐射剂量最低下，获取最高图像质量；（2）合理选择屏蔽方面并且积极研发新型屏蔽设备。下一阶段需要针对CT设备的全面辐射防护展开更加深入的研究。

参考文献:

[1] WAYMAN E N，SCLAVOUNOS P D，BUTTERFIELD S，et al.Coupled dynamic modeling of floating wind turbine systems:preprint[J].Wear.2006.

[2] YU M A，HU Z，XIAO L.Wind-wave induced dynamic response analysis for motions and mooring loads of a spar-type offshore floating wind turbine[J].Journal of Hydrodynamics， Ser.B.2015， 26(6):865-874.

[3] MATHA D.Model development and loads analysis of an offshore wind turbine on a tension leg platform with a comparison to other floating turbine concepts:April 2009[R].National Renewable Energy Laboratory(NREL)， Golden， CO.， 2010.

[4] NEMATBAKHSH A， BACHYNSKI E E， GAO Z， et al.Comparison of wave load effects on a TLP wind turbine by using computational fluid dynamics and potential flow theory approaches[J].Applied Ocean Research， 2015， 53:142-154.

[5] THOMSON W.Theory of vibration with applications[M].CRC Press，1996.

[6] 刘格梁.海上浮式风机系统动力响应的模拟与分析[D].上海:上海交通大学，2016.(LIU Geliang.Simulation and analysis for dynamic response of offshore floating wind turbine system[D].Shanghai:Shanghai Jiao Tong University，2016.(in Chinese))

[7] HANSEN M O.Aerodynamics of wind turbines[M].Routledge，2015.

除了事件本身，还可以通过其他细节让采访变得更加丰满。可以对采访对象所处的周边环境进行捕捉，给观众传递更多的信息；还可以依靠镜头语言，更好地展示采访对象的精神面貌，使观众提高对采访对象的认识，比如通过给采访对象的衣着、生活用品几个特写，观众就能够清楚知道采访对象的家庭等情况，提高观众的观看感受。

重型颅脑损伤是严重致残致死性疾病，国外报道占44岁以下成人死因首位，其治疗充满挑战，包括一系列集束化治疗方案，其中渗透性治疗占有举足轻重的地位[1]，临床工作中，静脉注射高渗药物成为治疗颅高压和脑疝综合征的常规治疗手段，然而，无论国内还是国外相关指南，均未对具体的药物选择作出推荐，理想的给药方式和剂量以及作用机制仍在探索之中，本研究旨在比较高渗盐水和甘露醇在重型颅脑损伤术后疗效、安全性方面的特点，为求最大程度减少数据结果偏倚，本研究剔除了合并原发脑干损伤病例。

[8] LEE C.WAMIT theory manual[M].Massachusetts Institute of Technology， Department of Ocean Engineering， 1995.

[9] KANE T R，LEVINSON D A.The use of Kane's dynamical equations in robotics[J].The International Journal of Robotics Research， 1983， 2(3):3-21.

[10]　JONKMAN J， BUTTERFIELD S， MUSIAL W， et al.Definition of a 5-MW reference wind turbine for offshore system development[J].National Renewable Energy Laboratory， Golden， CO， Technical Report No.NREL/TP-500-38060.2009.

[11]KOO B J，GOUPEE A J，KIMBALL R W，et al.Model tests for a floating wind turbine on three different floaters[J].Journal of Offshore Mechanics and Arctic Engineering， 2014， 136(2):20907.

[12]ROALD L， JONKMAN J， ROBERTSON A， et al.The effect of second-order hydrodynamics on floating offshore wind turbines[J].Energy Procedia， 2013， 35:253-264.

[13]STEWART G M，LACKNER M A，ROBERTSON A，et al.Calibration and validation of a FAST floating wind turbine model of the DeepCwind scaled tension-leg platform[C].International Society of Offshore and Polar Engineers.2012.

传统滨水动态人文景观分布的空间格局是以河涌水系为交通及景观廊道，以交通集汇、经济发达，休闲娱乐业兴盛的西关、长堤、许地、荔湾等为核心集聚区（面），构成了两轴、两核的整体格局（图8）：两条主轴一条是沿珠江前航道东西向的滨水景观轴线，一条是明清旧城墙范围内的传统城市中轴线；两核为传统滨水动态人文景观最为集中的两个区域：以祭祀水神的仁威神庙为精神核心节点的呈圆形扩散的西关功能核心区和以长堤滨江段为中线的呈椭圆形扩散的珠江两岸功能核心区。