bet222体育投注

 找回密码
 立即注册
QQ登录 只需一步,快速开始
搜索
查看: 544|bet222体育投注: 0
打印 上一主题 下一主题

ansys边坡程序

[复制链接]
20主题50帖子258积分
机械厂工人 Rank: 6Rank: 6
UID
894
威望
0
机械币
86
活力
0
  1. 发消息
跳转到指定楼层
楼主
发表于 2016-11-7 09:34:36 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
         /prep7
et,1,plane42
keypot,1,3,2
mp,ex,1,20e9             !岩体参数
mp,prxy,1,0.22
mp,dens,1,2600
mp,ex,2,3e9                !滑体参数
mp,prxy,2,0.35
mp,dens,2,2140
mp,ex,3,2.7e9             !35%饱和度
mp,prxy,3,0.35
mp,dens,3,2180
mp,ex,4,2.4e9             !50%饱和度
mp,prxy,4,0.35
mp,dens,4,2210
mp,ex,5,2.1e9             !65%饱和度
mp,prxy,5,0.35
mp,dens,5,2240
mp,ex,6,1.8e9            !80%饱和度
mp,prxy,6,0.35
mp,dens,6,2270
mp,ex,7,1.5e9             !100%饱和度
mp,prxy,7,0.35
mp,dens,7,2300
save
k,1,0,0,0
k,2,0,505,0
k,3,562,505,0
k,,621,459,0
k,,658,433,0
k,,693,393,0
k,,802,313,0
k,,850,303,0
k,,892,273,0
k,,913,253,0
k,,930,243,0
k,,1034,233,0
k,,1186,228,0
k,,1216,223,0
k,,1232,208,0
*do,i,1,14
l,i,i+1
*enddo
k,16,660,387,0
k,,770,248,0
k,,827,230,0
k,,888,217,0
k,,930,216,0
k,,1034,213,0
k,,1186,209,0
l,4,16
l,16,17
l,17,18
l,18,19
l,19,20
l,20,21
l,21,22
l,22,15
nummrg,all
numcmp,all
k,23,1243,198,0
k,,612,380,0
k,,745,214,0
k,,884,181,0
k,,1242,174,0
l,3,24
l,24,25
l,25,26
l,26,27
l,27,23
l,23,15
nummrg,all
numcmp,all
k,28,1610,198,0
k,,1610,0,0
k,,1243,0,0
k,,884,0,0
k,,745,0,0
k,,612,0,0
l,6,16
l,7,17
l,8,18
l,9,19
l,11,20
l,13,22
l,16,24
l,17,25
l,19,26
l,24,33
l,25,32
l,26,31
l,27,30
l,1,33
l,30,31
l,31,32
l,32,33
l,23,28
l,28,29
l,29,30
save
al,4,5,15,29
al,29,6,16,30
al,7,17,30,31
al,8,18,31,32
al,9,10,19,32,33
al,11,12,20,21,33,34
al,13,14,22,34
al,3,15,35,23
al,35,16,24,36
al,36,17,18,37,25
al,19,20,21,22,26,27,28,37
al,1,2,23,38,42
al,38,39,24,45
al,39,40,25,44
al,40,41,26,43
al,46,47,48,27,41
allsel,all                    !选择所有元素
type,1
mat,1
mshape,1,2d
mshkey,0                 !设置划分方式,自由划分
esize,10,0
amesh,1,11,1
esize,,,
allsel,all
amesh,12,16,1
save
finish
!施加边界条件,并进行求解设置
/solu
eplot
save
nsubst,20           !指定荷载子步
nlgeom,on
nropt,full
pred,on
lnsrch,on
outres,all,last
!施加约束条件
nsel,s,loc,x,0
d,all,ux
nsel,s,loc,x,1610
d,all,ux
nsel,s,loc,y,0
d,all,all
allsel,all
acel,0,10,0
allsel,all
save
finish
!求解自重应力场
/solu
time,1
asel,s,,,1,7,1
esla,s
mpchg,2,all
allsel,all
solve
save
finish
!求解不同饱和度下的应力场
!35%饱和度
/solu
antype,,rest
time,2
asel,s,,,1,7,1
esla,s
mpchg,3,all
allsel,all
solve
save
finish
!50%饱和度
/solu
antype,,rest
time,3
asel,s,,,1,7,1
esla,s
mpchg,4,all
allsel,all
solve
save
finish
!65%饱和度
/solu
antype,,rest
time,4
asel,s,,,1,7,1
esla,s
mpchg,5,all
allsel,all
solve
save
finish
!80%饱和度
/solu
antype,,rest
time,5
asel,s,,,1,7,1
esla,s
mpchg,6,all
allsel,all
solve
save
finish
!100%饱和度
/solu
antype,,rest
time,6
asel,s,,,1,7,1
esla,s
mpchg,7,all
allsel,all
solve
save
finish
!初始状态下边坡稳定分析
/post1
set,1
plnsol,s,y,0,1          !初始状态下自重应力场
plnsol,u,y,0,1         !初始状态下位移场
                     plnsol,s,xy,0,1          !初始状态下剪应力
plnsol,epto,xy,0,1 !初始状态下剪应变
!荷载工况
/post1
lcdef,1,1,          !将第一个计算步定义为工况1
lcdef,2,6,          !将第六个计算步定义为工况2
lcase,2,              !读入荷载工况2
lcoper,sub,1,,,        !将荷载工况2减去工况1
!绘制位移、剪力、剪应变增量图
plnsol,u,y,0,1         !初始状态下位移场
plnsol,s,xy,0,1        !初始状态下剪应力
plnsol,epto,xy,0,1 !初始状态下剪应变
!采用mohr-column准则对饱水状态下的塑性区分布进行计算,C=50kPa,内摩擦角30度
/post1
asel,s,,,1,7,1
esla,s
etable,s1_1,s,1     !将所选单元的第1,3主应力定义为单元表
etable,s3_1,s,3
sadd,s1,s3_1,,-1.00,1
sadd,s3,s1_1,,-1.00,1
sadd,smf_over,s1,s3,1.00,-1.00
cos_fai=cos(30/180*3.1415926)       !计算内摩擦角的余弦
sin_fai=sin(30/180*3.1415926)
constant=2*50000*cos_fai
sadd,s1_s3,s1,s3,1,1
sexp,s1_s3,s1_s3,,1
sadd,smf_down,s1_s3,,sin_fai,,constant
sexp,s_m_f,smf_over,smf_down,1,-1
esel,all
/contour
pletab,s_m_f,avg  !绘制屈服单元的云图
分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏 分享淘帖 鲜花鲜花 板砖板砖
bet222体育投注 使用道具 举报
您需要登录后才可以回帖 登录 | 立即注册
本版积分规则
QQ|手机版|小黑屋|bet222体育投注 ( 渝ICP备10201771号-3 ) GMT+8, 2018-1-22 08:47 , Processed in 0.030956 second(s), 30 queries , Gzip On.
Powered by bet222体育投注 © 2013- Jixiecun.com
快速bet222体育投注 返回顶部 返回列表
bet222