导缆孔的主要作用是限制缆绳导出位置,其广泛应用于船舶、海洋工程的系泊．船舶配备的拖索力要远大于其系泊索作用力[1],因此应急拖带中的导缆孔一般直接承受巨大的缆绳负荷.在导缆孔的使用时,发现导缆孔发生塑形变形甚至断裂,导致船舶无法正常工作甚至造成严重的安全事故．因此,针对导缆孔进行应力、应变、疲劳分析具有十分重要的意义．目前,在导缆孔新产品开发过程中,文献[2]通过经典理论力学计算,给出一种巴拿马导缆孔的设计计算方法,但此方法适用范围窄、分析精度低、获取信息量小、计算过于保守导致结构笨重；文献[3]主要针对巴拿马导缆孔数量和位置要求进行分析；文献[4]针对带缆桩、导缆孔等主要系泊拖带属具标准进行梳理；文献[5]分析了导缆孔位置对系统稳定性的影响；文献[6]研究巴拿马导缆孔选型问题，但没有提供优化结构方案和设计思想．目前较多研究人员和厂家通常采用经典的经验公式对导缆孔进行设计,只能采用较高的安全系数以保障其安全性,使得导缆孔尺寸和重量越来越大,理论的承载能力和实际的承载能力存在较大差异,但对于结构复杂、新材料或新工艺的导缆孔,此方法则可能存在较大误差．

文中首先分析导缆孔受力过程,运用Abaqus对其进行静力学分析；然后以HyperWorks 为平台建立导缆孔有限元模型,对导缆孔进行拓扑优化设计．最后分别以Abaqus和Femfat验证力学性能和疲劳强度,综合验证优化后结构的正确性．

1 有限元建模

1.1 结构及工况分析

文中主要针对BC450 GB11586-89型号的巴拿马导缆孔进行分析和优化,其结构关于水平和竖直方向对称,为简化模型,在导缆孔四分之一内,取a,b,c 3个方向作为分析的3种工况,如图1．

式中，各符号意义同式(1)；下标i分别代表各产业或者行业；*代表各产业(或行业)生产函数里的产出弹性参数。通过将L作为分母转变为人均形式的生产函数，式(1)和式(2)可改写为：

  

图1 导缆孔受力Fig.1 Chock′s force

《巴拿马运河导缆孔》(GB11586-89)是目前国内最为广泛的导缆孔标准,但该标准没有给出相关导缆孔的安全载荷[7],根据BC450 GB11586-89巴拿马导缆孔一般使用情况,文中选取缆绳拉力为320 kN．

根据拓扑优化结果,建立新的网格模型,重新对导缆孔进行应力应变分析,应力应变分析结果分别如图7、8．

目前，从全国总体来看，“纲要”课的课堂教学质量和教学效果质量有明显的改善和进步，学生的获得感不断增强。但是各校的水平还是层次不齐的。一所高校的调查显示，大学生对学习“纲要”课缺乏主动性，64.7%的学生回答对此课的学习积极性不高或一般，46.9%的人只是把老师讲的重点内容划在书上；12.7%的人只是听一听而已。很多学生来上课，手里光拿着一本书，纸笔全不带。学过“纲要”课程后，还有近1/3的学生对中国为什么选择了马克思主义不太明白，[2]上述现象绝非孤立，而是有代表性的。

文中以通用CAE 软件HyperWorks 为平台,对导缆孔进行优化设计,其中Optistract 用来作为优化求解器对导缆孔进行拓扑优化设计．

1.2 导缆孔有限元模型

影响因子方法用来评估零件工艺上的初期破裂以预测零件的疲劳寿命．经过有限元计算得到材料在各节点的S-N曲线．局部零件的材料属性、普遍的应力,尤其是应力梯度影响着S-N曲线的参数变化．如图9零件局部S-N曲线,主要通过以下3个参数来描述:疲劳强度极限σaf,C;曲线斜率Kc;极循环次数Ncf,C．通过这3个参数的确定可以得出局部应力状态,进而完成疲劳强度的计算．3个参数受到很多因素的影响,主要包括:凹槽效应、平均应力、表面粗糙度、工艺参数、回火工艺参数、表面处理方法、几何尺寸、温度、统计学等．在对巴拿马导缆孔进行疲劳分析时,主要考虑应力影响、表面粗糙度影响、温度影响、统计学等影响,表面处理选择退火．

  

图2 导缆孔模型Fig.2 Model of chock

2 优化模型

根据应急拖带样本[8],其基本模式为拖索一端与拖轮或者拖带船舶连接,另一端通过位于被拖带船舶上的导缆孔系于带缆桩或者系泊绞车之上．缆绳作用于导缆孔时,与导缆孔孔体表面圆弧段产生接触,从而产生压力．缆绳直径取60 mm,宽度取缆绳直径的60%～75%,假定该压力沿缆绳方向在接触面上均匀分布,接触面长度依照缆绳角度确定，从而计算出接触面积取18 286 mm2左右,此时对应压强为35 MPa．

(2) 对导缆孔进行拓扑优化后,优化后导缆孔强度比原始结构受力性能方面,在a,b,c 3个方向分别提高46%、20%、13%,整体受力情况显著改善．

E=ηαE0,

(1)

式中: E 为当前材料弹性模量; η为材料变密度; E0为材料的真实弹性模量;α 为离散参数．

寻找输入：在需求分析和目标确定以后，寻找输入成为设计 《汽车英语》课程交际任务的出发点。需求分析的结果对于交际任务中输入的选择具有重要的作用，也就是说教师要尽量选择汽车行业内使用的真实文本作为输入，主要通过互联网和去进口品牌汽车的4S店中搜集。这种输入对学习者来说是最实用的，

优化时以材料的变密度为拓扑设计变量,这样结构拓扑优化问题就转换为材料的最优分布问题,优化目标为柔度最小,约束是导缆孔的体积分数,则导缆孔拓扑优化的数学模型为:

 

3 结果分析与讨论

3.1 导缆孔应力应变计算分析

导缆孔使用的材料是G230-450，为铸铁类,故采用第一强度理论，即最大拉应力理论进行校核．

图3、4分别给出了导缆孔应力云图和位移云图．

  

图3 导缆孔应力云图(单位：MPa)Fig.3 Stress profile of chock(unit：MPa)

由图3可知,应力最大值分别出现在加强筋处、孔体表面处、加强筋处,3种工况对应的最大应力值分别为217、196、244 MPa,均已接近或超出材料的屈服极限230 MPa,因此,导缆孔结构需进行优化设计．由图4可知,导缆孔最大位移处在受力方向的孔体表面,且最大位移为分别为0.047 5、0.382、0.115 mm．此3种位移较小,满足刚度要求．

  

图4 导缆孔位移云图(单位：mm)Fig.4 Displacement profile of chock(unit：mm)

3.2 导缆孔优化结果

图5(a)、(b)分别为导缆孔拓扑优化结果和重建的三维模型．其中,优化后导缆孔形状为内、外双层薄壳,同时,外层壳采用薄壁斜支撑方式,取消了传统的加强筋．

  

图5 拓扑优化结果Fig.5 Result of topology optimization

根据优化后的导缆孔模型,采用origin对内层壳截面曲线进行拟合．结果表明,内层壳截面曲线为正高斯曲线,在结构受到相对均匀的力时,正高斯曲线具有很好的稳定性[10]．图6给出了优化后内层壳截面曲线图,且内层壳截面曲线方程为：

 

(1)

  

图6 优化后内层壳截面曲线Fig.6 Optimizated inner shell section curve

“截至7月30日，今年我们的吞吐量已有72399标箱。”陈淑卫介绍说，为了更好地融入长江经济带，他们有的放矢，针对航线不足、江海联运等现状进行破题。

根据Abaqus静力学分析结果,使用疲劳分析软件Femfat对导缆孔进行疲劳分析．其中,材料强度极限为450 MPa,屈服极限为230 MPa,进而得出S-N曲线、σ-ε曲线和Haigh图．

  

图7 优化后导缆孔应力云图Fig.7 Stress profile of the optimized chock

  

图8 优化后导缆孔位移云图Fig.8 Stress distribution graph of optimized chock

3.3 导缆孔优化前后对比分析

3.3.1 静力学分析对比

通过图3和图7,图4和图8对比可以看出,优化后导缆孔的应力云图和位移云图的分布情况和优化前(图3、4)基本一致．应力最大值分别出现在3处受力方向的孔体表面处,且3种工况对应的最大应力值分别为116、157、212 MPa,较优化前分别降低了46%、20%、13%,3种工况中,最大应力值相对于优化前均明显减小,说明导缆孔强度提高．3种工况对应的最大位移值分别为0.084 7、0.142、0.227 mm，较优化前,a方向增加0.037 2 mm,b方向减少0.24 mm,c方向增加0.112 mm,位移尺寸远远小于整体尺寸,因此可以忽略位移尺寸对整体结构的影响．

在中国丧葬文化中，“丧”指对死者的哀悼礼仪，“葬”指对死者的处置方式。丧仪有着一整套从初终、入殓、居丧、出殡到做七的礼制细节，涉及报丧、吊祭、小敛、大敛、丧服、哭丧、卜时、下葬等许多习俗规定。葬礼包含葬制葬式、葬地坟高、葬器葬具、封树碑志等许多习俗规范，在接下来我们要探讨的汉魏六朝人的丧葬主张中都有涉及。

3.3.2 疲劳分析对比

2.6 影响颈动脉粥样硬化斑块形成的多因素logistic回归分析 以颈动脉粥样硬化斑块形成状况为应变量，建立多因素logistic回归模型。分析可得，吸烟史、糖代谢异常、高血压病史、TG、TC以及LDL-C均是影响颈动脉粥样硬化斑块形成的危险因素（均P＜0.05）。见表6。

运用HyperMesh作为前处理器对导缆孔进行网格划分,用软件Abaqus进行材料和属性定义以及边界和工况施加．采用HyperMesh 对导缆孔进行网格划分,划分后的数值模型包含C3D8和C3D20两种单元类型．导缆孔的材料是G230-450,弹性模量、泊松比和屈服极限分别为210 000 MPa、0.3 和230 MPa; 根据3种工况,依次分别在a、b、c方向施加压强,压强方向垂直于内孔曲面，压强取35 MPa,约束底部侧面所有节点的6个自由度(图2)．

  

图9 材料S-N曲线Fig.9 S-N curve diagram of material

选择疲劳损伤值作为导缆孔疲劳分析计算的参数，应力循环次数与损伤值的比值即为的疲劳寿命[11]．其中,应力强度准则选择最大主应力,存活率取为60%．

图10为导缆孔疲劳分析损伤值云图,通过50 000次循环变应力的作用,导缆孔损伤值最大处与导缆孔应力最大处一致,最大损伤值分为8.6×10-10、2.8×10-4、0.11．疲劳寿命=应力循环次数/损伤值,分别为5.8×1013、1.75×108、4.44×105．

  

图10 导缆孔疲劳分析损伤值云图Fig.10 Fatigue analytical damage graph of the chock

图11为优化后导缆孔疲劳分析损伤值云图．优化后的模型添加的所有分析条件参数均与之前的相同,通过50 000次循环变应力的作用,导缆孔损伤值最大处与导缆孔应力最大处基本一致(图10、图3),最大损伤值分别为3.38×0-10、1.07×10-5、3.5×10-4．疲劳寿命分别为1.48×1014、4.68×109、1.43×108,分别是之前疲劳寿命2.55、26.7、322.4倍．

  

图11 优化后导缆孔疲劳分析损伤值云图Fig.11 Fatigue analytical damage graph of optimized chock

表1为优化前后导缆孔性能参数，由表1可知,优化后的导缆孔工作可靠性得到较大提高,符合工程实际要求．

 

表1 导缆孔优化前后性能对比Table 1 Performance of chock beforeand after optimization

  

参数原始模型优化模型工况1最大应力/MPa217116工况2最大应力/MPa196157工况3最大应力/MPa244212工况1最大位移/mm0.04750.0847工况2最大位移/mm0.3820.142工况3最大位移/mm0.1150.227工况1疲劳寿命/次5.8×10131.48×1014工况2疲劳寿命/次1.75×1084.68×109工况3疲劳寿命/次4.44×1051.43×108质量/kg385242

4 结论

(1) 导缆孔经拓扑优化后,其形状为内、外双层薄壳,无加强筋,且质量较优化前减少37%．

2.资产负债率。受到债务契约的影响，管理者往往通过降低资产负债率来避免违约行为的发生。由于研发支出强度较大的公司大多为高新技术企业，为了能够迅速发展，这类企业往往需要保持较高的增长速度，相应就会增加流动负债，以应对其大规模研发投入；同时这类企业还要提高偿债能力，以获得较低融资成本，开拓更为宽阔的融资渠道。因此，为了降低资产负债率，提高利润水平，企业往往会将研发支出资本化。据此，我们提出第二个假设：

拓扑优化采用变密度法对导缆孔进行优化设计,其基本思想是人为引入一种假想的密度可变的材料,材料特性参数(如许用应力、弹性模量等) 与材料当量密度间的关系也是人为假定的[9]

(3) 3种工况下,优化后的导缆孔疲劳寿命,分别为优化前的2.55、26.7、322.4倍．

(4) 文中所用的分析方法有效可靠,可以为导缆孔结构设计和性能分析提供理论指导．

参考文献(References)

[ 1 ] IACS.Common structural rules for bulk carriers[S]. America:International Association of Classification Societies, 2012:19-23

D级供应商：该级别的供应商绩效评价指标各方面均低于行业水平，属于最低等级的供应商，应慎重选用该类供应商的产品或者对其进行淘汰管理等。

[ 2 ] 刘庆, 曲德选.巴拿马导缆孔的设计计算[J]. 造船技术, 2010(5):20-21.

LIU Qing,QU Dexuan. Design and calculation of the panama chock [J]. Marine Technology,2010(5):20-21. (in Chinese)

[ 3 ] 杨乃青, 豆正伟, 华晶. 浅析巴拿马运河对散货船系泊布置的影响[J]. 黑龙江科技信息, 2013(16):132-132.

YANG Naiqing, DOU Zhengwei, HUA Jing. Analysis on the influence of the panama canal on the mooring arrangement of bulk carriers[J]. Heilongjiang Science and Technology Information, 2013(16):132-132. (in Chinese)

LWT-LSSVM建模相比于其他建模方法的显著优势在于它所需要的训练数据量相对较少但精度很高，为提高运算速度，隔4 min取一组数据，共取48组数据进行训练建模。

[ 4 ] 倪伟平, 江晓玲. 船舶系泊属具标准梳理分析[J]. 船舶设计通讯, 2014(2):82-86.

NI Weiping, JIANG Xiaoling. Study of the standard of ship mooring and towing fittings [J]. Journal of Ship Design, 2014 (2): 82-86. (in Chinese)

[ 5 ] 穆安乐, 张玉龙, 由艳萍,等. 系泊参数对漂浮式风力机稳定性的影响规律研究[J]. 中国电机工程学报, 2015(1):151-158.

MU Anle, ZHANG Yulong,YOU Yanping, et al. Study on the effects of mooring parameters on the stability of floating wind turbine [J]. Proceedings of The Chinese Society for Electrical Engineering, 2015(1):151-158. (in Chinese)

[ 6 ] 倪伟平, 翟帅帅. 用于应急拖带的巴拿马导缆孔选型问题研究[J]. 船海工程, 2015(6):45-49.

历史学习有着一定的时间限制，由于事件发生的时间相对较远，学生对当时的社会状况缺少一定的认识，不能对历史进行充分理解。但是科学技术的发展为教育教学提供了多媒体教学设备，教师可以通过多媒体信息技术为学生提供相关的历史视频与音频资料。例如在学习中国抗战胜利，中华人民共和国成立的课程内容时，教师可以通过视频资料的播放，让学生身临其境地感受祖国的伟大胜利。

NI Weiping,ZHAI Shuaishuai.Selection of Panama canal chock for emergency towing [J]. Ship & Ocean Engineering, 2015(6):45-49. (in Chinese)

素养考查分析：这是一道数学应用题，综合考查了组合学、概率论、导函数、函数的单调性、函数的最值、数学期望等基本知识，以及运用组合公式、求函数的导函数，判断函数单调性、计算数学期望等基本技能．

[ 7 ] 广州造船厂. GB11569-1989. 巴拿马运河导缆孔[S].北京:国家技术监督局,1989:1-7.

[ 8 ] IMO,An example of emergency towing procedures[S]. London:International Maritime Organization, 2008:1-2.

[ 9 ] KANG B S, PARK G J, ARORA J S. Optimization of flexible multibody dynamic systems using the equivalent static load method[J]. AIAA Journal, 2005, 43(4): 846-852. DOI:10.2514/1.4294.

[10] JASION P, MAGNUCKI K. Elastic buckling of clothoidal-spherical shells under external pressure-theoretical study[J]. Thin-Walled Structures, 2015,86:18-23.DOI:10.1016/j.tws.2014.10.001.

式中，ω为系统电压的角速度；Id、Iq为电流的d、q轴分量；Ud、Uq为电压的d、q轴分量；L为滤波电感感值，R为变压器和逆变器的等效电阻。

[11] 张建，李金泰，唐文献，等.霍尔锚与强度和疲劳特性研究[J].江苏科技大学(自然科学版)，2017，31(3)：278-281.DOI：10.3969/j.issn.1673-4807.2017.03.004.

ZHANG Jian, LI Jintai, TANG Wenxian, et al. Investigation on the strength and fatigue of Hall anchor[J]. Journal of Jiangsu Unversity of Science and Technology(Natural Science Edition),2017,31(3):278-281.DOI：10.3969/j.issn.1673-4807.2017.03.004.(in Chinese)