1、7-1Workbench Mechanical结构非线性结构非线性第七章第七章诊断诊断Workbench Mechanical-Diagnostics7-2Training Manual章节综述章节综述 因为某些非线性结构分析求解会十分困难，因此理解怎样诊断不收敛因为某些非线性结构分析求解会十分困难，因此理解怎样诊断不收敛问题对获得答案十分关键。问题对获得答案十分关键。本章将介绍以下内容本章将介绍以下内容:求解器输出求解器输出 监视求解监视求解 Newton-Raphson 残差残差 本章介绍的功能需要本章介绍的功能需要ANSYS Structural licenses 及其以上的授权。及其以
2、上的授权。Workbench Mechanical-Diagnostics7-3Training ManualA.求解信息求解信息 在第二章中，已经介绍过求解信息分支在第二章中，已经介绍过求解信息分支 回顾通过回顾通过 Solution Information 分支可以观察到详细分支可以观察到详细 ANSYS 求解器输出和收求解器输出和收敛曲线敛曲线,例如，绘制力敛曲线例如，绘制力敛曲线.Text OutputGraphical OutputWorkbench Mechanical-Diagnostics7-4Training Manual 信息窗口信息窗口“Messages”位于位于Solu
3、tion Information Worksheet正下方，提供了正下方，提供了一些警告错误的列表一些警告错误的列表.右键任何信息右键任何信息:Go to Object(高亮和这些信息相关的项目树对象）高亮和这些信息相关的项目树对象）Show Full Message从分离的弹出窗口显示了所有的信息从分离的弹出窗口显示了所有的信息 Copy 拷贝信息拷贝信息(到分开的文档文件到分开的文档文件)Delete 从列表中删除信息从列表中删除信息.求解信息求解信息Workbench Mechanical-Diagnostics7-5Training ManualSolver Output 可提供关于求
4、解的详细的文字输出可提供关于求解的详细的文字输出.熟悉怎样阅读这个文熟悉怎样阅读这个文件将会非常有用件将会非常有用.1)Solver Output开始简单显示使用的开始简单显示使用的ANSYS license(这里是这里是 ANSYS Multiphysics)和版本号和版本号.求解信息求解信息Workbench Mechanical-Diagnostics7-6Training Manual2)Solver Output 记录了根据指定的单元记录了根据指定的单元(中节点中节点)选择激活的单元技术和相选择激活的单元技术和相关的材料，参照附录关的材料，参照附录B中更多关于单元技术的细节中更多关于
5、单元技术的细节.低阶单元的低阶单元的2D平面应力平面应力弹性材料或金属塑性弹性材料或金属塑性高阶高阶/低阶单元的完全不低阶单元的完全不可压缩超弹性可压缩超弹性高阶单元的弹性材料或高阶单元的弹性材料或金属塑性金属塑性低阶单元的低阶单元的2D 平面应变平面应变弹性应变或金属塑性弹性应变或金属塑性Default URISimplified Enhanced StrainEnhanced StrainB-Bar with Mixed u-P.求解信息求解信息Workbench Mechanical-Diagnostics7-7Training Manual3)向下翻页直到一系列星号出现，这时就会向下翻
6、页直到一系列星号出现，这时就会 看到求解器读取有限单元数据的信息。看到求解器读取有限单元数据的信息。这个列表非常有用，它不仅提供了模型这个列表非常有用，它不仅提供了模型 中零件数量的信息，同样也列出了接触中零件数量的信息，同样也列出了接触 区域的区域的ID号。号。需要指出的是虽然在需要指出的是虽然在Simulation中接触中接触 区域可以使用任何名字，区域可以使用任何名字，ANSYS 求解器只根求解器只根 据接触区域的唯一编号来进行处理据接触区域的唯一编号来进行处理(ID编号编号).找到接触区域的找到接触区域的ID编号对于程序调试非常有用。编号对于程序调试非常有用。例如，在上面的片段中，接触
7、区域例如，在上面的片段中，接触区域“Teeth 3”的的 接触编号为接触编号为9 和和 10.求解信息求解信息Workbench Mechanical-Diagnostics7-8Training Manual4)当方程求解开始以后，输出窗口如右图所示当方程求解开始以后，输出窗口如右图所示这里有用的信息包括求解器的选用这里有用的信息包括求解器的选用(是是“Program Chosen”或手动设定或手动设定)、大变形效应是否打开、是否考虑材料非大变形效应是否打开、是否考虑材料非线性线性(如果出现塑性如果出现塑性)和所使用的子步数和所使用的子步数.这一部分的求解器输出并不是很关键，这一部分的求解器
8、输出并不是很关键，但它给出了矩阵什么时候被求解及在但它给出了矩阵什么时候被求解及在Simulation中进行了哪些求解设置中进行了哪些求解设置.求解信息求解信息Workbench Mechanical-Diagnostics7-9Training Manual5)随后将会列出接触单元的细节随后将会列出接触单元的细节.这里，将列出有关接触单元的各种选项包括这里，将列出有关接触单元的各种选项包括 接触接触Normal Stiffness 和和Pinball 半径半径 任何注意或警告信息都会对观察有用任何注意或警告信息都会对观察有用.例如，在这里会给出初始渗透或缝隙大小例如，在这里会给出初始渗透或缝
9、隙大小 (以当前的长度单位以当前的长度单位)。从从Solver Section 2 回顾列出的接触区回顾列出的接触区域名称和域名称和ID号号.从那个例子可知，以上从那个例子可知，以上的接触实常数号的接触实常数号9 是是“Teeth 3”区域部区域部分分.求解信息求解信息Workbench Mechanical-Diagnostics7-10Training Manual6)随着非线性求解的进行，平衡迭代的信息将会显示在窗口的底部随着非线性求解的进行，平衡迭代的信息将会显示在窗口的底部(下图下图)注意，在每一次平衡迭代中，残余力注意，在每一次平衡迭代中，残余力(FORCE CONVERGENCE
10、 VALUE)必须小于准必须小于准则则 一个理想的平衡系统，残差或不平衡力应该为零。但由于机器的精度和实际的考虑，一个理想的平衡系统，残差或不平衡力应该为零。但由于机器的精度和实际的考虑，Simulation将确定一个小到几乎可以忽略所造成的误差的残差值。这个值就是准则将确定一个小到几乎可以忽略所造成的误差的残差值。这个值就是准则,而力收敛准则值必须小于准则而力收敛准则值必须小于准则 以使这个子步收敛。以使这个子步收敛。在下面分析中，经过在下面分析中，经过3次平衡迭代，残余力小于准则次平衡迭代，残余力小于准则,因此，求解收敛因此，求解收敛.通知信息通知信息(例如收敛或对分例如收敛或对分)子输出
11、窗口中将以子输出窗口中将以“”和和“标识标识)可以找到对分的原因可以找到对分的原因.求解信息求解信息Workbench Mechanical-Diagnostics7-14Training ManualB.监视求解监视求解 当求解非线性模型时，因为有可能需要多次迭代，观察非线性求解趋势将会当求解非线性模型时，因为有可能需要多次迭代，观察非线性求解趋势将会非常有用非常有用 如果求解看来会发生意外情况，用户可以终止分析来研究问题而不是等到求解如果求解看来会发生意外情况，用户可以终止分析来研究问题而不是等到求解结束，这就大大地节约了时间结束，这就大大地节约了时间.监视求解同样可以使用户更深入地了解系
12、统的响应监视求解同样可以使用户更深入地了解系统的响应 在在WB-Mechanical中有两种方法可以供用户监视求解中有两种方法可以供用户监视求解:Solution Information 分支提供方程求解器行为分支提供方程求解器行为 Results Tracker 提供求解过程中系统的响应提供求解过程中系统的响应Workbench Mechanical-Diagnostics7-15Training Manual 评估求解收敛行为的最好方法是观察力收敛曲线评估求解收敛行为的最好方法是观察力收敛曲线:在前一节中我们提到，在前一节中我们提到，Force Convergence 曲线可以在曲线可以在
13、Solution Information 分支中得到分支中得到 力收敛力收敛回顾在回顾在 Newton-Raphson 方法里，就是找方法里，就是找寻力平衡。如果不平衡力寻力平衡。如果不平衡力(残差残差)力收敛准则力收敛准则，那么这个子步就认为是收敛的。求解过程，那么这个子步就认为是收敛的。求解过程中，可以观察力收敛曲线跟踪过程。上面的中，可以观察力收敛曲线跟踪过程。上面的力收敛曲线说明二分出现了。这意味着施加力收敛曲线说明二分出现了。这意味着施加了更小的载荷步长。了更小的载荷步长。底部的底部的“TIME”曲线代表了总载荷的分数。曲线代表了总载荷的分数。虽然在静态分析中时间并无明确意义，它用虽
14、然在静态分析中时间并无明确意义，它用于作为于作为TIME=1.0最终求解的计数器。如果时最终求解的计数器。如果时间目前是间目前是0.2，那意味着，那意味着20%的载荷正在施加的载荷正在施加。Workbench Mechanical-Diagnostics7-16Training Manual 力收敛力收敛 通常需要注意的是求解过程中什么时候力收敛曲线开始通常需要注意的是求解过程中什么时候力收敛曲线开始“平稳平稳”这通常意味着这通常意味着(a)需要减小载荷增量或需要减小载荷增量或(b)接触法向刚度太高接触法向刚度太高.对于情况对于情况(a)Simulation 将自动对分求解将自动对分求解 对于
15、情况对于情况(b),二分或手动降低接触二分或手动降低接触 法向刚度法向刚度Workbench Mechanical-Diagnostics7-17Training Manual 结果追踪器结果追踪器 除了监视不平衡力，在除了监视不平衡力，在“Solution Information”分支分支中还可以添加结果追踪器中还可以添加结果追踪器(Result Tracker)Results Tracker 使用户能够在求解过程中监视某点点变使用户能够在求解过程中监视某点点变形或接触区域信息形或接触区域信息.对对“Results Tracker Deformation,”选择感兴趣的一选择感兴趣的一点，指
16、定需要监视的点，指定需要监视的 x,y,或或 z 方向的变形方向的变形.对对“Results Tracker Contact,”在下拉菜单中选择接触在下拉菜单中选择接触区域就可以观察跟踪的数量区域就可以观察跟踪的数量(例如，接触单元的数量）例如，接触单元的数量）Workbench Mechanical-Diagnostics7-18Training Manual 结果追踪器结果追踪器 当定义了当定义了Results Tracker 项和开始求解后项和开始求解后,用户就可以在用户就可以在 求解求解 过程过程 中中“跟踪跟踪”变形或接触结果。变形或接触结果。本例中，特定接触区域接触单元的数本例中，
17、特定接触区域接触单元的数目被监控。右图明显显示目被监控。右图明显显示 在在Time=1.4 1.7,接触单元数目从接触单元数目从0跳到了跳到了29。既然既然“Time”在非线性分析中是在非线性分析中是“位位置定格置定格”，这意味着在第一个载荷步，这意味着在第一个载荷步（Time=1.0）之后）之后,在载荷的在载荷的40%和和70%之间，接触确立了。之间，接触确立了。Workbench Mechanical-Diagnostics7-19Training Manual 结果追踪器结果追踪器 在求解过程中可以同时选择和观察多个结果在求解过程中可以同时选择和观察多个结果(同一类的同一类的)本例中本例
18、中,显示了显示了8个接触区的穿透量。可以观察最个接触区的穿透量。可以观察最大穿透量是否太大，如果是这样，求解会过早停止；大穿透量是否太大，如果是这样，求解会过早停止；既然用户知道了哪个接触区域有太多的穿透，那么可既然用户知道了哪个接触区域有太多的穿透，那么可以增加该区域的法向刚度。以增加该区域的法向刚度。从这出发，如果用户期望接触区分开到接触，在从这出发，如果用户期望接触区分开到接触，在这点检测过程没得到反映，用户也想知道求解完成后这点检测过程没得到反映，用户也想知道求解完成后什么结果能观测。什么结果能观测。本例中，接触区本例中，接触区“Teeth 1-4”在第二个载荷步进入在第二个载荷步进入
19、接触，接触，“Teeth 7”似乎一直到第似乎一直到第4个接触步才接触个接触步才接触 上上，这些都给用户展示了，这些都给用户展示了 哪些接触区正在接触，什么时哪些接触区正在接触，什么时候接触。候接触。Workbench Mechanical-Diagnostics7-20Training ManualC.Newton-Raphson 残差残差 到目前为止，已经讨论过了获得详细的求解器信息到目前为止，已经讨论过了获得详细的求解器信息(节节 A)和监视求解和监视求解(节节 B).当出现不收敛问题时，当出现不收敛问题时，Newton-Raphson 残差对找到可能有问题的区域非残差对找到可能有问题的
20、区域非常有用常有用 如果求解不收敛是由于力平衡原因，这将会反映在求解器输出和力收敛曲线上如果求解不收敛是由于力平衡原因，这将会反映在求解器输出和力收敛曲线上.力收敛力收敛(残差残差)将会大于力判据将会大于力判据.在不收敛情况下，在不收敛情况下，Newton-Raphson残差残差(如需要如需要)将会指出力残差高的区域。将会指出力残差高的区域。这通常有利于帮助确定阻止力平衡的部位这通常有利于帮助确定阻止力平衡的部位,通常是由于这些部位的力和支撑或接通常是由于这些部位的力和支撑或接触区域的设置触区域的设置.Workbench Mechanical-Diagnostics7-21Training M
21、anual Newton-Raphson 残差残差 在在“Solution Information”详细窗口中，输入平衡迭代的数目来获取详细窗口中，输入平衡迭代的数目来获取Newton-Raphson 残差残差.例如，输入例如，输入“4”,当求解中断或不收敛时，最后四次迭代的残当求解中断或不收敛时，最后四次迭代的残余力将会给出余力将会给出.Workbench Mechanical-Diagnostics7-22Training Manual Newton-Raphson 残差残差 当求解停止或不能收敛，就可以从当求解停止或不能收敛，就可以从“Solution Information”分支中获得
22、残差，分支中获得残差，如下图所示如下图所示.如果求解收敛失败或用户中止如果求解收敛失败或用户中止,可得到残可得到残差的数量差的数量.看残差，可以检查出现不平衡看残差，可以检查出现不平衡力的位置力的位置.有助于用户识别可能的有问有助于用户识别可能的有问题位置题位置,以致采取矫正措施以致采取矫正措施.Workbench Mechanical-Diagnostics7-23Training Manual Newton-Raphson 残差残差 在下例中，力残差十分高在下例中，力残差十分高(不平衡不平衡),Newton-Raphson Residuals 使使用户可以看到哪些区域对造成高的不平衡力有贡
23、献用户可以看到哪些区域对造成高的不平衡力有贡献 在这个例子中在这个例子中,高残余力的高残余力的 一圈一圈 是接触区域的一部分是接触区域的一部分,因此用户就因此用户就知道从哪里开始检查。知道从哪里开始检查。Workbench Mechanical-Diagnostics7-24Training ManualD.例子例子 在本节中，针对不同的情况将会作简要的介绍。观察求解器输出，监在本节中，针对不同的情况将会作简要的介绍。观察求解器输出，监视求解，及进行非线性诊断都会根据不同的例子进行讨论视求解，及进行非线性诊断都会根据不同的例子进行讨论 包含所有的不收敛的原因是不现实的，因此，在此只对用户会遇到
24、的一包含所有的不收敛的原因是不现实的，因此，在此只对用户会遇到的一些常见的问题进行讨论些常见的问题进行讨论.Workbench Mechanical-Diagnostics7-25Training Manual 刚体运动刚体运动 常见的错误信息常见的错误信息“internal solution magnitude limit was exceeded.”本质上说本质上说,这意味着出现了刚体运动这意味着出现了刚体运动.不足的约束，即使是有支撑或接触也有可能会使零件不足的约束，即使是有支撑或接触也有可能会使零件飞向飞向 空间空间 Solver Output 也会提供刚体运动方向也会提供刚体运动方向
25、Workbench Mechanical-Diagnostics7-26Training Manual 刚体运动刚体运动 为了检测哪个零件发生了刚体位移，用户可以进行一次模态分析来找到接近于为了检测哪个零件发生了刚体位移，用户可以进行一次模态分析来找到接近于零频率的模态。这有可能会十分耗时。零频率的模态。这有可能会十分耗时。另一种方法就是查看哪个接触对是初始打开的。回顾在求解器输出的第一部分另一种方法就是查看哪个接触对是初始打开的。回顾在求解器输出的第一部分可以找到和接触可以找到和接触ID对应的接触区域名。对应的接触区域名。本例中本例中,“contact offset block-2 To c
26、ontact offset bolt-1”接触接触区域区域(set#10 和和11)有一个初始间隙。可以检查模型看间隙是否有一个初始间隙。可以检查模型看间隙是否存在，以及因为由于缺乏初始接触的建立是否正在引起刚体位移存在，以及因为由于缺乏初始接触的建立是否正在引起刚体位移。Workbench Mechanical-Diagnostics7-27Training Manual 接触刚度接触刚度 有时，接触有时，接触Normal Stiffness 会太高而造成求解不收敛会太高而造成求解不收敛。对于这个螺纹扣件，在第一子步中只求解了对于这个螺纹扣件，在第一子步中只求解了20%的载荷，求解就发散了的
27、载荷，求解就发散了Workbench Mechanical-Diagnostics7-28Training Manual 接触刚度接触刚度 通过观察通过观察Newton-Raphson 残差残差,可以发现在接触区域内可能不会达到力平衡可以发现在接触区域内可能不会达到力平衡.网格看起来很粗糙网格看起来很粗糙 通过降低接触法向刚度和细化接触区域网格，就可以解决这个问题通过降低接触法向刚度和细化接触区域网格，就可以解决这个问题不收敛的不收敛的Newton-Raphson残差残差调整后求解收敛调整后求解收敛Workbench Mechanical-Diagnostics7-29Training Man
28、ual 模型设置模型设置 错误信息错误信息“solver engine was unable to converge on a solution for the nonlinear problem as constrained”可能会有多种因素造可能会有多种因素造成，虽然当这个信息出现时，重复检查模型设置通常是个好的办法成，虽然当这个信息出现时，重复检查模型设置通常是个好的办法.这个模型中这个模型中,出现了出现了Bolt Loads,所以有两个载荷所以有两个载荷步步.第一载荷步第一载荷步(预加载预加载)收敛收敛.然而第二载荷步甚然而第二载荷步甚至没有一个子步收敛至没有一个子步收敛.Workbe
29、nch Mechanical-Diagnostics7-30Training Manual 模型设置模型设置 通过检查通过检查Results Tracker 和和Solver Output,用户可以确定不收敛的原因用户可以确定不收敛的原因:在第一载荷步中，只有一个接触区域保持零件接触在第一载荷步中，只有一个接触区域保持零件接触,如左边的如左边的Results Tracker 所示所示.从右边的从右边的Solver Output 中可以发现相互接触的零件分开了，这就意味着很可能中可以发现相互接触的零件分开了，这就意味着很可能载荷的施加而造成零件脱离接触载荷的施加而造成零件脱离接触.检查项检查项:
30、是否其他接触区假设进入了初始接触是否其他接触区假设进入了初始接触？加载幅度不正确？忽略了其他约束？加载幅度不正确？忽略了其他约束吗？吗？Workbench Mechanical-Diagnostics7-31Training ManualE.求解模型的技巧求解模型的技巧 同上一节中提到的一样，综合利用下列步骤可帮助确定可能的不收敛原因同上一节中提到的一样，综合利用下列步骤可帮助确定可能的不收敛原因:观察观察Force Convergence 曲线了解求解行为曲线了解求解行为 通过通过Solver Output 检查接触信息检查接触信息(接触区域是否初始接触？初始缝隙接触区域是否初始接触？初始缝
31、隙/穿透是多穿透是多少？少？)和求解过程中断警告或错误信息和求解过程中断警告或错误信息.利用利用Results Tracker 监视接触信息或某点处的变形。模型响应是否和预期的一监视接触信息或某点处的变形。模型响应是否和预期的一致致?如果发生了不收敛，通过检查如果发生了不收敛，通过检查Newton-Raphson 余量来确定高余量区域，这些余量来确定高余量区域，这些区域就可能是有问题的区域。在这些有问题的区域是否施加了载荷或约束区域就可能是有问题的区域。在这些有问题的区域是否施加了载荷或约束,或这或这些区域是否是零件的接触区域？再次检查模型设置些区域是否是零件的接触区域？再次检查模型设置.Wo
32、rkbench Mechanical-Diagnostics7-32Training Manual 求解模型的技巧求解模型的技巧 如果以上步骤仍然不能提供足够的信息来确定问题所在，就可以使用以下的如果以上步骤仍然不能提供足够的信息来确定问题所在，就可以使用以下的方法方法:如果存在塑性和过度的单元畸变，可以试着先去掉塑性来运行模型以确定是否如果存在塑性和过度的单元畸变，可以试着先去掉塑性来运行模型以确定是否是由于材料模型而造成的这个问题是由于材料模型而造成的这个问题.如果问题可以确定为材料塑性所造成的，检查塑性定义。是否是完全塑性？如果问题可以确定为材料塑性所造成的，检查塑性定义。是否是完全塑性
33、？如果是，如果是，是否塑性铰链是否塑性铰链(mechanism)形式？形式？是否有可能是由于单元的刚度太小从而使单元变是否有可能是由于单元的刚度太小从而使单元变形过大形过大?对于接触问题对于接触问题,复制分支利用绑定接触来求解模型以确定是否是由于接触定义的复制分支利用绑定接触来求解模型以确定是否是由于接触定义的原因而造成这个问题原因而造成这个问题.如果是接触造成的问题而且力余量很高，改变所有接触为如果是接触造成的问题而且力余量很高，改变所有接触为Pure Penalty 同时输入较小同时输入较小的的 Normal Stiffness(0.01 到到 0.1).试着再次求解模型，有可能会出现过度
34、穿透试着再次求解模型，有可能会出现过度穿透,但如但如果得到结果，用户就可以对零件间接触的相互作用有一定的了解果得到结果，用户就可以对零件间接触的相互作用有一定的了解.Workbench Mechanical-Diagnostics7-33Training Manual 求解模型的技巧求解模型的技巧 如果对于接触使用了如果对于接触使用了Normal Lagrange Formulation,就检查振荡状态就检查振荡状态.如果遇到不收敛，振荡信息会在求解过程中和结束后打印出来如果遇到不收敛，振荡信息会在求解过程中和结束后打印出来 用户还可以在用户还可以在 Results Tracker 中查询中查
35、询Contact Chattering.如果振荡太多，就可以选用如果振荡太多，就可以选用Augmented Lagrange 方法或使用方法或使用Preprocessing Commands 分支来控制分支来控制(FTOLN)Normal Lagrange 接触的振荡接触的振荡.Workbench Mechanical-Diagnostics7-34Training Manual 总结总结 WB-Mechanical 提供了许多工具来帮助用户追踪非线性分析和诊断问题提供了许多工具来帮助用户追踪非线性分析和诊断问题.通常，可以从简单问题开始然后再逐渐增加问题的复杂性，这样确定问题的原通常，可以从
36、简单问题开始然后再逐渐增加问题的复杂性，这样确定问题的原因更加有效。在初次分析就加入许多复杂性会在随后的分析中浪费许多时间因更加有效。在初次分析就加入许多复杂性会在随后的分析中浪费许多时间.不要随机改变设置。首先使用推荐的设置，然后如果有明确的原因才能改变接不要随机改变设置。首先使用推荐的设置，然后如果有明确的原因才能改变接触或求解器设置就象在触或求解器设置就象在Solver Output,Results Tracker,或或 Newton-Raphson Residuals所阐释的一样所阐释的一样.虽然接触区域是自动建立的，总是所有的接触虽然接触区域是自动建立的，总是所有的接触/目标表面和接
37、触设置来确保接触目标表面和接触设置来确保接触区域如所希望的一样定义区域如所希望的一样定义.在在 Solver Output 中观察详细的接触输出，检查接触中观察详细的接触输出，检查接触区域的初始接触状态，极可能存在初始穿透或缝隙，如果有，那么可找出它们区域的初始接触状态，极可能存在初始穿透或缝隙，如果有，那么可找出它们的大小的大小.Workbench Mechanical-Diagnostics7-35Training Manual 请参考作业中的补充说明请参考作业中的补充说明:作业作业7A 接触诊断工具接触诊断工具F.作业作业 非线性诊断非线性诊断Workbench Mechanical-Diagnostics7-36Training Manual 请参考作业中的补充说明请参考作业中的补充说明:作业作业7B 奇奇异性异性.作业作业 非线性诊断非线性诊断
