金属塑性成型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 set hx xlim 0 设置 x 轴的界限 if xlim0 0 set hx xlim xlim0 set hx xlim xlim0 xtick0 get h0 xtick 获取原 x 轴刻度 set hx xtick xtick0 设置 x 轴刻度 创建新的 y 轴 hy axes position pos 1 pos 3 2 pos 2 pos 3 1e 5 pos 4 hy axes position pos 1 pos 2 pos 3 1e 5 pos 4 ylim0 get h0 ylim 获取原 y 轴界限 set hy ylim ylim0 设置 y 轴的界限 ytick0 get h0 ytick 获取原 y 轴刻度 set hy ytick ytick0 设置 y 轴刻度 画 x y 轴箭头 如果 x y 轴不需要箭头 下面两句可省略 annotation arrow pos 1 pos 3 2 pos 1 pos 3 pos 2 pos 4 2 pos 2 pos 4 2 annotation arrow pos 1 pos 1 pos 3 pos 2 pos 4 2 pos 2 pos 4 2 annotation arrow pos 1 pos 3 2 pos 1 pos 3 2 pos 2 pos 4 2 pos 2 pos 4 annotation arrow pos 1 pos 1 pos 2 pos 2 pos 4 set h0 Visible off 使原坐标轴不可见 画应力椭球 x y 轴箭头 figure x y z ellipsoid 0 0 0 abs S123 1 abs S123 2 abs S123 3 150 surfl x y z axis equal shading interp camlight lighting gouraud spinmap 10 画平面椭圆 syms x y ezplot x 2 x y y 2 1 1 5 1 5 annotation1 annotation arrow 0 5179 0 5161 0 09524 0 9571 annotation1 annotation arrow 0 08036 0 9536 0 5167 0 519 axis off axis equal 画三维圆柱 f x y z x y 2 z y 2