已关闭

tozhou

已关闭

inkey

哈哈             做出这类随机性的好贴

FreddyMusic

有意思！

1。可生成球形颗粒、椭圆形颗粒、片状颗粒、不规则多面体颗粒和弥散型微观结构等不同的随机微观结构模型。

现有命令流只是圆 ？

2. 颗粒的最小间距（本文中设为MIN/10）的调整则可根据实际情况模拟出材料的不同微观结构形态。上述参数均可以根据材料的实际情况选定，以最大限度地模拟材料的实际情况。
控制参数确定后，各颗粒（包括夹杂相颗粒）在基体中的位置依给定的随机方法生成，不小于给定的最小间距。位置、直径和间距的随机分布规律可为平均分布、正态分布等，依实际情况选定。

我觉得这是价值所在，根据试验结论和仿真联系在了一起。

benjackxu

事实上现在已经有人做出多边形随机生长方式建立随机骨料的投放算法，见孙立国,杜成斌,戴春霞. 大体积混凝土随机骨料数值模拟[J].河海大学学报(自然科学版),2005, 33(3):291-295.
tozhou的方法是这样的吗：
取已知级配骨料，假定骨料颗粒为球状，其横断面颗粒直径可借助于Fuller曲线 确定，由该级配曲线浇筑的混凝土可产生优化的结构密度和强度。骨料形状和大小采用随机模拟方法，分布依据“互不靠”规则？

ilxy

还是用英语说吧。
the  strength of concrete depends on strength of cement paste and interfacial transition
zone between cement paste and aggregate. if you want to predict concrete strength well,
you must pay more attention to interfacial transition zone.
generally,the water cement ratio in interfacial transition zone is much bigger than cement
block paste. so the hydration process in this zone must be simulated carefully. in fact,
compared with simulation of cement hydration, that of aggregate is much easier.

ilxy

function huayuan
很早以前写的一个随机产生圆的程序，
输入参数的说明如下：
depth画图范围，
radi为圆的半径（数组）
number对应各种半径的圆的数量（数组）

clear
clc
format compact
%depth画图范围，radi为圆的半径，number对应各种半径的圆的数量
depth=10
radi=0.8:-0.1:0.1
number=8:-1:1
%产生画第一批最大的圆
xxold=[]
yyold=[]
[x,y]=huacirstart(depth,number(1),radi(1),xxold,yyold)
xxold=x
yyold=y
%产生第二批圆及以后的圆
for i=2:length(radi)
    radrawn=radi(1:i-1)
    numberdrawn=number(1:i-1)
    [x,y]=huacir(depth,number(i),radi(i),xxold,yyold,radrawn,numberdrawn)
    xxold=x
    yyold=y
end
%改变画圆的方法
aa=[]
for i=1:length(radi)
aa=[aa ones(1,number(i))*radi(i)]
end
for i=1:length(aa)
rectangle('position',[x(i)-aa(i),y(i)-aa(i),2*aa(i),2*aa(i)],'curvature',[1,1],'facecolor','b')
axis([0,depth,0,depth])
grid
hold on
end
%开始画第一批最大的圆的方法
function [xxnew,yynew]=huacirstart(depth,number,radi,xxold,yyold)
xxlen=length(xxold)
yylen=length(yyold)
xxold(xxlen+1)=radi+(depth-2*radi)*rand
yyold(yylen+1)=radi+(depth-2*radi)*rand
for i=2:number
     while 1
        xxold(xxlen+i)=radi+(depth-2*radi)*rand
        yyold(yylen+i)=radi+(depth-2*radi)*rand
        for j=1:xxlen+i-1
            distance(j)=((xxold(xxlen+i)-xxold(j))^2+(yyold(yylen+i)-yyold(j))^2)^0.5
        end
        judg=distance-2*radi
       if judg>0
       break
        end
    end
end
xxnew=xxold
yynew=yyold

%画第二批及以后的圆的方法
function [xxnew,yynew]=huacir(depth,number,radi,xxold,yyold,radrawn,numberdrawn)
xxlen=length(xxold)
yylen=length(yyold)
for i=1:number
     while 1
        xxold(xxlen+i)=radi+(depth-2*radi)*rand
        yyold(yylen+i)=radi+(depth-2*radi)*rand
        for j=1:xxlen+i-1
            distance(j)=((xxold(xxlen+i)-xxold(j))^2+(yyold(yylen+i)-yyold(j))^2)^0.5
        end
        %产生对比的矩阵
          aa=[]
          for I=1:length(radrawn)
          aa=[aa ones(1,numberdrawn(I))*radrawn(I)]
           end
          aa=[aa ones(1,i-1)*radi]
          aa=aa+radi
         judg=distance-aa
       if all(judg>0)
       break
        end
    end
end
xxnew=xxold
yynew=yyold

[ 本帖最后由 ilxy 于 2007-3-24 12:28 编辑 ]

wishangtian

个人觉得还是ilxy写的命令流好些。

tozhou 的命令流是GUI操作产生的log，看起来很费劲的。

xiaopohair

不错，可是如果我考虑增强颗粒密度很大的情况时候怎么办？颗粒之间假如有相互粘接甚至嵌入的情况，怎么定义接触？且没有基体相，只有相互搭接的骨架。

sybhit

有几个问题：
夹杂项的体积分数在模型中通过rate参数控制吗？如果是，体积分数怎么好象是大于0.4。
能够控制各种粒子直径的比例吗？
如何控制呢
我是新手，多多帮忙

sybhit

你的程序如何zaiansys 中实现亚，打下

sybhit

我是刚刚接触这方面的知识，
你的模型我有几个问题：
夹杂向体积控制问题，好象你的模型里面要大于0.4
能否生成已知直径的夹杂呢
，例如
半径为0.1与0.2的圆所占体积比为：1：2

yuxiaoyuh1138

这个好像
是matlab程序吧

sybhit

我想知道，flag是不是算作一个调用函数呢，我把命令中添加了不同粒径材料的体积比。但是运行不了，看看有什么毛病
/PREP7  
/COM,  Structural   
!*  
!*  
ET,1,PLANE42
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,4.2e11
MPDATA,PRXY,1,,0.17
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,2.07e11
MPDATA,PRXY,2,,0.32
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,3,,3.5e11
MPDATA,PRXY,3,,0.22
*DIM,MM,ARRAY,4000,3
*SET,XL,0.0
*SET,YB,0.0
*SET,XR,0.001   
*SET,YT,0.001   
*SET,B,XR-XL
*SET,H,YT-YB
*SET,TT,0.0
*SET,TOTALAREA,B*H  
*SET,AREA,0
*SET,AREA1,0
*SET,AREA2,0
*SET,RATE,0.4
*SET,RATE1,0.05
*SET,RATE2,0.1
*SET,FLAG,1
*SET,FLAGE,1
*SET,FLAGEE,1
*SET,FLAG1,0
*SET,i,1
*SET,Rmin,0.007*B
*SET,Rmax,0.049*B
BLC4,,,B,H  
*SET,MM(1,1),RAND(Rmin,Rmax)
*SET,d,Rmin/10  
*SET,MM(1,2),RAND(XL+MM(1,1)+d,XR-MM(1,1)-d)
*SET,MM(1,3),RAND(YB+MM(1,1)+d,YT-MM(1,1)-d)
CYL4,MM(1,2),MM(1,3),7E-6
*SET,AREA,3.14159265*MM(1,1)*MM(1,1)
*SET,MM(1,1),RAND(Rmin,Rmax)
*SET,MM(2,2),RAND(XL+MM(2,1)+d,XR-MM(2,1)-d)
*SET,MM(2,3),RAND(YB+MM(2,1)+d,YT-MM(2,1)-d)
CYL4,MM(2,2),MM(2,3),1E-6
*SET,AREA1,3.14159265*MM(2,1)*MM(2,1)
*SET,MM(3,1),RAND(Rmin,Rmax)
*SET,MM(3,2),RAND(XL+MM(2,1)+d,XR-MM(2,1)-d)
*SET,MM(3,3),RAND(YB+MM(2,1)+d,YT-MM(2,1)-d)
CYL4,MM(3,2),MM(3,3),0.5E-6
*SET,AREA2,3.14159265*MM(3,1)*MM(3,1)
*DOWHILE,FLAG
*SET,i,i+1  
*SET,MM(i,1),RAND(Rmin,Rmax)
*SET,MM(i,2),RAND(XL+MM(i,1)+d,XR-MM(i,1)-d)
*SET,MM(i,3),RAND(YB+MM(i,1)+d,YT-MM(i,1)-d)
*set,flag1,1
*Do,j,1,i-1
*SET,SHU1,(MM(i,2)-MM(j,2))*(MM(i,2)-MM(j,2))+(MM(i,3)-MM(j,3))*(MM(i,3)-MM(j,3))   
*SET,R2,(MM(i,1)+MM(j,1)+d)*(MM(i,1)+MM(j,1)+d)
*IF,SHU1,LT,R2,THEN
*SET,FLAG1,0
*ENDIF  
*ENDDO  
*IF,FLAG1,GT,0,THEN
CYL4,MM(i,2),MM(i,3),7E-6
*SET,AREA,3.14159265*MM(1,1)*MM(1,1)  
*ELSE   
*SET,i,i-1  
*ENDIF
*DOWHILE,FLAGE
*SET,i,i+1  
*SET,MM(i,1),RAND(Rmin,Rmax)
*SET,MM(i,2),RAND(XL+MM(i,1)+d,XR-MM(i,1)-d)
*SET,MM(i,3),RAND(YB+MM(i,1)+d,YT-MM(i,1)-d)
*set,flagE1,1
*Do,j,1,i-1
*SET,SHU2,(MM(i,2)-MM(j,2))*(MM(i,2)-MM(j,2))+(MM(i,3)-MM(j,3))*(MM(i,3)-MM(j,3))   
*SET,R3,(MM(i,1)+MM(j,1)+d)*(MM(i,1)+MM(j,1)+d)
*IF,SHU2,LT,R3,THEN
*SET,FLAGE1,0
*ENDIF  
*ENDDO  
*IF,FLAGE1,GT,0,THEN
CYL4,MM(i,2),MM(i,3),1E-6
*SET,AREA1,AREA1+3.14159265*MM(i,1)*MM(i,1)
*ELSE   
*SET,i,i-1  
*ENDIF  
*DOWHILE,FLAGEE
*SET,i,i+1  
*SET,MM(i,1),RAND(Rmin,Rmax)
*SET,MM(i,2),RAND(XL+MM(i,1)+d,XR-MM(i,1)-d)
*SET,MM(i,3),RAND(YB+MM(i,1)+d,YT-MM(i,1)-d)
*set,flagEE1,1
*Do,j,1,i-1
*SET,SHU3,(MM(i,2)-MM(j,2))*(MM(i,2)-MM(j,2))+(MM(i,3)-MM(j,3))*(MM(i,3)-MM(j,3))   
*SET,R4,(MM(i,1)+MM(j,1)+d)*(MM(i,1)+MM(j,1)+d)
*IF,SHU3,LT,R4,THEN
*SET,FLAGEE1,0
*ENDIF  
*ENDDO  
*IF,FLAGEE1,GT,0,THEN
CYL4,MM(i,2),MM(i,3),0.5E-6
*SET,AREA2,AREA2+3.14159265*MM(i,1)*MM(i,1)      
*ELSE   
*SET,i,i-1  
*ENDIF
*SET,TT,AREA/TOTALAREA  
*SET,TT1,AREA1/TOTALAREA
*SET,TT2,AREA2/TOTALAREA  
*IF,TT,LE,RATE,THEN
*SET,FLAG,1  
*ELSEIF,TT1,LE,RATE1,AND,TT,GE,RATE,THEN
*SET,FLAGE,1
*ELSEIF,TT2,LE,RATE2,AND,TT1,GE,RATE1,AND,TT,GE,RATE,THEN
*SET,FLAGEE,1
*ELSE
*SET,FLAG,0
*SET,FLAGE,O
*SET,FLAGEE,0
*ENDIF  
*ENDDO

cj2006n

请问ANSYS做微观分析需要看那学东西,我刚接触,希望大家指教,谢谢!

lakerkid

我现在在做堆石模型，在堆石体中除了颗粒就是空隙，因此颗粒占的rate相当大，达到80%，用楼主的方法算到50%就无法计算下去了。故请教有没有更好的方法，使rate较大。现在我很急，望大家帮忙。

chengweifeng

tozhou大侠ansys用到如此地步，小弟只有仰视！！满眼都是羡慕！！

椭圆（卵石），见：
http://www.simwe.com/forum/viewthread.php?tid=788042&highlight=%CB%E6%BB%FA%2B%CD%D6%D4%B2

混凝土细观结构模型，上个世纪八十年代，国外已经开始研究：
F.H. Wittmann,P.E.Roelfstra and H.Sadouki,Simulation and Analysis of Composite Structures,Materials Science and Engineering,68（1984-1985）239-248。

[ 本帖最后由 chengweifeng 于 2008-4-27 13:13 编辑 ]

chengweifeng

多边形：

[ 本帖最后由 chengweifeng 于 2007-7-29 20:49 编辑 ]