求ANSYS悬臂梁的实例分析

如题所述

给你个例子吧
FINISH
/CLEAR, NOSTART
/FileName, EX2.6
! (1) 定义工程选项,分析类型为静力(STATIC)分析,单元类型和实常数,材料参数等
/PREP7
/TITLE, EX2.6(1), BENDING OF A TAPERED PLATE (BEAM) WITH SHELL63 ELEMENT
! INTROD. TO STRESS ANALYSIS, HARRIS, 1ST PRINTING, PAGE 114, PROB. 61
ANTYPE,STATIC
ET,1,SHELL63,2
R,1,12.7
MP,EX,1,207E3
MP,NUXY,1,0
! (2) 通过节点和单元建模
N,1
N,8,508,-38.1
FILL
N,11
N,18, 508,38.1
FILL
E,1,2,12
E,2,3,12
E,13,12,3
E,3,4,14
E,14,13,3
EGEN,3,2,2,5
! (3) 定义耦合自由度,以保证变形对称
CP,1,UZ,2,12
CP,2,ROTY,2,12
CPSGEN,6,1,1,2
OUTPR,ALL,ALL
D,8,ALL,,,18,10
D, ALL,ROTX,0
F,1,FZ,-44.483
FINISH
! (4) 进入求解模块,开始求解
/SOLU
SOLVE
FINISH
! (5) 进入一般后处理模块,绘制变形图定义单元表,提取应力计算结果和节点位移值
/POST1
/VIEW,1,1,1,1
PLDISP,2
ETABLE,STRS,S,1
ESORT,STRS
*GET,SMAX1,SORT,,MAX
PRNSOL,DOF
LFT_NODE = NODE(0,0,0)
*GET,DEFL1,NODE,LFT_NODE,U,Z
! (6) 定义数组变量,并将理论值和计算值保存到变量
*DIM,LABEL,CHAR,2,2
*DIM,VALUE1,,2,3
LABEL(1,1) = 'DEFLECTI','MX_PRIN_'
LABEL(1,2) = 'ON(mm) ','STRS(MPa)'
*VFILL,VALUE1(1,1),DATA,-1.083,11.032
*VFILL,VALUE1(1,2),DATA,DEFL1,SMAX1
*VFILL,VALUE1(1,3),DATA,ABS(DEFL1/1.083),ABS(SMAX1/11.032)
SAVE, TABLE_1
FINISH
! (7) 利用三维梁单元BEAM44,重新计算前面的问题
/CLEAR, NOSTART
/PREP7
/TITLE, EX2.6(2), BENDING OF A TAPERED PLATE WITH TAPERED BEAM ELEMENTS BEAM44
! (8) 定义分析类型,单元类型和实常数
ANTYPE,STATIC
ET,1,BEAM44
*DO,I,1,10
R,I,1,1,(I-1)*76.2*12.7*12.7*12.7/120,6.35,1
RMORE,1,1,I*76.2*12.7*12.7*12.7/120,6.35,1
*ENDDO
RMOD,1,3,76.2*12.7*12.7*12.7/120/2
RLIST
MP,EX,1,207E3
MP,GXY,1,207E3/2.6
! (9) 定义节点位置
N,1
N,11,508
FILL
N,12,,,10
NGEN,10,1,12
! (10) 定义梁单元
E,1,2,12
EGEN,10,1,1,,,,,1
! (11) 施加位移约束条件
D,11,ALL
D,1,UY,,,10,,ROTX,ROTZ
F,1,FZ,-44.483
OUTPR,ALL,ALL
FINISH
! (12) 在求解模块中求解
/SOLU
SOLVE
FINISH
! (13) 在后处理模块中,绘制变形图,定义单元表,提取应力计算结果
/POST1
/VIEW,1,1,1,1
PLDISP,2
ETABLE,STRS,NMISC,1
ESORT,STRS
*GET,SMAX2,SORT,,MAX
PRNSOL,DOF
LFT_NODE = NODE (0,0,0)
*GET,DEFL2,NODE,LFT_NODE,U,Z
! (14) 定义数组,将计算结果及其与理论结果的比较输出到文件
*DIM,LABEL,CHAR,2,2
*DIM,VALUE2,,2,3
LABEL(1,1) = 'DEFLECTI','MXPRIN_'
LABEL(1,2) = 'ON(mm) ','STRS(MPa)'
*VFILL,VALUE2(1,1),DATA,-1.083,11.032
*VFILL,VALUE2(1,2),DATA,DEFL2,SMAX2
*VFILL,VALUE2(1,3),DATA,ABS(DEFL2/1.083),ABS(SMAX2/11.032)
SAVE, TABLE_2
温馨提示:答案为网友推荐,仅供参考
第1个回答  2012-06-17
一工字钢悬臂梁长10m,截面高0.42m,截面面积0.0072m2,截面惯性矩2.108×10-4m4,材料弹性模量为2×1011Pa, 泊松比0.3,在集中力P=10000N作用下求该梁各点的挠度。采用beam3单元(分10份)
第2个回答  2011-06-01
恩,这个我也想学学