MATLAB的答案

我似乎不能找出如何阴谋最终答案theta-x轴和一个,B、C在y方向上。3个单独的情节。

12个视图(30天)
费利佩•阿尔瓦拉多
费利佩•阿尔瓦拉多 2022年9月30日
编辑: dpb 2022年9月30日
我尝试着不同的方法如何阴谋θx轴的值从3 B和C不同的情节。甚至我已经麻烦弄清楚如何让不同θ,B和C,这是唯一的方法我可以使它工作。但是现在我不能似乎情节。我的方式我可以画出我所有的不同的θ所有不同的值,B和C,但情节不同图表而不是绘制在图1。如果任何人都可以帮忙。这是复合类的代码,我们使用常量职位减少合规和刚度矩阵,然后把他们从地方到全球和策划一些数据。
函数(年代,Q) = IsotropicComplianceStiffness (E1, E2、v12 G12)
%的输入:E(扬斯模量、GPa), v(泊松比),G剪切模
E1 = 40 * 10 ^ 3;
E2 = 8 * 10 ^ 3;
v12 = 0.25;
G12 = 4 * 10 ^ 3;
sigxx = 5;% MPa
sigyy = 4;% MPa
txy = 2;% MPa
%输出:C-stiffness矩阵,S-compliance矩阵
%初始化C和S矩阵,通过设置所有矩阵元素为零
S = 0 (3,3);
Q = 0 (3、3);
%计算方面的合规矩阵。
%记住:合规矩阵是对称的。
(1,1)= 1 / E1;
S (2, 2) = 1 / E2;
S = v12 / E1 (2, 1);
(3,3)= 1 / G12;
%的对称
(1、2)=年代(2,1);
%计算刚度矩阵。
%记住:刚度矩阵是对称的,逆的C。
Q =发票(年代);
为i = 0:15:90%做θ(15 0 - 90的增量)
θ=我
m = cos(θ);
n =罪(θ);
%的变换矩阵
T = [m ^ 2 n ^ 2 * m * n;n ^ 2 m ^ 2 2 * m * n;- m * n m * n ((m ^ 2) - (n ^ 2)));
%的倒数变换矩阵
Tinv =发票(T);
R = [1 0 0;0 1 0;0 0 2];
Rinv =发票(R);
条形= R S * T * Tinv * Rinv *;%将减少合规矩阵
一个Q = Tinv * * R * T * Rinv;%转换降低刚度矩阵
%全球应力应变计算
Stress_G = [sigxx; sigyy txy];%给压力在全球
Strain_G = S * Stress_G;%在全球找到应变(合规矩阵* stressglobal)
%局部应力和应变calcualtion(使用矩阵变换)
StressL = T * Stress_G% stressL = stresslocal (stress1 stress2,剪切应力(6))
StrainL = T * Strain_G;% strainL = strainlocal (strain1 strain2,σ(剪切应变(6))
%的比例和角度问题2.3合理——我们使用
一个= StressL (1,1) / sigxx% Stresslocal(1方向)/压力全球(1方向)
%我们使用stressL(1, 1)因为这是我们stressL在1方向)
B = StressL (2, 1) / sigyy% Stresslocal(2方向)/压力全球(两个方向)
C = StressL (3,1) / txy%剪应力本地/剪切应力全球
图
plot1 =情节(θ)
1评论
dpb
dpb 2022年9月30日
你的代码不是完整的…循环结束,但是
…
图
plot1 =情节(θ)
将创建一个新的人物和情节一行。没有其他可以做的只有θ,一个论点。
情节上的所有三个相同的情节,接电话 图 的循环,然后使用 抓住 添加额外的数据相同的情节。
…
hL =情节(θ,[A B C]);%我用hL作为通用的处理
注:然而,上面只会把一个点集的情节在每个θ值,它不会在这两点之间画一条线。你上面的构造,你覆盖 A, B, C 每个通过循环;您需要为每一个为数组分配空间和存储与索引数组每个通过循环来拯救他们。
也做,那么移动 情节 循环之前调用的一次性完成和情节。

登录置评。

在回答这个问题。

接受的答案

dpb
dpb 2022年9月30日
编辑:dpb 2022年9月30日
跑:
E1 = 40 * 10 ^ 3;
E2 = 8 * 10 ^ 3;
v12 = 0.25;
G12 = 4 * 10 ^ 3;
IsotropicComplianceStiffness (E1, E2、v12 G12);
函数(年代,Q) = IsotropicComplianceStiffness (E1, E2、v12 G12)
%的输入:E(扬斯模量、GPa), v(泊松比),G剪切模
sigxx = 5;% MPa
sigyy = 4;% MPa
txy = 2;% MPa
%输出:C-stiffness矩阵,S-compliance矩阵
%初始化C和S矩阵,通过设置所有矩阵元素为零
S = 0 (3,3);
Q = 0 (3、3);
%计算方面的合规矩阵。
%记住:合规矩阵是对称的。
(1,1)= 1 / E1;
S (2, 2) = 1 / E2;
S = v12 / E1 (2, 1);
(3,3)= 1 / G12;
%的对称
(1、2)=年代(2,1);
%计算刚度矩阵。
%记住:刚度矩阵是对称的,逆的C。
Q =发票(年代);
θ= [0:15:90]。”;%定义所需的角度
= 0(大小(θ));B =一个;C =;% preallocate输出
为i = 1:元素个数(θ)%迭代的角度
θ=θ(我);%设置迭代θ值
m = cos(θ);
n =罪(θ);
%的变换矩阵
T = [m ^ 2 n ^ 2 * m * n;n ^ 2 m ^ 2 2 * m * n;- m * n m * n ((m ^ 2) - (n ^ 2)));
%的倒数变换矩阵
Tinv =发票(T);
R = [1 0 0;0 1 0;0 0 2];
Rinv =发票(R);
条形= R S * T * Tinv * Rinv *;%将减少合规矩阵
一个Q = Tinv * * R * T * Rinv;%转换降低刚度矩阵
%全球应力应变计算
Stress_G = [sigxx; sigyy txy];%给压力在全球
Strain_G = S * Stress_G;%在全球找到应变(合规矩阵* stressglobal)
%局部应力和应变calcualtion(使用矩阵变换)
StressL = T * Stress_G;% stressL = stresslocal (stress1 stress2,剪切应力(6))
StrainL = T * Strain_G;% strainL = strainlocal (strain1 strain2,σ(剪切应变(6))
%的比例和角度问题2.3合理——我们使用
(我)= StressL (1,1) / sigxx;% Stresslocal(1方向)/压力全球(1方向)
%我们使用stressL(1, 1)因为这是我们stressL在1方向)
B (i) = StressL (2, 1) / sigyy;% Stresslocal(2方向)/压力全球(两个方向)
C (i) = StressL (3,1) / txy;%剪应力本地/剪切应力全球
结束
hL =情节(θ,(A B C),“x -”);
传奇(“一个”,“B”,“C”,“位置”,“eastoutside”)
结束
是一种解决——
注意:在原始函数的参数变量被设置为固定值在函数;没必要让他们在这种情况下参数。我最小的修复纠正,然后叫这些值的函数。
函数会更有用的一般如果它还在θ向量可以通过改变不改变代码本身。
返回变量的年代,问不像数组,保存;我只固定那些需要创建请求的块;相同的修正需要做这些。
如果返回其他压力,那么该函数不需要做任何绘图(如果没有更好的分解)但会通过它的计算,让调用者决定什么以及如何阴谋。
这些清理和修复左为“锻炼学生”

更多的答案(0)

类别

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

标签

s manbetx 845


释放

R2020b

社区寻宝

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

开始狩猎!