大水深离岸型海上测风塔整体有限元分析
- 格式:pdf
- 大小:311.40 KB
- 文档页数:5
- 15 -高 新 技 术0 引言对比陆上风力发电整体结构,海上风力发电整体结构有直径更大的叶片、高度更高的塔筒,相应风机整体结构也有更大的柔性。
海上风机的基础结构更复杂,对地基承载力要求高。
海上风力发电结构除了风荷载的作用外,还常受到复杂的外载激励,如海冰荷载、海流荷载和波浪荷载等海洋环境荷载的影响,导致风机结构振动的外载激励增多,所以对三桩基础的海上风力发电整体结构的抗震性能提出了更高的要求。
1 模态分析该文研究的海上风力发电结构的基础形式为三桩门架式,属于三脚架式基础。
采用大型通用有限元软件ANSYS,建立了三桩基础海上风机整体结构——“基础—塔筒—机舱—轮毂—叶片”的有限元模型,考虑土与结构相互作用的影响,然后根据m 法对其进行了模态分析,得到了三桩基础海上风力发电整体结构的前十阶自振频率和固有振型,结果如下。
第一阶自振频率:0.277 Hz,固有振型:塔筒、机舱、轮毂和叶片沿y 方向摆动。
第二阶自振频率:0.284 Hz,固有振型:塔筒、机舱、轮毂和叶片振动沿x 方向振动。
第三阶自振频率:0.395 Hz,固有振型:三叶片各自沿y 方向振动,左边的叶片振动最大。
第四阶自振频率:0.443 Hz,固有振型:塔三叶片各自沿y 方向振动,右下的叶片振动最大。
第五阶自振频率:0.536 Hz,固有振型:塔筒沿y 方向振动较小,三叶片同步地沿着y 方向振动。
第六阶自振频率:0.737 Hz,固有振型:塔筒沿x 方向振动较小,三叶片绕着轮毂旋转振动。
第七阶自振频率:1.004 Hz,固有振型:三叶片各自绕轮毂旋转振动,左边叶片振动最大。
第八阶自振频率:1.042 Hz,固有振型:塔筒振动很小,左边的叶片振动较小,右边两叶片三桩基础海上风机整体结构的共振分析李 益 凡 威(中国电建集团河南省电力勘测设计院有限公司,河南 郑州 450007)摘 要:该文对三桩基础的海上风力发电整体结构进行了共振分析。
风力发电机组有限元分析(FEA)Finite Element Analysis (FEA) for wind turbines2007年7月4日4th July 2007风机的有限元分析(FEA)–越来越重要FEA for wind turbines -increasingly important仅仅是几年前,有限元分析还仅仅用于相对简单的模型的计算。
Only a few years ago, FEA was used for relatively simple models目前:Presently:•越来越多的风机部件使用有限元进行分析Increasing number wind turbine components analysed using FEA•越来越多的细节在有限元模型中的体现,例如Increasing amount of detail in FE model for instance°变接触角轴承的表示法variable-contact-angle bearing representation°部件间具有有限摩擦的接触面的表示法(法兰表面的张开和移位)representation of contact between parts with finite friction (opening/ sliding of flange surfaces)°已装配螺栓的详细表示法detailed representation of bolts in assemblies•越来越多的复杂的非线性有限元Increasing number of complex, non-linear FEAFE analysis for wind turbinesFE analysis for wind turbines风机的FEA –在设计中的位置FEA for wind turbines -place in the design processGH 有限元分析过程GH FEA process规范regulations制造限制manufacturing limits 客户经验/偏好client experience/ preference质量quality人员通过Personnel access安装/安全assembly/safety与其他部件的接口interfaces with other components材料限制material limits最优/可接受的设计Optimal /acceptable design设计性能指标Design performance indicators:成本(材料,产品,运输,操作)Cost (material, production,transport, handling)重量Weight 强度StrengthFE analysis for wind turbinesMain bearing assemblyNacelle structureRotor headNacelle extenderBlade sectionPitch bearingPitch bearing plates(on top and bottom of inner race)同一组载荷下不同设计概念的比较comparison of different concepts under identical load sets最终设计尺寸强度的详细分析detailed analysis of strength of final design geometry风机的FEA –在设计中的位置FEA for wind turbines -place in the design process风机的FEA –最新的分析技术FEA for wind turbines -state of art technology德国劳埃德船级社2003版规范GL (Germanisher Lloyd) 2003 regulations有很大一章是对有限元分析的原则性和细节的要求Extensive chapter with general and detailed FEA requirements挪威船级社(丹麦规范)DNV (Danish regulations)对有限元分析的要求很少,认证更具一般性Fewer FEA requirements, more generic in certificationGarrad Hassan的工作符合GL的要求或更进一步Garrad Hassan work to GL requirements or furtherGL 规范GL regulationsGH 分析发展GH analysis development 某个部件的有限元分析—模型和实际之间的小差距FEA for certain components –narrow gap between model and realityGL接受新的或者改进的有限元分析方法acceptance and adoption by GLof new or improved FEAFE analysis for wind turbines叶片Blade•叶片blade•叶根连接blade root connection机械部件Machinery components•轮毂hub•主轴mainshaft•轮毂主轴连接hub to mainshaft connection •锁定销locking pin•其他螺栓连接other bolted connections 机舱结构Nacelle structure•主机架/底盘main frame/ bed plate•发电机架generator frame•轴承座bearing housings塔架Tower•塔顶(包括偏航轴承的影响)tower top (including influence of yaw bearing)•通道口(门,电缆,航空灯)access openings (door, cables, aviation light)•屈曲分析(仅限于非标准门)buckling analysis (only: non-standard doors)风机的有限元分析-模型总结FEA for wind turbines -FE models summarisedFE analysis for wind turbines风机的有限元分析-模型总结FEA for wind turbines -FE models summarised叶片bladeFE analysis for wind turbinesFE analysis for wind turbines风机的有限元分析-模型总结FEA for wind turbines -FE models summarisedNLR线性非线性包括其它结构分析中的刚度影响螺栓的疲劳和极限失效,包括法兰表面的张开/滑移(所以螺栓弯曲)和变桨轴承的影响叶片叶根叶片NLR非线性关系non-linear relationshipSCFBladeblade include stiffness effects in analysis of other structureslinear blade rootfatigue and ultimate failure in bolts, including effects of opening/ sliding of flange surfaces (hence: bolt bending) and pitch bearingnon-linearNLR除了热点应力(计算应力集中系数),还需要对照设计极限应力检查绝对应力σudas well as hot spotstresses (to calculate SCF),absolute stresses are checked against the design ultimatestress, σud应力集中系数stress concentrationfactor对于非线性分析,应力范围I 近似关系适合于用于计算SCF 的载荷循环for non-linear analyses,approximate relationship for stress ranges appropriate for applied fatigue cyclesto calculate SCFFE analysis for wind turbinesyellow = normal contact (friction coefficient 0.2)red = bonded contact ÆMPC algorithm风机的有限元分析-模型总结FEA for wind turbines -FE models summarised机械部件machinery componentsFE analysis for wind turbines风机的有限元分析-模型总结FEA for wind turbines -FE models summarisedNLR非线性非线性非线性螺栓的疲劳和极限失效,包括法兰表面的张开/滑移(所以螺栓弯曲)和任意附近轴承的影响螺栓的疲劳和极限失效,包括法兰表面的张开/滑移(所以螺栓弯曲)和任意附近轴承的影响包括其它结构分析中的刚度影响和非线性接触轮毂对主轴的连接其它螺栓连接主轴成SCF 非线性非线性线性疲劳和极限破坏的应力热点,包括变桨轴承的非线性影响轴上的热点,凹槽或螺纹处的应力集中,包括轴承的非线性影响销和盘的热点轮毂主轴锁定销和锁定盘机械部件(g)pgMachinery componentshubstress hot spots for fatigue and ultimate failure, including non-linear effects of pitch bearings and bolted connectionsnon-linearSCF mainshafthot spots in shaft, stress concentrations at grooves or threads, including non-lineareffects of bearingsnon-linear SCF locking pin and disk hot spots in pin and disk linear SCF hub to mainshaft connection fatigue and ultimate failure in bolts, including effects of opening/ sliding of flange surfaces (hence: bolt bending) and any nearby bearings non-linear NLR other bolted connections fatigue and ultimate failure in bolts, including effects of opening/ sliding of flange surfaces (hence: bolt bending) and any nearby bearings non-linear NLR main bearinginclude stiffness effects and non-linear contact in analysis of other structuresnon-linear风机的有限元分析-模型总结FEA for wind turbines -FE models summarised机舱结构nacellestructureFE analysis for wind turbinesFE analysis for wind turbinesNLR/SCF SCF SCF非线性线性线性线性焊缝疲劳破坏和其它热点,包括法兰张开和偏航轴承的非线性接触焊缝的疲劳破坏模态分析包括其它结构分析中的刚度影响主机架/底盘副属机架/发电机架副属机架/发电机架轴承座机舱结构Nacelle structuremain frame/ bed platefatigue damage at welds and other hot spots, including effects of flange opening and non-linear contact in yaw bearing non-linearNLR/ SCFauxiliary frame/ generator frame fatigue damage at welds linear SCF auxiliary frame/ generator frame modal analysislinear SCF bearing housingsinclude stiffness effects of bearing housing in analysis of other structureslinearTFE analysis for wind turbines塔架tower8-point contactbearing (slewing ring)three bolted connections风机的有限元分析-模型总结FEA for wind turbines -FE models summarisedFE analysis for wind turbines风机的有限元分析-模型总结FEA for wind turbines -FE models summarisedSCF非线性非线性塔筒底部法兰焊缝的疲劳破坏,包括法兰张开的影响焊缝疲劳破坏,洞周围的应力集中,包括法兰张开的效果塔底法兰塔筒地基嵌入环线性压曲抗力系数线性特征值分析,分析方法不限,非均匀或非锥形壳屈居分析SCF非线性线性塔头法兰焊缝疲劳破坏,包括法兰张开和偏航轴承的非线性接触焊缝的疲劳破坏塔头通道口塔架Towertower top fatigue damage at top flange weld including effects of flange opening and non-linear contact in yaw bearing non-linearSCFaccess openings fatigue damage at weldslinear SCFbuckling analysis eigen-analysis of doors outside limits of analytical method, non-uniform or non-conical shellslinearlinear bucklingresistance factortower base flange fatigue damage at bottom flange weld including effects of flange openingnon-linear SCFtower foundation insert fatigue damage at welds, stress concentrations around holes including effects of flange openingnon-linearSCF有限元建模讨论FE modelling Workshop FE modelling讲到的有限元分析有:The following FEA will be presented:•塔筒底部的入口Access door in tower base•塔头法兰焊缝分析Tower top flange weld analysis•主轴分析Mainshaft analysis•轮毂和主轴之间的螺栓连接Bolted connection between hub and mainshaftFE analysis for wind turbinesFE analysis for wind turbines主要有限元目标Main FE modelling goals-重新获得焊趾位置的应力(疲劳分析)Retrieving stresses at weld toe locations (fatigue analysis)-重新获得门位置的von Mises 应力最大值(检查屈服)Retrieving the maximum von Mises stress at door location (check on yield)有限元模型细节FE model specifics-使用体单元而不是壳单元(GL 的要求)solid elements not shells (required by GL)-塔筒截段的高度最小为2.5 ×Dtower section height of minimum 2.5 ×D有限元建模塔筒入口FE modellingTower access door1/12重新获得焊趾位置的应力Retrieving stress results from the FE model at weldlocation推荐EuroCode3/ 国际焊接协会(IIW)Recommendations from EuroCode3/ International Institute of Welding (IIW)reference pointsstructural stresshot spot Fcomputed total stressF0.4*t0.4*t1*tt有限元建模塔筒入口FE modelling Tower access door2/12FE analysis for wind turbinesFE analysis for wind turbinesStep 1•在CAD (SolidWorks).软件中创建3维门模型3D door model created in CAD (SolidWorks).•在CAD 软件中创建面(用于切割焊缝位置的体)surfaces (to be used for cutting the solid at weld locations) created in CAD.•导入CAE 软件imported into CAE.切割CAD 模型中创建的面Cutting surfaces created in CAD package有限元建模塔筒入口FE modellingTower access door3/12Step 2•用导入的面切割导入的体imported solid cut using imported surfaces•创建新体new solids created•用从塔筒向外复制小面(操作后会删除)的方法切割法兰顶部区域(法兰侧焊接位置)flange top area (flange side weld locations) cut by copying small area offset from tower wall (cut areas are deleted after operation)用导入的面切割体Solid is cut by imported areas 顶部区域用复制的面切割(6个)Top area is cut bycopied areas (6 areas)有限元建模塔筒入口FE modelling Tower access door4/12FE analysis for wind turbinesFE analysis for wind turbiness w e e p m e s hStep 3用工作面切割体以便得到扫略网格The solid is also cut by the workplane to make a sweep mesh possible用水平工作面切割Cut by ‘horizontal’workplane有限元建模塔筒入口FE modellingTower access door5/12有限元建模塔筒入口FE modelling Tower access door6/12Step 41 门区域之外的主要塔筒截面是扫略得到的solid95高阶单元The main tower sections outside the door area are sweep meshed with solid95 higher order elements2 门区域截面用面单元(Mesh200)细化The door area cross section is finely meshed with surface (Mesh200) elements3 焊缝附近的位置面网格细化In the vicinity of the weld the surface mesh is refined网格细化Mesh refinementFE analysis for wind turbines有限元建模塔筒入口FE modelling Tower access door7/12Step 5门和周围直接连接的壳用扫略划分网格door and shell immediatelyaround it are sweep meshed门框架区域和塔架主体之间的壳用四面体单元直接划分shell between door frame region and main towermeshed automatically with tetrahedra15200个Solid95单元15200 Solid95 elementsFE analysis for wind turbinesFE analysis for wind turbinesX T Pointing South.Z T Vertically upwards.Y TPointing East.OriginAt each tower station.塔架载荷和挠曲的坐标系协调Co-ordinate systems for tower loads and deflections在CAD 例子中指定的轴和原点axes and origin designated in CAD,example shown here改变坐标系change system of coordinates有限元建模塔筒入口FE modellingTower access door8/12FE analysis for wind turbines对称条件Symmetry condition所有自由度约束Fixed in all DOFM y =1NmStep 6施加单位载荷Mx=1Nm 求解模型model solved with a unit load M x =1Nm应力Stress: vonMises变形比例Scale deformation =1e9有限元建模塔筒入口FE modellingTower access door9/12FE analysis for wind turbines不对称条件Asymmetry condition所有自由度约束Fixed in all DOFM z =1NmStep 7施加单位载荷M y = 1Nm 求解模型model solved with a unit load M y = 1Nm应力Stress: vonMises变形比例Scale deformation =1e9有限元建模塔筒入口FE modellingTower access door10/12FE analysis for wind turbines应力结果用沿焊缝的“移动的”随动局部坐标系Stress results are read in “moving”local coordinate system aligned with weldStep 8焊趾处的3个应力分量three stress components are read at weld toe:1.与焊趾平行的应力stress parallel with weld 2.与焊趾垂直的应力stress perpendicular with weld 3.面内剪应力in plane shear stress有限元建模塔筒入口FE modellingTower access door11/12FE analysis for wind turbinesS tre s s e s a t d o o r w e ld , o u ts id e to w e r w a llAA-6-4-2024681000.511.522.53T o w er h eig h t [m ]S t r e s s [P a ]S x_Mx S y_Mx S xy_Mx S x_My S y_My S xy_My P 1_Mx P 3_MxPosition Name Node X Pos Y Pos MaxPS_Mx Sx_Mx Sx_My Sy_Mx Sy_My [-][-][no][m][m][N/m 2][N/m 2][N/m 2][N/m 2][N/m 2]Door tower wall, outside A 137500.40 2.08-6.31-5.61-2.96-0.188.01Door tower wall, inside B 136260.37 2.19-9.22-6.08-12.04-2.93-13.92Door flange, outside C 199780.40 2.02-7.51-7.38-7.11-6.39-5.92Door flange, inside D 205820.37 2.13-8.09-6.40-8.93-5.08-7.86Base flange, outside E 623090.820.08-8.75-2.55 1.21-8.74 3.58Base flange, insideF622180.550.09-3.93-0.980.36-3.931.20Step 91计算M x 引起的主应力(PS )Principal Stresses (PS) due to M x load calculated 2 画应力结果图Graphs plotted of stress results3 Mx 引起的主应力(PS )的最大值决定热点max. (abs) PS due to the M x load determines hotspot position4 生成影响线矩阵Influence matrix generated有限元建模塔筒入口FE modellingTower access door12/12有限元建模塔顶法兰焊缝分析FE modelling Tower top flange weld analysis1/15Main FE modelling goal决定关联载荷和焊缝的应力(疲劳和极限分析)Determination of the relationship load and stress in weld (fatigue andextreme analysis)FE model specifics有限元模型细节-使用体单元Solid elements will be used-模型中包括部分塔架,法兰,偏航轴承Part of tower, flange, yaw bearing is modelled-非线性分析ÆNon linear analysisÆ-可能有(局部)连接的张开possible (partial) opening of connection-活动支座滚珠接触shifting bearing ball contact-变化的轴承滚珠刚度changing bearing ball stiffnessFE analysis for wind turbines有限元建模塔顶法兰焊缝分析FE modelling Tower top flange weld analysis2/15Step 1导入SolidWorks2维面模型SolidWorks2D model of areas is imported Step 2在轴承座圈位置,模型是由多线组成的面(90根线组成圆)At bearing race positions, areas formed from poly lines (90 lines forming circle) are modelled部分机舱底座part of nacelle bed plate偏航轴承yaw bearing刹车盘brake disk塔头法兰tower top flange塔筒tower sectionFE analysis for wind turbines有限元建模塔顶法兰焊缝分析FE modelling Tower top flange weld analysis3/15 Step 31 用圆来切割模型The circles are used to cut the model2 螺栓的位置需要切割(工作面切割)Cuts are made for bolt position (work plane cut)3法兰焊缝的位置需要切割(工作面切割)Cuts are made at flange weld location (work plane cut)螺栓位置bolt location多线组成的轴承座圈bearing race formed by poly lines螺栓位置bolt location焊缝位置收缩法兰weld locationnecked flangeFE analysis for wind turbinesFE analysis for wind turbinesStep 41 用表面单元划分面的网格areas meshed with surface elements 2 焊缝位置网格细化(关心的区域)mesh refined in weld area (area of interest)网格细化M e sh r e fin e m e nt有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis4/15有限元建模塔顶法兰焊缝分析FE modelling Tower top flange weld analysis5/15Step 5网格沿轴向旋转mesh rotated around axis28000个Solid45 单元28000 Solid45 elementsFE analysis for wind turbinesFE analysis for wind turbinesd ir e cti o n o f i n i t i a l c on t a c td ir e c t i o n o f i n i t i a lc o n t a c t 使用Li n k 10单元的初始应变特性表示间隙R e p r e s e n t a t i o n o f g a p m o d e l l e d u s i n g i n i t i a l s t r ai n p r o p e r t y o fL I N K 10 e le m e nt sStep 7用Link10单元(单元仅受压)表示轴承滚珠Link10 (compression only elements) are modelled to represent the bearing balls有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis6/15FE analysis for wind turbinesStep 8用Beam 4单元表示螺栓Beam 4 elements are modelled to represent the boltsStep 91 模型中包括接触Contact is modelled2 模型中包括“载荷伞”Load umbrella is modelled塔头=载荷施加Tower top =load applicationcontact有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis7/15FE analysis for wind turbinesStep 10载荷步1:通过给螺栓设置温度(收缩梁)实现螺栓预加载Load step 1: Bolt preload step is solved by applying temperature on bolts (contracting beams)结果显示:仅螺栓预载荷(M z =0kNm)Results shown: Only bolt pre load (M z =0kNm)变形比例= 500Scale deformation = 500有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis8/15Step 11载荷步3:分5个子步施加实际Mz弯矩(7000kNm)求解模型Load step 3: Model is solved with realistic Mz moment (7000kNm) in 5 sub steps结果显示:Mz =7000kNmResults shown: Mz =7000kNm变形比例= 100Scale deformation = 100失去接触Loss of contact有限元建模塔顶法兰焊缝分析FE modelling Tower top flange weld analysis9/15FE analysis for wind turbinesFE analysis for wind turbines受压侧Compression side M z =7000kNm Sy stress 受拉侧Tensile side M z =7000kNm Sy stressLoc1Loc2焊趾位置应力最大weld toe location highest stressed有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis10/15FE analysis for wind turbinesNLR between stress and moment at 4 tower top locations-150.00-100.00-50.000.0050.00100.00150.00200.00-8000000-6000000-4000000-20000002000000400000060000008000000moment [Nm]S t r e s s [M P a ]Outside 1Inside 1Outside 2Inside 2-螺栓应力和弯矩的非线性关系Non linear relationship bolt stress and applied moment -M y ,hub 或M z ,hub 载荷时程(旋转坐标系)Load time history M y ,hub or M z ,hub (rotating coordinate sys.)Step 11生成应力/ Nm 焊趾非线性影响关系non linear influence relationship weld toe stress / Nm is generated螺栓应力时程(使用Bladed )bolt stress time history (using Bladed)极限螺栓应力extreme bolt stress寿命周期内疲劳损伤lifetime fatigue damage有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis11/15FE analysis for wind turbinesR e s u lt a n t c o n t a c t a n g le-80.00-60.00-40.00-20.000.0020.0040.0060.0080.0020406080100120140160180200R a d ia l p o s it io n [d e g ]C o n t a c t a n g l e [d e g ]T e n s i le s i d eC o m p r e s s i o n s i d e对于每一个载荷步,可以计算得到合成的接触角For each load step, a resultant contact angle can be calculated接触结果与前面使用固定45度接触角建模时得到的结果十分不同Æ对焊缝位置应力结果影响很大(大概减少20%)Contact results are very different compared to previous modelling with a fixed contact angle of 45 degrees Ælarge influence on stress results weld location (around 20%reduction)Loc1, 受拉侧tensile side Loc2, 受压侧compression side有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis12/15FE analysis for wind turbines检查轴承滚珠接触Check on realism of bearing ball contact轴承面外所有自由度完全约束Outside bearing areais fully constraint in all DOF刚度很高的载荷伞Very stiff load umbrellaF radial = 8000 kN对称条件Symmetry conditions有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis13/15FE analysis for wind turbines无量纲力无量纲位移有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis14/15FE analysis for wind turbines8点接触轴承(回转支承)8-point contactbearing (slewing ring)3个螺栓连接three bolted connections有限元建模塔顶法兰焊缝分析FE modellingTower top flange weld analysis15/15有限元建模主轴FE modelling Mainshaft 1/11主要有限元建模目的Main FE modelling goal指定主轴表面“热点”的位置Allocate hotspot locations on shaft surface产生“热点”周围的影响线矩阵Produce influence matrix for all the foundhotspots有限元模型细节FE model specifics-静态分析(“热点”分析)Static analysis (hotspot analysis)-使用体单元Solid elements will be applied-模型中包括部分轮毂Part of hub is modelled-模型中包括轴承状态(双球滚子形式)Bearing behaviour (double spherical rollerlay-out) is modelledFE analysis for wind turbinesStep 1主轴的SolidWorks3维模型,轮毂和内轴承座圈很重要SolidWorks3D model of mainshaft, hub and inner bearing races is imported Step 2体模型被切开,沿主轴轴向的截面和部分轮毂被留下Solid model is cut and cross section over the length of the shaft and part of hub is left有限元建模主轴FE modelling Mainshaft2/11FE analysis for wind turbinesFE analysis for wind turbinesStep 3Mesh200单元划分面的网格area meshed with Mesh200 elementsStep 4在关键区域首次细化网格initial mesh refined in critical areasStep 5在关键区域二次细化网格(轮毂侧和齿轮箱侧轴承位置)second mesh refined in critical areas (hub side and gearbox side bearings positions)网格细化1Mesh refinement 1网格细化2Mesh refinement 2有限元建模主轴FE modellingMainshaft3/11有限元建模主轴FE modelling Mainshaft4/11 Step 6面网格绕主轴轴线旋转生成solid 95(高阶)体单元模型area mesh revolved around shaft axis to generate asolid model with solid 95 (higher order) elements36200个Solid95单元36200 Solid95 elementsFE analysis for wind turbines有限元建模主轴FE modelling Mainshaft5/11 Step 6用非常细(没有刚度)的shell 93单元划分体模型网格The solid model is meshed with very thin (nostiffness) shell 93 elements.这些单元会在后处理中用于获取模型表面的应力These elements will be used in the postprocessing phase to extract stresses from thesurface of the model6200个shell 93 单元6200 Shell93 elementsFE analysis for wind turbinesFE analysis for wind turbinesStep 7用非线性Link10单元的“蛛网”表示轴承滚子(刚度建模)Spider webs of non linear Link10 elements are modelled to represent the bearing rollers (stiffness is modelled)Link10单元仅用于承受拉伸载荷的建模The Link10 elements are modelled only to take loading under tensile432个非线性Link10单元(仅受拉伸)432 non linear Link10 elements (tensile only option)轴承中心位置的中心节点Central node at bearing centre position把Link 单元连接到内轴承座圈的外表面Link elements are connected to the outside surface of the inner bearing race有限元建模主轴FE modellingMainshaft6/11FE analysis for wind turbinesStep 81 载荷施加点和“载荷伞”建模1 Load application node and load umbrella is modelled2 齿轮箱安装刚度建模(具有恰当刚度的梁单元)2 Gearbox mounting stiffness is modelled (Beam elements with correct stiffness)3 施加边界条件3 Boundary conditions are applied边界条件Boundary conditions齿轮箱安装刚度Gearbox mounting stiffness施加齿轮箱重量Gearbox weight application施加载荷Load application边界条件Boundary conditions有限元建模主轴FE modellingMainshaft7/11FE analysis for wind turbinesStep 9模型上施加了下列载荷:The following loads have been applied to the model 11个单位载荷:11 unit loads:+F x , -F x , +F y , -F y , +F z , -F z +M x , +M y , -M y , +M z , -M z4个变速箱载荷(实际重力载荷):4 Gearbox loads (real gravity loads):+F y , -F y , +F z , -F zShell 93 单元M y = 1NmvonMises 应力Solid 95 单元M y = 1NmvonMises 应力Shell 93 elements M y = 1NmStress=vonMisesSolid 95 elements M y = 1NmStress=vonMises有限元建模主轴FE modellingMainshaft8/11。
第46卷第4期2017年8月船海工程SHIP&OCEAN ENGINEERINGVol.46 No.4Aug.2017D01:10. 3963/j. issn. 1671-7953. 2017.04.037新型深海风机运输安装一体化半潜式平台整体强度分析和校核于雯S李曙生S袁宇波2(1.泰州职业技术学院,江苏泰州225300;2.南通中远船务有限公司,江苏南通226000)摘要:为了满足海上风力开发的深海化、专业化、大型化的需求,解决因海上气候恶劣,有效工作时间短 而导致的工作效率低的问题,设计一种满足风车整体运输、安装一体化的新型半潜式风车安装平台。
根据风 车整体运输、安装一体化的设计指标,提出新型半潜式风力发电安装平台船型方案设计,确定平台的总布置、主尺度和结构形式;采用s esam软件的设计波进行平台整体强度分析,在满足结构强度以及保证平台性能的 基础上进行方案优化,设计出满足结构强度以及性能要求的最优船型。
关键词:风车安装;半潜式平台;整体强度评估;运输安装一体化中图分类号:U656.6 文献标志码:A风电机组的安装是海上风电场建设的关键技 术之一。
现有的近海风力发电安装作业船,都是 针对近海风机安装,不能满足风机大型化、深海化 的发展趋势[13]。
深海风电安装平台要适合在海 况恶劣的深海工作,具备自航、运输、起重、避风等 复合功能。
为满足以上要求,结合半潜式平台优 良的耐波性能[45],设计平台采用半潜式结构,采 用流线型浮筒减小航行阻力,动力方面采用喷水 推进方式,更能适应全球未来海上风电大型化和 深海化的发展趋势。
采用有限元方法对整船结构 强度进行分析和研究[68]。
在提出一种新型半潜 式风力发电安装平台船型方案设计的基础上,根 据平台的主尺度、总布置特点,采用D N V S E S A M 系列软件根据其受力情况建立平台有限元模型,根据计算仿真结果科学地设计船体结构,对高应力区域做适当加强,进而对船体结构进行优化。
离岸深水全直桩码头的地震损伤分析
地震是自然灾害中常见的一种,对码头等基础设施造成的损害往往是不可忽视的。
离
岸深水全直桩码头是一种常见的码头结构,其地震损伤分析对于码头的设计和加固具有重
要意义。
1. 码头结构的受力特点:离岸深水全直桩码头主要由桩、连梁和甲板组成,其整体
刚度较大,具有相对较好的抗震性能。
在地震作用下,桩、连梁和甲板之间的连接处会受
到较大的力矩和剪力作用,造成结构的损伤。
2. 地震作用的参数分析:地震力是地震作用下码头结构的重要载荷,其大小和分布
形式对结构的损伤具有重要影响。
通过地震工程研究和地震观测数据分析,可以确定地震
的烈度、震源距离和地震频率等参数,进而计算出地震力的大小和分布情况。
3. 码头结构的地震响应分析:地震作用下,离岸深水全直桩码头会产生振动和变形。
通过有限元分析和结构动力学方法,可以计算出码头结构在地震作用下的位移、加速度和
内力等响应,进而评估结构的抗震性能和确定可能的损伤位置。
4. 码头结构的地震加固措施:根据地震损伤分析的结果,可以确定码头结构的薄弱
部位和损伤可能性较大的位置。
对于这些位置,可以采取合理的加固措施,如增加桩的直
径和数量、加固连梁和甲板的连接等,提高码头结构的抗震性能。
大跨度桥梁抗震分析中的整体有限元法及其应用目录一、内容概要 (2)1. 桥梁工程的重要性 (2)2. 抗震分析的意义与挑战 (3)二、有限元法概述及其在桥梁抗震分析中的应用 (4)1. 有限元法基本概念与原理 (6)1.1 有限元法定义与发展历程 (7)1.2 基本原理与计算步骤 (8)2. 有限元法在桥梁抗震分析中的应用现状 (9)2.1 应用范围及优势 (10)2.2 存在的问题与挑战 (11)三、大跨度桥梁整体有限元建模与分析方法 (13)1. 整体有限元建模流程 (14)1.1 模型建立前的准备工作 (15)1.2 模型建立过程及参数设置 (16)1.3 模型验证与校准 (17)2. 大跨度桥梁整体分析方法 (19)2.1 静力分析方法 (21)2.2 动力分析方法 (22)2.3 抗震性能评估指标 (23)四、大跨度桥梁抗震分析中的关键技术与策略 (25)1. 地震波输入与选择 (27)1.1 地震波特性分析 (28)1.2 地震波输入方法比较与选择 (29)2. 结构损伤评估与修复策略 (30)2.1 结构损伤识别技术 (32)2.2 损伤程度评估方法 (34)2.3 修复策略与建议 (35)一、内容概要本文档主要介绍了大跨度桥梁抗震分析中的整体有限元法及其应用。
整体有限元法是一种将结构划分为多个单元,通过离散化的方法对整个结构进行建模和求解的方法。
在大跨度桥梁抗震分析中,整体有限元法具有较高的计算精度和效率,能够有效地模拟桥梁在地震作用下的响应过程,为桥梁的抗震设计提供有力的支持。
本文档首先介绍了大跨度桥梁的基本结构特点和抗震要求,然后详细阐述了整体有限元法的基本原理、方法和步骤,包括单元划分、刚度矩阵和边界条件设置等。
通过实例分析,展示了如何运用整体有限元法对大跨度桥梁进行抗震分析,以及如何根据分析结果优化结构设计,提高桥梁的抗震性能。
对整体有限元法在大跨度桥梁抗震分析中的应用前景和技术发展趋势进行了展望。
第50卷第1期2021年2月船海工程SHIP&OCEAN ENGINEERINGVol.50No.1Feb.2021DOI:10.3963/j.issn.1671-7953-2021-01.024海上风电大直径单桩基频率因素的数分析,马,,吴,陈帮(明阳智慧能源集团股份公司,广东中山528400)摘要:针对海上风电整机结构其低阶固有频率的控制问题,基于有限元分析技术开展小直径桩土相互作用的数值模拟分析,与试验结果、API P算结果进行对比,验证了数值模型的适用性,分析大直径单桩结构桩内土芯、冲度、表层施对整频率的,结表明,基于元实体单元计算所得一、整频率要明显大于传统方法计算结果;桩内土芯对一、整频率小;在一定范围内,随着冲度的,一、频率明显,认为传统方法对于大直径单桩的适用;在实际工程中,桩内土芯的忽略;对于大直径单桩,在一围内可以忽略冲刷的影响。
关键词:大直径单桩;固有频率;土芯;冲刷深度中图分类号:U661.43文献标志码:A文章编号:1671最953(2021)01-0099-05海上风机基础的制作和安装费用约占风电机组工程总成本的20%〜30%(门$,内海上风电场的风电基础结构类型主要力式基础、单桩基础、基、负压桶式基础、基等结构(2)$中,大直径单桩具有结构型简单,方法成熟以及施工维护较为方便等优点,是海上风电场近浅海水域开发过程中应用多的基础型式(3)$由于风、波浪、海等外部环境荷载作用,使风电结构相当大的水平力和弯矩荷载。
整的频率开叶轮转动频率区间,以证整体生,满足 安全使用的要求。
为此,重点分析大直径桩的桩相互作用以及桩内土芯、冲度等因素对整频率的$1桩土作用数值模型验证1.1有限元模型验证的模型是出自Maw lock(4)在Sabine地区的试桩,荷载作用位置距0.3m,桩的度为12.81m$三轴试验中最大主应力差一半时的应收稿日期:2020-07-14修回日期:2020-08-21第一作者:孙肖菲&1991-),男,硕士,助理工程师研究方向:海上基础结构变值.5。
第46卷第9期 2 0 1 5年5月 人 民 长 江
Yangtze River Vo1.46,No.9
May. 2Ol5
文章编号:1001—4179(2015)09—0009一O4
大水深离岸型海上测风塔整体有限元分析
周 颖,吴海明,齐志诚 (中国电建集团北京勘测设计研究院有限公司,北京100024) 摘要:我国目前没有可直接运用于大水深、工程环境复杂的海上测风塔的设计技术标准。根据我国《港口工 程地基规范》及相关工程技术标准,利用通用有限元结构软件分析了测风塔结构受力;同时依据国际通用的 美国石油行业《海上固定平台入级与建造规范》(即API RP 2A规范),利用海工结构有限元分析软件对测风 塔进行了结构分析。计算分析内容包括:测风塔多种设计荷载、多种设计工况、关键节点位移、杆件内力和疲 劳分析等。两种计算结果基本一致,均能满足相应规范要求,说明所选取的设计参数较为合理,相关的工程设 计经验可为类似工程及海上风电场建设提供参考借鉴。 关键词:测风塔;风场;疲劳分析;海上风电 中图法分类号:TM614 文献标志码:A DOI:10.16232/j.enki.1001—4179.2015.09.003
相较于陆上风力发电,海上风力发电具有年发电 量高,占用土地资源少等优点,是战略性新兴产业的重 要内容,也是风电产业发展的技术制高点。建设海上 风电场,首先要获得不同高度处的准确风资源参数,需 要在海上建造测风塔。 海洋环境的复杂性和测风设备的运行要求,使测 风塔的设计比常规的陆地工程复杂很多。目前,我国 已建成的海上测风塔很少,且主要位于较浅水域。由 于我国尚无这方面专门的技术标准,在设计中多参照 《港口工程地基规范》 及陆上相关技术标准进行。 某海上测风塔工程距离海岸约27 km,水域深度 约22.5 m,是大水深离岸型海洋工程项目。为确保测 风塔安全可靠、经济适用,设计者参照《港口工程地基 规范》和美国石油行业规范《海上固定平台入级与建 造规范》 (即API RP 2A规范)及相关工程技术标 准,利用通用有限元结构分析软件和海工结构有限元 分析软件分别建立整体结构模型,从测风塔多种设计 荷载、多种设计工况、关键节点位移和杆件内力等几个 方面进行计算分析和研究,所取得的工程设计经验可 为类似工程及海上风电场建设提供参考借鉴。 1工程概况 1.1工程地质 拟建测风塔塔址处60.00 m深度范围内的地层按 其成因可分为9层:①层为粉、细砂,分布于海域表层, 厚4.80 m;②层为粉质黏土,厚约4.50 m;③层为黏 土,厚约3.10 m;④层为粉质黏土,厚约4.60 m;⑤层 为粉质黏土,厚约7.70 m;⑥层为粉砂,厚约6.90 m; ⑦层为粉质黏土、黏土互层,厚约13.60 m;⑧层为粉 质黏土,厚约7.00 m;⑨层为粉质黏土,厚约7.80 m。 1.2 水文条件 工点海区内潮流一般为不规则半日潮流,以往复 流为主,涨落潮流方向基本与海流一致,实测潮流流速 范围为0~0.79 m/s。余流方向以西南流(SW)为主, 最大流速为0.5 m/s,表层余流大于底层余流。波浪 玫瑰花图见图1,波浪要素统计见表1。 1.3测风塔结构形式 测风塔上部结构采用自立式,塔高100 m,塔架截 面采用四边形桁架结构形式。测风塔下部基础采用钢
收稿日期:2014—10—10 作者简介:周 颖,女,工程师,主要从事海工结构设计工作。E—mail:zhouying@bhidi.corn 10 人 民 长 江 筋混凝土高桩承台,承台为8.0 m x8.0 m见方,厚1.7 m的C40混凝土结构,4根桩径为I.4 m的钢管桩作 为基桩,桩长为70 In,斜桩6:1。
频率 0 2% 4% .
 ̄9
1 ̄l_40.
28:
昌 . ̄63 ̄_1.59am 亡
>2.0m
图1波浪玫瑰图
表l测风塔处波浪要素统计(最不利方向sE)
(6)海生物。海生物密度取值1.3 t/m 。 (7)地震力。抗震设防烈度为Ⅶ度。 (8)风速仪、风向标等自重取0.5 kN,集中力作 用于外伸横梁末端。 (9)飞溅区范围和腐蚀余量。飞溅区范围为高程 一1.6—3.8 m,腐蚀余量取0.5 mm/a。
3控制标准 3.1变形控制标准 由于没有测风塔相关的基础设计专业标准,参考 《建筑地基基础设计规范》(GB50007—2002)和《高耸 结构设计规范》(GBJ50135—2006)的有关规定确定变 形控制标准 。测风塔结构承受风荷载为动荷载, 且承受水平荷载较大,考虑到海上测风塔水平位移对 其正常测风功能影响有限,基础的最大倾斜可按0.5% 控制。高耸结构水平位移限值为H/75,即1 300 mm。 3.2结构应力控制标准 该工程设计年限为20 a,桩基础钢筋混凝土承台 应力应满足《混凝土结构设计规范》(GB50010— 2002)等相关规范要求。桩基础钢管桩采用Q345B型 钢材,结构应力应满足《钢结构设计规范》(GB50017) 等相关规范要求。 2设计荷载 4结构受力分析
测风塔设计考虑的荷载主要包括自重、风荷载、波 浪力、水流力、冰压力、船舶或漂浮物撞击力及地震力 等。 (1)风荷载。10 m高度处的平均风速计算结果见 表2。 (2)波浪力。根据《海港水文规范》 ,按波浪要 素计算波流对桩基的作用,波浪力计算的水位分别采 用极端高水位3.62 nl和极端低水位一1.53 In,波高取 设计波高5.81 m。
表2平均风速计算结果 重现期/a 至旦堕 堕 竺:!::1 1 rain 10 rain
(3)水流力。根据《港口工程荷载规范》计算 , 水流力计算的设计流速为0.84 m/s。 (4)船舶(或漂浮物)撞击力按200 kN计算。 (5)冰压力。经调研,本海域等深线10 in外没有 冰,只有浮冰,冰温0oC,极限冰压750 kPa,冰温系数 取1.0。
结构受力分析首先采用通用有限元软件进行。计 算的荷载工况组合按照承载能力极限状态组合和正常 使用极限状态标准组合确定,然后,参照API RP 2A规 范相关内容,按对结构产生最不利影响的相互荷载条 件重新定义了荷载组合,并利用海工结构有限元分析 软件进行计算。最后,将两次计算结果进行比较,以确 定合理可行的设计方案。 4.1结构计算假定 计算的关键在于水平荷载作用下桩土计算模型的 选取及结构整体的空间分析。计算时采用假想嵌固点 法假定,桩入土一定深度处视为假想嵌固点,结构分析 时不再考虑土的影响。刚性固定端位于设计泥面垂直 向下 处。设计泥面的位置在自然泥面下的距离要按 照地质条件来决定。 值按照下列经验公式来确定:对 淤泥T=(7~8.5)D,对硬黏土T=(3.5~4.5)D, 当缺乏土壤资料时T:6D,其中,D为桩腿的外径。 在结构的总体静力分析中,桩基础用等效桩模拟,等效 桩的长度取5倍桩径。 4.2通用有限元分析 采用通用有限元结构分析与设计软件进行结构整 第9期 周 颖,等:大水深离岸型海上测风塔整体有限元分析 体计算。上部结构塔架用空间梁单元模拟,承台模拟 成块单元,节点按刚性节点处理。整体计算模型见图 2。
图2测风塔整体计算模型 荷载工况组合主要采用了以下组合。
4.2.1 承载能力极限状态持久组合 (1)1.2G+1.4WF(NOICE)+1.5×0.7(,+ WAVE+BOAT+FLOATINGICE) (2)1.2G+1.4WF(ICE)+1.5×0.7(,+WAVE +BOAT+FLOATINGICE) (3)1.2G+1.4W45(NOICE)+1.5×0.7(,+ WAVE+BOAT+FLOATINGICE) (4)1.2G+1.4W45(ICE)+1.5×0.7(,+WAVE +BOAT+FLoATINGICE、 4.2.2 正常使用极限状态标准组合 (1)G+ F(NOICE) (2)G+ F(ICE) (3)G+W45(NOICE) (4)G+W45(ICE) (5)G+ F(NOICE)+,+WAVE+BOAT+ FLoATINGICE (6)G+WF(ICE)+,+ +BOAT+ FLoATINGICE (7)G+W45(NOICE)+,+WAVE+BOAT+ FLoATlNGlCE (8)G+W45(ICE)+,+WAVE+BOAT+ FLoATlNGlCE 其中:G为重力,,为裹冰的竖向荷载, F为正 面风荷载,W45为45。风向风荷载,ICE/NOICE为裹 冰/无冰状态,WAVE为波流力,BOAT为船舶撞击荷 载,FLOATINGICE为浮冰撞击荷载;字母前面的系数 为各荷载的分项系数。 塔顶代表节点最大位移发生在工况组合(7)时,
为571.998 mm。桩顶代表节点最大位移发生在工况 组合(1)时,为146.14 mm。海底泥面代表节点最大 位移发生在工况组合(1)时,为19.651 mm。 由于代表单元应力数据较多,此文不再赘列,其中 混凝土单元应力计算结果满足《混凝土结构设计规 范》(GB50010—2002)等相关规范要求,上部塔架及钢 管桩应力均满足《钢结构设计规范》(GB50017)等相 关规范要求。 4.3海工结构有限元分析 海工结构有限元分析软件精于海上平台失效的非 线性结构分析以及风、海流、海浪和地震引起的动态响 应分析。在利用该软件进行结构整体分析计算时,塔 架用空间梁单元模拟,承台按板单元模拟,节点按刚性 节点处理。主要分析在各种极端环境工况下,上部塔 桅和下部腿桩基础的节点位移、杆件内力(应力)、节 点冲剪等内容 。 对于极端风暴条件、操作条件,风力的作用方向与 波浪、流的作用方向一致。本工程极端风暴条件波流 选择0。,45。,90。,135。,180。,225。,270。,315。八个方 向,其方向定义见图3。
/ 。 口一
+ I 90*
图3风力作用方向 根据API RP 2A 21 st的规定,平台应按对结构产 生最不利影响的相互荷载条件进行设计。荷载条件应 包括环境荷载中相应的固定荷载和活荷载。本次设计 中极端风暴荷载条件下的荷载组合系数见表3。表 中,APP表示附属结构的重量(塔桅爬梯重量);LIVE 表示活荷载,包括设备中储备的液体、消耗品、人员活 动等,设备未占区域按2.5 kN/m 考虑;BOAT表示靠 船力;DEDS表示平台结构自重;工况SW00~SW31 是极端风速荷载,工况SC00~SC31是极端环境荷载 条件下波流作用。波流荷载根据API RP 2A 21 st计 算。 桩顶及桩腿泥面高程处3个结构关键节点最大位 移见表4,代表梁最大内力见表5,构件长细比较核见 表6。 各桩顶节点、梁单元分组冲剪应力UC值均小于