|
- % 输出内点坐标,位移,应力
- fid=fopen('output_interCoorUP.txt','wt');
- fprintf(fid,'输出内点坐标、内点位移、内点应力\n\n');
- fprintf(fid, 'point X Y UX UY sigmaX tauXY sigmaY\n');
- oc(1:NIT,1)=1:NIT;
- oc(1:NIT,2)=CX;% 内点坐标
- oc(1:NIT,3)=CY;
- DSOL_X=zeros(1,NIT);% 内点位移
- DSOL_Y=zeros(1,NIT);
- SSOL_X=zeros(1,NIT);% 内点应力
- SSOL_Y=zeros(1,NIT);
- SSOL_XY=zeros(1,NIT);
- for i=1:NIT
- DSOL_X(i)=DSOL(2*i-1);
- DSOL_Y(i)=DSOL(2*i);
- SSOL_X(i)=SSOL(3*i-2);
- SSOL_XY(i)=SSOL(3*i-1);
- SSOL_Y(i)=SSOL(3*i);
- end
- oc(1:NIT,4)=DSOL_X;% 内点位移
- oc(1:NIT,5)=DSOL_Y;
- oc(1:NIT,6)=SSOL_X;% 内点应力
- oc(1:NIT,7)=SSOL_XY;
- oc(1:NIT,8)=SSOL_Y;
- oc=oc';
- fprintf(fid, '% 1.0f %8.2e %5.2e %18.6e %18.6e %18.6e %18.6e %18.6e\n',oc);
- fclose(fid)
- clear fid oc;
复制代码
我是这样搞的 |
|