如果你拿一根铁丝，使劲掰一下，它可能只是弯了（塑性变形）；但如果你来回反复折 50 次，它就会"嘎嘣"一声断掉。这就是疲劳（Fatigue）。
疲劳分析的物理本质：
疲劳是材料在循环载荷作用下，即使应力水平低于屈服强度，也会逐渐产生损伤并最终断裂的现象。与静载破坏不同，疲劳破坏具有以下特征：
在仿真界，模拟这种反复折磨最痛苦的地方在于：前 1000 次循环看起来都差不多。如果按部就班地算，99% 的算力都浪费在重复的"过日子"上了。在 Abaqus 的 API 宇宙里，DirectCyclicStep 就像是一个能够看穿轮回的智者，它不关心你每一次是怎么折的，它直接去寻找那个稳定后的状态（Stabilized Cycle）。
根据官方脚本指南，DirectCyclicStep 是一种专门用于计算受周期性载荷作用下结构响应的高级分析步。
在创建一个 DirectCyclicStep 时，最重要的任务是确定结构是否进入了"安定状态"。
当材料在反复加载下，其应力 - 应变曲线不再移动，而是围成一个固定的圈时，我们就说它达到了"稳定循环"。
这就像是你去健身房练哑铃。刚开始肌肉会酸痛变形，但练了几天后，你的动作和肌肉反应会固定在一个节奏上。DirectCyclic 脚本直接跳过最初的酸痛期，一键传送你到那个"稳定健身"的状态。
安定状态 (Shakedown) 的三种情况：
1. 弹性安定 (Elastic Shakedown)：
经过若干循环后，结构完全弹性响应
无塑性应变累积
寿命理论上无限
2. 塑性安定 (Plastic Shakedown)：
达到稳定的塑性循环
每个循环的塑性应变相同
存在有限的疲劳寿命
3. 棘轮效应 (Ratchetting)：
2. 时间的"乐谱"：傅里叶级数（Fourier Series）API 参数中有一个关键项 `numTermFourier`（傅里叶项数）。
DirectCyclic 不用时间步增量，而是把载荷和响应拆解成一堆正弦波和余弦波的叠加。
如果普通分析是"拍视频"，那 DirectCyclic 就是在"写乐谱"。它把复杂的来回运动拆解成高低不同的音符。你给出的傅里叶项数越多，这首"轮回舞曲"就还原得越真实。
傅里叶级数的数学原理：
任何周期函数都可以表示为傅里叶级数：
在 DirectCyclic 中：
载荷被分解为傅里叶级数
响应也用傅里叶级数表示
求解的是傅里叶系数，而非时域值
傅里叶项数的选择：
正弦波，1-3项，基波即可描述
方波，10-20项，需要高次谐波
锯齿波，15-30项，需要更多谐波
复杂波形，30-50项，捕捉细节
含突变的波形，50+项，需要高频成分
配合 Damage 模型，DirectCyclic 可以直接预测裂纹什么时候萌生。
它是代码界的"预言家"。通过计算稳定循环中的能量耗散，它能告诉你：这根铁丝在第 15432 次折弯时，心脏（材料内部）将彻底破碎。
低周疲劳损伤模型：
Abaqus 提供多种损伤演化模型：
1. Coffin-Manson 模型：
2. 能量法 (Energy-based)：
3. 连续损伤力学 (CDM)：
在 Python 中调用这个强大的"时间机器"，你的脚本需要这样写：
from abaqus import *from abaqusConstants import *# 1. 创建直接循环分析步# timePeriod: 每一个循环周期的时间（比如 1 秒折一次）# maxNumIterations: 寻找稳定状态的最大尝试次数# numTermFourier: 用 50 个傅里叶项来模拟这个"轮回"mdb.models['Model-1'].DirectCyclicStep( name='Samsara_Analysis', previous='Initial', timePeriod=1.0, maxNumIterations=200, numTermFourier=50, initialPredictor=ON # 开启初始预测，加速收敛)
别忘了配套定义 Amplitude（幅值曲线）。 因为是循环分析，Abaqus 需要知道你的载荷是怎么"荡秋千"的。
如果你的模型在 200 次迭代后还没稳定，记得检查是不是产生了严重的塑性累积（Ratchetting）。
完整的直接循环分析流程：
from abaqus import mdbfrom abaqusConstants import *def setup_complete_direct_cyclic(model_name): if model_name not in mdb.models: print(f"错误：模型 {model_name} 不存在！") return model = mdb.models[model_name] # --- 步骤 1：周期性幅值曲线 --- # 定义一个完整的周期 (Time=1.0) amp_name = 'Cyclic-Amplitude' model.TabularAmplitude( name=amp_name, timeSpan=STEP, # 注意：直接循环通常配合 STEP 时间跨度 data=((0.0, 0.0), (0.25, 1.0), (0.5, 0.0), (0.75, -1.0), (1.0, 0.0)) ) # --- 步骤 2：直接循环分析步 --- # initialPredictor=ON 有助于加快首次循环的收敛 step_name = 'Direct-Cyclic-Step' model.DirectCyclicStep( name=step_name, previous='Initial', timePeriod=1.0, maxNumIterations=100, # 允许的最大迭代次数 numTermFourier=20, # 傅里叶级数项数 initialPredictor=ON, fatigue=OFF # 如果要做低周疲劳，这里需设为 ON ) # --- 步骤 3：循环载荷 --- # 确保 Load-Point 这个 Set 真的存在 try: region = model.rootAssembly.sets['Load-Point'] model.ConcentratedForce( name='Cyclic-Load', createStepName=step_name, region=region, cf2=1000.0, amplitude=amp_name ) except KeyError: print("警告：未找到 'Load-Point' 集，载荷未施加！") # --- 步骤 4：精细化输出请求 --- # 默认场输出通常已经存在，建议用 setValues 更新 if 'F-Output-1' in model.fieldOutputRequests.keys(): model.fieldOutputRequests['F-Output-1'].setValues( variables=('S', 'E', 'PE', 'PEEQ', 'U', 'RF'), frequency=1 # 直接循环中，1 表示输出最后一个稳定循环 ) # 历史输出：监控收敛是分析的关键 model.HistoryOutputRequest( name='Convergence-Check', createStepName=step_name, variables=('ALLIE', 'ALLPD'), # 监控内能和塑性耗散能量 frequency=EVERY_INCREMENT ) print(f"模型 {model_name} 已完成直接循环分析设置！") print("科趣提示：提交作业前，请确保单元类型支持温度/力耦合（如适用）。")# 使用示例setup_complete_direct_cyclic('Cyclic-Model')
在 Abaqus 脚本建模中，DirectCyclic 是对"低效重复"的一次降维打击。
最佳实践建议：
1. 从简单开始：先用少量傅里叶项测试
2. 逐步增加：根据收敛情况调整项数
3. 检查安定：确认是否达到稳定循环
4. 验证载荷：确保幅值曲线正确
5. 材料数据：使用准确的循环应力-应变数据
很多时候，事物的本质不在于它经历了多少次磨难，而在于它在磨难中最终呈现出的那个"稳定姿态"。作为工程师，学会用代码跳过冗长的过程，直击物理规律的终点，这才是最高级的职业优雅。
👉互动话题：你见过最离奇的疲劳破坏是什么？是刚买的耳机线断了，还是飞机的发动机裂了？在评论区聊聊你对"折磨"的理解，咱们一起用代码消灭它们！
