用MATLAB绘制三向应力圆

1. % 定义主应力值  
2. sigma1 = 100; % MPa  
3. sigma2 = 50;  % MPa  
4. sigma3 = -33;   % MPa  
6. sigma_m1=(sigma1 + sigma3)/2;
7. sigma_m2=(sigma1 + sigma2)/2;
8. sigma_m3=(sigma2 + sigma3)/2;
9. % 计算半径  
10. r1 = (sigma1 - sigma3) / 2;  
11. r2 = (sigma1 - sigma2) / 2;  
12. r3 = (sigma2 - sigma3) / 2;  
13.   
14. % 角度范围（从0到2π）  
15. theta = linspace(0, 2*pi, 1000);  
16.   
17. % 计算圆上的点（以sigma_m为中心）  
18. % 第一个圆  
19. x1 = sigma_m1 + r1 * cos(theta);  
20. y1 = r1 * sin(theta);  
21.   
22. % 第二个圆  
23. x2 = sigma_m2 + r2 * cos(theta);  
24. y2 = r2 * sin(theta);  
25.   
26. % 第三个圆  
27. x3 = sigma_m3 + r3 * cos(theta);  
28. y3 = r3 * sin(theta);  
29.   
30. % 绘制应力圆  
31. figure;  
32. hold on;  
33. plot(x1, y1, 'b', 'LineWidth', 2); % 蓝色表示第一个圆  
34. plot(x2, y2, 'r', 'LineWidth', 2); % 红色表示第二个圆  
35. plot(x3, y3, 'g', 'LineWidth', 2); % 绿色表示第三个圆  
36.   
37. % 添加主应力点   
38. plot(sigma1, 0, 'kx', 'MarkerSize', 10, 'MarkerFaceColor', 'k'); % σ1点  
39. plot(sigma2, 0, 'kx', 'MarkerSize', 10, 'MarkerFaceColor', 'k'); % σ2点  
40. plot(sigma3, 0, 'kx', 'MarkerSize', 10, 'MarkerFaceColor', 'k');

复制代码