简介「PSSE自定义模型学习资料」面向电力系统动态仿真工程师、高校电力专业研究者及需要开展PSS/E二次开发的用户聚焦用户自定义模型从编写到接入仿真的关键问题。压缩包内共1个doc文档约27KB篇幅紧凑便于按主题查阅。文档围绕动态仿真标志位MODE、KPAUSE、MSTATE、MIDTRM展开说明STRT等行为如何通过连接子程序调用模型子程序并梳理SUBROUTINE modelname(MC,ISLOT)的书写规范、STRTIN数组位置表入口含义、TBLCNC自动调用机制及模型调用顺序。读者可据此理解自定义模型在动态数据输入文件中的记录格式掌握Fortran模型子程序的逻辑书写步骤与地址分配规则减少联调中的接口错误。目前已有680人学习适合具备一定PSS/E基础、希望深入自定义模型开发与动态仿真排错的读者参考。1. 标志位驱动为什么自定义模型不是被调用一次而是被反复“唤醒”刚上手 PSSE 自定义模型的人常以为这件事就是“写个 Fortran 子程序挂上去”这么简单——直到仿真跑起来状态量一直不动或者一拍就发散。卡住人的往往不是方程写错而是没搞懂底层那条线主框架从不“调用你的模型”它调用的是一串行为你的子程序只是这些行为共同的入口。这就解释了为什么模型子程序里第一件事永远是判断标志位。主框架在 COMMON 内存里放着一组标志量每一拍按行为逐个改写它们再通过连接子程序回调到你的代码里。同一个子程序这一拍进的是初始化分支下一拍进的是微分方程分支。不理解这套调度方式模型写对了也调不通。2. 行为串联MODE、KPAUSE、MSTATE、MIDTRM 如何调度一个仿真步2.1 主框架眼里没有“模型”只有“行为”动态仿真在 PSS/E 主体框架中被切成一个个行为这些行为首尾相接串成一次计算。行为之间不是靠函数调用栈嵌套的而是靠“写标志量 → 回调 → 读标志量 → 再写标志量”这条链条推进。主框架在一拍里经历的行为序列初始化、网络方程求解、状态量积分、输出刷新等是固定的变化的是当前走到哪一步。模型子程序正是靠这个机制被反复进入。你写的SUBROUTINE不是一次性的构造函数而是像一张按标志位分岔的路由表第一个分支对应初始化第二个分支对应状态量更新第三个分支对应输出。每进入一次只执行其中一个分支其余全部RETURN掉。这种设计的直接后果是模型子程序必须是“无记忆”的。你不能在子程序里用SAVE保存跨拍状态也不能假设两次进入之间局部变量还活着。所有需要跨拍保留的量必须写进主框架分配的 STATE 或 VAR 区。很多新手在第一次写模型时把积分中间量放在局部变量里结果每一拍都被清零波形永远是一条直线——这就是没接受“反复唤醒”这个前提。另一个后果是调试时你不能靠“函数被调用了几次”来判断逻辑而要打印标志位。行为序列跑了多少轮你的子程序就被进入多少次但只有标志位匹配的分支才真正干活。2.2 四个标志位在 COMMON 内存里的语义与生命周期MODE、KPAUSE、MSTATE、MIDTRM 这几个量常驻 PSS/E 的 COMMON 内存模型子程序INCLUDE相应的公共块后就能直接读到。它们的生命周期比任何一次子程序调用都长是主框架与模型之间唯一的“会话上下文”。从命名和用法上看可以这样理解确切取值表以本地版本手册为准不同版本有细微差别MODE是最主要的行为开关决定当前这一拍该走哪个分支。模型里几乎所有的IF分支都挂在它上面。原文举的例子就是STRT行为与MODE1内联初始化语句因此被执行。KPAUSE与暂停、步进相关交互式调试和分步推进时会变化。它不改变模型方程但会改变行为推进的节奏。MSTATE作用于状态量层面涉及初始化与复位这类动作。MIDTRM与中间量、迭代过程相关在网络方程与状态量更新交替进行的阶段起作用。真正的要点不是背准每个名字而是理解它们的共同用法模型子程序在开头读取它们然后决定执行哪一段。标志位由主框架写、模型只读这条单向约定一旦被打破比如模型自己去改 MODE行为序列立刻错乱而且错得非常隐蔽——仿真不会报错只是结果慢慢偏掉。提示想知道某个标志位在某一拍到底取什么值最省事的办法是在模型里无条件打印一行MODE, KPAUSE, MSTATE, MIDTRM跑一个最短的仿真时长把日志对着行为序列看一遍。这比翻手册快得多。2.3 TBLCNC、CONEC、CONET 的分工谁负责自动挂载谁要你手写连接子程序是主框架与模型之间的“接线员”但它们不是同一类东西分工差别很大。TBLCNC负责把与机组相关的自定义模型联入 PSS/E 主模块它对用户不可见。模型按规则写完之后编译链接通过地址由 PSS/E 自动分配调用语句也由它自动生成。这意味着你不需要在某个连接子程序里手写CALL MYGOV(X, Y)也不用手工指定数组偏移——这两件事最容易出错的部分被主框架接管了。CONEC走的是另一条路径。与机组无关的模型比如挂在母线上、支路上的自定义模型需要用户在连接子程序里显式写入调用语句并指定地址。原文明确说了机组相关模型“不是通过子程序 CONEC 调用的”反过来读就是CONEC 这条路径上接线是你的责任。CONET这一类则更多承担数据面的搬运把 CON 区、STATE 区里的值按位置取出来供给计算。分工带来的实际影响是机组类模型的调试重点在“切片偏移算没算对”非机组类模型的调试重点在“连接子程序里那行 CALL 写没写、参数顺序对不对”。两类问题表现完全不同混在一起排查会浪费大量时间。3. 两参数签名SUBROUTINE modelname(MC, ISLOT) 为什么不多传一个参数3.1 MC机组内部码为什么不能直接拿母线号当索引与机组相关的自定义模型子程序形式被固定成SUBROUTINE MYGOV(MC, ISLOT)只有两个参数。MC是该模型所连机组的 PSS/E 内部机组码。为什么不直接传母线号因为母线号是用户侧标识会随拓扑调整、母线分裂、等值简化而变化而内部机组码是 DYRE 在建立动态数据时一次性分配的稳定索引整个仿真过程中不变。主框架内部所有按机组组织的数组都拿它当下标。这个选择很实用模型里如果要读取本机组相关的其他量比如机端电压、转子角用MC去索引内部数组就行不需要再做一次“母线号 → 机组”的查表。代价是MC对用户没有语义光看数字你不知道是哪台机调试时得配合输出文件对照。ISLOT是另一回事它是该模型在内部数组位置表里的入口指针。注意这里“指针”是 Fortran 77 意义上的整数下标不是现代语言的引用。3.2 ISLOT 与 STRTIN 四行入口NC/NS/NV/NI 的切片规则DYRE 活动结束时每个与机组相关的模型都会在位置表里分到一块连续区域。你在模型里通过STRTIN这个二维数组拿到自己那块区域的首地址入口对应声明含义数量为 0 时STRTIN(1, ISLOT)NC本模型第一个 CON 的位置返回 0STRTIN(2, ISLOT)NS本模型第一个 STATE 的位置返回 0STRTIN(3, ISLOT)NV本模型第一个 VAR 的位置返回 0STRTIN(4, ISLOT)NI本模型第一个 ICON 的位置返回 0这套机制的本质是把全局的CON(1)、STATE(1)、VAR(1)、ICON(1)四个大数组按模型切片。主框架只维护一张“谁从哪开始、占多少个”的表模型只拿到基址靠约定的偏移访问自己的那一段。SUBROUTINE MYGOV(MC, ISLOT) C C MC : 本模型所连机组的 PSS/E 内部机组码 C ISLOT : 本模型在数组位置表中的入口 C INCLUDE COMON4.INS C INTEGER MC, ISLOT INTEGER STRTIN REAL CON, STATE, VAR COMMON /CONS/ CON(1) COMMON /STATE/ STATE(1) COMMON /VARS/ VAR(1) COMMON /ICON/ ICON(1) COMMON /STRTIN/ STRTIN(4, 1) C INTEGER JC, JS, JV, JI C C ---- 先把本模型的四个切片基址取出来 C 任何一个数量声明为 0对应基址就是 0不能当合法下标用 JC STRTIN(1, ISLOT) JS STRTIN(2, ISLOT) JV STRTIN(3, ISLOT) JI STRTIN(4, ISLOT) C IF (MODE .EQ. 1) THEN C ---- 初始化行为用 CON 区里声明的常数给 STATE 赋初值 IF (JS .NE. 0 .AND. JC .NE. 0) STATE(JS) CON(JC) RETURN END IF C C ---- 其余 MODE 分支按对应行为的语义补齐 C 每个分支只做一件事做完立刻 RETURN RETURN END代码里最关键的两行不是赋值而是JC、JS这几个基址的获取以及后面的NE. 0判断。STRTIN返回的 0 是“本模型没有这一类量”的信号不是“第 0 号元素”。直接拿 0 去索引读到的会是别人的数据或数组越界而这种错误在仿真里通常不报错只表现为数值莫名其妙。参数含义再明确一遍ISLOT的取值由主框架分配你的模型里不负责计算它STRTIN的第二个维度就是ISLOT的上界偏移量必须严格落在[0, 声明数量-1]区间内越界就是踩到邻居模型的数据。3.3 MODE1 初始化段的写法指针为 0 时该做什么初始化段是整个模型里最容易写错的地方因为它同时承担三件事校验切片、给状态量赋初值、把常量区的值搬进变量区。三件事混在一起出问题就很难定位。稳妥的写法是分三步走。第一步只做基址获取和零值判断把不可用的切片挡在外面第二步按模型声明的数量逐个赋初值不在这一步做任何计算第三步做一致性检查比如时间常数是否为正、增益是否落在合理区间——这一步不做也不会崩但能在数据卡片填错时提前暴露。零值判断尤其不能省。一个模型可能只声明了 CON 和 STATE没声明 VAR 和 ICON那么STRTIN(3, ISLOT)和STRTIN(4, ISLOT)都是 0。如果初始化段里有一段代码要往 VAR 写东西就必须先判断JV .NE. 0。很多人是在仿真跑了一半发现输出通道全是零回头才发现是这里少了一个判断。还有一个隐性陷阱初始化段只应该写 STATE 和 VAR不应该读它们。DYRE 阶段刚分配完STATE 里可能是上一轮残留或者未定义值拿它参与计算等于把不确定性带进初值。4. USRMDL 记录与调用顺序从一张数据卡片到主循环挂载4.1 记录逐字段一次对主框架的“资源申报”自定义模型在动态数据输入文件里以USRMDL记录的形式出现格式为BUSID USRMDL IM modelname IC IT NI NC NS NV data list一行实际的数据卡片大概长这样字段间分隔按本地版本的卡片格式走101 USRMDL 1 MYGOV 0 0 0 4 1 0 1.0 0.05 0.2 5.0 0.0各字段的含义BUSID标识模型挂在哪台机组上是模型与物理设备的绑定关系。USRMDL是固定关键字DYRE 靠它识别这是一条自定义模型记录。IM是模型序号同一台机组上挂多个同类模型时用来区分。modelname是模型名必须和 Fortran 子程序名对应长度受限制。IC、IT、NI、NC、NS、NV是六个尺寸声明本质是模型对主框架的一次资源申报我要几个整数、几个常量、几个状态、几个变量。DYRE 据此在位置表里划出对应大小的切片把首地址写进STRTIN。data list是按声明顺序排下来的初值常量和状态初值都从这里进。NI、NC、NS、NV与STRTIN四个入口一一对应这条映射是整个机制的骨架。前两个字段IC、IT在不同版本手册里的表述有细微差别实际写模型时以本地版本手册为准别跨版本抄卡片——字段错位一位后面所有初值都会跟着错而且 DYRE 不一定报错。注意data list的个数必须等于声明的总量多一个少一个都可能让解析器把下一行当成数据。手写卡片时建议先数一遍再提交。4.2 调用顺序的因果链发电机→补偿器→稳定器→励磁→调速器→混合模型库里的调用顺序是固定的发电机模型电流补偿器模型励磁稳定器模型励磁模型汽轮机/调速器模型混合模型这个顺序不是按字母排的也不是按重要性排的而是按信号流向排的。发电机模型先算出机端电气量补偿器在此基础上修正电流注入稳定器需要发电机的转速和功率偏差作为输入所以排在发电机之后励磁模型接收稳定器给出的附加电压信号因此又排在稳定器之后调速器处理的是机械侧的转速与功率与电气侧相对独立但在同一拍内要和电气量对齐。顺序错了会怎样在同一拍里后调用的模型读到的是前一个模型刚算出来的值。把励磁提到稳定器前面励磁这一拍拿到的就是上一拍的稳定器输出相当于凭空多了一拍延迟。对慢过程影响不大对快速励磁系统就可能在阶跃响应里看出相位偏移严重时出现本不该有的振荡。这也解释了一条实践中的规矩不要为了让自己的模型“拿到最新值”去调整调用顺序。真正需要跨模型协同的场合应该用标志位分阶段而不是靠顺序碰运气。4.3 TBLCNC 不可见编译即挂载背后的取舍机组相关模型由 TBLCNC 自动联入主模块、自动分配地址用户不写CALL、不指定偏移。这条设计省掉了大量手工接线也把最容易出错的部分锁进了主框架。好处很直接模型与主框架之间的接口由签名和位置表约定死了用户只要按SUBROUTINE modelname(MC, ISLOT)写、按声明填卡片就能挂上去。换个模型不需要动连接代码也不存在“忘了在连接子程序里加一行 CALL”这种事故。代价是黑匣子感很强。挂载过程不可见地址分配不可见一旦出问题你只能从现象反推。比如模型明明编译进去了但完全没生效可能是名字没对上、可能是卡片位置错了、可能是链接到了旧库——三种原因在外部表现上几乎一样。这套机制把复杂度从“每次都要写”变成了“出问题时要猜”是一笔明确的取舍。理解了这一点排查思路就清楚了先确认挂载成功模型被进入过没有再确认切片正确基址落在哪最后才看方程本身。5. 排查清单自定义模型接入时最容易翻车的五处5.1 状态量恒为 0把指针 0 当成了合法下标现象模型编译通过DYRE 也顺利读完卡片仿真能跑完不报错但所有输出通道画出来都是水平线状态量从头到尾是初值。原因STRTIN某一行返回 0表示本模型没有声明这一类量。代码里没做判断直接拿 0 去索引VAR(0)、ICON(0)读到的要么是数组头部的哨兵值要么是邻居模型的数据。写进去的值也一样落到了别人的地盘自己的状态量自然永远不变。解决在基址获取之后立刻做一次零值判断把所有可能为 0 的入口挡在计算之外。更稳妥的做法是在初始化段里加一段自检把四个基址和声明的数量打印出来人工确认切片区间互不重叠、落在合理范围内。这条自检代码在模型稳定之后可以留着代价只是一行输出。5.2 模型名写了却找不到长度与大小写的隐形约束现象动态数据文件里明明写了MYGOVERNOR编译也过了仿真时这个模型就像不存在一样没有任何报错。原因模型名有长度限制超过限制会被截断截断后的名字和 Fortran 子程序名对不上挂载就失败了。另一种情况是大小写处理不一致卡片里写小写、编译出来的符号是大写某些平台上链接阶段就匹配不到。这类问题不产生错误码只产生“静默忽略”。解决模型名控制在限制长度以内统一用大写卡片和源码里的名字逐字符对齐。改完名字之后重新跑一次最短仿真确认模型确实被进入过——最直接的验证是在子程序入口无条件打印一行看日志里有没有。5.3 改了代码结果一点没变链接到了旧产物现象改了一行增益重新编译仿真结果和改之前一模一样反复几次都一样开始怀疑是不是改错了文件。原因编译产物和运行环境里的库不是同一份。常见情况是编译输出到了 A 目录运行时加载的是 B 目录下的旧版本或者上一次编译失败但没注意实际跑的还是更早的成功产物。这类问题在多人协作的环境里出现频率极高。解决编译后核对产物时间戳确认运行环境指向的是刚生成的那一份。养成“改完先故意改错一个显眼的值确认结果变了再改回来”的习惯——这一步花不了一分钟能省掉半天自我怀疑。5.4 调整调用顺序后结果变了一拍延迟的等价效应现象为了让自己写的模型拿到“更新鲜”的输入把某个模型的调用位置往前挪了一位结果稳态没变暂态波形却出现了相位偏移阶跃响应的第一个峰值明显变小。原因同一拍内后调用的模型读到的是前面模型刚算出的值。顺序一变某个环节的输入就从“本拍值”退化成“上一拍值”等价于在信号链里插入了一拍延迟。对慢速环节无所谓对快速环节就是实打实的动态特性改变。解决不要为了拿最新值去改顺序。需要跨模型先后协同的场合用标志位分阶段处理在同一个 MODE 分支里按固定次序完成多个模型的更新而不是靠整体顺序碰运气。5.5 别人的模型突然行为异常越界写坏了邻居切片现象自己新加的模型上线后原本运行正常的另一个模型输出开始漂移幅度不大但持续发散。把新模型摘掉一切恢复正常。原因新模型在写 STATE 或 VAR 时偏移算多了写到了相邻模型的切片里。位置表是连续分配的你的切片末尾紧接着就是下一个模型的起始地址多写一个元素就污染了别人。这种污染不会触发任何检查只会让邻居的初值莫名其妙地变。解决严格用STRTIN给的基址加声明范围内的偏移绝不用硬编码的绝对下标。模型里所有循环的上界都写成声明数量不写魔法数字。排查时把每个模型的切片区间打印出来检查有没有重叠——重叠就是越界的直接证据。6. 设计取舍为什么这套“标志位 指针表”的老结构还值得学把这几层拆开看整套机制其实只有三个约定全局标志位决定行为分支两参数签名固定调用接口位置表切片分配存储。用今天的眼光看这三条都是反模式全局状态难测、隐式契约多、越界无法自动检测、并发场景下完全不成立。但如果把时间放回它被设计的年代取舍就很清楚。Fortran 77 没有动态内存分配、没有结构体、没有模块命名空间一个能容纳成千上万个模型、跨厂商扩展、还要保证 ABI 稳定的框架几乎只能这么写。一个大数组加手工切片是当时唯一既省内存又保证访问速度的方案用标志位而不是参数传上下文是因为每增加一种行为就要改一次所有模型的签名代价高到无法接受让 TBLCNC 自动挂载是因为手工接线一旦出错排查成本远超框架本身。它的真正优势在于增量扩展新增一个模型不需要改主框架的任何一行代码只要按签名写、按声明填卡片、编译通过就能挂上去。这个性质在几十年的生命周期里被反复验证过——模型库能长到今天这个规模靠的就是这条。代价同样明确所有错误都推迟到运行期才暴露而且往往以数值异常而非报错的形式出现。这就决定了验证方法必须是“先小后大”。我一般会先写一个只声明一个 CON、一个 STATE 的最小模型初始化段里把STRTIN四个基址和ISLOT全部打印出来跑零秒仿真确认切片落在预期区间、不与其他模型重叠然后再往里填第一条微分方程用阶跃输入看响应方向和量级对不对最后才接真实数据。三步走下来绝大部分接线问题在第一步就暴露了。从那以后我每接一个新模型都强制先走一遍这个最小验证流程哪怕自认为对机制已经很熟。希望帮到你。本文还有配套的精品资源点击获取
