【正文】 Pel=0。 Pek=0。 t=Pj*R/(*Limitstress)。%根據(jù)K2=1求出殼板的厚度%開始進(jìn)入循環(huán)。 while (stressf*Limitstress)||(stress1Limitstress*)||(stress2*Limitstress)||(PekPj)||(Pel*Pj) t=t+。 u=+*E*t^2*(R^2*Pj)。 l=sqrt(R*t)*u/。%根據(jù)殼板的厚度和穩(wěn)定性求出肋骨間距。 %前期系數(shù)r,u1,u2,f1,f2,f4,f5的計(jì)算。 r=^2*R^2*Pj/(10000*t^2)。 u1=u*sqrt(1r)。 u2=u*sqrt(1+r)。 f5=u2*sinh(2*u1)+u1*sin(2*u2)。 f1=sqrt(1r^2)*(cosh(2*u1)cos(2*u2))/f5。 f2=*(u2*sinh(2*u1)u1*sin(2*u2))/(sqrt(3*(1bs^2))*f5)。 f4=*(u1*cosh(u1)*sin(u2)+u2*sinh(u1)*cos(u2))/f5。 %開始帶入肋骨的特征值進(jìn)行循環(huán)計(jì)算。 for n=1:22 A=ribbinga(n)。 I0=ribbingi(n)。 y0=ribbingy(n)。 %定義肋骨對(duì)殼板的影響系數(shù)為B。 B=l*t/A。 k2=1f4/(1+B*f1)。 k1=+f2/(1+B*f1)。 kf=*B*f1/(1+B*f1)。 stress2=k2*Pj*R/t。 stress1=k1*Pj*R/t。 stressf=kf*Pj*R/t。 %定義艙段內(nèi)L的總體穩(wěn)定中的系數(shù) %艙段長(zhǎng)與半徑比為coe1。 coe1=pi*R/L。 %計(jì)及帶板的肋骨慣性矩為I。 I=I0+l*t^3/12+(y0+t/2)^2*l*t*A/(l*t+A)。 %系數(shù)coe2。 coe2=l*t*(R/100)*(R/100)/I。 %開始尋求最小的Pel的值,并將其定為tem。 tem=3*E*I*((10^4*coe1^4*coe2/(2^2+coe1^2)^2+(2^21)^2)/(3*(2^21+*coe1^2)))/(R^3*l)。 for m=3:100 X=(10^4*coe1^4*coe2/(m^2+coe1^2)^2+(m^21)^2)/(3*(m^21+*coe1^2))。 p=3*E*I/(R^3*l)*X。 if p=tem tem=p。 end end %計(jì)算肋骨和殼板的穩(wěn)定性。 stresse=Pj*R/t。 xcs=stresse/Limitstress。 if xcs= cs=1。 elseif xcs cs=+**xcs^2+*xcs^*xcs^4。 if u= cg=1。 elseif ucg=*u+*u^*u^3+*u^4。 end Pel=cg*cs*tem。 Pek=cg*cs*196000*t^2*(R^2*())。 if (stressf=*Limitstress)amp。amp。(stress1=Limitstress*)amp。amp。(stress2=*Limitstress)amp。amp。(Pek=Pj)amp。amp。(Pel=*Pj) t0=t+A/l。 w=2*pi*R*t0。 break else continue end end end disp(t)。 disp(l)。 disp(n)。 disp(Pel)。 disp(Pek)。 disp(stressf)。 disp(stress1)。 disp(stress2)。 disp(w)。 dlmwrite(39。39。,t)。 dlmwrite(39。39。,l)。 dlmwrite(39。39。,n)。 dlmwrite(39。39。,w)。 dlmwrite(39。39。,stress1)。 dlmwrite(39。39。,stress2)。 dlmwrite(39。39。,stressf)。 dlmwrite(39。39。,Pek)。 dlmwrite(39。39。,Pel)。 部分環(huán)肋圓錐殼結(jié)構(gòu)設(shè)計(jì)MATLAB語(yǔ)言程序 %請(qǐng)輸入變量值。 fid1=fopen(39。39。,39。rt39。)。 fid2=fopen(39。39。,39。rt39。)。 fid3=fopen(39。39。,39。rt39。)。 fid4=fopen(39。39。,39。rt39。)。 fid5=fopen(39。39。,39。rt39。)。 fid6=fopen(39。39。,39。rt39。)。 Fid7=fopen(39。39。,39。rt39。) Pj=fscanf(fid1,39。%f39。)。 R=fscanf(fid2,39。%f39。)。 L=fscanf(fid3,39。%f39。)。 Limitstress=fscanf(fid4,39。%f39。)。 E=fscanf(fid5,39。%f39。)。 bs=fscanf(fid6,39。%f39。)。 gama=fscanf(fid7,39。%f39。)。 %首先定義出肋骨的參數(shù)數(shù)組。 ribbinga=[,,,,,,]。 ribbingi=[,,,,1003,1180,1503,1673,2080,2373,3047,3269]。 ribbingy=[,,,,,,]。 %定義循環(huán)變量的初始值。 stress1=1000000。 stressf=1000000。 stress2=1000000。 Pel=0。 Pek=0。 t=Pj*R/(*Limitstress*cos(gama))。%根據(jù)K2=1求出殼板的厚度。 %開始進(jìn)入循環(huán)。 while (stressf*Limitstress)||(stress1Limitstress*)||(stress2*Limitstress)||(PekPj)||(Pel*Pj) t=t+。 u=+*E*t^2*cos(gama).^2*(R^2*Pj)。 l=sqrt(R*t*cos(gama))*u/。%根據(jù)殼板的厚度和穩(wěn)定性求出肋骨間距。 %前期系數(shù)r,u1,u2,f1,f2,f4,f5的計(jì)算。 r=^2*R^2*Pj/(10000*t^2*cos(gama).^2)。 u1=u*sqrt(1r)。 u2=u*sqrt(1+r)。 f5=u2*sinh(2*u1)+u1*sin(2*u2)。 f1=sqrt(1r^2)*(cosh(2*u1)cos(2*u2))/f5。 f2=*(u2*sinh(2*u1)u1*sin(2*u2))/(sqrt(3*(1bs^2))*f5)。 f4=*(u1*cosh(u1)*sin(u2)+u2*sinh(u1)*cos(u2))/f5。 %開始帶入肋骨的特征值進(jìn)行循環(huán)計(jì)算。 for n=1:22 A=ribbinga(n)。 I0=ribbingi(n)。 y0=ribbingy(n)。 %定義肋骨對(duì)殼板的影響系數(shù)為B。 B=l*t/(A*cos(gama))。 k2=1f4/(1+B*f1)。 k1=+f2/(1+B*f1)。 kf=*B*f1/(1+B*f1)。 stress2=k2*Pj*R/(t*cos(gama))。 stress1=k1*Pj*R/(t*cos(gama))。 stressf=kf*Pj*R/(t*cos(gama))。 %定義系數(shù)Coe1。 Coe1=pi*R*cos(gama)/L。 %計(jì)及帶板的肋骨慣性矩為I。 I=I0+l*t^3/(12*cos(gama))+(y0+t/2)^2*l*t*A/(l*t+A*cos(gama))。 %開始尋求最小的Pel的值,并將其定為tem。 tem=(E*t*coe1^4*cos(gama).^3/(R*(2^2+coe1^2)^2)+E*I *(2^21)^2)/(R^3*l)) /(2^2+*coe1^2cos(gama))。for m=3:100 p=(E*t*coe1^4*cos(gama).^3/(R*(m^2+coe1^2)^2)+E*I *(m^21)^2)/(R^3*l)) /(m^2+*coe1^2cos(gama))。 if p=tem tem=p。 end end %計(jì)算肋骨和殼板的穩(wěn)定性。 stresse=Pj*R/(t*cos(gama))。 xcs=stresse/Limitstress。 if xcs= cs=1。 elseif xcs cs=+**xcs^2+*xcs^*xcs^4。 if u= cg=1。 elseif ucg=*u+*u^*u^3+*u^4。 end Pel=cg*cs*tem。 Pek=*cs*196000*t^2**cos(gama).^2/(R^2*())。 if (stressf=*Limitstress)amp。amp。(stress1=Limitstress*)amp。amp。(stress2=*Limitstress)amp。amp。(Pek=Pj)amp。amp。(Pel=*Pj) t0=t+A/l。 w=2*pi*R*t0/cos(gama)。 break else continue end end end disp(t)。 disp(l)。 disp(n)。 disp(Pel)。 disp(Pek)。 disp(stressf)。 disp(stress1)。 disp(stress2)。 disp(w)。 dlmwrite(39。39。,t)。 dlmwrite(39。39。,l)。 dlmwrite(39。39。,n)。 dlmwrite(39。39。,w)。 dlmwrite(39。39。,stress1)。 dlmwrite(39。39。,stress2)。 dlmwrite(39。39。,stressf)。 dlmwrite(39。39。,Pek)。 dlmwrite(39。39。,Pel)。部分球殼結(jié)構(gòu)設(shè)計(jì)MATLAB語(yǔ)言程序 %請(qǐng)輸入變量值。 fid1=fopen(39。39。,39。rt39。)。 fid2=fopen(39。39。,39。rt39。)。 fid3=fopen(39。39。,39。rt39。)。 fid5=fopen(39。39。,39。rt39。)。 Pj=fscanf(fid1,39。%f39。)。 R=fscanf(fid2,39。%f39。)。 Limitstress=fscanf(fid3,39。%f39。)。 E=fscanf(fid4,39。%f39。)。 %定義循環(huán)變量的初始值。 stress1=1000000。Pek=0。 t=Pj*R/(*Limitstress)。%求出殼板的厚度%開始進(jìn)入循環(huán)。 while (stress1*Limitstress)||(PekPj)t=t+。 %計(jì)算殼板的穩(wěn)定性。 Pek=*E*(t/R)^。if (stress1=*Limitstress)amp。amp。(Pek=Pj)w=2*pi*R^2*t。 break else continue end enddisp(t)。 disp(Pek)。 disp(stress1)。 disp(w)。 dlmwrite(39。39。,t)。 dlmwrite(39。39。,w)。 dlmwrite(39。39。,stress1)。 dlmwrite(39。39。,Pek)。