★ FEATURED ARTICLE 仿真这件事做静力学的人多能把动力学跑通、并且跑出来的数敢拿去支撑设计决策的人少。SolidWorks Simulation 这套工具的门槛其实不高点几下就能出彩图难的是知道什么时候该用哪个算例类型、参数凭什么这么填、结果为什么和实测对不上。这篇内容就是想把这层窗户纸捅破从最基础的线性静力学讲起说清楚它的能力边界在哪然后一步步过渡到模态、谐响应、瞬态动力学这些真正吃功夫的模块。核心关键词是 SolidWorks Simulation、静力学、动力学、模态分析、谐响应、瞬态仿真。适合已经能独立完成一个静力学算例、但一碰到振动和冲击就心里没底的朋友也适合带项目、需要判断仿真结论可不可信的技术负责人。1. 从静力学跨到动力学思路该怎么搭1.1 静力学算得准的部分和它算不了的部分很多人对仿真失去信任是从静力学算出来很漂亮装到设备上一跑就裂开始的。问题往往不在软件而在于把静力学当成了万能工具。线性静力学的底层假设有四条材料线弹性、变形足够小小变形假设、载荷恒定不随时间变化、忽略惯性力和阻尼。这四条只要破一条结果的参考价值就打折扣。我通常用一句话跟新人解释静力学回答的是慢慢压上去它扛不扛得住动力学回答的是突然来一下或者一直在抖它受不受得了。一个承重支架、一个静态放置的机架用静力学加安全系数完全够用。但同样是这个支架如果装在往复运动的机构上、或者旁边有台电机以固定转速旋转那静力学给出的应力值就只是整个故事的第一页。举一个典型的误判场景某设备的悬臂支架静力学算下来最大应力 60 MPa材料屈服 235 MPa安全系数接近 4看着很放心。实际运行三个月焊缝处出现裂纹。回头做模态分析才发现支架一阶固有频率 27 Hz而现场电机的激励频率在 25 Hz 附近长期共振导致交变应力幅值虽然只有十几兆帕但循环次数上去了疲劳自然找上门。这个例子里静力学没有算错它只是没被问对问题。所以判断该不该上动力学我一般看三个信号载荷是否随时间快速变化、结构是否存在明显的周期性激励源、失效模式是否表现为疲劳或松动而非一次性屈服。命中任意一条就该考虑往动力学走。1.2 动力学几个分支各自解决什么问题SolidWorks Simulation 里能用的动力学相关算例有好几种新手最容易犯的错是抓到哪个用哪个。下面这张表是我自己整理的选择逻辑基本上覆盖了日常八成以上的需求。算例类型输入需要什么输出给什么典型场景相对耗时线性模态材料、约束、网格固有频率、振型、参与系数避开共振、判断刚度薄弱环节低预应力模态模态基础上加载荷受力状态下的频率受拉薄膜、受压结构、旋转件中低谐响应模态 激励幅值/相位 阻尼频域响应曲线、峰值响应电机、风机、路面激励中瞬态动力学模态或直接积分 时域载荷曲线时域应力/位移历程冲击、跌落、启停过程高随机振动功率谱密度 PSD统计意义下的响应 RMS运输振动、宽频激励中高响应谱加速度谱峰值响应包络抗震、冲击规范校核低选型的核心是一句话能用频域解决的别用时域能用模态叠加的别用直接积分。谐响应背后跑的还是模态叠加计算成本比瞬态低一个数量级只要激励是稳态周期性的就没必要上瞬态。反过来如果载荷是一次性的脉冲比如跌落、碰撞、点火冲击频域方法从原理上就不适用只能老老实实做时域积分。1.3 建模降阶别一上来就全尺寸精细模型我见过太多人第一次做动力学直接拿整机装配体开算几十万个单元、几百个接触对跑一晚上出一堆警告。动力学对自由度的敏感度比静力学高得多——模态分析要做特征值求解计算量随自由度呈平方甚至立方增长而瞬态还要在每个时间步上重复求解。我的处理顺序一般是先确定分析目标关注哪个部件把这个部件从装配体中摘出来用远程质量、质量点、等效刚度边界代替周边结构。拿不准边界刚度的时候宁可做两套极端假设全刚性约束 / 加一个已知刚度的弹性支撑看频率变化的区间有多大。如果两套假设下频率差不到 5%说明边界不敏感可以放心简化如果差出 30%那说明这个连接本身就是模型的关键必须老老实实建出来。对称性也是常被忽略的降阶手段。轴对称、平面对称、周期对称的结构用对称边界能砍掉一半甚至四分之三的自由度而且对模态结果几乎没有损失。唯一要注意的是对称面只能约束法向位移和两个面内转动别把面内位移也锁死否则频率会虚高。2. 网格、连接、约束动力学模型真正的地基2.1 网格策略不能照搬静力学那一套静力学里大家习惯看应力集中所以在圆角、孔边做局部加密。动力学完全不同它的网格需求是全局性的振型要被准确表达结构在每一个波长范围内至少得有几个单元。经验做法是按关注的最高阶模态反推。假设你最关心到 500 Hz结构里的等效声速对钢大约是 5000 m/s 量级的纵波速度弯曲波更低那么 500 Hz 对应的弯曲波长在几十毫米量级单元尺寸控制在波长的 1/6 到 1/10也就是几毫米。这个尺寸往往比静力学用的粗网格还细所以动力学模型通常更重这也是为什么降阶特别重要。单元类型上实体结构默认用二阶四面体抛物线单元问题不大但壳单元在薄壁件上效率高得多。一块 2 mm 厚的钣金用实体建要三层单元才能表达弯曲用壳单元一层就够自由度差出好几倍而弯曲模态频率的误差通常在 2% 以内。我一般的判断标准是厚度与最小特征尺寸之比小于 1:10就用壳。有一点必须强调用模态频率做网格收敛性检查比看应力可靠得多。同一模型粗中细三套网格一阶频率变化小于 1%就可以停手了。应力收敛慢、还容易在奇异点卡住不适合当作网格够不够的判据。2.2 连接和接触是频率误差的最大来源装配体的动力学精度八成死在连接上。SolidWorks Simulation 的接合类型主要分两类 bonded bonded 接合、接点和接触无穿透、允许分离。这个区别在模态分析里尤其致命。模态分析本质是线性特征值问题接触非线性没法直接处理软件的做法是把接触状态在初始时刻冻结——也就是按当前几何位置判断哪里贴合之后就一直保持这个状态。如果实际工况中零件会分离、会滑移那模态结果就得打问号。螺栓连接是最典型的坑。全 bonded 把两个零件焊死局部刚度会明显偏大高阶模态尤其失真。我的做法是如果螺栓预紧足够、界面不会分离可以用 bonded 加一个等效的接触面积取螺栓头承压面而不是整个贴合面全 bond如果关心的是整体低频模态前几阶连接细节影响有限简化可接受如果要算到几百赫兹以上就得考虑用梁单元模拟螺栓杆、或者干脆引入接触非线性做预应力模态。提示模态分析算完一定要回头看一眼振型图。如果某一阶振型表现为两个零件之间相对错动而模型里它们又是全 bonded 的这一阶基本可以判定是虚假模态现实里不会以那个频率出现。2.3 约束不要过定位静力学里多加点约束顶多让结果偏保守动力学里约束过定位会直接把频率抬高。判断方法很朴素看看约束是不是限制了本该自由的方向。自由模态不加任何约束是另一条路很多人不知道它的价值。做自由模态前六阶应该是接近零的刚体模态三个平动、三个转动如果算出来是 0.3 Hz、0.5 Hz 这种量级说明模型有异常——要么有零件悬空要么网格里有断开。刚体模态不接近零说明模型散了。这个检查几乎是零成本的我每次建完装配体都会先跑一次自由模态验一验。约束位置的选取也有讲究。把约束加在一个小面上会在约束附近制造人工的局部刚度产生远高于实际的高频模态。我通常把约束面稍微扩大一点或者用虚拟壁配合弹性刚度更贴近真实的连接状态。3. 模态分析实操先拿到系统的身份证3.1 完整设置流程与关键参数模态算例的设置流程本身不复杂但几个参数必须想清楚第一是提取阶数。软件默认提取 5 阶这远远不够。原则是提取到关注频率上限的 2 倍。比如你关心 200 Hz 以内的响应那就设置提取所有 400 Hz 以下的模态。原因是谐响应和瞬态做模态叠加时高于激励频率的模态虽然不共振但会提供静刚度贡献截断太狠会让低频响应偏大。第二是频率范围 vs 阶数。SolidWorks Simulation 里可以按阶数提取也可以设频率上限。我倾向设频率上限配合足够的阶数上限避免在小模型上提取过多无意义的局部模态。第三是求解器选择。默认的 FFEPlus 迭代求解器在自由度大时速度快但对密集模态频率接近容易漏根直接稀疏求解器更稳但吃内存。如果发现两阶频率非常接近、或者结果和另一套网格差异大换直接求解器复算一次是值得的。完整流程我一般这么走新建算例选频率材料按实际填写不确定的话做敏感性分析模量 ±10% 看频率变化连接关系逐条过一遍bonded 的地方确认接触面积合理夹具按实际安装状态加能不加过约束就不加网格先粗后细用一阶频率验证收敛运行看频率列表、振型动画、参与系数三样东西。3.2 模态结果到底该怎么读频率列表只是结果的三分之一另外两样才是判断依据。振型动画要一阶一阶看判断是整体模态还是局部模态。局部模态通常表现为某个薄板、某根细杆单独在抖频率往往很高对整体响应贡献有限。如果你发现前几阶全是某个小支架在振说明这个零件刚度偏弱但也可能是网格或连接造成的假象需要甄别。有效质量参与系数是判断提取够不够的硬指标。它告诉你这一阶模态在某个方向上承担了多少质量。把各方向X/Y/Z的参与质量累加达到总质量的 80% 以上说明主要的结构响应已经被捕捉到了低于 60% 说明还有大量质量在更高频的模态里运动需要提高提取阶数。这个数据在结果里能直接查到但很多人从来没点开过。还有一个容易被忽略的点同一阶模态可能在不同方向上都有参与。如果一个结构某个方向的累积参与质量始终上不去往往意味着那个方向缺少约束结构在那个方向是飘的这本身就是个有价值的发现。3.3 预应力模态什么时候必须做结构在承载状态下刚度会变化这叫应力刚化或者应力软化。一根拉紧的钢丝横向振动频率比松弛状态高得多一块受压的薄板压力接近屈曲载荷时频率趋近于零。这两种情况下不加载荷直接算模态误差是量级级别的。判断标准结构中的应力水平与材料模量之比超过 1% 到 5%且应力以轴向拉压或薄膜应力为主时就该做预应力模态。常见的必须做的情况包括旋转件离心力引起的刚化、承受内压的薄壁容器、张紧的索和膜结构、细长受压杆件。SolidWorks Simulation 里做预应力模态的操作是先建一个静力学算例把载荷加上算完再建频率算例在属性里把静力学算例设为预应力来源。注意两边的网格和连接必须一致否则预应力传递不过去。3.4 模态和实测对不上怎么查频率算出来偏高是最常见的现象原因按可能性排序大概是约束过刚把弹性连接做成了刚性固定连接全部 bonded界面刚度虚高材料模量取值偏大尤其是塑料、复合材料忽略了附加质量线缆、涂层、紧固件网格偏粗导致刚度虚高这个方向也会偏高。算出来偏低通常是质量偏大密度、附加质量重复计入或者约束不足结构处于自由状态。我自己的经验是金属焊接结构模态仿真和实测的偏差控制在 5% 以内是合理的螺栓连接较多的装配体10% 到 15% 也属正常含橡胶件、含大量非结构附件的30% 偏差都不稀奇。拿这个量级去判断结果可不可信比追求和实测一模一样现实得多。4. 谐响应与瞬态载荷怎么加、步长怎么定4.1 谐响应分析的关键三件事谐响应回答的是结构在某个频率的正弦激励下响应多大它的设置就三件事频段、扫频密度、阻尼。频段要覆盖你关心的范围同时把主要模态包进去。我一般取 0 到 2 倍关注频率或者至少覆盖前 3 到 5 阶模态。频段设太窄容易漏掉低频的整体响应设太宽计算点浪费在高频无响应区。扫频密度是最容易出问题的地方。谐响应的曲线是逐个频率点算出来的如果点太稀峰值会被跳过——你看到的曲线很平滑实际峰值可能是它的三倍。经验做法是在每个共振峰附近加密先做一次宽带、稀疏的扫描比如 20 到 50 个点找出峰值位置再在峰值附近 ±20% 的范围内加密到半功率带宽内至少 5 到 10 个点。半功率带宽的估算带宽 ≈ 2ζf其中 ζ 是阻尼比。假设一阶频率 45 Hz、阻尼比 2%带宽约 1.8 Hz那么在 44 到 46 Hz 之间至少要有 5 个点步长 0.4 Hz 左右。这个计算不难但能省掉很多曲线看起来不对劲的返工。阻尼的选取在第 4.3 节单独说这里先记住一句话谐响应的峰值高度几乎完全由阻尼决定阻尼填错结果就没有意义。载荷施加方面谐响应需要的是幅值和相位。如果激励是旋转机械的不平衡力幅值随转速平方增长F m·e·ω²这意味着不同频率下的激励幅值不是常数必须按频率给一条载荷曲线而不是填一个固定值。这一点很多人第一次做会忽略导致高频段的响应被严重低估。4.2 瞬态动力学时间步长是生死线瞬态的精度七成取决于时间步长。步长太大高频响应被过滤掉、能量凭空消失结果看着很稳但其实是假的步长太小计算时间飞涨。理论依据是采样定理的工程化版本时间步长应小于最高关注频率对应周期的 1/20。最高关注频率怎么定不是看激励频率而是看结构的响应频率——一般是激励频率和前几阶固有频率中较高的那个。假设结构一阶频率 120 Hz激励是 20 Hz 的脉冲那么响应里会激起 120 Hz 的振荡周期 8.3 ms步长取 0.4 ms 以内。如果做的是冲击、跌落这类含接触的瞬态还要额外考虑接触的收敛步长往往要更小同时开启自动时间步。总时长要覆盖完整的过程并且留出余量。冲击结束后结构还会自由振荡几个周期如果你在振荡衰减完之前就截断能量没有耗散完结果不完整。我一般取激励持续时间 3 到 5 个主要振荡周期。载荷曲线用表格定义可以导入实测数据。这里有个细节载荷曲线的时间分辨率不要比求解步长更粗否则求解器在两点之间线性插值会人为削掉峰值。实测数据采样率低的话宁可做一次平滑重采样。另外瞬态的初始条件不能忽视。默认从静止开始如果你的实际工况是设备已经在稳态运行突然受到一次冲击那就得先算一个稳态作为初始状态否则起步的瞬态会污染结果。4.3 阻尼设置最常被瞎填、也最能决定成败的参数阻尼是动力仿真里唯一一个没有可靠理论值、只能靠经验区间的参数但它对结果的影响又极大。我见过太多人在这里随手填个 0.05 就开算然后拿结果去说服别人。常用阻尼比的参考区间结构类型与连接方式阻尼比参考范围全焊接钢结构0.5% – 2%螺栓连接、铆接结构2% – 5%含少量橡胶隔振件的结构5% – 10%大量橡胶、聚合物连接10% – 20%有阻尼器、摩擦副的结构需专项测定如果做的是模态叠加法SolidWorks Simulation 里常用的是模态阻尼或 Rayleigh 阻尼。Rayleigh 阻尼的两个系数 α质量比例项和 β刚度比例项可以按两个已知频率反算α 2ζ·ω₁·ω₂ / (ω₁ ω₂) β 2ζ / (ω₁ ω₂)举个例子结构一阶频率 45 Hz、二阶 180 Hz取 ζ 2%。先换算角频率ω₁ 2π×45 282.7 rad/sω₂ 2π×180 1131 rad/s。α 2 × 0.02 × 282.7 × 1131 / (282.7 1131) 12789 / 1413.7 ≈ 9.05 β 2 × 0.02 / (282.7 1131) 0.04 / 1413.7 ≈ 2.83×10⁻⁵这两个系数代入后在 45 到 180 Hz 的区间内实际阻尼比大致维持在 2% 附近超出这个区间后阻尼比会偏离高频段可能偏大、低频段偏小。所以用 Rayleigh 阻尼时两个参考频率一定要选在你最关心的频段两端不要随便取。注意阻尼比和响应幅值是反比关系。阻尼从 2% 改到 1%共振峰值理论上翻倍。所以做谐响应时阻尼的不确定性要作为结论的一部分说出来别只报一个数。5. 常见问题与排查速查表5.1 异常现象与对应原因下面这张表基本是我这些年踩坑的浓缩版遇到问题时从上到下对照命中率不低。现象最可能的原因处理方式一阶频率明显高于实测约束过刚、全 bonded、模量偏大改用弹性支撑、按接触面积 bond、做模量敏感性分析一阶频率明显低于实测密度偏大、附加质量重复计入核对材料密度检查是否有重复的质量点前六阶频率不趋近于零自由模态模型中有零件悬空或网格断裂检查接触对与网格连通性振型里有零件飞出去该零件未连接补接合或接触谐响应曲线过于平滑、无明显峰值扫频点太稀跳过共振峰峰值附近加密步长取半功率带宽的 1/5 以内谐响应峰值异常高阻尼填得过小如填成 0.0005按结构类型重取阻尼比瞬态结果出现高频毛刺或发散时间步长过大或过小步长控制在 1/(20f_max) 量级开启自动步长瞬态能量不守恒、响应衰减异常数值阻尼过大积分参数检查积分格式HHT 参数取默认或调小应力最大值出现在约束点应力奇异点约束、尖角忽略该点取附近节点的收敛值或用圆角/扩大约束面求解时间过长接触对过多、直接积分、网格过密先用模态叠加法非必要不引入接触结果和静力学差异很大惯性效应主导或载荷施加方式不同对比载荷幅值与频率确认静力等效是否成立5.2 一份可以照着走的判断清单算完动力学结果我固定会问自己四个问题任何一个答不上来就先别拿去汇报。第一个问题这个频率合理吗拿手算做一个粗略校核。一个悬臂梁的一阶频率可以按 f ≈ 0.56·√(EI/(mL⁴))除以 2π 之后的工程近似形式估个量级或者更简单——用静力学算出的柔度反推给结构加单位力得到位移 δ则一阶频率 f ≈ 1/(2π√(δm))。这个估算只要一次静力学计算五分钟能出结果用来判断模态结果是否离谱足够。第二个问题参与质量够不够前面说过80% 是门槛。第三个问题边界条件在真实结构里对应什么模型里那块被完全固定的面现场到底是怎么装的垫了几个垫片螺栓拧多紧这些细节决定了你算出来的频率是 27 Hz 还是 45 Hz。第四个问题结论对参数有多敏感把模量、阻尼、约束刚度各动 ±20%看结论变不变。如果某个参数一动结论就翻盘那这个结论的可靠性就取决于那个参数测得准不准必须说明清楚。6. 一些从实际项目里攒下来的体会做动力学这几年最大的感受是建模的取巧之处最终都会以结果不可信的形式还回来。静力学里偷的懒影响往往局限在局部应力动力学里偷的懒会直接体现在频率和响应的量级上而且是系统性的偏移不是局部误差。另一个体会是关于精度的期待管理。有限元动力学的绝对精度其实不高模态频率能到 5% 以内已经算优秀响应幅值的误差动辄 30% 以上。所以它的正确用法不是预测一个精确数字而是比较两个方案谁更好判断结构离共振有多远找出最薄弱的环节。把仿真当作趋势工具和筛选工具而不是验收依据心态会稳很多。最后分享一个我自己一直在用的小技巧任何动力学模型先跑自由模态看刚体频率再跑约束模态看前几阶振型最后才去加载荷做响应。这三步的顺序不要颠倒因为前两步是免费的体检能筛掉八成的建模错误。跳过体检直接上瞬态往往跑了几小时才发现模型是散的那种感觉试过一次就记住了。
阅读完成 · 觉得有帮助？
