0 引 言

大型集装箱船具有主尺度大和甲板开口大的特点，相较于其他船型需承受更大的总纵弯曲和扭转变形，而高强度钢的使用会导致板厚较小，从而使得船体结构的整体刚度和局部刚度偏“软”，进而易引起严重的结构疲劳破坏，因此对船体结构疲劳强度进行准确预报和评估研究非常重要。对于船舶结构疲劳问题，目前常用的解决方法有规范简化算法和直接计算法2种，其中：规范简化算法是在大量假定和经验公式的基础上得到的，在前期设计时较为简便，但无法反映真实海况和船体结构特点的影响；直接计算法是在疲劳载荷概率模型的基础上得到的，有谱分析法和设计波法2种。谱分析法分别对波浪频率范围内的多个规则波和浪向角进行计算，计算精度较高，计算量较大；设计波法则是通过对船体波浪载荷进行长期分析，寻找一种典型的规则波，使之产生与长期预报值相等的波浪载荷，从而计算在该设计波载荷循环作用下的疲劳强度，其优点是计算量小且能反映结构细节，缺点是仅仅以一种主控载荷来表达多种载荷组合的作用，在精度上稍差。

英国劳氏船级社（Lloyd′s Register，LR）船体疲劳计算评估（FDA）分为 3个等级水平[1]。水平 1[2]是基于结构构件和综合已有的航行经验得出的形式的比较。水平2 [3]和水平3[4]是简化及完全谱分析的计算程序，其中：水平2适用于纵骨节点疲劳分析；水平3采用波浪载荷直接计算和基于精细有限元模型的应力计算，可针对特定结构部位给出较为准确的疲劳评估。运用谱分析方法、Palmgren-Miner 线性累计损伤原理和S-N（应力-循环次数）曲线，在计算各短期海况的波浪诱导应力和疲劳损伤响应的基础上，评估船体结构在设计寿命期内的长期累计损伤和疲劳寿命。

1 疲劳谱分析理论基础

在某个给定的浪向、波频、航速及装载工况下，根据线性理论，在单位波幅规则正弦波下应力时域响应可表示为

式(1)中：σ(t)为应力响应；Ci为 ith加载过程结构影响系数； hi( t)为ith加载过程应力响应；n为加载总数。

在某个给定的工况、航速、浪向下，海况以有义波高Hs和过零平均周期Tz表示，短期频域应力统计可表示为

式(2)中：h i (ω)为应力响应； S∈(ω)为波浪谱（ISSC波浪谱）；Ci Cj为 ith 和jth加载过程结构影响系数；h i(ω)为ith加载过程载荷传递函数；hj(ω)为 jth加载过程变化复杂的载荷传递函数。

成联方：我不喜欢用“丑书”一词来形容这类风格的书家，我喜欢用“乡土主义”来形容他们。“乡土主义”虽然有抄袭文学、美术之嫌，但是，从这些书家成熟的时间节点以及取法对象等等方面来看，我觉得“乡土主义”更准确一些，也更好理解。我认为“美”与“丑”是相对的，这一点中西艺术史中都可以找到大量的例证。晋代郭璞在《尔雅》中就说“美恶不嫌同名”。所以，在书法上存在这种现象，你认为“丑”的，别人却认为是“美”的，你认为是“雅”的，别人却认为是“俗”的。至于对当下这些所谓的“丑书”进行评论，最好不用“丑”字，要去寻找一些新的、真正贴切的评价方式也许更好。

用于计算响应谱标准差、带宽及跨零频率的谱距计算式为

根据表7显示：标准舞练习组女生在上肢快速反应测试中使用时间平均减少了7.35s，拉丁舞练习组女生平均减少了8.17s，两组之间没有显著性差异。触杆反应测试和十字绕杆跑测试中，标准舞练习组女生使用时间分别平均减少了3.06s和3.09s，拉丁舞练习组女生分别平均减少了2.88s和2.67s，两组之间差异不显著。全身协调测试中两组之间差异显著，因标准舞练习组女生比拉丁舞练习组女生前测成绩相对较好，且无训练基础但协调能力较好的人群在经历训练后协调能力能够得到提升，但相对协调能力较差的人群，其提高程度一般较低[19]。

目前，百安居通过第三方物流对门店的配送形式是，由三大区域物流中心配送至全国的38家门店，配送线路总计达到120条，非常烦琐。简单来说，也就是现在北京的门店、深圳等门店的配送有时也会需要上海物流中心进行配送。

式(4)中：p( S)为应力范围概率分布函数；0σ为应力范围标准差；S为应力范围。

在某个给定的工况、航速、浪向及海况下，根据Palmgren-Miner理论，热点的短期疲劳损伤可表示为

1.语言激趣：“兴趣是最好的老师”，教师入情入境的语言往往能沟通学生的心灵，调动学生的积极性，活跃作文形式，达到较好的作文教学效果。如教完《伙伴》一课时，我想让学生续写一篇《我们是伙伴》。我这样用语言导入：自从上回乌鸦把猎狗救了之后，马哈多和猎狗还有乌鸦一直幸福的生活着，可是有一天马哈多去打猎的时候遇到了危险，猎狗和乌鸦是如何帮助马哈多脱险的，接下来会发生什么故事呢？在这样的语言设计中，激发了学生创作的欲望，能使学生产生强烈的表达思想的欲望。这样学生写出来的东西也一定是色彩纷呈的。

式(5)中：D为短期疲劳损伤；i为Si循环应力下达到破坏所需的循环次数； n( Si)为Si循环应力实际循环次数；nt为给定海况下总应力循环次数。

2 波浪频域计算和疲劳累计损伤组合计算的方法

LR的疲劳分析采用的是全概率谱分析法，由航程预报、水动力计算、有限元分析及累计疲劳损伤计算等4部分组成。FDA水平3计算流程见图1[5]。

基因编辑技术对农作物基因改良的速度在加快，一些消费者和环保组织开始担心该领域的发展已经大大超过监管机构的监管范畴。该项技术的倡导者和反对者都认为，具有30年历史的转基因作物审查法律的框架未能跟上CRISPR和TALENs等基因编辑技术的创新步伐。

2.1 航程预报

热点应力S-N曲线参数见表3。

Bone resorption was measured as the percentage of trabecular resorption surface(TRS%),which represents the percentage of the irregular,uneven surface of the trabecular bone relative to the total trabecular bone surface.TRS%can also be used to determine osteoclast activity.

图1 FDA水平3计算流程

图2 北大西洋有义波高频率概率分布图

表1 LR箱船疲劳分析工况（＞5000TEU集装箱船）

FDA 工况号 重箱重量/(t/箱) 空箱重量/(t/箱) 重箱装载率/% 空箱装载率/%FDA 装载工况FDA LC5.1 12.5 1.8 70 23 FDA LC5.2 10.0 1.8 78 8

表2 装载工况占有率（＞5000TEU集装箱船）

装载工况 整个服务寿命所占百分比/%FDA LC5.1 FDA LC5.2休息28 57 15

2.2 水动力频域计算

船舶运动及舷外水压力采用基于三维势流理论边界元法的线性刚体频域计算软件Wave-FD计算。计算的波浪频率为圆频率在0.2～1.2rad/s以每0.04rad/s为间隔，浪向角在0°～180°以每20°为间隔，0°为随浪，90°为横浪，180°为迎浪。根据2种疲劳计算工况，分别计算各工况、各浪向角及各波浪频率下的船舶六自由度运动，得到单位波幅下的幅值响应函数（RAO）。

2.3 有限元分析

采用单位力分析法对疲劳分析点进行精细网格有限元分析。即先将船体外壳按照一定的大小划分为片，在每块片上施加单位力，通过有限元计算得到每块片在施加力之后的应力响应结果，最终的响应以结合实际每块片上的水动力载荷及单位力响应结果得到。该方法可极大地减少计算量，有效地得到各海况、各工况下的应力响应谱。涉及液舱的计算可采用以下2种方法：

(1) 采用单位加速度分析法，在液舱重心位置施加单位加速度，根据实际船舶运动加速度叠加到最终应力响应谱中；

(2) 将液舱载荷嵌入到计算模型中，在计算单位力时将其包含进去，不单独考虑。

首先，每篇开场白第一段中的话题都跟本场辩论拟解决的问题相关。Yoni Cohen-Idov用五句话(五个话题和一个从属话题)阐述了本场辩论的总话题(禁止录像监控)拟解决的问题，即北京及其它城市都面临着扩大录像监控范围和规模的问题。

经过实际的计算对比发现：对于不在液舱范围内的疲劳点，采用2种方法所得计算结果误差很小；而对于存在于液舱边界中的疲劳点，采用第一种方法所得计算结果比采用第二种方法所得结果更保守，但二者相差≤15%。

2.4 疲劳损伤及年限计算

疲劳强度由S-N曲线表达式表示为

结合S-N曲线及Palmgren-Miner累计损伤原理，同时考虑疲劳点的厚度效应及焊接集中效应，计算疲劳点的累计损伤及疲劳年限。

假设响应过程是窄带分布，应力范围可用Rayleigh分布表示，即

本文所研究的船要求满足在北大西洋运营20a的疲劳寿命。根据北大西洋海域的航程预报，对该海域的有义波高进行统计分析，得到北大西洋有义波高概率分布图见图 2[6]。基于该船所在航道的信息为美国东海岸到挪威的航道信息，波浪环境为各浪向同概率，波浪以余弦平方的形式传播，航速为75%最大航速，计算工况见表1，装载占有率见表2。

表3 热点应力S-N曲线参数

参数 钢级l o g K hs m′N≤1 0 7 m′N＞1 0 7平均 设计A,A H 3 2, A H 3 6,A H 4 0, B 1 4.0 3 1 3.6 4 3.5 5.5自由边D, D H 3 2, D H 3 6,D H 4 0, E, E H 2 7 1 4.0 3 1 3.6 4 3.5 5.5 E H 3 2和F H 3 2 1 5.1 8 1 4.7 8 4 6 E H 3 6和F H 3 6 1 5.2 8 1 4.8 8 4 6 E H 4 0和F H 4 0 1 5.3 3 1 4.9 3 4 6焊接 所有钢级 1 2.6 4 1 2.1 9 3 5 1 0 7循环次数下应力范围Δ s 0/(N/m m)平均 设计2 l o g N标准差1 0 2.0 0 7 8.9 1 0.2 2 1 8 1 0 2.0 0 7 8.9 1 0.2 2 1 8 1 1 1.1 2 8 8.2 7 0.2 0 0 0 1 1 7.7 1 9 3.5 0 0.2 0 0 0 1 2 1.1 4 9 6.2 3 0.2 0 0 0 7 5.6 3 5 3.8 0 0.2 2 1 8

式(7)中：KSPR为弹振折减系数；Tf为不包含弹振影响的刚体疲劳年限（FDA 水平3结果）。

考虑到超大型集装箱船的航速较高，艏部外飘显著，且一阶固有频率比常规船型低，容易引起颤振和弹振。相对于颤振，弹振响应持续的时间更长，对船体结构的疲劳影响较大。该船要求满足LR FDA SPR符号，最终的疲劳年限计算结果考虑弹振的影响。

pH值参数统计与上述方法不同，应按下列步骤进行，先将土壤pH值换算成[H+]平均浓度进行统计计算，然后再换算成pH值，其公式为：[H+]=10-pH，[H+]平均浓度=∑10-pHn/n,pH=-lg[H+]平均浓度。

3 实船疲劳强度评估

全船采用粗网格建模，疲劳区域采用板厚乘以板厚的精细网格建模。全船有限元模型为自由体的静力计算，避免六自由度约束不当造成应力发生变化，边界条件采用惯性释放技术，充分模拟计算船舶实际情况。

舱口角隅可根据其与箱角的宽度选用圆形、椭圆形和负角隅钥匙孔等3种形式，三者应力集中系数为圆形大于椭圆形、椭圆形大于负角隅钥匙孔。但是，负角隅形式施工比较麻烦，且需在主甲板开口处加上挡雨板，防止雨水进入货舱，因此在应力集中系数要求不高的位置多采用椭圆形及圆形角隅。对于双岛型船，机舱和燃油舱上甲板处多采用负角隅设计，因为处于扭转边界的应力水平较高。纵向舱口围板及主甲板位置距离中和轴最远，纵向应力最大，本文选取机舱前端壁、燃油舱后端壁、船中水密舱壁、船中支撑舱壁、纵向舱口围及主甲板角隅为评估对象。

3.1 水动力响应

船舶最大航速为23.6kn，选取17.7kn（75%最大航速）为计算航速，分别计算出2种工况下各浪向对应的频域幅值响应函数RAO及对应的舷外水压力。由于计算结果繁多，这里仅列举LC52工况的横摇RAO（见图3）及浪向角60°和波频0.36rad/s情况下的舷外水压力（见图4）。实际计算的横摇阻尼系数结合LR推荐的范围及试算结果选取。由图3和图4可知：当浪向角为60°时，单位幅值下的横摇最大幅度为2.1°，认为计算结果可靠，可用于下一步的结构匹配。

图3 LC52工况的横摇 RAO

图4 航速17.7kn，浪向角60°，波频0.36rad/s情况下的舷外水压力

3.2 疲劳分析

采用LR的计算软件ShipRight FDA3计算模块，划分外板PANEL，加载单位舷外水压力及单位燃油舱纵摇加速度和横摇加速度。通过有限元计算，得到每块PANEL单独受力情况下船体的响应谱。根据有限元外板PANEL尺度，对水动力结果进行相应的网格划分，对边界位置进行插值，得到对应的每块PANEL的计算结果。液舱加速度根据水动力计算所得的RAO的两阶导数得到，根据转动加速度与平动加速度的关系，划分到对应的有限元网格上。同时，结合航程预报的概率系数、线性放大单位力及单位加速度下的应力RAO，根据S-N曲线及Palmgren-Miner累计损伤原理，得到最终的刚体疲劳总损伤值，其倒数为刚体疲劳年限。本文分析的4个热点的疲劳强度评估结果见表4。

表4 4个热点的疲劳强度评估结果

疲劳强度计算点 位置描述 FDA3疲劳寿命/a SRP 疲劳寿命/a 修改 修改后疲劳寿命/a纵向舱口围 18 0.67 12 增加中间支撑 21 Spot1 机舱前端壁FR79处（见图5） 横舱壁（20AH R900+R900）上甲板 41 0.66 27 — —4.8 — 4.8 40EH40 R650+R600面内及边缘打磨21 FR221纵向舱口围(R800+R450) 52 0.80 41 — —Spot2水密舱壁FR221及FR223处（见图6和图7）FR221主甲板（R1100+R560） 28 0.83 23 — —FR223纵向舱口围(R800+R450) 29 0.80 23 — —FR223主甲板（R1100+R560） 22 0.83 18 R800+R300 32 FR239纵向舱口围(R800+R450) 37 0.80 30 — —Spot3支撑舱壁FR239及FR241处（见图8和图9）FR239主甲板（R1100+R560） 23 0.83 19 R800+R300 38 FR241纵向舱口围(R800+R450) 31 0.80 25 — —FR241主甲板(R1100+R560) 42 0.83 35 — —Spot4纵向舱口围 47 0.49 23 — —燃油舱后端壁FR329处上甲板（见图10） 横舱壁(R900+R800) 4.2 — 4.2 R650+R600面内及边缘打磨 20(R4200+R970) 12 0.51 6 R4400+R1000细节见图11 22

图5 Spot1热点细化图

图6 Spot2热点细化图（FR221）

图7 Spot2热点细化图(FR223)

图8 Spot3热点细化图（FR239）

图9 Spot3热点细化图（FR241）

图10 Spot4热点细化图

由计算结果可知，大部分结果不满足疲劳寿命要求（20a），这是由于该船与以往设计的超大型集装箱船（20000TEU）相比，双岛间增加了一个货舱，机舱与燃油舱之间的距离增加，扭转应力显著增大，燃油舱后端壁尤为明显。具体的修改措施及修改之后的疲劳年限如表 4所示，负角隅形式开孔，在横舱壁与主甲板连接处应力集中明显，更改为采用增加间距引入大圆弧设计，减少应力集中；中舱段水密舱壁前档及支撑舱壁后档主甲板处修改形状，将椭圆形的尺寸改为（R800+300），纵向过渡内收满足20a疲劳年限。燃油舱壁处作为扭转边界应力集中尤为严重，燃油舱处主甲板具体的修改细节见图11，通过增加纵向舱口围处的中间横向支撑来减小变形，并在建造过程中采用表面打磨，控制建造监控，最终满足疲劳年限要求。

图11 燃油舱后端壁上甲板修改后细节图

4 结 语

通常在一个波长范围内需用 5个网格来精确表达船体运动和舷外水压力，无限水深波，波频越大波长越短。由北大西洋有义波高概率分布图可知，概率分布集中在中段海况，考虑到水动力计算耗时较长，可根据需要适当放大水动力网格，精确计算出概率分布占比较大部分的应力谱即可。

收集到2016～2017年本院门急诊药房成功干预的不合理处方共512张。依据《处方管理办法》、《临床用药须知》、《中国国家处方集》、药品说明书、疾病诊疗指南及专家共识等，对处方规范性、用药适宜性等方面进行分析和总结。运用Excel 14.1.0软件进行分类统计和数据处理。

由计算结果可知，弹振折减系数在扭转边界位置达到50%，对于超大型集装箱船而言，在进行疲劳强度校核时弹振现象需合理考虑。

在扭转边界、机舱前后舱壁、燃油舱前后舱壁及前端艏货舱舱壁处，建议采用减小应力集中效果最好的负角隅结构形式，横舱壁开孔与主甲板连接处采用大间距圆弧过渡。双岛之间的纵向舱口围及主甲板角隅建议采用椭圆形式，在保证角隅宽度的情况下纵向圆弧尽量内收，达到减小应力集中的效果。

对于超大型集装箱船的疲劳谱分析，考虑到不同装载工况和不同海况对疲劳损伤的累计作用，计算结果全面可靠，可为同类集装箱船的研发和设计提供参考。

【 参 考 文 献 】

[1] LR.Rules and regulations for the classification of ships[S].2017.

[2] LR.FDA Level 1 procedure structural detail design guide[S].2009.

[3] LR.FDA Level 2 procedure[S].2002.

[4] LR.FDA Level 3 procedure[S].2009.

[5] LR.Ship right analysis the FEM model user guide[S].2015.

[6] IACS.Rec 34 standard wave data[S].2001.