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