《ANSYS12.0结构分析工程应用实例解析第3版》连载14
发表时间： && 作者: 张朝辉&&来源: 机械工业出版社
关键字: &&&
本文是《ANSYS12.0结构分析工程应用实例解析 第3版》连载。由机械工业出版社独家授权e-works转载，任何人不得复制、转载、摘编等任何方式进行使用。如需联系出版相关书籍，请联系机械工业出版社张淑谦先生，电话：010-
8.2& 结构断裂分析实例详解――二维断裂问题 8.2.1& 问题描述 & 图8.5所示为一断裂试样结构示意图，厚度为5mm，试计算其应力强度因子。
& 试样材料参数：弹性模量E=220GPa；泊松比n=0.25；载荷P=0.12MPa
8.2.2& 问题分析 & 由于长度和宽度方向的尺寸远大于厚度方向的尺寸，且所承受的载荷位于长宽方向所构成的平面内，所以该问题满足平面应力问题的条件，可以简化为平面应力问题进行求解。 根据对称性，取整体模型的1/2建立几何模型；选择六节点三角形单元PLANE183模拟加载过程； 先进行普通结构分析求解，再采用特殊的后处理命令计算断裂参数。 8.2.3& 求解步骤 & &&& 1.定义工作文件名和工作标题 1）选择Utility&MenuFileChange&Jobname命令，出现Change&Jobname对话框，在［/FILNAM］&Enter&new&jobname文本框中输入工作文件名EXERCISE1，单击OK按钮关闭该对话框。 2）选择Utility&MenuFileChange&Title命令，出现Change&Title对话框，在文本框中输入ANALYSIS&OF&THE&STRESS&INTENSITY&FACTOR，单击OK按钮关闭该对话框。 & &&& 2.定义单元类型 1）选择Main&MenuPreprocessorElement&TypeAdd/Edit/Delete命令，出现Element&Types对话框，单击Add按钮，出现Library&of&Element&Types对话框。 2）在Library&of&Element&Types列表框中选择Structural&SolidQuad&8node&183，在Element&type&reference&number文本框中输入1，如图8.6所示，单击OK按钮关闭该对话框。
& 3）单击Element&Types对话框上的Options按钮，出现PLANE183&element&type&options对话框，在Element&behavior&K3下拉列表中选择Plane&strs&w/thk，其余选项采用默认设置，如图8.7所示，单击OK按钮关闭该对话框。 4）单击Element&Types对话框上的Close按钮，关闭该对话框。
& 5）选择Main&MenuPreprocessorReal&ConstantsAdd/Edit/Delete命令，出现Real&Constants对话框，单击Add按钮，出现Element&Type&for&Real&Constants对话框，单击OK按钮，出现Real&Constants&Set&Number&2，for&PLANE183对话框，在Real&Constant&Set&No.&文本框中输入1，在Thickness&THK文本框中输入5，如图8.8所示，单击OK按钮关闭该对话框。
& 6）单击Real&Constants对话框上的Close按钮关闭该对话框， & &&& 3.定义材料性能参数 1）选择Main&MenuPreprocessorMaterial&PropsMaterial&Models命令，出现Define&Material&Model&Behavior对话框。 2）在Material&Models&Available一栏中依次单击Structural、Linear、Elastic、Isotropic选项，出现Linear&Isotropic&Propeties&for&Material&Number&1对话框，在EX文本框中输入2.2E5，在PRXY文本框中输入0.25，单击OK按钮关闭该对话框。 3）在Define&Material&Model&Behavior对话框上选择MaterialExit命令，关闭该对话框。 & &&& 4.创建几何模型、划分网格 1）选择Utility&MenuPlotCtrlsNumbering命令，出现Plot&Numbering&Controls对话框，选择KP&Keypoint&numbers、LINE&Line&numbers和AREA&Area&numbers选项，使其状态从Off变为On，单击OK按钮关闭该对话框。
2）选择Main&MenuPreprocessorModelingCreateKeypointsIn&Active&CS命令，出现Create&Keypoints&In&Active&Coordinate&System对话框，在NPT&Keypoint&nuber文本框中输入1，在X,Y,Z&Location&in&active&CS文本框中分别输入50，0，0，如图8.9所示，单击OK按钮关闭该对话框。
& 3）参照上一步的操作过程，依次在ANSYS显示窗口生成以下关键点编号及坐标： 2（100，0，0）；3（100，60，0）；4（-25，60，0）；5（-25，0，0）。 & 分页 & 4）选择Main&MenuPreprocessorModelingCreateLinesLinesStraight&Line命令，出现Create&Straight拾取菜单，用鼠标在ANSYS显示窗口依次选择编号为1、2，2、3，3、4，4、5，5、1的关键点生成5条线段，单击OK按钮关闭该菜单。
5）&选择Main&MenuPreprocessorModelingCreateAreasArbitraryBy&Lines命令，出现Create&Area&by&L拾取菜单，在文本框中输入5，1，2，3，4，单击OK按钮关闭该菜单。
& 6）选择Main&MenuPreprocessorModelingCreateAreasCircleSolid&Circle命令，出现Solid&Circle&Area对话框，在WP&X文本框中输入0，在WP&Y文本框中输入27.5，在Radius文本框中输入12.5，如图8.10所示，单击OK按钮关闭该对话框。 7）选择Main&MenuPreprocessorModelingOperateBooleansSubtractAreas命令，出现Subtract&Areas拾取菜单，在文本框中输入1，单击OK按钮；在文本框中输入2，单击OK按钮关闭该菜单。 8）选择Main&MenuPreprocessorNumbering&CtrlsCompress&Numbers命令，出现Compress&Numbers对话框，在Label&Item&to&be&compressed下拉列表中选择Areas，单击OK按钮关闭该对话框。
& 9）选择Main&MenuPreprocessorMeshingSize&CntrlsConcentrat&KPsCreate命令，出现Concentration&Keypoint拾取菜单，在文本框中输入1，单击OK按钮，出现Concentration&Keypoint对话框，参照图8.11对其进行设置，单击OK按钮关闭该对话框。
& 10）选择Main&MenuPreprocessorMeshingSize&CntrlsManualSizeGlobalSize命令，出现Global&Element&Sizes对话框，在SIZE&Element&edge&length文本框中输入3，其余选项采用默认设置，如图8.12所示，单击OK按钮关闭该对话框。
& 11）选择Main&MenuPreprocessorMeshingMeshAreasFree命令，出现Mesh&Areas拾取菜单，单击Pick&All按钮关闭该菜单。 12）选择Utility&MenuPlotElements命令，ANSYS显示窗口将显示网格划分结果，如图8.13所示。
& 13）选择Utility&MenuFileSave&as命令，出现Save&Database&对话框，在Save&Database&to文本框中输入EXERCISE11.db，保存上述操作过程，单击OK按钮关闭该对话框。 & &&& 5.加载求解 1）选择Main&MenuSolutionAnalysis&TypeNew&Analysis命令，出现New&Analysis对话框，选择分析类型为Static，单击OK按钮关闭该对话框。 2）选择Utility&MenuPlotLines命令，显示所有线段。 3）选择Utility&MenuSelectEntities命令，出现Select&Entities对话框，在第1个下拉列表中选择Lines，在第2个下拉列表中选择By&Num/Pick，在第3栏中选择From&Full单选项，如图8.14a所示。单击OK按钮，出现Select&Lines拾取菜单，在文本框中输入1，单击OK按钮关闭该菜单。
4）选择Utility&MenuSelectEntities命令，出现Select&Entities对话框，在第1个下拉列表中选择Nodes，在第2个下拉列表中选择Attached&to，在第3栏中选择Line,all单选项，在第4栏中选择From&Full单选项，如图8.14b所示，单击OK按钮关闭该对话框。
& 分页 & 5）选择Main&MenuSolutionDefine&LoadsApplyStructuralDisplacementOn&Nodes命令，出现Apply&U，ROT&on&Nodes拾取菜单，单击Pick&All按钮，出现Apply&U,ROT&on&Nodes对话框，在Lab2&DOFs&to&be&contrained列表框中选择UY，在Apply&as下拉列表中选择Constant&value，在VALUE&Displacement&value文本框中输入0，如图8.15所示，单击OK按钮关闭该对话框。 6）选择Utility&MenuSelectEverything命令，选择所有实体。 7）选择Main&MenuSolutionDefine&LoadsApplyStructuralPressureOn&Lines命令，出现Apply&PRES&on&Lines拾取菜单，在文本框中输入6，单击OK按钮，出现Apply&PRES&on&Lines对话框，在［SFL］&Apply&PRES&on&lines&as&a下拉菜单中选择Constant&value，在VALUE&&Load&PRES&value文本框中输入0，在Value文本框中输入0.12，如图8.16所示，单击OK按钮关闭该对话框。
& 8）选择Main&MenuSolutionDefine&LoadsApplyStructuralPressureOn&Lines命令，出现Apply&PRES&on&Lines拾取菜单，在文本框中输入7，单击OK按钮，出现Apply&PRES&on&Lines对话框，在［SFL］&Apply&PRES&on&lines&as&a下拉列表中选择Constant&value，在VALUE&Load&PRES&value文本框中输入0.12，在Value文本框中输入0，单击OK按钮关闭该对话框。 9）选择Utility&MenuWorkPlaneLocal&Coordinate&SystemsCreate&Local&CSAt&Specified&Location命令，出现Create&CS&at&Location拾取菜单，在文本框中输入50，0，0，单击OK按钮，出现Create&Local&CS&at&Specified&Location对话框，参照图8.17对其进行设置，单击OK按钮关闭该对话框。 10）选择Main&MenuSolutionSolve&Current&LS命令，出现Solve&Current&Load&Step对话框，单击OK按钮，ANSYS开始求解计算。
& 11）求解结束时，出现Note提示框，单击Close按钮关闭该对话框。 12）选择Utility&MenuFileSave&as命令，出现Save&Database&对话框，在Save&Database&to文本框中输入EXERCISE12.db，保存求解结果，单击OK按钮关闭该对话框。 & &&& 6.查看求解结果 1）选择Main&MenuGeneral&PostprocPlot&ResultsContour&PlotNodal&Solu命令，出现Contour&Nodal&Solution&Data对话框，在Item&to&be&contoured列表框中选择Nodal&SolutionDOF&SolutionDisplacement&vector&sum，单击OK按钮，ANSYS窗口将显示如图8.18所示的合位移等值线图。 2）选择Main&MenuGeneral&PostprocPlot&ResultsContour&PlotNodal&Solu命令，出现Contour&Nodal&Solution&Data对话框，在Item&to&be&contoured列表框中选择Nodal&SolutionStressvon&Mises&stress，单击OK按钮，ANSYS窗口将显示如图8.19所示的等效应力等值线图。
& 3）选择Main&MenuGeneral&PostprocList&ResultsReaction&Solu命令，出现List&Reaction&Solution对话框，在Lab&&Item&to&be&listed列表框中选择All&items，单击OK按钮，ANSYS显示窗口将列表显示支反力结果，如图8.20所示。
& 分页 & 4）选择Main&MenuGeneral&PostprocPath&OperationsDefine&PathBy&Nodes命令，出现By&Nodes拾取菜单，在文本框中输入2，35，33，单击OK按钮，出现By&Nodes对话框，在Name&Define&Path&Name:&文本框中输入DF，其余选项采用默认设置，如图8.21所示，单击OK按钮关闭该对话框。 5）选择Main&MenuGeneral&PostprocPath&OperationsDefine&PathPath&StatusCurrent&Path命令，ANSYS显示窗口将显示当前路径，如图8.22所示。
& 6）选择Main&MenuGeneral&PostprocNodal&CalcsStress&Int&Factr命令，出现Stress&Intensity&Factor对话框，参照图8.23对其进行设置，单击OK按钮，ANSYS显示窗口将显示应力强度因子计算结果，如图8.24所示。
& 7）选择Utility&MenuFileExit命令，出现Exit&from&ANSYS对话框，选择Quit－No&Save!选项，单击OK按钮，关闭ANSYS。
& 8.2.4& 命令流 & /FILNAME，EXERCISE1 ！定义工作文件名 /TITLE，ANALYSIS&OF&THE&STRESS&INTENSITY&FACTOR &&&&&&&&&&&&&&&&&&&&&&&&& ！定义工作标题 /PREP7&&&&&&&&&&&&&&&&&&&& ！进入前处理器 ET，1，PLANE183&&&&&&&&&&&&&&&&&&&& ！指定单元类型 KEYOPT，1，3，3&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！设置单元关键字 R，1，5&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！定义单元实常数 MP，EX，1，2.2E5&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！输入弹性模量 MP，PRXY，1，0.25&&&&&&&&&&&&&&&&&&&&&&&&&&& ！输入泊松比
/PNUM，KP，1&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！显示关键点编号 /PNUM，LINE，1&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！显示线段编号 /PNUM，AREA，1&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！显示面编号 K，1，50&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！生成关键点 K，2，100 K，3，100，60 K，4，-25，60 K，5，-25 L，1，2&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！生成线段 L，2，3 L，3，4 L，4，5 L，5，1 AL，5，1，2，3，4&&&&&&&&&&&&&&&&&&&&&&& ！由线段生成面 CYL4，0，27.5，12.5&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！生成圆面 ASBA，1，2&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！面相减操作 NUMCMP，AREA&&&&&&&&&&&&&&&&&&&&&&&&& ！压缩面编号 KSCON，1，5，1，10，0.75&&&&&&&&&&&&&&&&&&&&&& ！定义裂纹尖端 ESIZE，3，0&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！定义单元尺寸 MSHKEY，0 AMESH，1&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！对面进行网格划分 /TITLE，ELEMENTS&IN&MODEL EPLOT&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！显示单元 FINISH
/SOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！进入求解器 ANTYPE，STATIC LPLOT&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！显示线段 LSEL，S，，，1&&&&&&&&&&&&&&&&&&&&&&&&& ！选择线段 NSLL，S，1&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！选择线段上的所有节点 D，ALL，UY&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！施加位移载荷 ALLSEL&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！选择所有实体 SFL，6，PRES，0，0.12&&&&&&&&&&&&&&&&&&&&&&&&& ！在线段上施加压力载荷 SFL，7，PRES，0.12，0 LOCAL，11，0，50，0，0，，，，1，1&&&&&&&& ！创建局部坐标系 SOLVE&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！开始求解计算 FINISH
/POST1&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！进入POST1后处理器 PLNSOL，U，SUM&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！绘制合位移等值线图 PLNSOL，S，EQV&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！绘制等效应力等值线图 PRRSOL&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！列表显示支反力求解结果 PATH，DF，3，30，20&&&&&&&&&&&&&&&&&&&&&&&&&& ！定义路径 PPATH，1，2&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！设置路径点 PPATH，2，35 PPATH，3，33 PDEF，STAT&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！显示当前路径 KCALC，1，1，0，0&&&&&&&&&&&&&&&&&&&&&&&&&& ！计算应力强度因子 FINISH /EXIT，ALL&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& ！退出ANSYS
责任编辑：赵栋此页面上的内容需要较新版本的 Adobe Flash Player。
当前位置：>>广州有道有限元>>培训服务成功案例
悬臂梁在循环加载作用下的弹塑性计算(GUI)
10:46&&by:有限元&&来源:广州有道有限元
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 悬臂梁在循环加载作用下的弹塑性计算(GUI)
为考察悬臂梁根部P点的应力-应变历程，采用2D的计算模型，使用平面单元PLANE42，材料采用多线性弹塑性模型(mkin)，进行循环加载过程的分析。
建模的要点：
􀁘设置几何以及材料参数，
􀁙输入材料的多线性弹塑性模型(包括：弹性模量、屈服极限)，见图9-8；
􀁰通过设置time来给出加载历程，每次加载都输入当时的状态载荷值，不是增量加载，每次加载后，必须进行计算，再进入下一步的计算；
􀁱在时间后处理中，通过设置几何位置来查询对应的P观测点的节点编号，并设置观测点的应力显示变量(2号变量)以及塑性应变为显示变量(3号变量)，最后将3号变量设置为横轴，画出2号变量随3号变量的变化曲线见图9-9，可以看出，该材料具有非常明显的Bauschinger效应(即正向屈服与反向屈服之和是单拉实验屈服极限的2倍)。
给出的基于图形界面(GUI)的交互式操作(step by step)过程如下。
(1) 进入ANSYS（设定工作目录和工作文件）
程序 → ANSYS → ANSYS Interactive → Working directory（设置工作目录）→ Initial jobname（设置工作文件名）: Beams → Run → OK
(2) 设置计算类型
ANSYS Main Menu：Preferences… → Structural → OK
(3) 设定不显示时间
ANSYS Utility Menu：PlotCtrls → Window Controls → Window Options… → DATE：No Date or Time → OK
(4) 定义单元类型
ANSYS Main Menu：Preprocessor → Element Type → Add/Edit/Delete... → Add…→ Solid: Quad 4node 42 → OK（返回到Element Types窗口）→ Close
(5) 定义材料参数
ANSYS Main Menu: Preprocessor → Material Props → Material Models → Structural → Linear → Elastic → Isotropic → 输入 EX: 2E4, PRXY: 0.3 (定义弹性模量及泊松比) → OK → 返回Define Material Model Behavior 窗口Structural → NonLinear → Inelastic → Rate Independent → Kinematic Hardening Plasticity → Mises Plasticity → Multilinear (Fixed table) → 在Strain一行中对应1至4号点输入0.004、0.015、0.03、0.08 → 在Curve1中对应1至4号点输入80、160、210、280 → 点击右下角Graph→ OK → Close（关闭材料定义窗口），，见图9-8，观察窗口中的多线性弹塑性模型
(6) 构造模型
生成关键点
ANSYS Main Menu：Preprocessor → Modeling → Create → Keypoints → In Active CS → Keypoints number：1，X，Y，Z Location in active CS：0，0，0 → Apply → 同样依次输入其他三个关键点（100，0，0）、（100，10，0）与（0，10，0）→ OK
ANSYS Main Menu：Preprocessor → Modeling → Create → Areas → Arbitrary → Through KPs → 用鼠标依次点击1、2、3、4关键点，生成面单元
(7) 网格划分
ANSYS Main Menu：Preprocessor → Meshing → Mesher Opts → Mesher Type : Mapped → OK → 2D Shape Key : Quad → OK
ANSYS Main Menu：Preprocessor → Meshing → size contrls → ManualSize → Lines → Picked Lines → 选择上下两条横边线，Ok → NDIV 设置为20 → Apply → 选择两条竖边线 → Ok → NDIV设置为8 → OK
ANSYS Main Menu：Preprocessor → Meshing → Mesh → Areas → Target Surf → 点击生成面几何体的位置，显示矩形面被选中 → OK
(8) 模型加约束
ANSYS Main Menu: Solution → Define Loads →Apply → Structural → Displacement On Lines →选取左侧边线（L4）→ OK → select Lab2: All DOF(施加全部约束) → OK
(9)求解设置
ANSYS Main Menu : Solution → Analysis Type → Sol’n Controls → 在Basic标签下设置Analysis Options 为 Large Displacement Satic，Number of substeps: 8, Max no. of substeps :25 Min no. Of substeps:2, Frequency 设置为Write N number of substeps Where N = 10 → OK
(10)按照时间步施加循环载荷
ANSYS Main Menu : Solution → Analysis Type → Sol’n Controls → 在Basic标签下设置Time at end of loadstep:1 → OK
ANSYS Main Menu : Solution → Define Loads → Apply → Structural → Force/Moment → On Nodes → 选择右侧边缘中点（26号节点）→ OK → Lab：Fy，Value：-40 → OK
ANSYS Main Menu：Solution → Solve → Current LS → OK
ANSYS Utility Menu : Plot → Replot
ANSYS Main Menu : Solution → Analysis Type → Sol’n Controls → 在Basic标签下设置Time at end of loadstep: 2 → OK
ANSYS Main Menu : Solution → Define Loads → Apply → Structural → Force/Moment → On Nodes → 选择右侧边缘中点（26号节点）→ OK→ Lab：Fy，Value：0 → OK
ANSYS Main Menu：Solution → Solve → Current LS → OK
ANSYS Utility Menu : Plot → Replot
ANSYS Main Menu : Solution → Analysis Type → Sol’n Controls → 在Basic标签下设置Time at end of loadstep: 3 → OK
ANSYS Main Menu : Solution → Define Loads → Apply → Structural →Force/Moment → On Nodes → 选择右侧边缘中点（26号节点）→ OK → Lab：Fy，Value：40 → OK
ANSYS Main Menu：Solution → Solve → Current LS → OK
ANSYS Utility Menu : Plot → Replot
ANSYS Main Menu : Solution → Analysis Type → Sol’n Controls → 在Basic标签下设置Time at end of loadstep: 4 → OK
ANSYS Main Menu : Solution → Define Loads → Apply → Structural → Force/Moment → On Nodes → 选择右侧边缘中点（26号节点）→ Lab：Fy，Value：0 → OK
ANSYS Main Menu：Solution → Solve → Current LS → OK
ANSYS Utility Menu : Plot → Replot
ANSYS Main Menu : Solution → Analysis Type → Sol’n Controls → 在Basic标签下设置Time at end of loadstep: 5 → OK
ANSYS Main Menu : Solution → Define Loads → Apply → Structural →Force/Moment → On Nodes → 选择右侧边缘中点（26号节点）→ Lab：Fy，Value：-40 → OK
ANSYS Main Menu：Solution → Solve → Current LS → OK
ANSYS Utility Menu : Plot → Replot
ANSYS Main Menu : Solution → Analysis Type → Sol’n Controls → 在Basic标签下设置Time at end of loadstep: 6 → OK
ANSYS Main Menu : Solution → Define Loads → Apply → Structural →Force/Moment → On Nodes → 选择右侧边缘中点（26号节点）→ Lab：Fy，Value：0 → OK
ANSYS Main Menu：Solution → Solve → Current LS → OK
(11) 计算结果
ANSYS Main Menu：General Postproc → Read Results → Last Set
ANSYS Main Menu：General Postproc → Plot Results → Deformed Shape → Def + Undeformed →
OK（观察最后变形情况）
ANSYS Main Menu：General Postproc → Plot Results → Contour Plot → Element solu → Plastic
Strain → Equivalent plastic strain → OK（观察累计的等效塑性应变）
ANSYS Main Menu：TimeHist Postpro → 关闭弹出窗口 → Define Variables → Add… → Element Results → OK 在方框中输入2 → OK 在方框中输入4 → OK → 在Item，Comp Data item 中选择 Stress, X-direction SX → OK返回Define Time-History Variables →Add… → Element Results → OK 在方框中输入2 → OK 在方框中输入4 → OK → 在Item，Comp Data item 中选择 Strain-plastic, X-dir’n EPPL X → OK → Close
ANSYS Main Menu：TimeHist Postpro → 关闭弹出窗口→ Settings → Graph → Single Variable No. 输入3 → OK
ANSYS Main Menu：TimeHist Postpro → 关闭弹出窗口→ Graph Variables → Nvar1中输入2 → OK观察观测点P上的应力应变历程(SX)，见图9-9。
ANSYS Utility Menu：File → Exit → Save Everything → OK
悬臂梁在循环加载作用下的弹塑性计算(命令流)
给出的命令流如下。
/PREP7 !进入前处理
/PLOPTS,DATE,0 !设置不显示日期和时间
!=====设置几何以及材料参数
FORCE=40 !设置力参数
LTH=100 !设置悬臂梁的长度
THK=10 !设置梁的高度
EXX=20000 !设置弹模
sts1=80 !设置屈服应力
!=====以下为多线性弹塑性模型的应力-应变数据点
stn1=sts1/exx
stn2=0.015
mp,ex,1,EXX
mp,nuxy,1,0.3
tb,mkin,1 !设置多线性随动硬化模型
tbtemp,,strain
tbdata,,stn1,stn2,stn3,stn4 !多线性弹塑性模型的应变数据点
tbdata,,sts1,sts2,sts3,sts4 !多线性弹塑性模型的应力数据点
TBPLOT, mkin,1 !图形显示多线性弹塑性模型实际输入曲线
!=====设置单元以及几何
ANTYPE,STATIC !设置为静力结构分析
ET,1,PLANE42 !选取单元类型1(平面单元PLANE42)
K,1,0,0,0 !生成几何点No.1
K,2,LTH,0,0 !生成几何点No.2
K,3,LTH,THK,0 !生成几何点No.3
K,4,0,THK,0 !生成几何点No.4
A,1,2,3,4 !由4个几何点连成几何面No.1
LESIZE,1, , ,20, , , , ,1 !对线1，画分20段(单元划分前的设置)
LESIZE,3, , ,20, , , , ,1 !对线3，画分20段(单元划分前的设置)
LESIZE,2, , ,8, , , , ,1 !对线2，画分8段(单元划分前的设置)
LESIZE,4, , ,8, , , , ,1 !对线4，画分8段(单元划分前的设置)
MSHAPE,0,2D !设置四边形单元
MSHKEY,1 !设置映射划分
AMESH,1 !对面No.1进行网格划分
DL,4, ,ALL, !对几何线4进行完全约束
FINISH !结束前处理
!=====在求解模块中，按照时间历程施加外力，进行求解
/SOLU !进入求解模块
nlgeom,on !打开几何大变形的开关
outres,all,-10 !每10步存储一次计算结果
nsubst,8,25,2 !每一个步长设置8个子步，最大25步，最小2步
NP=NODE(100,5,0) !获取几何位置为(100,5,0)所对应的节点号码，赋值给NP
!以下按照时间步来施加循环载荷
time,1 !时刻1(第1个时间步)
F,NP,FY,-FORCE !在节点NP上施加力FY=-FORCE
solve !求解
time,2 !时刻2(第1个时间步)
F,NP,FY,0 !在节点NP上施加力FY=0
solve !求解
time,3 !时刻3(第1个时间步)
F,NP,FY, FORCE !在节点NP上施加力FY=FORCE
solve !求解
time,4 !时刻4(第1个时间步)
F,NP,FY,0 !在节点NP上施加力FY=0
solve !求解
time,5 !时刻5(第1个时间步)
F,NP,FY,- FORCE !在节点NP上施加力FY=-FORCE
solve !求解
time,6 !时刻6(第1个时间步)
F,NP,FY,0 !在节点NP上施加力FY=0
solve !求解
!=====进入一般的后处理模块
/post1 !进入一般的后处理
set,last !读取最后一步的计算结果 .
PLDISP,2 !图形显示计算结果的变形后与前的形貌
PLESOL,nl,epeq !图形显示累计的等效塑性应变
FINISH !结束后处理
!=====进入时间历程的后处理模块
/post26 !进入时间历程的后处理
KEY= NODE (10,0,0) !获取位置为(10,0,0)的节点编号，并赋值给KEY
ELM = ENEARN (KEY) !获取最接近节点为KEY的单元号，并赋值给ELM
ESOL,2,ELM,KEY,S,X, !将ELM号单元的KEY号节点的应力SX设置为2号变量
ESOL,3,ELM,KEY,EPPL,X, !将ELM号单元的KEY号节点的弹塑性应变EPPL设置为3号变量
XVAR,3 !将3号变量设置为横坐标
PLVAR,2 !图形显示2号变量随时间变化的曲线
FINISH !结束后处理
相关标签搜索：&&&&&&&&&&&&&&&&&&&&&&
姓 &&&& 名:
请留下您的有效联系方式，以便我们能联系到您 ^v^
客服在线请直接联系我们的客服，您也可以通过上面的方式进行在线报名，我们会及时给您回复电话，谢谢！
工程项目服务
有限元培训课程
产品设计培训课程
广州市东风东路808号 华宫大厦17层
电话：020-
培训QQ咨询：
项目QQ咨询：
点击图标或添加QQ：
深圳龙华新区和平路 28号福轩大厦1708
培训QQ咨询：
项目QQ咨询：
点击图标或添加QQ：
东莞市莞城区莞太路34号创意产业园11号楼5楼
培训QQ咨询：
项目QQ咨询：
点击图标或添加QQ：
佛山市禅城区弼塘东二街20号4楼（季华四路创意产业园东北门出口左侧）
培训QQ咨询：
项目QQ咨询：
点击图标或添加QQ：}