二参数的weibull分布随机数序列

lccbaby

如何利用matlab生成一个二参数的weibull分布随机数序列？

jyliao

n2 = wblrnd(1/2,1/2,[1 6])
n2 =
  0.0046  1.7214  2.2108  0.0367  0.0531  0.0917

lccbaby

回复 2# jyliao
哈哈 谢谢。这个我知道怎么弄，我就是不懂别人怎么把这个转换成这样的图的。

xukyle

A=wblrnd(100,2,[1 10000]);
plot(A)

xukyle

A=wblrnd(100,2,[1 10000]);
plot(A)

lccbaby

xukyle 发表于 2011-8-30 22:06
A=wblrnd(100,2,[1 10000]);
plot(A)

真是非常感谢你，哈哈。你画的图的格式和我在别人论文里截图的格式一模一样的。画图的纵坐标是随机数，横坐标是按照间隔为一对应的数。

jyliao

% A=? & Charateristic Stress Range
% B=?  = 1 for EXP  (shape factor)
% freq. (Hz)
% T = total duration
% Ni = Fatigue cycles at stress Xi
dP = 0.01;
Pi = dP:dP:1 ;  %
Xi = weibinv(Pi,A,B);
CumDamage = sum(X./Ni)*0.01*T*freq;

jyliao

A=? ;% Charateristic Stress Range
B=?;%  (shape factor)
freq=?; % (Hz)
T =?; % total duration
Ni =?;%  Fatigue cycles at stress Xi

lccbaby

jyliao 发表于 2011-9-2 21:07
A=? ;% Charateristic Stress Range
B=?;%  (shape factor)
freq=?; % (Hz)

CumDamage = sum(X./Ni)*0.01*T*freq; 这个公式不明白是什么意思。x.是什么？

jyliao

lccbaby 发表于 2011-9-3 12:36
CumDamage = sum(X./Ni)*0.01*T*freq; 这个公式不明白是什么意思。x.是什么？

CumDamage = sum(Xi./Ni)*0.01*T*freq;

何不先说明你的完整问题，例如A、B参数值，作用时间与发生频率等，可提供MATLAB计算程式

lccbaby

我的应力范围服从二参数的weibull分布，A=119.1705，B=1.1266，作用时间为20年，应力范围循环次数为0.6e8。非常感谢你。。。

lccbaby

jyliao 发表于 2011-9-4 07:00
CumDamage = sum(Xi./Ni)*0.01*T*freq;

何不先说明你的完整问题，例如A、B参数值，作用时间与发生频率等 ...

我的应力范围服从二参数的weibull分布，A=119.1705，B=1.1266，作用时间为20年，应力范围循环次数为0.6e8。非常感谢你。。。

lccbaby

jyliao 发表于 2011-9-4 07:00
CumDamage = sum(Xi./Ni)*0.01*T*freq;

何不先说明你的完整问题，例如A、B参数值，作用时间与发生频率等 ...

另一种情况就是交变应力为一零均值的窄带平稳正态随机过程，应力峰值服从瑞利分布。

jyliao

Material S-N Curve?

lccbaby

jyliao 发表于 2011-9-4 19:26
Material S-N Curve?

SN曲线采用国际焊接学会的或者各国船级社的提供的SN曲线，他们给的是应力范围和循环次数之间的关系曲线，具体的参数都有。

lccbaby

jyliao 发表于 2011-9-4 19:26
Material S-N Curve?

材料的sn曲线用国际焊接学会提供的或者各国船级社提供的焊接结构的SN曲线。

jyliao

% A=? & Charateristic Stress Range
% B=?  = 1 for EXP  (shape factor)
% m = fatigue Exponent
% C = Material constant
% Ntotal = total Range cycles
% Ni = Fatigue cycles at stress Xi
A=119.1705;
B=1.1266;
Ntotal = 0.6e8;
dP = 0.01;
Pi = 0.005:dP:1 ;  %
Xi = weibinv(Pi,A,B);
% S-N Curve Ni = C*Xi^m;
C = exp(10.45);
m = -3;
Ni = C*(Xi.^m);
CumDamage = sum(Xi./Ni)*dP*Ntotal % >1 then fatigue

lccbaby

jyliao 发表于 2011-9-5 12:57
% A=? & Charateristic Stress Range
% B=?  = 1 for EXP  (shape factor)
% m = fatigue Exponent

不明白这个Pi在这里代表什么意思，还有后面这个算累积损伤的公式。累积损伤按照Miner线性累积损伤理论公式，你这个公式我理解不了啊。解释一下？

jyliao

Percentile of Stress-range

jyliao

check by:
plot(Xi, Pi)