MATLAB的答案

为什么情节空白?

2视图(30天)
济Zakhem
济Zakhem 2020年12月11日
评论道: 明星黾 2020年12月11日
clc
关闭 所有
清晰的
φ= 0;
f = 0.995;
k = 985;
n = 0.3;
e1B = 0 (21);
e1A = 0 (21);
rhoA = 0.01;
为 i = 21
ea2 (i) = rhoA。* e1A(我); % 2-direction undefected地区的压力
alpha_a = (2 * rhoA + 1) / (2 + rhoA); %在undefected地区应力比
eA (i) = e1A(我)。* sqrt (4/3 * (1 + rhoA + rhoA ^ 2) ^ 2); %在undefected地区有效应变
sigmaA (i) = k * eA(我)^ n; % undefected区域等效应力
sigma1A (i) = sigmaA(我)。/√(alpha_a ^ 2-alpha_a + 1); % 1-direction undefected地区的压力
sigma2A (i) = alpha_a。* sigma1A(我); % 2-direction undefected地区的压力
e2B (i) = ea2 (i);
q_b (i) = e2B (i)。/ e1B(我); %应变率在缺陷区域
alpha_b (i) = (2 * q_b (i) + 1)。/ (2 + q_b(我)); %在缺陷区域应力比
eB (i) = e1B(我)。* sqrt (4/3 * (1 + q_b(我)+ q_b(我)^ 2)^ 2); %有效应变在缺陷区域
sigmaB (i) = k * eB(我)。^ n; %在缺陷区域等效应力
sigma1B (i) = sigmaB(我)。/√(alpha_b(我)。^ 2-alpha_b (i) + 1); %的压力在1-direction缺陷区域
sigma2B (i) = alpha_b(我)。* sigma1B(我); %的压力在2-direction缺陷区域
e1B (i + 1) = e1B(我)+ 0.05;
rhoB (i) = e2B (i)。/ e1B(我);
e3A (i) = 0-e1A (i) ea2(我);
e3B (i) = 0-e1B (i) -e2B(我);
deA (i) = (((f * exp (e3B (i) -e3A(我))* (1 + sigmaB (i)。/ sigma1B(我))* (0.05 + eB (i)) ^ n) / (1 + sigmaA (i)。/ sigma1A(我)))^ (1 / n)) - ea(我);
e1A (i + 1) = e1A(我)+ deA(我);
结束
%的阴谋(e1A (21)、ea2 ' b ');
%抓住
%的阴谋(e1B (21)、e2B ' r ');
%抓住

在回答这个问题。

接受的答案

明星黾
明星黾 2020年12月11日
除了第一个元素,所有的向量 南 。那些不阴谋。(如果你只使用一个标记,你只能看到第一个元素)。
如果你设置这些以外的东西 0 ,如:
e1B = 0 (21) + 1 e-8;
e1A = 0 (21) + 1 e-8;
然后得到一个情节:
图
情节(e1A (21)、ea2“b”);
持有在
情节(e1B (21)、e2B“r”);
持有从
网格
我没有检查你的代码密切确定产生问题( 南 值的结果 0/0 , 正/正 或现有的 南 条目)所以我离开你。在此期间,这些向量的平方根是多少 每股收益 至少产生 的东西 。
2的评论
显示1年长的评论
明星黾
明星黾 2020年12月11日
像往常一样,这是我的荣幸!

登录置评。

更多的答案(0)

类别

找到更多的在应力和应变在帮助中心和文件交换

标签

s manbetx 845

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!