版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
//四参数随机生长法，生成三维多孔介质//voidQuartetStructureGenerationSet(){ doubled1_6=0.02; doubled7_18=d1_6/2.; doubled19_26=d1_6/8.; doublen=0.8; doublecdd=0.001; doublenumtotal_need=(1-n)*NX*NY*NZ; intnumsoild=0; intTnumsoild; intsoild[M+1][3]; srand((unsigned)time(NULL)); //第一步，遍历所有网格，生成固相内核 for(i=0;i<=NX;i++) { for(j=0;j<=NY;j++) { for(k=0;k<=NZ;k++) { arrgrid[i][j][k]=0; soild[numsoild][0]=0; soild[numsoild][1]=0; soild[numsoild][2]=0; if(((rand()%1000)/1000.0)<cdd) { arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; numsoild=numsoild+1; } } } } Tnumsoild=numsoild; //第二步，从固相核向周围26个方向生长 while(Tnumsoild<numtotal_need) { for(intindex_soild=0;index_soild<Tnumsoild;index_soild++) { intindex_i=soild[index_soild][0]; intindex_j=soild[index_soild][1]; intindex_k=soild[index_soild][2]; //1向右方向生长 if(index_j<NY-1) { i=index_i; j=index_j+1; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d1_6) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //2向后方向生长 if(index_i>0) { i=index_i-1; j=index_j; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d1_6) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //3向左方向生长 if(index_j>0) { i=index_i; j=index_j-1; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d1_6) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //4向前方向生长 if(index_i<NX-1) { i=index_i+1; j=index_j; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d1_6) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //5向上方向生长 if(index_k<NZ-1) { i=index_i; j=index_j; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d1_6) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //6向下方向生长 if(index_k>0) { i=index_i; j=index_j; k=index_k-1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d1_6) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //7向水平方向右前生长 if(index_i<NX-1&&index_j<NY-1) { i=index_i+1; j=index_j+1; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //8向水平方向左前生长 if(index_i<NX-1&&index_j>0) { i=index_i+1; j=index_j-1; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //9向水平方向右后生长 if(index_i>0&&index_j<NY-1) { i=index_i-1; j=index_j+1; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //10向水平方向左后生长 if(index_i>0&&index_j>0) { i=index_i-1; j=index_j-1; k=index_k; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //11向右上方向生长 if(index_j<NY-1&&index_k<NZ-1) { i=index_i; j=index_j+1; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //12向右下方向生长 if(index_j<NY-1&&index_k>0) { i=index_i; j=index_j+1; k=index_k-1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //13向左上方向生长 if(index_j>0&&index_k<NZ-1) { i=index_i; j=index_j-1; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //14向左下方向生长 if(index_j>0&&index_k>0) { i=index_i; j=index_j-1; k=index_k-1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //15向前上方向生长 if(index_i<NX-1&&index_k<NZ-1) { i=index_i+1; j=index_j; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //16向前下方向生长 if(index_i<NX-1&&index_k>0) { i=index_i+1; j=index_j; k=index_k-1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //17向后上方向生长 if(index_i>0&&index_k<NZ-1) { i=index_i-1; j=index_j; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //18向后下方向生长 if(index_i>0&&index_k>0) { i=index_i-1; j=index_j; k=index_k-1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d7_18) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //19向右前上对角线方向生长 if(index_i<NX-1&&index_j<NY-1&&index_k<NZ-1) { i=index_i+1; j=index_j+1; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d19_26) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //20向右后上对角线方向生长 if(index_i>0&&index_j<NY-1&&index_k<NZ-1) { i=index_i-1; j=index_j+1; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d19_26) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //21向左后上对角线方向生长 if(index_i>0&&index_j>0&&index_k<NZ-1) { i=index_i-1; j=index_j-1; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d19_26) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //22向左前上对角线方向生长 if(index_i<NX-1&&index_j>0&&index_k<NZ-1) { i=index_i+1; j=index_j-1; k=index_k+1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d19_26) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //23向右前下对角线方向生长 if(index_i<NX-1&&index_j<NY-1&&index_k>0) { i=index_i+1; j=index_j+1; k=index_k-1; if(arrgrid[i][j][k]==0&&((rand()%1000)/1000.0)<d19_26) { numsoild=numsoild+1; arrgrid[i][j][k]=1; soild[numsoild][0]=i; soild[numsoild][1]=j; soild[numsoild][2]=k; } } //24向右后下对角线方向生长 if(index_i>0&&index_j<NY-1&&index_k>0) { i=index_i-1; j=index_j+1; k
温馨提示
- 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年度年福建省高校教师资格证之高等教育心理学模拟考试试卷B卷含答案
- 2024年度山西省高校教师资格证之高等教育法规考前练习题及答案
- 历史教师培训心得体会
- 2024年度茶叶批发销售协议范本
- 2024年私人贷款协议样式
- 房产买卖居间服务协议2024全攻略
- 2024年家庭装修协议
- 2024游乐场设施租赁协议模板
- 2024年居间合作项目协议精简
- 2024年跨境资本贷款协议示例
- 《光伏发电工程安全预评价规程》（NBT 32039-2017）
- 是否混声四部合唱谱
- 从局部到整体：5G系统观-完整版
- 《白酒基础知识培训》课件
- 第4课西汉与东汉统一多民族封建国家的巩固教学设计-高中历史必修中外历史纲要上册
- 团队建设创造和谐的人际关系与文化
- 新时代中小学思政课一体化建设探究
- 陪诊项目商业计划书
- 第3.2课《简单相信傻傻坚持》（课件）-【中职专用】高二语文同步课件（高教版2023·职业模块）
- 小学科学项目式学习教学设计与实践研究
- 家纺面料未来发展趋势报告
评论
0/150
提交评论