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

user: 1. 引言 随着计算机技术，微处理器与电力电子技术的发展，逆变器控制系统也开始向数字化，智能化方向发展，并成为该领域的研究热点 [ 1 ]。SVPWM是近年发展的一种比较新颖的控制方法，是由三相功率逆变器的六个功率开关元件组成的特定开关模式产生的脉宽调制波，能够使输出电流波形尽可能接近于理想的正弦波形。空间电压矢量PWM与传统的正弦PWM不同，它是从三相输出电压的整体效果出发，着眼于如何使电机获得理想圆形磁链轨迹。SVPWM技术与SPWM相比较，绕组电流波形的谐波成分小，使得电机转矩脉动降低，旋转磁场更逼近圆形，而且使直流母线电压的利用率有了很大提高，且更易于实现数字化 [ 2 ]。 国内外学者对环流抑制提出了很多方法，硬件方面主要是从阻断环流通路考虑，如加入隔离变压器、增加硬件辅助电路等，但因提高了成本和并联系统的体积，已经鲜有人研究。软件方面，文献 [ 3 ] 提出了谐波消除PWM方法可以抑制零序环流的高频和低频分量，但会导致较高的开关损耗。文献 [ 4 ] 提出多载波PWM控制方案，通过废除零矢量的使用可以减小零序环流，但该方法对环流的抑制能力有限。文献 [ 5 ] 提出了非线性控制方法来抑制环流，但算法过于复杂，实际工程中应用困难。本文推导了T型三电平逆变器并联系统的平均模型，从零序环流的产生机理出发，对传统的SVPWM控制策略进行改进，设计出一种稳定的开环算法，并在SIMULINK中建立仿真模型验证了其可行性。
system: [SYSTEM]

user: 2. 逆变器并联系统数学模型 图1为两台T型三电平逆变器采用共直流母线和交流电网的结构并联组成的系统，选择直流负侧为参考点，则逆变器1的三相回路方程为(1)，其中 V a N x ， V b N x ， V c N x 是逆变器三相每个周期的平均电压；是电感电压；是电网中性点的电压。根据电感电压与电流的关系并考虑电感及线路内阻r，(1)可以改写为(2) { V a N x = V L a x + e a + V o N V b N x = V L b x + e b + V o N V c N x = V L c x + e c + V o N (1) ⇒ { L a x ⋅ d i a x d t = V a N x − r x ⋅ i a x − e a − V o N L b x ⋅ d i b x d t = V b N x − r x ⋅ i b x − e b − V o N L c x ⋅ d i a x d t = V c N x − r x ⋅ i c x − e c − V o N (2) 第x台逆变器的环流 i o x 为(3)，假定逆变器三相的电感值是相等的，令(2)的三个回路相加，并将(3)代入得 (4) i o x = 1 3 ∑ k = a ， b ， c i k x (3) L o x ⋅ d i o z d t = V o x − r x ⋅ i o z − 3 ⋅ V o N (4) 则两台逆变器零序电压的差值为 Δ V = V o 1 − V o 2 = ( L o 1 − L o 2 ) ⋅ d i o z d t + ( r 1 − r 2 ) ⋅ i o z (5) 如(5)所示，尽管零序环流的产生许多影响因素如硬件或者软件的不完美对称，并联逆变器的独立控制等，但最重要的还是并联逆变器的零序电压差异。由(5)可知零序环流的值与电网或者负载无关，主要取决于逆变器回路的零序阻抗，该阻抗为逆变器的电感内阻和线路内阻，因此非常小的零序电压差异将导致很大的环流。 图1. 并联T型三电平逆变器系统
system: [SYSTEM]

user: 3. 基于SVPWM的逆变器并联模型分析 3.1. SVPWM算法原理 三电平的空间电压矢量分布图如图2(a)所示。分为六大扇区，每个大扇区又分为四个小扇区 [ 6 ]。根据空间电压矢量的幅度，可以将其分为表1中的四组。小矢量是多余的，成对小矢量产生相同的相电压。 图2. (a)基本矢量分布图；(b)矢量对应零序电压图 Table 1 Vector classificatio 零矢量 PPP/OOO/NNN 小矢量 N型 OON/ONN/NON/NOO/NNO/ONO P型 POO/PPO/OPO/OPP/OOP/POP 中矢量 PON/OPN/NPO/NOP/ONP/PNO 大矢量 PNN/PPN/NPN/NPP/NNP/PNP 表1. 矢量分类 由于对称性，仅考虑扇区I就足够了，其他扇区都可以通过扇区I旋转一定角度得到。(6)所示的参考矢量Vref表示的交流侧电压可以在每个采样时刻由三个最近的矢量合成。 V r e f = 2 3 ( V a + V b ⋅ e j ⋅ 2 π 3 + V c ⋅ e j ⋅ 4 π 3 ) (6) 图3. (a) 参考矢量合成图；(b) 落入子扇区4的切换系列 当参考矢量落入I大扇区的4小区如图3(a)所示，根据SVPWM的矢量合成原则，参考矢量由距离最近的三个矢量合成，且由一个开关状态切换到另一个开关状态要求其开关管切换最少 [ 7 ]。则有(7) { d 0 V S 1 + d 1 V M + d 2 V L 1 = V r e f d 0 + d 1 + d 2 = 1 (7) 其中 d 0 ， d 1 ， d 2 分别是空间矢量 V S 1 ， V M ， V L 1 的占空比，一个典型的开关矢量分配图如图3(b)所示，本文均采用N型小矢量为起始小矢量，采用P型小矢量会得到相同的结果。 3.2. 由SVPWM引起的环流 在实际应用中，由于硬件或软件的不对称性，每个转换器d轴和q轴的电流输出控制不能保持一致，导致每个转换器的参考矢量有幅值或相位的差异 [ 8 ]。如图4所示，由于参考矢量相位或幅值的差异，参考矢量将会以不同时间、不同位置从子扇区2跨越到子扇区4，由于当参考矢量在子扇区2和4内，逆变器的零序电压值不同，当两台逆变器跨越扇区的时间不同，则两逆变器零序电压不同，导致环流。 假定逆变器的直流电容环节是平衡的，则逆变器的零序电压可以通过式(8)求出： V o = V a + V b + V c 6 (8) 图4. (a) 参考矢量同相位不同幅值；(b) 参考矢量同幅值不同相位 根据图2(b)所示的不同开关状态对应的零序电压值可以计算出当参考矢量落入I大扇区的子扇区4时，零序电压值为(9)同理，当参考电压落入I扇区的子扇区2时，其零序电压为(10) V 0 = V d c 3 ⋅ d 0 2 + V d c 2 ⋅ d 1 + 2 3 ⋅ d 2 + 5 V d c 6 ⋅ d 0 2 = 7 V d c 12 ⋅ d 0 − V d c 12 ⋅ d 1 + V d c 12 ⋅ d 2 (9) V 0 = V d c 6 ⋅ d 0 2 + V d c 3 ⋅ d 1 + 2 2 ⋅ d 2 + 2 V d c 3 ⋅ d 0 2 = 5 V d c 12 ⋅ d 0 − V d c 12 ⋅ d 1 + V d c 12 ⋅ d 2 (10) 参考矢量在各个扇区的零序电压值在每个大扇区内子扇区1-3的零序电压一致，且它们的起始小矢量相同，而4的零序电压与其他三个扇区不同，起始小矢量也不同于其他三个扇区。起始小矢量决定了参考矢量的合成序列，相同的起始小矢量在一个周期内产生相同的零序电压，因此结合之前的分析，零序电压的跳变主要是起始小矢量的异步变化引起的。 3.3. NP电压偏差导致的环流 上述分析都是基于直流电容环节是平衡的，但实际中，即使是加入了中点平衡策略，中点电位也是处于波动状态。这也导致了并联运行的三电平逆变器中点电位电压由偏差，导致环流产生 [ 9 ]。控制中点平衡的策略有很多，本文采用的是加入调节因子y调节成对正、负小矢量的持续时间来实现。并引入新的变量e，e定义为 { e = ( V c 1 − V c 2 ) / V d c − 1 ≤ e ≤ 1 (11) 则I扇区的四个子扇区的值如12-15所示，其他扇区的表达式是对称的。 V 01 = ( 5 12 − d 1 12 + d 2 12 + y ) V d c + e ( − 1 4 − d 1 12 − d 2 4 − y 3 ) V d c (12) (13) (14) (15)
system: [SYSTEM]

user: 4. 改进后的空间矢量调制算法 4.1. 抑制环流的周期性跳变 如3.2节分析所述，零序环流的周期性跳变主要是并联三电平逆变器之间起始小矢量的异步变化，参考电压之间幅值和相位的差异将引起该问题。为了抑制幅值差异导致的零序电压差异，可以通过将每个大扇区的子扇区1和2进一步分为两个区域：1a、1b和2a、2b，如图5(a)所示，当Vref从a区域移动到b区域时，起始小矢量发生改变，即零序电压发生变化，越过后则不再改变，这样可以抑制幅值差异引起的环流。为了消除相位差引起的环流，可以重新分配矢量图的扇区号，如图5(b)所示，I扇区变为[−30˚， 30˚]。改进后的扇区分布有如下好处： 1) 将每个扇区划为6个三角形区域，这样幅值差异将不会引起零序电压的周期性跳变。 2) 当有相位差异的并联逆变器的参考矢量落入同一扇区的不同子扇区时，零序电压差值不发生变化，但当参考电压由于相位差异异步地越过两相邻扇区的边界时，零序电压差值会发生跳变，导致环流的跳变。 3) 参考矢量扇区的判断一般是通过其角度和幅值来实现，角度决定了它所落入的大扇区，而通过大扇区的值结合幅值可以判断小扇区。本文通过将并联逆变器的大扇区判断值同步来实现相位差异引起的零序环流的跳变，大扇区同步之后，结合相应的幅值，再进行小扇区的判断，在跨越不同大扇区时候不再引入环流，但这种方法也有缺陷，当并联逆变器参考矢量相位差异过大时该方法将使影响从逆变器的期望输出。 图5. (a) 每个扇区划分6个扇区消除参考电压幅值差异引起的环流；(b) 新扇区分布，以消除由不同阶段引起的定期ZSCC跳跃 4.2. 抑制环流的恒定直流成分 环流中恒定的直流成分是由于并联逆变器的中点电位不一致导致的，为了消除该成分，提出了如图6所示中性点互联结构，使两逆变器的中点电位保持一致，消除了激励源，抑制环流。 图6. 模型控制框图 4.3. 抑制剩余的低频环流 如前所述，成对的正、负小矢量产生相同的线电压，但产生不同的零序电压，因此可以通过调整正负小矢量持续时间来抑制剩余的低频环流分量，但由于调整小矢量对中点电位也有影响，环流问题和中点平衡问题耦合使控制更加复杂。本文提出了一种开环的方法，无需解耦，就能同时实现二者的控制。假定有两台逆变器并联，选择其中一台为主逆变器，引入变量y 1 调整正、负小矢量持续时间来控制并联逆变器的中点电位平衡，获得y 1 的值后可以根据式(16)计算主逆变器的零序电压，并将其发送给逆变器2，逆变器2引入变量y 2 调整正、负小矢量持续时间进行零序电压的跟踪，逆变器2产生的零序电压同样根据式(16)计算，y 2 看为0，得到从逆变器的零序电压，这样计算两逆变器的零序电压差值，为了令该差值为零可以根据式(17)计算y 2 的所需值。这样就达到从逆变器对主逆变器的零序电压跟踪，从而抑制环流。综合以上三种方法即可对逆变器并联系统中的环流起到良好的抑制效果，且作为开环的算法在实际工程中更加稳定，无需添加任何硬件辅助电路，算法也只是对于传统SVPWM的改进，并不复杂，相较于一些复杂算法更易于实现。 (16) (17)
system: [SYSTEM]

user: 5. SIMULINK建模及仿真分析 5.1. 矢量控制系统仿真 Simulink是MATLAB中的一种可视化仿真工具，基于框图设计环境，可以实现系统的动态仿真。在此环境中，不必书写大量代码，而是通过简单连线操作即可实现复杂的仿真模型，图6为文章采用算法的控制模型框架。 5.2. 矢量控制系统仿真 为了确定理论分析的正确性和提出方法的有效性，在MATLAB中搭建了两台T型三电平系统的仿真，对上述所采取的环流控制策略进行验证，实验参数：并联逆变器台数n = 2，桥臂电感L = 1.37 mH，滤波电容Cf = 2 mF，开关频率f = 10 kHz，阻尼电阻R = 1 Ω。直流侧电压为650 V。 图7(a)为两逆变器参考电压幅值差20 V时逆变器的输出，其中ia1为逆变器1的a相输出电流，ia2为逆变器2的a相输出电流，iz为两逆变器间的环流，可以看到环流幅值高达50 A且使输出电流波形畸变。图7(b)是逆变器相位差异5˚时，a相输出电流及环流。可以看到环流幅值达到了40 A，且运行时间越长，输出电流波形畸变越严重。 图7. (a) 幅值差20 V时逆变器的输出及环流；(b) 相位差5˚时逆变器的输出及环流 图8(a)是两逆变器参考电压幅值相差20 V相位差5˚时，采取文章介绍控制算法，逆变器的a相输出电流及环流，可以看到环流相对于输出电流为极小的值，在图8(b)中放大显示，环流在±0.1 A附近波动，验证了本文控制算法的可行性。 图8. (a) 幅值差异20 V相位差5˚时的控制效果图；(b) 放大显示的环流
system: [SYSTEM]

user: 6. 结论 本文通过建立基于SVPWM的T型三电平逆变器并联系统的数学模型，分析了零序环流的产生机理，对传统的SVPWM控制策略进行改进，并在MATLAB SIMULINK中仿真验证了该策略的可行性，为实际数字控制系统的设计提供了理论依据。研究表明两台逆变器硬件参数和脉冲的不一致都会导致环流产生，在实际应用中是无法避免的，环流不仅增加系统损耗，还会使逆变器输出波形畸变，无环流控制的逆变器并联系统不能直接使用。而采取本文提出的环流抑制策略后，可以保证系统的稳定运行。相较于现有的环流抑制策略，本文文采用的方法无需增加硬件辅助电路，成本低，且是一种开环的控制方法，相较与闭环算法更加稳定，算法仅是从传统SVPWM改进而来，相较于一些复杂算法更易于实现，在实际工程应用中优势明显。
system: [SYSTEM]

user: 摘 要 为应对“互联网+”时代对房地产经纪行业的转型推动，本文在分析APP应用对房地产经纪行业的作用及影响基础上，从政府主管部门、房地产经纪行业协会、经纪机构及房地产经纪人几方面出发提出应对“互联网+”时代的对策。 关键词 :“互联网+”，线上交易平台，应用，影响，对策 Copyright © 2019 by author(s) 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. 前言 O2O (Online to Off-Line)模式，即线上到线下，互联网作为线下交易的基础，把线下的商务机会与互联网结合在一起。线下实体店的服务通过互联网在线上吸引消费者；消费者也可以通过互联网在线上筛选商家并支付选中的商品或服务，最终在实体店进行线下消费 [ 1 ]。 线下交易到O2O模式，是房地产经纪行业面对“互联网+”时代的及时反应，真实、快捷的信息发布能够降低二手房交易成本，但是平台发布的信息内容的完整性、现时性及平台沟通的通畅性关系到平台作用是否能够很好发挥，线上不是线下的补充，线下线上应平行工作并相互配合。目前杭州市中介行业协会登记在册会员单位共66家，笔者通过对比杭州房地产经纪市场占有量排名前三的我爱我家、链家和豪世华邦三家房地产经纪公司的APP应用，分析线上交易平台对房地产经纪行业的作用及影响。
system: [SYSTEM]

user: 2. 线上交易平台对房地产经纪行业的作用 2.1. 拓宽房源信息发布渠道 信息发布量及客户可及性是经纪人获取客户的关键，APP应用前，经纪人每月需要付费在公司官网外的专业网站发布房源信息，被动等待客户联系，信息发布带来的有效客户量不高。O2O模式下，APP成为最有效的房源信息发布端口，卖家委托房源，经纪人第一件事情就是在APP上登记房源基本信息并对外发布。我爱我家、链家和豪世华邦三家房地产经纪公司的APP近两年建设逐渐完善，功能增加、房源真实、经纪人线上反应迅速，极大地提高了APP的普及度，成为信息发布首选端口。 2.2. 提高房源情况咨询便利度 APP的应用改变了经纪人和客户沟通渠道，房源情况咨询从过去的电话或门店咨询变成了APP上直接沟通。我爱我家、链家和豪世华邦三家房地产经纪公司的APP上每个房源信息都对应有专门经纪人，客户可以针对感兴趣房源询问房源对应经纪人，经纪人也可以通过观察房源访客情况主动联系客户。经纪人甚至可以通过观察客户访问APP的痕迹分析客户是买家还是卖家，还可以通过客户浏览房源情况对比分析购房需求，从而提高房源匹配成功率和带看成功率。 2.3. 更新看房方式 VR看房是“互联网+”时代下应运而生的一项技术，改变了看房方式，也提高了带看成功率。调研的3家经纪公司APP中，除了华邦是视频看房外，链家和我爱我家均具备VR看房功能，其中链家VR看房有经纪人同步语音讲解和在线答疑，可以最大限度减少无效的房源带看。预约看房是APP相较于之前网站端口新增功能，客户根据自身合适时间预约看房，所有看房经历在平台上形成看房记录，看房记录中可以显示已看房源状态，便于经纪人分析客户购房需求，也便于客户跟踪感兴趣房源状态。
system: [SYSTEM]

user: 3. 线上交易平台对房地产经纪行业的影响 APP的使用给经纪行业带来了改变也提高了挑战，信息化方便了顾客，但对经纪人带来了更高的要求。 3.1. 工作模式变更 APP出现前，经纪行业和互联网已有紧密合作，经纪公司利用官网或者专业网站发布房源信息，但网站是单向的，仅是房源发布端口，顾客搜索到感兴趣的房源通过电话联系经纪人，从线上走到线下，网站房源更多的是吸引客户走到门店的工具，所以出现了假房源。 APP应用后，除了发布房源，VR看房、初步沟通、预约看房、已看房源状态变更等都可以在手机端完成，假房源少了，门店不再是高质量访客的首选，电话不再是经纪人与客户的唯一沟通渠道，经纪人的工作模式从网站发布房源、电话联系客户和门店接待客户，变成APP发布房源和接待客户、房源看好门店洽谈。 3.2. 工作内容增加 互联网最大的特点就是信息流通速度快、信息可达性高，APP的应用给客户带来便利的同时给经纪人带来了更多的工作内容。 APP使用前，新房源在网站端口发布后，经纪人主要工作是通过大量电话或者门店访客找到合适的客户进行带看，带看记录是否需要填写各公司要求不一，客户对带看房源感兴趣后经纪人撮合直到交易成功。 APP应用后，除了正常的线下接待、带看和撮合以外，经纪人增加了线上工作内容，在房源信息发布后要不断更新房源状态，例如房源带看情况；要时刻关注访客量，负责房源一旦有访客浏览要主动联系，并通过访客浏览痕迹分析访客需求，做好线上客户接待；若有客户留言要及时跟进，把潜在客户变成真实客户。 3.3. 工作质量要求提升 经纪行业的产品是服务，服务质量影响着行业地位，APP使用前，经纪人的服务对象主要是客户，服务通过接待体现，而接待方式主要是电话和门店。 APP使用后，特别是像贝壳找房这类面向全行业的APP应用后，经纪人的服务面除了客户以外，还有非同公司的同行，工作面变宽，从公司内部合作变成全行业合作；经纪人的接待方式除了电话和门店，更多的是线上接待，线上接待要求经纪人的服务具有更高的主动性、及时性和连贯性，因为互联网时代，反应速度是赢取客户的关键。
system: [SYSTEM]

user: 4. 房地产经纪行业应对“互联网+”时代的对策 4.1. 政府监管部门及行业协会 杭州市房地产经纪行业的政府监管部门是杭州市住房保障和房产管理局，行业协会是杭州市房地产中介行业协会(以下简称协会)，同时注册房地产经纪人又由浙江省房地产估价师与经纪人学会(以下简称学会)统一管理，为促进行业健康有序发展，交叉管理的同时更需要明确分工。 4.1.1. 牵头建立统一线上交易平台 20世纪30年代美国诞生了MLS (Multiple Listing System)制度，实质是由众多房地产经纪人联合起来共同建立一个信息资源共享和相互合作的系统，使客户委托的房地产买卖、租赁业务能够在较短的时间内完成，在MLS模式下，房地产经纪人真正实现了行业内的信息资源共享和合作共赢，极大缩短了房地产交易时间，节约了交易成本，提升了交易效率 [ 2 ]。 我国近些年也出现了很多类似于MLS的平台第三方，例如安居客、搜房、58同城等，但更多的作用是房源信息发布端口，且房源信息不在所有经纪机构间共享，因为对于经纪机构而言，房源信息是生存的根本。要借鉴MLS制度，实现行业内共享，必须要有一个各经纪机构都信任的第三方牵头建立平台，在交易限购的市场前提下，可以由杭州市房地产中介行业协会，或协会引导有意向和有能力企业建立平台，指导各经纪机构挖掘自身优势，在平台服务、房源寻找和客源寻找等方面加强合作，共同推进行业健康发展。 4.1.2. 明确房地产经纪人管理方 当前杭州市房地产经纪人出现双重管理现象，由中国房地产估价师与经纪人学会组织考试，通过的经纪人由学会管理并组织线上继续教育，主要考试内容是交易制度政策、职业导论、专业基础和业务操作，这类登记注册的经纪人实际在岗率不高，真正在岗的经纪人参加的是协会组织的从业人员职业培训和继续教育，并颁发专业培训证书，培训的内容主要是行业基础知识、交易基本法律常识、行为规范、职业素养、交易基本流程、合同条款解读及政策解读等。在岗的注册房地产经纪人既受学会的管理也受协会的管理，为避免重复学习，可以将经纪人培训及继续教育下放到协会，协会牵头会员单位建立行业内合作章程并引导业内合作，在原有培训内容基础上，增加合作章程学习、合作模式探讨，以实现合作方共赢及行业利益最大化。 4.1.3. 扩大房地产经纪人培训范围 目前协会组织的考证培训因培训力量有限导致覆盖面较窄，从业人员依据入职时间排序参加，培训考证不是上岗必须条件，经纪人参加公司组织入职培训即可上岗，这就造成经纪人因缺乏系统学习导致基础知识不扎实、法律意识和风险防范意识不强、职业素养不高且行为欠规范现象出现，为避免此类问题出现，协会可借助国家大力发展职业教育契机，在各职业院校设立培训点以扩大培训范围，将考证作为上岗前提，对每一位入行人员实施岗前培训。 4.1.4. 建立准入及退出机制 O2O模式在杭州乃至全国都还处在起步阶段，为保证行业健康发展，在借鉴国外运营模式的同时，需要行业政府监管部门，即杭州市住房保障和房产管理局加大监管力度，制定行业进入标准，公开投诉渠道，对于投诉率高的企业实施强制退出机制，特别是当下较热门的以租赁为主营业务的经纪机构，例如蛋壳公寓，没有实体店，线上交易线下消费，一旦出现问题，除了客服电话无任何渠道可以解决问题，在行业转型升级阶段，不能让客户产生投诉无门的感觉，这样将降低O2O模式的市场信任度，阻碍行业有序发展。 4.2. 房地产经纪机构 4.2.1. 加大APP投入，完善APP应用功能 从杭州市经纪行业来看，近两年各大经纪机构都加大了APP投入，功能逐渐完善。链家到贝壳找房体现了链家地产作为一家数据驱动的房产服务平台的初心；华邦从最开始将APP作为单纯房源发布端口，甚至房源信息都无法保证完整的情况下，发展成客户一旦浏览经纪人就主动跟进，反映了企业和经纪人对APP的重视和转型的决心。在接下来阶段，已有APP的企业为巩固发展成果，应加大对客户体验感受的关注，做好客户满意度测评和调查，从技术和服务两方面完善APP应用功能，例如华邦将视频看房更新为VR看房。不具备APP的企业可以考虑以加盟或者购买平台的方式实现O2O转型。 4.2.2. 合理安排在岗经纪人培训内容 各经纪机构在经纪人上岗前有对应岗位培训，从业过程中通过各类考试实现技能提升，例如链家有贝壳博学大考。在O2O模式下，经纪机构应重点加强经纪人线上服务能力，除了公司集中培训外，各门店可以依据自身情况顶定期开展具体化业务技能分享，例如房源撰写、卖点提炼、买家浏览痕迹数据分析及购房需求分析、房源客源匹配技巧、推介技巧、疑难案例分享等，由业绩好的经验丰富的经纪人分享从业经历和交易案例，形成良好工作氛围，实现共同进步。 4.2.3. 调整佣金及考核制度 佣金及考核制度是激励经纪人提升服务质量的重要手段，除了常规的依据交易价定佣金外，经纪机构可以借鉴杭州市房地产经纪行业管理服务平台客户点评功能，在APP上增加对经纪人评分功能，以服务阶段为划分依据，客户评价纳入月度考核，形成奖励性绩效，以鼓励经纪人提高服务主动性。 4.3. 房地产经纪人 4.3.1. 加强学习，提升业务技能 行业转型升级大背景下，市场、公司及客户对经纪人的要求都会提高，为了跟上行业发展步伐，经纪人应加强自身学习，除了协会及公司组织学习的基本职业素养外，要重点加强房源介绍能力，例如APP上房源卖点提炼；其次要往数据分析和数据处理方向努力，时刻关注APP上所负责房源动态，养成及时跟进客户习惯，依据客户浏览痕迹分析客户购房需求及偏好，以此加强房源客源匹配能力，提高带看成功率。 4.3.2. 提高服务意识，增强主动性 服务质量是经纪人在市场中取胜的关键，目前经纪人对买家服务质量相较于卖家而言更好，普遍认为卖家将房源委托给公司，公司将房源在合理价位以较快速度交易既是对卖家最大的服务，其实不然，经纪人服务质量提升空间还很大，例如在委托时为房东分析房源优缺点，指导定价；每天带看后及时向房东反馈房源推介和带看情况，让房东随时掌握销售进程；及时将市场行情变化反馈房东，给出调整价格或其他便于快速交易的建议等。
system: [SYSTEM]

user: 摘 要 随着我国对海洋渔业政策的不断调整，渔港的业态升级与空间的改造越来越显得必要及突出。而以渔港为基本依托的休闲渔业业态呈现出快速增长的需求，渔港的“游港化”改造成为客观必要。本文从渔港改造的当前外在背景以及“游港化”改造的内在需求出发，充分阐述并论证渔港的“游港化”改造的必要性，提出“游港”是一项综合性的全新业态，是渔区实施近海渔业减船减产、渔港推进渔业“双转”的替代性业态，更是地方渔业转型过程必由之路的观点。 关键词 :渔港“游港化”改造、滨海旅游、休闲渔业 Copyright © 2019 by author(s) 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.1. 关于渔港“游港化”定义 目前，国内对渔港的“游港化”改造的研究还处于起步阶段，对“游港化”的定义业界尚未形成统一的认识，笔者通过学习和借鉴国外及台湾地区渔港改造经历的成功案例，对渔港“游港化”的定义理解为：以市场需求为导向，为既有渔港引入休闲业态元素，并对其业态、空间、生态和文化的局部或整体进行针对性旅游利用的改造，通过兼用或专用港手段达到渔港功能多元化的目的，从而使经过改造的渔港更好发挥产业效益及社会服务效用。 1.2. 渔港“游港化”基本内涵 渔港“游港化”主要包括三个方面基本特征：第一，业态升级。渔港原有的传统业态，将当地渔业生产依托的渔港原有的码头修补、加工作业、船艇修造等生产方式通过节目化改造转化为渔业休闲体验项目，使区域局部或整体逐步转向休闲娱乐型的渔港业态；第二，元素转换。建设游艇公共泊位及休闲渔船专用及兼用港区，实现游客集散及始发管理与传统渔港的分区管理，把渔港打造成旅游消费的集聚地；第三，生态优化。在改造增设配套服务设施的同时，按照游港化要求改善渔港岸线景观，美化海岛渔村道路，活化当地的渔家生活方式及游艺方式，加深游客的深度体验 [ 1 ]。 1.3. 目前渔港存在的发展瓶颈 现阶段，一方面受到区位、信息、交通等外部环境的制约局限，另一方面也受渔港管理主体不明、社会资本难以介入、专门管理人才紧缺等内在因素困扰，渔港改造普遍存在发展瓶颈，具体可概括为以下几点： (1) 渔港现状业态过于单一，渔港配套产业存在起点不高等问题。仅以舟山市沈家门渔港为例，作为世界三大著名群众渔港之一的沈家门渔港中心港区业态空心化明显，渔港两岸约10公里岸线除了50余家海鲜排挡，公共专用区停泊功能缺失，岸线的游艇、休闲渔船的始发集散及文化创意服务延伸没有得到充分利用，全年500余万的中转游客依旧停留在渔港赏景，难以产生逗留消费。至于其他区域不同等级的渔港码头处于闲置及低水平利用的状况就更加普遍。 (2) 渔港改造受到现存渔港管理体制的局限，未能根据市场需求得以同步发展。由于我国目前的渔港责任主体与投资主体无法给以进一步明确，渔港岸线的分段设施利用、分项产品建设、监管及收益的分享等问题，很大程度受到产权界定的约束。看似谁都在管，结果往往是谁都可以不管，体制障碍严重阻滞了产业的健康发展 [ 2 ]。 (3) 渔港建设资金来源单一，拒绝社会资本进入的现状制约了渔港的快速发展。据了解，我国大中型渔港的改造建设资金通常是由中央、省、市各级政府共同出资，协商解决，这种资金配置方式一旦地方政府支付比例遇到资金紧张问题，就会严重制约渔港的整体建设完成及分块改造项目的如期完结。这种不可控的动力机制，往往缺失市场需求为导向的资本配置，会导致资源浪费及工程项目延期。
system: [SYSTEM]

user: 2. “游港化”改造外在背景 从渔港外在条件和背景出发，渔港的“游港化”改造是海洋资源环境压力的必然选择。主要表现在以下三个方面： 2.1. 过度捕捞导致海洋渔业资源枯竭 我国海域辽阔，渔业资源相对丰富。但上世纪八十年代中后期以来，一方面随着沿海渔船数量规模扩大，捕捞能力猛增；另一方面，不科学的捕捞方式也对海洋环境造成不可逆的破坏，导致渔业资源的消耗速度远远大于海洋资源的自然增长速度，海洋渔业资源趋于枯竭。据2018年全国渔业经济统计公报资料反映，目前国内仍拥有39.14万艘海洋捕捞渔船，国内近海捕捞产量达到490.48万吨，尽管以每年大于5%及以上的幅度缩减捕捞总量，但总的渔业资源萎缩趋势仍然没有得到有效遏制。 2.2. 渔业结构发生根本改变 据资料反映，2018年我国全社会渔业经济总产值达到25864.47亿元，其中渔业产值12815.41亿元，占比49.6%，捕捞及养殖业总趋势逐年下降；渔业经济三次产业产值的比例为49.6:21.9:28.5，渔业转方式调结构取得积极进展，三产占比快速上升 [ 3 ]。“十二五”以来，我国休闲渔业年均增长率达18.82%。2018年全国休闲渔业产值已经达到839.53亿元，占我国渔业经济总产值的3.25%。休闲渔业经营主体数量增长迅速，总量已达12.39万个，全国休闲渔业从业人员共80.49万人，同比增长17.86%，人均产值达10.43万元；全国共创建休闲渔业示范基地500余个，2018年全国共接待休闲渔业游客2.59亿人次，同比增长17.82% [ 4 ]。渔业五大产业均衡发展格局初步形成，渔港的游港化改造作为促进休闲渔业旅游发展的强劲动力，已然成为一种必然的理性选择。 2.3. 滨海旅游需求稳步增长 这是渔港改造的市场需求客观反映，也是“改造”行动的直接动力。海洋旅游正在逐步成为旅游业届的“宠儿”，备受旅游界的关注。以舟山市为例，2018年舟山市国民经济和社会发展公报数据显示，舟山市全年实现旅游总收入942.2亿元，比上年增长16.8%。旅游接待人数6321.4万人次，增长14.8%。数据表明全市已连续十年保持稳定的两位数增长，见图1。以休闲渔业为重要支撑的海洋旅游已然呈现出产品多选化、体验亲水化、消费大众化的倾向，以渔港为中心，以旅游消费为依托的发展潜力空间巨大。近年来，舟山市范围内海洋旅游企业及其他社会资本投资热情高涨，给逐步兴起的海洋旅游业发展注入了新鲜血液。据《舟山市休闲渔船发展规划》项目投资计划显示，到规划末期(2025年)休闲渔业及相关产业基础设施投资估算总额为16.43亿元，其中近期投入7.24亿元，占总投资的44.1%；中远期投入9.19亿元，占总投资的55.9%。休闲渔业综合体建设将引领产业先导 [ 5 ]。
system: [SYSTEM]

user: 3. “游港化”改造的内在需求 3.1. 休闲渔业业态发展必要 渔港是休闲渔业业态发展的重要基础。渔港的“游港化”改造主要基于渔港自身功能结构转型的主观需要。如充分利用渔港及周边水域原有的海水养殖筏排暂养设施、岛礁海钓资源、贝藻类养殖采集、休闲船艇运动等水域基础，通过产业整合、产品梳理引入“游港”休闲元素，推进培育“渔人码头”综合体、游艇游钓俱乐部、休闲船队等平台型等运营主体，重点转化原有的渔家采集、渔民捕捞、生产性垂钓等生产方式为休闲渔业旅游方式，有效缩减生产性捕捞能力，并使渔民进一步实现转业不离海，转岗不离船，减轻当地近海海洋资源环境压力，进一步促进渔区社会、经济的均衡、和谐发展。 图1. 2009~2018年旅游接待人数及其增速 3.2. 渔港产业空间调整需要 渔港的“游港化”改造可促使渔港有限空间得到更好利用。通过始发码头、集散中心、休闲渔业综合体、专吃一条街、专用港区等项目建设，使原渔港有限空间得到充分释放，提高渔港集聚功能效率。在当前推进渔业捕捞、加工和服务产业融合发展的势头下，一方面是传统渔港码头岸线及周边海域缺乏产业合理分工，空间布局不适宜，水域与岸线严重闲置；另一方面是目前很多好的旅游项目拿不到优质岸线，更申请不到水上运动休闲的专用水域或兼用水域，造成休闲渔业旅游功能与渔业水域的对立，导致渔港固有的海鲜夜排档、渔港文化体验、渔港夜游等延伸业态未能发挥最大的效用及效应。未来我国休闲渔业的发展，非常需要通过加强规范管理、提升产业发展层次和激发内生动力，促进休闲渔业高质量发展，而渔港的游港化改造过程正好能客观上满足这种发展的推进需求。 3.3. 城市文化及旅游品牌建设需要 渔港的“游港化”改造为渔港城市的市民游憩娱乐、旅游品质提升、文化传播与传承将奠定一定的发展基础。一方面通过完善休闲船艇码头活动中心、海钓俱乐部、赛事会展中心等活动基地建设，扩展了渔港为当地市民水上运动、传统渔家游艺赛事的活动空间，给当地城市人的日常休闲增加了新的服务设施提升了生活品质；另一方面也改变了之前因渔港设施简陋、服务低劣给游客造成的不满足、不满意的诸多困惑与难题，并极大丰富了旅游客源对垂钓、采集、环岛观光游览、水上运动休闲的节目选择；更重要的是以上的功能与展现，进一步改善了滨海旅游城市的对外形象，使传统渔港通过整合集聚资源要素，构建功能完善、产业融合的渔港生态系统，向全景渔港、人文渔港、休闲渔港、活力渔港转变，渔港及特色小镇的形象更加鲜明，这将在全国范围内对探索滨海及海岛地区构建渔港旅游样板具有重要创新意义。
system: [SYSTEM]

user: 4. 结论 目前我国休闲渔业总产值占渔业经济总产值的比重不到5%，与国际上休闲渔业发展水平较高的国家相比差距较大。根据近年来我国休闲渔业发展趋势，预测到2025年休闲渔业增加值达到2000亿元 [ 5 ]。我国的休闲渔业发展前景广阔，将拥有更加美好的未来。综上所述，“游港”是一项综合性的全新业态，是渔区实施近海渔业减船减产、渔港推进“双转”的替代性业态，也是地方渔业转型过程的必由之路，更是沿海地区传统渔港海洋旅游产品提升、渔港城市服务消费升级的载体与重要运营平台。以传统渔港引入“游港”元素，并在业态结构调整、空间形态重构、生态环境改善、传统文化彰显等方面展开全面的融合性改造，在满足传统渔业生产需要的同时，又能最大程度发挥休闲一类船艇及其附加服务的集聚及延伸功能。
system: [SYSTEM]

user: 摘 要 采用二维非定常不可压缩雷诺平均N-S方程结合RNG k-ε湍流模型，模拟分析了TSI (Technical Specifications for Interoperability)阵风以−5˚~+5˚攻角作用下箱梁断面的气动力及流场特性。TSI阵风作用过程中，箱梁断面的气动力总体上围绕准定常假设的气动力做小幅波动，但较大风攻角时箱梁断面的气动力波动幅度较大。分析探讨了箱梁断面的阻力系数跳跃现象、升力系数变化滞后于TSI阵风风速变化、流场静压分布的顺压和逆压梯度特性等。 关键词 :TSI阵风，箱梁断面，气动特性，数值模拟 Copyright © 2020 by author(s) 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. 引言 非定常风速的模拟通常有以下3种方法：1) 阵风系数法将非定常风速简化为定常风速，采用风速峰值计算风荷载。对于动力响应不敏感的刚性桥梁结构可采用该方法。如，公路桥梁抗风设计中采用静阵风系数计算±3˚风攻角内桥梁主梁的横向最大静阵风荷载 [ 1 ] [ 2 ]。2) 频域方法基于风谱表达非定常的风速，然后在时域进行动力响应分析。该方法能够充分考虑结构与动力响应间的相关性 [ 3 ]。3) 理想阵风法采用时间相关确定阵风时程，采用积分法计算系统的动力响应。该方法可以快速获得系统的最大响应，不需要进行复杂的分析计算。TSI (Technical Specifications for Interoperability)阵风模型 [ 4 ] 就是其中一种，可以有效近似风速最大值邻域的随机过程，常用于横向风作用下列车的脱轨分析 [ 5 ] [ 6 ]、带风屏障的高架桥上的高速列车的气动特性分析 [ 7 ]。在这些应用研究中，均基于准定常假设，直接应用静力三分力系数计算TSI阵风作用下列车的气动力。本文基于冻结湍流假设(frozen turbulence hypothesis) [ 8 ] [ 9 ]，应用FLUENT软件，采用CFD (Computational Fluid Dynamics)方法研究TSI阵风作用过程中箱梁断面的气动力及流场特性。
system: [SYSTEM]

user: 2. 物理模型 2.1. 阵风模型 本文采用仅考虑确定性阵风部分TSI阵风模型 [ 4 ] [ 7 ]。阵风风速 V ( t ) 定义： V ( t ) = β ( t ) U (1) 式中，U为相应的均匀来流风速，取 U = 10   m / s ； β ( t ) 为考虑阵风效应的系数： β ( t ) = { A + ( A − 1 ) sin ( π t + 0.5 π )                                   1 ≤ t ≤ 1.5 A + ( A − 1 ) sin ( π t − 0.5 π )                                   1.5 ≤ t ≤ 2 2 − A − ( A − 1 ) sin ( π t − 1.5 π )                             3 ≤ t ≤ 3.5 2 − A + ( A − 1 ) sin ( π t + 2.5 π )                         3.5 ≤ t ≤ 4 1.0                                                                                             t   取 其 它 值 (2) 式中，A为静阵风风速振幅，根据文献 [ 4 ] [ 7 ] 取1.7。 2.2. 箱梁模型及计算域 参考文献 [ 10 ] 中的桥梁模型，箱梁断面如图1所示。风攻角采用 α = − 5 ∘ ~ + 5 ∘ 。CFD计算时取缩尺比1:40。为了防止箱梁断面附近的分离旋涡由外边界反射回流场、外边界附近的流场参数分布与外边界条件不相容，计算域(如图2)取20B × 30B的矩形区域(B为箱梁断面的计算宽度，B = 0.30 m)，以确保计算结果收敛。其中，入口及上/下侧边界距离箱梁断面形心10B、出口边界距离箱梁断面形心20B。计算域边界定义：入口边界为速度入口(velocity-inlet)、上下边界为对称边界(symmetry)、出口边界为压力出口(pressure-outlet)、箱梁表面为无滑移壁面(wall)。自定义函数确定入口TSI阵风风速。 图1. 箱梁断面(单位：cm) 图2. 数值模拟计算域 2.3. 气动力 定义风轴坐标系下均匀来流风速U以风攻角 α 作用时的箱梁断面的阻力系数 C D ( α ) 、升力系数 C L ( α ) 和力矩系数 C M ( α ) ： C D ( α ) = F D ( α ) / ( 0.5 ρ U 2 B L ) (3a) C L ( α ) = F L ( α ) / ( 0.5 ρ U 2 B L ) (3b) C M ( α ) = M ( α ) / ( 0.5 ρ U 2 B 2 L ) (3c) 式中，B为箱梁断面的计算宽度，取 B = 0.30   m ；L为箱梁的节段长度，此处为二维断面，取 L = 1.0   m ； ρ 为空气密度，取 ρ = 1.225   kg / m 3 ； F D ( α ) 、 F L ( α ) 和 M ( α ) 分别为风轴坐标系下作用在箱梁断面上的阻力、升力和力矩(箱梁断面形心为扭转中心) (如图1)，均取非定常计算结果的平均值。 定义风轴坐标系下TSI阵风以风攻角 α 作用时箱梁断面的阻力系数 C D ， G ( α ， t ) 、升力系数 C L ， G ( α ， t ) 和力矩系数 C M ， G ( α ， t ) ： C D ， G ( α ， t ) = F D ， G ( α ， t ) / ( 0.5 ρ U 2 B L ) (4a) C L ， G ( α ， t ) = F L ， G ( α ， t ) / ( 0.5 ρ U 2 B L ) (4b) C M ， G ( α ， t ) = M G ( α ， t ) / ( 0.5 ρ U 2 B 2 L ) (4c) 式中， F D ， G ( α ， t ) 、 F L ， G ( α ， t ) 和 M G ( α ， t ) 分别为t时刻箱梁断面的阻力、升力和力矩，其正方向同图1所示的均匀来流作用下的静力三分力。 根据准定常假设，TSI阵风作用下箱梁断面的阻力系数、升力系数和力矩系数的参考值 C D ， G 0 ( α ， t ) 、 C L ， G 0 ( α ， t ) 和 C M ， G 0 ( α ， t ) 与均匀来流作用下的静力三分力系数存在如下关系： C D ， G 0 ( α ， t ) = β 2 ( t ) C D ( α ) (5a) C L ， G 0 ( α ， t ) = β 2 ( t ) C L ( α ) (5b) C M ， G 0 ( α ， t ) = β 2 ( t ) C M ( α ) (5c)
system: [SYSTEM]

user: 3. 数学模型 3.1. 控制方程 采用有限体积法求解基于Boussineq涡黏性模型的二维非定常不可压缩雷诺平均N-S方程： ∂ u ∂ x + ∂ v ∂ y = 0 (6a) ρ ( ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y ) = − ∂ p ˜ ∂ x + ∂ ∂ x ( μ ef ∂ u ∂ x ) + ∂ ∂ y ( μ ef ∂ u ∂ y ) (6b) ρ ( ∂ v ∂ t + u ∂ v ∂ x + v ∂ v ∂ y ) = − ∂ p ˜ ∂ x + ∂ ∂ x ( μ ef ∂ v ∂ x ) + ∂ ∂ y ( μ ef ∂ v ∂ y ) (6c) 式中， ρ 为空气密度，取 ρ = 1.225   kg / m 3 ；u和v分别为流体沿x和y坐标轴方向的速度；t为计算时间； p ˜ = p + 2 ρ k / 3 ，p为空气静压，k为湍流动能； μ ef = μ + μ t ， μ 为空气的动力黏性系数，取 μ = 1.7894 × 10 − 5 kg / ( m ⋅ s ) ， μ t 为空气的湍流黏性系数。 k和 μ t 由RNG k-ε湍流模型的运输方程 [ 11 ] 确定： ρ ( ∂ k ∂ t + u ∂ k ∂ x + v ∂ k ∂ y ) = ∂ ∂ x ( α k μ ef ∂ k ∂ x ) + ∂ ∂ y ( α k μ ef ∂ k ∂ y ) + S k (7a) ρ ( ∂ ε ∂ t + u ∂ ε ∂ x + v ∂ ε ∂ y ) = ∂ ∂ x ( α ε μ ef ∂ ε ∂ x ) + ∂ ∂ y ( α ε μ ef ∂ ε ∂ y ) + S ε − R ε (7c) 式中， α k 和 α ε 分别代表k和 ε 的逆有效Prandtl常数； R ε 代表平均应变率对湍流耗散率 ε 的影响； S k 、 S ε 和 R ε 的详细定义参见文献 [ 11 ] [ 12 ]。高雷诺数( μ ≪ μ t )下， μ t = ρ C μ k 2 / ε 。本文采用的湍流模型常数： C 1 ε = 1.42 ， C 2 ε = 1.68 ， C υ = 100 ， α 0 = 1.0 ， η 0 = 4.38 ， β ≈ 0.012 ， C μ = 0.0845 。 3.2. 计算网格及时间步长 以风攻角 α = + 4 ∘ 为例，由图3所示的3种网格计算均匀来流风速 U = 10   m / s 作用下箱梁断面的静力三分力系数，对比分析计算结果后确定TSI阵风分析的计算网格。采用四边形结构化网格对计算域进行网格划分，其中贴体网格高度根据文献 [ 12 ] [ 13 ] 确定。网格1的总网格数量为62，828，贴体网格高度为1.0 mm，计算得到 y + ≤ 13 。网格2在网格1基础上对尾流区域进行了网格加密，总网格数量为64，351。网格3进一步将贴体网格高度减小到0.5 mm，并对箱梁表面附近区域进行了网格加密，总网格数量为109，060。采用迭代时间步长0.0001 s计算得到的3种网格下的箱梁断面的静力三分力系数差别不大，与风洞试验结果 [ 10 ] 对比如表1所示。网格1可以满足计算精度要求。因此，采用网格1的方式对其余风攻角下的计算域进行网格划分。 图3. 3种网格局部图 Table 1 Comparison of the CFD results and wind tunnel test data of static coefficients of aerodynamic force 三分力系数 试验结果 [ 10 ] 网格1 网格2 网格3 CFD计算 相对误差 CFD计算 相对误差 CFD计算 相对误差 C D 0.4127 0.4793 16.14% 0.4822 16.84% 0.4793 16.14% C L 0.4365 0.4814 10.29% 0.5025 15.12% 0.5266 20.64% C M 0.1157 0.1400 21.00% 0.1388 19.97% 0.1393 20.40% 表1. CFD静力三分力系数与风洞试验结果对比 进一步采用网格1，分别以迭代时间步长0.00005 s和0.0005 s计算箱梁断面的静力三分力系数。计算结果及相对于风洞试验结果的误差如表2所示。可以看到，将迭代时间步长增大到0.0005 s后，CFD数值计算的精度并没有受到明显的影响。因此，本文的CFD计算均采用迭代时间步长0.0005 s。 Table 2 Numerical results of static coefficients of aerodynamic forces at different iteration time step 迭代时间步长 C D C L C M CFD计算 相对误差 CFD计算 相对误差 CFD计算 相对误差 0.00005 s 0.4791 16.09% 0.4691 7.47% 0.1394 20.48% 0.0001 s 0.4793 16.14% 0.4814 10.29% 0.1400 21.00% 0.0005 s 0.4796 16.21% 0.4823 10.49% 0.1401 21.09% 表2. 不同迭代时间步长的静力三分力系数计算结果
system: [SYSTEM]

user: 4. 箱梁断面的静力三分力系数验证 均匀来流风速 U = 10   m / s 以风攻角 α = − 5 ∘ ~ + 5 ∘ 作用，计算得到箱梁断面的静力三分力系数总体上大于风洞试验结果 [ 10 ]，相对误差20%左右，如表3所示。计算误差可能来源于箱梁表面附近的气流分离和再附，以及旋涡的脱落使得流场比较复杂，影响了CFD数值计算的精度。RNG k-ε湍流模型及计算条件也限制了计算精度。CFD计算结果总体上满足箱梁断面气动力的分析要求。 Table 3 Comparison of CFD results and wind tunnel data of static coefficients of aerodynamic force 风攻角(度) C D C L C M CFD计算 风洞试验 [ 10 ] 相对误差 CFD计算 风洞试验 [ 10 ] 相对误差 CFD计算 风洞试验 [ 10 ] 相对误差 −5 0.4206 0.3290 27.84% 0.4839 0.4914 −1.53% 0.2271 0.1798 26.31% −4 0.3913 0.3279 19.34% 0.6406 0.5820 10.07% 0.2259 0.1807 25.01% −3 0.3931 0.3298 19.19% 0.8009 0.6438 24.40% 0.2193 0.1691 29.69% −2 0.3778 0.3413 10.69% 0.8139 0.6366 27.85% 0.1733 0.1543 12.31% −1 0.4047 0.3574 13.23% 0.7659 0.6041 26.78% 0.1625 0.1410 15.25% 0 0.4231 0.3709 14.07% 0.7060 0.5584 26.43% 0.1611 0.1298 24.11% 1 0.4340 0.3803 14.12% 0.6487 0.5032 28.91% 0.1461 0.1211 20.64% 2 0.4553 0.3920 16.15% 0.5882 0.4712 24.83% 0.1423 0.1173 21.31% 3 0.4714 0.4042 16.63% 0.5221 0.4513 15.69% 0.1378 0.1163 18.49% 4 0.4796 0.4127 16.21% 0.4823 0.4365 10.49% 0.1401 0.1157 21.09% 5 0.4873 0.4251 14.63% 0.4607 0.4109 12.12% 0.1386 0.1153 20.21% 表3. CFD计算的静力三分力系数与风洞试验结果的对比
system: [SYSTEM]

user: 5. TSI阵风作用下箱梁断面的气动力及流场特性 5.1. 箱梁断面的气动力特性 当TSI阵风以风攻角 α = − 5 ∘ ~ + 5 ∘ 作用时，CFD计算得到的阻力系数 C D ， G ( α ， t ) 、升力系数 C L ， G ( α ， t ) 和力矩系数 C M ， G ( α ， t ) 总体上分别围绕参考值 C D ， G 0 ( α ， t ) 、 C L ， G 0 ( α ， t ) 和 C M ， G 0 ( α ， t ) 小幅波动，但较大风攻角时三分力系数的波动幅度较大。以 α = 0 ∘ ， ± 2 ∘ ， − 4 ∘ ， ± 5 ∘ 风攻角为例，CFD计算得到的阻力系数、升力系数和力矩系数时程曲线分别如图4~6所示。TSI阵风作用下箱梁断面的气动力有以下基本特性。1) TSI阵风作用下的阻力与准定常假设的阻力基本一致，但负攻角较大时阻力系数波动幅度较大，如图4(c)所示。在风速的峰/谷值时刻 t = 1.5   s ， 3.5   s 后，阻力系数时程曲线 C D ， G ( α ， t ) 存在明显的跳跃现象，如图4所示。其原因在于 t = 1.5   s ， 3.5   s 时刻风速的导数 d V / d t 的理论值不存在。理论上，CFD迭代计算的时间步长越小，这种跳跃现象越明显。当风速 V ( t ) < U 且风速减小到风速谷值附近时，存在轻微的负阻力现象，如图4(b)和图4(d)所示。2) 正/负风攻角较大时，阻力和升力系数时程的波动幅度较大。在风速 V ( t ) > U 的区段，即 t ∈ [ 1 ， 2 ] s ，攻角 α = − 4 ∘ ， − 5 ∘ 的阻力系数 C D ， G ( α ， t ) 和攻角 α = − 4 ∘ ， ± 5 ∘ 的升力系数 的波动幅度较大。但在风速 的区段，即 ，三分力系数时程曲线在风速谷值附近均有波动幅度减小趋势。3) TSI阵风作用下升力系数 变化滞后于阵风风速变化。在 的区段，风速由谷值逐步增大的初始阶段，升力系数 明显小于 ；在风速上升的末段，升力系数 相对于 有一定程度的偏离，如图5(b)、图5(d)。4) 风攻角较小时，在风速 区段的风速上升末段力矩系数 明显大于 ，如图6(b)、图6(d)。 图4. TSI阵风作用下的阻力系数时程曲线 图5. TSI阵风作用下的升力系数时程曲线 5.2. 箱梁断面的绕流流场特性 以风攻角 为例，TSI阵风作用下计算得到的流场静压分布如图7 (仅显示箱梁附近计算区域 ， )所示。箱梁断面的气动力计算仅考虑箱梁表面的静压的相对值，故流 场入口的初始静压设置为零，出口设置为压力出口并且固定静压为0 Pa。在风速大于均匀来流的时段 ( )，箱梁表面的远场静压随计算时间增大而增大。当入口风速由10 m/s逐步增大到最大值17 m/s时( )，CFD计算得到的入口静压由0 Pa逐渐增大到约245 Pa。当入口风速达到最大值 17 m/s后，再计算1个迭代时间步，入口静压急剧下跌到约−480 Pa；然后，入口风速逐步减小到10 m/s 时( )，入口静压由约−238 Pa逐步增大到0 Pa。入口静压的跳跃正是前述箱梁断面的阻力系数在 跳跃的原因。在阵风风速小于均匀来流的时段( )，箱梁表面的远场静压随计算时 间增大而减小。当入口风速由10 m/s逐步减小到最小值3 m/s时，CFD计算得到的入口静压由0 Pa逐渐减小到约−242 Pa。当入口风速达到最小值 3 m/s后，再计算1个迭代时间步，入口静压急剧增加到约482 Pa； 然后，入口风速由最小值3 m/s逐步增大到10 m/s时( )时，入口静压由约242 Pa逐步减小到 0 Pa。此处的入口静压跳跃也是前述箱梁断面的阻力系数在 跳跃的原因。 图6. TSI阵风作用下的力矩系数曲线 图7. 0˚风攻角静压分布(单位：Pa) 在风速增大时段( ， )，从入口到出口的流场静压整体上由正压逐渐减少到0 Pa，形成顺压梯度 ，造成箱梁表面的阻力增加；在风速减小时段( ， )， 从入口到出口的流场静压整体上由负压逐渐增加0 Pa，形成逆压梯度 ，造成箱梁表面的阻力减小。流场静压的总体特点与N-S方程的分析结果一致，但入口静压在TSI阵风风速峰/谷值时刻发生跳跃 可能不真实。在风速大于均匀来流风速的时段( )，TSI阵风使得箱梁的尾流旋涡增强了；而在阵风风速小于均匀来流风速的时段( )，TSI阵风使得箱梁的尾流旋涡减弱了。其余风攻角 下的静压分布的计算结果类似，除风攻角绝对值较大时，箱梁的尾流旋涡更强外。 箱梁尾流旋涡分布特点可以进一步从图8所示的0˚和+5˚风攻角下的速度分布对比图中看到。 图8. 0˚和+5˚攻角在 时的速度分布对比(单位：m/s)
system: [SYSTEM]

user: 6. 结论 本文采用RNG k-ε湍流模型模拟分析了TSI阵风以不同风攻角作用过程中箱梁断面的气动力及流场特性。分析得到以下主要结论：1) TSI阵风作用过程中箱梁断面的气动力总体上围绕准定常假设的气动力作小幅波动，但风攻角较大时气动力波动幅度较大。2) CFD计算得到的TSI阵风作用过程中箱梁断面的阻力系数在TSI阵风风速峰/谷值时刻存在跳跃、升力系数变化滞后于TSI阵风风速变化、流场整体静压存在顺压和逆压梯度等。TSI阵风模型应用于箱梁的气动力数值模拟存在局限性，将在后续的CFD数值模拟中改进。