版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
//四参数随机生长法，生成三维多孔介质//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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 新疆阿克苏市阿瓦提县第四中学2025届物理高二第二学期期末学业质量监测试题含解析
- 二零二五年度建筑材料采购合同
- 二零二五年度个人车辆抵押贷款简易合同
- 2025年度车辆个人抵押权登记及管理合同
- 二零二五年度安置点配套设施建设及施工一体化合同
- 2025版生物科技孵化中心入驻及研发合作合同
- 2025版厂房物业安全防范系统维护合同
- 二零二五版矿山尾矿处理居间服务合同
- 2025年度餐饮外卖配送服务合同范本
- 2025版车库购置及车位分配合同示范文本
- 胸痛单元建设课件介绍
- 超市消防安全管理制度制度
- 2025年广西中考语文试题卷（含答案）
- 2025年南京市中考数学真题试卷
- 2025年呼伦贝尔农垦集团有限公司招聘笔试参考题库含答案解析
- 2025年重庆市中考数学试卷真题（含标准答案）
- 2024广西专业技术人员继续教育公需科目参考答案（97分）
- 湘少版小学全部英语单词
- 电子秤cal标定方法
- 艺术品及品牌项目合作协议书
- 管道工字钢支护施工方案
评论
0/150
提交评论