波浪力计算表格
- 格式:xlsx
- 大小:259.98 KB
- 文档页数:11
设计高潮位2.87△h= 2.26风向组NE-ENE 42.65风向组N-NNE 2.35风向组Vc=13.3Cv=0.31Vc=15Cv=0.32Vc=Kp= 1.85Kp= 1.88Kp=Vap=30.24df=20.09Vap=34.44df=Vap=极限风区长实测风区长5.2极限风区长实测风区长6极限风区长风向组E-ESE 87.65风向组Vc=12.5Cv=0.35Vc=Kp= 1.99Kp=Vap=30.55df=Vap=极限风区长实测风区长5极限风区长L塘(假设)240塘前涂面高-0.550251030d前 3.372015dF/Lo 0.8480Ks=14130.5H'o=0.972030.5d前/Lo0.1133Ks=0.91353025H'前=0.882025L/Lo 0.8541115.5L=25.401010.5H*=0.262099.504.5H1% 1.7620.09H2% 1.66H13% 1.30H5% 1.50d前/Lo0.1133Hb/DbHb=5.26查表5.2.1-2Kv= 1.285Lo/H1%16.90查附图十一Ro= 1.32当m=0.4时海塘外护坡面取用砼K△0.8Ro=1.32R 2%=2.39Kf=1R 13%=1.36Kf=0.77假定压载宽=15高=2B/L=0.591L/H=19.491d1/H=0.779Ky=0.87R 2%(实)=2.08R 13%(实)=1.18R 2%(不)= 5.65 6.1R 2%(允)=4.451.15T L=25.40H塘/d前0.3H塘/L塘0.0348A=0.0132B=27.429查表5.2.1-3,F1%查表5.2.1-3,F13%三江查表5.2.9-1塘前有压载修正塘顶高程查表6.2.2无风越浪实际采用风区长=主风向左45度5.2实际采用风区长=波浪爬高V/(d前*9.81)^0.56实际采用风区右45度实际采用风区长=5实际采用风区越顶水量q=0.00032Hc=3.23Wf=2θ=68.1986K'= 5.590越顶水量q=0.0018风速大于26.8有风越浪N-NNE 12.65风向组N-NNE 27.6515Cv=0.33Vc=15Cv=0.331.92Kp= 1.9235.15df=Vap=35.15df=实测风区长5.25极限风区长实测风区长 5.25NE-ENE 72.65风向组NE-ENE 57.6513.5Cv=0.3Vc=13.3Cv=0.31.82Kp= 1.8230.20df=Vap=29.77df=实测风区长6极限风区长实测风区长65.13180S-SSW实际采用风区长=6左15度实际采用风区长= 5.25左30度 5.25用风区长=风区水深Df=右30度右15度用风区长=6。
波浪对桩柱作用力的计算分析港口码头、跨海桥梁、海洋工程中经常遇到波浪力对桩柱的作用,水中的桩柱结构所受的水平力主要来自波浪的作用。
对于一般的桩柱结构,当桩柱的直径D与波浪的波长L的比值D/L<0.2时,称此结构为小尺度结构物;当D/L>0.2时,则称为大尺度结构物[1]。
港口码头、跨海桥梁、海洋工程常用的桩柱结构,一般为小尺度结构物,因此本文重点探究小尺度结构物的波浪力计算。
1 波浪力的计算方法国内外对波浪力的计算理论相差不大,都是基于Morison公式[2]內对波浪力的计算主要有《海港水文规范》规定的方法,国外对波浪力的计算有美国的API、挪威的DNV等标准提出的方法。
基于Morison理论的波浪力计算公式如下:(公式-1)(公式-2)(公式-3)上式中:PD是波浪力的速度分力(也叫拖曳力);PI是波浪力的惯性分力;CD是速度力系数;CM是质量力系数;D是桩柱直径;dz是桩柱上每一小分段的长度;A是桩柱的断面积。
u、u/t——分别为水质点轨道运动的水平速度和水平加速度;是圆频率;t是时间,当波峰通过柱体中心线时t=0;H、T分别是波高和波周期。
美国的API规范给出拖曳力系数取0.6~1.0,惯性力系数取2.0。
国内的《海港水文规范》给出拖曳力系数取1.2,惯性力系数取2.0。
《海港水文规范》提出的拖曳力系数和惯性力系数比较适合我国的情况,因此本文以《海港水文规范》给出拖曳力系数1.2、惯性力系数2.0进行波浪力计算。
2 波浪力随高度和时间的变化趋势本文模拟了一个我国东部沿海的项目的桩基,该项目水深10m,采用垂直的钢管桩,桩体直径2.5m。
计算海浪的波长72m,波高3.9m,波周期7.8s。
采用Excel对周期T内的各时间点、沿桩柱各高度的拖曳力和惯性力的分力分别进行计算,然后用数值积分的方法将各分力沿桩柱高度进行积分,得出各时间点作用在桩柱上的波浪力的合力。
为考察波浪力随桩体高度变化趋势,按照公式1~3,将桩体在平均海平面以上分为10等份,在平均海平面以上分为2等份,进行拖曳力和惯性力计算。
共 19 页平台波浪力计算书二、平台基本数据和环境条件平台主尺度:长57.75 m,宽34.5 m,型深5 m桩腿尺度:Φ2.3*45 m桩腿间距:横向28.5 m 纵向44.55 m基线距海底:28.5 m设计水深:d = 20 m最大天文潮高:d t= 4 m最大波高:H = 6.5 m波浪周期:T = 7.0 ~ 13.0S最大流速:U = 1.5 m/S最大风速:V = 41.15 m/S作业工况波高:H1= 4 m作业工况风速:V1 = 25.8 m/S迁航波高:H C = 4 m迁航航速:V C = 4 knot迁航时桩腿在基线下:0.75m (Min.)图1 荷载方向三、自存工况(一)波浪荷载和海流荷载1.计算原理因为桩腿的直径和波长之比小于0.2,桩腿所受的波浪与海流荷载按Morison公式进行计算,单根桩腿单位长度所受的波流力为:F W=F D+F I (1)=C d×ρ×D×| U|×U/2+C m×ρ×π×D2/4×a=C d×ρ×D×( | u+v| )×( u+v )/2+C m×ρ×π×D2/4×a式中: F W =波流力,N。
F D=阻力,N。
F I=惯性力,N。
ρ=海水密度。
ρ=1. 025*103 kg/m3。
C d=垂直于构件轴线的阻力系数。
按照规范规定取值为1.0。
C m =惯性力系数。
按照规范规定取值为2.0 。
D =构件的直径。
D =2.3m 。
U =垂直于构件轴线的水质点相对于构件的总速度分量,m/s 。
u =垂直于构件轴线的波浪引起的水质点相对于构件的速度分量,m/s 。
v =垂直于构件轴线的海流引起的水质点速度分量。
计算中海流的方 向取和波浪相同的方向,v = 1.5m/s 。
a =垂直于构件轴线的水质点相对于构件的加速度分量,m/s 2。
附录C 波浪计算C.1 波浪要素确定C.1.1 计算风浪的风速、风向、风区长度、风时与水域水深的确定,应符合下列规定:1 风速应采用水面以上10m 高度处的自记10min平均风速。
2 风向宜按水域计算点的主风向及左右22.5°、45°的方位角确定。
3 当计算风向两侧较宽广、水域周界比较规则时,风区长度可采用由计算点逆风向量到对岸的距离;当水域周界不规则、水域中有岛屿时,或在河道的转弯、汊道处,风区长度可采用等效风区长度Fe,Fe可按下式计算确定:式中ri——在主风向两侧各45°范围内,每隔Δα角由计算点引到对岸的射线长度(m);αi——射线ri与主风向上射线r0之间的夹角(度),αi=i×Δα。
计算时可取Δα=7.5°(i=0,±1,±2,…,±6),初步计算也可取Δα=15°(i=0,±1,±2,±3),(图C.1.1)。
图C.1.1 等效风区长度计算4 当风区长度F小于或等于100km 时,可不计入风时的影响。
5 水深可按风区内水域平均深度确定。
当风区内水域的水深变化较小时,水域平均深度可按计算风向的水下地形剖面图确定。
C.1.2 风浪要素可按下列公式计算确定:式中——平均波高(m);——平均波周期(s);V——计算风速(m/s);F——风区长度(m);d——水域的平均水深(m);g——重力加速度(9.81m/s2);tmin——风浪达到稳定状态的最小风时(s)。
C.1.3 不规则波的不同累积频率波高Hp与平均图C.1.1 等效风区长度计算波高之比值Hp/可按表C.1.3-1确定。
表C.1.3.1 不同累积频率波高换算不规则波的波周期可采用平均波周期表示,按平均波周期计算的波长L 可按下式计算,也可直接按表C.1.3-2确定。
表C.1.3.2 波长~周期~水深关系表L=f(T,d)续表 C.1.3.2C.1.4 设计波浪推算应符合下列规定:1 对河、湖堤防,设计波浪要素可采用风速推算的方法,并按本附录第C.1.2条计算确定。
大尺度桩柱波浪力计算一览表大尺度桩柱波浪力1、极端高水位水深d = -(底面高程) + 极端高水位=+=m波长L =ma 、胸墙, 折算半径D 为:4**m D / L=/=查表惯性力系数 C M =Z 1 =+=m H / d=/=查表ηmax / H=ωt=270oZ 2 =d+ηmax-H/2=+*=m ⑴、最大总波浪力 P 、最大总波浪力矩 M max**sh()-sh()ch()-=***KN·M 29.610.13 6.96.90.3481=8914.13045.7126.929.61886.9025.0 4.6125.00.0015.052π880.17129.610.23320.6325.002.1141=1.9510.25sh=sh-88.088.0= 4.3586 2.89594.202= 2.1782 1.785126.96.9*ch2π*29.6188.02π*25.00=1.954π88.0*0.0727=0.348115.052 1.95mKN30.51=(14.0012.71)1/2=2π(-)ch()=P max =P Imax =M max =M Imax=群墩系数 K x 和附着生物修正系数 nP max ==**=M max ==**=⑵、考虑波浪浮托力作用由于两者结合紧密,故不考虑浮托力作用。
b 、整体, 折算半径D 为:4**m D / L=/=查表惯性力系数 C M =Z 1 =m⑴、计算沉箱部分的最大总波浪力 P max 、最大总波浪力矩 M max取 Z =m**sh()-sh()ch()-=0.0025.000.6892ch2π**sh88.0-1.785 2.1141= 2.895904.20288.0=KN0.6892==1.9610.25126.926.9=0=12.1141=ch2π*29.6188.014.47688 1.968914.1KNKN·M 3045.7KN KN·M11.158914.12π*2.178288.05.51-(ch = 4.20161[-(ch(0.0727(4.4718-0.3932* 4.3586-)]))-ch-ch( 1.785)88.0) 2.178288.088.0)]25.0030.5130.512π*=(12.9512.712π*88.0sh([)1/2=14.476m3.0637]1.153045.73502.60.164529.6188.00.00102516061sh2π*1KnP max KnM max π***2π(-)ch()=P max =P Imax =M max =M Imax=群墩系数 K x 和附着生物修正系数 nP max ==**=M max ==**=⑵、整体的最大总波浪力 P max 、最大总波浪力矩 M maxP max =沉箱Pmax + 胸墙 Pmax =+=M max =沉箱Mmax + 胸墙换算 Mmax=+ (/+)*=⑶、考虑波浪浮托力作用浮托力浮托力矩ωt=o z =mD / L=/=f 0 =f 1 =f 2 =f 3 =**2ch(/)ch(/)*[sin()-cos()]=**3ch(/)ch(/)*[sin()+cos()]27010473KN 104686102513502.625.003502.6104686KN·M)]1.156061(3.0637-1.785=[ 1.785*14.2016=1[2π*25.00sh(=ch2π*29.6188.0-(88.02π*KnM max KnP max 1125.00-)6970.12.8959KN·M -ch2π*3502.6910316970.1KN )]88.02025011.15KN·M9103188.02.114188.0=1.96126.9ch ]0.0088.025.002π*0.7391))1(ch( 1.785)-ch(-KN·M6.988.0*0.73914π=91031270*[1.32270-0.25](0.000)-1.32-0.000.3614.47629.61-0.251.45880.164514.4761.32KN P u =10.25 6.946061880.0088*2π*2π*KN 0.36=3705.4*14.2016-116414.476*1.45270M u =10.25 6.927088322π*29.61882π*0.00=P max =M max =2、设计高水位水深d = -(底面高程) + 设计高水位=+=m波长L =ma 、胸墙, 折算半径D 为:4**m D / L=/=查表惯性力系数 C M =Z 1 =+=m H / d=/=查表ηmax / H=ωt=270oZ 2 =d+ηmax-H/2=+*=m ⑴、最大总波浪力 P max 、最大总波浪力矩 M max**sh()-sh()ch()-=***2π(-)ch()π15.052 1.94850.17710.640.14 6.929.488514.00+KN-574.5KN·M *[(0.000)m15.052]-0.36)1/2KN·M12.71=6705.2*14.201625.00.0025.0025.0 3.5128.51=1.9410.25126.9==(6.9028.510.242=2641.428.51KN26.9*0.3035-sh2π*=sh1.8482.107525.0085.085.0ch2π*28.5185.0= 4.3615 3.09484.174= 2.1789=1.94126.94π0.3035=6225.3KN·M 6.985.0*0.0529=ch2π*28.5185.029.4885.085.0)]85.085.029.48-ch2π*2π*=1[2π*2.107525.00-(ch 4.48sh( 2.1789)85.0=P max =P Imax =M max =M Imax=群墩系数 K x 和附着生物修正系数 nP max ==**=M max ==**=⑵、考虑波浪浮托力作用由于两者结合紧密,故不考虑浮托力作用。
附录C 波浪计算C.1 波浪要素确定C.1.1 计算风浪的风速、风向、风区长度、风时与水域水深的确定,应符合下列规定:1 风速应采用水面以上10m 高度处的自记10min平均风速。
2 风向宜按水域计算点的主风向及左右22.5°、45°的方位角确定。
3 当计算风向两侧较宽广、水域周界比较规则时,风区长度可采用由计算点逆风向量到对岸的距离;当水域周界不规则、水域中有岛屿时,或在河道的转弯、汊道处,风区长度可采用等效风区长度Fe,Fe可按下式计算确定:式中ri——在主风向两侧各45°范围内,每隔Δα角由计算点引到对岸的射线长度(m);αi——射线ri与主风向上射线r0之间的夹角(度),αi=i×Δα。
计算时可取Δα=7.5°(i=0,±1,±2,…,±6),初步计算也可取Δα=15°(i=0,±1,±2,±3),(图C.1.1)。
图C.1.1 等效风区长度计算4 当风区长度F小于或等于100km 时,可不计入风时的影响。
5 水深可按风区内水域平均深度确定。
当风区内水域的水深变化较小时,水域平均深度可按计算风向的水下地形剖面图确定。
C.1.2 风浪要素可按下列公式计算确定:式中——平均波高(m);——平均波周期(s);V——计算风速(m/s);F——风区长度(m);d——水域的平均水深(m);g——重力加速度(9.81m/s2);tmin——风浪达到稳定状态的最小风时(s)。
C.1.3 不规则波的不同累积频率波高Hp与平均图C.1.1 等效风区长度计算波高之比值Hp/可按表C.1.3-1确定。
表C.1.3.1 不同累积频率波高换算不规则波的波周期可采用平均波周期表示,按平均波周期计算的波长L 可按下式计算,也可直接按表C.1.3-2确定。
表C.1.3.2 波长~周期~水深关系表L=f(T,d)续表 C.1.3.2C.1.4 设计波浪推算应符合下列规定:1 对河、湖堤防,设计波浪要素可采用风速推算的方法,并按本附录第C.1.2条计算确定。
附录C 波浪计算时间:2007-01-26 来源:作者:C.1波浪要素确定C.1.1计算风浪的风速、风向、风区长度、风时与水域水深的确定,应符合下列规定:1风速应采用水面以上10m高度处的自记10m i n平均风速。
2风向宜按水域计算点的主风向及左右22.5°、45°的方位角确定。
3当计算风向两侧较宽广、水域周界比较规则时,风区长度可采用由计算点逆风向量到对岸的距离;当水域周界不规则、水域中有岛屿时,或在河道的转弯、汊道处,风区长度可采用等效风区长度F e,F e可按下式计算确定:式中r i——在主风向两侧各45°范围内,每隔Δα角由计算点引到对岸的射线长度(m);αi——射线r i与主风向上射线r0之间的夹角(度),αi=i×Δα。
计算时可取Δα=7.5°(i=0,±1,±2,…,±6),初步计算也可取Δα=15°(i=0,±1,±2,±3),(图 C.1.1)。
图 C.1.1等效风区长度计算4当风区长度F小于或等于100k m时,可不计入风时的影响。
5水深可按风区内水域平均深度确定。
当风区内水域的水深变化较小时,水域平均深度可按计算风向的水下地形剖面图确定。
C.1.2风浪要素可按下列公式计算确定:式中——平均波高(m);——平均波周期(s);V——计算风速(m/s);F——风区长度(m);d——水域的平均水深(m);g——重力加速度(9.81m/s2);t m i n——风浪达到稳定状态的最小风时(s)。
C.1.3不规则波的不同累积频率波高H p与平均图 C.1.1等效风区长度计算波高之比值H p/可按表 C.1.3-1确定。
表 C.1.3.1不同累积频率波高换算不规则波的波周期可采用平均波周期表示,按平均波周期计算的波长L可按下式计算,也可直接按表 C.1.3-2确定。
表 C.1.3.2波长~周期~水深关系表L=f(T,d)续表 C.1.3.2C.1.4设计波浪推算应符合下列规定:1对河、湖堤防,设计波浪要素可采用风速推算的方法,并按本附录第 C.1.2条计算确定。