有限元分析基础教程-ANSYS算例 曾攀
force_2=-1000 ! 设置外力F2的数值 !设置参数-----end---
/PREP7 !进入前处理
ANTYPE,STATIC !设置分析类型(静力分析)
ET,1,LINK1 ! 设置单元类型(第1号:杆单元LINK1) R,1,bar_a !设置实常数(第1号:横截面积为bar_a) MP,EX,1, bar_el !设置材料常数(弹性模量)(第1号:bar_el)
N,1 $N,2,, bar_1 $N,3,, bar_2 $N,4,, bar_L !根据坐标建立1,2,3,4号节点 E,1,2 !建立一个单元(通过连接第1号节点和第2号节点) EGEN,3,1,1 !再重复产生3个单元,每次的节点号都加1
D,1,ALL,,,4,3 !施加完全约束,节点起始号为1,终止号为4,增量为3 F,2,FY, force_1 !施加节点力(对2节点,FY= force_1) F,3,FY, force_2 !施加节点力(对3节点,FY= force_2) FINISH !结束以上流程 /SOLU !进入求解状态 OUTPR,BASIC,1 !输出方式的设置 SOLVE !求解
FINISH !结束以上流程 /POST1 !进入后处理
/ESHAPE,0.5 !对于线单元(如杆、梁)按实体效果进行显示(以0.5倍实常数的比例) !将变形图输出到JPEG
JPGPRF,1500,300,1 !设置输出文件格式 /SHOW,JPEG ! 设置输出文件格式
PLNSOL,U,Y,0,1,0 !显示变形图,并将输出file000.jpg ETABLE,bar_I,SMISC,1 !定义线型单元I节点的轴力 ETABLE,bar_J,SMISC,1 !定义线型单元J节点的轴力 PLLS, bar_I, bar_J !画出线型单元的受力图
PRESOL,SMISC,1 ! 显示单元的轴向力(可以看到支反力),并将输出file001.jpg !(续)有时最后一个图形文件将在退出ansys才能得到 PRRSOL !打印单元的支反力结果 FINISH !结束以上流程
!%%%%%%%%parameterized log file: bar_para.log %%%% end %%%%%%
58