采用ANSYS软件建立混凝土试件三维随机骨料模型

ANSYS模拟仿真软件apdl语言命令流-——混凝土试件的三维随机骨料建模FINISH/clear,nostart!-----------三维混凝土150试件，单轴拉伸模拟!-------映射剖分，以最小骨料半径剖分网格，判断属性，不细分!=================================================xmin=1$xmax=149$ymin=1$ymax=149$zmin=1$zmax=149!定义坐标范围anum=120$bnum=230$cnum=510!定义三种组分骨料个数num=anum+bnum+cnumrmin=8$rmax=10$brmin=5$brmax=8$crmin=2.5$crmax=5!定义骨料半径范围*dim,agv,array,num,4!存放骨料位置及半径的数组cum=0icum=0*do,i,1,20000*if,cum,eq,num,then*exit*endif*if,icum,eq,anum,thenrmin=brmin$rmax=brmax*elseif,icum,eq,(anum+bnum),thenrmin=crmin$rmax=crmax*endifx=rand(xmin,xmax)$y=rand(ymin,ymax)$z=rand(zmin,zmax)$r=rand(rmin,rmax)*if,x-r,gt,xmin,and,x+r,lt,xmax,then*if,y-r,gt,ymin,and,y+r,lt,ymax,then*if,z-r,gt,zmin,and,z+r,lt,zmax,then*if,cum,eq,0,thencum=cum+1icum=icum+1agv(cum,1)=x$agv(cum,2)=y$agv(cum,3)=z$agv(cum,4)=r*elsesum=0*do,j,1,cumdist=sqrt((agv(j,1)-x)**2+(agv(j,2)-y)**2+(agv(j,3)-z)**2)!1.05为骨料影响范围系数*if,dist,lt,1.1*(agv(j,4)+r),then*exit*elsesum=sum+1*endif*enddo*if,sum,eq,cum,thencum=cum+1icum=icum+1agv(cum,1)=x$agv(cum,2)=y$agv(cum,3)=z$agv(cum,4)=r*endif*endif*endif*endif*endif*enddo!/TRLCY,volu,1,221!/TRLCY,area,1,221!=====================先划分网格，然后每个单元分给一个材料号====================!=====================按weibull概率分布赋予属性====================/prep7ET,1,SOLID45BLOCK,0,150,0,150,0,150!生成投放区域LSEL,ALLLESIZE,ALL,,,100VSEL,ALLmshape,0,3d$mshkey,1!自由网格，划分砂浆vmesh,allEMODIF,ALL,MAT,2,GULIAO=0JIEDIANX=0$JIEDIANY=0$JIEDIANZ=0!判断骨料砂浆属性!和所有的骨料比较*DO,i,1,NUMALLSEL,ALLNSEL,S,LOC,X,agv(i,1)-agv(i,4),agv(i,1)+agv(i,4)NSEL,R,LOC,Y,agv(i,2)-agv(i,4),agv(i,2)+agv(i,4)NSEL,R,LOC,Z,agv(i,3)-agv(i,4),agv(i,3)+agv(i,4)ESLN,S!选择集*GET,ENUM,ELEM,0,COUNT*GET,EMIN1,ELEM,0,NUM,MIN*DIM,ELEM,ARRAY,ENUMES=EMIN1-1*DO,J,1,ENUMES=ELNEXT(ES)ELEM(J)=ES!单元的编号*ENDDO!遍历单元*DO,K,1,ENUMGULIAO=0*IF,ELMIQR(ELEM(K),-1),EQ,2,THEN!!循环比较8个节点*DO,L,1,6JIEDIANX=NX(NELEM(ELEM(K),L))JIEDIANY=NY(NELEM(ELEM(K),L))JIEDIANZ=NZ(NELEM(ELEM(K),L))distA=sqrt((agv(i,1)-JIEDIANX)**2+(agv(i,2)-JIEDIANY)**2+(agv(i,3)-JIEDIANZ)**2)*IF,distA,LT,AGV(i,4),THENGULIAO=GULIAO+1*ENDIF*ENDDO*IF,GULIAO,EQ,6,THENEMODIF,ELEM(K),MAT,1*ELSEIF,GULIAO,LT,6,AND,GULIAO,GT,0,THENEMODIF,ELEM(K),MAT,3*ENDIF*ENDIF*ENDDO*SET,ELEM(1),*ENDDO!==========================建模完成====================