道路的基层、垫层模量与应力之间的关系

发布日期:[08-12-04 12:42:17] 浏览人次:[]

320*K+20*(J-1)+I+1,320*(K-1)+20*(J-1)+I, 320*(K-1)+20*J+I, 320*(K-1)+20*J+I+1, 320*(K-1)+20*(J-1)+I+1
 
*ENDDO
 *ENDDO
*ENDDO
!修改单元材料
*DO,I, 1,192, 1
EMODIF,I,mat,99
*ENDDO
*DO,I, 1153,2112, 1
EMODIF,I,mat,100
*ENDDO
*DO,K, 1,5, 1
 *DO,J, 1,16, 1
 *DO,I, 1,12, 1
EMODIF,192*K+12*(J-1)+I,mat,16*(K-1)+J
 *ENDDO
 *ENDDO
*ENDDO

 ESEL, s, ELEM, , 1, 12, 1, 

*DO,I, 2,11, 1 !2-11 layers
 ESEL, a, ELEM, ,192*(I-1)+1, 192*(I-1)+12, 1, 
*ENDDO
 
EPLOT 
/WAIT, 3
EMODIF,ALL,TYPE,2, 
/WAIT, 3
NUMSTR, node, 4001
EMID,ADD,ALL
ALLSEL,ALL
 

!1.4边界条件
/SOLU 
!固结底面节点
csys,0
NSEL,S,LOC,Y,-3.05 
NPLOT 
D,all, , , , , ,ALL, , , , , 
ALLSEL,ALL 
!固结圆周面节点
csys,5
NSEL,S,LOC,x,3.10 
NPLOT 
D,all, , , , , ,all, , , , , 
ALLSEL,ALL 
!约束中心线节点,使之成为对称面,关于X轴和Z轴。
csys,0
 *DO,I, 1,12, 1
D, I, UX, 0, , , , UZ, 
*ENDDO 
csys,0
!约束X=0平面节点,使之成为对称面
*DO,J, 1,11, 1 
*DO,I, 2,16, 1
D, 320*(J-1)+20*(I-1)+1, UX, 0
 *ENDDO 
*ENDDO 
csys,0
!约束Z=0平面节点,使之成为对称面
*DO,J, 1,11, 1 
*DO,I, 2,16, 1
D, 320*(J-1)+20*(I-1)+13, UZ, 0
 *ENDDO 
*ENDDO 


/SOLU 
!4.3施加荷载,在半径10.65cm圆面上
*DO,I, 1,12, 1
SFE,I,6,PRES, ,700, , , 
*ENDDO

!1.5.求解
FINISH 
/SOLU 
ANTYPE,0
NLGEOM,0
NROPT,AUTO, , 
EQSLV, , ,0,
PIVCHECK,1 
SOLVE 
FINISH 
!!未完待续,继续执行iteration.ans文件
!2.0 跌代计算,以获得弹性模量的收敛值
NMAX=100 !指定最大循环次数100
nnum=1 !循环次数计算器nnum,1为初始值
eps=0.5e2 !跌代收敛精度 0.05mpa
!2.1 说明

!2.2 循环跌代计算
*DO,m, 1,NMAX, 1

/PREP7 
!2.2.1 定义参数数组
*DIM,ss_nodex,ARRAY,6,17 !存放ss_nodex 节点的径向应力,6个高度线(自上至下),17个节点(自内至外)
*DIM,ss_nodez,ARRAY,6,17 !存放ss_nodez 节点的垂直向应力,6个高度线(自上至下),17个节点(自内至外)
*DIM,ss_elex,ARRAY,5,16 !存放ss_elex 单元的径向应力,位置:(2/3单元边长处,1/2单元高度处),5层(自上至下),16圈(自内至外)
*DIM,ss_elez,ARRAY,5,16 !存放ss_elez 单元的垂直向应力,位置:(2/3单元边长处,1/2单元高度处),5层(自上至下),16圈(自内至外)
*DIM,sigma1,ARRAY,5,16 !存放sigma1 5x16 三轴试验中的单元的σ1=ss_elex +ss_elez 
*DIM,sigma3,ARRAY,5,16 !存放sigma3 5x16 三轴试验中的单元的σ3=ss_elex
*DIM,excal,ARRAY,5,16 !存放excal 5x16 本次计算得到的弹性模量值excal
*DIM,exsub,ARRAY,5,16 !存放exsub 5x16 EX-excal=EXSUB(kPa) 

/POST1 
RSYS,5
!2.2.2 获得水平应力ss_nodeX
 *DO,I, 1,6, 1
*GET, ss_nodex(I,1), NODE, I+1, S, X !ss_nodex X of NO.0 round''s nodes
 *ENDDO
*DO,I, 1,6, 1
 *DO,J, 2,17, 1 
*GET, ss_nodex(I,J), NODE, 320*I+20*(J-1)+7, S, X !ss_nodex X of NO.1-16 round''s nodes 
*ENDDO
 *ENDDO 
/POST1 
RSYS,5
!2.2.3获得垂直应力ss_nodeZ
 *DO,I, 1,6, 1
*GET, ss_nodez(I,1), NODE, I+1, S, Z !ss_nodez X of NO.0 round''s nodes
 *ENDDO
*DO,I, 1,6, 1
 *DO,J, 2,17, 1 
*GET, ss_nodez(I,J), NODE, 320*I+20*(J-1)+7, S, Z !ss_nodez X of NO.1-16 round''s nodes 
*ENDDO
 *ENDDO 
/POST1 
RSYS,0

!2.2.4计算得到单元的ss_eleX
*DO,I, 1,5, 1
 *DO,J, 1,16, 1 
s
|<< << < 1 2 3 4 5 6 7 8 9 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:道路  基层  垫层  模量  应力  关系
相关资讯
热点资讯
推荐资讯

关于我们 | 站点导航 | 使用帮助 | 友情链接 | 广告服务 | 免责声明 | 新手上路
设为首页 | 加入收藏 | 在线留言 | 马棚网QQ群:{92562572}{102901272}{333259257} | 交流QQ: 客户服务 客户服务 客户服务