这个问题还不清楚。你能给一个例子,证明你的意思完全和最终的结果是什么?
我怎么能翻译(移动)应力场给定的距离?
2视图(30天)
显示旧的评论
我造型应力场的位错(额外的一半的原子在晶体中)在二维平面上
位错线可以在不同的角度。我可以任何角度旋转错位;但在relaity角由晶体的结构是有限的:θ0°,54.7°或125.2°。如位错必须躺在一个晶面允许这些角度。
我的代码设置材料属性u;νb
他们建立了一个网格。这给x和y的值在二维平面上的应力计算
第二行定义强调Qxx的大小,数量和Qxy x和y的值。
现在如果我想旋转错位,说54.7°;x和y是重新定义为x2, y2强调重新计算QxxPRIME等。
问题的下一步是将错位的距离,年代,以及它坐在飞机。是说,如果θ0°,这将方便与x轴,位错就应该沿着x轴的距离。
如果θ是90°,它沿着y轴的距离
如果θ= 54.7°,它沿着一个平面倾斜于54.7点,距离。
我的问题是正确的位错移动。原点(0,0)必须保持它在哪里。位错是集中在原点。那么x = y = 0时的压力是无限的——这是压力中心的位错是无限的。
所以我猜我不能做的就是添加值x和y,因为x和y不会是零在混乱的中心。
(运动的目的是有几个(最后,成千上万的)混乱沿着1 d飞机在2 d区域并计算整个压力在任何时候在任何时间)。
如果一个代码运行两次,角为54.7 0,一个可以看到旋转的影响,我很高兴。
谢谢你！
u = 80000;%剪切模量,测量在Mpa
ν= 0.3;% Poission材料的比例
b = .00025;%汉堡矢量测量微米
(x, y) = meshgrid (100:1:100 100:1:100);%的区域应力场映射网格
角= 54.7;%的滑移面度
θ=角*π/ 180;%角转换为弧度
Qxx = - * b / (2 (u。*π。* (1-nu))) * (3 * ((y。* x。^ 2 + y ^ 2) / ((x)。^ 2 + y ^ 2) ^ 2));% Qxx压力
Qyy = (u。* b / (2。*π。* (1-nu))) * ((y。* (x。^ 2 - y ^ 2) / ((x)。^ 2 + y ^ 2) ^ 2));% Qyy压力
Qxy = (u。* b / (2。*π。* (1-nu))) * ((x) * (x。^ 2 - y ^ 2) / ((x)。^ 2 + y ^ 2) ^ 2));% Qxy压力
x2 = x。* cos(θ)+ y。* sin(θ);%旋转错位
y2 = y。* cos(θ)- x。* sin(θ);%旋转错位
QxxPRIME = - * b / (2 (u。*π。* (1-nu))) * ((y2。* (3。* x2。^ 2 + y2。^ 2) / ((x2)。^ 2 + y2。^ 2) ^ 2));% Qyy压力'
QyyPRIME = (u。* b / (2。*π。* (1-nu))) * ((y2。* (x2。^ 2 ^ 2)y2)。/ ((x2。^ 2 + y2。^ 2) ^ 2));% Qyy压力'
QxyPRIME = (u。* b / (2。*π。* (1-nu))) * ((x2。* (x2。^ 2 ^ 2)y2)。/ ((x2。^ 2 + y2。^ 2) ^ 2));% Qxy压力'
contourf (x (1:), y (: 1) ', QxxPRIME)
轴平等的
