ansys5.m

利用ansys计算的应力结果

clc;
clear;
[data]=distilldata('PRNSOL1.txt');
data1=data(:,3);
data2=data1';
fid=fopen('sy0.txt','w');
fprintf(fid,'%10.1f\n',data2);
fclose(fid);


%绘制坝基面上垂直正应力与距坝踵距离的关系图 
%坝踵长158.4
%坐标原点在坝踵
load xxx.txt;
load sy0.txt;

x=xxx(:);
y=sy0(:);
%排序
for i=1:35
    if x(i)>x(i+1)
        t=x(i);
        x(i)=x(i+1);
        x(i+1)=t;
        tt=y(i);
        y(i)=y(i+1);
        y(i+1)=tt;
    else continue
    end
end
grid on
plot(x,y)

%合成坝基面上的轴力N,弯矩M
N=0;
M1=0;
M2=0;
M3=0;
b=1;
h=158.4432;
S=b*h;                                %截面面积 
I=b*h^3/12;                           %截面转动惯量 
for j=1:35
    N=N+(y(j)+y(j+1))*(x(j+1)-x(j))/2;
   %M1=M1+(y(j)+y(j+1))*(x(j+1)-x(j))*(x(j+1)+x(j))/2/2;   %以坐标原点为转动中心
   M1=M1+(y(j)*x(j)+y(j+1)*x(j+1))*(x(j+1)-x(j))/2;       %以坐标原点为转动中心
   %M2=M2+(y(j)+y(j+1))*(x(j+1)-x(j))*((x(j+1)+x(j))/2-h/2)/2;   %以截面中心为转动中心
   M2=M2+(y(j)*(x(j)-h/2)+y(j+1)*(x(j+1)-h/2))*(x(j+1)-x(j))/2;       %以截面中心为转动中心
   M3=M3+(y(j)*(x(j)-15)+y(j+1)*(x(j+1)-15))*(x(j+1)-x(j))/2;       %以任一点为转动中心
end

%计算等效应力eq

for n=1:36
eq1(n)=N/S+M1*x(n)/I;
eq2(n)=N/S+M2*(x(n)-h/2)/I; 
eq3(n)=N/S+M2*(x(n)-15)/I;
%eq1(n)=N/S+M1*abs(x(n)-h/2)/I;
%eq2(n)=N/S+M2*abs(x(n)-h/2)/I;
end
hold on
grid on
plot(x,eq1,'r')

hold on
grid on
plot(x,eq2,':')

hold on
grid on
plot(x,eq3,'k')