金属塑性成型matlab.docVIP

下载本文档

8
0
约3.88千字
约 4页
2018-03-14 发布于河南
举报
版权申诉

金属塑性成型matlab.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

金属塑性成型matlab

%function[J1,J2,J3,r,x1,x2,x3]=ZZ8DPQ(A) %计算主、8面体、等效、偏、球应力，输入：A-应力矩阵 %A=[10 0 3;2 6 1;3 1 5] A=input(A=?) for i=1:3 for j=1:3 str=strcat(请输入A(,num2str(i),,,num2str(j),):); A(i,j)=input(str); end end %A=[10 0 3;2 6 1;3 1 5] %A=[5 0 -5;0 -5 0;-5 0 5] %A=[4 2 3;2 6 1;3 1 5] %A=[50 50 80;50 0 -75;80 -75 -30] %计算第一不变量 J1=A(1,1)+A(2,2)+A(3,3); %计算第二不变量 J2=det(A(1:2,1:2))+det(A(2:3,2:3))+det(A([1,3],[1,3])); %计算第三不变量 J3=det(A); %求主应力赋值r p = [1 -J1 +J2 -J3] r = [roots(p)] S123=sort(r,descend ) %%%%%%% %应力矩阵主对角元素减主应力 B=A; for i=1:3 for j=1:3 B(j,j)=A(j,j)-r(i); %pause end %求主应力方向余弦 [x1 x2 x3]=solve((B(1,1)*x1+B(1,2)*x2+B(1,3)*x3)=0,... (B(2,1)*x1+B(2,2)*x2+B(2,3)*x3)=0,... (x1*x1+x2*x2+x3*x3)=1,... x1,x2, x3); %numeric(x1) lmn=[eval(x1);eval(x2);eval(x3)] end %求主切应力值 t121323=[(abs(r(1)-r(2)))/2 (abs(r(1)-r(3)))/2 (abs(r(3)-r(2)))/2] %求最大切应力值 tmax=max(max(t121323)) %求8面体应力值 FX=[1/sqrt(3);1/sqrt(3);1/sqrt(3)] C=[S123(1) 0 0;0 S123(2) 0;0 0 S123(3)] %计算等倾面上全应力分量 S=C*FX; %求全应力分量的平方赋值为S1 S1=S.^2; %求等倾面全应力赋值为SQ SQ=sqrt(sum(S1)); %求等倾面正应力赋值为SZ8 Sz8=S*FX %求剪应力赋值为SJ Sj8=sqrt(SQ^2-Sz8^2) %求等效应力赋值为SJ Sdx=sqrt(((S123(1)-S123(2))^2+(S123(3)-S123(2))^2+(S123(1)-S123(3))^2)/2) %求应力球张量 Sm=[J1/3 0 0;0 J1/3 0;0 0 J1/3] %求应力偏张量 Sp=A-Sm %画莫尔圆 alpha=0:pi/20:2*pi;%角度[0,2*pi] R=(S123(1)-S123(3))/2;%半径 x=(S123(1)+S123(3))/2+R*cos(alpha); y=R*sin(alpha); plot(x,y,r-) fill(x,y,r); %计算莫尔圆图面范围 if S123(3)=0,XMIN=0,XMAX=S123(1)+1,YMIN=-R-1,YMAX=+R+1; elseif S123(1)=0,XMIN=S123(3)-1,XMAX=0,YMIN=-R-1,YMAX=+R+1; else XMIN=S123(3)-1,XMAX=S123(1)+1,YMIN=-R-1,YMAX=+R+1; end axis equal hold on R=(S123(1)-S123(2))/2;%半径 x=(S123(1)+S123(2))/2+R*cos(alpha); y=R*sin(alpha); plot(x,y,k-) fill(x,y,k); axis equal hold on R=(S123(2)-S123(3))/2;%半径 x=(S123(2)+S123(3))/2+R*cos(alpha); y=R*sin(alpha); plot(x,y,-) fill(x,y,b); %line([0,0],[YMIN,YMAX]) axis equal axis([XMIN XMAX YMIN YMAX])%设定轴取值范围 %quiver(0,YMIN,0,YMAX,r,filled,LineWidth,2); %line([0,0],[YMIN,YMAX],Color,k,LineWid