TA的每日心情 | 无聊 2020-11-4 11:03 |
---|
签到天数: 3 天 连续签到: 1 天 [LV.2]偶尔看看I 累计签到:3 天 连续签到:1 天
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
本帖最后由 yuanzhch 于 2011-6-1 20:23 编辑
8 m' |; {, s" ^1 g" p
- X/ H8 a6 u4 L& a/ ?在pscad中要实现新的算法,首先要建立一个自定义元件。自定义元件分为两种,一种是电气元件,一种是信号处理元件。本帖实现的内容是一个非线性电感,其电感值随着电流大小而变化,属于电气元件。写控制算法的同学可能主要用到后者,将输入信号经过一系列的处理计算得到一个输出信号。 " F1 S7 i/ m( \% ?
问题1:将程序写在哪里?
7 ^# `: \9 o. _% e2 c. {# N0 ` 首先要将算法程序分成两个部分:* m H7 U$ Q: `7 C6 H7 p5 ]% }
一部分是PSCAD可以实现的,如处理仿真参数、读入仿真数据、输出仿真结果以及积分、微分等简单函数。这部分内容写在自定义元件definition的Script中。可以写在Fortran中也可以写在DSDYN中,两者略有不同但不会影响简单应用,具体内容请参阅帮助文件。
; ]( V" N/ k! @& }8 ?6 a- h 另一部分是我们自己的C语言算法,要将C算法总结成y=f(x,p)的输入输出形式的形式写在一个*.c文件中。并添加到PSCAD的File Reference中。 C函数的命名形式一定是 void XX_XXX__()的形式,注意名字中的下划线,本人的试验结果是必须有这三个下划线,哪位同学有其他发现可以一起探讨一下。+ w; l. A* T. z6 p; E
问题2:如何传数据?
/ }4 ?3 z; K* W3 d8 i 由于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程序中的变量,后面参数前面加“$”是引用参数的意思,具体内容请参阅帮助文件。+ t& G/ l% z* n5 @/ A, _" x
问题3:如何保存前一时刻的变量并在本时刻调用?2 q$ ^% V4 K' r) H J& H
在时域仿真的时候我们不可避免的要用到上一时刻的某些计算结果。而PSCAD中的Fortran程序和C程序每一时刻都重新执行一次,变量都要清零,所以就要单独存放需要的数据。PSCAD提供了一个STORF数组,这里直接给出例子:
( J0 C$ |# c9 ?/ e. s! ]) {, v#STORAGE REAL:4) `1 m* _ O7 x
MY_NSTORF=NSTORF Z4 q5 }8 k6 R5 W: p7 w
NSTORF=NSTORF+4* ]( ?% b* E/ S4 f+ h4 v
BT=STORF(MY_NSTORF)
, @) s& J& r4 h6 H8 }, D9 @; r. M0 A HT=STORF(MY_NSTORF+1)
* H$ S; S$ h$ `2 c2 ^: Q! {1 V MT=STORF(MY_NSTORF+2)7 U a8 h9 B2 }3 @
H=STORF(MY_NSTORF+3)*$N/$L
( G# P5 C4 A9 F# ]$ A这是读上一时刻值。
6 R* i4 K2 g7 O STORF(MY_NSTORF)=B
' g! s m5 Z* d: H- T) n" X0 @ STORF(MY_NSTORF+2)=M
1 w, c Q! G4 ]这是保存存本时刻值。
7 K4 |: B$ R# ]4 {) |解决了上面3个问题就可以完成信号处理元件的建模。
. t- x) _, S4 m* G( K问题4:如何建立电气元件?2 P: y6 N1 ]! y9 Z$ K/ V
我们对一个电气元件建模就是要描述其伏安特性。如一个双端口元件,测它两个结点电压,然后计算出支路电流就可以了。PSCAD提供了一个基于结点的电气接口函数,说白了就是一个受控电流源叫GGIN和CCIN。帮助里可以搜到,下面给出例子:
2 ?" c3 f2 `5 {# [8 h CCIN($NA,$SS)=CCIN($NA,$SS)-STORF(MY_NSTORF+3)/1000.0
0 t3 w! l% u+ j CCIN($NB,$SS)=CCIN($NB,$SS)+STORF(MY_NSTORF+3)/1000.0
' J7 [2 u9 \. ]& l GGIN($NA,$SS)=GGIN($NA,$SS)+G
. L- m! t* B: E# O7 P9 R+ Z, e+ N GGIN($NB,$SS)=GGIN($NB,$SS)+G
: T7 q. s1 f" s z& E我觉得解决以上4个问题就能满足一般应用了。这里只是给大家提供一个思路,还有一些细节没有说到还需要看帮助文件,希望能帮到大家!附件中是个非线性电感,J-A模型是个微分方程在C函数中用龙格库塔法解之。
1 M2 `3 ~% E& g, ]& i
基于JA理论的非线性电感.rar
(1.03 KB, 下载次数: 241)
|
评分
-
查看全部评分
|