简介：这份文档面向结构工程与有限元分析人员，聚焦ANSYS中预应力混凝土（有粘结）的建模方法，帮助读者在等效荷载法与实体力筋法之间做出合理选择。内容围绕两种方法的原理与优缺点展开：等效荷载法建模简单、整体效应易求，但难以反映力筋分布方向与应力损失；实体力筋法则用solid模拟混凝土、link模拟力筋，通过降温法或初应变法施加预应力，精度更高。文档还介绍了体分割法与独立建模耦合法的实现思路，并给出预应力简支梁弹性分析的APDL示例，涵盖材料属性定义、单元类型设置、力筋线建立、节点耦合及边界条件施加等关键环节。资源包为1个doc文件，约59KB，篇幅精炼、要点集中，适合作为预应力分析入门与建模思路参考。目前已有347人学习下载。
1. 预应力分析为什么总在“加完预紧力”这一步翻车
做结构仿真的人迟早会撞上预应力这个坎。螺栓要预紧、混凝土要张拉、轴承要压装、弹簧要预压缩——只要结构在承受外载之前内部已经存在一套自平衡的应力场，你就绕不开预应力分析。很多人第一次在 ANSYS 里做预应力，流程看着都跑通了，结果一看变形云图，要么位移大得离谱，要么应力分布完全不符合直觉，回头查边界条件又找不出毛病。问题往往不在求解器，而在“预应力到底怎么加进去”这一步被想当然了。
这篇笔记面向的是已经会用 ANSYS 做静力分析、但一碰预应力就心里没底的工程师，也适合需要把预紧结构仿真结果写进报告、被评审追问过“你这个预紧力怎么等效的”的人。我会按“预应力在求解器里到底以什么形式存在 → 怎么建模和加载 → 怎么验证结果可信 → 哪些参数一改就翻车”的顺序讲，中间给可直接抄的命令流和参数设置。读完你应该能独立搭出一个能收敛、结果能自圆其说的预应力模型，并且知道哪些地方是玄学、哪些地方是硬约束。
2. 预应力在 ANSYS 里到底是怎么被求解器“记住”的
2.1 预应力不是一种载荷，而是一种初始状态
很多人把预应力当成一个“力”去理解，于是在加载步里到处找“施加预应力”的按钮，找不到就慌了。实际上在 ANSYS 里，预应力从来不是一个独立的载荷类型，它是结构在某个参考构型下已经存在的应力状态。求解器要做的，是把这个初始应力场作为后续分析的起点。
这就带来一个关键区分：预应力分析分两大类。第一类是“预应力模态”或“预应力谐响应”，也就是先做一次静力分析把应力场算出来，再把刚度矩阵里由应力引起的几何刚度项带上，去做后续的模态或频响。第二类是“带预应力的静力分析”，比如螺栓预紧后继续加载，这时候预应力本身就是静力分析的一部分，只是分步施加。
理解这一点，后面所有操作就顺了：你要么是在“制造一个初始应力场”，要么是在“把已有应力场传递到下一步”。ANSYS 里对应的手段无非三种——初应力法、预紧单元法、分步加载法。选哪种，取决于你的结构是螺栓、是张拉索、还是压装配合。
2.2 三种常见实现路径的选型对比
| 方法 | 适用场景 | 关键命令/单元 | 收敛性 | 结果可信度 |
|---|---|---|---|---|
| 初应力法 | 已知初始应力场，如混凝土张拉 | INISTATE | 好 | 依赖输入应力准确性 |
| 预紧单元法 | 螺栓、拉杆预紧 | PRETS179 | 中等 | 高，能模拟预紧过程 |
| 分步加载法 | 压装、过盈配合 | 多载荷步 + 接触 | 较差 | 高，但难收敛 |
初应力法最直接，你把算好的或测得的应力值通过 INISTATE 写进单元，求解器直接把它当初始状态。优点是简单，缺点是你要先知道应力是多少，而且它不模拟“怎么加上去的”，只模拟“加完之后”。预紧单元法专门为螺栓类结构设计，PRETS179 单元可以施加一个预紧力，求解器会自动把它转换成螺栓内部的拉力，还能在后续载荷步里锁定这个预紧位移。分步加载法最贴近物理，但要把过盈量、接触、摩擦都设对，收敛难度直线上升。
我一般会这样选：螺栓连接优先用预紧单元，混凝土或预应力索用初应力法，压装配合如果接触定义得清楚就用分步加载，否则退而求其次用初应力等效。
2.3 一个最小可复现的预紧螺栓模型
下面这段命令流是一个最简化的预紧螺栓模型，两块板夹一根螺栓，用 PRETS179 施加预紧力。你可以直接改尺寸和力值跑通。
! 简化预紧螺栓模型 - 命令流 /PREP7 ET,1,SOLID185 ! 板单元 ET,2,SOLID185 ! 螺栓单元 ET,3,PRETS179 ! 预紧单元 MP,EX,1,2.1E11 ! 钢弹性模量 MP,PRXY,1,0.3 ! 建两块板，尺寸 50x50x10，中间留螺栓孔 BLOCK,0,50,0,50,0,10 BLOCK,0,50,0,50,10,20 VGLUE,ALL ! 建螺栓杆，直径 10，贯穿两块板 CYLIND,0,5,0,20,0,360 ! 网格划分（此处省略尺寸控制，实际需细化孔边） VSEL,S,LOC,Z,0,20 VATT,1,1,1 VMESH,ALL ! 定义预紧截面：在螺栓中间切一刀 WPOFFS,0,0,10 VSBW,ALL ! 在切面上生成预紧单元 VSEL,S,LOC,Z,10 ASEL,R,LOC,Z,10 NSLA,S,1 ESEL,S,MAT,,1 ! 此处需根据实际节点生成 PRETS179，简化示意 ! 实际操作用 EINTF 或手动定义预紧节点这段代码的重点不在网格，而在最后几步：预紧单元必须建在螺栓的横截面上，而且这个截面要垂直于螺栓轴线。PRETS179 单元由一个预紧节点和一组被预紧的节点组成，求解器通过这个单元施加力，并自动计算需要多少位移才能达到这个力。参数上，预紧力大小通过 SLOAD 命令施加，方向沿单元坐标系。
逻辑说明：先建几何、划网格，然后在螺栓中间切出一个截面，在这个截面上定义预紧单元。预紧单元不是实体单元，它是一层“虚拟”的连接，把截面两侧的节点耦合起来，同时允许施加一个拉力。求解时，ANSYS 会先施加预紧力，算出位移，然后在后续载荷步里保持这个位移或力。
参数说明：预紧力大小根据工程要求定，比如 M10 螺栓常见预紧力在 20~40 kN。弹性模量和泊松比按材料填。接触如果考虑，需要在板与板、螺栓头与板之间定义 CONTACT 单元，摩擦系数一般取 0.1~0.2。
2.4 初应力法的 INISTATE 怎么写才不出错
如果你做的是混凝土张拉或预应力索，初应力法更顺手。核心命令是 INISTATE。
! 初应力法示例 - 给指定单元组施加初始应力 ESEL,S,MAT,,2 ! 选择预应力索单元 INISTATE,SET,CSYS,0 ! 使用全局坐标系 INISTATE,SET,DTYP,STRESS INISTATE,DEFINE,,,,,1.2E8,0,0,0,0,0 ! 施加 X 向 120 MPa 初应力这里最容易翻车的是坐标系。INISTATE 默认用全局坐标，但你的索方向可能不在全局 X 向。如果方向错了，应力就加到了错误的方向上，结果自然离谱。另一个坑是初应力必须自平衡，如果你只给索加应力而不给混凝土加对应的压应力，第一步求解就会因为不平衡力过大而发散。常见做法是给索和混凝土分别施加大小相等、方向相反的初应力，或者用 ISFILE 读入一个自平衡的应力场。
3. 从建模到求解：预应力分析的完整操作链
3.1 几何处理与网格划分的硬约束
预应力分析的几何处理比普通静力分析苛刻。螺栓孔周围、张拉端锚固区、压装配合面，这些地方的网格质量直接决定收敛性。我一般会这样做：孔边至少两层单元，配合面附近网格尺寸不超过过盈量的三分之一。如果用过盈配合，网格太粗会导致接触穿透，求解器要么不收敛，要么给出错误的接触压力。
另一个硬约束是预紧截面的位置。预紧单元必须建在螺栓的直线段上，不能建在螺纹或倒角处。截面要垂直于轴线，否则预紧力方向会偏。实际操作中，我会在螺栓中间切一刀，确保截面干净、节点整齐。
3.2 载荷步设置：先预紧，再锁定，后加载
预应力分析的载荷步不能一步到位。正确顺序是：第一步施加预紧力，第二步锁定预紧位移，第三步施加外载。用 PRETS179 时，锁定通过 LOCK 命令实现。
/SOLU ANTYPE,STATIC ! 第一步：施加预紧力 SLOAD,1,PL01,LOCK,FORC,30000 ! 预紧力 30 kN LSWRITE,1 ! 第二步：锁定预紧位移 SLOAD,1,PL01,LOCK,ON LSWRITE,2 ! 第三步：施加外载 F,节点号,FX,5000 LSWRITE,3 LSSOLVE,1,3,1逻辑说明：第一步 SLOAD 施加力，求解器算出达到这个力所需的位移。第二步把预紧状态锁定，意味着后续载荷不会改变预紧力的大小，但预紧力会作为内力参与平衡。第三步加外载，求解器在预紧基础上继续求解。
参数说明：LOCK,FORC 后面的数值是预紧力，单位跟你的模型单位制一致。如果模型用毫米-牛顿-秒，力就是牛顿。LOCK,ON 表示锁定当前位移。外载大小按实际工况给。
3.3 接触与摩擦：预紧分析里最容易忽略的变量
螺栓头与板、板与板之间的接触，如果不定义，预紧力就传不出去，结果会显示螺栓被拉长但板没被压紧。接触定义要注意三点：接触刚度、摩擦系数、接触算法。
接触刚度 FKN 一般取 0.1~1.0，太大收敛难，太小穿透多。摩擦系数按实际表面处理取，机加工面 0.1~0.15，喷砂面 0.3~0.5。算法用增广拉格朗日法比罚函数法更稳，但计算量大。我一般先用罚函数跑，不收敛再换增广拉格朗日。
提示：接触定义后，一定要检查接触状态。求解前用 CNCHECK 命令查看初始接触是否有间隙或穿透。初始穿透会导致第一步就发散。
3.4 求解器选择与收敛调试
预应力分析建议用稀疏矩阵求解器，默认就行。如果模型大，可以用 PCG 求解器省内存。收敛调试的核心是看力残差和位移残差。如果残差震荡不降，先查接触，再查网格，最后查预紧力是否过大导致局部屈服。
常见收敛问题及对策：
| 现象 | 可能原因 | 对策 |
|---|---|---|
| 第一步就发散 | 初应力不自平衡 | 检查 INISTATE 方向，补平衡应力 |
| 预紧步收敛慢 | 接触刚度太大 | 降低 FKN 到 0.1 |
| 锁定步报错 | 预紧单元方向反了 | 检查 PRETS179 节点顺序 |
| 外载步不收敛 | 材料屈服或接触滑移 | 加子步，开自动时间步 |
4. 预应力分析避坑：5 个血泪教训
4.1 现象：预紧力加完，螺栓应力只有理论值的一半
原因：预紧单元建在了螺栓的螺纹段或倒角处，有效承载面积比光杆段小，求解器按实际截面算应力，自然偏低。解决：把预紧截面移到光杆段，或者手动修正预紧力值，按光杆截面积换算。
4.2 现象：锁定后加外载，预紧力自己变了
原因：LOCK,ON 锁定的是位移，不是力。如果外载导致螺栓进一步伸长，预紧力会增大；如果外载导致板分离，预紧力会减小。这是物理事实，不是 bug。解决：如果你希望预紧力严格不变，用 LOCK,FORC 而不是 LOCK,ON，但这样不符合实际。一般工程上接受预紧力在小范围内波动。
4.3 现象：初应力法加完，第一步求解就报“不平衡力过大”
原因：初应力场不自平衡。你给索加了拉应力，但混凝土没有对应的压应力，整个系统合力不为零。解决：要么给混凝土也加初应力，要么用 ISFILE 读入一个自平衡的应力场，要么先做一次无外载的静力分析让求解器自己找平衡。
4.4 现象：过盈配合分析，接触压力远小于理论值
原因：网格太粗，接触面穿透严重，求解器把过盈量“吃掉”了。解决：细化接触面网格，至少三层单元覆盖过盈区，同时降低接触刚度，让接触更容易建立。
4.5 现象：预应力模态分析，频率比预期低很多
原因：预应力模态里，几何刚度项可能被忽略，或者预紧力方向搞反了，导致刚度被削弱而不是增强。解决：检查 PSTART 命令是否打开，确认预紧力是拉力而不是压力。如果是压力，结构可能已经失稳，频率自然会降。
5. 怎么验证你的预应力结果不是“看起来对”
5.1 用反力平衡做第一道校验
求解完成后，第一件事是查反力。预紧力是内力，在预紧步，支座反力应该为零（如果模型没有其他外载）。如果反力不为零，说明预紧力没有自平衡，或者边界条件有问题。在锁定步和外载步，反力应该等于外载。这个校验能筛掉大部分低级错误。
5.2 用解析解对比螺栓预紧
对于单螺栓连接，理论预紧力 F、螺栓应力 σ = F/A，A 是光杆截面积。仿真结果应该在这个值附近，误差不超过 10%。如果差太多，查预紧截面位置和网格。另外，板间接触压力应该均匀分布，如果边缘压力远大于中心，说明接触刚度或网格有问题。
5.3 用能量法判断收敛质量
ANSYS 输出里有应变能和动能。静力分析中动能应该接近零。如果动能不为零，说明有刚体位移或接触滑移。应变能应该随载荷步平稳增加，如果突然跳变，说明有单元畸变或接触状态突变。我一般会看 ETABLE 里的 SENE 和 KENE，对比一下量级。
5.4 一个我常犯的错：忘了关掉大变形
预应力分析里，如果预紧力导致结构刚度变化明显，必须打开大变形（NLGEOM,ON）。我早期做压装分析时忘了开，结果过盈量被线性近似，接触压力算出来偏小 30%。打开大变形后，结果和理论值对上了。这个开关在 /SOLU 里设，NLGEOM,ON，然后重新求解。
5.5 最后一步：把预紧力写成报告里的可追溯参数
预应力分析的报告不能只写“施加了预紧力”，要写清楚预紧力大小、施加方式、锁定方式、接触参数。我习惯在报告里附一个参数表，包括螺栓规格、预紧力、摩擦系数、接触刚度、网格尺寸。这样评审问起来，你能直接指给他看。这个习惯帮我省了很多返工。
希望帮到你。
本文还有配套的精品资源，点击获取
)
