TA的每日心情 | 无聊
2020-11-4 11:03 |
|---|
签到天数: 3 天 连续签到: 1 天 [LV.2]偶尔看看I 累计签到：3 天
连续签到：1 天
|
马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！
您需要 登录 才可以下载或查看，没有账号？立即加入
×
本帖最后由 yuanzhch 于 2011-6-1 20:23 编辑
; p. O& b& n7 f; J$ _2 S% B; G* v' X' B# C% F: r) e0 O! P" Q; `& a
在pscad中要实现新的算法，首先要建立一个自定义元件。自定义元件分为两种，一种是电气元件，一种是信号处理元件。本帖实现的内容是一个非线性电感，其电感值随着电流大小而变化，属于电气元件。写控制算法的同学可能主要用到后者，将输入信号经过一系列的处理计算得到一个输出信号。 % @" a2 p% M, V, z3 I# H# Z( Y
问题1：将程序写在哪里？
; z7 K- F+ U4 ]1 ^& d 首先要将算法程序分成两个部分：9 B& G' m4 L2 `' L! X8 [" b6 j) |
一部分是PSCAD可以实现的，如处理仿真参数、读入仿真数据、输出仿真结果以及积分、微分等简单函数。这部分内容写在自定义元件definition的Script中。可以写在Fortran中也可以写在DSDYN中，两者略有不同但不会影响简单应用，具体内容请参阅帮助文件。
6 J: N9 f2 u. q6 U* }3 w: B 另一部分是我们自己的C语言算法，要将C算法总结成y=f（x，p）的输入输出形式的形式写在一个*.c文件中。并添加到PSCAD的File Reference中。 C函数的命名形式一定是 void XX_XXX__（）的形式，注意名字中的下划线，本人的试验结果是必须有这三个下划线，哪位同学有其他发现可以一起探讨一下。" \7 j' O+ |7 n' r7 l' ^, W& Q
问题2：如何传数据？8 F, Q. R' S4 p" J
由于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程序中的变量，后面参数前面加“$”是引用参数的意思，具体内容请参阅帮助文件。/ v% ?8 o3 U. C
问题3：如何保存前一时刻的变量并在本时刻调用？( u) k& \# z; ]2 u$ p) @
在时域仿真的时候我们不可避免的要用到上一时刻的某些计算结果。而PSCAD中的Fortran程序和C程序每一时刻都重新执行一次，变量都要清零，所以就要单独存放需要的数据。PSCAD提供了一个STORF数组，这里直接给出例子：
) {- n: i# U! V#STORAGE REAL:45 p1 p f; W( }/ [) I6 P+ J
MY_NSTORF=NSTORF% q8 `* B4 u ~! Z5 i9 X U
NSTORF=NSTORF+4
% A& j4 S3 U" d3 S& D9 _ BT=STORF(MY_NSTORF)
, _; ], {% c% s0 U HT=STORF(MY_NSTORF+1)
' f/ v9 H) j" h/ T MT=STORF(MY_NSTORF+2): y9 }# L5 J$ k& P; ~. H/ g# t
H=STORF(MY_NSTORF+3)*$N/$L" w% S! w \6 t- B
这是读上一时刻值。
$ ~ o- i. k6 b) z$ p% ?( l STORF(MY_NSTORF)=B
6 L ]; l3 e3 t+ G& P STORF(MY_NSTORF+2)=M; `4 q7 L. h4 ^% i! e. U# z
这是保存存本时刻值。
8 w1 [& w& @* U解决了上面3个问题就可以完成信号处理元件的建模。
- P5 T$ n% [: [/ E8 [& a7 r4 B问题4：如何建立电气元件？
" t% p5 d9 G8 u$ @: h 我们对一个电气元件建模就是要描述其伏安特性。如一个双端口元件，测它两个结点电压，然后计算出支路电流就可以了。PSCAD提供了一个基于结点的电气接口函数，说白了就是一个受控电流源叫GGIN和CCIN。帮助里可以搜到，下面给出例子：. M, x# s2 N; F6 Y3 L2 j: d
CCIN($NA,$SS)=CCIN($NA,$SS)-STORF(MY_NSTORF+3)/1000.0
/ R; @4 h; I& N8 S CCIN($NB,$SS)=CCIN($NB,$SS)+STORF(MY_NSTORF+3)/1000.0' O$ ]* W" \$ Q0 g' ^# Q
GGIN($NA,$SS)=GGIN($NA,$SS)+G
$ r2 k4 \* e$ H/ L GGIN($NB,$SS)=GGIN($NB,$SS)+G
. M8 O) {: l6 M% t" K* H' K6 a我觉得解决以上4个问题就能满足一般应用了。这里只是给大家提供一个思路，还有一些细节没有说到还需要看帮助文件，希望能帮到大家！附件中是个非线性电感，J-A模型是个微分方程在C函数中用龙格库塔法解之。* D" f. s: d$ J9 C9 B* a
基于JA理论的非线性电感.rar
(1.03 KB, 下载次数: 241)
|
评分
-
查看全部评分
|
赣公网安备36040302000210号 )|网站地图
狗仔卡
发表于 2011-6-1 20:19:02
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡
发表于 2011-6-2 01:41:43
