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

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

s_elex(I,J)=1/6*(ss_nodex(I,J)+ss_nodex(I+1,J))+2/6*(ss_nodex(I,J+1)+ss_nodex(I+1,J+1))
 *ENDDO
 *ENDDO 

!2.2.5计算得到单元的ss_eleZ
*DO,I, 1,5, 1
 *DO,J, 1,16, 1 
ss_elez(I,J)=1/6*(ss_nodez(I,J)+ss_nodez(I+1,J))+2/6*(ss_nodez(I,J+1)+ss_nodez(I+1,J+1))
 *ENDDO
 *ENDDO 
!2.2.6 计算得到单元的sigma1=ss_elex+ss_elez
*DO,I, 1,5, 1
 *DO,J, 1,16, 1 
*VOPER, sigma1(I,J),ss_elex(I,J),add,ss_elez(I,J)
*VOPER, sigma3(I,J),ss_elex(I,J),MULT,1
 *ENDDO
 *ENDDO 


!2.2.7 求得本次计算得到的弹性模量值excal(I,J)
*DO,I, 1,5, 1
 *DO,J, 1,16, 1 
!excal(I,J)= (45.481+ 1.105*(abs(sigma1(I,J)))**1.0)*(abs((sigma3(I,J)))**1.0)+183368.8 !!!unit kpa 参数组1
!excal(I,J)=(8617.70+ 3064.0*(abs(sigma1(I,J)))**0.421)*(abs((sigma3(I,J)))**0.33)+30000 !!!unit kpa 参数组2
excal(I,J)=(10341.10+5021.90*(abs(sigma1(I,J)))**0.4)*(abs((sigma3(I,J)))**0.33)+38000 !!!unit kpa 参数组3
 *ENDDO
 *ENDDO 

!!!E=(45.481+1.105*σ1^1.0)*σ3^1.0+183368.8 kPa 参数组1
!!!E=(8617.7+3064.0*σ1^0.421)*σ3^0.33+30000 kPa 参数组2
!!!E=(10341.1+5021.9*σ1^0.4)*σ3^0.33+38000 kPa 参数组3

!2.2.8两次的弹性模量值相减,求RMS(根方差),并用RMS
*DO,I, 1,5, 1
 *DO,J, 1,16, 1 
*VOPER,exsub(I,J),ex(I,J),SUB,excal(I,J)
 *ENDDO
 *ENDDO 

*VSCFUN,exrms,RMS,exsub
*IF, exrms,ge,eps, THEN
 *MFUN, ex,copy,excal
/PREP7 
!2.2.9 重新定义材料属性,循环
*DO,I, 1,2, 1
 *DO,J, 1,16, 1
MP, EX,16*(I-1)+J, EX(I,J) 
MP,PRXY,16*(I-1)+J, PRXY2
MP,DENS,16*(I-1)+J, 0
 *ENDDO
*ENDDO
*DO,I, 3,5, 1
 *DO,J, 1,16, 1
MP, EX,16*(I-1)+J, EX(I,J) 
MP,PRXY,16*(I-1)+J, PRXY3
MP,DENS,16*(I-1)+J, 0
 *ENDDO
*ENDDO
!2.2.10 修改单元材料,循环
*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
nnum=nnum+1 
*IF,nnum,ge,100, THEN
 *MSG,UI
 Sorry, the step number is greater than 100, Exit! BYE......
 *EXIT
*ENDIF 
!2.2.11 循环求解
FINISH
/SOLU
ANTYPE,0
NLGEOM,0
NROPT,AUTO, ,
EQSLV, , ,0,
PIVCHECK,1
SOLVE
FINISH

*ELSE
 *MSG,UI
 Ok,Congratulation,XiaoMei! EX has been found!! Very Good!!!
 *EXIT
*ENDIF
*ENDDO
!2.2.12 后处理
/post1
!!未完待续,继续执行post.ans文件
!3.0后处理 一共输出12项内容,内容12涵盖了内容1-11

/PREP7 
!3.1 定义参数数组
*DIM,LABEX,CHAR,1,16 
LABEX(1, 1)=''EX1''
LABEX(1, 2)=''EX2''
LABEX(1, 3)=''EX3''
LABEX(1, 4)=''EX4''
LABEX(1, 5)=''EX5''
LABEX(1, 6)=''EX6''
LABEX(1, 7)=''EX7''
LABEX(1, 8)=''EX8''
LABEX(1, 9)=''EX9''
LABEX(1,10)=''EX10''
LABEX(1,11)=''EX11''
LABEX(1,12)=''EX12''
LABEX(1,13)=''EX13''
LABEX(1,14)=''EX14''
LABEX(1,15)=''EX15''
LABEX(1,16)=''EX16''
*DIM,LABSS,CHAR,6,17
LABSS(1, 1)=''SSNX 1''
LABSS(1, 2)=''SSNX 2''
LABSS(1, 3)=''SSNX 3''
LABSS(1, 4)=''SSNX 4''
LABSS(1, 5)=''SSNX 5''
LABSS(1, 6)=''SSNX 6''
LABSS(1, 7)=''SSNX 7''
LABSS(1, 8)=''SSNX 8''
LABSS(1, 9)=''SSNX 9''
LABSS(1,10)=''SSNX10''
LABSS(1,11)=''SSNX11''
LABSS(1,12)=''SSNX12''
LABSS(1,13)=''SSNX13''
LABSS(1,14)=''SSNX14''
LABSS(1,15)=''SSNX
|<< << < 1 2 3 4 5 6 7 8 9 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:道路  基层  垫层  模量  应力  关系
相关资讯
热点资讯
推荐资讯

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