6.1 研究类型详解:给求解器下任务书 本节摘要:研究类型回答的是"你想问方程组什么问题":稳态问平衡、瞬态问过程、特征频率问天性、频域问响应。本节讲清四类研究的适用场景与容差约定,并用热时间常数给铜排补一个上电六千秒的瞬态仿真。 问不同的问题,下不同的任务书 前五章你一直用稳态研究,它回答"系统最终停在哪个平衡态"。但工程问题里的另一半是过程:上电后多久到达热平衡?停机后多久凉下来?振动平台的激励扫过固有频率时响应放大几倍?研究类型就是给求解器下的任务书,任务书写错,算得再努力也是答非所问。四类常见研究的分工一句话各自说清:稳态(Stationary)求不随时间变化的平衡解;瞬态(Time Dependent)求物理量随时间的演化轨迹;
本节摘要:研究类型回答的是"你想问方程组什么问题":稳态问平衡、瞬态问过程、特征频率问天性、频域问响应。本节讲清四类研究的适用场景与容差约定,并用热时间常数给铜排补一个上电六千秒的瞬态仿真。
前五章你一直用稳态研究,它回答"系统最终停在哪个平衡态"。但工程问题里的另一半是过程:上电后多久到达热平衡?停机后多久凉下来?振动平台的激励扫过固有频率时响应放大几倍?研究类型就是给求解器下的任务书,任务书写错,算得再努力也是答非所问。四类常见研究的分工一句话各自说清:稳态(Stationary)求不随时间变化的平衡解;瞬态(Time Dependent)求物理量随时间的演化轨迹;特征频率(Eigenfrequency)求系统天然的振动模式与频率,不需要任何载荷;频域(Frequency Domain)求系统在单频谐波激励下的稳态响应,扫频即得频响曲线。
选择时的判断口诀按问题来:只关心最终温升,稳态;关心到达温升的时间与过程峰值,瞬态;结构会不会共振,特征频率;交流激励下的损耗与响应,频域。注意频域研究内含谐波假设——所有量按正弦振荡且只保留基频,开关电源的方波激励要先做傅里叶分解再对关心的谐波逐个算,这是新手高频损耗算错的常见源头。
稳态研究默认相对容差 0.001,含义是非线性迭代的残差要降到初始残差的千分之一,第 1 章日志里那列一路下降的误差数字干的就是这件事。瞬态研究的默认相对容差是 0.01,且多了一层时间步误差控制:求解器在每一步里估算时间离散误差,超了就自动减半步长,富余就放大——所以瞬态仿真的"步长"是算出来的不是设出来的,你在设置里看到的只是输出时步。输出时步决定存哪些时刻的解供后处理,内部步长由误差控制机制自主调节,两者各管各的:把输出步设密不会让计算更精确,只会让文件更大;想提精度要动的是相对容差与步进严格程度。
瞬态设置的具体操作走一遍:右键"研究"→ 添加"瞬态"研究或直接改研究步骤;时间栏填 range(0,300,6000)——零到六千秒、每三百秒输出一个解。步进方式选"自由"是默认,求解器按误差自调;"中间"强制在每个输出间隔至少再算若干内部步;"严格"最重,输出时刻本身也作为精确解算点,只在关心特定时刻精确值时用。
瞬态仿真开多长的窗口?用集总参数的热时间常数估:τ = ρ·Cp·V/(h·A),热容除以散热能力。铜排代入:质量 ρ·V = 8960 × 1.5e-4 ≈ 1.34 kg,热容 385 J/(kg·K) 乘上去约 516 J/K;散热能力 h·A = 10 × 0.037 ≈ 0.37 W/K;τ ≈ 516/0.37 ≈ 1400 s,约二十三分钟。上电后温度按指数逼近平衡,三到五个时间常数即基本到位,所以六千秒窗口足够看到铜排爬到终温——这就是 6.1 开头那个 range 里 6000 的来历。算瞬态前先花一分钟算 τ,窗口开得过大浪费步数、过小看不到平台,这把口袋尺子对所有一阶热问题都适用。
特征频率研究在 4.1 已经操练过:给约束、不给载荷、填模态数,返回一张频率表与振型库。它有个兄弟用法常被忽略——电容电感的特征模态分析,电磁腔体的谐振频率同样用特征值求解器找,思路与结构模态如出一辙。
研究节点支持多步骤串联,这让"组合技"成为日常。其一是稳态初始化:瞬态计算前先跑一步稳态(比如先把流场算到充分发展),瞬态步骤的初值设置里引用它的解,起步就站在平衡点附近,前几百秒的冗长爬坡直接跳过。其二是频域前挂稳态:先算直流偏置下的温度场,再在热变形的模型上扫频算响应,偏置与交流一锅出。设置位置都在瞬态或频域步骤的"因变量初值"里,把"解"字段从默认的零换成指定研究的解——两分钟的操作,常常省下半小时的无效计算。
💡 瞬态算完先画关键量随时间的曲线再看云图:曲线没到平台,说明窗口太短;曲线抖成锯齿,说明容差太松或输出步太密,先修任务书再谈结论。
| 研究类型 | 求解什么 | 典型场景 | 注意事项 |
|---|---|---|---|
| 稳态 | 不随时间变化的解 | 散热器稳态温度、结构静力 | 无需初始值 |
| 瞬态 | 随时间变化的解 | 加热过程、脉冲响应 | 需要初始值和时间步控制 |
| 特征频率 | 固有频率与模态 | 结构振动、谐振腔 | 线性化问题 |
| 频域 | 正弦稳态响应 | 电磁波、交流阻抗 | 需要频率列表 |
| 参数化扫描 | 扫描参数批量求解 | 尺寸优化、材料对比 | 每个参数值完整求解一次 |
选择速查的用法:先确定你的物理问题是稳态还是瞬态(是否随时间变化)、是线性还是非线性(方程是否含未知量的乘积或指数)、是单参数还是多参数(是否需要扫描)——三个回答确定了研究类型,COMSOL 的研究节点就选好了。选错研究类型的典型症状:稳态问题选了瞬态导致算不出终值(只看了中间过程)、瞬态问题选了稳态导致初始瞬态被忽略。
以保险丝的热熔断为例走一遍瞬态求解。物理:电流通过保险丝产生焦耳热,温度升高直到达到熔点。COMSOL 设置:固体传热接口 + 电流接口,研究选瞬态,时间范围从 0 到 10 秒、步长 0.01 秒。关键设置:材料电阻率的温度系数让电阻随温度上升——这就引入了热-电正反馈,瞬态求解器需要捕捉这个加速升温过程。求解器的设置:时间步进方法用 BDF(向后差分公式,COMSOL 默认),最大步长限制在 0.01 秒以避免跳过熔断时刻,容差设为物理量的合理精度。后处理:画温度 vs 时间的曲线看升温速率何时加速(正反馈生效的标志)、最大温度何时达到熔点(熔断时刻)。这个案例的深层教学点是:瞬态求解的难度不在求解器设置,在于物理的正反馈是否被正确建模——物理对了,COMSOL 的自适应时间步会自动处理好数值问题。