clear all;close all
a=0.005;%定义盲孔半径为a=1;
sigama_1=120;
sigama_2=100;%假设实际的残余应力值;
u=0.28;
E=2.1e11;%定义材料为合金钢,给出材料的泊松比及弹性模量;
r=0.010;%选择参考点为半径为r=2圆上的点;
theta=2*pi*rand(1,100);%选择参考点的角度参数;
A=0.5*(1+u)*a^2/r^2;
B=0.5*(3*(1+u)*a^4/r^4-4*a^2/r^2);%推导公式中的A、B的代数式;
H=zeros(100,2);%定义矩阵H;
Y=zeros(100,1);%定义矩阵Y;
S=zeros(100,1);%定义矩阵S;
for k=1:1:100
H(k,:)=1/E*[A+B*cos(2*theta(k)),A-B*cos(2*theta(k))];
S(k)=1/E*[A+B*cos(2*theta(k)),A-B*cos(2*theta(k))]*[sigama_1;sigama_2];%for循环,实际应变推导公式;
end
coeff=0.01;%定义参考系数;
Y=S+coeff*randn(size(S));%测量得出的应变值;
X=inv(H'*H)*H'*Y
a=0.005;%定义盲孔半径为a=1;
sigama_1=120;
sigama_2=100;%假设实际的残余应力值;
u=0.28;
E=2.1e11;%定义材料为合金钢,给出材料的泊松比及弹性模量;
r=0.010;%选择参考点为半径为r=2圆上的点;
theta=2*pi*rand(1,100);%选择参考点的角度参数;
A=0.5*(1+u)*a^2/r^2;
B=0.5*(3*(1+u)*a^4/r^4-4*a^2/r^2);%推导公式中的A、B的代数式;
H=zeros(100,2);%定义矩阵H;
Y=zeros(100,1);%定义矩阵Y;
S=zeros(100,1);%定义矩阵S;
for k=1:1:100
H(k,:)=1/E*[A+B*cos(2*theta(k)),A-B*cos(2*theta(k))];
S(k)=1/E*[A+B*cos(2*theta(k)),A-B*cos(2*theta(k))]*[sigama_1;sigama_2];%for循环,实际应变推导公式;
end
coeff=0.01;%定义参考系数;
Y=S+coeff*randn(size(S));%测量得出的应变值;
X=inv(H'*H)*H'*Y