计算机在输电线路导线、避雷线弧垂应力计算中的应用.doc

计算机在输电线路导线、避雷线弧垂应力计算中的应用摘 要架空输电导线、避雷线的弧垂，在实际工程中是非常重要的，这是因为架空线常年暴置于大自然中，受气象条件的影响较大，当温度变化或外部荷载大小变化时，架空线的长度会发生改变，虽然与其实际长度相比较其所占的百分比很小，但架空线长度的很小变化，会引起导线弧垂和应力产生相当大的变化。在线路设计中，若设计弧垂过小，架空线收得过紧，则会引起应力过大，振动加剧，同时杆塔荷载增大因而要求强度提高。反之，若弧垂过大，架空线较松弛，则为了满足导线对地的安全距离要求，需增加杆塔的高度和增大杆塔头部尺寸，使线路投资增加。因此，弧垂的理论研究和计算，已成为输电线路设计和施工的重要理论计算工作之一。特别是近几年来，我国超高压交、直流输电线路的迅速发展，越来越多地使用大载面导线或分裂导线，档距也越来越大，所以更应该进一步去研究弧垂计算理论，计算出更精确的弧垂数值来。本文从阐述在线路设计过程中的所涉及的架空线弧垂应力、不平衡张力及荷载计算等方面的理论、概念及计算方法等入手，其中重点是导线、避雷线弧垂应力计算方面。并通过结合计算机技术的应用进行算例，从中体现出计算机技术在线路设计繁杂的计算中应用的优越性。关键字：计算机，线路计算，应用目录引言11.导线、避雷线弧垂应力1.1 概述21.2 设计用气象条件21.3 导线、避雷线的比载计算21.4 架空线的最大使用应力和年平均应力41.5 临界档距41.5.1临界档距的概念. 51.5.2临界档距的计算公式.61.5.3临界档距的判定. 71.6 最大弧垂的判定.71.7 弧垂应力曲线计算及制作 81.7.1设计线路概况. 81.7.2导线弧垂应力曲线计算及制作. 91.7.3避雷线弧垂应力曲线计算及制作.171.8 弧垂应力计算语言、自设程序 222.避雷线最大使用应力的确定2.1 简述252.2 避雷线最大使用应力的选择253.导线、避雷线不平衡张力计算3.1 概述263.2 断线张力的确定263.3 线路正常运行中的不平衡张力计算273.4 避雷线的支持力274.杆塔荷载的计算4.1荷载的分类和荷载组合274.2荷载的计算条件284.3荷载计算用档距的确定294.4荷载的计算29结论30参考文献30附件302引 言架空线常年暴置于大自然中，经受各种恶劣气候的影响。为了保证其安全运行，必须使之在结构强度和电气性能等各方面能够很好地适应自然界的气象变化，以保证送电线路的安全运行，这有赖于线路设计环节，而设计环节中各项数据计算的准确性是保证基础，所以应该进一步去研究线路设计计算理论，计算出更精确的数值来。在架空输电线路的设计中，需要确定导线、避雷线在多种气象条件下的应力和弧垂，一般的计算方法是：1求出给定条件下的临界档距，判断控制架空线应力的档距范围；2求出耐张段的代表档距；3根据代表档距和临界档距确定导线、避雷线最大应力的控制条件；4利用架空线状态方程式求出给定气象条件下导线、避雷线的应力；5求出导线、避雷线在给定气象条件下的弧垂。但导线、避雷线应力孤垂的计算项目有时多达十种，即最大风速、覆冰情况、安装情况、事故断线、最低气温、最高气温、外过电压、内过电压、平均气温，这十种情况又对应十种气象条件，要进行大量反复的定值计算、比较和筛选，工作量很大；而传统的手工计算导线、避雷线弧垂应力的方法不仅计算量大、效率低,且容易出错。怎样把线路设计人员从繁杂的计算中解放出来，已成为许多专家学者和技术人员追求的目标。随着计算机技术的发展，计算机应用给各行各业的生产和管理带来了极大的方便，计算机的普及使计算机的应用渗透到社会生活的各个方面。同样的，在电力系统的线路设计、线路施工及线路运行和管理等方面计算机应用已发生巨大的作用。因此，计算机技术的迅速发展使实现这个目标有了技术支持，要解决这些问题，就要用到计算机将计算经常需用到一系列的数据资料，如设计手册或设计规范的数表、图表这些数据资料采用一定的方法形式存于计算机内或以数据库的形式存放于数据文件中，建立可修改的规则库，建立一定的逻辑关系，逻辑关系一旦建立，无论的其他参数如何变化，整定计算机都能自动完成弧垂应力计算、分析、保存和输出等工作，大大地提高弧垂计算精确度和设计效率,降低工作强度。本文从阐述在线路设计过程中的所涉及的架空线弧垂应力、不平衡张力及荷载计算等方面的理论、概念及计算方法入手，重点是导线、避雷线弧垂应力计算方面。其后，通过实际线路设计计算过程 (主要是导线、避雷线弧垂应力计算)中如何应用计算机技术，运用了自编程设计的弧垂应力计算通用程序和利用Microsoft公司推出的office套装软件中的Excel电子表格软件编辑进行计算两种方法，初步将现代新技术和传统方法进行了结合，从中体现出计算机技术应用在线路设计繁杂的计算中的优越性，以达到把线路设计人员从繁杂的计算中解放出来，提高计算精度、设计效率及降低人员劳动强度的目的，即本文的主要研究内容及应解决的问题。1 导线、避雷线弧垂应力1.1概述架空输电导线、避雷线的弧垂，在实际工程中是非常重要的，这是因为架空线常年暴置于大自然中，受气象条件的影响较大，当温度变化或外部荷载大小变化时，导线的长度会发生改变，虽然与导线的实际长度相比较其所占的百分比很小，但导线长度的很小变化，会引起导线弧垂和应力产生非常大的变化。在线路设计和施工中，若弧垂过小，导线收得过紧，则会引起应力过大，振动加剧，同时杆塔荷载增大因而要求强度提高。反之，若弧垂过大，导线较松弛，则为了满足导线对地的安全距离要求，需增加杆塔的高度和增大杆塔头部尺寸，使线路投资增加。因此，弧垂的理论研究和计算，已成为输电线路设计和施工的重要理论计算工作之一。由于架空输电线路的档距比其线径大得多，同时架空线又多采用多股细金属线铰合而成，所以架空线的刚性对其悬挂空间曲线形状影响很小，可忽略不计。可将导线看作一根理想的、柔软的，而且荷载是均匀分布的线段。故悬挂在两基杆塔间的架空线是悬链线形状。积分通式为：（1-1）式中C1，C2均为积分常数，其值取决于坐标系原点位置的选择； 0导线各点的水平应力（即最低点的应力），N/mm2； g导线的比载，MPa/m。1.2 设计用气象条件所谓设计用气象条件，是指那些与线路的电气强度和机械强度有关的气象参数，如风速、温度、湿度、和雷电情况等。而对架空线机械强度设计起主导作用是风速、温度和覆冰厚度这三个气象参数，所以在设计过程中，必须对沿线的气象情况进行全面的了解。风速、大气温度和覆冰厚度是架空线路设计中要考虑的主要因素，它们的取值方法就是线路设计中气象条件的组合。气象条件的组合，除了应合理反映自然规律的变化外，还要考虑技术经济上的合理性及设计计算上的方便性。高压架空输电线路，计算用气象条件组合主要有：线路正常运行情况下的气象条件组合、线路安装和检修情况下的气象条件组合、外过电压时的气象条件组合、线路断线事故情况下的气象条件组合、线路耐振计算用的气象条件组合、内过电压的气象条件组合。1.3导线、避雷线的比载计算1.3.1导线、避雷线的比载种类和计算公式。架空线的比载是单位长度架空线上所受的荷载折算到单位截面积上的数值，其常用单位是N/（m.mm2）或MPa/m。作用在架空线上的荷载主要有自重、冰重、风荷载以及覆冰厚的风荷载等，相应的有自重比载、冰重比载、风压比载和覆冰比载等。根据作用方向的不同，比载又可分为垂直比载（自重、冰重比载）、水平比载（无冰风压、覆冰风压比载）和综合比载。下列是各种比载的计算公式，为了清楚表示，用符号gx（b，v），其中x=17，表示七种不同的比载，b表示计算比载时覆冰的厚度，mm：v表示计算风速m/s。（1）自重比载g1（0.0）是架空线自身的质量引起的比载，其大小可认为不受气象条件变化的影响。其计算公式为（1-2）式中 G0架空线单位长度质量，kg/km； A架空线的截面积，mm2。（2）冰重比载g2（.0）覆冰时的冰重由架空线承受。其计算公式为（1-3）式中 架空线的计算直径，mm。（3）垂直总比载g3（b.0）是架空线自重比载与冰重比载之和。其计算公式为（1-4）（4）导线无冰时的风压比载g4（0.v）。架空线的风压是由作用于架空线的空气动能所引起的。其计算公式为 （1-5）式中 风速不均匀系数，按设计风速不同取值； C风载体型系数：线径d17mm,取1.2；线径17mm,取1.1;覆冰时无论线径大小,取1.2。 风向与线路方向夹角（度）。（5）导线覆冰时的风压比载g5（b.v）。受风方向的投影面积应为导线的直径d及覆冰厚度2b之和，并且体型系数c均取1.2。其计算公式为（1-6）（6）导线无冰有风时的综合比载g6（0.V）是架空线自重比载与水平风压比载的几荷和。其计算公式为（1-7）（7）导线覆冰有风时的综合比载g7（b.V）是架空线垂直总比载与覆冰时风压比载的几荷和。其计算公式为（1-8）1.3.2应用Microsoft Excel进行比载计算架空线比载计算是一项繁琐的计算，要进行十多种气象情况下的比载计算，如风压比载g4（0，v）及无冰有风综合比载g6（0，v）计算要分别计算外过电压(安装有风)、内过电压、最大风速三种气象条件下的比载。人工手算费时且容易出错，而比载算数的正确是弧垂应力算数正确的基础。Microsoft公司推出的office套装软件是应用十分广泛的办公软件，在国内拥有很多用户。现运用计算机，应用Microsoft Excel电子表格软件编辑进行比载的计算，由于Excel电子表格具有可编程、即时计算、递归的功能，因此，通过输入架空线的截面积、重量、直径、覆冰厚度、风速等，再定制出比载计算公式，从而完成比载的自动计算，使原本繁琐的计算变得简单、快捷、准确。1.3.2.1应用Microsoft Excel进行比载计算公式编辑及原理(1)在进行计算前先确定设计线路所采用架空线的规格、物理性能及气象区等数据。(2)打开Microsoft Excel电子表格，根据比载计算涉及项制定表格格式，具体见表1-1:表1-1 架空线比载计算电子表格(3) 单元格D3K13是用来输入比载计算类别、线路气象条件、采用架空线规格及物理性能等的数据，其中风速不均匀系数的取值对应设计风速、风载体型系数C(线径g1所以垂直比载愈大，与地面垂直的弧垂也愈大，但此时31，则又将使弧垂f减少，由以上两种气象情况分析，很难判断最大弧垂究竟是在最高气温或最大垂直比载时出现。为了求得最大弧垂，可以利用状态方程式分别求得两种气象条件下的应力，然后再运用弧垂公式计算各自的弧垂，加以比较而得到。也可以先判定出现最大弧垂的是哪一种气象条件，再计算该种气象条件下架空线的弧垂即为最大弧垂。常用的判定方法有临界温度法和临界比载法两种。(1)临界温度判定法若在某一温度下，架空线在自重比载(最高气温时的比载)作用下产生的弧垂与覆冰无风(最大垂直比载)时产生的弧垂相等，则此温度称为临界温度。以覆冰无风为第一种状态，临界温度为第二种状态，从状态方程式中可得临界温度tk计算公式为 ( ) (1-14)式中 ：tK临界温度，；tb架空线在最大垂直比载时温度，；3架空线在最大垂直比载时的应力，N/mm2。g3架空线的最大垂直比载，MPa/m；g1架空线的自重比载，MPa/m。将计算出的临界温度tb与最高温度tmax相比较，温度高者为出现最大弧垂的控制条件。若tbtmax，则最大弧垂发生在最大垂直比载气象条件，反之最大弧垂发生在最高气温气象条件。(2)临界比载判定法若架空线在气温tb下，某一垂直比载使其产生的弧垂于最高气温条件下的弧垂相等，则此比载称为临界比载。以最高气温为第一种状态，临界比载为第二种状态，从状态方程式中可得冰层临界比载g2k计算公式为 (MPa/m) (1-15)式中：tmax架空线的最高气温，； 对应最高气温时的应力，N/mm2。 其它符号的意义与式（1-23）相同。当计算的冰层临界比载g2k大于某一气象区的冰层比载g2时，导线最大弧垂出现在最高温度，反之最大弧垂出现在覆冰无风时。1.7 架空线弧垂应力曲线的计算和制作通过计算架空线应力弧垂后制成的应力弧垂曲线图是输电线路设计图纸中的一部分，是用坐标纸以横坐标为代表档距纵坐标为应力，根据状态方程式算出的应力，按合适比例作出=f()之间的函数关系曲线。它反映了某地区各种气象情况下应力、弧垂与代表档距之间的变化关系，计算它们的主要目的是为了提供导线、避雷线的力学特性，从而能进行线路杆塔设计的进一步计算。前面已把架空线应力弧垂计算的基本理论和公式应用阐述清楚，下面将通过实际线路设计计算过程中的导线、避雷线弧垂应力计算过程怎样应用计算机技术，来进一步从实际上说明在线路导线、避雷线弧垂应力计算中应用计算机的目的是把线路工程设计人员从繁杂的计算中解放出来和提高计算精度、设计效率及降低人员的劳动强度。1.7.1 线路概述(1)线路基本情况：线路穿越的地区属于西北地区第三气象区，设计所采用的技术参数为：西北级气象区，最大设计风速取30m/s。导线采用LGJ-95/20、避雷线采用GJ-35。(2)气象条件：表1-2 气象组合条件的相关数据表(西北级)计算情况最高气温最低气温最大风速覆冰年平均气温外过电压(无风)外过电压(有风)内过电压安装有风事故情况温度40-30-5-5515155-150风速（m/s）003010001015100冰厚（mm）0005000000(3)线路采用的架空线规格及物理性能见表1-3、1-4：表1-3 LGJ-95/20规格及物理性能表总截面积A（mm2）直径d（mm）线膨胀系数（1/）弹性模数（mm2/N）单位长度重量G（kg/km）瞬时破坏应力p（N/mm2）安全系数K113.9613.8718.510-613.15710-6408.9326.42.5表1-4 GJ-35规格及物理性能表总截面积A（mm2）直径d（mm）线膨胀系数（1/）弹性模数（mm2/N）单位长度重量G（kg/km）瞬时破坏应力p（N/mm2）安全系数K37.157.811.510-65.51110-6318.21224.013.21.7.2 LGJ-95/20导线弧垂应力曲线的计算和制作1.7.2.1 应用表1-1，按对应单元格项输入LGJ-95/20导线的截面积、重量、直径及线路所经过气象区的覆冰厚度、风速等数据，其中风速不均匀系数对应设计风速取值、风载体型系数C取1.2；自动生成的比载计算结果见表1-5： 表1-5 LGJ-95/20导线的比载计算结果从上的计算中可以看出，定制好Excel电子表格后，只要相对应地输入相关数据，结果“L”栏中就会自动进行比载的计算并显示结果，避免了手工计算的繁琐使计算变得简单、快捷、准确，减轻了设计人员的劳动强度。1.7.2.2 计算LGJ-95/20导线最大使用应力及年平均运行应力：最大使用应力： （）年平均使用应力： （）1.7.2.3 控制LGJ-95/20导线应力的档距值计算及判断，计算的相关数据见表1-6：表1-6判断控制导线LGJ-95/20应力的档距范围的相关数据表项 目最大风速年平均气温最低气温最大覆冰最大使用应力N/mm2130.681.6130.6130.6比载10-3（MPa/m）69.87535.16335.16360.109气 温（）-55-30-5g/ 10-45.3504.3092.6924.603顺 序 编 号DBAC计算导线的临界档距值：将表1-6中数据分别代入公式（1-13）计算各控制气象下导线的临界档距值：;；。有效临界档距判别如表1-7所示，判别结果如图1-8所示： A控制 B控制 D控制24.2 330.98 图1-8 有效临界档距判别结果表1-7 有效临界档距判别表（m）ABC=24.2=648.7=0=282=330.98=227.7即 LGJ-95/20导线的临界档距： 当24.2m时，控制气象条件为最低气温=130.6Mpa; 当24.2m330.98m, 控制气象条件为年平均气温=81.6Mpa; 当330.98时，控制气象条件为最大风速=130.6Mpa;1.7.2.4 最大弧垂出现气象情况的判定由(1-15)公式 当档距0时，根据公式（1-11）最高气温时导线的应力为： Mpam又22.941 Mpam=68.47=22.941所以导线最大弧垂发生在最高气温气象条件下。1.7.2.5 各种气象情况下LGJ-95/20导线应力和弧垂值的计算在架空线状态方程式(1-10)中令 则公式(1-10)可简化为 (1-25)对于上式，在设计给定参数后，则A、E均为已知量，解此方程可用牛顿迭代法求解，具体过程(略)：由公式(1-25)式中求出了该方程的通式为 (1-26)解上式时，先任意假定一个代入，求出值来，判断是否成立(为设定迭代精度，取10-6)。若不成立，再将代入，求出值来，判断是否成立。如此循环，直到求得，且，则则为公式(1-26)的解。求出后，代入公式： (1-27)则可以求出与之对应的弧垂值了。由于弧垂应力计算过程繁琐，稍不小心就会出错。现继续应用office套装软件中的Excel电子表格软件编辑进行计算，因为Excel电子表格具有可编程、即时计算、递归的功能，而且Microsoft公司推出的office套装软件是应用十分广泛的办公软件，在国内拥有很多用户，技术大众化。前面比载计算时已应用Excel电子表格软件编辑进行计算，领略了其完成繁琐的计算的简单、快捷及准确性且其特点是随着设计数据的变化而自动进行计算，结果一目了然，因此现继续利用其来计算导线的应力和弧垂，编辑原理跟比载计算一样，只不过表格的行、列数根据内容的不同而有所区别，其计算过程及结果如下：一、A值的计算：插入Excel电子表格，按计算要求建立表格(1)在单元格G3G11进行编辑公式工作。以单元格G3为例说明编辑公式步骤：单击单元格G3，使之成为当前编辑对象。单击编辑栏上的“=”，打开公式编辑对话框。根据A值的计算通式，输入“C32/24/D3”。当公式编辑对话框下显示出计算结果后，单击“确定”，退出编辑状态，完成公式编辑。(2) 在对应气象情况单元格中输入 (表1-5)的比载值g和(表1-3)LGJ-95/20导线的弹性模数，相应的结果就会自动显示在“G”栏，表格格式和计算结果见表1-7：二、E值的计算：插入Excel电子表格，按要求建立表格后(1)在单元格I4I12、J4J12进行编辑公式工作。以单元格J4为例说明编辑公式步骤：单击单元格J4，使之成为当前编辑对象。单击编辑栏上的“=”，打开公式编辑对话框，根据其计算公式变形，输入“D4-E4/F4*(G4-H4)”，当公式编辑对话框下显示出计算结果后，单击“确定”，退出编辑状态，完成公式编辑。(2) 根据公式及，单元格I4I12、J4J12各对应的计算公式如下：I4=C42/24/F4； J4 =D4-E4/F4*(G4-H4)；I5=C52/24/F5； J5 = D5-E5/F5*(G5-H5)；I6=C62/24/F6； J6 = D6-E6/F6*(G6-H6)；I7=C72/24/F7； J7 = D7-E7/F7*(G7-H7)；I8=C82/24/F8； J8 = D8-E8/F8*(G8-H8)；I9=C92/24/F9； J9 = D9-E9/F9*(G9-H9)；I10=C102/24/F10； J10 = D10-E10/F10*(G10-H10)；I11=C112/24/F11； J11 = D11-E11/F11*(G11-H11)；I12=C122/24/F12； J12 = D12-E12/F12*(G12-H12)。(3) 编辑公式工作完成后，本表格即可用于E值计算，表中各参数的取值是根据判断出控制导线应力的档距范围的气象条件来定的。在本设计线路：当在292.93(m)的范围内由年平均气温B(t=20、cp=71.05 N/mm2、g=34.01510-3MPa/m）控制；在292.93(m)的范围内由最大风速A(t=10、max=113.68N/mm2、g=65.17310-3MPa/m）控制。故分两种情况进行计算，表格格式及计算结果见表1-8、1-9：表1-7 LGJ-95/20导线A值的计算表表1-8 LGJ-95/20导线E值的计算表1 表1-9 LGJ-95/20导线E值的计算表2表1-10 LGJ-95/20导线E值的计算表3 三、各种气象情况下LGJ-95/20导线应力、弧垂f的计算：点击插入、建立Excel电子表格并按对应单元格项输入相应线路参数，其表格格式及计算结果见表1-10，具体编辑方法、应用过程说明如下：(1)建立电子表格表1-10后，在单元格D3F119、I3P119进行编辑公式工作。以单元格P3为例说明编辑公式步骤：单击单元格P4，使之成为当前编辑对象。单击编辑栏上的“=”，打开公式编辑对话框，根据公式（1-27），对应输入“C2*G/8/O”，当公式编辑对话框下显示出计算结果后，单击“确定”，退出编辑状态，完成了弧垂“f”公式编辑，即可进行弧垂计算。(2) 应力 “”是利用牛顿迭代法求，在单元格“I3O119”中对应公式(1-26)编辑好各应力计算公式，其计算过程：假定一个代入，求出值来，判断是否成立(为设定迭代精度，取10-6)。若不成立，再将代入，求出值来，判断是否成立。如此循环，直到求得，且，则为所求导线的应力，取=n，完成求解(其中表中隐藏了“”的求解过程，只要点击打开电子表格就可以了解全部过程)。(3)单元格中对应公式说明(以I2、P3为例，其它单元格对应公式是同样，只是对应单元不同)：对应计算公式如下：I2=(2*H22-$E2*H2+$C2/H2)/(3*H2-2*$E2)； P3= =B22*F2/(8*1000*N2)；其余单元格的数据，按鼠标左鍵不动托动当前单元格的右下角移动就可以得到对应档距的应力和弧垂。(4) 编辑公式工作完成后，本表格即可用于架空线应力弧垂计算。 表 1-10 各种气象情况下 LGJ-95/20导线的弧垂应力计算表气象情况(）ADEg1(10-3MPa/m)max1(Mpa)(Mpa)(m)最 高 气 温 4000.00.00032.173 35.163 130.60 100.462 32.181 0.0000 597.90.00632.167 35.163 130.60 100.462 32.269 0.0034 10391.60.02332.150 35.163 130.60 100.462 32.527 0.0135 15881.10.05232.121 35.163 130.60 100.462 32.938 0.0300 201566.40.09232.081 35.163 130.60 100.462 33.482 0.0525 24.22293.40.13532.038 35.163 130.60 100.462 34.022 0.0757 509790.01.47030.917 35.163 130.60 100.462 37.777 0.2909 10039160.05.88026.507 35.163 130.60 100.462 45.458 0.9669 15088110.013.23019.157 35.163 130.60 100.462 51.886 1.9060 200156640.023.5208.867 35.163 130.60 100.462 57.029 3.0829 250244750.036.750-4.363 35.163 130.60 100.462 61.131 4.4938 300352440.052.920-20.533 35.163 130.60 100.462 64.413 6.1414 330.98428989.064.414-32.027 35.163 130.60 100.462 66.114 7.2829 350479710.0111.108-43.782 35.163 130.60 100.462 66.080 8.1483 400626560.0145.120-77.794 35.163 130.60 100.462 66.008 10.6541 450792990.0183.668-116.342 35.163 130.60 100.462 65.955 13.4951 500979000.0226.750-159.424 35.163 130.60 100.462 65.914 16.6710 6001409760.0326.520-259.194 35.163 130.60 100.462 65.856 24.0271 6501654510.0383.208-315.882 35.163 130.60 100.462 65.836 28.2071 最 低 气 温 -3000.00.000130.600 35.163 130.60 100.462 130.600 0.0000 597.90.006130.594 35.163 130.60 100.462 130.600 0.0008 10391.60.023130.577 35.163 130.60 100.462 130.600 0.0034 15881.10.052130.548 35.163 130.60 100.462 130.600 0.0076 201566.40.092130.508 35.163 130.60 100.462 130.600 0.0135 24.22293.40.135130.465 35.163 130.60 100.462 130.600 0.0197 509790.01.470129.343 35.163 130.60 100.462 129.923 0.0846 10039160.05.880124.933 35.163 130.60 100.462 127.348 0.3451 15088110.013.230117.583 35.163 130.60 100.462 123.372 0.8016 200156640.023.520107.293 35.163 130.60 100.462 118.456 1.4842 250244750.036.75094.063 35.163 130.60 100.462 113.172 2.4274 300352440.052.92077.893 35.163 130.60 100.462 108.070 3.6604 330.98428989.064.41466.399 35.163 130.60 100.462 105.178 4.5780 350479710.0111.10854.645 35.163 130.60 100.462 101.348 5.3127 400626560.0145.12020.632 35.163 130.60 100.462 93.029 7.5595 450792990.0183.668-17.916 35.163 130.60 100.462 86.957 10.2357 500979000.0226.750-60.998 35.163 130.60 100.462 82.576 13.3071 6001409760.0326.520-160.768 35.163 130.60 100.462 77.001 20.5496 6501654510