ARTICLE · 1084862
开箱:软件的零件柜——《仿真理论手册》逐章精读
仿 真 精 读 · 第 六 站（完结篇）
开箱：软件的零件柜
——《仿真理论手册》逐章精读：从 TRUSS3D 到 SHELL3，六件套各自的脾气
五站走完，理论全齐：一个点的力学、一串方程的流水线、一套振动、一次翻脸、一本热账。但按钮背后还有最后一块拼图——单元库。手册第六章是个「货架清单」：六个单元，每个都配好四件套——刚度矩阵、变分泛函、形函数、数值积分。
这一章读起来像逛五金店：不起眼，但每一件都是前五章理论的具体化身。今天带各位把六个货架逐个开箱——看完这期，你在软件里点选「杆单元」「梁单元」「壳单元」时，就知道自己把谁的脾气请进了模型。
🛒 选单元像选工具：扳手拧不了十字螺丝。单元没有「最好的」，只有「最合适的」——选错的代价，软件不会当场告诉你，它会在结果里慢慢告诉你。
本 章 路 线 图（终 点 站）
第一站第一章 基本关系（已更新）
第二站第二章 有限元法（已更新）
第三站第三章 结构振动频率（已更新）
第四站第四章 结构屈曲（已更新）
第五站第五章 热传导（已更新）
第六站第六章 单元库（本次 · 完结）
壹 · 一根杆：TRUSS3D
柜子里最简单的家伙。两节点、每节点三个平移自由度，只传轴力——像铰接的连杆，不认弯矩。形函数是第二章见过的线性插值（每方向 u = (1−x/ΔL)u₁ + (x/ΔL)u₂），刚度、质量、重力载荷全部「精确、零积分点」——线性函数被高斯积分秒杀，都不用上积分点。
图 6-1原书 Figure 15：TRUSS3D 线性三维桁架单元。图源：SOLIDWORKS Simulation 理论手册
小家伙也不含糊：应力公式 σx = E(Δux/ΔL − αΔT) 连热载荷都算好了（杆子被热胀顶住会受压，这钱它认账）。方向变换矩阵 T 只装方向余弦 ΔX/ΔL、ΔY/ΔL、ΔZ/ΔL——杆往哪摆，刚度就往哪摆。
🔩 桁架的世界只有拉和压，像只会发「嗯」「啊」两种消息的聊天对象——简单，但大桥和铁塔全靠它们撑着。
贰 · 一根梁：BEAM3D，与它的「第三点」
梁单元是柜子里的全能选手：每节点 6 个自由度（3 平移＋3 转动），2 节点共 12 阶刚度阵。第二章的梁理论在这里完整落地：轴向、双轴弯曲、扭转、剪切，一个不落。考虑剪切变形时，插值函数里冒出两个参数 φy、φz = 24(1+ν)Ī/(k̄AL²)——梁又细又长时它们趋近于零，自动退化为经典欧拉梁；越短粗，剪切的存在感越强。一个参数，两代理论无缝衔接。
图 6-2原书 Figure 16：BEAM3D 线性三维弹性梁单元。图源：同上
真正的戏肉是定向：梁截面绕轴线怎么摆？手册的答案：再给第三个点——它和两端的连线定出主平面（图 6-3）。软件里选梁单元时那条「方向参照」，就是这个点在值班。没它，工字梁的腹板朝上还是朝侧，全凭缘分。
图 6-3原书 Figure 17：偏置节点与第三点定义的梁单元坐标系。图源：同上
从截面主轴到整体坐标，刚度阵要坐五趟换乘：T1 主轴转杆轴、T2 端截面偏置、T3 杆轴转单元轴、T4 单元转整体、T5 剪力中心换手（式 5 连乘，配合图 6-3 的 LCG、LA、LB、LN 几何）。连锥形梁都照顾到了：两端截面不等时，用等效特性 A = (1/3)(A₁+√(A₁A₂)+A₂) 等一组公式折算——老一辈工程师手算变截面杆的记忆，被浓缩成三行代数。
🧭 「第三点」就像自拍的参照角度：不给参照，截面朝向全靠缘分。在软件里给梁点方向时手别抖——它一抖，整根梁的弯曲轴跟着抖。
叁 · 硬汉与垫片：RBAR 与 SPRING
RBAR 刚性杆（图 6-4）：两节点之间焊死——平移、转动全部满足刚体关系（六条约束方程）。有限元里「焊死」不能真焊，手册用罚函数：给每条约束挂一个量级 10¹⁰ 的大数刚度，谁敢违反，就吃一张天价罚单。上一站热分析里那个「大导热率钉温度」的技巧，在这里正式报上学名。
图 6-4原书 Figure 18：RBAR 两节点刚性杆。图源：同上
SPRING 弹簧单元（图 6-5）：两节点之间只给两根「筋」——轴向刚度 kaxial 和扭转刚度 ktorsional，刚度阵直接手写，连积分都省了。它是仿真里的万能垫片：地脚螺栓的柔性、减振器的刚度、电机轴的扭振，都能拿它「垫」一下。
图 6-5原书 Figure 19：SPRING 弹簧单元。图源：同上
🪨 罚函数的哲学：不讲道理，只讲力气——10¹⁰ 的「罚单」一开，约束比真焊死还硬。弹簧的哲学正好相反：能屈能伸，刚度你说了算。
肆 · 两片三角壳：SHELL3T 与 SHELL3
壳单元是柜子里的压轴货，而且一柜两派。SHELL3T 厚壳（图 6-6）走 Mindlin 路线：三节点、每节点 6 自由度、18 阶刚度阵，转角 Ψ 独立出场，横向剪切留着余地。形函数用三角形的面积坐标 L₁、L₂、L₃（在某个顶点是 1、在对边是 0 的「三分天下」函数）；面内能量精确积分，弯曲能量用一个高斯点——第二章讲过的减缩积分，在这块单元上正好上岗。
图 6-6原书 Figure 20：SHELL3T 三角形厚壳单元。图源：同上
SHELL3 薄壳（图 6-7）走 Kirchhoff 路线：不允许横向剪切，法线变形后仍垂直中面。它的处理堪称漂亮：内部先多配六个中点转角自由度，再用 Kirchhoff 条件（转角绑死在挠度斜率上）逐个消除——手册里那三行 ∂uz/∂x + Ψy = 0 就是「裁员通知书」。瘦完身的单元只保留角点自由度，刚度阵却精确积分。这就是文献里常说的离散 Kirchhoff 三角（DKT）一族的做派——手册参考文献第 25 条（Batoz-Bathe-Ho 1980 年那篇三角板单元研究）是它的户口本。
图 6-7原书 Figure 21：SHELL3 三角形薄壳单元。图源：同上
🚙 厚壳像 SUV：不挑路（粗网格、翘曲、剪切都扛得住）；薄壳像轿车：高速优雅（薄板精确），烂路上（粗网格＋剪切敏感）容易托底。选哪个，先看路况。
伍 · 六件套点名册
——一张表放六站功夫。选型的功课，赛前做，别赛后做。
🖥️ 软 件 对话框里的零件柜
在 SOLIDWORKS Simulation 里，右键树里的「应用单元」就是把零件指派给这六个货架：细长焊架选桁架、框骨架选横梁、薄壁件抽中面用壳、螺栓垫圈垫弹簧、需要焊死的地方用刚性杆。指派之前多想一秒「它的脾气合不合」，比事后检查十遍结果都值。
📜 史 料 零件柜里住着人名
翻到手册参考文献页会发现：前三条全是铁木辛柯——《弹性理论》《板壳理论》《弹性稳定理论》。梁单元里那个剪切修正参数，源头是他 1921 年的梁理论；这位从俄国到密歇根、再到斯坦福的工程师，被后辈称作工程力学界的奠基人之一。软件的零件柜看着冷冰冰，柜板背面全是人名：柯西、胡克、泊松、铁木辛柯、柯朗、兰索斯——每个单元里，都住着一个人的一生。
六站到此收官。回头看这条线：第一章备好一个点的三本字典；第二章把零件切成方程，装入K·q=F；第三章让结构开口唱歌；第四章算出它翻脸的门槛；第五章把温度的账算清；第六章开箱，看见理论的肉身。全书 116 页，讲的其实是一件事：你点下的每一次「运行」，都是这一百年力学的一次回放。
完结不散场：六篇全部可回看，建议按顺序二刷——第一遍看热闹，第二遍看门道。哪篇没看懂，评论区喊我；想加更的题目（非线性、疲劳、拓扑优化），也欢迎点菜。
🏁 全系列完结撒花。手册合上，按钮还在——下次点「运行」时，你已经知道那几秒钟里发生了什么。这，就是这六站想给你的东西。
全 系 列 完 结
《仿真精读 · 六站》全六篇已更新完毕
一个点的基本关系 → 有限元的流水线 → 振动频率 → 屈曲 → 热传导 → 单元库。感谢一路追更的各位；回看请戳公众号历史消息，从「断口，是应力画好的」开始。
MC话机械，聊点机器背后的事。
仿 真 精 读 · 第 6 / 6 站 · 完
本文参考：Theoretical Manual — SOLIDWORKS Simulation（Jeffrey Wong 著）第六章「The Element Library」。文中图 6-1 至 6-7 均裁自原书 Figure 15 至 21。
MC话机械，聊点机器背后的事。