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

第9章 专题部分

break; }

bid[0]=b[a][0]+(b[a-1][0]-b[a][0])*(alf-jd[a])/(jd[a-1]-jd[a]); bid[1]=b[a][1]+(b[a-1][1]-b[a][1])*(alf-jd[a])/(jd[a-1]-jd[a]);

for(i=0;i<10;i++) {po[i][X]=0; po[i][Y]=0; po[i][Z]=0; }

fz1=0; fz2=0; fm1=0; fm2=0;

fz1=tan((alf+bid[0])*3.14159/180)*pa[X]+tan(bid[1]*3.14159/180)*pb[X]-pa[Y]+pb[Y]; fm1=tan((alf+bid[0])*3.14159/180)+tan(bid[1]*3.14159/180);

fz2=tan(bid[1]*3.14159/180)*pa[Y]+tan((alf+bid[0])*3.14159/180)*pb[Y]-tan((alf+bid[0])*3.14159/180)*tan(bid[1]*3.14159/180)*(pa[X]-pb[X]);

fm2=tan((alf+bid[0])*3.14159/180)+tan(bid[1]*3.14159/180); po[0][X]=fz1/fm1; po[0][Y]=fz2/fm2; po[0][Z]=0;

for(i=0;i<10;i++) r[i]=0;

r[0]=ads_distance(po[0],pa); for(i=1;i<10;i++) { po[i][X]=po[i-1][X]-0.008*H;

po[i][Y]=(po[0][Y]-pd[Y])/(po[0][X]-pd[X])*(po[i][X]-pd[X])+pd[Y]; po[i][Z]=0;

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

acutPrintf(\ %f \ }

for(i=0;i<10;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点坐标 //acedCommand(RTSTR,\); ///////滑动面与各层底交点 for(j=0;j

64

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

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

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

65

第9章 专题部分

{ xtz[j][p][X]=pcj[i][j][X]+p*d[j][0]+d[j][0]/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; }

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[200][200]与tg[200][200] 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; 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;

66

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

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; } xlt[j][p][X]=xtz[j][p][X]-d[j][1]/2;

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][1]/2;

67

第9章 专题部分

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; // }//圆弧端点 for(p=0;p1) { for(p=0;p=(pb[X]+lq)&&xrt[j][p][X]<=(pb[X]+lq+dq)) { zk=zk+zk1+(zk2+zk3)*zk4; mk=mk+(mk1+mk2)*mk3; } else{ zk=zk+zk1+zk3*zk4; mk=mk+mk2*mk3; } acutPrintf(\ } zk11=c[j]*al[j][p]; 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][p][Y])/r[i]*tan(fai[j]*3.14159/180); mk11=q*d[j][1]; mk22=gam[j]*d[j][1]*(tg[j][p]-gh)+d[j][1]*ctz; mk33=(xtz[j][p][X]-po[i][X])/r[i]; if(xlt[j][p][X]>=pb[X]+lq&&xrt[j][p][X]<=pb[X]+lq+dq) { zk=zk+zk11+(zk22+zk33)*zk44; mk=mk+(mk11+mk22)*mk33; } else{ zk=zk+zk11+zk33*zk44;

68