首页/文章/ 详情

快速学会一项分析-板材弯曲工况的非线性隐式分析-OS-T:1500

28天前浏览535
OptiStruct的非线性隐式分析能够处理材料非线性、几何非线性以及接触非线性等问题。通过增量迭代的方式逐步逼近真实路径,使用牛顿-拉夫森方法求解非线性平衡方程,通常能够提供快速的收敛速度。在处理诸如大变形、塑性变形、接触和约束等高度非线性问题时,这种方法表现出了强大的鲁棒性。此外,还提供了自动时间步长控制功能,能够根据当前增量的收敛难度来调整下一步的时间步长,如果迭代次数少于设定的最优迭代次数,时间步长会增加;反之,如果迭代次数过多,则会减小时间步长。非线性隐式分析通常用于那些需要考虑加载路径历史依赖性的工程应用,如金属成形、轮胎分析和医疗器械分析等领域。
本教学案例通过模拟平板在恒定压力下的弯曲,演示OptiStruct中的非线性大位移分析。
在开始之前,请将本教程中使用的文件复 制到您的工作目录。
http://majorv.help.altair.com/minorv/simulation/tutorials/hwsolvers/optistruct/OS-T-1500/plate.zip
此外,还提供了与线性分析结果的比较,以说明进行非线性大位移分析的要求。本教程考虑了几何非线性(大位移)和弹塑性材料。
对于static分析,能量方程简化为F=Ku,其中求解未知位移(u)。在非线性大位移分析中,刚度(K)可能是材料、几何体和边界条件的函数。因此,采用增量迭代方法来计算未知位移。
可以了解线性分析中考虑了哪些条件。首先,在分析过程中,材料是被认为是线性的,如图1 所示。其次,几何变形应该足够小,以确保初始形状的变化不会显著到足以影响刚度。最后,在分析过程中,边界条件不应发生变化。
图1.
该分析同时考虑了几何和材料非线性。考虑几何非线性是因为在几何结构中观察到的大位移。此外,合成应力超过屈服应力极限,这意味着材料不符合线性应力-应变行为,并且开始观察到塑性效应。
模型说明如下:
图2 说明了本教程中使用的结构模型:一个两端支撑的矩形长板,并在顶面上施加了分布式荷载。板的尺寸和材料参数,如下所示。    
图2.板上的约束和压力载荷的表示
Units
Length: mm
Time: s
Mass: ton
Force: N
Stress: MPa
Length:1000 mm
Width:200 mm
Thickness :4.0 mm
Material Steel, elastoplastic
Initial Density:7.86e-9 t/mm3
Young's modulus (E):200000 MPa
Poisson coefficient:0.29
Yield stress :100.0 MPa
Tangent modulus:20000.0 MPa
Imposed:压力0.02 MPa,垂直于板施加
本教程包括以下练习:
·创建plastic Material和相应的壳属性
·设置边界条件和施加的载荷
·设置非线性和线性分析
·提交作业并查看结果

一、启动HyperMesh并设置OptiStruct用户配置文件

1.启动HyperMesh。此时将打开User Profile对话框
2.选择OptiStruct,然后单击OK
这将加载用户配置文件。它包括相应的模板、宏菜单和导入阅读器,将HyperMesh的功能缩减为与生成OptiStruct模型相关的功能。

二、导入模型    

1.点击File>Import>Solver Deck
导入选项卡将添加到您的选项卡菜单中。
2.对于File type ,选择OptiStruct
3.选择文件图标。
此时将打开Select OptiStruct文件Browser。
4.选择保存到工作目录的plate.fem文件。
5.单击Open
6.单击Import,然后单击Close以关闭Import选项卡。
三、设置模型
3.1更新Material
1.在Model Browser中,点击Material MAT1_1
2.输入值,如下所示。有关更多信息,请参阅MAT1。
最初,MAT1 card image表示线性各向同性材料。
3.选中MATS1旁边的框,为非线性分析的双线性材料弹塑性定义其他材料属性。
LIM1 是材料屈服点,在本例中为屈服应力。
H是塑性模量(加工硬化斜率),它定义了塑性区域中应变和应力之间的线性关系。对于双线性材料曲线,它与Young's模量通过切线模量(ETET) 相关。有关更多信息,请参阅MATS1。
图3.应力-应变曲线中的杨氏模量和切线模量以及它们相关的方程    
图4.
3.2Update属性
1.在Model Browser中,点击属性PSHELL_1
2.输入值,如下所示。有关更多信息,请参阅PSHELL。   
图5.
3.3创建边界条件
1.在Model Browser中,右键单击并选择Create>Load Collector
2.对于Name ,输入LC_SPC。
3.单击Color并从调色板中选择一种颜色。
4.对于Card Image,选择None
5.在Analysis页面中,选择constraints,然后切换create
6.将entity selector切换到节点并选择节点。
图6.
7.选择dof1dof3。取消选择其他自由度
8.对于载荷类型,请选择SPC
9.单击create创建边界约束。
10.接下来,从对面的边缘选择节点。    
图7.
11.选择dof3。取消选择所有其他这会限制z方向的平移移动。
12.对于载荷类型,请选择SPC
13.单击create创建边界约束。
14.单击return返回主菜单。
3.4创建Uniform压力
1.在Model Browser中,右键单击并选择Create>Load Collector
2.对于Name ,输入LC_IMPLOAD。
3.单击Color并从调色板中选择一种颜色。
4.对于Card Image,选择None
5.在Analysis页面中,单击pressures并切换create
6.将entity selector切换到elems,然后选择all elements
7.点击magnitude=旁边的切换按钮,并将压力定义方法切换为常量。
8.对于magnitude,输入-0.02。
9.对于载荷类型,选择PLOAD作为载荷类型。
10.单击create>return
3.5定义非线性分析参数
1.在Model Browser中,右键单击并选择Create>Load Step Inputs
2.对于Name ,输入LC_NLPARM。
3.对于Config type ,选择Nonlinear Parameters
默认类型为NLPARM
4.选择DT并输入0.0 作为值。
5.选择CONV作为UPW,EPSUP和ESPS为0.001,EPSW为1e-7,这是默认值。    
CONV标志选择收敛标准。在这种情况下,是位移(U)、载荷(P)和功(W)准则的组合。EPSU、EPSP和EPSW分别是它们各自的误差容限,UPW
6.输入值,如下所示。有关更多信息,请参阅NLOUT。
图8.
3.6创建Output Parameters
1.在Model Browser中,右键单击并选择Create>Load Step Inputs
2.对于Name ,输入LC_NLOUT。
3.对于Config type ,选择Output Parameters
默认类型将为NLOUT
4.选中NINT旁边的框,然后输入10
此参数设置中间结果的输出中显示的区间数。
5.输入值,如下所示。有关更多信息,请参阅NLOUT。 
图9.
图10.
3.7定义输出控制参数
1.在Analysis页面中,选择control cards
2.单击GLOBAL_OUTPUT_REQUEST
3.对于DISPLACEMENT、ELFORCE、OLOAD、STRESS和STRAIN,将Option设置为Yes
4.单击return两次以转到主菜单。
3.8激活大位移非线性分析
1.在Analysis页面中,选择control cards
2.从控制卡中,选择PARAM
3.选择HASHASSMLGDISP
对于非线性分析,建议设置HASHASSM卡以减少内存需求,并设置LGDISP以激活大位移非线性static分析。
4.对于HASHASSM,将Option设置为YES
5.对于LGDISP,设置为1    
6.单击return返回主菜单。
3.9创建非线性static分析SUBCASE
在此步骤中,将创建线性和非线性Load Case,它们具有相同的边界条件和施加的载荷。
1.在Model Browser中,右键单击并选择Create>Load Step
2.对于Name ,输入nonlinear_lgdisp。
3.对于Analysis type ,选择Nonlinear static
4.对于SPC,请选择Unspecified>Loadcol
5.Select Loadcol对话框中,从Load Collector列表中选择LC_SPC,然后单击OK
这将设置上面创建的边界条件。
6.对于LOAD,选择Unspecified>Loadcol
7.Select Loadcol对话框中,从Load Collector列表中选择LC_IMPLOAD,然后单击OK
这将设置上面创建的施加载荷。
8.对于NLPARM,选择Unspecified>Load step inputs
9.Select Load step inputs对话框中,从Load step inputs列表中选择LC_NLPARM,然后单击OK
这将设置上面创建的分析参数。
10.对于NLOUT,选择Unspecified>Load step inputs
11.Select Load step inputs对话框中,从Load step inputs列表中选择LC_NLOUT,然后单击OK
这将设置上面创建的output control参数。非线性分析已设置完毕。
12.在Model Browser中,右键单击并选择Create>Load Step
13.对于Name ,输入linear。
14.对于Analysis type ,选择Linear Static
15.对于SPC,请选择Unspecified>Loadcol
16.Select Loadcol对话框中,从Load Collector列表中选择LC_SPC,然后单击OK
17.对于LOAD,选择Unspecified>Loadcol
18.Select Loadcol对话框中,从Load Collector列表中选择LC_IMPLOAD,然后单击OK
线性分析已设置。
四、提交作业    
1.在Analysis页面中,单击OptiStruct面板。
图11.访问OptiStruct面板
2.单击save as
3.Save As对话框中,指定写入OptiStruct模型文件的位置,并在文件名中输入plate
对于OptiStruct求解器模型,建议使用.fem扩展名。
4.单击Save
input file字段显示在Save As对话框中指定的文件名和位置。
5.将导出选项切换设置为all
6.将run options切换设置为analysis
7.将内存选项切换设置为memory default
8.单击OptiStruct启动OptiStruct作业。
如果作业成功,则新的结果文件应位于写入plate.fem的目录中。plate.out文件是查找错误消息的好地方,如果存在任何错误,这些错误消息可以帮助调试输入模型。
写入目录的默认文件为:
plate.html:HTML分析报告,提供问题表述和分析结果的摘要。
plate.out:OptiStruct输出文件,包含有关文件设置的特定信息、优化问题的设置、运行所需的RAM和磁盘空间量的估计值、每次优化迭代的信息以及计算时间信息。查看此文件是否有警告和错误。
plate.h3d:HyperView二进制结果文件。
plate.stat:摘要,提供分析过程中每个步骤的CPU信息。
在输出文件(plate.out) 中,可以遵循根据定义非线性分析参数中选择的标准参数计算的收敛迭代。此外,还计算了每个步骤的最大塑性应变。 
图12.
五、查看结果
1.使用HyperView打开结果文件(plate.out),并绘制SUBCASE1 和SUBCASE2 在100% 载荷(载荷因子= 1.0)下的位移和von Mises应力云图。
2.比较结果。
图13.位移
图14.应力(von Mises)    
图15.塑性应变
             


来源:TodayCAEer
OptiStructHyperMeshHyperView非线性材料控制Altair
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-10-18
最近编辑:28天前
TodayCAEer
本科 签名征集中
获赞 18粉丝 27文章 242课程 0
点赞
收藏
作者推荐

快速学会一项分析-支架的直接法瞬态动力分析OS-T:1310

直接瞬态动力分析是一种在物理空间中直接求解瞬态动力学方程的方法,它不依赖于模态分析。在 OptiStruct 中,这种分析通常使用 Newmark-β 逐步积分算法来实现。直接法的优点是计算精度高,因为它没有模态截断导致的误差。然而,这种方法的缺点是随着模型自由度的增加,计算复杂度会显著提高,尤其是在处理大规模有限元模型时。在本教程中,使用支架的现有有限元模型来演示如何使用 OptiStruct 进行直接瞬态动态分析。HyperGraph 用于对支架在瞬态动态载荷下的变形特性进行后处理。在开始之前,请将本教程中使用的文件复 制到您的工作目录。http://majorv.help.altair.com/minorv/simulation/tutorials/hwsolvers/optistruct/OS-T-1310/bracket_transient.zip图 1.支架的有限元模型支架两条腿的底部被约束。瞬态动态载荷将施加在支架顶部平面的网格点处,该平面围绕孔在负 z 方向上。加载的时间历程如图 2 所示。直接瞬态分析的总运行时间为 4 秒,时间分为 800 个增量(即时间步长为 0.005)。该模型考虑了结构阻尼。在Spider中心定义了一个集中质量单元,并在该孔中心的集中质量处监测 z轴位移。 图 2.施加载荷的时间历史一、启动 HyperMesh 并设置 OptiStruct 用户配置文件1.启动 HyperMesh。 此时将打开User Profile对话框。2.选择OptiStruct并单击OK。这将加载用户配置文件。它包括相应的模板、宏菜单和导入阅读器,将 HyperMesh 的功能缩减为与生成 OptiStruct 模型相关的功能。二、导入模型1.点击File>Import>Solver Deck。导入选项卡将添加到您的选项卡菜单中。2.对于 File type ,选择OptiStruct。3.选择文件图标 。此时将打开 Select OptiStruct 文件浏览器。4.选择保存到工作目录的bracket_transient.hm文件。5.单击Open。6.点击Import,然后点击Close关闭 Import 选项卡。三、搭建模型3.1 创建TABLED1曲线1.在模型浏览器中,右键单击并选择Create>Curve。2.对于 Name ,输入 tabled1。3.在 Curve Editor 窗口中,输入如图 3 所示的值。 图 3.显示加载时间历史的曲线4.关闭Curve 编辑器。5.在 Curves中,选择tabled 1。6.单击Color并从调色板中选择一种颜色。7.对于 Card Image,从下拉菜单中选择TABLED 1。8.单击Close。定义加载时间历程的 TABLED1 已创建。3.2 创建 TSTEP Load Collector1.在 Model Browser 中,右键单击并选择Create>Load Collector。2.对于 Name ,输入 tstep。用于定义生成和输出解决方案的时间步长间隔的瞬态时间步长。3.单击Color并从调色板中选择一种颜色。4.对于 Card Image,从下拉菜单中选择TSTEP。5.对于 TSTEP_NUM,输入 1 并按Enter。6.对于 N,输入时间步长数 800。7.对于 DT,输入 0.005 的时间增量。施加到载荷的总时间为:800 x 0.005 = 4 秒。这是请求 output 的时间步长。NO 的默认值为 1.0。8.单击Close。3.3 创建 DAREA Load Collector 定义支架顶面上的力。1.在 Model Browser 中,右键单击并选择Create>Load Collector。2.对于 Name ,输入 darea。3.单击Color并从调色板中选择一种颜色。4.对于 Card Image,选择NONE。5.单击BCs>Create>Constraints以打开约束面板。6.点击nodes>by sets。将显示两组set。7.选择force并单击select。属于set force的节点将被选中。8.单击除 dof3 之外的所有自由度 (dof3) 旁边的框,取消选中每个自由度 (dof3) 旁边的框,表示 dof3 是唯一活动的自由度。9.对于 dof3,输入值 -1500。10.对于 load types=,选择DAREA。11.单击create。这将创建一个 1500 个单位的力,该力应用于负 z 方向上的选定节点。12.点击return返回主菜单。3.4 创建 TLOAD Load Step 输入1.在模型浏览器中,右键单击并选择Create>Load Step Inputs。2.对于 Name ,输入 tload1。3.对于 Config type ,从下拉列表中选择Dynamic Load – Time Dependent4.对于 Type,从下拉菜单中选择TLOAD 1。5.对于 Exciteid ,单击Unspecified>Loadcol。6.在Select Loadcol对话框中,从Load Collector 列表中选择darea。7.单击OK以完成选择。8.同样,为 TID 字段选择tabled 1曲线(已定义加载的时间历程)。激励类型可以是施加的载荷 (力或力矩)、强制位移、速度或加速度。在Load Step inputs中的[TLOAD]字段定义了载荷的类型。默认情况下,该类型被设置为applied load。3.5 创建Load Step Load Step 选择 Direct Transient Analysis。1.在模型浏览器中,右键单击并选择Create>Load Step。默认load step模板现在显示在 Model Browser 下方的 Entity Editor 中。2.对于 Name ,输入 transient。3.对于 Analysis type,从下拉菜单中选择Transient(direct)。4.从Select Loadcol对话框中,选择spcs。5.对于 DLOAD,从 Load Step Inputs弹出窗口中选择tload。6.激活TSTEP(TIME)并选择之前创建的Load Collector tstep。将创建一个subcase,用于指定直接瞬态动态分析的载荷和边界条件。3.6 创建阻尼参数1.点击Setup>Create>Control Cards进入控制卡面板。2.单击next以查看更多卡片。3.单击PARAM以定义参数卡。4.向下滚动以激活G,单击G_V1,然后输入 0.2。此参数指定直接瞬态动态分析的均匀结构阻尼系数。5.向下滚动以激活W3,点击W3_V1,输入 300。此参数在瞬态分析中用于将结构阻尼转换为等效黏性阻尼。6.点击return。3.7 创建输出请求1.单击GLOBAL_OUTPUT_REQUESTS并选择DISPLACEMENT并将 FORMAT 下方的空白处留空。2.对于 FORM(1),请选择BOTH。3.对于 OPTION(1),选择SID。此时将显示一个标记为 SID 的黄色按钮。4.双击SID并选择center。5.选择center选项。此集表示位于连接到质量元素(即节点 395)的 spider 中心的节点。6.点击return>next。7.单击OUTPUT。8.在 number_of_outputs = 下,输入 2。9.对于 KEYWORD,选择H 3D和HGTRANS。 10.对于 FREQ,请同时选择ALL。11.单击return两次以退出 Control Cards 面板。四、保存数据库设置保存文件的目录。1.点击File>Save as>Model。2.对于 File name,输入bracket_transient_direct.hm。3.单击Save。五、提交工作1.在 Analysis 页面中,单击 OptiStruct 面板。 图 4.访问 OptiStruct 面板2.单击save As。3.在Save As对话框中,指定写入 OptiStruct 模型文件的位置,并在文件名中输入 bracket_transient_direct。对于 OptiStruct 输入模型,建议使用 .fem 扩展名。4.单击Save。input file 字段显示在Save As对话框中指定的文件名和位置。5.将导出选项切换设置为all。6.将 run options toggle 设置为analysis。7.将内存选项切换设置为memory default。8.单击OptiStruct启动 OptiStruct 作业。如果作业成功,则新的结果文件应位于写入 bracket_transient_direct.fem 的目录中。bracket_transient_direct.out 文件是查找错误消息的好地方,如果存在任何错误,这些消息可以帮助调试输入模型。写入目录的默认文件为:bracket_transient_direct.html:HTML 分析报告,提供问题表述和分析结果的摘要。 bracket_transient_direct.out:OptiStruct 输出文件,包含有关文件设置、优化问题设置、运行所需 RAM 和磁盘空间量的估计值、每次优化迭代的信息以及计算时间信息的特定信息。查看此文件是否有警告和错误。bracket_transient_direct.h3d:HyperView 二进制结果文件。bracket_transient_direct.res:HyperMesh 二进制结果文件。bracket_transient_direct.stat:Summary,提供分析过程中每个步骤的 CPU 信息。bracket_transient_direct.mvw:HyperView 会话文件。只有在执行瞬态分析时,才会创建此文件。此文件会自动为文件中包含的位移、速度和加速度结果创建绘图。六、后处理位移结果1.在 OptiStruct 面板中,单击 HyperView 以启动 HyperView。2.点击File>Open>Session。3. 从运行输入文件的目录中选择 HyperView 会话文件 bracket_transient_direct.mvw。将显示以下提示: 图 5.4.单击Yes关闭消息窗口。由于载荷仅施加在 z 方向上,因此您对节点 395 的 z 位移时间历程感兴趣。此文件会自动为文件中包含的位移结果创建绘图。5.点击Curve Attributes工具栏图标 并关闭曲线X Trans和Y Trans。这可以通过选择单个曲线(X Trans 和 Y Trans),然后单击线条属性Off来完成,如下所示: 图 6.6.单击 以适合 GUI 中节点 395 的 y 轴(即 Z 位移)。7.如果需要,您可以更改曲线的颜色和/或线条属性。 图 7.用于直接瞬态动力学分析的 Spider 中心集中质量的 Z 位移时程从上图中可以看出,节点 395 的位移在负 z 方向上,因为载荷也在 -z 方向上。由于模型中存在结构阻尼,位移最终会减弱。 来源:TodayCAEer

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈