边坡分析，ARX, CAD二次开发，程序源代码(3)

辽宁工程技术大学毕业设计（论文）

mk=mk+mk22*mk33; } } else{ zk11=c[j]*al[j][0]; zk22=q*d[j][1]; zk33=gam[j]*d[j][1]*(tg[j][p]-gh)+d[j][1]*ctz; zk44=(po[i][Y]-xtz[j][0][Y])/r[i]*tan(fai[j]*3.14159/180); mk11=zk22; mk22=zk33; mk33=(xtz[j][0][X]-po[i][X])/r[i]; if(xlt[j][0][X]>=pb[X]+lq&&xrt[j][0][X]<=pb[X]+lq+dq) { zk=zk+zk11+(zk22+zk33)*zk44; mk=mk+(mk11+mk22)*mk33; } else{ zk=zk+zk11+zk33*zk44; mk=mk+mk22*mk33; } } gh=gh+ch[j]; ctz=ctz+gam[j]*ch[j]; }

k[i]=zk/mk; }

/////////保存OD线上的计算结果/////////////////////////////////////////// K=0; KR=0; ke[X]=0; ke[Y]=0; ke[Z]=0;

for(i=0;i<10;i++) { K=k[0]; KR=r[0]; ko[X]=po[0][X]; ko[Y]=po[0][Y]; ko[Z]=0; ke[X]=pe[0][X]; ke[Y]=pe[0][Y];

69

第9章 专题部分

ke[Z]=0; if(k[i]

//////////////////////Dmo垂线/////////////////////////////// ads_point pt1,pt2,pt3,pt4; pt1[X]=ko[X]+1; pt1[Y]=-(ko[X]-pd[X])/(ko[Y]-pd[Y])*(pt1[X]-ko[X])+ko[Y];//od上的点 pt1[Z]=0; pt2[X]=(pa[X]+pb[X])/2; pt2[Y]=(pa[Y]+pb[Y])/2; pt2[Z]=0;//ab垂线上的点 pt3[X]=pt2[X]+1; pt3[Y]=-(pa[X]-pb[X])/(pa[Y]-pb[Y])*(pt3[X]-pt2[X])+pt2[Y]; pt3[Z]=0;//ab垂线上的点 ads_inters(pt1,ko,pt2,pt3,0,pt4); //求OD线上的圆心边界点 po[0][X]=ko[X]; po[0][Y]=ko[Y]; po[0][Z]=0; for(i=1;i<6;i++) { po[i][X]=ko[X]+0.008*i*H; po[i][Y]=-(ko[X]-pd[X])/(ko[Y]-pd[Y])*(po[i][X]-ko[X])+ko[Y]; po[i][Z]=0; } for(i=6;i<10;i++) { po[i][X]=ko[X]-0.008*i*H; if(po[i][X]

//////////////////////////////od垂线上的试算//////////////////////

70

辽宁工程技术大学毕业设计（论文）

for(i=0;i

r[i]=ads_distance(pa,po[i]); for(i=0;i

{ pe[i][Y]=pb[Y];

pe[i][X]=po[i][X]+sqrt(r[i]*r[i]-(pb[Y]-po[i][Y])*(pb[Y]-po[i][Y])); pe[i][Z]=0;//求E点坐标 for(j=0;j

pcj[i][j][X]=po[i][X]+sqrt(r[i]*r[i]-(pcj[i][j][Y]-po[i][Y])*(pcj[i][j][Y]-po[i][Y])); pcj[i][j][Z]=0; } ///初始化th[],n[],d[][] for(j=0;j<100;j++) { th[i][j]=0; n[j]=0; d[j][0]=0; d[j][1]=0; } //求个层滑动面宽 th[i][0]=pe[i][X]-pcj[i][0][X]; if(cs>1) for(j=1;j

////////////////////为各层分条，以0.01R为一条/////////////// for(j=0;j

acutPrintf(\ n[j]= %d \ } //初使化xtz[][] for(j=0;j

71

第9章 专题部分

/////各土条底端中点与滑动面的交点 for(j=0;j

xtz[j][0][Y]=po[i][Y]-sqrt(r[i]*r[i]-(xtz[j][0][X]-po[i][X])*(xtz[j][0][X]-po[i][X])); xtz[j][0][Z]=0;} else{

for(p=0;p

xtz[j][p][Y]=po[i][Y]-sqrt(r[i]*r[i]-(xtz[j][p][X]-po[i][X])*(xtz[j][p][X]-po[i][X])); xtz[j][p][Z]=0;

xtz[j][p][X]=pcj[i][j][X]+p*d[j][0]+d[j][1]/2; xtz[j][p][Y]=po[i][Y]-sqrt(r[i]*r[i]-(xtz[j][p][X]-po[i][X])*(xtz[j][p][X]-po[i][X])); xtz[j][p][Z]=0; } }

//初始化stz[100][100]与tg[100][100] for(j=0;j<100;j++) for(p=0;p<100;p++) { stz[j][p][X]=0; stz[j][p][Y]=0; stz[j][p][Z]=0; tg[j][p]=0; }

//////////求土条高度 for(j=0;j

///////////初始化zk,mk,ctz,gh,xlt[][],xrt[][],al[][] for(j=0;j<100;j++) for(p=0;p<100;p++) { xlt[j][p][X]=0;

72

辽宁工程技术大学毕业设计（论文）

xlt[j][p][Y]=0; xlt[j][p][Z]=0; xrt[j][p][X]=0; xrt[j][p][Y]=0; xrt[j][p][Z]=0; al[j][p]=0; } zk=0; mk=0; gh=0; ctz=0; zk1=0; zk2=0; zk3=0; zk4=0; mk1=0; mk2=0; mk3=0; zk11=0; zk22=0; zk33=0; zk44=0; mk11=0; mk22=0; mk33=0; //////11111///#####////////求安全系数///////////#######///////////////////////////////// for(j=0;j

xlt[j][p][Y]=po[i][Y]-sqrt(r[i]*r[i]-(xlt[j][p][X]-po[i][X])*(xlt[j][p][X]-po[i][X])); xlt[j][p][Z]=0; xrt[j][p][X]=xtz[j][p][X]+d[j][0]/2;

xrt[j][p][Y]=po[i][Y]-sqrt(r[i]*r[i]-(xrt[j][p][X]-po[i][X])*(xrt[j][p][X]-po[i][X])); xrt[j][p][Z]=0;

73