一类光伏电站架构体系的风荷载特性及折减分析

_类光伏电站架构体系的风荷载特性及折减分析周炜;何斌;蔡晶;邓霞;刘丹【摘要】Thesolarphotovoltaicpowergenerationhasbecomeanimportantpartofthegreenenergydevelopmentandutilizationintheworld.Thecablestructuralsystemofphotovoltaicpowerstationhasbeenfocusedbecauseofmanyadvantages.Windcanbringaboutthemainloadactingonthephotovoltaicpanelsofphotovoltaicpowerstation.ThetwodimensionalComputationalFluidDynamics(CFD)numericalsimulationtechnologyisadoptedtoanalyzethewindloadcharacteristicsofarrayedphotovoltaicpanels.Fortwokindsofextremeconditions,whicharetheloadedcasesofforwardwindandbackwardwind,thewindloadanditsreductioncharacteristicsareanalyzedintheinvestigationconsideringtheeffectsofthephotovoltaicpanelangle,distancebetweenpanelsandheightofpanels.Theanalysismethodandcorrelativeconclusionarehelpfultothestructuralstrengthdesignofphotovoltaicpowerstation.%太阳能光伏发电已成为全球绿色能源开发利用的重要组成部分，拉索式架构体系光伏电站由于诸多优点备受关注.风是光伏电站架构体系中光伏板的重要荷载因素，采用二维计算流体动力学数值模拟技术研究阵列排放光伏板的风荷载特性,针对正向风载、背向风载两种极限工况，分析了光伏板倾角、板间距和板高度等因素对光伏板的风荷载及其折减特性的影响.相关分析方法和结论可为光伏电站架构体系的结构强度设计提供载荷依据.【期刊名称】《结构工程师》【年(卷)，期】2018(034)002【总页数】9页(P86-94)【关键词】光伏电站;拉索式架构体系;计算流体动力学;风载荷;折减系数【作者】周炜;何斌;蔡晶;邓霞;刘丹【作者单位】南京工业大学工程力学系，南京210009;南京工业大学工程力学系，南京210009浦京工业大学工程力学系浦京210009;羲和太阳能电力有限公司浦京210012;羲和太阳能电力有限公司，南京210012【正文语种】中文0引言太阳能是一种可再生能源，具有清洁无污染，取之不尽，用之不竭等优点，因此相关新型工程结构引起工程人员的重视［1-2］。光伏发电已成为全球绿色能源开发利用组成中的重要部分［3-4］,杨金焕等根据天空散射辐射各向异性的Hay模型，在倾斜面上辐射量的基础上推导得到了冬半年朝向赤道倾斜面最佳倾角的数学表达式［5-6］。丁明等运用粒子群算法求解以最大发电量和最大收益为目标的光伏阵列自身阴影损失的改进模型［7］。孙英云等在光伏阵列的隐互补模型基础上提出效用函数的光伏阵列I-U特性求解算法［8］。王晶晶等基于BSRN3000太阳辐射观测系统近3年的全年度太阳能气象实测数据，分析了呼和浩特地区倾角可调式太阳能光伏板表面接收的太阳辐射，得出了半年调整是该地区的光伏板倾角调整的最适合方案［9］。光伏板是光伏电站架构体系中直接承受风载荷的主要组件，风荷载是光伏电站架构体系中作用力最大的系统载荷,其研究手段主要围绕风洞试验和计算流体动力学(简称CFD)仿真两方面有效展开，风洞试验受到研究周期长、费用高等问题的限制不易实施,CFD仿真倍受工程技术及科研人员的青睐［10］。AlyMousaad等针对安装在具有两个倾斜面的屋顶太阳能构件风载荷分析，发现屋顶不同区域风载荷特性存在明显差异［11］。Cao等针对屋顶阵列光伏构件风载荷进行了系列风洞试验，寻找合适的建筑高度、屋顶结构和光伏板安装间距等参数［12］。牛斌等采用CFD方法分析光伏阵列所受风压，得到第一排光伏板为主要受力结构，其余阵列受首排保护受力较小［13］。阮辉等针对4.2m安装间距的5块光伏阵列，进行了表面风荷载数值的CFD仿真，发现在电池板中间截面受到较大正负荷载作用，而靠近电池板边壁处受到较小正负荷载作用［14］。该研究在以上研究基础上采用CFD仿真分析了光伏阵列的风荷载折减规律，以及光伏板的高度、倾角和间隔对折减规律的影响。CFX是全球第一个通过ISO9001质量认证的大型商业CFD软件,是英国AEATechnology公司为解决工业实际问题而开发，已用于光伏发电系统分析［13-14］。Standardk-8模型自从被Launder和Spalding提出之后，就变成工程流场计算中主要的工具,它是个半经验的公式，是从实验现象中总结出来的。RNGk-s模型在Standardk-s模型的基础上进行了改进，考虑到了湍流漩涡，提高这方面的精度，目前有钝体绕流的低速风场，湍流模型一般选用RNGk-s模型［15］。本文针对拉索式架构体系光伏电站结构，基于CFX计算平台，采用RNGk-s湍流模型，进行光伏板风荷载仿真计算，可为拉索式架构体系光伏电站架构体系的结构强度设计提供载荷依据。1控制方程在结构风工程中，结构物处于各梯度低速风场中，因此采用不可压缩黏性流体模型［16］。若不考虑流场中温度变化，则可认为流体的动力学黏度p和密度p为常数。在连续介质假设下,二维流场可通过下列基本方程描述［16］:连续性方程(1)运动方程(2)本构方程⑶将式(3)带入式(2),并考虑式(1),得到Navier-Stokes方程⑷式中：vi,vj是流体速度矢量;xi,xj是平面坐标;fi为单位质量流受到的体积力矢■；pji,pij为流体压力的二阶张量;&j是克罗内克符号;前面变量的下标i,j是自由坐标,可取1,2,分别表示x,y;v为流体的运动黏度系数;p为流体密度;p为流体压力的非静力分量;t为时间矛为动力学黏度系数。常温常压下(20。^1个标准大气压)，取空气的黏性系数p=1.7894x10-5kg/(m・s),密度p=1.225kg/m3,式(4)和式(1)组成了求解vi和p的基本方程。大涡模拟方法的主要思想是通过求解用滤波函数处理过的瞬时N-S方程来模拟大尺度涡,小尺度涡则通过求解附加方程得到［17］。在二维不可压缩湍流中，对于瞬时变量中'有(5)式中：是大尺度平均分量冲是小尺度分量。(6)式中：D是流动区域;x1,x2是实际流动区域的平面坐标;是滤波后大尺度平面上的平面坐标;是高斯滤波函数(i=1,2)。⑺式中：是方向的网格尺寸。通过方程式(6)、式(7)将瞬时的N-S方程和连续方程过滤,得到大涡模拟所使用的控制方程：(8)⑼式中：是滤波后的速度矢量;是滤波后压力的非静力分量;是常数;前面变量的下标i,j是自由坐标，可取12分别表示x,y［18］。Standardk-8和RNGk-8湍流模型都属于基于雷诺平均法的涡黏模型［19］。Standardk-£湍流模型主要考虑了湍流动能方程和湍动扩散率方程,并假定流场完全是湍流。RNGk-8和Standardk-8湍流模型很相似相比Standardk-&它修正了湍动黏度从而可以更好地处理湍流涡旋［15］。2力学建模2.1研究模型本文采用的研究对象是如图1所示的拉索式光伏架构体系，图2是其几何三维模型。因为受风载的大部分是光伏板，所以可以将模型简化主要受风载荷的光伏板。选择—种具有代表性的光伏板，其尺寸为1.65mx0.992mx0.02m,倾角为28°,前后距离为3.187m,距地面高度为3m,—排共12列光伏板，如图2所示。组件坐标系见图3,模型均以实际尺寸建立，所有计算模型阻塞比满足数值模拟计算要求，阻塞比表示模型中风道截面上的投影面积与风道截面积之比［20-22］。图1光伏架构体系Fig.1Photovoltaicstructuralsystem图2光伏架构体系三维模型图Fig.23DmodelofPhotovoltaicstructuralsystem图3光伏板模型示意图Fig.3Diagramofphotovoltaicpanels2.2网格划分与边界条件网格划分时把流场分成属于一个部件的三个部分，中间部分,流体的变化比较复杂，特别是靠近阵列板的部分,该区域采用膨胀网格，在板的长方向设置60个节点,宽方向设置3个节点，如图4所示。中间部分一共存在70万个节点，网格250万。上下两区域采用multizone网格，网格尺寸设置为0.2m,节点21.6万网格38万如图5所示。在单元的一些特定点上计算出雅可比矩阵行列式，最大值跟最小值的比率就是雅克比值，值为1最好，值越大就说明单元越扭曲。基于一个给定单元的体积与边长间的比率就是网格质量值。其值处于0和1之间,0为最差,1为最好。这里流场网格的雅克比平均值1.1,网格质量平均值0.7,网格质量符合要求。图4光伏板局部网格划分Fig.4Localgridpartitionofphotovoltaicpanels图5流场网格划分Fig.5Gridpartitionofflowfield边界条件设置，流体材料设为空气，密度为1.225kg/m3,入口处平均风速根据建筑设计规范，取十级风速20m/s。出口处静压为1个大气压。设置光伏板为固壁壁面，对近壁面区采用标准壁面函数法。由于计算模型为钝体低速绕流,RNGk-8湍流模型较为合适［15］。动量采用二阶迎风离散格式，压力和速度耦合采用SIMPLE算法［18］。3数值模拟与结果分析建立光伏板离地高度3m/顷角28°,板间距3.187m的CFD模型，计算得到z=0轴截面的正向工况和背向工况速度云图，如图6、图7所示。分析发现前两列光伏板的板前风速变化大,压力变化也会相应比较大;这是由于均匀来流受到迎风面第一列光伏板阻挡后，流场发生适应性变化;在第三列光伏板位置处高低方向稍微远离组件的位置流场趋于均匀，于是后续光伏板周围流场相对类似、变化较小。图8、图9分别为背向离地高度3m、板间距3.187m时光伏板倾角为28°和14°的流线图。从图中可以看出,倾角为28°时在第一、二列板之间会形成明显的涡流，而在倾角为14°时，在第一、二列板之间没有形成涡流,结构参数对流场特性影响很大工程设计中应充分考虑结构参数变化对风载荷的影响。下面将从工程应用的角度出发，特别是从结构体系强度设计等方面考虑，重点分析风载荷对整体某一光伏板的作用，分析不同安装高度的光伏板风荷载折减特性、不同安装倾角的光伏板风荷载折减特性、不同光伏板间距的光伏板风荷载折减特性等的影响。图6正向速度云图(Z=0截面，速度单位:m/s)Fig.6Velocityfieldcloudpictureofforwardwind(Z=0section,velocityUnit:m/s)图7背向速度云图(Z=0截面，速度单位:m/s)Fig.7Velocityfieldcloudpictureofbackwardwind(Z=0section,velocityUnit:m/s)图8背向倾角为28°前三列光伏板周围流线图Fig.8Streamlinediagramaroundfirstthreecolumnsphotovoltaicpanelswithaninclinationof28°图9背向倾角为14°前三列光伏板周围流线图Fig.9Streamlinediagramaroundthefirstthreecolumnsphotovoltaicpanelswithaninclinationof14°为了提取光伏板所受风阻对结构体系的影响，将光伏板视为刚体，于是作用于光伏板上的复杂分布空气动力可简化为力系的主矢和主矩。光伏板所受的载荷主要是垂直于光伏板的正压荷载和所受的倾覆作用力矩，其对光伏结构体系的承载有重要影响。这里将分布风阻向光伏板上距离地面最近的点简化得倾覆作用力矩和阻力载荷。这两个参量将成为下面分析的重点，由CFD计算可得光伏阵列各个光伏板的正压荷载和作用力矩，如图10、图11所示。观察图10和图11发现无论正向还是背向光伏板所受正压荷载都随流场的深入逐步衰减，且前两列光伏板所受风荷载变化最大，正向正压载荷和作用力矩折减达到50%,背向正压载荷折减达到30%,作用力矩折减达到60%。这是因为第一、二列光伏板之间和第二、三列光伏板之间的涡流特性影响了第二列光伏板的流场压力分布，使其整体受力出现显著的变化，详细结果见图12所示的前三列光伏板周围压力云图。后续光伏板风荷载单调折减，且最后几列光伏板风荷载折减率比较小，风荷载趋于稳定。背向载荷总体比正向载荷大,光伏板背向载荷是结构强度设计需要重点关注。图10正向背向光伏板所受正压载荷对比Fig.10Comparisonofpositivepressureloadsbetweenforwardwindandbackwardwind图11正向背向光伏板所受作用力矩对比Fig.11Comparisonoftorquebetweenforwardwindandbackwardwind图12背向前三列光伏板周围压力云图Fig.12Pressurefieldcloudpicturearoundfirstthreecolumnsphotovoltaicpanelsofbackwardwind3.1不同高度的光伏板风荷载折减特性保持光伏板倾角28°、板间距3.187m,光伏板的安装高度一般根据使用条件来确定，比如不能被遮挡，不能被水淹，一般离地高度大于0.5m,有农业耕种的光伏板离地高度3m左右。本文选取常用的光伏板离地高度0.5m、1m、2m、3m,分别建立CFD模型，分析光伏板离地高度对风荷载的影响。图13和图14为正向、背向第一列光伏板的所受的正压载荷和作用力矩。观察图13和图14发现正压载荷和作用力矩都随着光伏板高度的增大而减小，且背向所受正压载荷和作用力矩都比正向时大，且随着高度增加，正向、背向的正压载荷差距急剧变小。图13第一列光伏板所受正压载荷Fig.13Positivepressureoffirstcolumnphotovoltaicpanels图14第一列光伏板所受作用力矩Fig.14Torqueoffirstcolumnphotovoltaicpanels图15、图16为不同离地高度光伏板阵列所受作用力矩和正压载荷对比，分析发现不同高度的光伏板阵列折减趋势基本一致。离地高度越大，第一列到第二列的折减率越小，末尾几列的风荷载也越大一些。图15不同离地高度光伏板所受正压载荷对比Fig.15Comparisonofpositivepressureloadsindifferentheights图16不同离地高度光伏板所受作用力矩对比Fig.16Comparisonoftorqueindifferentheights3.2不同倾角的光伏板风荷载折减特性在保持光伏板离地高度3m、板间距3.187m不变的情况下，由于光伏板倾角和地理位置相关，本文选取具有代表性的4个倾角14°、21°、28°、35°,建立4个CFD模型。根据模拟计算结果，可以得到不同倾角正压载荷和作用力矩，如图17、图18所示。随着倾角增大，第一列光伏板的风荷载单调增大。倾角为14°,21。时，第二列的风荷载小于第三列，倾角为28°,35°时，第二列风荷载大于第三列，倾角的变化会影响第二、三列的光伏板风荷载变化规律。倾角对末尾的光伏板风载荷影响较小。除第一列光伏板外，倾角变化对其他光伏板受力影响复杂，大范围内难以发现较强规律性。图17不同倾角光伏板所受正压载荷对比Fig.17Comparisonofpositivepressureloadsatdifferentangles图18不同倾角光伏板所受作用力矩对比Fig.18Comparisonoftorqueatdifferentangles3.3不同间距的光伏板风荷载折减特性保持光伏板离地高度3m、光伏板倾角28°不变，光伏板间距可以根据光伏发电站设计规范(GB50797—2012)计算得到。本文选取该尺寸光伏板常用的间距2m、2.5m、3.187m、3.5m、4m,建立5个CFD模型。仿真可得不同间距情况下的正压载荷和作用力矩，如图19、图20所示。分析发现由于后续光伏板对绕过来流方向第一列光伏板尾流的遮挡作用，使得不同间距的光伏列阵，来流方向第一列光伏板风载荷也不相同。并且由于涡流的影响，第二列光伏板的风荷载比第一列和第三列的风荷载小。随着间距的增大，第二列板的风载荷快速增大,后续光伏板的风荷载折减趋势一致。间距对末尾的光伏板风载荷影响较小。图19不同间距光伏板所受正压载荷对比Fig.19Comparisonofpositivepressureloads>withdifferentdistances图20不同间距光伏板所受作用力矩对比Fig.20Comparisonoftorquewithdifferentdistances4结论本文针对光伏电站拉索式构架体系风荷载特性及其折减性能进行数值仿真，分析了正向风载、背向风载工况下，光伏板高度、光伏板倾角和光伏板间距对风荷载以及其折减特性的影响，结论如下：在数值上，背向风载荷大于正向风载荷。正向、背向两种工况下,第一列光伏板所受的正压载荷和作用力矩都随着高度增大而减小;且随着高度增加，正向、背向的正压载荷差距变小，第一列到第二列的折减率变小，末尾几列的风荷载变大。⑶倾角增大，第一列光伏板的风荷载也会随之增大。光伏板倾角的变化会影响第二三列的光伏板风荷载的变化规律。倾角对末尾的光伏板风载荷影响较小。⑷光伏板间距变化时，第一列光伏板的风荷载值不同;第二列板的风荷载比第一列和第三列的风荷载小;第二列板的风载荷随着间距增大而增大;后续的光伏板风荷载折减趋势一致。参考文献【相关文献】章红梅，董金芝,师振华，等.一种扣件+卡槽组合式连接的带光伏板剪力墙拟静力试验研究[J].结构工程师,2015,31(5):127-134.ZhangHongmei,DongJinzhi,ShiZhenhua,etal.ExperimentalstudyonpseudostaticresponseofaRCshearwallwithabuild-upfastening-grooveconnectedphotovoltaicpanel[J].StructuralEngineers,2015,31(5):127-134.(inChinese)梁亚伟,章红梅，师振华.卡式连接的带光伏板剪力墙拟静力与发电性能试验研究[J].结构工程师,2016,32(3):141-148.LiangYawei,ZhangHongmei,ShiZhenhua.LowcyclicloadingandpowergenerationperformancestudyonPV-shearwallintegratedsystembyaplug-inconnector[J].StructuralEngineer,2016,32(3):141-148.(inChinese)崔容强，汪建强,孟凡英，等.太阳能光伏发电之未来[J].可再生能源,2008,26(3):96-101.CuiRongqiang,WangJianqiang,MengFanying,etal.ThefutureofPVpowergeneration[J].RenewableEnergyResources,2008,26(3):96-101.(inChinese)曾祥光，殷俊，冯鉴.动车组太阳能光伏发电技术应用研究[J].可再生能源,2012,30(12):15-20.ZengXiangguang,YinJun,FengJian.ThestudyofsolarphotovoltaictechnologyappliedonEMU[J].RenewableEnergyResources,2012,30(12):15-20.(inChinese)杨金焕.固定式光伏方阵最佳倾角的分析[J].太阳能学报,1992,13(1):86-92.YangJinhuan.Analysisofoptimumtiltedangleforfixedphotovoltaicarray[J].ActaEnergiaeSolarisSinica,1992,13(1):86-92.(inChinese)杨金焕，毛家俊,陈中华.不同方位倾斜面上太阳辐射量及最佳倾角的计算[J].上海交通大学学报,2002,36(7):1032-1036.YangJinhuan,MaoJiajun,ChenZhonghua.Calculationofsolarradiationonvariouslyorientedtiltedsurfaceandoptimumtiltangle[J].JournalofShanghaiJiaotongUniversity,2002,36(7):1032-1036.(inChinese)丁明，刘盛，徐志成.光伏阵列改进优化设计方法与应用[J].中国电机工程学报,2013,33(34):1-8.DingMing,LiuSheng,XuZhicheng.Animprovedoptimizationdesignmethodforsolarphotovoltaicarraysandengineeringapplication[J].ProceedingsoftheCSEE,2013,33(34):1-8.(inChinese)孙英云，侯建兰，李润，等.基于隐互补问题的光伏阵列模型及其求解算法[J].中国电机工程学报,2014,34(34):6066-6073.SunYingyun,HouJianlan,LiRun,etal.Aphotovoltaicarraymodelbasedonimplicitcomplementarityproblemanditsalgorithm[J].ProceedingsoftheCSEE,2014,34(34):6066-6073.(inChinese)王晶晶，王志敏,冯志诚，等.基于实测数据的可调式太阳能光伏板倾角计算分析[J].可再生能源,2014,32(3):265-269.WangJingjing,WangZhimin,FengZhicheng,etal.CalculationandanalysisofadjustablePVpanelstiltanglebasedonmeasureddata[J].RenewableEnergyResources,2014,32(3):265-269.(inChinese)日本太阳光电发电协会.太阳能光伏发电系统的设计与施工[M].北京:科学出版社,2006.JapanphotovoltaicpowergenerationAssociation.Thedesignandconstructionofsolarphotovoltaicpowergenerationsystem[M].Beijing:SciencePress,2006.(inChinese)AlyMousaadAly,GirmaBitsuamlak.Wind-Inducedpressuresonsolarpanelsmountedonresidentialhomes[J].JournalofArchitecturalEngineering,2014,20:1-12.JinxinCao,AkihitoYoshida,ProshitKumarSaha,etal.Windloadingcharacteristicsofsolararraysmountedonflatroofs[J].JournalofWindEngineering&IndustrialAerodynamics,2013,123:214-225.牛斌，张超，侯巍，等.基于CFD方法的地面光伏阵列风压时程特性研究[J].太阳能学报,2016,37(7):1774-1779.NiuBin,ZhangChao,HouWei,etal.TimehistoryanalysisofwindloadonarrayedsolarpanelsbasedonCFDsimulations[J].ActaEnergiaeSolarisSinica,2016,37(7):1774-1779.(inChinese)阮辉，廖伟丽,王康生，等.光伏阵列表面风荷载数值研究[J].太阳能学报,2015,36(4):871-877.RuanHui,LiaoWeili,WangKangsheng,etal.NumericalresearchonsurfacewindloadofPVarray[J].ActaEnergiaeSolarisSinica,2015,36(4):871-877.(inChinese)黄张裕，阎虹旭.太阳能光伏板风荷载体型系数群体遮挡效应数值模拟研究[J].特种结构,2015,32(3):18-22.HuangZhangyu,YanHongxu.NumericalsimulationofthewindshapecoefficientsandblockingeffectaboutPVpanels[J].SpecialStructures,2015,32(3):18-22.(inChinese)李华峰，马骏，周岱，等.空间结构风场风载的数值模拟[J].上海交通大学学报,2006,40(12):2112-2117.LiHuafeng,MaJun,ZhouDai,etal.Thenumericalsimulationforthewindfieldsandwindloadsaroundspatialstructures[J].JournalofShanghaiJiaotongUniversity,2006,40(12):2112-2117.(inChinese)王文龙，何斌，傅杰，