1、if( (k2*X2(n) + C2)=0 ) Y(n) = Y(n-1),并且置位OFW和Q; else Y(n) (k1*X1(n)+C1) / (k2*X2(n) + C2);Y(n)被限幅在有效数值之间,如果越限,就置位本功能模块的数值溢出状态位。其它情况下,本功能块按属性中参数所定义品质传递方式传递前续功能块的状态。4 开方器 SQRT ID=4图形表示:X SQRT Y Y=(k*X+C)1/2标记说明:开方计算输出输入 X参数 k输入X的增益X的偏置参数 DB零点切除值, 0.0算法说明:本功能块对输入浮点变量开方,输出一个浮点变量。if( k*X(n)+C DB) Y(n)=;
2、else Y(n)=0,并且置位OFW和Q;5 取绝对值 ABS ID=5 X ABS Y 绝对值计算输出参数 k, C1.0, 0.0输入X的增益和偏置算法描述:本功能块对输入浮点变量取绝对值,输出一个浮点变量。Y(n)=|k*X(n)+C|6 五次多项式 POLYNOM ID=6X POLYNOM Y多项式计算输出参数 C0C5多项式的系数本功能块对浮点变量进行五次多项式运算,输出一个浮点变量。Y(n) C0 + C1*X(n) + C2*X(n)*2 + C3*X(n)*3 +C4*X(n)*4 +C5*X(n)*57 8输入数学统计器 SUM8 ID=7X1 Y SUM8 X8 Cnt
3、数学统计输出输出 Cnt模拟输入量中绝对值DB的个数输入 X1,X2,X8Null参数 Modeword方式,0=累加和,1=取平均,2=取最大,3=取最小参数 k1,k2,k8输入X1X8的增益参数 C1,C2,C8输入X1X8的偏置判非零的死区本功能块对8个浮点变量加或减,输出一个浮点变量。输出Y(n)按方式取输入的累加和、平均、最大、最小。所有输入为k*X(n)+C。空脚不参与运算。输出Cnt(n)等于 fabs(k*X(n)+C)DB的个数。8 12段函数变换F(X) ID=8图形描述:X F(X) YYF(X)标记描述:函数变换输出参数 (x1,y1)0.0,0.0第1个坐标点参数
4、.第n个坐标点, 要求x1x2.x12参数 (x12,y12)第12个坐标点 F(X)由12个坐标点(X1,Y1)(X12,Y12)确定 如 上图所示:F(X)为一个十二段折线,如XX1,YY1;XX12,YY12。 X的值应递增填写,不递增曲线取到递增的转折点为止。9 贝塞尔函数 BESSEL, ID=9X Bessel Y函数计算输出函数类型:0第一类Bessel函数jn;1第二类Bessel函数 yn参数 Rankword 函数阶次,0n本功能块计算Bessel函数。10 指数/对数函数/求模 POW/LOG/MOD ID=10 X POW/LOG YY B指数或对数,以B.I方式表示的
5、指针,或是立即浮点数输入 B10.0底数,以B.I方式表示的指针,或是立即浮点数方式,0=为指数运算,1=对数运算,2=求模运算输入和底变量的增益输入和底变量的偏置本功能块对浮点变量进行指数或数值或求模,输出一个浮点变量。当Mode = 0时: Y(n) (k2*B(n)+C2) * (k1*X(n)+C1),底数位0.0时置OFW和Q标志位。当Mode=1时: Y(n) LOG (k2*B(n)+C2) (k1*X(n)+C1) 当(k2*B(n)+C2为立即浮点数且为0.0时,该函数为自然对数ln。当Mode=2时:Y(n)=(k1*X(n)+C1)MOD (k2*B(n)+C2) 当(k
6、2*B(n)+C2为0.0时,置OFW和Q标志位。求模运算对运算双方自动取整后再运算。11 三角和反三角函数 TRIANGLE ID=11X TRIANGLE Y Y=sin(k*X+C)三角函数计算输出以B.I方式表示的指针,或是立即浮点数,单位弧度函数类型, 0=sin, 1=cos, 2=tan, 3=asin, 4=acos, 5=atan 变量的增益变量的偏置本功能块按Mode对浮点变量进行三角或反三角运算,输出一个浮点变量。角度都以弧度为单位。例:Y(n) sin(k*X(n)+C), Y(n) asin(k*X(n)+C) 12 热力性质计算 PTCal ID=12P/S/H P
7、TCal YT/S/H热力性能计算输出输入 P/S/H压力Mpa/熵值S/焓值H输入 T/S/H温度或熵值S或焓值H计算方式,0=焓值,1=比容,2=熵值,3=温度本功能块用于热力性能计算。在Mode=0时,计算出给定压力(MPa)和温度()的水或汽的焓值(KJ/Kg)在Mode=1时,计算出给定压力(MPa)和温度()的水或汽的比容(M3/Kg)在Mode=2时,计算出给定压力(MPa)和温度()的水或汽的熵值(KJ/Kg*K)在Mode=3时,计算出给定压力(MPa)和熵值(KJ/Kg*K)的水或汽温度()在Mode=4时,计算出给定压力(MPa)和温度()的水的焓值(KJ/Kg)在Mod
8、e=5时,计算出给定压力(MPa)和温度()的水的比容(M3/Kg)在Mode=6时,计算出给定压力(MPa)和温度()的水的熵值(KJ/Kg*K)在Mode=7时,计算出给定熵值(KJ/Kg*K)和温度()的水或汽的压力(MPa)在Mode=8时,计算出给定压力(MPa)和温度()的汽的焓值(KJ/Kg)在Mode=9时,计算出给定压力(MPa)和温度()的汽的比容(M3/Kg)在Mode=10时,计算出给定压力(MPa)和温度()的汽的熵值(KJ/Kg*K)在Mode=11时,计算出给定压力(MPa)或温度()的水的沸点温度或压力 当P为空脚时计算压力;当T为空脚时计算温度.在Mode=1
9、2时,计算出给定压力(MPa)和焓值(KJ/Kg)的水或汽温度()在Mode=13时,计算出给定焓值(KJ/Kg)和温度()的水或汽的压力(MPa) Mode=02 分别计算出给定压力和温度的水或汽的焓值、比容和熵值;气态或液态由给定 的压力和温度自动识别。 Mode=46 分别计算出给定压力和温度的水的焓值、比容和熵值;当只有一个输入而另一 个输入为空脚时,该水缺省为饱和。 Mode=810 分别计算出给定压力和温度的 汽的焓值、比容和熵值;当只有一个输入而另 一个输入为空脚时,该气体缺省为饱和。13 温度场计算模块 FDM, ID=15边界温度TSrnd 换热系数ALFA Tnode20
10、节点温度 应力集中系数StrC FDM 初始温度系数Tinit3 Tavg 体积平均温度复位TReset 计算允许FdmEn Valid 数据有效输出 Tnode20节点温度1-20输出 Tavg体积平均温度输出 Validbool1输出数据是否有效:1有效,0无效输入 TSrnd20.0输入 ALFA输入 StrC输入 Tinit3 20,0,0初始温度系数02输入 TReset所有输出到初始值输入 FdmEn1允许计算;0禁止计算参数 Ro0.4外径参数 Ri0.1内径参数 iR15截面数参数 Rou37800,0,0密度系数02参数 Nmda341.413,-0.0085,-0.0000
11、1导热系数02参数 Cp3flloat554.04-0.0809比热系数02参数 WayWho算法选择:0-Robin;1-DirichletFDM模块用于计算汽轮机变动工况过程中危险截面的温度场,为ATC提供数值依据,算法基于标准的二维轴对称有限差分方法。14 区间代数基本性质 iProp, ID=16区间下限ivaliL M 区间中 W 区间宽区间上限ivaliU IProp S 区间长 Lx 新区间下限宽限Trim Ux 新区间上限输出 M区间中输出 W区间宽输出 S区间长输出 Lx新区间下限输出 Ux新区间上限输入 ival1L输入 ival1U输入 Trim区间的基本性质:中、宽、长
12、、新区间的下限和上限。 M=( ival1L + ival1U)/2; W=( ival1U - ival1L)/2; S=(abs(ival1L)+abs(ival1U)/2; Lx=M-Trim*W; Ux=M+Trim*W.15 区间代数四则运算 iBsc, ID=17区间1下限L1 区间1上限U1 L区间2下限L2 IBsc 区间2上限U2 U指数输入f 输出 L输出 U输入 L1输入 U1输入 L2输入 U2输入 f乘方时为指数,其他情况为比例系数参数 Oplode0: 加“+”; 1: 减“-”; 2: 乘“*”; 3: 除“/”; 4: 乘方“”改进算法; 1:标准算法; 2:改
13、进算法2完成区间的基本四则运算:加、减、乘、除、乘方。第2章 时间过程函数1 超前滞后模块 LEADLAG ID=20X YTR LEADLAG TS 超前滞后输出输入 TR以B.I方式表示的指针,或是立即浮点数。被跟踪量,空脚总不跟踪。输入 TSF以B.I方式表示的指针,或是立即布尔数。跟踪切换开关参数 H100.0输出高限参数 L输出低限参数 LD超前时间常数,如果 LD为0,则为滞后模块参数 LG滞后时间常数,LG必须大于等于1/2T(采样周期)参数 GN增益GNif( TS1 ) Y(n) =TR(n)*GN;else K1=GN(T+2LD)/(2LG+T); K2=GN(T2LD)
14、/(2LG+T); K3=(2LGT)/(2LG+T); Y(n) = K1*X(n)+ K2*X(n1)+ K3*Y(n1);if ( Y(n)H ) Y(n)=H; else if( Y(n)参数 LT惯性时间常数,单位秒,=0.0时为纯滞后本功能块对输入进行纯滞运算,计算公式的拉氏传递函数如下:Y(S) = e*(-DT*S) * k / ( LT*S + 1 )递推计算公式如下: pn=KG*ts*x+lt*tmp-PreAI; /k/(LS+1) pn/=(ts+lt);a = k*T/(T+L), b =L/(T+L), L= int(DT/T), T=计算周期P(n) = a*X
15、(n) + b*X(n-1), Y(n) = P(n-L);因缓冲区长度为120,所以L大于120时,其中有些输出为线性插值的结果。3 微分模块 DIFF ID=22X DIFF Y 微分输出输入参数 Kd增益参数 Td微分时间常数,必须0,单位秒对应的传递函数为: Y(s)Y(n) = ( Kd*Td*(X(n)-X(n-1) + Td*Y(n-1) ) / (T+Td);然后,将Y限制在YH和YL之间。4 时域统计模块 TSUM ID=23 TSUM Set Yj 输出 Yj前一次的统计结果值输入 SetOn/Off开关,以B.I方式表示的指针,或是立即布尔数方式,0=累加和,1=取平均,
16、2=取最大,3=取最小,4=梯形累加和参数 Y0每次On时,Y的初值参数 Cj记忆转换系数本功能块对输入模拟变量在Set指定的一段时间内进行累加,平均,或取最大、最小值、梯形累加和,并记录前次统计值。if( Set从0-1 ) 统计开始,Yj(n)=Y(n)*Cj,Y(n)=Y0, sum=1;else if(Set从1-0) 统计结束 else if( Mode = 0 ) Y(n) = Y(n-1) + X(n)*T; / 累加和 else if( Mode = 1 ) Y(n) = 累计和 / sum; sum=内部计数; / 取平均 else if( Mode = 2 ) Y(n) = Max ( Y(n-1), X(n) ); / 取最大 else if( Mode = 3 ) Y(n) = Min ( Y(n-1), X(n) ); / 取最小 esle if( Mode = 4 ) Y(n) = Y(n-1)+ (X(n)+X(n-1)*T/2; / 梯形累加和其中T为该页计算周期(以秒为单位),因此,累计方式时,输出Y为输入X在计时间内的积分值。5 数字滤波 FILTER ID=24X FILTER Y滤波后的数值参数 k1,.,k81/8滤波器系数,应1.0,八个系数和为1本功能块对输入模拟变量进行8阶数字滤波。Y(n) = k
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1