1Methodofcalculation-国际船舶网.doc

1、锚链疲劳载荷计算方法 用有限元分析软件ANSYS进行静载荷, 计算模型选用实体单元. 因为单个锚链实体和所受到的载荷是关于X轴方向对称, 所以取半个锚链进行计算、 实体被划分为75421四角形单元(如图1，取平行于锚链长轴方向为X方向)。 根据计算需要对称体在截面处X方向加约束，并且截面上的一个节点在X-，Y-，Z-方向均加约束。载荷看作为均布载荷在直径为 25.52 mm的面上。 半径计算为： R==12.76mm 其中：E1=E2=210000N/mm2=2.1×1011Pa = =0.3 = + = 则=147/4； 仅做线性计算，既不考虑材料非线性也不考虑几何非线性计算。 疲劳寿命计算基于 Miner线性累积疲劳准则（即曼耐尔疲劳损伤积累假说）计算的指数方程为： 其中： N—循环次数 ---应力变化范围 m---S/N曲线关于对数横坐标下的斜率 ，m 取为3.36 a-----方程的实常数, 由实验获得, loga （以10为底a的对数）取为13.6451 寿命损伤率： 其中：ni----各 对应的积累循环次数（实际情况下的） Ni-----各 对应的材料发生疲劳破坏时极限循环次数 据Miner疲劳准则材料达到疲劳失效时，总疲劳寿命损伤率： =1 意味着, 如果 D1,已经失效, 如果 D1,失效还没发生（还没达到疲劳寿命极限）. 取0.3倍的最小破断载荷作为平均拉力, 12421×0.3=3726 KN. 取 372650 KN(或 3676 KN和3776 KN)作为计算 S/L值的载荷范围 ，用在3676 KN 和3776 KN下的锚链最大主应力的差值计算得到应力-载荷比值S/L, S/L==0.080508815 S/L乘以载荷变化范围得到应力变化范围, =×S/L。由得到各的值。由得到疲劳寿命损伤率。其中载荷加载频谱（实际的）由镇江锚链厂提供, 可以得到每年的累积损伤率 和200年的总的累积疲劳损伤率。 计算结果 静力有限元分析 图2.1—2.2 表示在试验载荷（proof-load 8695KN）下的体应力。 图3.1-3.2表示在最小断裂载荷（minimum breaking load 12421KN）下的体应力。 . 图 4表示在单位载荷（unit load 100KN）下的应力 在锚链的加载区域加载8695 KN 时最大的体应力26300 Mpa. 除了加载区域之外，最大的体应力在 圆弧的连接处的截面上为868Mpa. 在锚链的加载区域加载12421 KN 时最大的体应力37500Mpa. 除了加载区域之外，最大的体应力在 圆弧的连接处的截面上为1270Mpa. 这些结果表示在试验载荷（proof-load 8695KN）和最小断裂载荷（minimum breaking load 12421KN）下的体应力都超过材料的强度极限860 Mpa. 但是这些是线性计算的结果。事实上，高的应力将随着局部塑性变形而减小。根据同种材料相似结构的非线性分析，在最小破断载荷下，只有局部塑性变形而不会发生失效。 塑性变形的区域比 R3 钢的锚链大。 疲劳寿命分析 基于二次开发，计算疲劳极限值，并计算各应力水平的疲劳损伤度，得到累积疲劳损伤率小于1，从表2中可以看出，200年的累积疲劳损伤率小于1，因此，没有疲劳失效发生。 结论 在这篇报告中，锚链环 (SE)的应力和疲劳寿命已经被分析和计算出，锚链环可以承受最小断裂载荷（ minimum breaking load.）。 前面提及的锚链的疲劳寿命远远超过200年。 图1实体划分网格 (取平行于锚链长轴方向为X方向根据计算需要对称体在截面处X方向加约束，并且截面上的一个节点在X-，Y-，Z-方向均加约束) 图2.1 Proof load 8695 KN 图2.2 Proof load 8695 KN 图3.1 Minimum Breaking Load 12421 KN 图3.1 Minimum Breaking Load 12421 KN 图 4 Unit Load (100 KN) Table 1 Log(a) 13.6451 m 3.36 应力/载荷（Stress/Load S/L） 0.0805088(Mpa/KN) 设计寿命（Designed life ） 20(Year) 安全系数（Safety factor） 10 Table 2 载荷范围Load Range (KN) 疲劳