节点应力时间图是锯齿状的

chenyujie711

做了个简单的拉伸试验，输出节点的应力时间图，有时候会是锯齿状的。 这是什么原因？
会不会是因为力外推算法的原因？ 能不能只将最大值连成一平滑曲线。

附图见zip。

huaijuliu

如果单纯取最大值连成曲线 可以采用MATLAB
% Find upper and lower envelopes of a given signal
% The idea is from Envelope1.1 by Lei Wang, but here it works well when the signal contains
% successive equal samples and also includes first and last samples of the signal in the envelopes.
% inputs:
%   sig: vector of input signal
%   method: method of interpolation (defined as in interp1)
% outputs:
%   upperenv: upper envelope of the input signal
%   lowerenv: lower envelope of the input signal
function [upperenv lowerenv] = envelope(sig, method)
if nargin == 1
    method = 'linear';
end
upperind = find(diff(sign(diff(sig))) < 0) + 1;
lowerind = find(diff(sign(diff(sig))) > 0) + 1;
f = 1;
l = length(sig);
try
    upperind = [f upperind l];
    lowerind = [f lowerind l];
catch
    upperind = [f; upperind; l];
    lowerind = [f; lowerind; l];
end
xi = f : l;
upperenv = interp1(upperind, sig(upperind), xi, method, 'extrap');
lowerenv = interp1(lowerind, sig(lowerind), xi, method, 'extrap');
-----
用法如
sig=[1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11];
envelope(sig,'linear')

ans =

  Columns 1 through 12

    1.0000    2.0000    2.5000    3.0000    3.5000    4.0000    4.5000    5.0000    5.5000    6.0000    6.5000    7.0000

  Columns 13 through 20

    7.5000    8.0000    8.5000    9.0000    9.5000   10.0000   10.5000   11.0000

huaijuliu

[upperenv lowerenv] = envelope(sig, 'linear'); plot(sig); hold on; plot(upperenv);

chenyujie711

谢谢， 感激涕零呀~
能不能给解释一下，对于这种简单的金属板拉伸实验，节点的应力怎么回事锯齿状？

3# huaijuliu

huaijuliu

我自己算了一哈 不会啊

chenyujie711

附上cae file ， 最大应力点node3，输出mises-time的图的确是锯齿状的，s11和max principle是曲线。

btw，判断金属试件是否失效，是不是该用mises的值？
5# huaijuliu

huaijuliu

你如果是试件拉伸 判断准则应该是第一强度理论最大拉应力准则吧
单向拉伸

chenyujie711

恩，我找书看看哈，von mises 屈服准则和强度理论都啥时候用哈。俺初学者，谢谢ls大力支持哈。 7# huaijuliu

iambadman

金属件的拉伸压缩用VON.MISES的比较多。

一般塑料件塑性不好，用第一主应力和主应变的比较多



8# chenyujie711

chenyujie711

恩，但是mises-time图形是锯齿状是怎么回事，不知道能不能解释一下？
我用最大值的连线表示mises可以么？
9# iambadman

iambadman

没看到CAE文件啊？我也不知道你具体模型什么样子

Robert_Su

金属件的拉伸压缩用VON.MISES的比较多。

一般塑料件塑性不好，用第一主应力和主应变的比较多



8# chenyujie711
iambadman 发表于 2010-10-3 09:31

我的經驗中,

塑料件如PP,POM等靜態拉伸曲線的塑性區很大,

和iambadman的經驗不同.

請參考.

chenyujie711

附上cae file ~ 11# iambadman

chenyujie711

13# chenyujie711

chenyujie711

不好意思，六楼的时候忘了附cae了，13 14楼附上~
Node 3 的stress是锯齿状，难不成是step 里的frequency 选的不对？
11# iambadman

chenyujie711

我用的是钢，但node3 的 mises为什么会是锯齿状呢 12# Robert_Su

chenyujie711

这种拉伸试验，可以说当node3超过yield point的时候，试件就失效么？ 11# iambadman

chenyujie711

最后决定用node3 附近的点了~  唉唉  谁要是知道该怎么画平滑的最大应力点的mises 麻烦通知下，以后再改了~

iambadman

你的最大值点严重奇异，你选择的节点是有9个MISES应力数值的。如果你非要看一个单数曲线，把averaging threshod 修改为100%就可以了。

顺便给你说一下就是你这个网格根本不适合做你需要的结论。非要给别人看的话把网格隐藏掉吧。

iambadman

同意，POM等材料的塑性是非常好的。


但是即使是POM，我看到的资料用MISES的也非常少。

第二：POM材料虽然有20%，甚至30%的静应变极限，但是实际在工程中，很少这样用，一般都是5%左右，或者8%，9%的都有，特别是动载荷中，极端的设计甚至是用3%，不支持用20%这样的参数，这个问题我倒是正好曾经专门和宝理的工程师以及他们的总经理什么的聊过。


另外就是我也说了一般来说～～塑料是这样的，可能不严密吧，毕竟工程塑料种类还是挺多的，差别也不小。


12# Robert_Su