波浪荷载计算
- 格式:doc
- 大小:251.67 KB
- 文档页数:15
波浪力的计算需要两方面理论的支持:波浪运动理论及波浪荷载计算理论。
前者研究波浪的运动,后者在已知波浪运动的前提下计算波浪对水中物体的作用。
几种常用的波浪普: 1.P-M 谱Pierson 和Moskowitz适用于无限风速发在的波浪普。
国际船模水池会议(ITTC)推荐采用这一形式的波,故也称为ITTC波谱。
JONSWAP(Joint north sea wave project).是一种频谱。
3.应力范围的长期分布模型:1.离散型模型,2.分段连续型模型,3.连续模型。
1. 离散模型:用Hs作为波高,Tz为波浪周期,定义一个余弦波。
然后用规则波理论计算作用在结构上的波浪力。
并用准静定的方法计算结构呢I的应力。
缺陷:没有将波浪作为一个随机过程来处理。
每一海况的应力范围只有一个确的数值。
因此又称为确定性模型。
2.分段连续型模型每一短期海况中,交变应力过程是一个均值为0的平稳正态过程。
综合所有海况中应力范围的短期分布,并得出各个海况出现的疲劳,就得到应力范围的长期分布,它的形式是分段连续的。
应力范围的两种短期分布模型:1.Rayleigh分布和Rice分布。
在某一海况中交变应力均值为。
应力峰值服从Rayleigh分布。
通过计算得出应力范围也服从Rayleigh分布。
3.在船舶及海洋工程结构疲劳可靠性分析中,希望应力范围的长期分布能用一个连续的分布函数来描述。
这就是应力范围长期分布的连续模型.最常用的就是Weibull分布。
4.有义波高:(significant wave height)所有波浪中波高最大的三分之一波浪的平均高度。
用Hs表示。
5.Stokes五阶波给出了波陡的量度(H/L)H/L越大,波就越陡。
当波高与波长的比值大到一定程度时,波会破碎。
6.波速=波长与频率的乘积 C=λ/T或者C=λf,其中f是频率。
或者T=2π/ω7.圆频率1.圆频率即2π秒内振动的次数,又叫角频率,和角速度的ω没有任何关系。
起重船波浪载荷直接计算方法研究程正华【期刊名称】《《河南科技》》【年(卷),期】2019(000)017【总页数】3页(P97-99)【关键词】起重船; 波浪载荷; 直接计算; 单元应力响应函数【作者】程正华【作者单位】招商局重工(江苏)有限公司江苏海门 226100【正文语种】中文【中图分类】TH218起重船又被称为浮吊,翻译自英文crane ship或crane vessel,源于内河和港区水上作业,是为了在水上起吊重物而产生的一种工程船,起重能力跨度从几百吨到上万吨不等。
随着海洋工程的发展,近些年,起重船的作业区域逐渐由港口区域延伸到了近海(Offshore)区域[1]。
波浪荷载,通常也称为波浪力,是波浪对海洋中的结构物所产生的作用,也是海洋工程结构在海洋风浪流环境下地受载的总称。
波浪是一种随机性运动,很难在数学上精确描述,常用特征波法和谱分析法确定[2]。
波浪载荷计算对于起重船的性能有重要作用。
本文根据CCS发布的《钢质海船入级规范》(2015)(以下简称《钢规》)第二分册2.2.1.2中所指定的范围[3],对于L/B≤5和B/D≥2.5的船型,简述一种波浪载荷直接计算方法。
1 坐标系定义坐标系通常选用右手笛卡尔坐标系,X坐标沿船体纵向,向船艏方向为正;Y坐标沿船体横向,向左舷为正;Z轴沿船体垂向,向上为正。
坐标原点在Fr0中纵剖面基线处。
浪向角按逆时针方向定义,随浪(波浪从船艉向船艏传播)为0°,横浪(波浪从右舷向左舷传播)为90°,迎浪(波浪从船艏向船艉传播)为180°。
在水动力分析中,船体作为一个刚体,在波浪中的6个自由度的运动定义如下,沿X 方向的运动为纵荡(Surge),向船艏为正;沿Y方向的运动为横荡(Sway),向左舷为正;沿Z方向的运动为垂荡,向上为正;沿X、Y、Z三轴的旋转分别为横摇(Roll)、纵摇(Pitch)、艏摇(Yaw)3个旋转运动的正向,均按右手法则定义,如图1所示。
游艇码头结构荷载设计樊亮亮;吴国鑫【摘要】针对中国的游艇码头规范中关于结构荷载设计的规定不很详尽的问题, 结合国内外游艇码头设计规范, 对游艇码头的风荷载、波浪荷载、水流力、系揽力以及撞击速度进行研究.通过研究可知, 中国规范中对风荷载、波浪荷载以及水流力的计算较为合理; 而对于船舶系缆力的计算方法略欠妥, 建议参考美国或者澳大利亚规范进行设计计算; 对于撞击工况的计算, 建议参考澳大利亚规范中的撞击速度进行设计计算.【期刊名称】《水运工程》【年(卷),期】2018(000)011【总页数】4页(P64-67)【关键词】游艇码头;荷载;设计【作者】樊亮亮;吴国鑫【作者单位】中交第四航务工程勘察设计院有限公司, 广东广州510230;中交第四航务工程勘察设计院有限公司, 广东广州510230【正文语种】中文【中图分类】U656.1浮式游艇码头或者浮式防波堤兼顾单侧靠船时一般采用浮筒或者浮箱结构。
对于该结构及船舶所受荷载类型,各国规范的计算方法略有不同。
浮式结构设计荷载大致可分为风荷载、波浪荷载、水流力、船舶荷载、人群荷载以及流动机械荷载等。
对于整体结构的受力分析,只需要考虑风、浪、流和船舶荷载即可,竖向荷载只针对浮筒自身的稳定性分析。
浮筒自身的稳定性一般由浮筒自身的产品性能选定,因此本文只对结构整体的设计荷载进行计算分析。
本文基于中国规范[1]、美国规范[2]以及澳大利亚规范[3]等,通过对比各国游艇码头结构荷载计算,进一步完善我国游艇码头的结构设计。
1 荷载计算1.1 风荷载计算1.1.1 中国规范中国游艇码头的风荷载计算与常规码头的相同,计算公式按照《港口工程荷载规范》[4]中的E.0.1条进行计算。
(1)(2)式中:Fxw、Fyw分别为作用在船舶上的计算风压力的横向和纵向分力(kN);Axw、Ayw分别为船体水面以上横向和纵向受风面积(m2);vx、vy分别设计风速的横向和纵向分量(ms);ζ1为风压不均匀折减系数;ζ2为风压高度变化修正系数。
高 新 技 术0 引言由于大量海生物的出现,核电站经常发生海生物和漂浮物进入冷源入水口的事件,从而造成核电厂冷源系统堵塞[1-2]。
目前,核电站采用传统“网兜拦截+人工海上渔船驳运清理”的方式来防护冷源系统,这种方式效率低且夜间作业风险高,会直接对核电站安全运行造成影响[3]。
因此,辽宁某公司和大连某公司联合设计、制造了一种冷源取水口拦污网网兜海生物抽吸装置,将其安装在工作船上,将被拦截入网兜的海生物和漂浮物输送到海面,并运送至指定地点,代替传统人工海上渔船驳运清理冷源入水口拦截网兜中海生物的作业方式,该装置可以在设计工况下全天候、无间断地开展清理作业,彻底解决人工海上渔船驳运清理存在的安全隐患,提高经济效益,保障冷源安全。
冷源取水口拦污网网兜抽吸装置工作船在进行工作时,启动锚机系统,抛锚索固定,将海生物切割装置调整至工作状态,连接网兜,即可在操作台上按要求模式进行抽吸作业。
因此,需要根据具体工作的海洋环境条件对工作船的运动响应进行分析,获取关键位置的运动响应结果,评估其抽吸工作的安全性能。
1 工作船设计需求1.1 工作海况条件根据船体的对称性,考虑风、浪和流环境荷载入射方向与船长方向的夹角为0°~180°(风、浪和流共向),如图1所示。
选取2级海况波浪,3级风风速,具体的环境参数值见表1。
表1 作业工况环境参数环境参数数值不规则波有义波高/m0.5谱峰周期/s8均匀风风速/(m/s)4均匀流流速/(m/s)0.31.2 运动响应限制工作船结构示意图如图2所示。
工作船能够承载提升装置、输送泵系统、满负荷为1 m 3的过滤箱、锚机以及就地控制柜等设备和4名操作人员,工作台设计有4个锚绞机和导缆孔,以保证工作船水下连接网兜位置处的耦合装置的漂移范围不超过1.5 m。
2 运动响应分析模型根据工作船的运行情况,建立其对应的简化运动响应分析模型,如图3所示。
工作船会受到风、浪以及流环境的联合作用。
§ 2-1水工建筑物的荷载计算水工建筑物上的作用有:垂力、水作用、渗透作用力、风及波浪作用、冰及冰冻作用、温度、土及泥沙作用、地震作用等。
一、自重W=Vy一般素碗取23. 5'24kN/m‘,钢筋碗取24. 5"25kN/m3,浆砌石取21. 5"23kN/m3,对土石坝的材料重度应根据具体性能及不同部位,分别取湿重度、干重度、饱和重度、浮重度等几种情况计算。
水工建筑物上永久固定设备,如闸门、启闭机等,其口重标准值采用设备标牌重量作用分项系数:大体积混凝土、土石坝取1.0;对普通水工混凝土、金属结构(设备)取1. 05,当口重对结构有利时取0. 950地下工程的混凝土衬砌取1. 1,其对结构有利时取0.9。
二、水压力水体对各种水工结构均发生作用,作用结果是对结构产生水圧力,其可分为静水压力和动水压力。
1.静水压力水体静止状态下对某结构表面的作用力称为静水圧力(1)作用在坝、闸等结构面上的水压力r SI P、二|rJ(2)管道及地下结构上的水压力计算。
内水压力:作用在管道内壁上的静水圧力;外水压力:作用于管道或衬砌外侧的水斥力。
对内水圧力,为计算方便,常将其分解成均匀内水压力和非均匀内水压力两部分。
对有床隧洞的碗衬砌的外水床强标准值可按式(2-6)计算。
(2-6)式中:3 一一作用于衬砌上的外水压强标准值(KN/m:);a一一外水压力折减系数,可按表2-1采用:a 一一作用水头(m),按设计采用的地下水位线与隧洞中心线的高差确定。
同内水压力一样,外水压力也可分解成均匀外水压力和非均匀外水圧力。
非均匀外水压力的合力方向垂直向上,合力的大小应等于单位洞长排开水体的重2.动水压力(1) 渐变流时的时均压强:式中:H 一一过流面上计算点的时均圧强代表值(N/m‘);3 -一水的密度(kg/m 3):g 一一重力加速度(m/s 3);h 一一计算点A 的水深(m):0——结构物底面与平面的夹角。
水工建筑物的荷载计算水工建筑物上的作用有:重力、水作用、渗透作用力、风及波浪作用、冰及冰冻作用、温度、土及泥沙作用、地震作用等。 一、自重 W=Vγ 一般素砼取23.5~24kN/m3,钢筋砼取24.5~25kN/m3,浆砌石取21.5~23kN/m3,对土石坝的材料重度应根据具体性能及不同部位,分别取湿重度、干重度、饱和重度、浮重度等几种情况计算。 水工建筑物上永久固定设备,如闸门、启闭机等,其自重标准值采用设备标牌重量 作用分项系数:大体积混凝土、土石坝取1.0;对普通水工混凝土、金属结构(设备)取1.05,当自重对结构有利时取0.95。地下工程的混凝土衬砌取1.1,其对结构有利时取0.9。 二、水压力 水体对各种水工结构均发生作用,作用结果是对结构产生水压力,其可分为静水压力和动水压力。 1.静水压力 水体静止状态下对某结构表面的作用力称为静水压力 (1)作用在坝、闸等结构面上的水压力
PH=221Hw
PV=wwV (2)管道及地下结构上的水压力计算。 内水压力:作用在管道内壁上的静水压力; 外水压力:作用于管道或衬砌外侧的水压力。 对内水压力,为计算方便,常将其分解成均匀内水压力和非均匀内水压力两部分。
hpwwr'
)cos1(''iwwrrp
对有压隧洞的砼衬砌的外水压强标准值可按式(2-6)计算。 eeekHp (2-6)
式中:ekp——作用于衬砌上的外水压强标准值(KN/m2); e——外水压力折减系数,可按表2-1采用;
eH——作用水头(m),按设计采用的地下水位线与隧洞中心线的
高差确定。 同内水压力一样,外水压力也可分解成均匀外水压力和非均匀外水压力。 非均匀外水压力的合力方向垂直向上,合力的大小应等于单位洞长排开水体的重量。
2.动水压力 (1)渐变流时的时均压强: cosghpwtr
式中:trp——过流面上计算点的时均压强代表值(N/m2); w——水的密度(kg/m3);
1ˎ某水库库区多年平均最大风速为15m/s, 吹程2000m, 水库平均水深30m,水库正常蓄水位为1426m,设计洪水位为1427.7m, 校核洪水位为1428.8m, 工程等级为4等, 上游坝坡为混凝土护坡, 坡度为1:2。
气象站位于河岸广阔的平原上,距建筑物及树木100米以外。
计算浪压力及坝顶超高。
解析:由题意可知:1,坝顶超高的计算当浪压力参与荷载的基本组合时,计算风速采用重现期为50年的最大风速。
当风速测量高度m h 00.10c =时,依据《碾压式混凝土设计规范》查得:风速高度修正系数00.1z =k ,则计算风速s m 15.0015.001.00v 0=⨯=按《碾压式混凝土设计规范》可知平原地区按莆田实验公式确定:⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛v h g v gD v h g v h m th th m th m 07.013.0020018.007.013.0g 7.045.07.00 []⎭⎬⎫⎩⎨⎧⨯⨯=0.097.470018.07.013.0g 1.0380.70th th v h m=0.09×0.15=0.0135 m h 31.081.90135.0152m =⨯= 对于丘陵,平原地区水库,当26.5m v 0≤,7500m ≤D 时可用鹤地公式计算: 04.037.400625.002000625.00g 156131612%2=⨯⨯==⎪⎪⎭⎫ ⎝⎛v gD v v h 0.989.810.0415h 22%=⨯=36.034.90386.0020386.00g 212=⨯=⨯=⎪⎪⎭⎫ ⎝⎛v gD v l m 8.279.810.3615L 2m =⨯= 风壅水水面高度:m g D k H v h m1015103262z 45.289.03081.9220006.3cos 20--⨯=⨯⨯⨯⨯⨯⨯==α 由规范可知: 23.242.2h 2%1%=h ,则m 2.1h 1%= 又因为该混凝土坝体的工程等级为4等,则由规范可知m h 4.0c = 所以坝顶超高h h h c z H ++=∆%1=1.2+0.00245+0.4=1.6m2ˎ浪压力计算 因为2m L H ≥时,波浪运动不受库底约束,这样的条件下的波浪成为深水波,则由题意和图形可知三角形分布的面积之差纪委浪压力,则浪压力为:()m h h L KN P z 2%1m 39.242.128.204=⨯=+=γ。
关于水库风浪高度计算公式的几个问题1 关于鹤地水库波高计算公式现行规范《水工建筑物荷载设计规范》(SL744-2016),《水工建筑物荷载设计规范》(DL5077-1997),推荐的波浪要素计算方法分别是莆田试验站公式、鹤地水库公式及官厅水库公式。
莆田试验站公式如下。
式中-平均波高,m;-平均波周期,s;-计算风速,m/s;D-风区长度,m;-水域平均水深,m;g-重力加速度,取9.81m/s2。
注意。
在SL744-2016与DL5077-1997中,平均波长计算公式中的水深符号与平均波高计算公式中的水深符号不同,但没有相应说明内容。
在《碾压式土石坝设计规范》(SL274-2001)、《碾压式土石坝设计规范》(DLT5395-2007)及《小型水利水电工程碾压式土石坝设计规范》(SL189-2013)中,平均波长计算公式中水深采用坝迎水面前水深。
在《碾压式土石坝设计规范》(SDJ218-84,作废)、《堤防工程设计规范》(GB50286-2013)、《滩涂治理工程技术规范》(SL389-2008)、《海堤工程设计规范》(SL435-2008)中及《广东省海堤工程设计导则(试行)》(DB44/T182-2004),平均波长计算公式中水深采用水域平均水深。
《小型水利水电工程碾压式土石坝设计导则》(SL189-1996,作废)第6.1.7条,波高可采用莆田试验站公式或官厅-鹤地公式等算出。
官厅-鹤地公式是指波高按官厅水库公式计算、波长按鹤地公式计算。
新版SL189-2013仅推荐采用莆田试验站公式计算波浪要素,但没有说明原因。
在SDJ218-84附录一中,推荐的波浪要素计算公式有莆田试验站公式、安德烈扬诺夫公式及官厅-鹤地公式。
安德烈扬诺夫公式如下。
,并在注中说明,原公式作者未规定计算波高的累积频率,经比较,当时,可取为,当时,可取为。
官厅水库波高、波长计算公式如下。
并说明,经比较,官厅水库波高公式的波高累积概率当时,可取为,当时,可取为,现规范中为。
整理后:波浪荷载的计算理论波浪是发生在海洋表面的一种波动现象,其波动性质因受浅水区域海底地形影响和水深的变浅,发生波浪破碎现象,成为影响海岸侵蚀和变形以及海岸带污染物迁移与扩散的最主要的水动力环境之一。
破浪破碎与冲击现象对海上工程设施的安全也十分重要。
由于波浪破碎及冲击作用的机理极其复杂,至今仍然是海岸工程领域没有解决的困难课题之一。
因此,开展近海波浪破碎与冲击过程数值模型的研究,就有着重要的理论意义和工程意义。
波浪荷载,也称波浪力,是波浪对港口码头和海洋平台等结构所产生的作用。
目前按绕射理论进行分析。
波浪对结构物的作用由四部分组成:水流粘性所引起的摩阻力(与水质点速度平方成正比);不恒定水流的惯性或结构物在水流中作变速运动所产生的附加质量力(与波浪中水质点加速度成正比);结构物的存在对入射波浪流动场的辐射作用所产生的压力和结构物运动对入射波浪流动场的辐射作用所引起的压力。
包括上述全部作用影响的波浪力理论称为绕射理论。
在目前实际工作中,常用只考虑了结构受到波浪摩阻力和质量力影响的半经验半理论的莫里森(Mrison)方程分析波浪力。
波浪荷载是由波浪水质点与结构间的相对运动所引起的。
波浪是一随机性运动,很难在数学上精确描述。
当结构构件(部件)的直径小于波长的20%时,波浪荷载的计算通常用半经验半理论的美国莫里森方程;大于波长的20%时,应考虑结构对入射波场的影响,考虑入射波的绕射,计算时用绕射理论求解。
影响波浪荷载大小的因素很多,如波高、波浪周期、水深、结构尺寸和形状、群桩的相互干扰和遮蔽作用以及海生物附着等。
波浪荷载常用特征波法和谱分析法确定。
对一些特殊形状或特别重要的海洋工程结构,除了用上述的方法进行计算分析外,还应进行物理模型试验,以确定波浪力。
①特征波法。
选用某一特征波作为单一的规则波,并以它的参数(有效波高、波浪周期、水深)和结构的有关尺寸代入莫里森方程或绕射理论的公式,求出作用在结构上的波浪力。
此法简便易行,在海洋工程设计广泛应用。
②谱分析法。
利用海浪谱进行波浪荷载计算、结构疲劳和动力响应分析的一种方法。
把波浪作为随机性的、由许多不同波高和波周期的规则波线性迭加而成的不规则波,用概率论和数理统计的方法收集、分析处理波浪观测数据,由于它能较精确地反映波浪的能量分布规律,所以是一种比较理想的方法。
海洋工程结构设计中常用的有P-M和联合(JONSWAP)谱。
波力谱确定后,可求出波浪力分布函数中的统计特征值,进而得到某一累积概率的波浪力。
由于波浪具有明显的随机性,难用确定的函数表达,故在波浪的研究中常采用多个或无限个振幅、频率、方向、位相不同的简单波的叠加,并规定组成波的振幅或相位是随机量,从而叠加的结果为随机函数,以反映波浪的随机性。
实践证明这种方法是可行的,它以成为研究波浪要素的统计特性的分布来描述它,另一是用波浪要素的“谱”来表征其内部的频率结构。
当然,波浪外观上表现出来的性质和它的内部结构是有联系的。
水库波浪在风里直接作用下产生的运动,表面十分复杂,在统计过程中是把波浪当作准稳定的随即过程来处理,每次测量时间为10-5小时。
如果观测是段短、波数少时,为了提高精度,可父子俩偏差大的缺点,也可将各组中每种波的出现概率进行加权统计。
左图是根据1966、1967年密云水库和还有那个水库的原形观测资料绘出的波高小于H的概率曲线,它与三元海浪概率分布的克雷洛夫共识甚为符合(即图上所示公式)。
波浪荷载作用下土体的动力特性:在波浪荷载作用下,海床中的土单元也受到一系列循环荷载作用. 在某一时刻,当波峰作用在所研究的土单元正上方,则会产生正的竖向压力;当波谷作用在其上时,则产生负的竖向压力. 这样,在一个波长距离的波作用下,产生的应力是由三轴应力作用的圆形轨迹.在波高为零的瞬时,波作用在土单元上,产生水平剪应力并引起单剪模式的应变,这一水平剪应力分量也随着波的传播而改变其方向,引起剪应力的另一类型的循环交替. 值得注意的是,上述两种循环剪应力是交替作用而不是同时作用的. 三轴试验的剪切模式所产生的循环应力与单剪模式的水平剪应力相位差为90°. 因此,由波浪荷载产生的作用在海床土体的循环应力是沿主应力方向连续旋转的,其应力交替的性状可由图2 (b) 中的τvh与(σv - σh) / 2 之间的圆形关系来表示.实例分析:直立浮式截圆柱柱群的绕射问题假设流体为不可压、无粘性、均匀的理想流体,流场中运动处处无旋。
对静止于水深为d的水域中的N个相同吃水h、半径a的柱体,坐标系统如图1,单个柱体的圆心Oj 坐标为(xoj,yoj,z)(j=1,…,N),建立局部的柱坐标系(γj ,θj,z)。
这样第k个柱体中心Ok相对于j柱局部坐标系有(Rjk ,θjk,z)极坐标,(j,k=1,2,…,N)。
入射波采用线性微幅波理论。
1 坐标系统为了将每个单元柱体的不同散射波成分迭加并计入柱间水动力干扰,引入大间距假设,即认为柱间距Rjk,j,k=1,2,……,N与入射波波长相比足够大,k0Rjk1。
这样由柱群中任一圆柱在入射波作用下产生的绕射波对其它圆柱的作用可近似为非平面修正的等效平面波,即改进平面波法。
对每一个柱体而言,其速度势中的未知系数与其编号无关,即决定未知系数的方程对每一个柱体都是相同的。
因此,柱群情况下只需考虑第j柱附近的速度势即可。
对于第j柱,传播方向与x轴正向夹角为β的线性规则波速度势Φj2I 可记为:(1)式中,A为入射波波幅(m),ω为入射波频率(s-1);Jm(x)为第一类m阶Bessel函数,εm 为Neumann常数,且ε=1,εm=2(m≥1);k为波数,应满足色散关系:ω2=gk0thkd。
PHj =exp{i(kxjcos β+kyjsin β)}为相位项。
对应第j柱的内、外域速度势应有:(2)(3)由入射波速度势Φj2I 引起的第j柱的绕射速度势Φj2s为:(4)式中,考虑另一个柱体k,半径亦为a,到j柱的距离为Rjk。
由第j柱入射波引起的第k柱绕射势为:(5)式中:代表了等效平面波幅:(6)为非平面波修正项,其中:这样对N(N>2)个柱体所组成的群柱中,外部入射波在第j柱的入射势仍如(1)式。
与上述二柱情形类似可得,由其他N-1个柱的绕射波在j柱:产生的等效平面入射波势Φj2I1(7)对于j柱,等效平面入射波的一阶非平面修正项包含两部分:对其它N-1个柱体外部入射波的绕射波的修正;其它柱体的绕射波在j柱的反射后的等效平面波修正。
(8)由于等效平面波幅C包含两部分:k柱对外部入射波引起的绕射;由其它jk柱绕射波引起在k柱的二次绕射,这样得:(9)式中,j,k=1,2,…,N,j≠k。
由(9)式可确定未知的等效平面波幅C。
jk(i=1,2),可由伯努利方程确定流场中任一点确定了流场速度势Φi压力:(10)这样沿湿表面上积分可得波浪对直立浮式柱群的q模态线性干扰力或力矩:(11)其中:q=1,…,6,分别对应纵荡、横荡、垂荡、横摇、纵摇和首摇等运动模态。
3 数值计算结果分析本文应用上述理论和方法计算了不同波长、波频情况下两根直立浮式圆柱群的波浪荷载,并与现有的文献结果进行了比较,取得了良好的一致。
图中圆柱横截面半径a=10m,柱吃水h=5m,水深d=100m,波幅A=1m,ρ为为波数,ω(s-1)为入射波频率,两柱间距用海水密度,g为重力加速度,kR表示。
图中箭头(→)表示入射波的方向,①表示前柱,②表示后柱。
图中的点号(。
)和叉号(×)分别表示相应状态下由A N Williams &Z Demirbilek计算所得的双柱的波浪力[6]。
图2为单柱及双柱迎浪状态下前柱(柱1)和后柱(柱2)所受x轴方向波浪荷载的幅频变化曲线,此时入射波向角β=0°,h/a=0.5,d/a=10,R=3a。
图3与图2的情况相同,只是R=5a。
图4为单柱及双柱前、后柱所受y轴方向波浪荷载的幅频变化曲线,此时β=0°,h/a=0.5,d/a=10,R=3a。
图5与图4的情况相同,只是R=5a。
图6和图7分别R=3a和R=5a时,单柱及双柱迎浪状态下前柱和后柱所受波浪力矩荷载的幅频变化曲线,此时β=0°,h/a=0.5,d/a=10。
图2 双柱迎浪间距3a时x轴方向受力幅频曲线图3 双柱迎浪间距5a时x轴方向受力幅频曲线图4 双柱迎浪间距3a时z轴方向受力幅频曲线Fig.4 Nondimensional amplitude-frequency curves of load in zaxle of two cylinders R=3a图5 双柱迎浪间距5a时z轴方向受力幅频曲线图6 双柱迎浪间距3a时时纵摇干扰力矩幅频曲线图7 双柱迎浪间距5a时纵摇干扰力矩幅频曲线(1)从图中可见,柱群所受的波浪力对于柱间距和入射波频率十分敏感。
随着柱间距离的增加,前柱1所受荷载峰值明显变小。
如图3,R =5a时,显然柱1在x轴方向的波浪荷载最大幅值比R=3a(图2)时小。
但是它随k0a的变化明显地比R=3a时的复杂:图2中ka在区间(0,3.0)即频率ω在区间]0,1.7[内变化时荷载Fx出现了两个峰值而图3同样区间内Fx却出现了四个峰值。
图5和图7中柱1所受荷载和力矩变化也呈现出这样的趋势。
相比之下,后柱2的变化总是比较平缓,它所受的荷载和力矩并不随着柱间距的增加而有剧烈的变化。
(2)柱间水动力相互作用在入射波频率很低时并不明显。
如图3,当频率ω在区间(0,0.5)即较低频范围内变化时,单柱、柱1和柱2所受荷载幅频曲线基本重合,图7中低频时柱体所受力矩曲线也基本重合。
可见低频时柱间的水动力相互作用并不明显,因此此时计算荷载时将其忽略也是合理的。
(3)不同的波浪特性对柱群所受荷载有不同的影响。
随着入射波频率的变化,前柱1所受力和力矩荷载以单柱荷载为平均值交错变化,在某一频率出现的最大值要超过同频率下单柱所受荷载。
某种参数组合的条件下,浮式柱群所受的荷载比单柱所受荷载要大许多,即使柱间距较大也是如此:如图2中无因次的Fx 大约在ka=0.8即ω=0.89时出现最大值2.2,远远大于同频率下单柱的值1.4,而即使在R=5a的情况下,如图3,柱1的无因次Fx仍在某一频率时达到最大值1.95,大于同频率下单柱的值1.45。
因此,在平台设计中,设计人员只有充分考虑到这种由于柱间的相互作用所产生的荷载增加,才能有效地保证平台的安全。
相比之下,后柱2所受荷载都比同频率下单柱所受荷载小,而且二者的荷载幅频曲线的变化趋势大体相当。
出现这种现象主要是由柱间水动力相互作用引起的。
柱间水动力相互作用可分为干扰效应和遮蔽效应,干扰效应使波浪荷载变化剧烈且使幅值增大,而遮蔽效应也使波浪荷载变化剧烈但却使其幅值减少。