用ansys做非稳定渗流及边坡稳定分析第8、10页均为动画表示的渗流场随时步变化的情况-下载后就能看到了.pptVIP

 ANSYS使用经验交流　　　　 作者邮箱：erep@ 汇报安排 一、简介使用ANSYS对课题的研究情况 基本理论 基本理论 非稳定渗流程序结构框图 非稳定渗流程序结构框图 库水位变动时的渗流场模拟 库水位变动时的渗流场模拟 降雨条件下渗流场模拟 降雨条件下渗流场模拟 强度折减有限元进行边坡稳定分析——塑性区图 降雨条件下渗流场模拟 边坡稳定分析结果 二、浅谈ANSYS学习方法 二、浅谈ANSYS学习方法—心得体会 二、浅谈ANSYS学习方法 * * 一、简介使用ANSYS对课题的研究情况 二、浅谈ANSYS学习方法 三、相互交流学习 基于ANSYS,开发出了土石坝渗流与边坡稳定分析程序. 1、库水位变动时的土石坝渗流场模拟； 2、降雨条件下土石坝渗流场模拟； 3、以上两种非稳定渗流条件下土石坝边坡稳定分析。 非稳定饱和—非饱和渗流控制方程 非饱和土抗剪强度 ！调整单元渗透系数的语句 /PREP7 MATNUM=2 *DO,I,1,E_MAX ePRESS=BBB(i)-ccc(i) *IF,ePRESS,GT,0,THEN MPCHG,1,I *ELSEIF,ePRESS,LT,-13.812,THEN MPCHG,2,I *ELSE MP,KXX,MATNUM+1,kxxbh*hsl(ePRESS) MPCHG,MATNUM+1,I MATNUM=MATNUM+1 *ENDIF *ENDDO 上面模拟的关键点在于：逸出边界的确定 逸出边界需要在定水头边界转化为不透水边界的动态变化过程中确定。 *get,q_day1,fsum,0,item,heat *IF,q_day1,lt,0,THEN cmsel,s,nn2 nsel,a,,,node1 cm,nn2,node *ENDIF cmsel,s,nn1 node1=ndnext(node1) *ENDDO /post1 cmsel,s,nn1 *GET,LN_MAX,NODE,,NUM,MAX *GET,LN_MIN,NODE,,NUM,MIN *GET,Num1,NODE,0,count node1=LN_MIN *DO,I,1,Num1 a1=ny(node1) nsel,s,,,node1 prrsol,heat fsum,heat 上面模拟的关键点在于：降雨边界的确定 降雨边界在定流量边界与定水头边界之间相互转化。 *endif cmsel,s,nn4 node1=ndnext(node1) *ENDDO !退出条件 cmsel,s,nn4 *GET,Num2,NODE,0,count *if,num1,eq,num2,exit allsel cmsel,s,nn4 *GET,LN_MAX,NODE,,NUM,MAX *GET,LN_MIN,NODE,,NUM,MIN! *GET,Num1,NODE,0,count node1=LN_MIN *DO,I,1,Num1 *if,hc1%k%(node1),ge,0,then cmsel,s,nn1 nsel,a,,,node1 cm,nn1,node 二分法编制的宏文件 zhuangtai=_status *if,zhuangtai,gt,1,then fmax=fs *else fmin=fs *endif jingdu=abs(fss-fs) *if,jingdu,le,0.01,then fs=(fmin+fmax)/2 *exit *endif *enddo *end!稳定宏结束 *creat,wending,mac allsel fmin=0 fmax=2.5 fs=0 /prep7 /solu *do,j,1,30 f ss=fs fs=(f