★ FEATURED ARTICLE 双稳态折纸立方体这种东西玩实物的时候最直观的感受就是那两个“咔嗒”停靠点摊开来是方方正正的立方体沿着折痕一压哗啦一下就塌成另一形态中间总有一股明显的“别扭感”要翻过去。很多人第一次摸到都会问一句“它为什么能停住为什么到中间要使劲”——这正是双稳态结构的力学本质系统势能曲线存在两个极小值点中间隔着一个能量高峰外力得先把结构“推过”那个峰才能让它自己滑进另一个阱。但手上把玩归把玩真正要用于产品设计、机构验证、或是判断这个东西能不能可靠工作上万次时靠手感就不够了。你需要知道折叠到底需要多大的驱动力、哪片板材在哪个瞬间应力最大、如果能垒选得太高会不会导致折叠机构卡死——这些问题的标准答案都是在ABAQUS里把折叠过程完整跑一遍之后才拿得到的。这篇文章就按我实际做项目的顺序把ABAQUS模拟双稳态折纸立方体的完整思路讲清楚从双稳态的力学本质、建模前的关键决定到参数取值、后处理判据再到动真格时踩过的坑。适合正在做可展结构、折纸机器人、柔性机构仿真的研究生和工程师也适合想把“多稳态机构”用有限元讲清楚的初学者。本文不会贴特别复杂的公式重点放在怎么把这些概念翻译成ABAQUS里的具体操作。1. 双稳态折纸立方体的力学本质能量曲线才是主角很多人一上来就在ABAQUS里画几何体、分网格结果跑出来一坨东西既不像折叠也不像展开根本不知道问题出在哪。我建议第一步先别碰软件把双稳态的力学模型在脑子里立起来因为后面所有参数设置本质上都是在复现这张能量曲线。1.1 两个凹坑之间的翻越用势能曲线理解“咔嗒”感把整个折纸立方体看作一个力学系统它的状态可以用一个广义坐标来描述最常用的就是“折叠角”——也就是立方体某一对相邻面板之间的夹角。当这个角度从180°完全展开的立方体变化到接近0°完全折叠收拢时系统的总势能不会单调上升或下降而是呈现一条W形的双凹曲线。展开状态和折叠收拢状态分别对应曲线上的两个极小值点所以结构能稳定停在那个位置。而两者之间必然存在一个局部极大值点对应折叠角约在四五十度附近的位置这个点就是能量壁垒也叫能垒。折叠时手感觉到的那股“别扭感”实际上就是你通过外力做功把系统从第一个势能阱里拉出来、推上山头、再放它滚进第二个势能阱。这个模型虽然简单但直接回答了三个核心问题能垒的高度决定了折叠需要的峰值驱动力能垒越高越难折但一旦折过去就越不容易自己弹回来。两个阱的深浅差异决定了哪个状态更“占优”如果展开态的阱明显更深长期存放时结构倾向于摊开。峰值附近曲线的斜率决定了你推到最后一下的时候是“突然自己冲进去”还是“还要顶住慢慢压进去”这就是手感来源。1.2 为什么偏偏是折纸结构最容易做出双稳态不是所有多稳态结构都好做但折纸结构天生适合原因在于它把变形集中到了折痕上。面板可以近似看作刚体所有弯曲变形都发生在狭窄的折痕区域这使得总势能的主要贡献是折痕的弯曲能而面板自身的拉伸能和弯曲能几乎可以忽略。也正因为面板近似刚体我们可以放心地用壳单元来模拟网格也不用做得特别细。更关键的一点是折纸结构的几何非线性极其强烈。折叠过程中面板会发生大幅转动和位移相邻面板之间还会发生自接触——这些恰恰是ABAQUS这类大变形非线性求解器最擅长的场景。拿折纸立方体来说它比单条折痕的折叠复杂得多因为多个折痕的弯曲是联动的折叠路径不是唯一确定的。同一个立方体可能从这组折痕先折也可能从另一组先折不同的折痕变形顺序会导致完全不同的受力情况。仿真里如果不做任何路径约束静力学求解经常会在这类分支点附近徘徊不定这也是为什么后面我要建议用显式动力学来跑。1.3 仿真到底要回答哪几个问题明确了能量曲线的概念之后我每次做这种仿真都会在项目开始前把“要回答的问题”列成清单。这不是为了写文档而是为了决定后处理看什么极限驱动力折叠过程中载荷-位移曲线的峰值对应的力是多少这直接决定电机、气缸或手动机构的选型。稳态保持力在两个稳态位置结构抵抗外界扰动的刚度有多大如果保持力太弱稍微一颠簸结构就会自己跳状态。应力集中位置折痕根部或面板圆角处的峰值Mises应力有没有超过材料屈服点折纸机构的失效几乎都发生在折痕附近见过太多次仿真里明明峰过了、结果折痕区应力已经远超材料极限的情况。折叠路径稳定性如果初始条件有微小偏差结构是沿着设计路径走还是会扭成一个不对称的畸形状态这些问题在ABAQUS里都可以得到定量回答但前提是建模方式选对了。2. 建模前的四个决定折痕处理、单元类型、单位制和求解器有了力学模型接下来落地到ABAQUS。这一步有四个关键决定每一个都会直接影响仿真成败。我建议在打开CAE之前就全部想好不然改起来非常痛苦。2.1 折痕建模这是整个仿真的灵魂折痕是双稳态折纸结构里变形最集中的地方模拟得好不好几乎决定了整个仿真可信不可信。ABAQUS里做折痕主要有三种方案我放在一起对比方案做法优点缺点适用场景连续薄带法把折痕区域划分出一个窄条降低其弹性模量或厚度变形连续、接触自然、应力结果直观网格需要足够细分计算量增加关心折痕局部应力、需要把实物压痕对应的刚度退化考虑进来时连接器铰链法在面板间建立HINGE连接器定义扭转刚度曲线计算稳定、参数直接对应折痕弯矩-角度关系无法得到折痕区域的应力分布建模略繁琐只关心整体机构行为、不关心折痕应力时薄壳带法折痕用非常薄的壳单元处理面板用实体精度较高不同类型的单元连接麻烦折痕区域厚度变化明显、面板很厚时我绝大多数项目用的是连续薄带法。原因是折纸机构的失败往往发生在折痕区域如果折痕只是铰链连接器应力集中信息就全丢了。具体做法是面板厚度0.2mm折痕区同样厚度但降低弹性模量到面板的1/30左右——这模拟的是实物中压痕或折线预处理造成的刚度退化同时也让变形自然地集中在折痕区。折痕的宽度通常取0.5mm到1mm取决于实际折痕的压痕宽度。提示折痕区如果太窄网格尺寸必须跟着缩小否则计算量会爆炸。一个40mm的立方体折痕宽度0.5mm网格至少要在宽度方向剖两三个单元全局网格尺寸4mm的前提下折痕区网格要细化到0.25mm左右。2.2 单元类型用S4R但别迷信S4R折纸结构用壳单元是公认的标准做法首选S4R减缩积分壳元。它在面内变形不剧烈的情况下效率极高自带沙漏控制。但要注意两个坑第一减缩积分单元在发生面内严重翘曲或单元被压扁时容易出现沙漏模式——即单元看起来有变形实际上没有应变能。折纸折叠过程中单元的大转动很常见所以Section Controls里一定要打开增强沙漏控制并且在后处理里检查沙漏能ALLAE占总内能的比例。第二面板与折痕的单元尺寸最好不要差太多。如果面板区4mm、折痕区0.25mm两侧的过渡区网格会非常难看还会产生不规则的翘曲单元。我的习惯是把折痕两侧向外扩展一小段过渡区让网格尺寸从0.25mm渐变到4mm既保证精度又控制计算量。用一个40mm见方、折痕宽0.5mm的立方体来说全局网格4mm折痕区0.25mm总单元数通常在5000到12000之间这是很舒适的规模普通工作站几分钟就能跑完一折。2.3 单位制陷阱比热容、导热率和热膨胀系数的单位查一次少走一星期弯路ABAQUS没有内置单位制你输入什么数字它就按什么单位算。新手在这上面吃亏的太多了。最常用的两种自洽单位制是SI制长度m、质量kg、时间s、力N、应力Pamm制长度mm、质量tonne吨、时间s、力N、应力MPa在mm制下密度的单位是tonne/mm³。比如PET材料密度1.38×10³ kg/m³换算成mm制就是1.38×10⁻⁹ tonne/mm³别写成1.38×10⁻⁶ kg/mm³那差了整整1000倍。很多做热-力耦合的人问过ABAQUS里比热容、导热率、热膨胀系数到底用啥单位。答案其实一句话跟长度单位保持自洽就行。具体查起来是这样物理量SI制m-kg-smm制mm-tonne-s说明比热容J/(kg·K)J/(tonne·K)数值上比SI制大1000倍导热率W/(m·K)W/(mm·K)用灵敏单位时要换算不要直接抄热膨胀系数1/K1/K无量纲比值两种制下数值相同用钢材在mm制下举例比热容约4.6×10² J/(kg·K)在ABAQUS mm制里要填4.6×10⁵ J/(tonne·K)导热率约45 W/(m·K)在mm制里要填4.5×10⁻² W/(mm·K)。折纸模拟本身一般不做热分析但如果你的工况涉及热定型、热压工艺单位没换算对算出来的温度场整个都是错的。2.4 求解器用Explicit还是Standard这是做折纸仿真必问的问题。折叠过程要克服能垒本质上要经历一个由稳定到失稳再到稳定的不稳定阶段隐式求解在失稳点附近会遇到刚度矩阵奇异收敛非常困难。我的经验是Explicit显式动力学跑折叠过程的首选。它不做隐式迭代天然能通过失稳点处理自接触也远比Standard鲁棒。代价是必须检查准静态性——折叠速度不能太快动能不能占主导。做准静态判断的核心指标就是动能与内能之比后面细说。Standard隐式适合分别提取两个稳态附近的刚度特性和应力场但不适合直接跨过能垒。如果必须用Standard跨峰得开自动稳定Automatic Stabilization或者一点点加阻尼非常折磨。所以我的推荐路线是用Explicit跑完整折叠-展开循环拿到关键载荷-位移曲线再用Standard单独算两个稳态点做精确的模态和刚度校核。大多数项目其实只跑Explicit就够了。3. 折叠仿真的实际参数怎么定以40mm立方体为例理论讲完了这部分是我实际建摸时的具体参数取值。每个数值都有出处不是随手填的我会写清背后的理由。3.1 几何与网格划分把折痕分区清清爽爽地剖出来在CAE里导入立方体几何后第一件事就是把每个面板的折痕区域partition出来。以40mm的立方体为例每个面是40×40的正方形面板每个面板的四条边各有一条0.5mm宽的折痕带相邻面板共享同一条折痕带。这样整个立方体就分成了24块面板区若干条折痕区再加上折痕向面板内的过渡区。分区的意义是让不同区域能分配不同的材料属性。面板区用PET全刚度折痕区用折痕刚度这就避免了用tie去连接不同材料的部件——整块几何用同一个part再按区域分配截面属性接触和网格都干净很多。网格划分的时候我建议这么处理全局种子4mm折痕区种子0.25mm过渡区让ABAQUS自动渐变。单元类型S4RSection Controls里开启增强沙漏控制。壳厚方向默认5个积分点就够了因为面板主要是弯曲变形不是接触挤压。前一版有个项目图省事直接把折痕区也划成4mm网格结果折痕方向的弯曲刚度被严重高估折叠力算出来比实测大了近3倍。这个案例很典型折痕宽度0.5mm却只用4mm的单元去表达相当于把0.5mm宽的柔性铰链糊成一整片4mm宽的软板能垒当然整体漂移。3.2 材料与接触PET面板、折痕退化参数和自接触材料参数按PET聚酯薄膜来设置这是折纸原型最常见的材料。面板区直接给PET弹性参数。折痕区的折纸刚度我一般按面板弹性模量的1/30到1/50来取前提是实物折痕真的做过压痕处理。如果只是光板直接折折痕刚度可能和面板只差3到5倍那就不是折纸了是硬掰塑料板。材料卡用inp片段来说明最直观** 面板区材料 *Material, namePET_Panel *Density 1.38e-9 *Elastic 3000.0, 0.36 ** 折痕区材料 *Material, namePET_Crease *Density 1.38e-9 *Elastic 100.0, 0.36接触设置是整个仿真里最容易被低估的一步。折纸折叠到后半程必然出现面板压到面板、面板压到折痕的情况所以接触必须包含自接触。我的标准配置是*Contact *Contact Inclusions, ALL EXTERIOR, SELF *Contact Property Assignment, , , FOLD_CONTACT ** 接触属性 *Surface Interaction, nameFOLD_CONTACT *Friction 0.15 *Surface Behavior, pressure-overclosureHARD摩擦系数取0.15是为了模拟PET表面光滑的自润滑特性。如果你做的是纸张材料摩擦系数起码给到0.3以上纸张的接触行为对摩擦非常敏感。如果模拟结果中折叠力偏高先别急着怀疑材料刚度先想想摩擦是不是给大了。3.3 边界条件、加载方式与分析步设置别“手捏”要“机压”边界条件的设计直接决定结果可重复性。用手捏的加载方式是无法在仿真里稳定复现的因为手指接触面、施力方向、捏的先后顺序都有高度不确定性仿真的结果就会带上一堆你自己都没意识到的边界噪声。我推荐的做法是在立方体底面的某个稳定接触面上创建一个刚性平板约束刚性平板的所有自由度让立方体底面刚好架在上面。顶面再建一个刚性平板作为加载板通过参考点施加位移载荷。这样加载过程是平面对平面的匀质压入几何对称性和接触稳定性都极好后处理提取支反力也简单。加载的幅值曲线用Smooth StepABAQUS自带的光滑阶跃幅值在一个分析步内平滑地增加位移避免启动时刻的冲击载荷。分析步时长我给0.03秒左右网格0.25mm时Explicit稳定时间增量大约是1e-8秒量级加固定质量缩放把目标增量控制在1e-7到2e-7秒整个折叠过程计算量是可控的。smooth step幅值、0.03秒的分析步长、固定质量缩放——这三个搭配对我的40mm模型来说折叠力曲线平稳无震荡。如果你发现自己的模型在载荷加载初期出现高频振荡往往就是幅值曲线太陡、加载时间太短把0.03秒改成0.1秒重跑一次对比一下通常能直接消除振荡。4. 后处理判据不是变形好看就叫双稳态跑完计算只是第一步怎么从海量的输出数据里判断“这个模型到底是不是合格的双稳态结构”是需要一套明确判据的。只看变形动画是远远不够的我见过太多人拿着一个看起来很丝滑的折叠动画就宣布仿真成功但能量曲线一拉出来全是一毛一样的单阱结构——那根本就不是双稳态。4.1 能量校验先查动能再谈准静态Explicit算出来的结果要可信前提是整个过程接近准静态。怎么判断看能量历史曲线。在历史变量输出里把ALLIE内能、ALLKE动能、ALLSE应变能和ALLAE沙漏能全部输出*Output, history, frequency100 *Energy Output ALLIE, ALLKE, ALLSE, ALLAE, ETOTAL判断标准有两条ALLKE动能要远小于ALLIE内能。业界常用的经验阈值是动能占内能的5%以下。动能占比超过10%说明惯性力已经在主导折叠过程折叠速度太快算出来的载荷峰值会明显偏大。ALLAE沙漏能占ALLIE的比例应低于5%。沙漏能过高说明网格里的零能模式在主导变形结果已经没有物理意义了。如果检查出动能偏高降低加载速度或者减弱质量缩放重跑一遍即可。如果沙漏能偏高回到Section Controls里增强沙漏刚度或者把折痕区网格再细化。4.2 载荷-位移曲线能垒高度和负刚度段才是核心证据能量检查通过后接下来是提取载荷-位移曲线。在加载板参考点上输出反力RF和位移U画出来的RF-U曲线就是整个双稳态行为的核心证据。双稳态的RF-U曲线应包含以下特征曲线从零开始力随位移上升到达第一个峰值。这个峰值就是你克服能垒需要的力对应手感的“最难推的那一下”。峰值之后曲线应有一段负刚度区力随位移增大而下降甚至变成负值。负刚度意味着结构过了失稳点后不再需要外力驱动自己在“往阱里滚”。这段负刚度区就是“咔嗒”感的数学表达。到达第二稳态后力重新回到零点附近。这时候结构已经停进第二个势能阱不需要外力维持稳定。如果总势能两个阱相等RF-U曲线正反两个方向的峰值和过零位置几乎对称如果不对称说明两个稳态一个更稳一个较亚稳长期存放会倾向于停在深阱侧。我每次都会把峰值力和峰位角标出来最大峰值力直接拿去做驱动机构选型——这是我做双稳态折纸仿真最实用的一步。4.3 应力与变形折痕区的Mises应力不能超过材料极限能量曲线和载荷曲线都符合预期后最后看应力场。双稳态折纸的应力极值几乎总是出现在折痕区尤其是折痕与面板交界的地方。折叠过程中折痕内半径一侧会出现明显压应力集中外半径一侧出现拉应力集中。PET的屈服强度大约在50到70MPa断裂延伸率很高。如果仿真结果显示折痕区的Mises应力远超这个范围说明单纯的弹性模拟已经不足以描述真实行为了需要考虑材料的塑性项或者修改折痕宽度、面板厚度等几何参数让变形更分散。另一个常被忽略的点是应力云的对称性。理想对称折叠时立方体的应力云图在相同结构位置上应呈现对称分布。如果左右两侧应力分布明显不对称而模型本身是对称的说明加载路径或者初始接触产生了非对称扰动。这时候要想清楚这是仿真自身的数值缺陷还是实物中真实存在的制造公差敏感性。后者在设计中反而更有价值。5. 排障与落地那些让折叠卡死、错乱、和实物对不上的坑最后这部分是真正动手以后才会暴露出来的问题。每个问题我都标注了典型的判断标志和排查方法你可以直接拿去做排除索引。5.1 隐式求解跑不动先想想是不是换Explicit折纸结构用Standard计算最容易看到的现象就是时间增量步一减再减最后停在一个极小的时间步长上原地颤抖。底层原因是系统刚度过零或者变负——也就是结构进入了失稳段。ABAQUS报错信息里高频出现的“负特征值”警告通常就是这个阶段产生的。遇到这种情况第一反应不要想着再怎么调网格、调阻尼先想想这一步的物理过程本来就包含失稳—重稳。失稳点附近隐式求解失败是“正确地反映了物理过程”不是你设错了参数。正确做法是切换到Explicit用显式时间积分走过这段失稳区。当然Explicit也有自己的翻车方式质量缩放设置过大导致惯性影响显著折叠过程变成慢动作撞击测试接触定义遗漏自接触导致面板互相穿透却发现不了。5.2 穿透、穿透、还是穿透折叠到中后期大量的面板会相互贴近。最常见的问题就是折叠闭合瞬态两片面板相互穿越看起来像科幻片里的穿模效果。这可能由几个原因导致接触定义里没有包含自接触或接触对只覆盖了外表面EXTERIOR而没有选中所有可能出现接触的面。壳单元的厚度方向相反导致ABAQUS不知道面板的“哪一面”是外面。这时候最简单的方法是进Mesh模块检查法向指向把面板外法向统一朝外。接触刚度设置过低罚接触太软导致允许的穿透量过大。对于这种薄壁结构接触的罚刚度可以适当提高只要不发散就行。排查时最容易漏掉的是“初始穿透”即模型初始状态就有单元互相嵌入。在Interactions模块里运行Interference Check先把初始穿透清零否则后续折叠过程会从错误的位置开始。5.3 从仿真“过拟合”回到实物可靠性的最后一公里仿真和实物对不上几乎每个做折纸机构的人都经历过。我的经验是偏差来源通常集中在三个地方一是折痕刚度的标定偏差。仿真里折痕区给的是一个折减后的弹性模量但实物折痕由压痕深度、板材厚度、材料韧性共同决定。更可靠的做法是做一个简单的单折痕实验用万能试验机给一条单独的折痕做弯折加载测出弯矩-角度曲线把这条曲线作为折痕区的非线性材料属性或者连接器属性反填进模型。这样一整套流程跑下来驱动力预测误差能压到15%以内。二是摩擦系数的“虚高”。仿真里默认一个0.15或0.3的恒定摩擦系数但实际折叠过程中接触面会发生相对滑移甚至微黏着摩擦系数并不是恒定常数。折叠力偏大时第一嫌疑就是摩擦力。对比实测值把摩擦系数往低修正通常能拉回差距。三是温度与时间效应。PET这类材料有粘弹性快折和慢折的力-位移曲线有明显差异。如果产品实际使用是慢慢压过去的0.03秒级的Explicit分析结果会高估峰值力。此时可以跑两个不同加载速度做对照看峰值力是否稳定如果差异过大就得考虑引入粘弹性本构来做速度敏感性评估。我自己的项目流程一般是先用统一参数的Explicit模型跑通折叠-展开循环→做能量和载荷判断→立刻做一个最小化的实物样机测单条折痕和整机折叠力→拿实测数据反过来校核折痕刚度与摩擦系数→用修正模型做最终的驱动机构和疲劳估算。这一圈下来仿真的价值才算真正兑现。双稳态折纸立方体的麻烦之处在于它的性能几乎完全由折痕区和接触行为支配这两个恰恰都是有限元里最考验经验的地方。但反过来说一旦你用能量曲线、载荷-位移曲线和应力分布这把尺子把整个折叠过程量出来这个结构的性能边界就非常清楚了——能垒多高、驱动力多大、哪里最脆弱、多久会疲劳全在数据里不在手感里。
阅读完成 · 觉得有帮助？
