搜索
7-1Workbench–Mechanical结构非线性第七章诊断WorkbenchMechanical-Diagnostics7-2TrainingManual章节综述•因为某些非线性结构分析求解会十分困难，因此理解怎样诊断不收敛问题对获得答案十分关键。•本章将介绍以下内容:–求解器输出–监视求解–Newton-Raphson残差•本章介绍的功能需要ANSYSStructurallicenses及其以上的授权。WorkbenchMechanical-Diagnostics7-3TrainingManualA.求解信息•在第二章中，已经介绍过求解信息分支–回顾通过SolutionInformation分支可以观察到详细ANSYS求解器输出和收敛曲线,例如，绘制力敛曲线.TextOutputGraphicalOutputWorkbenchMechanical-Diagnostics7-4TrainingManual•信息窗口“Messages”位于SolutionInformationWorksheet正下方，提供了一些警告错误的列表.–右键任何信息:•GotoObject(高亮和这些信息相关的项目树对象）•ShowFullMessage从分离的弹出窗口显示了所有的信息•Copy拷贝信息(到分开的文档文件)•Delete从列表中删除信息...求解信息WorkbenchMechanical-Diagnostics7-5TrainingManual•SolverOutput可提供关于求解的详细的文字输出.熟悉怎样阅读这个文件将会非常有用.1)SolverOutput开始简单显示使用的ANSYSlicense(这里是ANSYSMultiphysics)和版本号....求解信息WorkbenchMechanical-Diagnostics7-6TrainingManual2)SolverOutput记录了根据指定的单元(中节点)选择激活的单元技术和相关的材料，参照附录B中更多关于单元技术的细节.低阶单元的2D平面应力弹性材料或金属塑性高阶/低阶单元的完全不可压缩超弹性高阶单元的弹性材料或金属塑性低阶单元的2D平面应变弹性应变或金属塑性DefaultURISimplifiedEnhancedStrainEnhancedStrainB-BarwithMixedu-P...求解信息WorkbenchMechanical-Diagnostics7-7TrainingManual3)向下翻页直到一系列星号出现，这时就会看到求解器读取有限单元数据的信息。这个列表非常有用，它不仅提供了模型中零件数量的信息，同样也列出了接触区域的ID号。需要指出的是虽然在Simulation中接触区域可以使用任何名字，ANSYS求解器只根据接触区域的唯一编号来进行处理(ID编号).找到接触区域的ID编号对于程序调试非常有用。例如，在上面的片段中，接触区域“Teeth3”的接触编号为9和10....求解信息WorkbenchMechanical-Diagnostics7-8TrainingManual4)当方程求解开始以后，输出窗口如右图所示这里有用的信息包括求解器的选用(是“ProgramChosen”或手动设定)、大变形效应是否打开、是否考虑材料非线性(如果出现塑性)和所使用的子步数.这一部分的求解器输出并不是很关键，但它给出了矩阵什么时候被求解及在Simulation中进行了哪些求解设置....求解信息WorkbenchMechanical-Diagnostics7-9TrainingManual5)随后将会列出接触单元的细节.这里，将列出有关接触单元的各种选项包括接触NormalStiffness和Pinball半径任何注意或警告信息都会对观察有用.例如，在这里会给出初始渗透或缝隙大小(以当前的长度单位)。从SolverSection2回顾列出的接触区域名称和ID号.从那个例子可知，以上的接触实常数号9是“Teeth3”区域部分....求解信息WorkbenchMechanical-Diagnostics7-10TrainingManual6)随着非线性求解的进行，平衡迭代的信息将会显示在窗口的底部(下图)•注意，在每一次平衡迭代中，残余力(FORCECONVERGENCEVALUE)必须小于准则•一个理想的平衡系统，残差或不平衡力应该为零。但由于机器的精度和实际的考虑，Simulation将确定一个小到几乎可以忽略所造成的误差的残差值。这个值就是准则,而力收敛准则值必须小于准则以使这个子步收敛。•在下面分析中，经过3次平衡迭代，残余力小于准则,因此，求解收敛.•通知信息(例如收敛或对分)子输出窗口中将以“”和““标识....求解信息WorkbenchMechanical-Diagnostics7-11TrainingManual警告和错误信息同样会显示在输出窗口中•当接触状态突变时，就会出现一个警告信息指出接触单元剧烈地进出‘pinball区域’.这可能是由於载荷太大造成零件间剧烈地相互滑动或分离。如果需要Simulation会自动二分求解。•单元畸变通常是由于载荷过大或过度约束所造成的严重问题。对分载荷会自动进行，但有时会需要一些修正措施来解决这个问题。...求解信息WorkbenchMechanical-Diagnostics7-12TrainingManual7)最后,当求解完成以后，在求解输出窗口的最后会提供一些本次分析的统计数据。求解时间与计算其它单元的时间的百分比以及整个方程的求解时间.在最后会显示整个Elapsedtime.•对于单个处理器，如果walltime远远大于CPU时间(例如2x),那就表示慢的磁盘I/O大大影响了整个计算时间.•对于双处理器，CPU时间将不会精确，因为它是两个处理器所用时间之和....求解信息WorkbenchMechanical-Diagnostics7-13TrainingManual•熟悉求解器输出的内容及结构对理解有限元求解十分有用–用户可以通过在求解过程中观察求解器输出来了解当前平衡迭代的详细信息。–在求解后观察求解器输出可以确定不收敛的原因(如果求解失败),获得求解器性能信息,和找到接触单元的设置信息–详细的警告或错误信息(由***标识)也将会显示在求解器输出中。同样在求解过程中观察子步状态(由标识)可以找到对分的原因....求解信息WorkbenchMechanical-Diagnostics7-14TrainingManualB.监视求解•当求解非线性模型时，因为有可能需要多次迭代，观察非线性求解趋势将会非常有用–如果求解看来会发生意外情况，用户可以终止分析来研究问题而不是等到求解结束，这就大大地节约了时间.–监视求解同样可以使用户更深入地了解系统的响应•在WB-Mechanical中有两种方法可以供用户监视求解:–SolutionInformation分支提供方程求解器行为–ResultsTracker提供求解过程中系统的响应WorkbenchMechanical-Diagnostics7-15TrainingManual•评估求解收敛行为的最好方法是观察力收敛曲线:–在前一节中我们提到，ForceConvergence曲线可以在SolutionInformation分支中得到…力收敛回顾在Newton-Raphson方法里，就是找寻力平衡。如果不平衡力(残差)力收敛准则，那么这个子步就认为是收敛的。求解过程中，可以观察力收敛曲线跟踪过程。上面的力收敛曲线说明二分出现了。这意味着施加了更小的载荷步长。底部的“TIME”曲线代表了总载荷的分数。虽然在静态分析中时间并无明确意义，它用于作为TIME=1.0最终求解的计数器。如果时间目前是0.2，那意味着20%的载荷正在施加。WorkbenchMechanical-Diagnostics7-16TrainingManual…力收敛•通常需要注意的是求解过程中什么时候力收敛曲线开始“平稳”–这通常意味着(a)需要减小载荷增量或(b)接触法向刚度太高.–对于情况(a)Simulation将自动对分求解–对于情况(b),二分或手动降低接触法向刚度WorkbenchMechanical-Diagnostics7-17TrainingManual…结果追踪器•除了监视不平衡力，在“SolutionInformation”分支中还可以添加结果追踪器(ResultTracker)–ResultsTracker使用户能够在求解过程中监视某点点变形或接触区域信息.–对“ResultsTrackerDeformation,”选择感兴趣的一点，指定需要监视的x,y,或z方向的变形.–对“ResultsTrackerContact,”在下拉菜单中选择接触区域就可以观察跟踪的数量(例如，接触单元的数量）WorkbenchMechanical-Diagnostics7-18TrainingManual…结果追踪器•当定义了ResultsTracker项和开始求解后,用户就可以在求解过程中“跟踪”变形或接触结果。本例中，特定接触区域接触单元的数目被监控。右图明显显示在Time=1.4~1.7,接触单元数目从0跳到了29。既然“Time”在非线性分析中是“位置定格”，这意味着在第一个载荷步（Time=1.0）之后,在载荷的40%和70%之间，接触确立了。WorkbenchMechanical-Diagnostics7-19TrainingManual…结果追踪器•在求解过程中可以同时选择和观察多个结果(同一类的)本例中,显示了8个接触区的穿透量。可以观察最大穿透量是否太大，如果是这样，求解会过早停止；既然用户知道了哪个接触区域有太多的穿透，那么可以增加该区域的法向刚度。从这出发，如果用户期望接触区分开到接触，在这点检测过程没得到反映，用户也想知道求解完成后什么结果能观测。本例中，接触区“Teeth1-4”在第二个载荷步进入接触，“Teeth7”似乎一直到第4个接触步才接触上，这些都给用户展示了哪些接触区正在接触，什么时候接触。WorkbenchMechanical-Diagnostics7-20TrainingManualC.Newton-Raphson残差•到目前为止，已经讨论过了获得详细的求解器信息(节A)和监视求解(节B)..•当出现不收敛问题时，Newton-Raphson残差对找到可能有问题的区域非常有用–如果求解不收敛是由于力平衡原因，这将会反映在求解器输出和力收敛曲线上.力收敛(残差)将会大于力判据.–在不收敛情况下，Newton-Raphson残差(如需要)将会指出力残差高的区域。这通常有利于帮助确定阻止力平衡的部位,通常是由于这些部位的力和支撑或接触区域的设置.WorkbenchMechanical-Diagnostics7-21TrainingManual…Newton-Raphson残差–在“SolutionInformation”详细窗口中，输入平衡迭代的数目来获取Newton-Raphson残差.例如，输入“4”,当求解中断或不收敛时，最后四次迭代的残余力将会给出.WorkbenchMechanical-Diagnostics7-22TrainingManual…Newton-Raphson残差–当求解停止或不能收敛，就可以从“SolutionInformation”分支中获得残差，如下图所示.如果求解收敛失败或用户中止,可得到残差的数量.看残差，可以检查出现不平衡力的位置.有助于用户识别可能的有问题位置,以致采取矫正措施.WorkbenchMechanical-Diagnostics7-23TrainingManual…Newton-Raphson残差•在下例中，力残差十分高(不平衡),Newton-RaphsonResiduals使用户可以看到哪些区域对造成高的不平衡力有贡献–在这个例子中,高残余力的‘一圈’是接触区域的