kw四冲程汽油机曲轴设计电子.doc

学 号： 0120707210133题 目100kw四冲程汽油机曲轴设计学 院汽车工程学院专 业热能与动力工程班 级0701姓 名刘玲玲指导教师罗马吉2010年11月23日目 录0前言-41汽油机的结构参数-4 1.1初始条件-4 1.2发动机类型-4 1.2.1冲程数的选择-4 1.2.2冷却方式-4 1.2.3汽缸数与气缸布置方式-4 1.3基本参数-4 1.3.1形成缸径比的选择-4 1.3.2气缸工作容积V，缸径D的选择-42热力学计算-6 2.1热力循环基本参数的确定-6 2.2P-V图的绘制-6 2.3P-V图的调整-73运动学计算-8 3.1曲柄连杆机构的类型-8 3.2连杆比的选择-8 3.3活塞运动规律-83.4连杆运动规律-104动力学计算-10 4.1质量转换-104.2作用在活塞上的力 -11 4.2.1作用在活塞销上的力-11 4.2.2惯性力-144.3输出合成转矩-155曲轴零件的结构设计-15 5.1曲轴的工作条件，结构形式和材料的选择-15 5.1.1曲轴的工作条件和设计要求-15 5.1.2曲轴的结构形式-16 5.1.3曲轴的材料-16 5.2曲轴主要尺寸的确定和结构设计细节-16 5.2.1主要尺寸-16 5.2.2一些细节设计-175.2.2.1油道布置-175.2.2.2曲轴两端的结构-175.2.2.3曲轴的止推-185.2.2.4过渡圆角-185.2.2.5平衡分析-186曲轴强度的校核-196.1曲轴的弯曲强度校核-196.2曲轴的扭转强度校核-217小结-238参考文献-24附录-25附表1-25附录2-25附录3-29100kw四行程汽油机曲轴设计0前言曲轴是内燃机的关键零件，其尺寸大小在很大程度上影响着内燃机的外形尺寸和重量以及工作可靠和寿命。曲轴的破坏事故可能引起其他零件的严重损坏。为了更深入地了解曲轴乃至发动机，学校为我们安排了汽车发动机设计课程设计。本设计的主要任务是100kw四行程汽油机曲轴的设计。1. 汽油机的结构参数11初始条件额定功率：P100KW 平均有效压力:=0.81.2MPa活塞平均速度：Vm18ms12发动机类型1.2.1冲程数的选择四冲程（=4）1.2.2冷却方式水冷 1.2.3气缸数与气缸布置方式直列式四缸机（i=4）1.3基本参数1.3.1行程缸径比SD的选择 初步选择行程缸径比为SD=1.01.3.2气缸工作容积V，缸径D的选择根据内燃机学的基本计算公式：P=Pme*Vs*i*n（30）Vm=S*n30000Vs=D*S4000000 其中Pe为发动机的有效功率，依题100kw Pme为发动机的平均有效压力，依题为 Vs 为气缸的工作容积，单位为L i为发动机的气缸数目 ，依题为4 n 为发动机的转速，单位为r/min Vm为活塞的平均速度，依题取为16m/s S为发动机活塞的行程,单位为mm D为发动机汽缸直径,单位为mm为发动机的行程数，依题为4根据以上的条件代入以上公式得：D=S=89.2mm , n=5380.8 r/min , Vs=0.557L将额定转速n取整为5200r/minVs=0.577LD=92.3mm再由S/D=1.0计算并同D一样圆整得S=92mm，D=92mm根据计算得出的D和S计算实际的单缸排量、行程缸径比、平均有效压力和活塞平均速度得Vs=0.61LS/D=1.0Pme=0.946MPaVm=15.95m/s利用以上算出的实际Pme、Vm值与给定值比较，他们应与同用途类型先进发动机水平相当。本设计曲轴一些基本参数见表一。2热力学计算 汽油机利用化油器对供给气缸的混合气进行定量调节，以适应汽油机负荷变化的需要。汽油机的工作过程包括进气、压缩、做功和排气四个过程。在本设计过程中，先确定热力循环基本参数然后重点针对压缩和膨胀过程进行计算，绘制P-V图并校核。2.1热力循环基本参数的确定 根据参考文献【汽车拖拉机发动机】压缩过程绝热指数n=1.321.38，初步取n=1.35 膨胀过程绝热指数n=1.231.28，初步取n=1.25 根据参考文献【内燃机学】汽油机压缩比=612，初取=9 根据参考文献【发动机设计】p=69，初取p=62.2P-V图的绘制通常情况下，压缩始点的压强在Pa=（0.80.9）（为当地大气压力值），假定外界=0.10MPa，选定Pa=0.09 MPa，将压缩过程近似看作绝热过程，由上一步n=1.35，并利用PV=const，利用前面计算所得单缸容积，可以在excel中绘出压缩过程线。混合气体在气缸中压缩后，经等容加热，利用值可得最大爆发压力值，膨胀过程类似于压缩过程，n由上一步取得1.25，绘出膨胀线。最后连接膨胀终点和压缩始点。得出理论的P-V图。简化的条件为：假设工质是理想气体，其物理常数与标准状态下的空气物理常数相同。假设工质是在闭口系统中作封闭循环。假设工质的压缩及膨胀是绝热等熵过程。假设燃烧过程为等容加热过程，工质放热为定容放热。Pa=Pc, Pa=0.09MPa, Va=Vc+Vs, Vs=0.61L, =9,得Va=0.69L,Vc=0.076L Pc=1.75MPaPz=p* Pc=12.23MPa。Pz=Pb，Pb=0.785MPa。理论的P-V图如下图所示： 图12.3 P-V图的调整 实际的P-V图和利用多变过程状态方程绘制的P-V图还存在一些差别，主要是点火提前角和配气相位的原因。对上图作以下调整：最大爆发压力：P取理论水平的2/3，具体值为8.15MPa8.5MPa，以此值与原图形相交，水平线以上的部分去掉，余下部分作些调整。考虑到实际过程与理论过程的差异，最大爆发压力发生在上止点之后1215，选择最高爆发压力出现在上止点后12.点火提前角：据资料得常用的范围是2030，经调整后取26。排气提前角：常使用的范围是4080，经调整后取60。调整后的P-V图如下： 图2由热力学计算所绘制的示功图为理论循环的示功图，其围成的面积表示的是汽油机所做的指示功数值由对示功图积分后求得的面积来表示。由理论P-V图可求解出该示功图下的指示功的大小。经过计算后得： =N*Pm*0.02*1000=688JN为理论P-V图中的格子数，由图知为86，Pm为每小格纵坐标表示的压强值为0.4MPa，0.02L表示每小格横坐标表示的容积值。所以汽油机的平均有效压力： =*/=688*0.85/0.00061=0.959MPa代入公式得P=101.38KW其中为汽油机的机械效率，与前面计算的结果大致一样，在4%以内，故上面选取的参数和以后的相关计算在满足制造的同时能够前后一致。3.运动学计算3.1曲柄连杆机构的类型在往复活塞式内燃机中基本上采用三种曲柄连杆机构：中心曲柄连杆机构，偏心曲柄连杆机构和关节曲柄连杆机构。其中中心曲柄连杆机构应用最广泛。本次设计选择中心曲柄连杆机构。3.2连杆比的选择据【发动机设计】知，=1/31/5，车用发动机多用小连杆，初选=1/4.3.3活塞运动规律 活塞位移：=0.25，为曲轴半径，r=S/2=92/2=46mm 经计算后X-图如下： 图3活塞速度：V=r（Sin+Sin2/2） =n/30 n=5200r/min r=0.046m经计算后V-图如下： 图4活塞加速度：a=r（Cos+ Cos2）经计算后a-图如下： 图53.4连杆运动规律 连杆式做复合平面运动，即其运动是由随活塞的往复运动以及绕活塞销的摆动合成。连杆相对于气缸中心的摆角为=arcsin（sin）经计算后-图如下： 图64.动力学计算4.1 质量转换沿气缸轴线作直线运动的活塞组零件，可以按质量不变的原则简单相加，并集中在活塞销中心。mp=mpi匀速旋转的曲拐质量，可以按产生离心力不变的原则换算，并集中在曲柄销的中心mc=miri/rmi曲拐单元的质量，ri各单元的旋转半径 做平面运动的连杆组，根据动力学等效性的质量，质心和转动惯量守恒三原则进行质量换算。3个条件决定三个未知数，可用位于比较方便的位置上即连杆小头，大头和质心处三个质量来代替连杆。实际结果表明m与m、m相比很小，为简化受力分析，常用集中在连杆小头和大头的2个质量代替连杆。 0.253kg=0.202kg往复质量：m= + m=1.8=0.455kg旋转质量：mr=m+ m=2.276kg其中D=92mm，l为活塞厚度l=8mm，活塞材料为共晶铝合金：=2.7g/cm3，H为活塞高度H=（0.81.0）D=83.7mm4.2作用在活塞上的力作用在活塞销中心的力，是P和的合力，P为气体作用力，为往复惯性力4.2.1作用在活塞销上的力气体作用力在曲柄连杆中的传递情况如图7，F对连杆产生拉力，还对气缸产生侧向力F=1/4D（P-P） P活塞顶上的压力，-活塞背压 侧压力 =Ptan 连杆力=P/cos 切向力= sin(+)=P 径向力 =cos(+)=P 综合应用上述等式，求解出各力随曲轴转角的变化曲线。把计算过程的相关数据收录在附录3中。引入的角是连杆和气缸中心线方向的夹角，它和曲轴转角满足下列关系式：lsin=rsin 图7 =-0.491() 图8 P= + P 图9 =Ptan 图10 =P/cos 图11=P 图12 =P 图134.2.2惯性力往复惯性力Fj在机构中的传递情况与Fg很相似，Fj也使机构受负荷，也产生转矩和倾覆力矩，由于Fj对汽缸盖没有作用，所以它不能在机内自行抵消，是向外表现的力，需要由支承承受。 图14旋转惯性力Fr=mrr2，当曲轴角速度不变时，Fr大小不变，其方向总是沿着曲轴半径向外。如果不用结构措施（如平衡块）消除，它也是自由力，使曲轴轴承和内燃机承受支反力，它不产生转矩和倾覆力矩。4.3合成扭矩 在讨论输出合成扭矩时，要把气体力和往复惯性力合成。虽然往复惯性力对输出转矩的平均值没有影响，但对输出转矩瞬时值的影响很大。多缸内燃机的总转矩等于不同相位的各缸转矩的叠加，对于发火间隔角均匀的内燃机来说，总转矩曲线是将各缸转矩曲线错开一个相当于发火间隔角（=720/i）的距离，然后叠加的结果，其变化周期就是。5.曲轴零件的结构的设计5.1曲轴的工作条件，结构形式和材料的选择5.1.1曲轴的工作条件和设计要求 曲轴的基本工作载荷是弯曲载荷和扭转载荷；对内不平衡的发动机，曲轴还承受内弯矩和剪力；未采取扭转振动减振措施的曲轴还可能作用着幅值较大的扭转振动弹性力矩。这些载荷都是交变性的，可能引起曲轴疲劳失效。实践表明，弯曲载荷具有决定性作用，弯曲疲劳失效是主要的破坏形式。因此，曲轴结构强度的研究重点是弯曲疲劳强度，曲轴设计上要致力于提高曲轴的疲劳强度。 曲轴形状复杂应力集中相当严重，特别是在连杆轴颈与曲柄臂的过渡圆角处和润滑油孔出口附近的应力集中尤为突出。曲轴的设计制造要特别注意，设法缓和应力集中现象，强化应力集中的部位。曲轴各轴颈在很高的比压下，以很大的相对速度在轴承中发生滑动摩擦。为了提高曲轴的使用寿命，应使各摩擦表面耐磨，并匹配好适当条件的轴瓦。曲轴的刚度亦很重要，所以设计曲轴时应保证曲轴有尽可能高的弯曲刚度和扭转刚度。曲轴是内燃机的基础件，其主要结构的尺寸选择，还要照顾到内燃机强化，适当留出发展余地。5.1.2曲轴的结构形式曲轴从整体结构上看可以分为整体式和组合式，从支承方式看，曲轴有全支承结构和浮动支承结构，为了提高曲轴的弯曲强度和刚度，现代多缸内燃机的曲轴都采用全支承结构。5.1.3曲轴的材料 曲轴材料一般使用45，40Cr，35Mn2等中碳钢和中碳合金钢。轴颈表面经高频淬火或氮化处理，最后进行精加工。目前球磨铸铁由于性能优良，加工方便，价格便宜广泛的用于曲轴材料。本设计采用QT800.5.2曲轴主要尺寸的确定和结构设计细节5.2.1主要尺寸综合以上考虑，确定主要尺寸如下：主轴颈：D1=（0.65-0.75）D=64.4mm L1=（0.3-0.45）D=39.56mm曲柄销：D2=（0.6-0.65）D=57.96mm L2=（0.35-0.45）D=35.88mm曲柄臂：h=（0.18-0.25）D=20.24mm b=（0.75-1.2）D=90.16mm根据主轴颈长度和曲柄销长度以及曲柄臂的厚度，确定缸心距为 L=2h+L1+L2=115.92mm5.2.2一些细节设计5.2.2.1油道布置 在确定主轴颈上油道入口和曲柄销上油道出口的位置时，既要考虑到有利于供油又要考虑到油孔对轴颈强度的影响最小。曲柄销油孔多数选择在曲拐平面运转前方=的范围内。实验表明，若油孔偏向主轴颈圆角侧，则主轴颈圆角发生裂纹的时间提前，甚至导致破坏。又曲轴的扭转疲劳破坏大多数是从应力集中的油孔边缘或过渡圆角处开始的，因此应尽量减少油孔处的应力集中，油孔边缘的过渡圆角半径取得较大并进行抛光。5.2.2.2曲轴两端的结构 消除扭转振动的减振器装在曲轴前端，因为这里的振幅最大。曲轴后端设有法兰或加粗的轴颈，飞轮与后端螺栓和定位销连接。定位销用来保证重装飞轮时保持飞轮与曲轴的装配位置。故定位销的布置是不对称的或只有一个。这种连接方式结构简单、工作可靠。 另外发动机的各种辅助装置的驱动齿轮或皮带轮及配气正时齿轮也装在曲轴前端。5.2.2.3曲轴的止推曲轴由于受热膨胀而伸长或受斜齿轮及离合器等的轴向力会产生轴向移动，为防止曲轴的轴向位移，在曲轴与机体之间设置止推轴承。止推轴承只能设置一个，以使曲轴相对于机体能自由地沿轴向作热膨胀。曲轴轴向间隙应保持=0.050.2mm，其他各主轴承端面间隙应保证曲轴受热伸长时能自由延伸。5.2.2.4过渡圆角主轴颈到曲柄臂的国度圆角半径R对于曲轴弯曲疲劳强度影响很大，增加圆角对于提高曲轴疲劳强度非常有利，但对于表面耐磨性有不利影响，在保证耐磨条件下取最大圆角。一般R取为2mm。5.2.2.5平衡分析直列四缸机的旋转惯性力，一阶往复惯性力，旋转惯性力矩，一阶往复惯性力矩，二阶往复惯性力矩都可以平衡。具体分析如下：1）旋转惯性力R Rrx=m*r*sin=0 Rry=m*r*Cos=0 Rr= (Rrx)+( Rry)1/2=02)一阶往复惯性力P=-m*r*Cos=0 3)二阶往复惯性力P =-m*r*Cos2=-4 m*r*Cos2=04）旋转惯性力矩 M=Rrx*r*Cos=2 m*r*sin2=0 M= Rry r*sin=2 m*r*sin2=0 M=M-M=05）一阶往复惯性力矩 M=P*r=-m*r*Cos=06）二阶往复惯性力矩 M=P*r=-m*r*Cos=06.曲轴强度的校核6.1 曲轴的弯曲强度校核 圆角弯曲名义应力= = (b:曲柄的宽度，h：曲柄厚度) =90.16*/6=6155.79 =*a其中a为曲柄中心到主轴颈中心的距离。a=29.9 =其中:曲柄销承受的来自连杆力中沿曲柄销半径方向的力。 ：连杆旋转质量的离心力。根据连杆大小头质量和长度的关系式，一般情况下 ：=9:2，所以=4.5*0.202=0.909kg ：曲柄销离心力。估算 =0.743kg :曲柄离心力。=0.56=0.56*0.743=0.416kg ：平衡重离心力。=0.5=0.5*0.416=0.208kg由于知道 = 其中 ，=6644.2*=6644.2*7.366=48941.2N(见附录3) =6644.2*=6644.2*(-1.04)=-6910.0N(见附录3) =（+2-2）r =(0.909+0.743+2*0.416-2*0.208)*0.046* =28182.9N =10379.2N =-17546.45N =*a=10379.2*29.9=310338.1mm =*a=-17546.45*29.9=-524638.9mm则=310338.1/6155.79=50.4MPa =-17546.45/6155.79=-8.52MPa :圆角弯曲名义应力幅值 =29.46MPa ：平均应力 =20.94MPa =)其中：有效总不均匀度系数，=1+*（-1） ：应力集中敏感系数，参考内燃机手册知=0.7 则 =1+0.7*（1.46-1）=1.322 ：曲轴材料对称循环弯曲疲劳极限。参考内燃机手册知=350MPa ：强化系数，由于设计时选取的加工工艺为滚压圆角，参考内燃机手册知=1.6 ：绝对尺寸影响系数，参考内燃机手册知=0.91 ：材料对应应力循环不对称的敏感系数，参考内燃机手册知=0.10则=350/(1.322*29.46/(1.6*0.91)+0.1*24.54)=11.86=350/11.86=29.5MPa :圆角弯曲形状系数，=/= :标准曲轴的弯曲形状系数，参考内燃机手册知=2.2 :轴颈重叠度影响系数，参考内燃机手册知=0.94 : 曲柄宽度影响系数，参考内燃机手册知=0.84 : 曲柄销空心度影响系数， 参考内燃机手册知=1 : 轴颈减重孔偏心距离影响系数，参考内燃机手册知=1 : 与圆角连接的曲柄销中减重孔至主轴颈的距离l的影响系数，参考内燃机手册知=1.17所以=2.2*0.94*0.84*1*1*1.17=2.03 =/=50.4/2.03=24.8MPa29.5MPa满足，达到设计要求。6.2曲轴的扭转强度校核 ：圆角名义扭转应力 =/其中：曲柄销抗弯断面系数 =/32=3.14*/32=19105.8() ：曲柄销圆角所受扭矩 =-r=()r 其中：曲柄销承受的来自连杆力中沿曲柄销半径切向的力由前所述，可知=/4 其中 =/4=3.96*=26311.2N=/4=-0.962*=-6391.8N:主轴颈受到的支座支反力，对于对称结构的曲轴，=/2所以，=26311.2*46/2=605157.6Nmm =-6391.8*46/2=-204537.6Nmm =/=605157.6/19105.8=31.67MPa =/=-204537.6/19105.8=-10.7MPa =(31.67+10.7)/2=21.2MPa =(31.67-10.7)/2=10.49MPa =) 其中：材料扭转疲劳极限，参考内燃机手册知=200MPa ：有效总不均匀度系数，1+*（-1）则1+0.7*（1.32-1）=1.73 ：应力集中敏感系数，参考内燃机手册知 ：因圆角半径过小引起的固有应力集中系数，参考内燃机手册知=1.32 ：强化系数，参考内燃机手册知=1.6 ：尺寸系数，参考内燃机手册知=0.89 ：材料对应力循环不对称的敏感系数，参考内燃机手册知=0.05 =200/(1.73*29.46)/(1.6*0.89)+0.05*10.49=5.51 :圆角扭转形状系数，=/其中：轴线对称之阶梯轴的扭转形状系数，参考内燃机手册知=1.30 ：曲柄宽度影响系数，参考内燃机手册知=1.57 ：曲柄厚度影响系数，参考内燃机手册知=1.26 : 轴颈重叠度影响系数，参考内燃机手册知=0.89 ：轴颈空心偏心距影响系数，参考内燃机手册知=1.0 =1.30*1.57*1.26*0.89*1.0=2.23 =/=31.67/2.23=14.2MPa =200/5.51=36.3MPa14.2MPa 满足， 该曲轴能达到设计使用的性能要求。 7小结通过这次汽车发动机设计课程设计，我对专业课的基础知识进行了深一层次的学习、整理，了解到随着内燃机的发展与强化，曲轴的工作条件愈加严酷，所以设计上正确选择曲轴的结构形式和尺寸、材料与工艺变得越来越重要，以期获得满意的技术经济效果。通过查阅各种设计手册和教材，参考网络资源，我终于完成了课程设计，收获了很多。8参考文献1 杨连生内燃机设计北京：中国农业机械出版社，19812 徐兀汽车发动机现代设计人民交通出版社3 唐增宝，何永然，刘安俊机械设计课程设计武汉：华中科技大学出版社，19994 周龙保内燃机学北京：机械工业出版社，20055 陈家瑞.汽车构造. 人民交通出版社，20026 吴兆汉内燃机设计北京：北京理工大学出版社，19907 沈维道工程热力学北京：高等教育出版社，2002附 录附表1 一些需要的参数活塞行程S（mm）活塞直径D(mm)转速n（r/min）往复惯性质量mj(kg)压缩比余隙容积Vc（L)929252000.45590.076曲柄半径r(mm)连杆长度l(mm)半径/长度比曲柄角速度（rad/s)进气压力Pa（MPa)压力升高比461840.25544.30.096多变指数n1多变指数n2气缸总容积Va(L)排量（L)旋转运动质量mr(kg)活塞面积Fh（cm2)1.351.250.690.612.27666.5 附表2 运动学计算图表对应弧度X（m）V(m/s)a（m/s)气缸体积（L）Pg（MPa）00.0000 0.0000 0.0000 17033.0068 0.0000 0.0760 0.0900 50.0873 0.0002 2.7255 16929.4002 -0.2421 0.0775 0.0900 100.1745 0.0009 5.4179 16620.5478 -0.1364 0.0818 0.0900 150.2618 0.0020 8.0446 16112.3002 0.1633 0.0890 0.0900 200.3491 0.0034 10.5745 15414.2406 0.2303 0.0989 0.0900 250.4363 0.0053 12.9781 14539.4385 -0.0331 0.1115 0.0900 300.5236 0.0076 15.2284 13504.1139 -0.2496 0.1265 0.0900 350.6109 0.0102 17.3010 12327.2241 -0.1073 0.1438 0.0900 400.6981 0.0131 19.1750 11029.9823 0.1874 0.1633 0.0900 450.7854 0.0163 20.8328 9635.3237 0.2144 0.1846 0.0900 500.8727 0.0198 22.2609 8167.3344 -0.0656 0.2076 0.0900 550.9599 0.0235 23.4493 6650.6588 -0.2526 0.2320 0.0900 601.0472 0.0273 24.3923 5109.9020 -0.0763 0.2575 0.0900 651.1345 0.0313 25.0879 3569.0466 0.2082 0.2838 0.0900 701.2217 0.0353 25.5380 2050.8971 0.1947 0.3108 0.0900 751.3090 0.0395 25.7479 576.5699 -0.0971 0.3382 0.0900 801.3963 0.0436 25.7263 -834.9577 -0.2511 0.3656 0.0900 851.4835 0.0477 25.4844 -2167.2279 -0.0440 0.3929 0.0900 901.5708 0.0518 25.0363 -3406.6014 0.2254 0.4198 0.0900 951.6581 0.0557 24.3976 -4542.4669 0.1717 0.4462 0.0900 1001.7453 0.0596 23.5855 -5567.3586 -0.1269 0.4718 0.0900 1051.8326 0.0633 22.6184 -6476.9765 -0.2451 0.4964 0.0900 1101.9199 0.0668 21.5148 -7270.1132 -0.0111 0.5199 0.0900 1152.0071 0.0702 20.2932 -7948.4889 0.2386 0.5422 0.0900 1202.0944 0.0733 18.9718 -8516.5034 0.1457 0.5631 0.0900 1252.1817 0.0762 17.5677 -8980.9113 -0.1546 0.5826 0.0900 1302.2689 0.0789 16.0969 -9350.4347 -0.2347 0.6005 0.0900 1352.3562 0.0814 14.5738 -9635.3237 0.0221 0.6169 0.0900 1402.4435 0.0836 13.0110 -9846.8820 0.2476 0.6316 0.0900 1452.5307 0.0856 11.4194 -9996.9716 0.1172 0.6446 0.0900 1502.6180 0.0873 9.8079 -10097.5126 -0.1797 0.6559 0.0900 1552.7053 0.0887 8.1834 -10159.9962 -0.2201 0.6655 0.0900 1602.7925 0.0899 6.5513 -10195.0246 0.0549 0.6733 0.0900 1652.8798 0.0908 4.9151 -10211.8936 0.2521 0.6794 0.0900 1702.9671 0.0915 3.2771 -10218.2315 0.0868 0.6838 0.0900 1753.0543 0.0919 1.6386 -10219.7054 -0.2016 0.6864 0.0900 1803.1416 0.0920 0.0000 -10219.8041 -0.2017 0.6873 0.0900 1853.2289 0.0919 -1.6386 -10219.7054 0.0868 0.6864 0.0902 1903.3161 0.0915 -3.2771 -10218.2315 0.2521 0.6838 0.0906 1953.4034 0.0908 -4.9151 -10211.8936 0.0549 0.6794 0.0914 2003.4907 0.0899 -6.5513 -10195.0246 -0.2201 0.6733 0.0925 2053.5779 0.0887 -8.1834 -10159.9962 -0.1797 0.6655 0.0940 2103.6652 0.0873 -9.8079 -10097.5126 0.1172 0.6559 0.0959 2153.7525 0.0856 -11.4194 -9996.9716 0.2476 0.6446 0.0981 2203.8397 0.0836 -13.0110 -9846.8820 0.0221 0.6316 0.1009 2253.9270 0.0814 -14.5738 -9635.3237 -0.2347 0.6169 0.1041 2304.0143 0.0789 -16.0969 -9350.4347 -0.1546 0.6005 0.1080 2354.1015 0.0762 -17.5677 -8980.9113 0.1457 0.5826 0.1125 2404.1888 0.0733 -18.9718 -8516.5034 0.2386 0.5631 0.1178 2454.2761 0.0702 -20.2932 -7948.4889 -0.0111 0.5422 0.1240 2504.3633 0.0668 -21.5148 -7270.1132 -0.2451 0.5199 0.1312 2554.4506 0.0633 -22.6184 -6476.9765 -0.1269 0.4964 0.1396 2604.5379 0.0596 -23.5855 -5567.3586 0.1717 0.4718 0.1496 2654.6251 0.0557 -24.3976 -4542.4669 0.2254 0.4462 0.1613 2704.7124 0.0518 -25.0363 -3406.6014 -0.0440 0.4198 0.1751 2754.7997 0.0477 -25.4844 -2167.2279 -0.2511 0.3929 0.1915 2804.8869 0.0436 -25.7263 -834.9577 -0.0971 0.3656 0.2110 2854.9742 0.0395 -25.7479 576.5699 0.1947 0.3382 0.2344 2905.0615 0.0353 -25.5380 2050.8971 0.2082 0.3108 0.2627 2955.1487 0.0313 -25.0879 3569.0466 -0.0763 0.2838 0.2970 3005.2360 0.0273 -24.3923 5109.9020 -0.2526 0.2575 0.3388 3055.3233 0.0235 -23.4493 6650.6588 -0.0656 0.2320 0.3900 3105.4105 0.0198 -22.2609 8167.3344 0.2144 0.2076 0.4530 3155.4978 0.0163 -20.8328 9635.3237 0.1874 0.1846 0.5307 3205.5851 0.0131 -19.1750 11029.9823 -0.1072 0.1633 0.6264 3255.6723 0.0102 -17.3010 12327.2241 -0.2496 0.1438 0.7434