TA的每日心情 | 无聊
2020-11-4 11:03 |
|---|
签到天数: 3 天 连续签到: 1 天 [LV.2]偶尔看看I 累计签到：3 天
连续签到：1 天
|
马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！
您需要 登录 才可以下载或查看，没有账号？立即加入
×
本帖最后由 yuanzhch 于 2011-6-1 20:23 编辑 / k' C0 f/ i" K* ^ u4 a
* s: z+ G4 d6 r5 X9 t6 y在pscad中要实现新的算法，首先要建立一个自定义元件。自定义元件分为两种，一种是电气元件，一种是信号处理元件。本帖实现的内容是一个非线性电感，其电感值随着电流大小而变化，属于电气元件。写控制算法的同学可能主要用到后者，将输入信号经过一系列的处理计算得到一个输出信号。
$ d/ i" s) j" F. `7 p1 d% G, p问题1：将程序写在哪里？# l3 U4 n% [8 t& U. l
首先要将算法程序分成两个部分：2 v8 f! F. y u' o2 `
一部分是PSCAD可以实现的，如处理仿真参数、读入仿真数据、输出仿真结果以及积分、微分等简单函数。这部分内容写在自定义元件definition的Script中。可以写在Fortran中也可以写在DSDYN中，两者略有不同但不会影响简单应用，具体内容请参阅帮助文件。 _) L" M+ B( [
另一部分是我们自己的C语言算法，要将C算法总结成y=f（x，p）的输入输出形式的形式写在一个*.c文件中。并添加到PSCAD的File Reference中。 C函数的命名形式一定是 void XX_XXX__（）的形式，注意名字中的下划线，本人的试验结果是必须有这三个下划线，哪位同学有其他发现可以一起探讨一下。
. V# Q( O, r5 g/ f问题2：如何传数据？3 `: ?, e* ~; e& A* g! i" 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程序中的变量，后面参数前面加“$”是引用参数的意思，具体内容请参阅帮助文件。
! R4 U. ~% }+ }$ D) h R' _6 t问题3：如何保存前一时刻的变量并在本时刻调用？8 n2 Z: P+ s- m+ W/ k; i! R6 f
在时域仿真的时候我们不可避免的要用到上一时刻的某些计算结果。而PSCAD中的Fortran程序和C程序每一时刻都重新执行一次，变量都要清零，所以就要单独存放需要的数据。PSCAD提供了一个STORF数组，这里直接给出例子：
$ S6 {/ Q0 S5 u- j+ R: L#STORAGE REAL:49 s. b; S, j/ Q) i
MY_NSTORF=NSTORF
: W9 w- D q$ Y+ C1 X* w6 _ NSTORF=NSTORF+4
3 M: E' _7 j/ T3 n- f BT=STORF(MY_NSTORF)
+ Q- C2 m6 m0 ? HT=STORF(MY_NSTORF+1)
! ?0 I2 h+ ]1 Z2 D! m MT=STORF(MY_NSTORF+2)9 Y2 D, s) _4 K" A$ N) }* f% q4 ~
H=STORF(MY_NSTORF+3)*$N/$L+ E) h5 m. n- ?0 w8 `
这是读上一时刻值。1 }, Z0 F6 S& j2 k) v0 ^4 u
STORF(MY_NSTORF)=B- `* b- x s7 Z0 x( @
STORF(MY_NSTORF+2)=M' }/ U M5 D6 W2 @' h; n u
这是保存存本时刻值。
& B2 s X1 y- {1 V+ [解决了上面3个问题就可以完成信号处理元件的建模。
4 O1 z* p1 t, f, A$ x, ?8 E问题4：如何建立电气元件？) W! }5 \/ K1 }* D" g
我们对一个电气元件建模就是要描述其伏安特性。如一个双端口元件，测它两个结点电压，然后计算出支路电流就可以了。PSCAD提供了一个基于结点的电气接口函数，说白了就是一个受控电流源叫GGIN和CCIN。帮助里可以搜到，下面给出例子：
* J3 L2 |/ a" K0 E& W+ V/ d5 { CCIN($NA,$SS)=CCIN($NA,$SS)-STORF(MY_NSTORF+3)/1000.0
$ B- @' ^) ]5 L9 v! J' [ CCIN($NB,$SS)=CCIN($NB,$SS)+STORF(MY_NSTORF+3)/1000.0
8 J1 C$ y8 P0 p4 R. x GGIN($NA,$SS)=GGIN($NA,$SS)+G
; z( O- p7 E" P. R4 y: Z GGIN($NB,$SS)=GGIN($NB,$SS)+G
3 g9 `) a* m1 {3 O, d$ r我觉得解决以上4个问题就能满足一般应用了。这里只是给大家提供一个思路，还有一些细节没有说到还需要看帮助文件，希望能帮到大家！附件中是个非线性电感，J-A模型是个微分方程在C函数中用龙格库塔法解之。
' V4 X$ u3 `* ~3 w- w6 q$ w: l
基于JA理论的非线性电感.rar
(1.03 KB, 下载次数: 238)
|
评分
-
查看全部评分
|
赣公网安备36040302000210号 )|网站地图
狗仔卡
发表于 2011-6-1 20:19:02
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡
发表于 2011-6-2 01:41:43
