《船舶疲劳ppt课件.ppt》由会员分享,可在线阅读,更多相关《船舶疲劳ppt课件.ppt(217页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第一章 绪论绪论 第一节 什么是疲劳? 1 现象:铁丝反复折 断 2 定义:美国试验与材料协会(ASTM)在“疲劳试验及数据统计分析之有关术语的标准定义”(ASTM E206-72)中规定: 在某点或某些点承受交变应力且在足够多的循环扰动作用之后形成裂纹或完全断裂的材料中发生的局部的永久结构变化的发展过程,称为疲劳。 3 疲劳问题的特点 只有在承受交变应力作用的条件下,疲劳才会发生。 随时间交替变化的应力,随时间交替变化的应力,也可称为扰动载荷。也可称为扰动载荷。 变化可以是有规则的,变化可以是有规则的,也可以是不规则的也可以是不规则的 ,甚至,甚至是随机的是随机的。 疲劳破坏起源于高应力或高
2、应变的,一般是几何形状变化或材料缺陷等引起应力集中的局部细节; 疲劳破坏是在足够多次的交变载荷作用后,形成裂纹或完全断裂的现象。 疲劳是一个发展的过程。 划分成三个阶段: 裂纹萌生(起始) 扩展 断裂 注:在疲劳分析中经常用到的参数 应力范围: 应力幅值: 平均应力: 应力比: 反映循环特性当 时, 对称循环当 时, 脉动循环当 时, 静载荷minmaxSSS2/ )(2/minmaxSSSSaminmaxmaxminSSSSR或minmaxSSminmaxSS2/ )(minmaxSSSm1R0minS0R1R 4 疲劳研究的目的 1)疲劳寿命:从结构开始使用到裂纹萌生,扩展并最后断裂,这个
3、过程所经历的时间或交变载荷作用次数,称为“寿命”。 它取决于载荷水平作用次数或时间及原材料抵抗疲劳破坏的能力。 2)目的:研究寿命预报的方法 Ntotal=Ninitiation+Npropagation Ninitiation:起始寿命,由应力寿命关系,应变寿命关系 Npropagation:扩展寿命,由断裂力学方法第二节 疲劳破坏机理 1 断口的宏观特征 有裂纹源、疲劳裂纹扩展区和最后断裂区三个部分; 裂纹扩展区断面较光滑、平整,通常可见“海滩条带”,有腐蚀痕迹; 裂纹源通常在高应力局部或材料缺陷处; 与静载破坏相比,即使是延性材料也没有的明显的塑性变形; 工程实际中的表面裂纹一般称半椭圆
4、形。 2 疲劳裂纹萌生机理 金属大多是多晶体,各晶粒有各自不同的排列方位。在高应力作用下,材料晶粒中易滑移平面的方位若与最大作用剪应力一致,则将发生滑移。 滑移可以在单调载荷下发生,也可以在循环载荷下发生。 在较大载荷作用下发生粗滑移和在较小的循环载荷作用下发生细滑移。 在循环载荷作用下,材料表面发生滑移带“挤出”和“凹入”,进一步形成应力集中,导致微裂纹产生。应当注意,滑移主要是在晶粒内进行的。深度大于几个微米的少数几条滑移带穿过晶粒,称为“持久滑移带”,微裂纹正是由这些持久滑移带发展而成的。滑移带的发展过程与施加的载荷及循环次数有关,随着循环次数的增加,滑移线越来越密集,越来越粗。 3 疲
5、劳裂纹扩展机理 疲劳裂纹在高应力处由持久滑移带萌生,是由最大剪应力控制的。形成的微裂纹与最大剪应力方向一致。 在循环载荷作用下,由持久滑移带形成的微裂纹沿45度最大剪应力作用面继续扩展或相互连接。此后,有少数几条裂纹达到几十微米的长度,逐步汇聚成一条主裂纹,并由沿最大剪应力面扩展逐步转向沿垂直于载荷作用线的最大拉应力面扩展。裂纹沿45度最大剪应力面的扩展是第1阶段的扩展,在最大拉应力面内的扩展是第2阶段的扩展。 从第1阶段向第2阶段转变所对应的裂纹尺寸主要取决于材料和作用应力的大小,但通常都在0.05mm内, 只有几个晶粒的尺寸。 第1阶段裂纹扩展的 尺寸虽小,对寿命的 贡献却很大,对于高 强
6、度材料,尤其如此。 与第1阶段相比,第2阶段的裂纹扩展较便于观察。Laird直接观察了循环应力作用下延性材料中裂纹尖端几何形状的改变,提出了描述疲劳裂纹扩展的“塑性钝化模型”,如图所示。 (a)给出了循环开始时的裂纹尖端形状; 随着循环应力的增加,裂纹逐步张开,裂尖材料由于应力集中而沿最大剪应力方向滑移(b); 应力进一步增大,裂纹充分张开,裂尖钝化成半圆形,开创出新的表面(c); 卸载时已张开的裂纹要收缩,但新开创的裂纹面却不能消失,将在卸载引入的压应力作用下失稳而在裂纹尖端形成凹槽形(d); 最后,在最大循环压应力作用下,又成为尖裂纹,但其长度已增加了一个小长度(e)。 下一个循环,裂纹又
7、张开、钝化、扩展、锐化,重复上述过程。这样,每一个应力循环,将在裂纹面上留下一条痕迹。 第三节 疲劳断裂研究方法 疲劳断裂问题需要研究载荷谱、裂纹萌生及扩展规律、构件细节应力分析,疲劳寿命预测和抗疲劳设计方法,等等。一方面,由于涉及因素多,问题复杂,难以找到解析的、普遍的寿命预报方法;另一方面,工程应用的需求迫切。 因此,研究问题时必须抓住主要因素,建立简化模型,逐步深化认识。 例如,对于载荷谱,先研究恒幅循环载荷的最简单情况,再考虑变幅载荷下的累积损伤,最后考虑随机载荷。对于裂纹萌生及扩展规律,则先研究不含有缺陷的光滑材料在恒幅循环载荷作用下的裂纹萌生规律,给出应力-寿命、应变-寿命以及不引
8、发裂纹的疲劳极限等的基本关系,在讨论应用于构件时所需要进行的必要的修正,建立裂纹萌生寿命估算方法,满足无限寿命设计、安全寿命设计的要求。再讨论含有裂纹材料的断裂和疲劳裂纹扩展规律,研究断裂判据,研究在不同载荷谱作用下裂纹扩展寿命的预测,建立损伤容限设计方法。 关于寿命预测和抗疲劳设计方法,是依赖于对问题的认识水平,从不考虑裂纹向考虑裂纹;从确定性分析向可靠性分析,从控制构件和结构的安全向综合控制设计-制造-使用-维修,以安全和经济为目标,逐步发展、丰富的。 此外,还应研究疲劳破坏的基本机理,不断积累、深化对于疲劳断裂破坏的更本质的认识,不断提高抗疲劳设计能力。疲劳断裂研究的基本思路如下,第二章
9、第二章 疲劳强度模型疲劳强度模型S-N曲线曲线 1、S-N曲线 材料的疲劳性能用作用的应力范围S与到破坏时的寿命N之间的关系描述,即S-N曲线。 寿命N定义为在给定应力比R下,恒幅载荷作用下循环到破坏的循环次数。问题:如何得到S-N曲线? 实验得到! 疲劳破坏有裂纹萌生,扩展至断裂三个阶段,这里破坏指的是裂纹萌生寿命。因此,破坏可以定义为: 1)标准小尺寸试件断裂。对于高、中强度钢等脆性材料,从裂纹萌生到扩展至小尺寸圆截面试件断裂的时间很短,对整个寿命的影响很小,考虑到裂纹萌生时尺度小,观察困难,故这样定义是合理的。 2)出现可见小裂纹,或有515应变降。对于延性较好的材料,裂纹萌生后有相当长
10、的一段扩展阶段,不应当计入裂纹萌生寿命。小尺寸裂纹观察困难时,可以监测恒幅循环应力作用下的应变变化。当试件出现裂纹后,刚度改变,应变也随之变化,故可用应变变化量来确定是否萌生了裂纹。 材料疲劳性能试验所用标准试件,(通常为710件),在给定的应力比R下,施加不同的应力范围S,进行疲劳试验,记录相应的寿命N,即可得到图示S-N曲线。NS 由图可知,在给定的应力比下,应力范围S越小,寿命越长。当应力范围S小于某极限值时,试件不发生破坏,寿命趋于无限长。 由S-N曲线确定的,对应于寿命N的应力范围 ,称为寿命为N循环的疲劳强度。寿命N趋于无穷大时所对应的应力范围S,称为材料的疲劳极限。 由于疲劳极限
11、是由试验确定的,试验又不可能一直做下去,故在许多试验研究的基础上,所谓的无穷大一般被定义为: 钢材,107次循环,焊接件:2*106。2、S-N曲线的数学表达式 NSm=A两边取对数, LogN +mLogS=LogA 选取几个不同的应力范围平 , ,进行n组疲劳试验,对各组实验数据1S2SnS1SS14S13S12S11N,N,N,NS1iN2SS24S23S22S21N,N,N,NS2iN3SS34S33S32S31N,N,N,NS3iN应力范围循环次数两个参数:m,A假定 为某一概率分布 (一般为Weibull分布)存活率 则可求得存活率为p的,分别对应于 , , 的试验次数多 少S1i
12、N, S2iN Nf pdNNfNp1S2SnSpnp3p2p1N.N,N,NpnnP33p22p11N,S,.,S,N,N,S,N,S 假定应力范围水平下疲劳寿命N的分布为对数正态分布时,采用极大似然法拟合得到P-S-N曲线为 其中m定值, 表示存活率为p时的 正态分布 标准差 个 mlgSlgAlgNpplgAplgA mN,S,.,S,N,N,S,N,Snn332211lgAlgAppulgAlgAn,lgAm,N,Siiin1iilgAn1lgAn1i22ilgAlgAnlgA1-n1对于船海工程,一般构件 )0 . 2u(72.97pp 00mlgS2lgAlgNlgA主要构件 )0
13、 . 3u(99.87pp 00mlgS3lgAlgNlgA 在实际设计或计算中,为了得到适合的S-N曲线,需要做实验吗? 可以查阅相关规范或资料,得到S-N曲线F2F2FF2总结: S-N曲线表征结构的抗疲劳能力,由实验得到。 实验中根据结构形式和载荷类型选取S-N曲线,此时S-N曲线都是对应于一定的概率水平的!3、平均应力的影响 材料的疲劳性能,用作用应力S与到破坏时的寿命N之间的关系描述。在疲劳载荷作用下,最简单的载荷谱是恒幅循环应力。 R=-1时,对称恒幅循环载荷控制下,试验给出的应力寿命关系,是材料的基本疲劳性能曲线。 本节讨论应力比R变化对疲劳性能的影响。 如图所示,应力比R增大,
14、表示循环平均应力Sm增大。且应力幅Sa给定时有 Sm=(1+R)Sa/(1-R) 一般趋势 当Sa给定时,R增大,平均应力Sm也增大。循环载荷中的拉伸部分增大,这对于疲劳裂纹的萌生和扩展都是不利的,将使得疲劳寿命降低。 平均应力对S-N曲 线影响的一般趋势 如图所示。 平均应力Sm=0时的S-N曲线是基本S-N曲线。当Sm0,即拉伸平均应力作用时,S-N曲线下移,表示同样应力幅作用下的寿命下降,或者说在同样寿命下的疲劳强度降低,对疲劳有不利的影响。SmS(拉)S(扭) 假定作用应力水平相同,拉压时高应力区体积等于试件整个试验段的体积;弯曲情形下的高应力区体积则要小得多。我们知道疲劳破坏主要取决
15、于作用应力的大小(外因)和材料抵抗疲劳破坏的能力(内因)二者,即疲劳破坏通常发生在高应力区或材料缺陷处。假如图中的作用的循环最大应力Smax相等,因为拉压循环时高应力区域的材料体积较大,存在缺陷并由此引发裂纹萌生的可能性也大。 所以,同样的应力水平作用下,拉压循环载荷作用时的寿命比弯曲时短;或者说,同样寿命下,拉压循环时的疲劳强度比弯曲时低。 扭转时疲劳寿命降低,体积的影响不大,需由不同应力状态下的破坏判据解释,在此不作进一步讨论。 2)尺寸效应 不同试件尺寸对疲劳性能的影响,也可以用高应力区体积的不同来解释。应力水平相同时,试件尺寸越大,高应力区域材料体积就越大。疲劳发生在高应力区材料最薄弱
16、处,体积越大,存在缺陷或薄弱处的可能就越大,故大尺寸构件的疲劳抗力低于小尺寸试件。或者说,在给定寿命N下,大尺寸构件的疲劳强度下降;在给定的应力水平下,大尺寸构件的疲劳寿命降低。 3)表面光洁度 由疲劳的局部性显然可知,若试件表面粗糙,将使局部应力集中的程度加大,裂纹萌生寿命缩短。材料的基本S-N曲线是由精磨后光洁度良好的标准试件测得的。 4) 表面处理 一般来说,疲劳裂纹总是起源于表面。为了提高疲劳性能,除前述改善光洁度外,常常采用各种方法在构件的高应力表面引入压缩残余应力,以达到提高疲劳寿命的目的。 若循环应力如图中1-2-3-4所示,平均应力为Sm,则当引入压缩残余应力Sres后,实际循
17、环应力水平是原1-2-3-4各应力与-Sres的叠加,成为1-2-3-4,平均应力降为Sm,疲劳性能将得到改善。 表面喷丸处理;零件冷挤压加工;在构件表面引入残余压应力,都是提高疲劳寿命的常用方法。材料强度越高,循环应力水平越低,寿命越长,延寿效果越好。在有应力梯度或缺口应力集中处采用喷丸,效果更好。 表面渗氮或渗碳处理,可以提高表面材料的强度并在材料表面引入压缩残余应力,这两种作用对于提高材料疲劳性能都是有利的。试验表明,渗氮或渗碳处理可使钢材疲劳极限提高一倍。对于缺口试件,效果更好。 5) 环境和温度的影响 材料的S-N 曲线一般是在室温、空气环境下得到的。在诸如海水、酸碱溶液等腐蚀介质环
18、境下的疲劳称为腐蚀疲劳。腐蚀介质的作用对疲劳是不利的。腐蚀疲劳过程是力学作用与化学作用的综合过程,其破坏机理十分复杂。影响腐蚀疲劳的因素很多,一般有如下趋势: a)载荷循环频率的影响显著 无腐蚀环境作用时,在相当宽的频率范围内(如200Hz以内),频率对材料S-N曲线的影响不大。但在腐蚀环境中,随着频率的降低,同样循环次数经历的时间增长,腐蚀的不利作用有较充分的时间显示,使疲劳性能下降的影响明显。 b)在腐蚀介质(如海水)中,半浸入状态(或海水飞溅区)比完全浸入更不利。 c)耐腐蚀钢材,抗腐蚀疲劳的性能较好;许多普通碳钢的疲劳极限则下降较多,甚至因腐蚀环境而消失。 d)金属材料的疲劳极限一般是
19、随温度的降低而增加的。但随着温度的下降,材料的断裂韧性也下降,表现出低温脆性。一旦出现裂纹,则易于发生失稳断裂。高温将降低材料的强度,可能引起蠕变,对疲劳也是不利的。同时还应注意,为改善疲劳性能而引入的残余压应力,也会因温度升高而消失。 第三章第三章 线性疲劳累积损伤理论线性疲劳累积损伤理论 若构件在某恒幅交变应力范围S作用下,循环破坏的寿命为N,则可以定义其在经受n次循环时的损伤为Dn/N n=0则D0,n=N则D1时,破坏。 构件在应力范围Si作用下经受ni次循环的损伤为Di=ni/Ni.则在K个应力范围Si作用下,各经受ni次循环则可定义其总损伤为 kikiNiniDiD11 破坏准则为
20、D=ni/Ni=1 若设计寿命为时间Td,时间Td内的损伤为D,则疲劳寿命为TfTd/D例例1:构件的:构件的S-N曲线为曲线为S2N=2.5*1010,设计寿命期间内的载设计寿命期间内的载荷谱如表中前二栏所列。试求累积损伤。假定对于荷谱如表中前二栏所列。试求累积损伤。假定对于100%P时的应力为时的应力为150Mpa.设计载荷Pi循环ni(106)Si(Mpa)Ni(106)Di=ni/NiP0.051500.8p0.10.6p0.50.4p5 1.111 0.045 120 1.736 0.058 90 3.086 0.162 60 6.944 0.72答案: D=0.985例例2 2:构
21、件的:构件的S-NS-N曲线为曲线为S S2 2N=2.5N=2.5* *10101010, ,若其一年所承受若其一年所承受的典型应力谱如表中前二栏所列。试估计其寿命。的典型应力谱如表中前二栏所列。试估计其寿命。Si(Mpa)循环数ni(106)Ni(106)Di=ni/Ni1500.011200.05900.1600.35D= _ 1.111 0.009 1.736 0.029 3.086 0.033 6.944 0.05 0.121答案 : T=1/D=1/0.121=8.27年当疲劳载荷谱不是用若干级应力范围水平的组合表示,而是用相应于一定时间期间的连续概率密度函数表示时,疲劳累积损伤度
22、的计算可表示为, 其中,S表示应力范围, 是应力范围分布的概率密度函数; N是应力范围为S的单一循环载荷作用下达到破坏所需的循环次数; LSLSLdSNSfNNdSSfNNdnD00)()()(SfS 是所考虑的整个时间期间内应力范围的总循环次数; 是在落在区间S, S+dS内的应力范围循环次数; 表示是对所考虑的整个时间期间积分。 LNdSSfNdnSL)(L应力范围长期分布为应力范围长期分布为Weibull分布的疲劳累积损伤计分布的疲劳累积损伤计算算Weibull 分布分布 应力范围在结构整个寿命期间的分布称为应力范围的长期分布。但是,在进行疲劳评估时,其疲劳寿命事先并不知道,因此,通常将
23、应力范围在一个适当的确定时间长度内有代表性的分布看作是应力范围的长期分布。这一时间长度称为疲劳载荷谱的回复期。 在船舶与海洋工程结构疲劳分析中,经常用两参数的Weibull分布表示应力范围S的长期分布,其概率密度和分布函数分别为SSSfSexp)(1 S0SSFSexp1)(式中, 称为尺度参数,称为形状参数。 S0疲劳累积损伤计算疲劳累积损伤计算 以回复期作为考虑的时间期间,将应力范围长期分布的表达式代入,得到相应的疲劳损伤计算式式中, 为伽玛函数。01expdSSSSANDmL)1(mANmL) (关于关于Weibull分布的两个参数分布的两个参数 在S-N曲线的两个疲劳实验参数A、m给定
24、的情况下,结构在一定回复期内的疲劳损伤取决于应力范围Weibull分布的形状参数 和尺度参数 。 形状参数 一般是根据结构所处的海洋环境、结构类型及响应特性以及构件在整个结构中的位置等因素来确定。到目前的研究结果表明,形状参数的数值一般是在0.7到1.3之间。通常是用某一海况资料对一批船舶进行疲劳载荷的长期分析,然后用Weibull分布对结果进行拟合,从而得到的 值。 在 值已知的情况下,尺度参数 可用回复期内疲劳载荷长期分析得到的对应某一超越概率的应力范围表示。通常是用超越概率为 的应力范围 来表示。其含义是,在该回复期内的全部 次应力范围循环中,大于 的应力范围仅可能出现一次。可得LLNS
25、SP/1)(LSLNLSLLLLSSFSSPSSPexp)(1)(1)( 再根据超越概率的定义,可得到尺度参数和形状参数的关系如下:/1lnLLNS应力范围为分段连续型分布的疲劳累积损伤应力范围为分段连续型分布的疲劳累积损伤计算计算分段连续分布模型分段连续分布模型 在船舶与海洋工程中,海洋波浪的长期状态通常看成是由许多短期海况的序列所组成。每一海况由表征波浪特性的参数以及该海况出现的频率来描述。对每一短期海况,通常是把波浪作为一个平稳正态随机过程来研究。相应地,船舶结构因波浪引起的交变应力过程也可以看成是由许多短期海况的序列所组成。对于航行在海洋中的船 舶而言,还应按航向进一步划分航行工况。对
26、每一海况和给定航向和航速,交变应力过程是一个均值为零的平稳正态过程,其相应的应力范围分布称为短期分布。根据平稳正态交变应力过程的统计特征,应力范围的短期分布可用连续的理论概率密度函数来描述。实际分析时,航速一般取为一个定值。综合所有海况和航向的应力范围短期分布以及各海况和航向出现的频率,就得到了应力范围的长期分布,其形式是分段连续的。 实际应用中,一般是将某一海况中在给定航向下的交变应力过程作为均值为零的窄带平稳随机过程,则根据随机过程理论可知,其应力峰值服从Rayleigh分布,概率密度函数为 式中,y表示应力峰值; 为交变应力过程的标准差。2222exp)(XXYyyyf y0X 设上述应
27、力交变过程的功率谱密度为 , 这一功率谱密度通常是用谱分析方法得到。记 , 分别为功率谱密度的0次矩和2次矩,则有 (n= 0, 2) 交变应力过程的标准差可由功率谱密度得到,根据随机过程理论表示为 )(XXG0m2m0)(dGmXXnn00)(mdGXXX 由随机过程理论,前述两个平稳随机过程的功率谱密度之间有下列关系:式中, 称为线性动力系统的传递函数或频率响应函数, 则称为响应幅值算子(RAO)。 GHGXX2)(H 2H 传递函数的物理意义是,它是在线性动力系统做圆频率为 的简单谐振时,响应过程的振幅与输入过程的振幅之比。当输入过程为波浪,响应过程为交变应力时,传递函数就是结构在圆频率
28、为 的规则余弦波作用下,应力幅值与波幅之比。 对于船舶结构疲劳评估的谱分析方法,波浪的功率谱密度可采用两参数的Pierson-Moskowitz谱(简称P-M谱),也就是国际船舶结构会议(ISSC)推荐使用的波浪谱。表达式为: 4454221exp24zzsTTHG 另外,为得到给定时间内的应力循环次数,要用到交变应力过程的跨零率 ,即单位时间内以正斜率跨越零均值的平均次数,其表达式为 当交变应力过程为窄带时,应力每跨越零均值一次就出现一个峰值,可以假设应力范围S和应力峰值y之间有以下关系: 或 0f02021mmfyS22/Sy 进而,利用概率论中随机变量函数的概率密度的计算方法,可得应力范
29、围的概率密度函数为2228exp4)(XXSSSSf S0疲劳累积损伤计算疲劳累积损伤计算 设所考虑的船舶在第i海况和第j航向中航行时间为 ,并用 表示在 期间的累积损伤度,则有 式中, 为该应力交变过程的跨零率; 为该航行状态期间内应力范围循环次数; 为该期间的短期应力范围分布。ijTijDijT00)(dSSfSAfTDSijmijijijijf0ijijfT0)(SfSij 将相应的短期应力范围分布 的表达式带入上式,可得 = 式中, 和 分别为该应力交变过程的标准差和功率谱密度的零次矩; 为伽玛函数。)(SfSij)21 (220mAfTDmXijijijij)21 (2200mmAf
30、TmijijijXijijm0) ( 设所考虑的长期时间期间为 ,相应的应力范围长期分布由 个海况组成,各海况出现的概率为 ,划分的航向数为 个,各航向出现的概率为 ,则 。该期间内总的疲劳累积损伤度 应为 = = TSnipHnjpjiijppTTTD SHninjijTDD11mijninjijjimfppmATSH011022)21 (随机谱与循环计数法 恒幅载荷作用下的疲劳寿命估算,可直接利用S-N曲线。变幅载荷谱下的寿命预测可用Miner理论加以解决。 现在进一步研究随机载荷的处理。 将不规则的、随机的载荷时间历程转化为一系列循环的方法,称为“循环计数法”。计数法有很多种,本课只讨论
31、简单、实用且与变幅循环载荷下的应力应变响应一致的简化雨流计数法。 简化雨流计数法,适用于以典型载荷谱段为基础的重复历程。既然载荷是某典型段的重复,则取最大峰或谷处的起至段作为典型段。 简化雨流计数方法如下: 1)有随机载荷谱中选取适合雨流计数的、最大峰或谷起止的典型段,作为计数典型段。如图中1-1或2-2段。 2)将谱历程曲线旋转90度放置。将载荷历程看做多层屋顶,假想有雨滴沿最大峰或谷处开始往下流。若无屋顶阻挡,则雨滴反向,继续流至端点。 3)记下雨滴经过的最大峰谷值,作为一个循环。 4)从载荷历程中删 除雨滴流过的部分, 对各剩余历程段,重 复上述雨流计数。直 至再无剩余历程为止。 第四章
32、 船舶结构疲劳载荷 ?你能想到哪些船体所受的载荷船体梁载荷总体载荷局部压力局部载荷垂 水 垂 水 扭向 平 向 平剪 剪 弯 弯力 力 矩 矩 矩海水压力内部货物压力问题1:如何来计算上面的载荷?问题2:上面的载荷在疲劳评估中都需要考虑吗? 载荷 第一部分 DNV规范的载荷一、船体梁载荷1.1船体梁垂向波浪矩式中 中拱波浪弯矩幅值; 中垂波浪弯矩幅值;BwwmrhwoBCLCKfM2,19.0mKN )7.0(11.02,bwwmrswoCBLCKfMmKN hwoM,swoM,wmKL m mBBCrf015 .0h8104100h Lhlog54.021.20 弯矩分布因子,从距尾垂线0.
33、40L至0.65L之间等于1,首尾垂线处等于0,其余部分进行线性插值; 理论船长 ; 水线处船的最大型宽 ; 方型系数; ,由超越概率 转化为 的因子, Weibull分布的基本形状参数, ; 波浪系数 wC1.2船体梁水平波浪弯矩LxCBTfMBactH2cos130. 022. 049mKN actT实际吃水 mx从尾垂线到所考虑横剖面的距离 m1.3大开口船扭矩 对于集装箱船这类具有甲板大开口的船舶来说,由扭转引起的翘曲应力在结构疲劳评估中必须加以考虑。二、局部载荷2.1 外部海水动压力载荷 船体外部水动压力按下式计算:dpeprpactwfBdrWactldpdTzkyCypzTByp
34、pp27 .016210)(2 .175135maxfwslkCkpLVkCkfws15. 08 . 0 尾垂线及其之后的部分距尾垂线0.2L至0.7L之间的部分首垂线及其之前的部分注:其它部分线性插值BBCC5 . 23BC3BBCC0 .43BC3skwz mTacty myy mB 4/ fwC8 . 0 mfkactTf radV kn 从基线到装载点的垂向距离,最大取值为 从中心线到装载点水平距离 ,最小为 从水线到舷侧顶部的距离,最大值取为 和 的最小值 横摇角 设计航速pr1 . 0wlwlzzzTact20 . 0wlzTactzwlwlzTactzzTactwlzTactzg
35、pzdTwl43 mddTppTactzw025. 13mt2.2由船舶运动引起的内部装载物的动压力由船舶运动引起的内部装载物的动压力船舶内部压载水或内部液体货物动压力的计算: 式中 垂向加速度引起的压力 横向加速度引起的压力 纵向加速度引起的压力slstsvaixapyaphapfp321max2/mKN1p2p3psxsy m m 舱内液面到考虑点的垂直距离shhaf15.0 Lhhlog54. 026. 205. 00025. 13mt202sinryytagaa2/sm03 . 0gaay2/sm舱内液体自由面中心到考虑点的纵向距离 舱内液体自由面中心到考虑点的横向距离 LVCLCaV
36、w/302.050/LCV 设计航速 , 同前VLwC75/50BcradKTcR25.025.1K见Notes30.7第28页说明 sGMkTrR30/2rk , 见Notes30.7第28页说明GMRZRryRTa2/2RZR横摇轴到舱中心的垂直距离 m202sinpxxlagaa2/smBxCgaa02 .0 方型系数 同前BC0aBCa /25. 00Notes30.7第29页PZPpxRTa2/22/smgLTP/80. 1 sNotes30.7第29页PZR纵摇轴到舱中心的垂直距离m2222maxzpzzrzvaaaaaRYRrzRTa2/2RYR横摇轴到舱中心的横向距离 mPXP
37、pzRTa2/2PXR纵摇轴到舱中心的纵向距离 mBzCgaa/7.00第二部分 JTP规范的载荷符号 :普通加速度参数 :由于纵荡产生的纵向加速度 :由于纵摇产生的纵向加速度 :由于纵摇产生的垂直加速度 :由于横荡和首摇产生的横向加速度 :由于横摇产生的横向加速度 :由于横摇产生的垂直加速度 :由于垂荡产生的垂直加速度0asurgeaxpitchazpitchaswayayrollazrollaheavea :水线处船的最大型宽 :特定装载工况吃水下水线位置的船宽,m :方形系数 :特定装载工况的 :型深 :1.2为不具备舭龙骨的船舶 1.0为具备舭龙骨的船 舶 :重力加速度,9.81 m/
38、s2 :初稳性高 :静水水线的波浪动压头 :规范船长,mBlocalBbCLCbCDbkfgGMWLhL :回转半径 :特定装载工况下的吃水,m :满足关于船舶尺度要求下的最大设计吃水 :纵摇固有周期 :横摇固有周期 :最大营运前进航速 :船舶速度(节),应取:0,为尺度要求和 强度评估 0.75,为疲劳强度 :纵向坐标,m :横向坐标,m :垂直坐标,mgyrrollrLCTscTpitchUrollUV0Vxyz :参考点的纵坐标,取液舱顶部的长度方向的中点,m :参考点的横坐标,取液舱顶部的宽度方向的中点,m :参考点的垂直坐标,取液舱的最高点,m :横摇角 :纵摇角0 x0y0z2.1
39、 运动运动 2.1.1 概述概述 船舶运动的特性包络线值为10-8可能性水平。2.1.2 横摇运动横摇运动2.1.2.1 自然横摇周期 由下式求得: 2.1.2.2 横摇角 由下式求得: rollUGMrUgyrrollroll30.2sbkrollfUB)025. 025. 1 (7550弧度2.1.3 纵摇运动纵摇运动2.1.3.1 特定纵摇周期 由下式求得: pitchULfgfUTVpitch)1 (26 . 0s式中: )67. 0525(0 . 10LVVfVscLCTTTf2.1.3.2 纵摇角 由下式求得: :取 值,但不应小于10节1801)(96025. 01LCVb弧度式
40、中:1VV2.2 船舶加速度船舶加速度2.2.1 概述概述 在六个自由度方向上的运动产生的平移加速度特性包络线值可计算得到。横向和纵向加速度包括由于横摇和纵摇引起的重力加速度。2.2.2 普通加速度参数普通加速度参数2.2.2.1 普通加速度参数由下式求得: )600344 . 2)(47. 058. 1 (20LLLCab2.2.3 垂向加速度垂向加速度2.2.3.1 任何位置的垂向加速度 由下式求得: m/s2 m/s2 :参见2.2.3.2 :参见2.3.3.22.2.3.2 对于疲劳强度: 应取0.45 va222zrollzpitchheaveprobvaaafaheaveagafV
41、0zpitchaLxULpitch45. 0)2()3253 . 0(2m/s2zrollayUroll2)2(2 . 1m/s2式中:probfVfprobf)10002 . 1 ()(2LCCfbLCbV2.2.4 横向加速度横向加速度所有位置的横向加速度包络线 由下式求得: 式中: m/s2 或 ,取其大者 m 应取0.5ta22)sin(yrollswayprobtagafam/s2swaya03 . 0 gayrollarollrollRU2)2()24(LCrollTDzR)2(Dz probfm/s2 对于疲劳强度:2.2.5 纵向加速度纵向加速度所有位置的纵向加速度包络线 由下
42、式求得: : 或 ,取其大者 m 对于疲劳强度: 应取0.5 应取1.7galn22ln)sin(3257 . 0 xpitchsurgeprobgagLafam/s2式中:surgea02 . 0gaxpitchapitchpitchVRUf2)/2(pitchR)24(LCTDz)2(Dz probfVf2.3 船体梁载荷2.3.1 垂直波浪弯矩中拱和中垂的垂直波浪弯矩包络线 和 垂直波浪弯矩沿船舶长度的分布因数 对于疲劳强度 应取: hogwvMsagwvM由下式求得: bwvvwvprobhogwvBCLCffM219. 0)7 . 0(11. 02bwvvwvprobsagwvCBL
43、CffMkNm kNm 式中: vwvfvwvf 0.0 在尾尖 0.1 在尾尖开始0.1L处 1.0 在尾尖开始0.4L到0.65L之间 0.1 在尾尖开始0.9L处 0.0 在首尖 通过线性内插法获得中间值,参见下图分布分布 尾尖 首尖vwvf 波浪系数 wCprobf:应取0.5 2.3.2 水平波浪弯矩水平波浪弯矩 沿船舶长度的水平波浪弯矩的分布因数,取值同2.3.1中 取值 应取0.5vwvf水平波浪弯矩包络线 hwvM由下式求得: bLCwvhwvprobhwvCTLCfLfM2)20003 . 0( kNm式中: hwvfwvC参见2.3.1probf2.3.3 波浪动压力波浪动
44、压力波浪动压力 应取 和 的大者。 和 由下式求得: kN/m2 dynexP1P2P1P2P2111111)75(4135)(2 . 1)75(41352fBBfzTBBPffPlocalLClocalPnprob21212)27 . 0(1425. 08)27 . 0(148 . 025. 0826fTzBCfBfTzCBCfBffPLClocalbTlocalLCwvlocalbTlocalPnprobkN/m2式中: localB:特定装载工况吃水下水线位置的船宽,m 当 当 , 在尾尖后 , 在尾尖开始0.2L到0.7L , 在首尖前 的中间值通过线性内插法获得wvsCfP)8 .
45、03(1122lnln1ffffffVgg) 14(25. 0localVByf) 14(localVByf2flocalBy25. 0localBy25. 0scLCTTTfbbCC33. 1bCbbCC33. 1sfsf 1.0,在尾尖后 = 0.7,在尾尖开始0.2L到0.7L 1.0,在首尖前 的中间值通过线性内插法获得 对于疲劳强度,波浪动压力的振幅 (如后图)由下式求得: kN/m2 当 或 之间小者 kN/m2 水线处 kN/m2 当 或0之间大者 当位于水线和 之间时,由线性内插法获得。式中:gflngflnampexP0WLP5 . 0dynexPWLLChTzDWLLChT
46、zWLLChTzampexPWLh10/WLP 取静水线处的 ,kN/m2 取 和 的大者,kN/m2 :参见2.3.3 ,kN/m2 , 在0.7向后 , 在首尖前 的中间值可以由线性插值法求得 :参见2.3.3,kN/m2 WLPmaxexPmaxexP1P2P1P5 . 0probf0 . 111Pnf0 . 15 . 1VfVf2P5 . 0probf0 . 121Pnf对于疲劳强度的波浪动压力振幅的对于疲劳强度的波浪动压力振幅的 横向分布横向分布ampexP2.4 货物产生的压力对于疲劳强度,在一个邻近舱室为空舱的舱壁上,液舱内部动压力振幅 可由下式求得: kN/m2 式中: kN/
47、m2 kN/m2 kN/m2 ampinPgingtintvinvampinPfPfPfPlnlnzzaPvvin0yyaPttin0 xxafPggullgin0lnlnln :参考点的纵坐标,取液舱顶部的长度方向的中点,m :参考点的横坐标,取液舱顶部的宽度方向的中点,m :参考点的垂直坐标,取液舱的最高点,m 的取值如下表0 x0y0zvftfgfln第五章 船舶结构疲劳应力的计算 思考 船舶结构应力的分类?船体梁应力(1) 板架应力(2) 骨材弯曲应力(3) 板弯曲应力(4) 应力分量定义应力分量定义 问题:每一种应力成分的计算方法?一、船体梁弯曲应力船体梁弯曲应力 垂向弯矩和水平弯矩
48、将产生下面的应力范围垂向弯矩和水平弯矩将产生下面的应力范围Ns ,woh ,woaxialvInzMMK0310CHaxialhgIyMK3102你 能 想 到的 理 论 基础?对船体梁弯矩作用下的应力进行合成,hgvvhhgvg222式中: hgv,分别为垂向整体应力范围和水平整体应力范围); 2mm/N( s ,woh ,woM,M中拱弯矩,中垂弯矩(kNm); NI4mm关于水平中和轴的剖面惯性矩(序计算得到;),利用程CI关于垂向中和轴的剖面惯性矩( 4mm),利用程序计算得到;y 计算点到垂向中和轴的距离( m); 两方向弯矩间的平均相关系数,取 10.vh; 提问:物理意义?二、板
49、架弯曲应力的计算板架弯曲应力的计算a.计算双壳体板格与横舱壁相交处的纵向第2弯曲应力方程为 当板架的长大于宽时,沿纵向的应力为baaebiirbpKK224abiiba当板架的宽大于长时,沿纵向的应力为aaebirapKK224baiiabb.计算双壳体板格与纵舱壁相交处的横向第2弯曲应力方程为 当板架的长大于宽时,沿横向的应力为bbebirbpKK224abiibac当板架的宽大于长时,沿横向的应力为 依赖于外形比例的系数,=1.0,初始情况对照规范查表11.1babebiirapKK224baiiab式中: K轴向载荷引起的应力集中系数,即 axialK; bK 板架宽度( ); 板架长度
50、( ); 计算点到板架中和轴的距离( );骨材弯曲引起的应力baarmmmmmmrZlEImKZsMKASaxiallateralA222SlateralAZMK2prpslM12201662.)lx()lx(rplx 0式中: 横向载荷作用下引起弯曲时用的应力集中系数; 计算点处的弯矩( ) 横向动压力,分 和 。 为海水动压力值, 为内部动压力值。这里计算时都取为1 ; 骨材间距( ); 加强筋的有效跨距,具体计算参见规范第23页图示3.4( ); 纵骨连带板/加强筋连带板的在翼板顶部的剖面模数; 力矩插值系数; 到热点的距离();lateralKMmmN 310pepipepip2m/K