!2.2 工字形截面外伸梁的平面弯曲
! 本程序来源于邢静忠等编著的《ANSYS分析实例与工程应用》，机械工业出版社（2004年）
Finish
/CLEAR, NOSTART
/PREP7
/FileName, EX2.2
/TITLE, EX2.2, BEAM STRESSES AND DEFLECTIONS
C*** STR. OF MATL., TIMOSHENKO, PART 1, 3RD ED., PAGE 98, PROB. 4
! (1) 定义求解类型，实常数和材料参数
ANTYPE,STATIC
ET,1,BEAM3
KEYOPT,1,9,9
R,1,35e-4,4.5E-5,0.20
MP,EX,1,210E9
MP,PRXY,1,0.3
! (2) 定义节点位置
N,1
N,5,4
FILL
! (3) 定义单元连接关系
E,1,2
EGEN,4,1,1
! (4) 定义位移约束和荷载条件
D,2,UX,,,,,UY
D,4,UY
SFBEAM,1,1,PRES,5000
SFBEAM,4,1,PRES,5000
FINISH
! (5) 进入求解模块，定义求解选项并求解
/SOLU
OUTPR,BASIC,1
SOLVE
FINISH
! (6) 进入后处理模块列表显示节点位移计算结果和图形显示变形图
/POST1
SET,1,1
PRNSOL,U,COMP
PRNSOL,ROT,COMP
PLDISP,1
! (7) 提取指定位置的节点和单元，定义单元表并获取弯曲应力
MID_NODE = NODE (2,,, )
*GET,DISP,NODE,MID_NODE,U,Y
MID_ELM = ENEARN (MID_NODE)
ETABLE,STRS,LS,3
*GET,STRSS,ELEM,MID_ELM,ETAB,STRS
! (8) 将计算结果写入到文本文件
*DIM,LABEL,CHAR,2
*DIM,VALUE,,2,3
LABEL(1) = 'STRS_MPa','DEF_mm'
*VFILL,VALUE(1,1),DATA,-5.5556,0.1323
*VFILL,VALUE(1,2),DATA,STRSS/1e6,DISP*1000
*VFILL,VALUE(1,3),DATA,ABS(STRSS/1e6/5.5556) ,ABS(DISP/0.1323e-3)
/OUT, EX2_2, out
/COM
/COM,------------------- EX2.2 RESULTS COMPARISON ---------------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)
/COM,---------------------------------------------------------------
/OUT
FINISH
*LIST, EX2_2, out
