[4,5] 由于硬件计算能力和软件计算模型的限制, 目前对结构的流固耦合数值分析大部分还是二维的研究 [4] 计算机和软件技术的发展使得全三维的流固耦合数值模拟成为可能, 本文基于通用有限元软件 ANSYS, 对 CFRP 索的风致振动进行了三维流固耦合数值仿真, 对计算中的相关参数进行了细致讨论,

Size: px

Start display at page:

Download "[4,5] 由于硬件计算能力和软件计算模型的限制, 目前对结构的流固耦合数值分析大部分还是二维的研究 [4] 计算机和软件技术的发展使得全三维的流固耦合数值模拟成为可能, 本文基于通用有限元软件 ANSYS, 对 CFRP 索的风致振动进行了三维流固耦合数值仿真, 对计算中的相关参数进行了细致讨论,"

辣焚厍
5 years ago
Views:

1 CFRP 索流固耦合风振研究陆新征叶列平马千里江见鲸北京清华大学土木工程系, 摘要 : 碳纤维增强复合材料 (CFRP) 由于轻质高强, 在大跨领域有着很好的应用前景 CFRP 索材由于重量很轻, 其风振行为需要给予更多关注本文利用通用有限元软件 ANSYS 提供的流固耦合分析功能, 对 CFRP 索的风振行为进行了研究, 讨论了不同的湍流模型, 以及 CFRP 索在均匀风场和脉动风场下的风振行为分析结果表明,CFRP 索由于质量很轻, 气弹阻尼效果明显, 如果设计合理, 其风振作用可能会优于同尺寸钢索关键词 :CFRP 索, 大跨, 风振, 流固耦合中图分类号 :TU328 文献表示码 :A Study on the Fluid-Solid Coupled Wind-induced Vibration of CFRP Cable Lu Xin-zheng, Ye Lie-ping, Ma Qian-li, Jiang Jianjing Department of Civil Engineering, Tsinghua University, Beijing, Abstract: Due to the lightweight and high strength of Carbon Fiber Reinforced Polymer (CFRP), it has good application future in large span structures. Since the self-weight of CFRP cable is very light, more study should be given on its wind-induced vibration. In this paper, the wind-induced vibration of CFRP cable is studied with the fluid-solid coupled simulation function of ANSYS, a common finite element package. Different turbulence models are discussed. And the responses of CFRP cable in stable and fluctuant wind fields are simulated. The numerical results show that due to the lightweight of CFRP cable, the aerodynamic damping effect is more obvious. So with proper design, the wind-induced vibration of CFRP cables may have better performance than steel cables with the same size. Keyword: CFRP cable, large span, wind-induced vibration, fluid-solid coupled 1 引言碳纤维增强复合材料 (Carbon Fiber Reinforced Polymer, 简称 CFRP) 由于具有轻质高强耐腐蚀等诸多优点, 因而在大跨结构领域有着广阔的应用前景 [1-3] 而用 CFRP 索代替钢索用于索支撑结构, 是目前 CFRP 在大跨领域应用中一个比较可行的方案然而, 由于索材质量很轻, 侧向刚度较小, 因此由风引起的风振效应不可忽略国内外已经发生过多起钢索由于风振而导致严重损伤甚至破坏的事故由于 CFRP 比钢索的质量还要轻, 可能应用的跨度会更大, 因此其风致振动问题需要进行专门研究尽管目前基于风洞实验及随机振动理论, 结构风工程已经得到了长足的发展但是, 由于建筑结构体形很大, 风洞实验尺寸效应问题突出, 且如果牵涉到结构和空气流场的耦合振动问题, 则实验模型相似比设计将变得更为复杂, 且风洞实验的代价也较高因此, 基于现代计算流体动力学 (Computational Fluid Dynamics, 简称 CFD ) 和计算固体动力学 (Computational Structural Dynamics, 简称 CSD) 的流固耦合数值仿真有着很好的应用前景基金项目 : 国家自然科学香港澳门青年学者合作研究基金 ( 项目编号 : ) 与国家自然科 1 学重点基金 ( 项目编号 : )

2 [4,5] 由于硬件计算能力和软件计算模型的限制, 目前对结构的流固耦合数值分析大部分还是二维的研究 [4] 计算机和软件技术的发展使得全三维的流固耦合数值模拟成为可能, 本文基于通用有限元软件 ANSYS, 对 CFRP 索的风致振动进行了三维流固耦合数值仿真, 对计算中的相关参数进行了细致讨论, 并和钢索的风致振动进行了对比 2 ANSYS 的流固耦合计算功能 2.1 计算方案一般流固耦合问题中, 包含了三个不同的分析模块 [4] : (1) 流体域 (CFD): 由于风速一般比较低, 可以将空气看作不可压缩流体 ; (2) 结构域 (CSD): 由于结构一般不允许在风荷载下屈服, 因此结构一般为考虑大位移带来的几何非线性的固体 ; (3) 流体 - 固体界面 : 由于流体计算一般基于欧拉网格描述, 而固体计算一般基于拉格朗日网格描述, 因而流体和固体之间的界面是流固耦合分析的一个关键问题目前比较成熟的分析方法是任意拉格朗日网格法 (Arbitrary Lagrangian-Eulerian, 简称 ALE) 计算的流程一般为 : (1) 以未变形的结构作为流体的计算边界, 计算流体流场, 得到作用在结构表面的压力 ; (2) 将该压力作用于结构, 使结构产生变形, 进而根据结构的变形修改流体边界及流体网格, 重新进行流体计算分析 ; (3) 重复迭代步骤 (1) 和 (2) 直到流体和固体界面两次计算的压力差小于误差容限由于结构域 (CSD) 计算目前已经非常成熟, 且土木结构即使是大位移往往变形也有限, 在有合适的网格划分技术的情况下, 界面网格调整相对也并不复杂, 因此流固耦合计算的精度和难点主要表现在对于土木结构这样的钝体外形, 如何能够准确计算其周边流场及压力场 2.2 流体分析模型在 ANSYS 中, 提供的流体分析模型为雷诺平均法 [5,6] 其核心是不直接求解瞬时 Navier-Stokes 方程, 而是求解平均化的雷诺方程这样可以有效减少计算量, 因而得到了大量的工程应用, 也相对比较成熟雷诺平均法由于需要对瞬时湍流进行平均化处理, 因此如何构建合适的平均化的湍流模型是该方法的关键问题,ANSYS 程序中提供了 7 种不同的湍流模型, 分别是 : (1) Standard k-ε model ( 本文简称 SKE 模型 ) (2) Re-Normalized Group Turbulence model( 简称 RNG 模型 ) (3) k-ε Model due to Shih ( 简称 NKE 模型 ) (4) Nonlinear Model of Girimaji ( 简称 GIR 模型 ) (5) Shih, Zhu, Lumley Model ( 简称 SZL 模型 ) (6) k-ω Turbulence Model( 本文简称 SKW 模型 ) (7) Shear Stress Transport Turbulence Model ( 本文简称 SST 模型 ) 已有研究表明, 对于同样的问题, 采用不同的湍流模型可能得到差异悬殊的结果 [5] 尤其是对于土木结构这样的钝体, 差距可能更大 Standard k-ε 模型已经被大量研究证明在分析钝体流场时是不合适的 [4,5,7], 而近年来计算风工程中采用得比较多的是 RNG 模型 [8] 因此, 为了对上述模型有一个全面的认识, 本文首先对不同湍流模型的计算结果进行了对比分析对比计算包含以下三个不同的基本分析对象 : 2

3 (a) 正方形房屋的风压分布, 如图 1(a) 所示 ; (b) 正方形柱体绕流流速分布, 如图 1(b) 所示 ; (c) 圆柱体绕流流场分析, 如图 1(c) 所示 ; H α H (a) 正方形房屋风压分布 (b) 正方形柱体绕流 (c) 圆柱体绕流 (a) Wind Pressure on the square building (b) Wind pass a square cylinder (c) Wind pass a round cylinder 图 1 流场基本分析对象 Fig. 1 Fundamental objects in fluid analysis 2.3 基本模型结果讨论各工况分析结果简要介绍如下 : (a) 对于正方形房屋的风压分布, 已有研究表明 [5,7,9],Standard k-ε 模型会过高估计迎风面檐口的风压, 其他模型均在此基础上有所改进, 不同计算模型得到的风压分布如图 2 所示 ( 实验资料来源 : 文献 [7]) 可见, 数值计算结果和风洞实验结果总体上吻合较好, 但是在迎风面檐口差异较大, 相对而言,RNG GIR SZL SKW 模型结果较好,SKE 模型结果较差 (b) 对于正方形柱体绕流流速分布, 计算结果和实验结果对比如图 3 所示 ( 实验资料来源 : 文献 [10]) 从对比中可以看出, 由于钝体绕流分析, 目前各湍流模型分析都有较大难度, 相对而言,SZL 模型和 GIR 模型和实验结果更为接近一些, 而其他模型差距更大 (c) 对于圆柱体绕流流场分析, 计算结果和实验结果对比如图 4 所示 ( 实验资料来源 : 文献 [5]) 从对比中可以看出, 正风向各模型和实验结果基本一致, 横风向 SST SKW 和 SKE 模型和实验结果较为接近, 其他模型在计算结果偏小如果比较不同模型的流速云图, 可以发现,RNG 模型 GIR 模型和 SZL 模型 ( 图 5) 的计算结果和实验结果比较相似, 而其他模型计算结果则差异较大 SKE NKE RNG GIR 相对速度 (V/V0) SZL SST SKW 风洞试验相对距离 (D/H) 图 2 房屋表面风压分布对比 Fig. 2 Comparison for wind pressure on the building 3 图 3 正方形柱体后部流速分布对比 Fig.3 Comparison for wind speed along the centerline of square cylinder

4 相对压力 (P/P0) SKE RNG NKE GIR SZL SKW SST 实验角度 α 图 4 圆柱体表面风压对比 Fig. 4 Wind pressure on the round cylinder 图 5 SZL 模型流速分布图 Fig. 5 Wind speed distribution with SZL model 为了进一步论证不同湍流模型在耦合分析中的差异, 需要进行考虑流固耦合效应的模型对比分析因为目前还很难找到考虑流固耦合效果的标准实验, 因此本文以文献 [11] 中提供的带隔板的二维流场问题作为对象进行了对比分析 ( 图 6) 不同湍流模型分析结果和文献 [11] 中结果的差异如图 7 所示其中,GIR 模型由于计算结果发散, 未能得到其振动频率从对比中可以看出,RNG 和 SZL 模型的计算结果和文献 [11] 中结果最为接近 10 U? d 8 6 最大位移误差振动频率误差 L E=3.5GPa ρ s =1200kg/m 3 L=1m d=10mm U? =10m/s 相对误差 SKE RNG NKE GIR SZL SKW SSL -10 模型名称图 6 带隔板二维流场问题 Fig. 6 2D flow problem with an isolate plate 图 7 不同湍流模型结果与文献 [11] 误差对比 Fig. 7. Comparison for errors between different models with Ref [11] 综合上述对比分析, 各模型计算效果对比如表 1 所示表 1 不同湍流模型计算结果对比 Table 1 Comparison for the results with different turbulent models SKE RNG NKE GIR SZL SKW SST 正方形房屋风压分布正方形柱体绕流圆柱体绕流 ( 压力分布 ) 圆柱体绕流 ( 流场分布 ) 带隔板流固耦合符号说明 : : 较好 ; : 一般 ; : 较差从表中可以看出, 对于本文研究的 CFRP 索流固耦合振动问题而言,RNG 模型和 SZL 模型相对比较合适, 而其他模型的问题相对要大一些, 该结论也和一些理论分析相一致 4

5 3 CFRP 索风振计算 3.1 计算模型取 CFRP 索长度为 100m, 直径为 10cm, 弹性模量为 230GPa, 密度为 1400kg/m 3 设定索内初始应力为 200MPa 或 500MPa 为了考虑 CFRP 索和周围空气的流固耦合作用, 在索周围取 4m 10m 100m 的立方体空气柱整个计算模型如图 8 所示空气柱索 (a) 模型尺寸 (b) 有限元网格 ( 局部 ) 图 8 CFRP 索风振计算模型 Fig. 8 Computational model for the wind-induced vibration of CFRP cable. 分别进行了一个平稳风场工况和三个脉动风场工况的计算, 平均风速均为 30m/s, 并和钢索进行了对比 3.2 平稳风场计算结果当初始张拉应力分别为 200MPa 和 500MPa 时, 分别采用 RNG 模型和 SZL 模型, 对 CFRP 索和钢索的流固耦合风振进行了分析分析发现 RNG 模型和 SZL 模型计算结果非常接近, 以 SZL 模型计算结果为例,CFRP 索和钢索在平稳风场作用下最大变形的时程如图 9 所示从图中可以得到以下结论 : (1) 周围空气对索的气弹阻尼效果非常明显, 虽然计算模型中未设置任何阻尼, 但是由于索周面空气的影响, 无论是钢索还是 CFRP 索的振动都表现出明显的随时间衰减以 CFRP 索为例, 将 CFRP 索迎风面 ( 正面风压 ) 和索的振动时程放在一起, 如图 10 所示, 就可以明显看出, 周围空气流场对 CFRP 索振动有一个明显的反向抑制作用并且, 索的振幅越大振动频率越快, 周围空气流场的抑制作用也就越明显考虑周围空气的气弹阻尼作用, 从图 9 索的振动中可以估算出 : 对于 CFRP 索, 当初始应力为 200MPa 时, 阻尼比大约为 0.053, 而当初始应力为 500MPa 时, 由于索振动频率更快, 气弹阻尼也就更大, 阻尼比达到了约 0.11 对于钢索也有类似现象当初始应力为 200MPa 时, 钢索阻尼比约为 0.015, 而初始应力为 500MPa 时, 阻尼比达到了 (2) 数值计算得到的索的一阶振动频率和理论频率对比如表 2 所示, 从表中可以看出, 计算频率均小于理论频率值, 说明考虑流固耦合作用后, 周围空气和索一起振动, 气弹质量降低了索的振动频率且由于 CFRP 索更轻, 周围空气质量影响也就更大, 相应的数值计算得到的频率和理论值的差异也就更大一些但是总的说来, 无论是钢索还是 CFRP 索, 频率和理论值的差别都小于 5%, 说 5

6 明对于索的振动而言, 气弹质量和气弹刚度的影响相对比较小 (3) 在本次分析中, 索主要为顺风向振动, 横风向的振动始终很小, 本文认为, 这主要是由于目前的湍流模型对侧风向的风压计算精度不高所致 ( 图 4) Steel CFRP % 3 Steel 25% CFRP 2 15% 1 5% (a) 初始应力 =200MPa (b) 初始应力 =500MPa 图 9 平稳风场作用下索最大相对变形时程曲线 (dx 为索最大变形,D 为索直径 ) Fig. 9 The time-history curves of maximal relative deformation of cables in stable wind field 压力 (Pa) 正面风压 8 相对变形相对位移 (dx/d) 压力 (Pa) % 2 15% 1 正面风压 5% 相对变形相对位移 (dx/d) (a) 初始应力 =200MPa (b) 初始应力 =500MPa 图 10 CFRP 索迎风面 ( 正面风压 ) 及变形时程曲线 Fig. 10 The time-history curves of wind pressure and deformation of CFRP cables 表 2 索的一阶自振频率 Table 2 Free vibration frequencies of cables (First model) 数值计算理论频率 ( 不考虑阻尼 ) 1 F = 2L T M 理论频率 ( 考虑阻尼 ) 2 ( 1 ) 1 T ξ F = 2L M CFRP 索钢索 CFRP 索钢索 CFRP 索钢索 200MPa MPa 注 :L 为索长度,T 为索拉力,M 为索单位长度质量, ξ 为阻尼比 3.3 脉动风场计算结果采用谐波合成法 (WAWS), 根据 Davenport 谱, 生成了 3 条随机脉动风时程, 计算 CFRP 索和钢索在脉动风荷载作用下的振动需要说明的是, 本文使用的脉动风未考虑风压的高度变化 6

7 计算得到的索在脉动风场作用下的变形时程曲线如图 11 所示, 由于脉动风场本身就很复杂, 因此从变形时程曲线上难以直接判断钢索和 CFRP 索的优劣因此, 对图 11 中的 6 组 12 条曲线的进行统计, 并用正态分布曲线去拟合, 得到各条时程曲线的均值和方差如表 3 所示从表 3 中可以看出, 当初始应力 =200MPa 时, 无论是钢索还是 CFRP 索振动的离散性都比较大, 二者对比没有明显优劣之分而当初始应力达到 500MPa 时,CFRP 索振动的方差均明显小于钢索分析其原因包括 :1CFRP 索质量轻频率高, 吸收脉动风中低频能量较少 ;2CFRP 索气弹阻尼高于钢索, 进而振动中能量耗散更多因此, 虽然从概念上,CFRP 索因为质量更轻, 风致振动问题可能会更严重, 但是在实际工程中, 如果设计得当的话, CFRP 索风振性能完全可能比钢索更好脉动风时程 01CFRP -10 脉动风时程 01Steel 5-5 脉动风时程 02CFRP -10 脉动风时程 02Steel (a) 初始应力 =200MPa, 脉动风时程 01 (b) 初始应力 =200MPa, 脉动风时程脉动风时程 03CFRP -10 脉动风时程 03Steel 脉动风时程 01CFRP -2 脉动风时程 01Steel (c) 初始应力 =200MPa, 脉动风时程 03 (d) 初始应力 =500MPa, 脉动风时程脉动风时程 02CFRP -2 脉动风时程 02Steel 脉动风时程 03CFRP 脉动风时程 03Steel (e) 初始应力 =500MPa, 脉动风时程 02 (f) 初始应力 =500MPa, 脉动风时程 03 图 11 脉动风场下索最大相对变形时程曲线 Fig.11 Time-history curves of maximal relative deformation of cables in fluctuant wind field 7

8 表 3 脉动风场相对变形时程统计结果 Table 3 Statistic results of relative deformation of cables in fluctuant wind field 初始应力 200MPa 500MPa 脉动风时程 CFRP 索钢索均值方差方差 / 均值均值方差方差 / 均值脉动风时程脉动风时程脉动风时程脉动风时程脉动风时程脉动风时程结论本文利用通用有限元软件 ANSYS 提供的流固耦合计算功能, 对 CFRP 索风振问题进行了数值分析, 得到以下结论 : (1) ANSYS 提供的湍流模型中, 比较适合分析索风振问题的湍流模型是 RNG 模型和 SZL 模型 ; (2) 在均匀风场作用下,CFRP 索的气弹阻尼效果明显, 气弹质量和气弹刚度效果较小 ; (3) CFRP 索的气弹阻尼比要高于同尺寸的钢索, 且索的振幅越大频率越高, 气弹阻尼效果越明显 ; (4) 在脉动风场作用下,CFRP 索由于自振频率比同尺寸的钢索要高, 且气弹阻尼比要大, 如果设计合理, 其风振效果可能会优于同尺寸的钢索该结论对 CFRP 索在大跨结构中的应用有一定意义参考文献 1. 张元凯, 肖汝诚, FRP 材料在大跨度桥梁结构中的应用展望, 公路交通科技,2004, 21(4), (Zhang YK., Xiao RC., Application of FRP material in long span bridge, Journal of Highway and Transportation Research and Development, 2004, 21(4), ) 2. 汪胜义, 碳纤维复合材料 (CFRP) 在桥梁工程中的应用及前景, 重庆交通学院学报,2005, 24(1), (Wang SY., The application and future of carbon fiber reinforced polymer in bridge structure, journal of Chongqing Jiaotong University, 2005, 24(1), ) 3. 梅葵花, 吕志涛,CFRP 斜拉索的静力特性分析, 中国公路学报,2004,17(2), (Mei KH., Lu ZT., Static characteristic analysis of CFRP cables, China Journal of Highway and Transport, 2004, 17(2), ) 4. 武岳, 考虑流固耦合作用的索膜结构风致振动相应研究, 哈尔滨工业大学博士学位论文, 黄本才, 结构抗风分析原理及应用, 上海 : 同济大学出版社, 王福军, 计算流体动力学分析 CFD 软件原理与应用, 北京 : 清华大学出版社, Tsuchiya M., Murakami S., Mochida A., Kondo K., Ishida Y., Development of a new k-ε model for flow and pressure fields around bluff body, Journal of Wind Engineering and Industrial Aerodynamics, 1997, 67&68: Murakami S., Overview of turbulence models applied in CWE-1997, Journal of Wind Engineering and Industrial Aerodynamics, 1998, 74-76: Jeong UY., Koh HM., Lee HS., Finite element formulation for the analysis of turbulent wind flow passing bluff structures using the RNG k-ε model. Journal of Wind Engineering and Industrial Aerodynamics. 2002, 8

9 90: Murakami S., Mochida A., On turbulent vortex shedding flow past 2D square cylinder predicted by CFD, Journal of Wind Engineering and Industrial Aerodynamics, 1995, 54&55: Gluck M., Breuer M., Durst F., Halfmann A., Rank E., Computation of fluid-structure interaction on lightweight structures, Journal of Wind Engineering and Industrial Aerodynamics, 2001, 89:

2 ( 自然科学版 ) 第 20 卷波 ). 这种压缩波空气必然有一部分要绕流到车身两端的环状空间中, 形成与列车运行方向相反的空气流动. 在列车尾部, 会产生低于大气压的空气流

2 ( 自然科学版 ) 第 20 卷波 ). 这种压缩波空气必然有一部分要绕流到车身两端的环状空间中, 形成与列车运行方向相反的空气流动. 在列车尾部, 会产生低于大气压的空气流第 20 卷第 3 期 2014 年 6 月 ( 自然科学版 ) JOURNAL OF SHANGHAI UNIVERSITY (NATURAL SCIENCE) Vol. 20 No. 3 June 2014 DOI: 10.3969/j.issn.1007-2861.2013.07.031 基于 FLUENT 测轨道交通隧道中电波折射率结构常数张永