check by:
plot(Xi, Pi)
理解起来还是有困难。Xi = weibinv(Pi,A,B)，这个函数是什么意思？是代表应力范围服从参数为AB的weibull分布当所占百分比为pi是对应的应力范围值吗？ Read the New one
% S-N curve....N(x) = C/(x^m)
m = 3;
C =10^10.45;
% Weibull Stress-Range distribution
A = 119;
B = 1.226;
%B = 1 for a Exp. distribution
%B = 2 for a Rayleigh distribution
Ntot = 0.6E8;% Total stress-cycles
Damage = @(x) Ntot/C*(x.^m).*wblpdf(x, A, B);
% Define integration upper bound
xup = wblinv(0.9999999, A, B);
Q = quadl(Damage, 0, xup);
% Exact solution, ref: Lalanne, 2Ed, V4, p.244, Eq.[**0].
Qex = Ntot/C*(A^m)*gamma(1+m/B);
jyliao 发表于 2011-9-7 18:50 static/image/common/back.gif
Read the New one
% S-N curve....N(x) = C/(x^m)
非常感谢。以后以问题再向学习。。。
页:
1
[2]