企业官网建设流程全解析

同步发电机突然短路这段暂态过程，是所有电机学课程里最刺激、也最容易被公式劝退的内容。说刺激，是因为它把“磁链不能突变”“电流不能突变”这些抽象概念一次性甩到你面前——短路瞬间定子电流能冲到额定电流的十倍以上，绕组端部要承受巨大的电磁力，转轴要承受振荡扭矩；说劝退，是因为整个过程同时包含超瞬态、瞬态、稳态三种分量，衰减快慢由不同的时间常数控制，手推公式能推满一黑板。这篇内容想做的事情很具体：用MATLAB把同步发电机突然短路的暂态过程仿真出来，一步步看清短路后电流为什么先冲高、再快速回落、再慢慢稳定，并告诉你怎么搭模型、怎么填参数、怎么看波形、怎么排查问题。内容既适合正在做电机学课程设计或电力系统暂态作业的同学，也适合想验证自己对暂态过程理解的工程师。

1. 突然短路这个课题：为什么课程设计和工程里都绕不开

1.1 短路到底有多“猛”：从绕组到轴系的连锁反应

同步发电机正常运行时机端电压稳定，电流在额定范围内，电磁过程处于一种平衡态。突然三相短路等于把机端三相电压瞬间压到零，这个变化比正常运行时的任何扰动都剧烈。结果是什么呢？

定子电流瞬间出现很大的冲击分量，峰值可以达到额定电流的十几倍。电流在线圈里会产生与电流平方成正比的电磁力，端部绕组受到的机械应力最大，历史上不少发电机端部损坏事故就发生在出口短路之后。与此同时，电磁转矩也会发生强烈振荡，这个转矩通过联轴器传到汽轮机或水轮机的轴系，相当于转轴被反复“猛拧”，轴系扭振累积下来可能导致疲劳损伤。所以突然短路不只是电机学的一个理论问题，它直接关联到发电机设计时的机械强度校核、断路器开断能力选择、继电保护整定计算。这也是为什么每个电气工程专业的学生都要碰一次这个课题。

我做这个课题的时期，最直观的感受是：书本上那些“超瞬态电抗”“暂态电抗”“同步电抗”的排序，只有在看到仿真波形那一刻才真正变成有物理意义的东西。短路电流先按超瞬态电抗决定的数量级冲上去，再按暂态时间常数回落，最后停在由同步电抗决定的稳态值上——整个过程在波形上一目了然。

1.2 用MATLAB仿真能解决什么问题

手推解析公式当然能算出短路电流的表达式，但算出来的结果要和动态过程对应上，难度不小。MATLAB仿真的价值在于把黑板上的公式变成可观察、可交互的曲线：

• 可以看到三相电流各自不同的直流偏置，理解为什么某相第一峰值特别大；
• 可以观察励磁电流在短路瞬间的冲击，以及阻尼绕组电流的衰减；
• 可以同时看转速和电磁转矩的变化，感受机电暂态和电磁暂态的时间尺度差异；
• 可以改参数，比如把阻尼绕组去掉、把短路时刻往后挪，波形立刻变样，这种“参数敏感性”是纯解析很难直观感受到的。

对课程设计来说，MATLAB配合Simulink能快速出图、快速交差；对想深入理解的人来说，自己写状态方程再求解，能逼你把每个方程、每个变量都搞清楚。两条路线我都走过，这篇内容会把两种做法都讲清楚，并给出我踩过的坑。

2. 暂态过程的底层物理：记住磁链守恒就够了

2.1 磁链守恒是总开关

突然短路分析的核心，概括起来就一句话：绕组磁链不能突变。这是电磁惯性定律在电路里的体现。物理上，一个无源绕组如果外部电压突然变化，为了维持自身磁链不变，必然感应出电流来“顶住”这个变化。

用生活里的场景类比：你坐在行驶的公交车里，司机突然急刹车，你的身体会向前冲，因为身体想保持原来的速度。电流也类似，短路前定子绕组里有一个稳态磁链，短路瞬间外部条件变了，但磁链这个“速度”不能跳变，于是绕组里必须产生新的电流分量来维持磁链连续。

同步发电机的转子上有励磁绕组，还可能有阻尼绕组。短路瞬间，定子电流要维持定子磁链，同时励磁绕组和阻尼绕组也要维持它们自己的磁链。这三个绕组之间通过互感耦合在一起，结果就是短路后的电流波形由多个衰减分量的叠加构成。磁链守恒是分析突然短路的总开关，所有分量都能从这个原理里推出来。

2.2 定子电流的三段式衰减从哪来

同步发电机突然三相短路后，定子电流里会出现几个性质完全不同的分量：

• 基波交流分量的起始部分，由超瞬态电抗 X''d 决定。短路刚开始，阻尼绕组和励磁绕组都努力维持磁链，相当于把电枢反应的磁通“挡”在漏磁路径上，电机呈现的阻抗很小，所以电流初始幅值非常大。
• 基波交流分量随着阻尼绕组电流衰减，过渡到由暂态电抗 X'd 决定的水平。阻尼绕组电阻小、时间常数短，大约几十毫秒内它的电流就衰减完了，这时只剩励磁绕组在维持磁链，短路电流回落到一个中间水平。
• 基波交流分量随着励磁电流衰减，最终进入由同步电抗 Xd 决定的稳态短路电流。励磁绕组时间常数长，可能要一两秒，所以这个阶段的衰减过程明显更“拖沓”。

除了这些交流分量，定子电流里还有一个衰减的直流分量。这个直流分量是为了维持短路瞬间定子磁链不突变而产生的，它的大小和短路发生的时刻直接相关：某相绕组磁链在短路瞬间越大，该相直流分量就越大。直流分量叠加在交流分量上，会让波形整体偏移，于是不同相的电流波形看起来“胖瘦”不同、第一峰值出现的时间也不同，这就是所谓的非对称短路电流。断路器开断能力校核时，必须按非对称峰值来考虑，而不是只看对称交流分量。

2.3 电抗和时间常数的排列规律

同步电机的几个电抗，按大小排序是 Xd > X'd > X''d。这个排序本身就是暂态过程的“目录”。为什么同步电抗最大？因为磁通走气隙主磁路，电枢反应全部参与。为什么暂态电抗小？因为励磁绕组被磁链守恒“短路”了，电枢反应磁通被逼到只走励磁绕组的漏磁路径，等效磁路变短，电抗变小。为什么超瞬态电抗更小？因为阻尼绕组也加入进来，磁通被进一步逼到更小的漏磁路径上。

与电抗配套的是时间常数，它们决定了各分量衰减的快慢：

看时间常数就知道，超瞬态过程快得像“啪”一下，暂态过程慢一些，稳态过程则要等很久。仿真时如果想完整看到这三段，仿真时长至少设到2~3秒，否则只能看到超瞬态和暂态的一截。

3. 仿真路线选型：Simscape模块还是自建状态方程

3.1 路线A：Simscape Electrical里的同步电机模块

Simscape Electrical（早期叫SimPowerSystems，再早叫Power System Blockset）里提供了现成的同步发电机模型，模块名一般是 Synchronous Machine SI Fundamental 或 Synchronous Machine pu Fundamental。这个模块内部已经实现了派克变换、转子绕组方程、运动方程，外部只留出三相端子、励磁电压输入 Vf、机械功率输入 Pm 和转速/角度输出端口。

用这条路线，你不需要自己推导和编写微分方程，核心工作在于：

• 正确填写电机参数（额定容量、电压、频率、电抗、时间常数、惯性常数）；
• 配置故障模块和测量模块；
• 设置初始条件让仿真起点处在空载额定电压状态；
• 从Scope里读取波形并解释。

如果目标是课程设计快速出图、验证教材结论、或者把发电机接入一个更大的电力系统做系统级分析，这条路线最合适。模块里的参数是按厂商数据手册习惯组织的，比如直接填 Xd、X'd、X''d 和 T'd0、T''d0，不需要自己去换算电感矩阵，这点很友好。

3.2 路线B：dq0标幺方程自己编程积分

如果你想真正钻进电机内部，Simscape模块多少像个黑箱。它内部怎么处理磁链耦合、怎么处理饱和、怎么处理坐标系，你只能通过外部特性去猜。这时可以自己写dq0坐标下的标幺方程，用MATLAB脚本配合 ode 求解器积分，把每个状态量的变化过程全部掌握在自己手里。

自建方程的核心工作包括：

• 列写磁链方程和电压方程；
• 求短路前的稳态初值；
• 在0到短路时刻先积分到稳态，然后切换故障条件继续积分；
• 把结果从dq0坐标系变换回abc三相坐标绘图。

这条路线比拖模块费时间，但做完之后你对同步电机的理解会上一个台阶。尤其是“初始稳态怎么求”“故障瞬间为什么状态变量不能跳变”这些问题，是一边写代码一边真正想明白的。

3.3 两条路线的取舍建议

我的建议很实际：如果时间紧、重点是看结果，走路线A；如果老师要求写推导过程，或者你想在文章里写“本文基于dq坐标建立同步发电机模型并编程实现”，那就走路线B。两条路线不矛盾，先用路线A验证参数，再用路线B加深理解，是很高效的学习路径。

4. Simscape实操：搭一个三相突然短路仿真的详细步骤

4.1 模型骨架与模块库定位

打开Simulink，新建一个模型，按下面的模块列表搭建骨架：

• powergui：路径在 Simscape > Electrical > Specialized Power Systems > powergui。这个模块必须放在模型里，Simscape电力系统才能求解，它相当于整个电气网络的“求解环境”。
• Synchronous Machine SI Fundamental：路径在 Simscape > Electrical > Specialized Power Systems > Fundamental Blocks > Machines。
• Three-Phase Fault：路径在 Simscape > Electrical > Specialized Power Systems > Fundamental Blocks > Elements。
• Three-Phase V-I Measurement：路径在 Simscape > Electrical > Specialized Power Systems > Fundamental Blocks > Measurements。
• Constant模块：两个，一个给 Pm 提供机械功率输入，一个给 Vf 提供励磁电压输入。
• Scope：观察电流、电压、转速、励磁电流等信号。

连接关系不复杂：发电机三相输出端 A、B、C 接三相连到 Three-Phase V-I Measurement 的一组输入，测量模块的输出再接 Three-Phase Fault 的输入端，故障模块的输出端接地；测量模块的电压、电流信号引到Scope。Pm 输入给机械功率，Vf 输入给励磁电压，转速端口 wm 也可以引到Scope。模型搭完就是这个拓扑，关键是后面的参数。

4.2 发电机参数怎么填

双击同步电机模块，参数页面按机组实际数据填写。以一台典型教学用隐极发电机为例，参数大致如下：

这里有个容易踩的坑：模块参数页里，如果选的是“Specify by manufacturer data”，填的是上面的电抗X和时间常数T'do系列；如果选“Specify by d/q inductances”，填的是Ld、L'd、L''d这些电感值。不要混着填，选一种方式到底。填完电抗后，模块内部会自动换算成磁链方程系数，前提是额定频率填对——因为有名值和标幺值之间的换算依赖频率，频率错了整个波形幅值都会跟着错。

4.3 故障、测量与初始条件配置

Three-Phase Fault 模块的设置，核心就几个点：

• Fault type选 Three-Phase Short Circuit，这是最标准的三相同时短路；
• Transition times填短路发生时刻，比如 [0.2 100]，意思是0.2秒时故障开始，100秒时故障结束——对观察完整暂态衰减过程来说，可以认为故障一直存在；
• Fault resistance填0.001左右。填0当然最理想，但纯0电阻在某些求解器里容易引起数值振荡，给一个极小的电阻更稳；
• Initial status设 open，让故障模块在短路前是断开的，发电机空载正常运行。

初始条件这一步很关键，处理不好模型会报“Initialization”错误。推荐做法：先把机械功率 Pm 设为0（空载），励磁电压 Vf 设一个额定空载励磁值，然后用 powergui 里的 Machine Initialization 工具自动完成稳态初始化。具体操作是双击 powergui，打开 Tools 里的 Machine Initialization，选中同步电机模块，软件会根据你填的 Vf、Pm 和发电机额定值自动算出一套满足稳态条件的初始状态。点 Apply 后，模块内部的初始转子角度、初始励磁电流、初始功角都会写进去，仿真起点就是发电机空载带额定电压的稳态。

如果你不想用初始化工具，也可以手动给模块填初始状态，但很容易出现“初始条件跟电路方程不匹配”的问题。我建议一律用 Machine Initialization，省时省事。

4.4 典型运行结果长什么样

搭好模型、参数填对之后，仿真时间设3秒，运行。波形大致会这样发展：

• 0到0.2秒：发电机空载，定子三相电流为0（严格说是很小的空载励磁电流对应的空载电流，对同步电机空载而言定子电流理论上为零），机端电压为标准正弦波；
• 0.2秒短路瞬间：A、B、C三相电流同时突然出现很大的冲击。由于三相电流的直流偏置不同，三根曲线不是简单的等比放大关系，而是各自带着不同的偏移量；某相如果在短路时刻磁链最大，那一相会出现特别突出的第一峰值；
• 短路后最初0.1秒：电流下降很快，这是超瞬态阶段，由阻尼绕组过程主导；
• 0.2秒到1秒左右：电流继续衰减，速度明显变慢，这是暂态阶段，励磁绕组过程主导；
• 1秒以后：电流逐渐逼近稳态短路电流，数值大约等于空载电动势除以同步电抗X d，对额定电压空载短路来说，即1/1.8≈0.56 p.u.。

同时看励磁电流 If，可以发现短路瞬间励磁电流也会发生一个冲击，然后慢慢回稳。这是因为定子突然短路产生的电枢反应磁通试图穿过励磁绕组，励磁绕组为了维持磁链，电流必然要调整。看电磁转矩Te，能看到强烈的振荡，甚至出现负值——发电机瞬时变成了拖着制动转矩的“电磁刹车”，转速会有一个小幅波动。

5. 自建状态方程的仿真实现：从公式到曲线

5.1 先把电机方程整理成能编程的形式

如果你选择自己写代码，第一步是把同步电机的dq0方程整理清楚。这里用标幺值、发电机惯例，定子电流方向按输出为正。

磁链方程可以写成：

• 定子d轴磁链：ψd = Xad(id + if + iD) + Xl·id
• 定子q轴磁链：ψq = Xaq(iq + iQ) + Xl·iq
• 励磁绕组磁链：ψf = Xf·if + Xad(id + if + iD)
• d轴阻尼绕组磁链：ψD = XD·iD + Xad(id + if + iD)
• q轴阻尼绕组磁链：ψQ = XQ·iQ + Xaq(iq + iQ)

其中Xad、Xaq是电枢反应电抗，Xl是定子漏抗，Xf、XD、XQ是各转子绕组的总自感电抗（含各自漏感）。

电压方程：

• 定子d轴：vd = -rs·id + ω·ψq - dψd/dt
• 定子q轴：vq = -rs·iq - ω·ψd - dψq/dt
• 励磁绕组：vf = rf·if + dψf/dt
• d轴阻尼绕组：0 = rD·iD + dψD/dt
• q轴阻尼绕组：0 = rQ·iQ + dψQ/dt

运动方程：

• 2H·dω/dt = Tm - Te，其中Te = ψd·iq - ψq·id

三相突然短路时，故障点的电压被钳到零，所以对发电机端而言 vd=0、vq=0（当然，这里忽略了故障点到机端之间的线路阻抗，如果做机端短路就是直接为零）。励磁电压 vf 保持不变，机械功率 Tm 在电磁暂态这个时间尺度内也可以当作不变，除非你特意接调速器模型。

这组方程的特点是：磁链方程是代数约束，电压方程是微分方程，两者耦合在一起。编程时通常把磁链作为中间变量，从当前电流状态反解出磁链，再结合电压方程求出电流导数，最后积分。d轴和q轴之间因为有 ω·ψq、ω·ψd 这样的交叉项，不能拆开独立求解。

5.2 两步法：先求稳态初值，再切换故障继续积分

写脚本仿真突然短路，最严谨的思路是两步法。

第一步，先求短路前的稳态。发电机空载运行，定子电流为零，机端电压为额定值，励磁电流 if0 由空载励磁决定。这个初值可以解析求出来：空载时 id0=iq0=0，iD0=iQ0=0，ψd0=Xad·if0，ψq0=0，机端电压全部由励磁磁链产生，于是 vf0 和 if0 的关系由励磁回路稳态方程决定。把这个稳态作为积分起点，先让仿真从 t=0 运行到 t=tshort，验证机端电压确实是额定正弦波，再准备切换故障。

第二步，在 t=tshort 时刻，把定子电压约束从“空载额定电压”换成“vd=vq=0”，也就是机端短接。记住一个关键点：切换时刻所有状态变量（各绕组电流、磁链、转速）必须保持连续，不能人为跳变。因为电感里的电流和磁链在物理上是不能突变的，如果你在代码里硬给某个电流赋了新值，波形就会出现不真实的尖峰。实现时可以用一个事件检测函数，让积分在 t=tshort 处停下，把当前状态读出来作为故障后的初值，再调用一次积分器继续算故障后的动态。这就是MATLAB里 ode 系列求解器和事件（Event）功能的典型用法。

5.3 一个可落地的MATLAB脚本骨架

下面的代码骨架展示了核心思路，你可以把它扩展到完整实现。这里省略了电感矩阵的具体数值，因为不同机组参数不一样，你需要按自己的数据填进去。

% 同步发电机突然三相短路仿真骨架 % 状态向量 x = [id; iq; i0; if; iD; iQ; omega; delta] % 先用解析公式求空载稳态初值 x0 % 阶段一：0~tshort 常态运行积分，检查空载电压 [t1, x1] = ode15s(@(t,x) gen_dynamics(t,x,0), [0 tshort], x0, ... odeset('Events', @(t,x) faultEvent(t,x,tshort))); % 阶段二：tshort 起短路，积分到 tfinal [t2, x2] = ode15s(@(t,x) gen_dynamics(t,x,1), [tshort tfinal], x1(end,:)); % 合并结果并变换到abc坐标绘图

核心的动态函数大致是这个结构：

function dx = gen_dynamics(t, x, fault_flag) % 从状态向量解出dq0电流和转子电流 % 先根据磁链方程反解电流：i = L \ psi % 再根据电压方程算导数 if fault_flag == 0 vd = 0; vq = Et; % 空载情况，Et为额定电压对应的内部量 else vd = 0; vq = 0; % 机端三相短路 end % 带入定子、励磁、阻尼绕组的电压方程求di/dt % 带入运动方程求domega/dt和ddelta/dt end

很多人会在这里卡住：为什么空载时 vq 不是零？因为你在dq旋转坐标系里看定子量，空载电压是一个旋转向量，在d、q轴上各有分量，具体怎么分配取决于转子位置和初始功角。所以你第一步求稳态初值时，除了要算出 if0，还要把初始转子位置角确定下来，否则后面的Park反变换就画不出正确的abc三相波形。这也是为什么用现成模块省心——模块把初始角度都替你管好了；自己编程则必须每一步都清楚。

5.4 数值求解要注意的刚性问题

同步电机方程是个典型刚性系统：阻尼绕组时间常数T''d是几十毫秒量级，机械转动的时间尺度是秒级，两者相差几十倍。用 ode45 这种非刚性求解器，要么算得很慢，要么干脆发散。我实测下来，ode15s 和 ode23t 这类刚性求解器最合适。建议设置 RelTol=1e-6，MaxStep不要太大，比如0.005秒，否则短路瞬间的陡变会被大步长糊弄过去，峰值都看不准。

还有一个细节：短路瞬间电流从零跳到十几倍额定值，虽然电流本身连续，但它的变化率极大。积分器在事件触发切换后，第一步计算很关键，最好在切换后的初始段强制用小步长。可以在第二次积分时给 odeset 里的 InitialStep 设置一个小值，比如 1e-5，让求解器先“看到”陡变过程，后面再自动放大步长。

6. 结果分析实战：从波形反推电机参数

6.1 三段式波形这么读

拿到仿真曲线，第一件事不是看细节，而是先找三段式结构。我用一张表总结怎么对应理论：

举个例子，如果 X''d=0.2 p.u.，X'd=0.3 p.u.，Xd=1.8 p.u.，那么短路电流的周期分量理论上会经历“5 p.u. → 3.3 p.u. → 0.56 p.u.”的下台阶过程。如果你的波形里看不到三个台阶，基本可以断定参数设置或者初始条件出了问题。

这里有个容易混淆的点：短路电流的瞬时值峰值不是1/X''d直接等于的数值，因为瞬时值里包含衰减直流分量，而且第一峰值出现在短路后约半个周波处，这段时间里超瞬态分量已经开始衰减。所以用表里的1/X''d去核对时，要留出一些余量，不要指望完全相等。

6.2 如何分离周期分量和非周期分量

三相短路电流里，周期分量和非周期分量混在一起。如果你想分别观察它们，有个很经典的技巧：利用相隔半个周期的两个瞬时值。

假设某一相电流在 t 时刻瞬时值为 i(t)，在 t+T/2 时刻为 i(t+T/2)，其中T是工频周期。周期分量在这两个时刻大小相等、方向相反，而非周期分量衰减很慢，近似认为这两个时刻的直流分量相等。于是：

• 周期分量瞬时值近似为 [i(t) - i(t+T/2)] / 2
• 非周期分量近似为 [i(t) + i(t+T/2)] / 2

在MATLAB里实现就是找到每个采样点往后推半个周期的数据，做加减运算。把周期分量单独画出来，你会看到一条清晰地从“大”衰减到“小”的交流包络；把非周期分量单独画出来，可以看到它从某个初始值指数衰减到零。这个过程如果手工写代码，还能顺便练习数据切片和索引操作，比直接看混合波形要更能锻炼人。

做完分离之后，还可以做进阶操作：把周期分量包络取对数，画在半对数坐标上。正常情况下包络会呈现两段斜率不同的直线段，前段陡、对应T''d，后段缓、对应T'd。用MATLAB的 polyfit 对这两段分别做线性拟合，就能从仿真波形里“反推”出时间常数，甚至和你填在参数卡片里的 T'd0、T''d0 对比。这个操作在课程设计里属于加分项，在论文里就是完整的参数辨识过程。

6.3 几个值得单独观察的量

不要盯着定子电流看个没完，同步发电机突然短路时，其他量同样值得体检：

• 励磁电流：短路瞬间会有一个明显冲击，可能冲到空载值的2~3倍，然后逐渐回落。这个冲击就是励磁绕组维持磁链的直接表现。如果励磁电流没有冲击，说明你用的模型可能把励磁回路简化成了恒定电流源，而不是真实的有惯性绕组。
• 电磁转矩：短路后转矩会出现高频振荡，甚至瞬时变负，平均值也不为零。观察转矩可以理解轴系受到的冲击从哪里来。把转矩曲线做FFT，能看到接近两倍工频的振荡分量，那是定子直流分量与转子磁场相互作用的结果。
• 转速：转速变化相对缓慢，波动幅度取决于惯性常数H。H越小，转速波动越明显。如果仿真里转速掉得离谱，先检查机械功率和初始状态是否一致，不要急着怀疑方程。
• 机端电压：如果你把测量模块放在故障点之前，短路瞬间电压会断崖式掉到接近零。这段电压跌落过程也是判断故障时序是否正确的重要标志。

这些量并不是孤立的，它们通过磁链、转矩、转速互相耦合。看波形时养成多量对比的习惯，比单看电流更有收获。

7. 高频问题与调试经验记录

7.1 仿真中常见的几个坑

这些坑我基本都踩过，整理成速查表，按现象排查会快很多：

7.2 三条能显著提高成功率的总原则

第一，先跑空载稳态，再做故障。不要一上来就把故障时间设成0，让模型从初始状态直接短路。这样既容易初始化失败，也看不清短路前和短路后的对照。我习惯的做法是短路前至少留0.2秒，这段时间用来确认空载电压正常、波形是一条干净的正弦波。

第二，短路电阻不要给0。理论上的金属性短路是0欧，但数值求解时完全0电阻会让电流变化率趋近无穷大，求解器非常难受。给0.001Ω甚至0.0001Ω，物理上等同于金属性短路，数值上却稳定得多。这条经验在Simscape和自建方程里都适用。

第三，输出数据，别只依赖Scope。Scope看趋势可以，做分析不行。把电流、励磁电流、转矩用 To Workspace 模块或者脚本里直接存到工作区，再统一用脚本重绘，加图例、加标注、把特写区间放大。仿真结果能不能写进报告，往往取决于你重绘的功夫，而不是原始Scope截图。

7.3 一点个人体会

同步发电机突然短路这个题目，看起来有标准答案，但真正自己动手做一遍，会发现理论和仿真之间隔着一层“参数映射”“初值处理”“数值稳定性”的皮。我做这个课题的经历是：前三天都在跟初始化错误和发散问题搏斗，等到波形终于正常了，反而觉得前面那些坑是经验加速器。如果你也卡在某个报错上，记住一个原则——大多数问题不是电机方程错了，而是参数单位、初始条件、求解器设置这三件事没对齐。按本文的顺序排查，大概率能把波形救回来。