哪位仁兄有一次二阶矩的程序呀？

zhangrh_1217

要是有matlab 和vb的最好都帮传一下啊．

zd83happy

zd83happy

我也很想看看这个程序~~~
:loveliness:

mjhzhjg

楼主可以参考《用Matlab编的计算结构可靠指标的改进一次二阶矩法程序（验算点法）》
http://littlebees.blog.hexun.com/14143949_d.html
或许对你有些帮助，呵呵

civilmanhuijun

要是有源程序就好了，哪位高手能提供呢？谢谢了

pabulum

同求~~~
上面贴的这个matlab的程序，我也搜到过。ms不通用。。。。。

johhan

原帖由 mjhzhjg 于 2008-2-24 18:10 发表

                               
登录/注册后可看大图

楼主可以参考《用Matlab编的计算结构可靠指标的改进一次二阶矩法程序（验算点法）》
http://littlebees.blog.hexun.com/14143949_d.html
或许对你有些帮助，呵呵

mjhzhjg 兄弟给的那个不错,谢谢了!

civilmanhuijun

这个是我自己编制的一个验算点法的算例（赵国藩编写的那本‘工程结构可靠性理论与应用’中的算例），想跟大家多交流交流，是用matlab编写的

format long
syms i f w niter m0 n0 meanf0 meanw0 deltaf deltaw thf thw beta g f1 w1 m1 n1 cosf0 cosw0 cosf1 cosw1 beta0 beta1 betaa;
meanf0=38;
meanw0=54;
deltaf=0.1;
deltaw=0.05;
thf=meanf0*deltaf;
thw=meanw0*deltaw;
g=f*w-1040;
m0=subs(-diff(g,f)*thf,w,meanw0);
n0=subs(-diff(g,w)*thw,f,meanf0);
cosf0=m0/sqrt(m0^2+n0^2);
cosw0=n0/sqrt(m0^2+n0^2);
f0=meanf0+beta0*thf*cosf0;
w0=meanw0+beta0*thw*cosw0;
beta0=subs(solve('(meanf0+beta0*thf*cosf0)*(meanw0+beta0*thw*cosw0)-1140=0',beta0));
beta0=beta0(2);
niter=0;
for i=1:1000000

f1=meanf0+beta0*thf*cosf0;
w1=meanw0+beta0*thw*cosw0;
m1=subs(-diff(g,f)*thf,w,w1);
n1=subs(-diff(g,w)*thw,f,f1);
cosf1=m1/sqrt(m1^2+n1^2);
cosw1=n1/sqrt(m1^2+n1^2);
beta1=subs(solve('(meanf0+betaa*thf*cosf1)*(meanw0+betaa*thw*cosw1)-1140=0',betaa));
beta1=beta1(2);

if abs(beta1-beta0)<10^(-6)
   break
end
niter=niter+1;
beta0=beta1;
end

beta1         %即为可靠度指标
niter           %为迭代次数