1、输出控制在线性分析中，ABAQUS仅求解一次平衡方程，并以此解答来计算结果。非线性 分析可以产生更多的输岀，因为在每一个收敛的增量步结朿时都可以要求输出结果。如 果你不注意选择输出要求，输出文件会成为非常之大，潜在地占满你的汁算机的磁盘空 间。如前所述，数据输出有四种不同的文件形式：输出数据库（.odb）文件，它包含以二进制格式存储的数据，需要应用 ABAQUS/CAE后处理结果：数据（.dat）文件，它包含了选左结果的数据报表（仅应用于ABAQUS/Standard）；重启动（.res）文件，应用于继续分析；结果（.fil）文件，由第三方后处理器使用的文件。这里只讨论输出数据库（.odb）文
2、件。如果注意选择，在模拟过程中可以经常存储 数据，而又不会过多地占用磁盘空间。从主菜单栏中，选择 0utput-Field Output Requests-Manager,打开 Field Output Requests Manager,在对话框的右边，点击Edit来打开场变量输出编借器。 在Output Variables （输出变量）域中,选择Preselected defaults.删除对线性分 析模型定义的场变量输出要求，并指泄默认的场变疑输岀要求。对于一般的静态过程, 这个输出变疑的预选设這是最经常应用的场变量输出设置。为了减小输出数据库文件的尺寸，选择在每第二个增量步写一次场变量输
3、出。如果 你是简直地感兴趣最终的结果，你也可以或者选择The last increment （最终增量步） 或者设置保存输出的频率等于一个大数。不论指泄什么值，在每个分析步结束时总会保 存结果；所以，使用一个大数会导致仅保存最终的结果。从前而的分析中，可以保留指泄在跨中盯点位移的历史输岀，我们将在Visualization 模块中应用X-K曲线图功能演示这些结果。运行及监控作业在Job模块中，为非线性（Nonlinear）模型创建一个作业,命名为NlSkewPlate, 并给出描述为Nonlinear Elastic Skew Platea记住将你的模型保存为一个新 的模型数据库文件。提交作业
4、进行分析并监控求解进程。如果遇到了任何错误，必须纠正它们：如果发 出了任何警告信息，必须调查它们的来源，并在必要时采取纠正的措施。对于这个非线性斜板例题，图8-12显示了 Job Monitor （作业监视器）的内容。 第一列显示了分析步序号，在本例中只有一个分析步。第二列给出了增疑步序号。第六 列显示了在每个增量步中为了得到收敛解，ABAQUS/Standard所需要的迭代步的数目； 例如，在增量步1中，ABAQUS/Standard需要3次迭代。第八列显示了已经完成的总 的分析步时间，第九列显示了增量步的大小（AT）。这个例子显示了 ABAQUS/Standard如何自动地控制增量步的大小
5、，即在每个增量 步中载荷施加的比例。在这个分析中，ABAQUS/Standard在第一个增量步中施加了总 载荷的10%：你指泄了初始增量”沏为0.1和分析步的总时间为1.0。在第一个增量 步.ABAQUS/Standard需要3次迭代才收敛到解答。在第二个增量步,ABAQUS/Standard 只需要2次迭代，因此，它自动地对下一个增量步的值增加了 50%,达到A7-=0.15o 在第四个和第五个增量步，ABAQUS/Standard也增加了 &T。它调整最后一个增量步的 值使得分析步刚好完成；任本例中，最后增量步的值为0.0875oJ（）: NISkewPlate Status: Compl
6、etedStepInGfcmcntAUSevere Dhcon lleiEquil ItcvTotal IterTotal Tirmc2FrcqSleo Timc/LPFTime 儿 PF Inc130.1ai20.2Q20.350350.15405750 575022550.91250.S1250.33751.J60.0875Log | Enocs Warnings | Output |Subrmitted: ThuApr G3 1035:45 2003 Slatted: Sokef Input Fie Processxx Completed: Solver Input File Procw
7、&or Started: ABAQUS/Stondocd Completed: ABAQU$/S landed Completed： Thu Apr 03 1 0:36:07 2003肋 Dimiss |图8-12 Job Monitor：非线性斜板分析8.4.2作业诊断ABAQUS/CAE不仅可以让你监控分析作业的过程，而且还提供了一个可视化的诊断工 具帮助你了解这个分析模型的收敛行为，以及在必要时对模型进行调试。ABAQUS/Standard 在输岀数拯库中存储了分析作业的每一个分析步、增量步、尝试计算和迭代的信息。当你运 行每一个作业时，将自动地存储诊断的信息。如果分析运算时间超岀了预先
8、估计的时间，或 者过早地被中断，你可以观察由ABAQUS/CAE提供的作业诊断信息，以帮助查找问题的原 因和修改模型的方法。进入Visulization模块，并打开输岀数据库NlSkewPlate . odb以检查收敛历史。从主 菜单栏中，选择Tools-Job Diagnostics打开Job Diagnostics （作业诊断）对话框。在 Job History （作业历史）列表中，点击“ + ”号以扩展列表，它包括了在分析作业中的分 析步、增量步、尝试计算和迭代列表。例如，在lncrement-1下，选择Attempt-1,如图 8-13所示。图8-13第一个増量步的第一次尝试计算的信息
9、摘要在对话框右侧的Attempt Summary （尝试计算信息摘要）中包含了基本信息，如增量 步大小和迭代尝试次数等。选择本次尝试计算的lteraction-1査看关于第一次迭代的详细信 息。在Summary （摘要）页中的信息表明在本次迭代并没有达到收敛，所以点击Residuals （残差）页以便查明原因。如图8-14所示，Residuals页显示了在模型中的平均力/和时间平均力犷的值。它 也显示了最大作用力残差，爲、最大位移增量山匚和最大位移修正值匚，以及发生这些值的节点和自由度。在对话框的底部,通过选择Highlight selection in viewpoint （在视图 窗高亮度
10、显示），可以在视图窗的模型中髙亮度显示发生这些节点和自由度的任何位置。诊 断标准的选择是实时跟踪的，所以你可以在对话框左边的迭代列表中快速浏览，以查看在迭 代过程中视图窗模型相对于判断准则位置的变化。如果你正在试图调试大型、复杂的模型, 这可能是非常有用的。类似的显示可用于查看转动自由度（在Variables （变量）列表中， 选择Rotation （转动）。图8-14第一次迭代的作用力残差信息在这个例题中，在分析步泄义中指泄了初始时间增量为O.lSo关于增量步的平均力为 30.29 N：由于这是第一个增量步，它与时间平均力的值相同。在这个模型中，最大残余 力唸t是-749.6 N,它明显地大
11、于0.005X qa 0 rax出现在节点编号167的自由度1上。由 于包含了壳单元，ABAQUS/Standard还必须检查在模型中力矩的平衡。力矩/转动场也未能 满足平衡检査。尽管不满足平衡检查就足以使ABAQUS/Standard尝试新一轮的迭代，但是你也应该检 査位移修正值。在第一个分析步的第一个增量步的第一次迭代中，位移的最大增量卜点“和 最大位移修正值c爲均为-5.587X 10-3 m：并且转动的最大增量和转动修正值都是-1.598X IO 2弧度。由于在第一个分析步的第一个增量步的第一次迭代中，增量值与修正值总是相等 的，所以关于节点变量的最大修正值是小于1%最大增量值的检验将
12、总是失败的。然而，如 果ABAQUS/Standard判左结果是线性的（基于残差量值的判断，猛 血访“），就会忽 略该准则。由于ABAQUS/Standard在首次迭代中未找到平衡解答，因此它尝试了第二次迭代。第 二次迭代的残差信息如图8-15所示。图8-15第二次迭代的作用力残差信息在第二次迭代中，在肖点167的自由度1上尸爲已降至-073N。然而，由于0.005X严 仍比唸小,其中?a=2.49N,在此次迭代中平衡尚未得到满足。最大位移修正准则也未能 满足，因为发生在节点5的自由度1上的位移修正值c爲二7.055XI0-5大于1%的最大位移 增量唸严5.584X10%在第二次迭代中，力矩残
13、差值检査和最大转动修正值检查都是满足的：然而，ABAQUS/Standard必须进行列一次迭代，因为解答未能通过作用力残差值检査（或最大位 移修正值准则）。图8-16显示了在第一个增量步中需要得到平衡解所做的又一次迭代的残差 信息。图8-16第三次迭代的作用力残差信息在第三次迭代后，严=2.476 N和在节点86的自由度2 =-5.855X 103Na这些值满足rLStep/Frame.弹出Step/Frame （分析步/画而）对话框。在分析时，ABAQUS/Standard根据要求在每第二个增疑步将场变量输出结果写入到输出数据库文件。ABAQUS/CAE显示画而列表，如图818所示。图 81
14、8 iuD ifri表中列出了储存场变量的分析步和增疑步。此分析中只包含一个分析步和6个增量 步，已经默认地保存了关于增量步0的结果（即分析步的初始状态），并按照要求 保存了第2、4和6增量步的结果。默认情况下，ABAQUS/CAE总是使用保存在输 出数据库文件中的最后一个增量步的数据。2.点击0K关闭分析Step/Frame对话框。显示变形前后的模型形状将未变形图叠加在变形图上，一起显示变形前后的模型形状。旋转视图得到类似于图 8J9所示的图形。图819斜板变形前和变形后的模型形状应用来自其它画面的结果从保存在输出数据库文件中的其它增量步数拯中，你可以选择适当的画而来评估结果。 选择一个新的
15、画面：1.从主菜单栏中，选择Result-显示Step/Frame对话框。2从Frame菜单中，选择Increment 4 （增虽:步4）。3.点击0K应用这些变化，并关闭Step/Frame对话框。现在，所需要的任何绘图将使用来自增量步4的结果。重复这个过程，应用所感兴 趣的增量步加以替换，自如地调用在输出数拯库文件中的数据。xy曲线图对于模拟中的每一个增量步，你保存了跨中肖点（节点集合Midspan）的位移作为输 出到数据库文件NlSkewPlate.odb中的历史变量部分，你可以使用这些结果来绘制XY 曲线图。特别是，你将绘制位于板跨中边界处节点的竖向位移历史。创建跨中位移的图形：1.首
16、先，创建一个显示组（display group）,它包括节点集Midspan中未变形的模型 图，显示出节点号以确定那些位于板跨中边界处的节点。2.从主菜单栏中，选择 Result-History Output.3.在弹出的History Output对话框中，选择（用Ctrl+点击）两个跨中边界节点的竖 向运动。曲线标注的形式为：Spatial displacement: U3 at Node xxx in NSET Midspan （用节点编号确泄你需要选择的曲线）4.点击 PlotoABAQUS从输出数据库文件中读出两条曲线的数据，并画出类似于图8-20所示的 曲线图（为了淸楚，第二条曲线
17、已变为虚线）。 U3 W: 1 HSBT MIDSPAN U3 tr： 2 NSET MIDSPAMTime图8-20在板跨中边界的位移历史从这些曲线中可以淸楚地看到该模拟的非线性性质：随着分析的进行，板会逐渐变 硬。几何非线性的效应意味着结构的刚度将随着变形而改变。在该模拟中，由于薄 膜效应使板当变形时变得刚硬。因此，所得到的位移唸值比线性分析预测的小，因 为在线性分析中没有包括这种效应。应用保存在输出数据库文件（.odb）中的历史变量数据或场变量数据，你可以创 建X-Y曲线图。X-Y曲线的数据也可从外部文件读入，或者交互地键入到Visulization 模块中。一旦创建了曲线，可以进一步利
18、用这些数据，并以图形的形式绘制到屏幕上。在第10章“材料”中将进一步讨论Visualization模块的X-Y曲线图功能。数据报表创建一个跨中位移的数据报表。应用丹点集合Midspan创建一个适当的显式组。 报表内容显示如下。Source 1ODB: NlSkewPlate. odbStep: Apply presBure11Frane： Increment 6: Step Time = 1.000Loc 1 : NodaL vaLues from source 1Output sorted by column Node LabelNode LabelU.U1 Loc 1U.U2 Loc 1U
19、U3 Loc 1-2.68697B-03-747.394E-06-49.4696E-03-2.27869E-03-806.331E-06-45.4817E-037-2.57405R-02-759.29SK.06-dfi.5985R-028-2.4908SB-03-775.16SE-06-47.7038E-039-2.4038E-03-793.355E-06-46.533E-0366-2.62603E-03-750.246E-06-d9.D086E-0370-2.53836B-03-762.497E-06-48.1876E-0373-2.45757B-03-778.207E-06-47.144E
20、-0376-2.36229E-03-794.069E-06-45.9613E-03Minimum-2.68697E-03AtNodeMaximum-2.27859B-03Nodo将这些位移值与在第5章应用壳单元”中应用线性分析得到的结果进行比较。该 模拟中的跨中最大位移比由线性分析预测的位移约小9%。在模拟中包括非线性几何效 应，减小了板跨中的竖向挠度（U3）。两种分析的另一个区别是在非线性模拟中沿1和2方向有非零挠度。在非线性分析 中，是什么效果使得而内位移U1和U2非零呢？为什么板的竖向挠度会小呢？板变形后成了弯曲形状：在非线性模拟中考虑了几何改变，作为结果，薄膜效应使 得部分载荷由薄膜作用来承受而不是仅由弯曲作用单独承受，这使得板更加刚硬。列外, 始终保持垂直于板而的压力载荷随着板的变形也开始具有沿1和2方向的分量。非线性 分析中考虑了这种刚性效应和压力方向的改变，而在线性分析中这两种效应均未考虑。在线性和非线性模拟之间的差别是相当大的，表明在这种特殊载荷条件下，对于该 板应用线性模拟是不合适的。对于5个自由度的壳单元，如在这个分析中应用的S8R5单元，ABAQUS/Standard 没有输出在节点处的所有转动。
