版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
//四参数随机生长法，生成三维多孔介质//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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2026湖北武汉市某国企写字楼会务岗位招聘1人考试备考题库及答案解析
- 2026江苏南京大学图书馆仓储人员招聘考试参考题库及答案解析
- 2026中国矿业大学（北京）招聘9人考试备考题库及答案解析
- 2026中国农业科学院农产品加工研究所食物营养与功能性食品创新团队招聘合同制科研助理考试备考题库及答案解析
- 2026江西九江市彭泽县红光港管理服务中心招聘海关协管员3人考试备考题库及答案解析
- 2026广东深圳九州光电子技术有限公司招聘运维工程师等岗位3人考试参考试题及答案解析
- 2025浙江杭州余杭水务有限公司招聘36人考试参考试题及答案解析
- 2026达州农商银行通川支行寒假实习生招聘考试参考试题及答案解析
- 2026广西桂林市七星区编外聘用人员招聘2人考试参考试题及答案解析
- 2026日照银行见习人员招聘10人考试备考题库及答案解析
- 2026年中国航空传媒有限责任公司市场化人才招聘备考题库有答案详解
- 2026年《全科》住院医师规范化培训结业理论考试题库及答案
- 2026北京大兴初二上学期期末语文试卷和答案
- 重力式挡土墙施工安全措施
- 2025年武汉大学专职管理人员和学生辅导员招聘真题
- 葫芦岛事业单位笔试真题2025年附答案
- 2023届广东省佛山市普通高中高三上学期教学质量检测（一模）物理试题含答案
- GB/T 9163-2001关节轴承向心关节轴承
- GB/T 36342-2018智慧校园总体框架
- GB/T 32264-2015气相色谱单四极质谱仪性能测定方法
- GB/T 25667.1-2010整体硬质合金直柄麻花钻第1部分：直柄麻花钻型式与尺寸
评论
0/150
提交评论