TA的每日心情 | 无聊 2020-11-4 11:03 |
---|
签到天数: 3 天 连续签到: 1 天 [LV.2]偶尔看看I 累计签到:3 天 连续签到:1 天
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
本帖最后由 yuanzhch 于 2011-6-1 20:23 编辑
5 B5 q4 N, V) T% Z z
& K/ E# V! @7 b在pscad中要实现新的算法,首先要建立一个自定义元件。自定义元件分为两种,一种是电气元件,一种是信号处理元件。本帖实现的内容是一个非线性电感,其电感值随着电流大小而变化,属于电气元件。写控制算法的同学可能主要用到后者,将输入信号经过一系列的处理计算得到一个输出信号。
. _+ T0 o& v- w$ B8 m' Y3 c问题1:将程序写在哪里?% T/ f+ _" h8 y( w( Z/ t4 B% L
首先要将算法程序分成两个部分:! h+ Q6 N7 I. e8 H% y
一部分是PSCAD可以实现的,如处理仿真参数、读入仿真数据、输出仿真结果以及积分、微分等简单函数。这部分内容写在自定义元件definition的Script中。可以写在Fortran中也可以写在DSDYN中,两者略有不同但不会影响简单应用,具体内容请参阅帮助文件。& ^! k# L! m0 |
另一部分是我们自己的C语言算法,要将C算法总结成y=f(x,p)的输入输出形式的形式写在一个*.c文件中。并添加到PSCAD的File Reference中。 C函数的命名形式一定是 void XX_XXX__()的形式,注意名字中的下划线,本人的试验结果是必须有这三个下划线,哪位同学有其他发现可以一起探讨一下。. S/ e5 A+ C" C
问题2:如何传数据?3 U8 s) B/ V6 F( o
由于PSCAD和C函数是互相传递指针,因此就需要将C函数写成f(y,x,p)的形式,例如:void ja_model__(real *h,real *ht,real *mt,real *ms,real *a,real *alpha,real *c,real *k,real *b,real *m)。其中h是输出,ht,mt是输入,后面的都是参数。在PSCAD中用call语句调用,例如:CALL JA_MODEL(H,HT,MT,$MS,$A,$ALPHA,$C,$K,B,M) 。注意函数名中的下划线,这里只有一个!其中H,HT,MT都是Fortran程序中的变量,后面参数前面加“$”是引用参数的意思,具体内容请参阅帮助文件。/ E3 c% y3 }% u# V( D
问题3:如何保存前一时刻的变量并在本时刻调用?
, |1 Y0 L' n+ }6 y/ [# j% [ 在时域仿真的时候我们不可避免的要用到上一时刻的某些计算结果。而PSCAD中的Fortran程序和C程序每一时刻都重新执行一次,变量都要清零,所以就要单独存放需要的数据。PSCAD提供了一个STORF数组,这里直接给出例子:# z' y4 T7 k: X: |+ R
#STORAGE REAL:43 x/ |7 I9 A* N2 P, X# Z2 T
MY_NSTORF=NSTORF. U, y# Q1 S4 ]4 o0 d# V9 `% F
NSTORF=NSTORF+4
, r7 L7 E! N& a" I BT=STORF(MY_NSTORF)
4 Y9 d7 F7 u- J$ \, y. w0 X" s HT=STORF(MY_NSTORF+1)% D9 c, }$ j0 N9 a0 ^: t
MT=STORF(MY_NSTORF+2). c. Q+ j1 h' c, }
H=STORF(MY_NSTORF+3)*$N/$L0 O1 I0 N( Q+ ^% T: l x! o ^. S
这是读上一时刻值。
! h) ]1 l/ B* ? g: ?- s7 l STORF(MY_NSTORF)=B
- J' D; c) D7 T5 ~5 K STORF(MY_NSTORF+2)=M! U( {0 ~7 ~) t
这是保存存本时刻值。
2 O5 Q+ m2 a" H1 g& j: H- Y解决了上面3个问题就可以完成信号处理元件的建模。9 `: u- G' R3 w# G2 n: U& N7 H }% w
问题4:如何建立电气元件?& T. q" F, F7 V: h$ C! \* m4 Q
我们对一个电气元件建模就是要描述其伏安特性。如一个双端口元件,测它两个结点电压,然后计算出支路电流就可以了。PSCAD提供了一个基于结点的电气接口函数,说白了就是一个受控电流源叫GGIN和CCIN。帮助里可以搜到,下面给出例子:
* q N. t, }, @6 M( N' Q+ C, ? CCIN($NA,$SS)=CCIN($NA,$SS)-STORF(MY_NSTORF+3)/1000.03 j. W! [9 W' L' L6 ?5 W
CCIN($NB,$SS)=CCIN($NB,$SS)+STORF(MY_NSTORF+3)/1000.09 x/ Q4 Q8 r0 h- C' Z- O
GGIN($NA,$SS)=GGIN($NA,$SS)+G6 C$ \* R9 h. i' {2 w$ u
GGIN($NB,$SS)=GGIN($NB,$SS)+G
& k' d# g- Y/ p6 {我觉得解决以上4个问题就能满足一般应用了。这里只是给大家提供一个思路,还有一些细节没有说到还需要看帮助文件,希望能帮到大家!附件中是个非线性电感,J-A模型是个微分方程在C函数中用龙格库塔法解之。
1 a7 Z7 U2 A) b- ?8 t* k( ?
基于JA理论的非线性电感.rar
(1.03 KB, 下载次数: 241)
|
评分
-
查看全部评分
|