l,5,6
l,6,7
l,7,8
l,8,5
l,9,5
l,9,6
l,9,7
l,9,8
lsel,,loc,y,0,119 ! 按位置(Y)选择立柱(线)
cm,lvert,line ! 定义部件Vertical lines
lsel,,loc,y,120 !选择水平横梁
cm,lhoriz,line ! Horizontal lines
lsel,,loc,y,121,180 !选择顶梁
cm,lslope,line ! Sloping lines
lsel,all !选择所有线
et,1,188 !定义188梁单元
mp,ex,1,2e6 !定义弹性模量
mp,nuxy,1,0.3
mp,dens,1,7800
sectype,1,beam,i,beam !定义第一种截面为工字梁
secdata,6.535,6.535,8.06,.465,.465,.285 !定义梁截面尺寸
sectype,2,beam,i,column !定义第二种截面为工字梁
secdata,12,12,12.12,.605,.605,.39 !定义梁截面尺寸
sectype,3,beam,hrec,peak !定义第三种种截面为口字梁
secdata,6,6,.25,.25,.25,.25 !定义梁截面尺寸
save,frame,db
cmsel,,lslope ! Sloping lines
latt,1,,1,,100,,3 ! 顶梁的截面号为3,方向关键点为100
cmsel,,lvert
lsel,r,loc,x,-90 ! Left vertical lines
latt,1,,1,,102,,2 ! 左侧的垂直梁(线1、4)截面号为2,方向关键点为102
cmsel,,lvert
lsel,r,loc,x,90 ! Right vertical lines
latt,1,,1,,101,,2 ! 右侧的垂直梁(线2、3)截面号为2,方向关键点为101
cmsel,,lhoriz
lsel,u,loc,z,-60
lsel,u,loc,x,90 ! Left and front horizontal lines
latt,1,,1,,1,,1 ! 左前水平梁(线5、8)截面号为1,方向关键点为1
cmsel,,lhoriz
lsel,u,loc,z,60
lsel,u,loc,x,-90 ! Right and back horizontal lines
latt,1,,1,,3,,1 ! 右后水平梁(线6、7)截面号为1,方向关键点为3
lsel,all
lesize,all,20 ! 所有的线分割数都为20
lmesh,all !对所有的线划分网格
nsel,,loc,y,0 !选择y=0处的所有节点
d,all,all ! 约束全部自由度Fix bottom nodes
ksel,all !选择所有关键点 !
fk,9,fy,-10000 !给关键点9加-10000的Y方向的集中力
alls !选择所有实体及有限元节点和单元
/solu !进入SOLUTION
solve !求解
后处理
打开单元形状显示,查看变形及应力结果。
- 利用前面的轴承座模型进行扩展后处理的讲解
具体体会路径操作,单元表的使用,等值图的画法,标注、切片图等技术
- 梁壳实体单元混合使用分析(命令流)具体讲解
问题描述:广告牌承受风载荷的模拟
采用单元:
SOLID45
SHELL181
BEAM188
约束为两立柱的底部节点全部约束
风载全部加在壳单元的面上
材料均为钢材
以下为该例题的命令流
fini
/cle
/title, mult-element in a model
/filname,ele_test,1
/prep7
ET,1,SHELL181
ET,2,SOLID45
ET,3,BEAM188
R,1,0.05
r,2,0
SECTYPE,1,BEAM,CSOLID,se1
SECOFFSET, CENT
SECDATA,0.2,8,0,0,0,0,0,0,0,0
mp,ex,1,2e11
mp,dens,1,7800
mp,mu,1,0.3
mp,nuxy,1,0.3
rect,0,4,,3
block,,4,-.5,0,,0.5
block,,4,-.5,0,,-0.5
k,,0,-5,0
k,,4,-5,0
k,,5,-5,0
nummrg,kp
l,6,21
l,8,22
wpoff,,-0.1
wprot,,90
vsbw,all
wpoff,,,0.3
vsel,s,,,3,5,2
vsbw,all
alls
nummrg,all
lsel,s,,,5,7
lsel,a,,,21
LATT,1,2,3,,23,,1
asel,s,,,1,18,17
AATT,1,1,1,0
vsel,all
VATT,1,1,2,0
alls
WPCSYS,-1,0
lsel,s,loc,x,0.01,4-0.01
lesi,all,,,12,1
lsel,s,loc,y,0.01,2-0.01
lesi,all,,,8,2.5
lsel,a,loc,y,-.5-0.01,-5+0.01
lesi,all,,,10,3
alls
lesi,4,,,8,1/2.5
esize,0.05
asel,s,loc,z,0
asel,r,loc,y,-0.1,2
amesh,all
alls
lsel,s,loc,x,0
lsel,a,loc,x,4
lsel,r,loc,z,0
lsel,r,loc,y,-0.4,-5
lmesh,all
alls
vsel,all
vmesh,all
alls
asel,s,loc,y,0+0.01,2
sfa,all,1,pres,100
alls
nsel,s,loc,y,-5
d,all,all
alls
save
/solu
solve
fini
/post1
PLNSOL,S,EQV,0,1
ANCNTR,10,0.5
第四日 练习主题:巩固APDL基础、热应力及模态分析内容
1.弹簧建模
!This example demo modeling a ring using APDL
fini
/cle
/filname,ring,1
/prep7
*AFUN,DEG
TOT=12
N=TOT*3
et,1,45
et,2,200
keyopt,2,1,6
DIV=360/TOT
DIS=0.6/TOT
rid=1
rid2=0.1
csys,1
*DO,I,1,N+1
k,I,rid,(I-1)*DIV,(I-1)*DIS
*ENDDO
flst,3,n+1,3
*do,i,1,n+1
fitem,3,i
*enddo
bsplin,,p51x
ldiv,1,,,3
csys,0
wprot,,90
wpoff,ridpcirc,rid2,,0,360
asel,all
lsla,s
lesi,all,,,4
type,2
mshkey,0
amesh,all
type,1
extopt,esize,60
vdrag,1,,,,,,1
vdrag,6,,,,,,2
vdrag,11,,,,,,3