叶片前缘参数的计算方法

叶片前缘参数的计算方法
叶片前缘参数的计算方法

题目:叶片前缘参数的计算方法

【摘要】航空发动机叶片及其截面线可以分为叶盆、叶背、前缘和后缘四个部分,其中前缘的形状对整个叶片的气动性能具有非常重要的作用。传统的方法是将前缘设计成一段圆弧但是近年来越来越多的实验和数值研究显示,采用非圆弧形前缘,例如椭圆弧形前缘可以明显改善叶片的气动性能。

针对问题1,我们建立了2个模型。

模型I :分析航天发动机叶片前缘参数,针对问题1中的数据分离,我们可建立模型,先设出椭圆的一般方程,再求解该方程的未知参数,从而可以分离出前缘数据。

模型Ⅱ:为了求出椭圆弧的中心坐标()00,x y 、长短半轴a 、b ,倾角,可以建立一个模型,设出椭圆的标准方程,从而求得叶片重要的前缘参数。

叶片前缘参数的计算方法

1.问题重述

航空发动机叶片及其截面线的形状如图1所示,它可以分为叶盆、叶背、前缘和后缘四个部分,其中前缘的形状对整个叶片的气动性能具有非常重要的作用。传统的方法是将前缘设计成一段圆弧,则这段圆弧的圆心坐标()00,x y 和半径(也称前缘半径)r 就是重要的前缘参数。但是近年来越来越多的实验和数值研究显示,采用非圆弧形前缘,例如椭圆弧形前缘可以明显改善叶片的气动性能。如果是椭圆弧形的前缘,则椭圆弧的中心坐标()00,x y 、长短半轴a 、b ,倾角β(椭圆长轴方向与y 轴的夹角)就是重要的前缘参数。

图1 航空发动机叶片及其截面线

图2. 测量得到的截面线前缘附近的离散数据

在某型航空发动机的仿制过程中,如图2所示测量了截面线前缘附近的一组离散数据。它们的坐标值记录在附件qianyuan_data.txt 中, 文件中的每一行代表一个数据点的坐标(第一个值为横坐标x ,第二个值为纵坐标y )。由于在测量时前缘数据与相邻的叶盆叶背数据无法直接分离,故其中混有相邻的叶盆叶背数据,请建立数学模型,解决以下问题:

1)

根据所给数据,通过计算分离出前缘数据,并判断出前缘的形状是一段圆弧还是一段椭圆弧,计算出相应的前缘参数。将整个由离散数据出发,经过分离前缘数据、判断前缘类型、计算前缘参数的过程用一个算法描述。

2)

评价你所给出模型的精确度和稳定性(指如果给数据添加均值较小的随机扰动,是否还能得到较高的计算结果),给出计算实例。

3) 能否对1)中计算前缘参数的方法进行改进提高计算精度,给出具体的方法并比较改进前后的计算结果。

2.模型的假设与符号说明 2.1.模型的假设:

实际测量所得到的前缘附近数据是离散的,在这里设各点的连线是一段 光滑曲线。 2.2.符号说明:

圆心坐标()00,x y , 半径r

A.B.C.D.E.F.分别为椭圆方程的系数;

(i y x , ) ( i , j = 1,2,3………) 为测量得到的前缘数据的坐标; (00,y x ) 为叶片前缘椭圆弧的中心坐标;

a ,

b 为叶片前缘椭圆弧的长短半轴;

β 叶片前缘椭圆长轴方向与y 轴方向的夹角。

3.模型的分析

从混有相邻的叶盆和叶背数据中,分离出前缘数据,必须找到只有前缘数据才满足的条件,如方程,将数据逐一代人,条件成立的数据,即为前缘数据,实现前缘数据的筛选。

因为叶片前缘是一段圆弧或椭圆弧,叶片前缘的数据都应满足圆或椭圆的方程,而叶盆叶背的数据则不满足此方程。就能做出判断数据是否属于叶片前缘。

求出叶片前缘数据所满足的方程是问题的关键。

假设叶片前缘是一段椭圆弧,求出椭圆一般方程需知道圆弧上的5个点的坐标。

要求出椭圆的参数(0,0y x ),a,b ,β值,须知道椭圆的标准方程,通过建立模型求解。

4.对问题一中模型的建立与计算 4.1模型一的建立:

由于圆是椭圆的特殊情况,所以可以直接设叶片前缘类型为椭圆弧,根据椭圆的一般方程建立模型:

Ax 2+By 2+Cxy+Dx+Ey+F=0

① 如果A=B 则前缘类型为圆; ② 如果A ≠B 则前缘类型为椭圆。 4.2.模型一的计算:

第一步:分析截面线前缘附近数据点,可得到相对准确的6个点在椭圆弧上,分

别是:

(-18.9572,-0.241143),(-18.9982,-0.0403025),(-19.1117,0.699986), (-19.0947,0.886326),(-18.7235,2.00039),(-18.6762,2.09696)

代入这6个点坐标,可得一个五元一次线性方程组:

?????

?

??

???=+++++=+++++=+++++=+++++=+++++=+++++000000666626265555252544442

42433332

32322222

22211112

121F Ey Dx y Cx By Ax F Ey Dx y Cx By Ax F Ey Dx y Cx By Ax F Ey Dx y Cx By Ax F Ey Dx y Cx By Ax F Ey Dx y Cx By Ax 由于线性方程组中常数项全为0,故此方程组为齐次线性方程组,而系数行 的列式为 :

1

111116

6

6

62

6

26

55552

52544442

42433332

32322222

222

11112

12

1y x y x y x y x y x y x y x y x y x y x y x y x y x y x y x y x y x y x F =

由此行列式可由MATLAB 软件编程可求得。 由于所求行列式F ≠0,故次方程组只有零解。

第二步:为了能求出精确且合理的系数,可设方程组的解F=-K (其中K ∈R ).则可得非齐次线性方程:Ax 2+By 2+Cxy+Dx+Ey=K 。 则在上述分析所得到的五组数据中任意找5点,即

(-19.1117,0.699986),(-19.0947,0.886326),(-18.7235,2.00039), (-18.6762,2.09696),(-18.5912,2.19285)

代入方程:Ax 2+By 2+Cxy+Dx+E=K ,得到非齐次线性方程组:

??

??

????

?=++++=++++=++++=++++=++++K Ey Dx y Cx By Ax K Ey Dx y Cx By Ax K Ey Dx y Cx By Ax K Ey Dx y Cx By Ax K Ey Dx y Cx By Ax 55552

525444424243333232322222

22211112

121

由此推出行列式为:

???????

?

??=55

5

52

5

25

44442

42

433332

32

322222

222

11112

12

1y x y x y x y x y x y x y x y x y x y x y x y x y x y x y x G ;

???????

? ??=E D C B A X

常数项行列式:???????

? ??=K

K K

K

K H 则GX=H.

该方程的解可由附录1中的程序在MALAB 软件上运行求得:

则椭圆的系数为A =-0.0030*k,B=-0.00070442*k,C=-0.00000044023,D=-0.1102,E=0.00061973;

把所求得系数代入方程Ax 2+By 2+Cxy+Dx+Ey=K 则可得:

(0.003x 2+-0.00070442y 2+-0.00000044023xy+-0.1102x+0.00061973y )k=k. 第三步:由此得k 为非0,则:

0.003x 2+-0.00070442y 2+-0.00000044023xy+-0.1102x+0.00061973y=1. ① 要分离出前缘数据,可把数据代入此方程,看是否满足条件: 0.003x 2

+-0.00070442y 2

+-0.00000044xy+-0.1102x+0.00061973y-1

≤0.08

② 由附录2中的C程序可判断分离出的前缘数据如下:

(-11.5725,-8.90186),(-12.6076,-8.26876),(-12.8985,-8.08197), (-13.257,-7.79105),(-15.0655,-6.25705), (-15.7753,-5.50894), (-15.8267,-5.45908),(-15.8516,-5.4325) ,(-16.0351,-5.21541), (-16.1321,-5.10414),(-16.1491,-5.082), (-16.279,-4.90552), (-16.7902,-4.22412),(-17.5836,-3.025), (-17.6692,-2.89498), (-17.7415,-2.7646),(-18.2115,-1.91654),(-18.8601,-0.48739), (-18.9239,-0.346534),(-18.9572,-0.241143),(-18.9982,-0.0403025), (-19.1117,0.699986) ,(-19.0947,0.886326) ,(-18.7235,2.00039), (-18.6762,2.09696),(-18.5912,2.19285),(-18.1736,2.65479), (-17.9496,2.79082), (-17.6239,2.98724), (-17.0294,3.14561), (-16.8289,3.1809), (-16.5749,3.15474),(-16.005 ,3.11765), (-15.6571,3.00622), (-15.2361,2.8083),(-14.7976,2.52742), (-14.1023,2.09416), (-13.8554,1.96404), 4.3模型二的建立。

求椭圆的参数是,必须得知道椭圆的标准方程,因此可以建立模型二,进行求解。

1

]

sin )(cos )([]

cos )(sin )[(2

2

002

2

00=-+--+

-+-b

y y x x a

y y x x ββββ根据模型一求得的椭圆一般方程:

0.003x 2+-0.00070442y 2+-0.00000044023xy+-0.1102x+0.00061973y=1 可求得具体的参数(0,0y x ),a,b ,β值。

5.模型的检验

5.1 模型的精确度分析与讨论

对建立的模型求解是,我们只是通过分析数据,找了相对比较准确的5组离散数据,带入方程组求得系数A,B,C,D,E,F 的解,但由于所测得数为离散的,在选点时必定会有误差,而且在MATLANB 软件上所求得数据虽为最精确值,但由于其数字太长,我们只能取其近似值,这必然会带来一定的误差,影响实验的数据。

在分离前缘数据时,我们是通过条件:

0.003x2+-0.00070442y2+-0.00000044xy+-0.1102x+0.00061973y-1 0.08 判断满足条件的前缘数据,而此时的精确度只到0.08,但如果增大其精确度,满足条件的数据自然会减少,缩小精确度,得到的数据太不精确,所以只有取适当的精确度求解。

5.2模型稳定性的分析与讨论

为了判断模型所得结果的稳定性,我们可给数据添加均值较小的随机扰动,方法同模型一中的计算完全一样,通过计算我们可得出较高的计算结果。

6.模型的进一步分析与讨论

为了提高前缘参数的精确度,我们可通过求均值的方法得更准确的值,或通过用最小二乘解求模型的方程。

同理分析得到6组数据,在这6组数据中任选5组代入方程组,有6种选法,这6种方法可求的6组不同的A,B,C,D,E.通过计算求得其平均值。

附录1:

>> syms k;A=sym([18.9572*18.9572 0.241143*0.241143 18.9572*(-0.241143) -18.9572 -0.241143; 18.9982*18.9982 0.0403025*0.0403025 18.9982*(-0.040325) -18.9982 -0.0403025;19.1117*19.1117 0.699986*0.699986 19.1117*0.699986 -19.1117 0.699986;19.0947*19.0947 0.886326*0.886326 19.0947*0.886326 -19.0947 0.886326;18.7352*18.7352 2.00039*2.00039 18.7352*2.00039 -18.7352 2.00039]); >> b=sym([k;k;k;k;k]);S1=A\b

S1 =

-1245788493499456970433983768238997765579718487697344827417316355058970 4994816/461292085021707399954519484748287329300964810893471597422687443 9622622529411555*k

4893338813900759464044798158947450594910180125564877514414357503617550 778368/2306460425108536999772597423741436646504824054467357987113437219 8113112647057775*k

-9058868471264516166183907064206905717477367760707889551402809253220304 551936/2306460425108536999772597423741436646504824054467357987113437219 8113112647057775*k

-9589661754370946199737396123319724709290558062331516938023517493919072

9113600/922584170043414799909038969496574658601929621786943194845374887 924524505882311*k 1662375975104813257260738822076791052957483485859607725150452661810145 72498944/23064604251085369997725974237414366465048240544673579871134372 198113112647057775*k

>>

A=-377875778385690354981326726967070890998012205084567414144598189 4439059062784/12469019793841652044787904264245784157156340

48871183618898565875274408901382079

A =-0.0030

>>

B=-878336522276960353681668402845807090680165106432104152002379575 444690173952/124690197938416520447879042642457841571563404

8871183618898565875274408901382079

B=-7.0442e-004

>>

C=-548927144593865254706058540533394401441932229521806274722588633 897893888/124690197938416520447879042642457841571563404887

1183618898565875274408901382079

C =-4.4023e-007

>>

D=-137455946265195794012784835013500721489179151363177007208220611 217949982720000/124690197938416520447879042642457841571563

4048871183618898565875274408901382079

D = -0.1102

>>

E=7727365118520173012333743710160626516055116434325818915401359289 22112000000/1246901979384165204478790426424578415715634048

871183618898565875274408901382079

E = 6.1973e-004

附录2:

#include

#include

main()

{

double a[58][2]= {{-6.53163,-10.4501},{-7.32463,-10.3246},{-8.09575,-10.2018},

{-8.19995,-10.1822},{-9.63034,-9.75588},{-10.1825,-9.55725},

{-10.61,-9.35925},{-11.5725,-8.90186},{-12.6076,-8.26876},

{-12.8985,-8.08197},{-13.257,-7.79105},{-15.0655,-6.25705},

{-15.7753,-5.50894},{-15.8267,-5.45908},{-15.8516,-5.4325},

{-16.0351,-5.21541},{-16.1321,-5.10414},{-16.1491,-5.082},

{-16.279,-4.90552},{-16.7902,-4.22412},{-17.5836,-3.025},

{-17.6692,-2.89498},{-17.7415,-2.7646},{-18.2115,-1.91654}

{-18.8601,-0.48739},{-18.9239,-0.346534},{-18.9572,-0.241143}

{-18.9982,-0.0403025},{-19.1117,0.699986},{-19.0947,0.886326},

{-18.7235,2.00039}, {-18.6762,2.09696},{-18.5912,2.19285},

{-18.1736,2.65479},{-17.9496,2.79082},{-17.6239,2.98724},

{-17.0294,3.14561},{-16.8289,3.1809},{-16.5749,3.15474},

{-16.005,3.11765},{-15.6571,3.00622},{-15.2361,2.8083},

{-14.7976,2.52742},{-14.1023,2.09416},{-13.8554,1.96404},

{-13.3816,1.74286},{-11.7274,0.883965}, {-11.3153,0.737037},

{-9.67702,0.282667}, {-9.65347,0.275782}, {-9.63495,0.272118},

{-9.09549,0.165523},{-8.05064,0.0817958},{-7.58706,0.0411919},

{-7.26465,0.0377957},{-6.22285,0.0596912},{-5.39057,0.159466},

{-4.98982,0.230775}};

int i,c=0;

for(i=0;i<58;i++)

if(fabs(a[i][0]*a[i][0]*(-0.003)+a[i][1]*a[i][1]*0.000000704442+a[i][0]*a[i][1]*(-0 .00044023)+a[i][0]* (-0.1102)+a[i][1]*(-0.00061973)-1)<=0.01)

{

c++;

printf("满足的数据:a1=%lf,a2=%lf\n",a[i][0],a[i][1]);

}

printf("总共有个%d数据",c);

}

附件qianyuan_data.txt

-6.53163 -10.4501

-7.32463 -10.3246

-8.09575 -10.2018

-8.19995 -10.1822

-9.63034 -9.75588

-10.1825 -9.55725 -10.61 -9.35925 -11.5725 -8.90186 -12.6076 -8.26876 -12.8985 -8.08197 -13.257 -7.79105 -15.0655 -6.25705 -15.7753 -5.50894 -15.8267 -5.45908 -15.8516 -5.4325

-16.0351 -5.21541 -16.1321 -5.10414 -16.1491 -5.082 -16.279 -4.90552 -16.7902 -4.22412 -17.5836 -3.025 -17.6692 -2.89498 -17.7415 -2.7646

-18.2115 -1.91654 -18.8601 -0.48739 -18.9239 -0.346534 -18.9572 -0.241143 -18.9982 -0.0403025 -19.1117 0.699986 -19.0947 0.886326 -18.7235 2.00039

-18.6762 2.09696

-18.5912 2.19285

-18.1736 2.65479

-17.9496 2.79082 -17.6239 2.98724 -17.0294 3.14561 -16.8289 3.1809 -16.5749 3.15474 -16.005 3.11765

-15.6571 3.00622 -15.2361 2.8083 -14.7976 2.52742 -14.1023 2.09416 -13.8554 1.96404 -13.3816 1.74286 -11.7274 0.883965 -11.3153 0.737037 -9.67702 0.282667 -9.65347 0.275782 -9.63495 0.272118 -9.09549 0.165523 -8.05064 0.0817958 -7.58706 0.0411919 -7.26465 0.0377957 -6.22285 0.0596912 -4.98982 0.230775

各地转换参数及转换参数的计算方法

坐标转换 一、中央经线(LONGITUDEORIGIN) 在坐标转换中,首先需要设置测区的中央经线,以下是新疆各地州的中央经线,仅供参考。 乌鲁木齐E87度 吐鲁番E87度 鄯善E93度 哈密E93度 阿勒泰E87度 塔城E81度 克拉玛依E87度 奎屯E87度 博乐E81度 伊犁E81度 阿克苏E81度 库尔勒E87度 喀什E75度 和田E81度 二、投影比例(SCALE) 系统一般默认值时+0.9996。将改值改为1 三、东西偏差(ALSEE) 系统一般默认值:+1000000.0m。将该值改为:+500000.0m 四、南北偏差(FALSEN) 系统一般默认值:+100000.0m。将该值改为:+0.0m 五、dx\dy\dz\da\df DX、DY、DZ是坐标在三个方向的平移量,原则上在不同的地区,值是不一样的。 六、下面用软件COORD 进行转换!!! 以下面这个实例来求解转换参数:某林内有一个北京-54坐标系下的已知点,中央经线E117°,属于3度带,其坐标为X=4426818.5,Y=456613.7,h=63.9,其对应的WGS84坐标系统下的坐标为B=39°58′27.120″N,L=116°29′32.874″E,H=58。.由这两套坐标进行系统坐标转换三参数Dx、Dy、Dz求解。 打开COORD转换软件,如图:

1、请按步骤操作,点击坐标转换,选择投影设置。测量地区属于高斯投影3度带的选择高斯投影3度带,测量地区属于高斯投影6度带的选择高斯投影6度带,中央子午线根据所在地区中央经线填入。由实例填入中央经线117度,高斯投影3度带。 图2 2、点击坐标转换,选择计算三参数。此时,需要到当地测绘部门去咨询当地的一个已知点的大地坐标和平面坐标。将大地坐标的三个参数和平面坐标的三个参数填入。左边椭球基准,选择WGS-84坐标系。右边根据用户要求可选择北京-54坐标系或者国家-80坐标系,点击确定。 如图3,由实例,我们填入大地坐标和平面坐标 图3

南方gps坐标转换参数设置

注:新版本已将"控制点坐标库"改为"求转换参数",实现的功能不变! 一、控制点坐标库的应用 GPS 接收机输出的数据是WGS-84 经纬度坐标,需要转化到施工测量坐标,这就需要软件进行坐标转换参数的计算和设置,控制点坐标库就是完成这一工作的主要工具。 控制点坐标库是计算四参数和高程拟合参数的工具,可以方便直观的编辑、查看、调用参与计算四参数和高程拟合参数的校正控制点。 利用控制点坐标库可以计算GPS 原始记录坐标到当地施工坐标的参数。在计算之前,需新建工程,输入当地的施工坐标系及中央子午线、投影高等。假设我们利用A、B 这两个已知点来求取参数,那么首先要有A、B 两点的GPS 原始记录坐标和测量施工坐标。 A、B 两点的GPS原始记录坐标的获取有两种方式: 一种是布设静态控制网,采用静态控制网布设时后处理软件的GPS 原始记录坐标; 另一种是GPS 移动站在没有任何校正参数起作用的Fixed(固定解)状态下记录的GPS 原始坐标。 1.1、校正参数 操作:工具→校正向导或设置→求转换参数(控制点坐标库) 所需已知点数:1个 校正参数是工程之星软件很特别的一个设计,它是结合国内的具体测量工作而设计的。校正参数实际上就是只用同一个公共控制点来计算两套坐标系的差异。根据坐标转换的理论,一个公共控制点计算两个坐标系误差是比较大的,除非两套坐标系之间不存在旋转或者控制的距离特别小。因此,校正参数的使用通常都是在已经使用了四参数或者七参数的基础上才使用的。

在工程之星新版本中,在校正向导中已经取消了两点校正功能,如果两个以上的已知点请使用控制点坐标库来求取参数。习惯使用校正向导的人请尽快学习新版本。 1.2 四参数 操作:设置→求转换参数(控制点坐标库) 四参数是同一个椭球内不同坐标系之间进行转换的参数。在工程之星软件中的四参数指的是在投影设置下选定的椭球内GPS 坐标系和施工测量坐标系之间的转换参数。工程之星提供的四参数的计算方式有两种,一种是利用“工具/参数计算/计算四参数”来计算,另一种是用“控制点坐标库”计算。。需要特别注意的是参予计算的控制点原则上至少要用两个或两个以上的点,控制点等级的高低和分布直接决定了四参数的控制范围。经验上四参数理想的控制范围一般都在5-7 公里以内。 四参数的四个基本项分别是:X 平移、Y 平移、旋转角和比例。 从参数来看,这里没有高程改正,所以建议采用“控制点坐标库”来求取参数,而根据已知点个数的不同所求取的参数也会不同,具体有以下几种。 1.2.1 四参数+校正参数 所需已知点个数:2个

电厂烟气环境监测常用计算公式

1.1.1 烟气流量的计算 s s V F Q ??=3600 (式 4-1) 式中:s Q -湿烟气排放量,m 3/h ; F -测定断面面积,m 2; s V -测定断面的平均烟气流速,m/s 。 1.1.2 标态下干烟气排放量的计算 )1() 273(101325273 sw s s a s m X t )P (B Q Q -?+??+?= (式4-2) 式中:m Q -标准状态下干烟气的排放量,Nm 3/h ; sw X -烟气中水分含量体积百份数,%; a B -大气压力,Pa ; s p -测点处烟气静压,Pa ; s t -烟气温度,℃。 1.1.3 采样体积的计算 s t P B V V s a m snd ++? =2730027.0 (式4-3) 式中:snd V -标准状态下的干烟气采样体积,L ; m V -实际工况下的干烟气采样体积,L ; s P -烟气静压,Pa ; s t -烟气温度,℃。 1.1.4 烟气含尘浓度计算 3 10?= snd V g C (式4-4) 式中:C -标准状态下干燥烟气的含尘浓度,mg/Nm 3; g -所采得的粉尘量,mg ;

21g g g -=; 1g -采样前滤筒质量,mg ; 2g -采样后滤筒质量,mg 。 1.1.5 烟尘排放量的计算 6 10m m Q C q ?= (式 4-5) 式中:m q -烟尘排放量 kg/h 。 1.1.6 漏风率的计算 % 100222?--= ?out in out O K O O α (式4-6) 式中:α?-除尘器漏风率,%; out O 2-除尘器出口断面烟气平均氧量,%; in O 2-除尘器入口断面烟气平均氧量,%; K -大气中的含氧量,%。 1.1.7 除尘效率的计算 % 100) 1(??+-= in out in C C C αη (式4-7) 式中:η-除尘效率,%; in C -进口烟尘浓度(标态干烟气),mg/m 3; out C -出口烟尘浓度(标态干烟气),mg/m 3。 1.1.8 除尘器本体压力降计算 H out in p p p p +-=? (式 4-8) 式中:p ?-除尘器压力降,Pa ; in p -除尘器入口全压平均值,Pa ; out p -除尘器出口全压平均值,Pa ; H p -高温气体浮力的校正值,Pa 。

计算转换参数方法

计算转换参数方法 一、知道基准站的北京54坐标x,y 如果是这种情况,那么直接可以用软件计算得到基准站的54经纬度,直接用这个坐标设置基准站就可以了,在软件里面四个转换参数都为默认值:dx,dy,scale,rotation 举例:54坐标x = 3391531.060 y = 408652.459 45.917 在菜单坐标转换里面选择投影设置,如下图1: 在这里选择投影方式和设置中央子午线,一般情况都是3度带和6度带。 然后如下图进行经纬度转换,如图2: 图中两边的椭球基准都选择北京-54坐标系,左边在“选择源坐标类型”里面选择“平面坐标”,右边的“选择目标坐标类型”中选择“大地坐标”。 这里计算的经纬度就是使用基准站的54坐标转换得到的。我们就可以把这个坐标设置到基准站的GPS中去。 在海测软件里面我们也只要设置中央子午线就可以了。4个地方转换参数都是默认值,不用

设置。 二、知道两个wgs84坐标和两个地方坐标,其中一个是基准站的坐标 这时候,我们可以很方便的设置基准站,但是在流动站的船上,我们必须设置相关的坐标转换参数。现在我们在软件里面设计了输入地方坐标转换参数的对话框。主要的工作就是要求出这四个转换参数。需要有几个步骤,但是对于一个工程来说,只要在工作前花半个小时就可以了。具体步骤和图示如下: 1、84经纬度转换成54投影坐标 如图1设置投影带和中央子午线 2、把两个点的经纬度转换成54投影坐标,如下图 如上图的设置,在左边的椭球基准选择WGS-84坐标系,右边也是相同的坐标系。举例已知的两组数据如下: 点一 B = 030:38:26.645 L = 122:02:49.556 地方坐标 x’ = -65839.283 y’ = 55680.371 点二 B = 030:37:59.928 L = 122:03:07.031 地方坐标 x’ = -66659.526 y’ = 56150.074 分别得到高斯平面投影坐标 点一 x = 3391469.448 y = 408651.927 点二 x = 3390704.294

手持GPS三参数计算及各地坐标转换经验参数

如何设置手持GPS相关参数及全国各地坐标转换参数一、如何设置手持GPS相关参数 (一)手持GPS的主要功能 手持GPS,指全球移动定位系统,是以移动互联网为支撑、以GPS智能手机为终端的GIS系统,是继桌面Gis、WebGis之后又一新的技术热点。目前功能最强的手持GPS,其集成GPRS通讯、蓝牙技术、数码相机、麦克风、海量数据存储、USB/RS232端口于一身,能全面满足您的使用需求。 主要功能:移动GIS数据采集、野外制图、航点存储坐标、计算长度、面积角度(测量经纬度,海拔高度)等各种野外数据测量;有些具有双坐标系一键转换功能;有些内置全国交通详图,配各地区地理详图,详细至乡镇村落,可升级细化。 (二)手持GPS的技术参数 因为GPS卫星星历是以WGS84大地坐标系为根据建立的,手持GPS单点定位的坐标属于WGS84大地坐标系。WGS84坐标系所采用的椭球基本常数为:地球长半轴a=6378137m;扁率F=1/298.257223563。 常用的北京54、西安80及国家2000公里网坐标系,属于平面高斯投影坐标系统。北京54坐标系,采用的参考椭球是克拉索夫斯基椭球,该椭球的参数为:地球长半轴a=6378245m;扁率F=1/298.2。西安80坐标系,其椭球的参数为:地球长半轴a=6378140m;扁率F=1/298.257。国家2000坐标系,其椭球的参数为:地球长半轴a=6378137m;扁率F=1/298.298.257222101。 (三)手持GPS的参数设置 要想测量点位的北京54、西安80及国家2000公里网高精度坐标数据,必须学习坐标转换的基础知识,并分别科学设置手持GPS的各项参数。 首先,在手持式GPS接收机应用的区域内(该区域不宜过大),从当地测绘部门收集1至两个已知点的北京54、西安80或国家2000坐标系统的坐标值;然后在

(完整版)烟气量计算公式

燃料空气需要量及燃烧产物量的计算 所有理论计算均按燃料中可燃物质化学当量反应式,在标准状态下进行,1kmol 反 应物质或生成物质的体积按22.4m 3计,空气中氧和氮的容积比为21:79,空气密度为 1.293kg/m 3。 理论计算中空气量按干空气计算。燃料按单位燃料量计算,即固体、液体燃料以1kg 计算,气体燃料以标准状态下的1m 3计算。 单位燃料燃烧需要理论干空气量表示为L 0 g ,实际燃烧过程中供应干空气量表示为 Ln g ; 单位燃料燃烧理论烟气量表示为V 0,实际燃烧过程中产生烟气量表示为Vn; 单位燃料燃烧理论干烟气量表示为V 0g ,实际燃烧过程中产生干烟气量表示为Vn g ; 一、通过已知燃料成分计算 1. 单位质量固体燃料和液体燃料的理论空气需要量(m 3/kg ) L 0=(8.89C +26.67H +3.33S -3.33O )×10﹣2式中的C 、H 、O 、S ——燃料中收到基 碳、氢、氧、硫的质量分数%。 2. 标态下单位体积气体燃料的理论空气需要量(m 3/m 3) L 0=4.76?? ????-+??? ??+++∑2222342121 O S H?CmHn n m H CO ×10﹣2式中CO 、H 2、H 2O 、H 2S 、CmHn 、O 2——燃料中气体相应成分体积分数(%). 3. 空气过剩系数及单位燃料实际空气供应量 空气消耗系数а=0 L 量单位燃料理论空气需要量单位燃料实际空气需要?L 在理想情况下,а=1即能达到完全燃烧,实际情况下,а必须大于1才能完全燃烧。а<1显然属不完全燃烧。 а值确定后,则单位实际空气需要量L а可由下式求得: L 0g =аgL 0 以上计算未考虑空气中所含水分 4. 燃烧产物量 a.单位质量固体和液体燃料理论燃烧产物量(m 3/kg) 当а=1时, V 0=0.7L 0+0.01(1.867C+11.2H+0.7S+1.244M+0.8N)式中 M ——燃料中水分(%)。 b.单位燃料实际燃烧产物量(m 3/kg ) 当a >1时,按下式计算: 干空气时,V a =V 0+(a-1)L 0 气体燃料 (2)单位燃料生成湿气量 ?V =1+α0L -[0.5H 2+0.5C O -(4 n -1) C m H n ] (标米3/公斤) (2-14) (3)单位干燃料生成气量 g V ?=1+α0L -[1.5H 2+0.5C O -( 4n -1) C m H n +2 n C m H n ) (标米3/公斤) (2-15)

坐标转换三参数计算器使用说明

坐标转换三参数计算器使用说明 4.0升级及使用说明: 1、增加了批量处理数据功能。 2、经纬度数据与直角坐标数据可混合输入(经纬度格式:DDD.MMSS,109度04分08.94343秒表示为109.040894343,直角坐标格式单位为米,如X为1234567.89,Y为123456.78,Y坐标无带号)。 3、批量处理数据文件为文本文件,格式为严格每行4个数据,以逗号或空格分开。 点号1,X坐标(或为纬度),Y坐标(或为经度),高程 4、输出文件为文本文件,格式为: 点号1,转换前的X坐标(或为纬度),Y坐标(或为经度),高程 > 转换后的X坐标(或为纬度),Y坐标(或为经度),高程 5、未注册软件无批量处理功能,部分参数隐形显示,但内部坐标转换仍可正常进行。 工作界面:

=========================================== 3.0使用说明 本软件分成上下二部分,上半部为在两个不同椭球体间求坐标转换的三参数,下半部为在两个不同椭球体间的坐标转换。在两个不同椭球体间进行坐标转换首先必需知道坐标转换参数,通常有三参数和七参数转换二种方式,本程序提供三参数转换方式。 例1:我要求手持GPS的北京54(或西安80)坐标转换参数。 向有关部门收集所在工作区内已知点(只要一个控制点)的WGS84坐标系中的经度、纬度、高程,以及同点的北京54(或西安80)坐标系中的直角坐标,即可进行本软件操作了。如某一个控制点的WGS84经度、纬度、高程为: 109度34分28.94343秒, 31度02分25.65526秒, 104.967米该控制点北京54坐标为:

X射线机暴光参数计算法

X射线机曝光参数计算法 基本参数确定 一、以透照厚度为准:单壁单影=T;双壁单影或双壁双影=2T 1、≤10mm时,1mm相当于5KV; 2、10~20mm时,1mm相当于6.2KV; 3、21~30 mm时,1mm相当于9KV; 4、31~40 mm时,1mm相当于12KV; 二、焦距 焦距每增加或者减少100mm,电压增大或者减少10KV。 三、时间 1分钟=25KV 三、X射线机曝光参数为(基数): 透照厚度T=8mm时,电压170KV,时间为1分钟。 四、X射线机焦点到窗口的距离 XXQ 2005 120 mm XXQ 2505 150 mm XXQ 3005 170 mm 五、计算方法 1、当透照厚度增加或者减少1 mm时,电压变化按(一)中各变化范围执行; 2、当焦距每增加或者减少100mm时,压变化按(二)中执行; 3、时间每增加或者减少1分钟,电压增加或者减少25KV; 例:计算φ219*14管焊口的曝光 第一步:确定所用X射线机型号,XXQ 2505或者XXQ 3005型; 第二步:计算焦距-----219+150=369 mm或者219+170=389 mm 第三步:确定焦距和电压变化量,我们一般以X射线机曝光正常基数为准,即600 mm;这里φ219*14的焦距为219+150=369 mm或者219+170=389 mm,比基数600 mm缩短231 mm或者211 mm,那么电压就应该减去23.1KV或者21.1KV。 第四步:计算透照厚度变化时,电压变化量,我们基本厚度是8 mm,现在透照厚度是 14×2=28 mm。这样比基本厚度8 mm增加20mm,根据(一)中4参照,电压补偿量为: 20 mm×8KV=160KV。因为基数是170KV,故正常曝光参数为:170KV+160KV-23.1KV=306.9KV 或者170KV+160KV-21.1KV=308.9KV,时间1分钟。 第五步:因为1分钟=25KV,在此基础上计算XXQ 2505或者XXQ 3005型的曝光参数: 1、XXQ 2505:用240KV拍片,其时间为(306.9 KV-240 KV)÷25KV/分钟=2.68 分钟;这里2.68分钟是在原来1分钟基础需要补偿的2.68分钟,故还应加上基础1分钟, 即正常曝光时间为2.68分钟+1分钟≈4分钟

烟气监测系统计算公式

烟气监测系统计算公式: 1. 流量 1.1原烟气流量(湿态) 【未用】 1.2净烟气流量 1.2.1工况下的湿烟气流量s Q : s s V F Q ??=3600 s Q ――工况下的湿烟气流量,h m 3; F ――监测孔处烟道截面积,2m ; s V ――监测孔处湿烟气平均流速,s m /。 1.2.2监测孔处湿烟气平均流速s V : s V = 流速仪输出值 1.2.3标准状态下干烟气流量sn Q : )1(273273101325sw s s a s sn X t P B Q Q -+?+?= sn Q ――标准状态下干烟气流量,m 3; sw X ――烟气湿度。 1.2.4烟气排放量 ∑=?=n i sni h Q n Q 1)1( ∑==24 1i hi d Q Q ∑==31 1i di m Q Q ∑==121i mi y Q Q 式中, Q h ——标准状况下干烟气小时排放量,m 3;

Q d ——标准状况下干烟气天排放量,m 3; Q m ——标准状况下干烟气月排放量,m 3; Q y ——标准状况下干烟气年排放量,m 3; Q sni ——标准状况下,第i 次采样测得的干烟气流量,m 3/h ; Q hi ——标准状况下,第i 个小时的干烟气小时排放量,m 3/h ; Q di ——标准状况下,第i 天的干烟气天排放量,m 3/h ; Q mi ——标准状况下,第i 个月的干烟气月排放量,m 3/h ; n ——每小时内的采样次数。 2.烟气湿度sw X : 222O O O sw X X X X '-'= 2O X ――湿烟气氧量,%; 2O X '――干烟气氧量,%。 3.过量空气系数α': 2 2121O X -='α 4.烟尘 4.1.1标准状态下干烟气的烟尘排放浓度 程截距烟尘方程斜率+烟尘方.dust dust C C ''=' 式中, dust C ''——实测的烟尘排放浓度,mg/m 3; dust C '——标准状态下干烟气烟尘排放浓度,mg/m 3。 4.1.2折算的烟尘排放浓度 α α'?'=dust dust C C 式中, dust C ——折算成过量空气系数为α时的烟尘排放浓度; dust C '——标准状态下干烟气烟尘排放浓度,mg/m 3; α' ——实测的过量空气系数;

坐标转换三参数计算器使用说明

坐标转换三参数计算器使用说明 一、软件功能 该软件可实现在北京54坐标系、西安80坐标系、WGS84坐标系(GPS通常采用WGS84坐标系)之间进行三参数条件下的高精度相互转换,求取手持GPS 的北京54(或西安80)DA、DF、DX、DY、DZ坐标转换的参数。 二、使用说明 软件分成上下二部分,上半部为在两个不同椭球体间求坐标转换的三参数DX、DY、DZ,下半部为在两个不同椭球体间的坐标转换(如下图)。 在两个不同椭球体间进行坐标转换首要条件是必需知道坐标转换参数,通常有三参数和七参数转换二种方式,本程序提供三参数转换方式。 实例1:我要求手持GPS的北京54(或西安80)坐标转换参数。 向有关部门收集所在工作区内已知点(只要一个控制点)的WGS84坐标系经纬度坐标,以及同点的北京54(或西安80)坐标系中的直角坐标,即可进行本软件操作了。如某一个控制点的WGS84经度、纬度、高程为: 109度34分28.94343秒, 31度02分25.65526秒, 104.967米,该控制点北京54坐标为:x=3436391.566m,y=37363926.964m(37为带号),h=108.717m ,将上述数据输入在软件上半部相应栏中,注意勾选前后坐标系正确(坐标系A,坐标系B),

输入中央经线(37带,输111),点击参数计算,计算结果为 DA=-108,DF=0.00000048,dx=32.284979,dy=-90.792978,dz=-57.993043, 此参数即为手持GPS北京54坐标参数。此三参数为不同椭球体间进行坐标转换奠定了基础。以上计算是精确算法,不存在漏洞。 如果收集控制点确实很困难,在不严谨的情况下,用手持GPS在工作区内某点上设置在WGS84状态下长时间观察读数,取平均值,获取WGS84经度、纬度、高程。北京54(或西安80)坐标你再想办法得到(因为你那已经有测量成果了就好说,如果还没开展测量的话,你就得在大比例尺图上读坐标,越精确越好),也能解决问题,但这个办法不推荐使用,你把求得的参数在其它地貌特征点上检验一下是否提高了定点精度,没提高的话,请重复几次,直到符合定点精度要求。 以上方法求得的坐标转换参数为北京54坐标系、西安80坐标系、WGS84坐标系之间相互转换提供了基础,请注意不同地区参数是不一样的。 实例2:如何将WGS84坐标转换为北京54坐标 已知某点WGS84坐标经纬度、高程(GPS通常采用WGS84坐标系)为: 113度12分34.5678秒, 34度56分12.3456秒, 123.888米,已知WGS84坐标转换为北京54坐标三参数为dx=32.284979,dy=-90.792978,dz=-57.993043。输入软件下半部相应栏中,中央经线111输入右上角相应栏中,点击单点转换,北京54坐标结果为X=3869865.711m, Y=19701880.461m(19带),H=127.052m

坐标转换问题

坐标转换问题 坐标转换问题的详细了解对于测量很重要,那么请和我一起来讨论这个问题。 首先,我们要弄清楚几种坐标表示方法。大致有三种坐标表示方法:经纬度和高程,空间直角坐标,平面坐标和高程。我们通常说的WGS-84坐标是经纬度和高程这一种,北京54坐标是平面坐标和高程着一种。 现在,再搞清楚转换的严密性问题,在同一个椭球里的转换都是严密的,而在不同的椭球之间的转换这时不严密的。举个例子,在WGS-84坐标和北京54坐标之间是不存在一套转换参数可以全国通用的,在每个地方会不一样,因为它们是两个不同的椭球基准。 那么,两个椭球间的坐标转换应该是怎样的呢?一般而言比较严密的是用七参数法,即X平移,Y平移,Z平移,X旋转,Y旋转,Z旋转,尺度变化K。要求得七参数就需要在一个地区需要3个以上的已知点,如果区域范围不大,最远点间的距离不大于30Km(经验值),这可以用三参数,即X平移,Y平移,Z平移,而将X旋转,Y旋转,Z旋转,尺度变化K 视为0,所以三参数只是七参数的一种特例。在本软件中提供了计算三参数、七参数的功能。 在一个椭球的不同坐标系中转换需要用到四参数转换,举个例子,在深圳既有北京54坐标又有深圳坐标,在这两种坐标之间转换就用到四参数,计算四参数需要两个已知点。本软件提供计算四参数的功能。 现在举个例子说明:在珠江有一个测区,需要完成WGS-84坐标到珠江坐标系(54椭球)的坐标转换,整个转换过程是这样的:

COORD软件使用说明: COORD软件采用文件化管理,用户可以将一种转换作为一个文件保存下来,下次使用时从文件菜单中选择打开这个文件来调用所有已有的转换参数。 实例一: 转换要求: 用户在一个佛山测区内使用RTK GPS接收机接受了一些点的WGS-84的坐标,现在希望将其转换为北京54和佛山坐标系下的坐标。用户有佛山测区的一些控制点,这些控制点有WGS-84坐标,也有北京-54坐标也有佛山坐标。 分析: WGS-84坐标和北京54坐标是不同两个椭球的坐标转换,所以要求得三参数或七参数,而北京54和佛山坐标都是同一个椭球,所以他们之间的转换是地方坐标转换,需要求得地

烟气流量计算公式

锅炉烟尘测试方法 1991—09—14发布1992—08—01实施 国家技术监督局 国家环境保护局发布 1、主题内容与适用范围 本标准规定了锅炉出口原始烟尘浓度、锅炉烟尘排放浓度、烟气黑度及有关参数的测试方法。 本标准适用于GBl3271有关参数的测试。 2、引用标准 GB l0180 工业锅炉热工测试规范 GB l327l 工业锅炉排放标准 3、测定的基本要求 3.1 新设计、研制的锅炉在按GBl0180标准进行热工试验的同时,测定锅炉出口原始烟尘浓度和锅炉烟尘排放浓度。 3.2 新锅炉安装后,锅炉出口原始烟尘浓度和烟尘排放浓度的验收测试,应在设计出力下进行。 3.3 在用锅炉烟尘排放浓度的测试,必须在锅炉设计出力70%以上的情况下进行,并按锅炉运行三年内和锅炉运行三年以上两种情况,将不同出力下实测的烟尘排放浓度乘以表l中所列出力影响系数K,作为该锅炉额定出力情况下的烟尘排放浓度,对于手烧炉应在不低于两个加煤周期的时间内测定。 表1 锅炉实测出力占锅炉设计出力的百分数,% 70-《75 75-《80 80-《85 85-《90 9 0-《95 》=95 运行三年内的出力影响系数K 1.6 1.4 1.2 1.1 1.05 1 运行三年以上的出力影响系数K 1.3 1.2 1.1 1 1 1 3.4 测定位置: 测定位置应尽量选择在垂直管段,并不宜靠近管道弯头及断面形状急剧变化的部位。测定位置应距弯头、接头、阀门和其他变径管的下游方向大于6倍直径处,和距上述部位的上游方向大于3倍直径处。 3.5 测孔规格: 在选定的测定位置上开测孔,在孔口接上直径dn为75mm,长度为30mm左右的短管,并装上丝堵。 3.6 测点位置、数目: 3.6.1 圆形断面:将管道断面划分为适当数量的等面积同心圆环,各测点均在环的等面积中心线上,所分的等面积圆环数由管道直径大小而定,并按表2确定环数和测点数。 表2 圆形管道分环及测点数的确定 管道直径D,mm 环数测点数 《200 1 2 200-400 1-2 2-4 400-600 2-3 4-6 600-800 3-4 6-8 800以上4-5 8-10

齿轮各参数计算方法

齿轮各参数计算方法 1、齿数Z 闭式齿轮传动一般转速较高,为了提高传动的平稳性,减小冲击振动,以齿数多一些为好,小一些为好,小齿轮的齿数可取为z1=20~40。开式(半开式)齿轮传动,由于轮齿主要为磨损失效,为使齿轮不致过小,故小齿轮不亦选用过多的齿数,一般可取z1=17~20。为使齿轮免于根切,对于α=20度的标准支持圆柱齿轮,应取z1≥17 2、模数m 齿距与齿数的乘积等于分度圆的周长,即pz=πd。为使d为有理数的条件是 p/π为有理数,称之为模数。即:m=p/π 模数m是决定齿轮尺寸的一个基本参数。齿数相同的齿轮模数大,则其尺寸也大。

3、分度圆直径d 齿轮的轮齿尺寸均以此圆为基准而加以确定,d=mz 4、齿顶圆直径da和齿根圆直径df 由齿顶高、齿根高计算公式可以推出齿顶圆直径和齿根圆直径的计算公式: da=d+2ha df=d-2hf =mz+2m=mz-2×1.25m =m(z+2)=m(z-2.5) 5、分度圆直径d 在齿轮计算中必须规定一个圆作为尺寸计算的基准圆,定义:直径为模数乘以齿数的乘积的圆。实际在齿轮中并不存在,只是一个定义上的圆。其直径和半径分别用d和r表示,值只和模数和齿数的乘积有关,模数为端面模数。与变位系数无关。标准齿轮中为槽宽和齿厚相等的那个圆(不考虑齿侧间隙)就为分度圆。标准齿轮传动中和节圆重合。但若是变位齿轮中,分度圆上齿槽和齿厚将不再相等。若为变位齿轮传动中高变位齿轮传动分度圆仍和节圆重合。但角变位的齿轮传动将分度圆和节圆分离。 6、压力角αrb=rcosα=1/2mzcosα 在两齿轮节圆相切点P处,两齿廓曲线的公法线(即齿廓的受力方向)与两节圆的公切线(即P点处的瞬时运动方向)所夹的锐角称为压力角,也称啮合角。对单个齿轮即为齿形角。标准齿轮的压力角一般为20”。在某些场合也有采用α=14.5°、15°、22.50°及25°等情况。

坐标系转换问题及转换参数的计算方法

坐标系转换问题及转换参数的计算方法 对于坐标系的转换,给很多GPS的使用者造成一些迷惑,尤其是对于刚刚接触的人,搞不明白到底是怎么一回事。我对坐标系的转换问题,也是一知半解,对于没学过测量专业的人来说,各种参数的搞来搞去实在让人迷糊。在我有限的理解范围内,我想在这里简单介绍一下,主要是抛砖引玉,希望能引出更多的高手来指点迷津。 我们常见的坐标转换问题,多数为WGS84转换成北京54或西安80坐标系。其中WGS84坐标系属于大地坐标,就是我们常说的经纬度坐标,而北京54或者西安80属于平面直角坐标。对于什么是大地坐标,什么是平面直角坐标,以及他们如何建立,我们可以另外讨论。这里不多罗嗦。 那么,为什么要做这样的坐标转换呢? 因为GPS卫星星历是以WGS84坐标系为根据而建立的,我国目前应用的地形图却属于1954年北京坐标系或1980年国家大地坐标系;因为不同坐标系之间存在着平移和旋转关系(WGS84坐标系与我国应用的坐标系之间的误差约为80),所以在我国应用GPS进行绝对定位必须进行坐标转换,转换后的绝对定位精度可由80提高到5-10米。简单的来说,就一句话,减小误差,提高精度。 下面要说到的,才是我们要讨论的根本问题:如何在WGS84坐标系和北京54坐标系之间进行转换。 说到坐标系转换,还要罗嗦两句,就是上面提到过的椭球模型。我们都知道,地球是一个近似的椭球体。因此为了研究方便,科学家们根据各自的理论建立了不同的椭球模型来模拟地球的形状。而且我们刚才讨论了半天的各种坐标系也是建立在这些椭球基准之上的。比如北京54坐标系采用的就是克拉索夫斯基椭球模型。而对应于WGS84坐标系有一个WGS84椭球,其常数采用IUGG第17届大会大地测量常数的推荐值。WGS84椭球两个最常用的几何常数:长半轴:6378137±2(m);扁率:1:298.257223563 之所以说到半长轴和扁率倒数是因为要在不同的坐标系之间转换,就需要转换不同的椭球基准。这就需要两个很重要的转换参数dA、dF。 dA的含义是两个椭球基准之间半长轴的差;dF的含义是两个椭球基准之间扁率倒数的差。 在进行坐标转换时,这两个转换参数是固定的,这里,我们给出在进行84—〉54,84—〉80坐标转换时候的这两个参数如下: WGS84>北京54:DA:-108;DF:0.0000005 WGS84>西安80:DA: -3 ;DF: 0 椭球的基准转换过来了,那么由于建立椭球的原点还是不一致的,还需要在dXdYdZ这三个空间平移参量,来将两个不同的椭球原点重合,这样一来才能使两个坐标系的椭球完全转换过来。而由于各地的地理位置不同,所以在各个地方的这三个坐标轴平移参量也是不同的,因此需要用当地的已知点来计算这三个参数。具体的计算方法是: 第一步:搜集应用区域内GPS“B”级网三个以上网点WGS84坐标系B、L、H值及我国坐标系(BJ54或西安80)B、L、h、x值。(注:B、L、H分别为大地坐标系中的大地纬度、大地经度及大地高,h、x分别为大地坐标系中的高程及高程异常。各参数可以通过各省级测绘局或测绘院具有“A”级、“B”级网的单位获得。) 第二步:计算不同坐标系三维直角坐标值。计算公式如下: X=(N+H)cosBcosL

实时动态(RTK)测量中坐标转换参数计算的几种方法

实时动态(RTK)测量中坐标转换参数计算的几种方法 摘要:RTK所接收到的数据是WGS-84坐标系下的数据,而我们使用的坐标系一般是1954北京坐标系、1980年国家大地坐标系以及一些城市工矿使用的独立坐标,因此,需要将RTK接收到的WGS-84坐标转换成我们工程所使用的坐标系坐标。为此,如何计算这些坐标系统转换参数成为RTK使用过程中的一个非常重要的环节。 关键词:GPS-RTK测量坐标转换 1、RTK技术概述 实时动态(RTK)测量系统,是GPS测量技术与数据传输技术的结合,是GPS测量技术中的一个新突破。GPS测量中,静态、快速静态、动态测量都需要事后进行解算处理才能获得待测点的坐标,而RTK测量实时差分定位是一种能够在野外实时得到厘米级精度的测点坐标。 RTK实时测量技术具有全天候、作业效率高、定位精度高、操作简便等优点,因而得到了广泛的应用,而且技术设备越来越先进与方便。RTK测量系统一般由以下三部分组成:GPS接收设备、数据传输设备、软件系统。数据传输系统由基准站的发射电台与流动站的接收电台组成,它是实现实时动态测量的关键设备。 2、RTK实时测量坐标参数转换 RTK所接收到的数据是WGS-84坐标系下的数据,而我们一般使用的坐标系是1954北京坐标系、1980年国家大地坐标系以及一些城市工矿使用的独立坐标,因此,需要将RTK接收到的WGS-84坐标转换成我们使用的1954北京坐标系坐标或1980年国家大地坐标系坐标或城市工矿使用的独立坐标系坐标。为此,如何计算坐标系统转换参数成为RTK使用过程中的很重要的一个环节。 根据RTK的原理,参考站和流动站直接采集的都为WGS84坐标,参考站一般以一个WGS84坐标作为起始值来发射,实时地计算点位误差并由电台发射出去,流动站同步接收WGS84坐标并通过电台来接收参考站的数据,条件满足后就可达到固定解,流动站就可实时得到高精度的相对于参考站的WGS84三维坐标,这样就保证了参考站与流动站之间的测量精度。如果要符合到已有的已知点上,需要把原坐标系统和现有坐标系统之间的转换参数求出。 3、三参数转换

第二章 烟气参数的测定

2.仪器 ①标准型皮托管。标准型皮托管的构造如图5-2-7所示。它是一个弯成90°的双层同心圆管,前端呈半圆形,正前方有一开孔,与内管相通,用来测定全压。在距前端6倍直径出外管壁上开有一圈孔径为1mm 的小孔,通至后端的侧出口,用于测定排气静压。 按照上述尺寸制作的皮托管其修正系数为1.99 ±0.01,如果未经标定,使用时可取修正系数K p 为0.99。 标准型皮托管的侧孔很小当烟道内颗粒物浓度大时,易被堵塞。它是用于测量较清洁的排气。 ②S 型皮托管。S 型皮托管的结构见图5-2-8.它是由两根相同的金属管并联组成。测量端有方向相反的两个开口,测定时,面向气流的开口测得的压力为全压,背向气流的开口测得的压力小于静压。按照图5-2-8设计要求制作的S 型皮托管,其修正系数K p 为0.84 ±0.01。制作尺寸与上述要求有差别S 型皮托管的修正系数需进行校正。其正,反方向的修正系数相差应不大于0.01。S 型皮托管的测压孔开口较大,不易被颗粒物堵塞额,且便于在厚壁烟道中使用。 S 型皮托管在使用前用标准皮托管在风洞中进行校正。S 型皮托管的速度校正系数按下式计算: PS K K = 式中:P S K 、P N K ——分别为标准皮托管和S 型皮托管的速度校正系数; d N P 、d S P ——分别为标准皮托管和S 型皮托管测得的动压值,Pa 。 ③U 形压力计。U 形压力计用于测定排气的全压和静压,其最小分度值应不大于10Pa 。压力计由U 形玻璃管制成,内装测压也挺i ,常用测压液体有水,乙醇和汞,视被测压力范围选用。压力P 按下式计算: P g h ρ=?? 式中:P ——压力,Pa ; h ——液柱差,mm ; ρ——液体密度,g/cm 3 ; 在实际工作中,常用mmH 2O 表示压力,这样压力P=ρ*h U 形压力计的误差较大,不适宜测量微小压力。 ④斜管微压计,斜管微压计用于测定排气的动压,测量范围0~2000Pa ,其精确度应不低于2%,最小分度值应不大于2Pa 。 斜管微压计,构造见示意图5-2-9.一端为截面积较大的容器,另一端为可调角度的玻璃管,管上刻度表示压力的读数。测压时,将微压计容器开口与测定系统中压力较高的一端相连,斜管与系统中压力较低的一端相连,作用于两个液面上的压力差,使液柱沿斜管上升,压力P 按下式计算: P=12sin S L g S αρ? ? ?+ ? ??? 令K=12sin S g S αρ? ? + ? ?? ?

手持GPS参数设置及全国各地坐标转换参数17597

如何设置手持GPS相关参数及全国各地坐标转换参数、如何设置手持GPS相关参数 (一)手持GPS的主要功能 手持GPS,指全球移动定位系统,是以移动互联网为支撑、以GPS 智能手机为终端的GIS系统,是继桌面Gis、WebGis之后又一新的技术热点。目前功能最强的手持GPS,其集成GPRS通讯、蓝牙技术、数码相机、麦克风、海量数据存储、USB/RS232端口于一身,能全面满足您的使用需求。 主要功能:移动GIS数据采集、野外制图、航点存储坐标、计算长度、面积角度(测量经纬度,海拔高度)等各种野外数据测量;有些具有双坐标系一键转换功能;有些内置全国交通详图,配各地区地理详图,详细至乡镇村落,可升级细化。 (二)手持GPS的技术参数因为GPS卫星星历是以WGS84大地坐标系为根据建立的,手 持GPS单点定位的坐标属于WGS84大地坐标系。WGS84坐标系 所采用的椭球基本常数为:地球长半轴a=6378137m;扁率F=1 / 298.257223563。 常用的北京54、西安80及国家2000公里网坐标系,属于平面 高斯投影坐标系统。北京54坐标系,采用的参考椭球是克拉索夫 斯基椭球,该椭球的参数为:地球长半轴a=6378245m;扁率F=1

/298.2。西安 80坐标系,其椭球的参数为:地球长半轴 a=6378140m ;扁率F=1 /298.257。国家2000坐标系,其椭球的参 数为:地球长半轴 a=6378137m ;扁率 F=1 /298.298.257222101。 (三)手持GPS 的参数设置 要想测量点位的北京 54、西安80及国家2000公里网高精度坐 标数据,必须学习坐标转换的基础知识,并分别科学设置手持 GPS 的各项参数。 首先,在手持式GPS 接收机应用的区域内(该区域不宜过大), 从当地测绘部门收集 1至两个已知点的北京 54、西安80或国家 2000坐标系统的坐标值;然后在对应的点位上读取 WGS84坐标 系的坐标值;之后采用《万能坐标转换》软件,可计算出 DY 、DZ 的值。 将计算出的DX 、DY 、DZ 三个参数与DA 、DF 、中 投影比例、东西偏差、南北偏差等六个常数值输入 GPS 接收机。 将GPS 接收机的网格转换为“ UserGrid ”格式,实际测量已知点 的公里网纵、横坐标值,并与对应的公里网纵、横坐标已知值进 行比较, 二者相差较大时要重新计算或查找出现问题的原因。 细过程可查看《万能坐标转换》软件的【手持 GPS 参数设置】界 面。 (四)自定义坐标系统(User )投影参数的确定 DX 、 央经线、

电机参数计算方法

我设定的自制马达规格如左:使用7.4V 1600mA锂电池,耗电在7A以内(马达功率约50W,电池放电系数约4.4C),采用直驱或减速皆可。 以上述条件,无刷马达应采用△接线铜损较小(因线电流=√3*相电流,故马达内线圈电流会较小,以相同的线径来说,铜损自然较小)。 我是采用AWG #28号线(直径0.32mm),每相每极绕21圈,采用△接线,使用7.4V 1600mA 锂电池。 以直驱测试,其数据如下: 螺旋桨测量转数(RPM) 测量电池电流(A) 测量马达线电流(A) 换算马达相电流(A) 计算功率(W) 4040 15000 6.2A 3.6A 2.1A 45W 5025 13000 7.4A 4.3A 2.5A 55W 以减速组测试(58/18=3.2),其数据如下: 螺旋桨测量螺旋桨转数(RPM) 换算马达转速(RPM) 测量电池电流(A) 计算功率(W) 7060 6250 20000 4.2A 31W 8060 5500 17600 6.2A 46W 9070 5000 16000 7.4A 55W 无刷马达/有碳刷马达效能计算 扭力常数: Kt=Kb x 1.345 Kt=1345 / kv 消耗电流: I = [V-(Kb x kRPM)] / Rm I = [V-(RPM / kv)] / Rm 输出扭力: J = (Kt x I) - (Kt x Inl) 每分钟转速: kRPM = (V - RmI) / Kb kRPM = (V - RmI) x kv / 1000 输出功率: Po = (J x RPM) / 1345 消耗功率: Pi = V x I 马达效率: Eff = (Po / Pi) x 100 最高效率电流: Ie max = Sqrt [(V x Inl) / Rm] 符号定义: Eff = 效率 I = 消耗电流值 Iemax=发挥最高效率之电流量 Inl = 无负载量测电流值 J = 扭力(oz-in) Kb = 电压常数(Volt / 1000 RPM) Kt = 扭力常数(oz-In / A) Pi = 消耗功率(Watts) Po = 机械输出功率(Watts) Rm = 马达内阻 RPM = 每分钟转速 V = 电压

烟气折算公式

烟气折算公式 流速 Vs = Kv * Vp 其中 Vs 为折算流速 Kv为速度场系数 Vp 为测量流速 粉尘 1 粉尘干基值 DustG = Dust / ( 1 – Xsw / 100 ) 其中 DustG 为粉尘干基值 Dust 为实测的粉尘浓度值 Xsw 为湿度 2 粉尘折算 DustZ = DustG * Coef 其中 DustZ 为折算的粉尘浓度值 DustG 为粉尘干基值 Coef 为折算系数,它的计算方式如下: Coef = 21 / ( 21 - O2 ) / Alphas 其中 O2 为实测的氧气体积百分比。 Alphas 为过量空气系数(燃煤锅炉小于等于45.5MW折算系数为1.8; 燃煤锅炉大于45.5MW折算系数为1.4; 燃气、燃油锅炉折算系数为1.2) 3粉尘排放率 DustP = DustG * Qs / 1000000 其中 DustP 为粉尘排放率 Dust 为粉尘干基值 Qs 为湿烟气流量,它的计算方式如下: Qs = 3600 * F * Vs 其中 Qs 为湿烟气流量 F 为测量断面面积 Vs 为折算流速 SO2 1 SO2干基值 SO2G = SO2 / ( 1 – Xsw / 100 ) 其中 SO2G 为SO2干基值 SO2 为实测SO2浓度值 Xsw 为湿度 2 SO2折算

SO2Z = SO2G * Coef 其中 SO2Z 为SO2折算率 SO2G 为SO2干基值 Coef 为折算系数,具体见粉尘折算 3 SO2排放率 SO2P = SO2G * Qsn / 1000000 其中 SO2P 为SO2排放率 SO2G 为SO2干基值 Qsn 为干烟气流量,它的计算方式如下: Qsn = Qs * 273 / ( 273 + Ts ) * ( Ba + Ps ) / 101325 * ( 1 – Xsw / 100 ) 其中 Qs 为湿烟气流量 Ts 为实测温度 Ba 为大气压力 Ps 为烟气压力 Xsw 为湿度 NO 1 NO干基值 NOG = NO / ( 1 – Xsw / 100 ) 其中 NOG 为NO干基值 NO 为实测NO浓度值 Xsw 为湿度 2 NO折算 NOZ = NOG * Coef 其中 NOZ 为NO折算率 NOG 为NO干基值 Coef 为折算系数,具体见粉尘折算 3 NO排放率 NOP = NOG * Qsn / 1000000 其中 NOP 为NO排放率 NOG 为NO干基值 Qsn 为干烟气流量,它的计算方式如下: Qsn = Qs * 273 / ( 273 + Ts ) * ( Ba + Ps ) / 101325 * ( 1 – Xsw / 100 ) 其中 Qs 为湿烟气流量 Ts 为实测温度 Ba 为大气压力 Ps 为烟气压力 Xsw 为湿度

相关文档
最新文档