ABAQUS有限元分析实例非线性斜板.docx

1、ABAQUS有限元分析实例非线性斜板Abaqus有限元分析实例 非线性斜板这个例子如图8-11所示。已经应用ABAQUS/Standard模拟了板的线性响应，现在你将 应用ABAQUS/Standard对它进行重新分析，包含几何非线性的影响。从线性模拟的结果表 明对于此问题非线性的效应可能是重要的，由此次分析的结果，你将判断这个结论是否正确。如果你愿意，可以根据本例题后而的指导，应用ABAQUS/Explicit将模拟扩展到动态 分析。在本手册的在线文档第A.6节“Norlincar skew plate提供了输入文件。当通过 ABAQUS/CAE运行这个输入文件时，将创建关于该问题的完整的分

2、析模型。根据下而给岀 的指导如果你遇到困难，或者如果你希望检查你的工作，则可以运行这个输入文件。在附录 A Example Files中，给出了如何提取和运行输入文件的指导。如果你没有进入ABAQUS/CAE或者其它的前处理器，可以人工创建关于这个问题的输 入文件，关于这方而的讨论，见 Getting Started with ABAQUS/Standard： Keywords Version, 第 7.4 iT “Example： norlincar skew plate。8.4.1修改模型打开模型数据库文件SkewPlate.cae,从主菜单栏中，选择Model-Copy Model-Li

3、near,将名字为Linear的模型复制成名字为Nonlinear的模型。对于非线性斜板模型，你将考虑包含几何非线性效应和改变输出要求。定义分析步进入分析步Step模块，从主菜单栏中，选择Step-Edit-Apply Pressure来编 辑分析步定义。在Edit Step对话框的Basic页中,选中NIgeom (注：几何非线性的 缩写)以考虑几何非线性的效应，并设垃分析步的时间周期为1.0。在Incrementation (增量步)页中，设置初始增量步的值(initial increment size)为0.1。默认的增量步 最大数目(maximum number of incremen

4、ts)为100； ABAQUS可能采用少于这个上限 的增量步数目，但是如果需要髙于这个上限的增量步数目，分析就会中止。你可能希望改变分析步的描述，以反映它现在是一个非线性分析步。输出控制在线性分析中，ABAQUS仅求解一次平衡方程，并以此解答来计算结果。非线性 分析可以产生更多的输岀，因为在每一个收敛的增量步结朿时都可以要求输出结果。如 果你不注意选择输出要求，输出文件会成为非常之大，潜在地占满你的汁算机的磁盘空 间。如前所述，数据输出有四种不同的文件形式：输出数据库（.odb）文件，它包含以二进制格式存储的数据，需要应用 ABAQUS/CAE后处理结果：数据（.dat）文件，它包含了选左结果

5、的数据报表（仅应用于ABAQUS/Standard）；重启动（.res）文件，应用于继续分析；结果（.fil）文件，由第三方后处理器使用的文件。这里只讨论输出数据库（.odb）文件。如果注意选择，在模拟过程中可以经常存储 数据，而又不会过多地占用磁盘空间。从主菜单栏中，选择 0utput-Field Output Requests-Manager,打开 Field Output Requests Manager,在对话框的右边，点击Edit来打开场变量输出编借器。 在Output Variables （输出变量）域中,选择Preselected defaults.删除对线性分 析模型定义的场变

6、量输出要求，并指泄默认的场变疑输岀要求。对于一般的静态过程, 这个输出变疑的预选设這是最经常应用的场变量输出设置。为了减小输出数据库文件的尺寸，选择在每第二个增量步写一次场变量输出。如果 你是简直地感兴趣最终的结果，你也可以或者选择The last increment （最终增量步） 或者设置保存输出的频率等于一个大数。不论指泄什么值，在每个分析步结束时总会保 存结果；所以，使用一个大数会导致仅保存最终的结果。从前而的分析中，可以保留指泄在跨中盯点位移的历史输岀，我们将在Visualization 模块中应用X-K曲线图功能演示这些结果。运行及监控作业在Job模块中，为非线性（Nonlinea

7、r）模型创建一个作业,命名为NlSkewPlate, 并给出描述为Nonlinear Elastic Skew Platea记住将你的模型保存为一个新 的模型数据库文件。提交作业进行分析并监控求解进程。如果遇到了任何错误，必须纠正它们：如果发 出了任何警告信息，必须调查它们的来源，并在必要时采取纠正的措施。对于这个非线性斜板例题，图8-12显示了 Job Monitor （作业监视器）的内容。 第一列显示了分析步序号，在本例中只有一个分析步。第二列给出了增疑步序号。第六 列显示了在每个增量步中为了得到收敛解，ABAQUS/Standard所需要的迭代步的数目； 例如，在增量步1中，ABAQUS

8、/Standard需要3次迭代。第八列显示了已经完成的总 的分析步时间，第九列显示了增量步的大小（AT）。这个例子显示了 ABAQUS/Standard如何自动地控制增量步的大小，即在每个增量 步中载荷施加的比例。在这个分析中，ABAQUS/Standard在第一个增量步中施加了总 载荷的10%：你指泄了初始增量”沏为0.1和分析步的总时间为1.0。在第一个增量 步.ABAQUS/Standard需要3次迭代才收敛到解答。在第二个增量步,ABAQUS/Standard 只需要2次迭代，因此，它自动地对下一个增量步的值增加了 50%,达到A7-=0.15o 在第四个和第五个增量步，ABAQUS/

9、Standard也增加了 &T。它调整最后一个增量步的 值使得分析步刚好完成；任本例中，最后增量步的值为0.0875oJ（）: NISkewPlate Status: CompletedStepInGfcmcntAUSevere Dhcon lleiEquil ItcvTotal IterTotal Tirmc2FrcqSleo Timc/LPFTime 儿 PF Inc1110330.10.1ai1210220.2Q2ai1310220.350350.1514102205750 57502251510330.91250.S1250.33751.J61022110.0875Log | Enoc

10、s Warnings | Output |Subrmitted: ThuApr G3 1035:45 2003 Slatted: Sokef Input Fie Processxx Completed: Solver Input File Procw&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不仅可以让你监控分析作业的过程，而且还提供了一个可视化的诊

11、断工 具帮助你了解这个分析模型的收敛行为，以及在必要时对模型进行调试。ABAQUS/Standard 在输岀数拯库中存储了分析作业的每一个分析步、增量步、尝试计算和迭代的信息。当你运 行每一个作业时，将自动地存储诊断的信息。如果分析运算时间超岀了预先估计的时间，或 者过早地被中断，你可以观察由ABAQUS/CAE提供的作业诊断信息，以帮助查找问题的原 因和修改模型的方法。进入Visulization模块，并打开输岀数据库NlSkewPlate . odb以检查收敛历史。从主 菜单栏中，选择Tools-Job Diagnostics打开Job Diagnostics （作业诊断）对话框。在 Jo

12、b History （作业历史）列表中，点击“ + ”号以扩展列表，它包括了在分析作业中的分 析步、增量步、尝试计算和迭代列表。例如，在lncrement-1下，选择Attempt-1,如图 8-13所示。图8-13第一个増量步的第一次尝试计算的信息摘要在对话框右侧的Attempt Summary （尝试计算信息摘要）中包含了基本信息，如增量 步大小和迭代尝试次数等。选择本次尝试计算的lteraction-1査看关于第一次迭代的详细信 息。在Summary （摘要）页中的信息表明在本次迭代并没有达到收敛，所以点击Residuals （残差）页以便查明原因。如图8-14所示，Residuals页

13、显示了在模型中的平均力/和时间平均力犷的值。它 也显示了最大作用力残差，爲、最大位移增量山匚和最大位移修正值匚，以及发生这些值的节点和自由度。在对话框的底部,通过选择Highlight selection in viewpoint （在视图 窗高亮度显示），可以在视图窗的模型中髙亮度显示发生这些节点和自由度的任何位置。诊 断标准的选择是实时跟踪的，所以你可以在对话框左边的迭代列表中快速浏览，以查看在迭 代过程中视图窗模型相对于判断准则位置的变化。如果你正在试图调试大型、复杂的模型, 这可能是非常有用的。类似的显示可用于查看转动自由度（在Variables （变量）列表中， 选择Rotation

14、 （转动）。图8-14第一次迭代的作用力残差信息在这个例题中，在分析步泄义中指泄了初始时间增量为O.lSo关于增量步的平均力为 30.29 N：由于这是第一个增量步，它与时间平均力的值相同。在这个模型中，最大残余 力唸t是-749.6 N,它明显地大于0.005X qa 0 rax出现在节点编号167的自由度1上。由 于包含了壳单元，ABAQUS/Standard还必须检查在模型中力矩的平衡。力矩/转动场也未能 满足平衡检査。尽管不满足平衡检查就足以使ABAQUS/Standard尝试新一轮的迭代，但是你也应该检 査位移修正值。在第一个分析步的第一个增量步的第一次迭代中，位移的最大增量卜点“和

15、 最大位移修正值c爲均为-5.587X 10-3 m：并且转动的最大增量和转动修正值都是-1.598X IO 2弧度。由于在第一个分析步的第一个增量步的第一次迭代中，增量值与修正值总是相等 的，所以关于节点变量的最大修正值是小于1%最大增量值的检验将总是失败的。然而，如 果ABAQUS/Standard判左结果是线性的（基于残差量值的判断，猛 血访“），就会忽 略该准则。由于ABAQUS/Standard在首次迭代中未找到平衡解答，因此它尝试了第二次迭代。第 二次迭代的残差信息如图8-15所示。图8-15第二次迭代的作用力残差信息在第二次迭代中，在肖点167的自由度1上尸爲已降至-073N。然

16、而，由于0.005X严 仍比唸小,其中?a=2.49N,在此次迭代中平衡尚未得到满足。最大位移修正准则也未能 满足，因为发生在节点5的自由度1上的位移修正值c爲二7.055XI0-5大于1%的最大位移 增量唸严5.584X10%在第二次迭代中，力矩残差值检査和最大转动修正值检查都是满足的：然而，ABAQUS/Standard必须进行列一次迭代，因为解答未能通过作用力残差值检査（或最大位 移修正值准则）。图8-16显示了在第一个增量步中需要得到平衡解所做的又一次迭代的残差 信息。图8-16第三次迭代的作用力残差信息在第三次迭代后，严=2.476 N和在节点86的自由度2 =-5.855X 103

17、Na这些值满足rLStep/Frame.弹出Step/Frame （分析步/画而）对话框。在分析时，ABAQUS/Standard根据要求在每第二个增疑步将场变量输出结果写入到输出数据库文件。ABAQUS/CAE显示画而列表，如图818所示。图 818 iuD ifri表中列出了储存场变量的分析步和增疑步。此分析中只包含一个分析步和6个增量 步，已经默认地保存了关于增量步0的结果（即分析步的初始状态），并按照要求 保存了第2、4和6增量步的结果。默认情况下，ABAQUS/CAE总是使用保存在输 出数据库文件中的最后一个增量步的数据。2.点击0K关闭分析Step/Frame对话框。显示变形前后的

18、模型形状将未变形图叠加在变形图上，一起显示变形前后的模型形状。旋转视图得到类似于图 8J9所示的图形。图819斜板变形前和变形后的模型形状应用来自其它画面的结果从保存在输出数据库文件中的其它增量步数拯中，你可以选择适当的画而来评估结果。 选择一个新的画面：1.从主菜单栏中，选择Result-Step/Frame.显示Step/Frame对话框。2从Frame菜单中，选择Increment 4 （增虽:步4）。3.点击0K应用这些变化，并关闭Step/Frame对话框。现在，所需要的任何绘图将使用来自增量步4的结果。重复这个过程，应用所感兴 趣的增量步加以替换，自如地调用在输出数拯库文件中的数据

19、。xy曲线图对于模拟中的每一个增量步，你保存了跨中肖点（节点集合Midspan）的位移作为输 出到数据库文件NlSkewPlate.odb中的历史变量部分，你可以使用这些结果来绘制XY 曲线图。特别是，你将绘制位于板跨中边界处节点的竖向位移历史。创建跨中位移的图形：1.首先，创建一个显示组（display group）,它包括节点集Midspan中未变形的模型 图，显示出节点号以确定那些位于板跨中边界处的节点。2.从主菜单栏中，选择 Result-History Output.3.在弹出的History Output对话框中，选择（用Ctrl+点击）两个跨中边界节点的竖 向运动。曲线标注的形式

20、为：Spatial displacement: U3 at Node xxx in NSET Midspan （用节点编号确泄你需要选择的曲线）4.点击 PlotoABAQUS从输出数据库文件中读出两条曲线的数据，并画出类似于图8-20所示的 曲线图（为了淸楚，第二条曲线已变为虚线）。 U3 W: 1 HSBT MIDSPAN U3 tr： 2 NSET MIDSPAMTime图8-20在板跨中边界的位移历史从这些曲线中可以淸楚地看到该模拟的非线性性质：随着分析的进行，板会逐渐变 硬。几何非线性的效应意味着结构的刚度将随着变形而改变。在该模拟中，由于薄 膜效应使板当变形时变得刚硬。因此，所得到

21、的位移唸值比线性分析预测的小，因 为在线性分析中没有包括这种效应。应用保存在输出数据库文件（.odb）中的历史变量数据或场变量数据，你可以创 建X-Y曲线图。X-Y曲线的数据也可从外部文件读入，或者交互地键入到Visulization 模块中。一旦创建了曲线，可以进一步利用这些数据，并以图形的形式绘制到屏幕上。在第10章“材料”中将进一步讨论Visualization模块的X-Y曲线图功能。数据报表创建一个跨中位移的数据报表。应用丹点集合Midspan创建一个适当的显式组。 报表内容显示如下。Source 1ODB: NlSkewPlate. odbStep: Apply presBure11

22、Frane： Increment 6: Step Time = 1.000Loc 1 : NodaL vaLues from source 1Output sorted by column Node LabelNode LabelU.U1 Loc 1U.U2 Loc 1UU3 Loc 11-2.68697B-03-747.394E-06-49.4696E-032-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.403

23、8E-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-0376-2.36229E-03-794.069E-06-45.9613E-03Minimum-2.68697E-03-806.331E-06-49.4696E-03AtNode121Maximum-2.27859B-03-747.394E-06-45.4817E-03AtNodo212将这些位移值与在第5章应用壳单元”

24、中应用线性分析得到的结果进行比较。该 模拟中的跨中最大位移比由线性分析预测的位移约小9%。在模拟中包括非线性几何效 应，减小了板跨中的竖向挠度（U3）。两种分析的另一个区别是在非线性模拟中沿1和2方向有非零挠度。在非线性分析 中，是什么效果使得而内位移U1和U2非零呢？为什么板的竖向挠度会小呢？板变形后成了弯曲形状：在非线性模拟中考虑了几何改变，作为结果，薄膜效应使 得部分载荷由薄膜作用来承受而不是仅由弯曲作用单独承受，这使得板更加刚硬。列外, 始终保持垂直于板而的压力载荷随着板的变形也开始具有沿1和2方向的分量。非线性 分析中考虑了这种刚性效应和压力方向的改变，而在线性分析中这两种效应均未考虑。在线性和非线性模拟之间的差别是相当大的，表明在这种特殊载荷条件下，对于该 板应用线性模拟是不合适的。对于5个自由度的壳单元，如在这个分析中应用的S8R5单元，ABAQUS/Standard 没有输出在节点处的所有转动。