CN104408295B - 一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 - Google Patents
一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 Download PDFInfo
- Publication number
- CN104408295B CN104408295B CN201410629276.6A CN201410629276A CN104408295B CN 104408295 B CN104408295 B CN 104408295B CN 201410629276 A CN201410629276 A CN 201410629276A CN 104408295 B CN104408295 B CN 104408295B
- Authority
- CN
- China
- Prior art keywords
- wave
- wind
- simulation
- wind speed
- coupling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 68
- 230000008878 coupling Effects 0.000 title claims abstract description 50
- 238000010168 coupling process Methods 0.000 title claims abstract description 50
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001228 spectrum Methods 0.000 claims abstract description 41
- 230000008569 process Effects 0.000 claims abstract description 13
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 238000005516 engineering process Methods 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000003595 spectral effect Effects 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 230000003993 interaction Effects 0.000 abstract description 4
- 230000001360 synchronised effect Effects 0.000 abstract description 3
- 241000282485 Vulpes vulpes Species 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 210000001367 artery Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Wind Motors (AREA)
Abstract
本发明提供了一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法,该方法解决了模拟过程中遇到的风与浪耦合、风与结构耦合和波浪与结构耦合三个方面的环境与结构相互作用多重耦合问题。首先,通过脉动风速谱和波浪频谱的耦合来实现风、浪两种随机过程的时域内同步耦合模拟;其次,通过雷诺应力张量修正结构周边风场的脉动风速谱来实现风与结构的耦合模拟;最后,通过MacCamy‑Fuchs绕射理论来实现波浪与结构的耦合模拟。该方法综合运用了计算流体力学(CFD)方法、快速傅里叶变换技术(FFT)、特征正交分解法(POD)和频率插值法,该方法计算效率高、模拟过程简单且物理意义明确,可推广到其他类似海洋结构的风‑浪耦合作用荷载数值模拟中。
Description
技术领域
本发明涉及一种大跨桥梁下部结构风-浪耦合作用荷载数值模拟方法,尤其是涉及数值模拟过程中风与浪耦合、风与结构耦合以及波浪与结构耦合三个方面的环境与结构相互作用多重耦合问题的解决方法。
背景技术
大跨桥梁下部结构在施工和运营过程中将承受复杂的、随时空变化的随机环境荷载,主要包括风荷载和波浪荷载。大量研究表明风、浪两种环境要素相互影响,存在较强的耦合特性。风、浪的相互耦合使得海面风场特性显著区别于陆地风场特性,主要表现在风场壁面(即波面)光滑但不断演变、近地面粗糙度减小、摩擦风速随时空变化、脉动风速谱的频率向低频偏移。同时,风是浪的重要动力来源,波浪形态受到风速控制,随机波浪频谱是风速和风区的函数。
风、浪耦合作用下结构动力响应的精确分析依赖于输入环境荷载数据的可靠性。一般而言,环境荷载作用数据来源于现场实测或数值模拟;前者可靠度好,但对设备仪器要求较高且耗时较长;后者适应性强、耗时短,但鲜有考虑风、浪两者的相互耦合。同时,置于风场和波浪中的大型结构本身对风、浪的干扰以及风-结构、波浪-结构之间的相互作用不可忽略,风-浪耦合作用荷载数值模拟需要解决这三方面的耦合问题。为实现大跨桥梁下部结构风-浪耦合作用荷载的高效模拟,发展一种以海面风场和随机波浪的概率统计及物理机制描述为特征的数值模拟程序十分必要。
发明内容
本发明的目的在于提供一种以海面风场和随机波浪的概率统计及物理机制描述为特征的风-浪耦合作用荷载数值模拟程序,在数值算法上解决模拟过程中遇到的风与浪耦合、风与结构耦合和波浪与结构耦合三种不同的相互作用耦合问题,为风-浪耦合作用下大跨桥梁下部结构动力响应分析及设计奠定了基础。
为此,本发明采用如下技术方案:
首先,通过脉动风速谱和波浪频谱的谱矩阵耦合来实现风、浪两种随机过程的时域内同步耦合模拟;其次,通过雷诺应力张量修正结构周边风场的脉动风速谱来实现风与结构的耦合模拟;最后,通过MacCamy-Fuchs绕射理论来实现波浪与结构的耦合模拟。
本发明提供一种大跨桥梁下部结构风-浪耦合作用荷载数值模拟方法,该方法包括以下步骤:
(1)确定大跨桥梁下部结构脉动风速模拟点和随机波浪模拟点的坐标;
(2)根据海面风场和随机波浪频谱特性,利用Ochi脉动风速谱和JONSWAP波浪频谱生成耦合谱矩阵S:
式中z为海面高度,为脉动风速u的谱分量,n为模拟点数,Sη为波浪频谱分量,η为随机波浪的波面。耦合谱矩阵S的建立为风、浪两种随机过程的时域内同步耦合模拟奠定了算法基础;
(3)建立大跨桥梁下部结构数值模拟风场,利用雷诺应力模型(RMS)完成该风场的计算流体力学(CFD)数值模拟,提取脉动风速模拟点处的雷诺应力张量来修正耦合谱矩阵S中的脉动风速谱分量修正后的自谱为:
互谱为:
式中τii为雷诺应力张量,X为莫宁坐标,f为频率,Uz为海面高度z处的平均风速,Coh(f)为雷诺应力张量修正后的点相干函数,Lu为湍流积分尺度。修正后的耦合谱矩阵S考虑了风与结构之间的相互耦合;
(4)通过对数等间距插值设置耦合谱矩阵S的圆频率插值基点
式中Δω=ωu/Nω,ωu为截止圆频率,Nω为频域离散点数,为圆频率插值基点数;
(5)在圆频率插值基点处对耦合谱矩阵S进行特征正交分解(POD),计算得到各基点处的特征值和特征向量通过插值公式获得任意频率点ωj(j=1,…,Nω)处的特征值和特征向量θj:
式中floor(sj)和ceil(sj)分别表示不大于sj的最大整数和不小于sj的最小整数,round(sj)为最接近sj的整数;
(6)合成由脉动风速和随机波浪组成的随机矢量过程V,引入快速傅里叶变换技术(FFT)提高合成效率:
式中V(tm)={u1(tm),…,un(tm),η(tm)},tm=mΔt(m=1,2,...,Nt),Δt为时间步长,Nt为时域离散点数,ωj=(j-1/2)Δω,Ns为频率点ωj处特征值和特征向量θj的截断数,为随机相位角,服从[0,2π]间的独立均匀分布,Δω为圆角频率步长;ηj,k为频率点ωj处的第k阶随机波浪的波面,θj,k为频率点ωj处的第k阶特征向量;{u1(tm),…,un(tm)}为tm时刻脉动风速模拟点1~n的风速取值,η(tm)为tm时刻随机波浪的波面取值;
(7)在{u1(tm),…,un(tm)}的基础上利用准定常假定计算脉动风速模拟点处大跨桥梁下部结构的脉动风压。在η(tm)的基础上利用MacCamy-Fuchs绕射理论计算大跨桥梁下部结构表面的随机波浪压力。通过绕射理论,波浪与结构之间的耦合作用得到充分考虑。
本方发明以海面风场和随机波浪两种环境要素的概率统计特性及物理机制为基础,结合计算流体力学(CFD)方法、特征正交分解法(POD)、快速傅里叶变换技术(FFT)和频率插值法,在数值算法上高效地解决了模拟过程中需要考虑的风与浪耦合、风与结构耦合和波浪与结构耦合三种耦合问题,可有效实现大跨桥梁下部结构和类似大型海洋结构风-浪耦合作用荷载数值模拟,该方法计算效率高、模拟过程简单且物理意义明确,可推广到其他类似海洋结构的风-浪耦合作用荷载数值模拟中。
附图说明
图1为本发明的流程图。
图2为某大跨桥梁下部结构示意图和脉动风速、随机波浪模拟点分布图。
图3为某大跨桥梁下部结构数值模拟风场图。
图4为模拟点1处的脉动风速模拟结果图。
图5为随机波浪的波面模拟结果图。
图6为模拟点2处的随机波浪压力模拟结果图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明,但该实例不应理解为对本发明的限制。
本发明阐述一种大跨桥梁下部结构风-浪耦合作用荷载的数值模拟方法,其流程如图1所示,该方法包括以下步骤:
(1)确定大跨桥梁下部结构脉动风速模拟点和随机波浪模拟点的坐标;
(2)根据海面风场特性和随机波浪特性,利用Ochi脉动风速谱和JONSWAP波浪频谱生成耦合谱矩阵S:
式中为脉动风速u的谱分量,n为模拟点数,Sη为波浪频谱分量,η为随机波浪的波面。耦合谱矩阵S的建立为风、浪两种随机过程的时域内同步耦合模拟奠定了基础;
(3)建立大跨桥梁下部结构数值模拟风场,利用雷诺应力模型(RSM)完成该风场的计算流体力学(CFD)数值模拟,提取脉动风速模拟点处的雷诺应力张量来修正耦合谱矩阵S中的脉动风速谱分量修正后的自谱为:
互谱为:
式中z为海面高度,τii为雷诺应力张量,X为莫宁坐标,f为频率,Uz为海面高度z处的平均风速,Coh(f)为雷诺应力张量修正后的点相干函数,Lu为湍流积分尺度。修正后的耦合谱矩阵S考虑了风与结构之间的相互耦合;
(4)通过对数等间距插值设置耦合谱矩阵S的圆频率插值基点
式中Δω=ωu/Nω,ωu为截止圆频率,Nω为频域离散点数,为圆频率插值基点数;
(5)在圆频率插值基点处对耦合谱矩阵S进行特征正交分解(POD),计算得到各基点处的特征值和特征向量通过插值公式获得任意频率点ωj(j=1,…,Nω)处的特征值和特征向量θj:
式中floor(sj)和ceil(sj)分别表示不大于sj的最大整数和不小于sj的最小整数,round(sj)为最接近sj的整数;
(6)合成由脉动风速和随机波浪组成的随机矢量过程V,引入快速傅里叶变换技术(FFT)提高合成效率:
式中V(tm)={u1(tm),…,un(tm),η(tm)},tm=mΔt(m=1,2,...,Nt),Δt为时间步长,Nt为时域离散点数,ωj=(j-1/2)Δω,Ns为频率点ωj处特征值和特征向量θj的截断数,为随机相位角,服从[0,2π]间的独立均匀分布,Δω为圆角频率步长;ηj,k为频率点ωj处的第k阶随机波浪的波面,θj,k为频率点ωj处的第k阶特征向量;{u1(tm),…,un(tm)}为tm时刻脉动风速模拟点1~n的风速取值,η(tm)为tm时刻随机波浪的波面取值;
(7)在{u1(tm),…,un(tm)}的基础上利用准定常假定计算脉动风速模拟点处大跨桥梁下部结构的脉动风压。在η(tm)的基础上利用MacCamy-Fuchs绕射理论计算大跨桥梁下部结构表面的随机波浪压力。通过绕射理论,波浪与结构之间的耦合作用得到充分考虑。
下面例举一个实施例。
图2为某大跨桥梁下部结构示意图,包括水面以上的桥塔部分和水面以下的大直径圆柱形沉井基础部分。
(1)确定该大跨桥梁下部结构脉动风速模拟点和随机波浪模拟点的坐标,其中脉动风速模拟点位于桥塔部分,每个模拟点的距离约为10m,共81个,随机波浪模拟点的位于水面,模拟点数为1,模拟点的具体分布见图2;
(2)根据海面风场特性和随机波浪特性,利用Ochi脉动风速谱和JONSWAP波浪频谱生成耦合矩阵S,其中n=81,S为82阶方阵;
(3)建立大跨桥梁下部结构数值模拟风场,如图3所示。沉井基础不是海面风场的组成部分,因此未包含于数值模拟风场。利用雷诺应力模型(RSM)完成该风场的计算流体力学(CFD)数值模拟,提取各脉动风速模拟点处的雷诺应力张量τii来修正耦合谱矩阵S中的脉动风速谱分量
(4)通过对数等间距插值设置耦合谱矩阵S的圆频率插值基点其中频域离散点数Nω=1024,圆频率插值基点数截止圆频率ωu=3π;
(5)在各频率插值基点处对耦合谱矩阵S进行特征正交分解(POD),计算得到各基点处的特征值和特征向量正交分解的频率截断数Ns=16。利用公式(16)和(17)得到任意频率点ωj处的特征值和特征向量θj;
(6)利用公式(18)合成由脉动风速和随机波浪组成的随机过程V,其中时间步长Δt=0.1s,时域离散点数Nt=6500。图4为模拟点1处的脉动风速模拟结果,图5为随机波浪的波面模拟结果;
(7)利用准定常假定计算脉动风速模拟点处的脉动风压,利用MacCamy-Fuchs绕射理论计算沉井基础表面的随机波浪压力。图6为沉井基础模拟点2处的随机波浪压力。
如上所述,尽管参照特定的优选实施例已经表示和表述了本发明,但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下,可对其在形式上和细节上作出各种变化。
Claims (1)
1.一种大跨桥梁下部结构风-浪耦合作用荷载数值模拟方法,其特征在于,该方法包括以下步骤:
第一步:确定大跨桥梁下部结构脉动风速模拟点和随机波浪模拟点的坐标;
第二步:根据海面风场特性和随机波浪特性,利用Ochi脉动风速谱和JONSWAP波浪频谱生成耦合谱矩阵S:
式中为脉动风速u的谱分量,n为模拟点数,Sη为波浪频谱分量,η为随机波浪的波面;
第三步:建立大跨桥梁下部结构周边数值模拟风场,利用雷诺应力模型(RSM)完成该风场的计算流体力学(CFD)数值模拟,提取脉动风速模拟点处的雷诺应力张量来修正耦合谱矩阵S中的脉动风速谱分量修正后的自谱为:
互谱为:
式中z为海面高度,τii为雷诺应力张量,X为莫宁坐标,f为频率,Uz为海面高度z处的平均风速,Coh(f)为雷诺应力张量修正后的点相干函数,Lu为湍流积分尺度,修正后的耦合谱矩阵S考虑了风与结构之间的相互耦合;
第四步:通过对数等间距插值设置耦合谱矩阵S的圆频率插值基点
式中Δω=ωu/Nω,ωu为截止圆频率,Nω为频域离散点数,为圆频率插值基点数;
第五步:在圆频率插值基点处对耦合谱矩阵S进行特征正交分解(POD),计算得到各基点处的特征值和特征向量通过插值公式获得任意频率点ωj,j=1,…,Nω,处的特征值和特征向量θj:
式中floor(sj)和ceil(sj)分别表示不大于sj的最大整数和不小于sj的最小整数,round(sj)为最接近sj的整数;
第六步:合成由脉动风速和随机波浪组成的随机矢量过程V,引入快速傅里叶变换技术(FFT)提高合成效率:
式中V(tm)={u1(tm),…,un(tm),η(tm)},tm=mΔt,m=1,2,…,Nt,Δt为时间步长,Nt为时域离散点数;ωj=(j-1/2)Δω,Ns为频率点ωj处特征值和特征向量θj的截断数,为随机相位角,服从[0,2π]间的独立均匀分布,Δω为圆角频率步长;ηj,k为频率点ωj处的第k阶随机波浪的波面,θj,k为频率点ωj处的第k阶特征向量;{u1(tm),…,un(tm)}为tm时刻脉动风速模拟点1~n的风速取值,η(tm)为tm时刻随机波浪的波面取值;
第七步:在{u1(tm),…,un(tm)}的基础上利用准定常假定计算脉动风速模拟点处大跨桥梁下部结构的脉动风压,在η(tm)的基础上利用MacCamy-Fuchs绕射理论计算大跨桥梁下部结构表面的随机波浪压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410629276.6A CN104408295B (zh) | 2014-11-10 | 2014-11-10 | 一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410629276.6A CN104408295B (zh) | 2014-11-10 | 2014-11-10 | 一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104408295A CN104408295A (zh) | 2015-03-11 |
CN104408295B true CN104408295B (zh) | 2017-07-07 |
Family
ID=52645926
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410629276.6A Active CN104408295B (zh) | 2014-11-10 | 2014-11-10 | 一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408295B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104778331B (zh) * | 2015-04-24 | 2017-12-05 | 浙江工业大学 | 一种大跨桥梁监测数据空间插值方法 |
CN105303055B (zh) * | 2015-11-16 | 2018-06-08 | 国网河南省电力公司郑州供电公司 | 风振激励作用下特高压输电线路非线性动力学分析方法 |
CN105653817A (zh) * | 2016-01-20 | 2016-06-08 | 吴本刚 | 一种具有精确识别功能的变电站 |
CN107247687B (zh) * | 2017-05-18 | 2021-06-08 | 西南交通大学 | 一种基于特征正交分解的非平稳随机过程快速模拟方法 |
CN109726422B (zh) * | 2018-07-17 | 2020-09-08 | 中国科学院力学研究所 | 基于绕射理论的圆柱阵列波浪力曲线计算步长确定方法 |
CN109635509A (zh) * | 2019-01-15 | 2019-04-16 | 西南交通大学 | 一种基于风浪流荷载组合的跨海桥梁动力响应计算方法 |
CN109657409B (zh) * | 2019-01-15 | 2022-04-29 | 西南交通大学 | 一种基于风浪联合分布极端响应的跨海桥结构优化方法 |
CN114166416B (zh) * | 2021-12-08 | 2023-05-30 | 中国船舶科学研究中心 | 一种湍流脉动压力的修正方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2300706A2 (en) * | 2008-04-24 | 2011-03-30 | Ocean Wave Rocker AS | Energy system |
CN103969010B (zh) * | 2013-01-24 | 2016-12-07 | 中交公路规划设计院有限公司 | 桥梁风浪流耦合场、弹性模型及动态响应试验测试系统 |
-
2014
- 2014-11-10 CN CN201410629276.6A patent/CN104408295B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104408295A (zh) | 2015-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408295B (zh) | 一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 | |
Chen et al. | A deterministic–stochastic hybrid model of downbursts and its impact on a cantilevered structure | |
CN110008509B (zh) | 一种考虑背景流场下的内孤立波作用力特性分析方法 | |
Gomes et al. | Wave power extraction of a heaving floating oscillating water column in a wave channel | |
Ke et al. | A new methodology for analysis of equivalent static wind loads on super-large cooling towers | |
CN105224715A (zh) | 一种山区地貌下强风三维脉动风场综合模拟方法 | |
de Leon et al. | Simulation of irregular waves in an offshore wind farm with a spectral wave model | |
Allaerts et al. | Sensitivity and feedback of wind-farm-induced gravity waves | |
Feng et al. | Development of an unstructured-grid wave-current coupled model and its application | |
Wang et al. | Assessment of wind energy and wave energy resources in Weifang sea area | |
Lock et al. | Demonstration of a cut-cell representation of 3D orography for studies of atmospheric flows over very steep hills | |
Song et al. | Investigation on wind-induced aero-elastic effects of tall buildings by wind tunnel test using a bi-axial forced vibration device | |
Wang | Efficient prediction of wave energy converters power output considering bottom effects | |
CN115470731A (zh) | 一种基于微气象的风场风电功率预测方法及系统 | |
Fan et al. | Optimized design of equivalent truncated mooring system based on similarity of static and damping characteristics | |
CN111651920B (zh) | 塔线体系下大跨越输电塔风致机械能的传递分析方法 | |
Weller et al. | Challenges facing adaptive mesh modeling of the atmosphere and ocean | |
Altaf et al. | Inverse shallow-water flow modeling using model reduction | |
CN103065044B (zh) | 在分形海面背景下畸形波的模拟方法 | |
CN114492220B (zh) | 基于OpenFOAM平台的海洋大气边界层流动特性的预测方法及系统 | |
Ferri et al. | Validation of a wave-body interaction model by experimental tests | |
Liu et al. | Dynamic elastic response testing method of bridge structure under wind-wave-current action | |
Zhou et al. | Numerical predictions of wave impacts on the supporting structures of Shanghai Donghai-Bridge offshore wind turbines | |
Nishi et al. | Growth and coexistence of structural and lift force modes in vortex-induced vibration of a flexible riser | |
Long et al. | Some experimental calculations for 3D currents in the strong upwelling region of Southern Central Vietnam using finite element method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |