埋地管道轴向力计算、屈曲判断程序

轴向压力为 P N, 安全系数为 n ,失稳临界力为 W N ,设该管道的圆
2/4
管罐结构设计
截面的惯性矩为 I m^4 ,管道内径为 d mm ,土壤的压缩抗力系数为 K N/m^3*/
{ float P, W, I ,d; float n=0.7, K=0.5E7; d=D-2*b; I=3.14*pow((D*1E-3),4)*(1-pow((d/D),4))/64; {
{
x=p*D/(4*b);
printf("管道自由段轴向应力：x=%0.4f MPa\",x);
}
{
e=(1-2*v)*x/(E*1000)+(a*1E-2)*(T-t);
printf("管道自由段最大轴向应变：e=%f\n",e);
}
/*土壤反力为 Q N ,嵌固段轴向应力为 y MPa ,管道横截面积为 A
油气储运 109 班
#include<stdio.h>
张玉林
学号：0915070133
#include<math.h>
main( )
/*设输油管道直径为 D mm,壁厚为 b mm,管线压力为 p MPa,覆土深度
为 H m，土壤密度为 c kg/m^3，管线安装温度为 t℃,操作温度为 T℃，
土壤与管道的摩擦系数为 u,内摩擦角为 r°,管材弹性模量为 E GPa,线
㎡ ,管道单位长度的摩擦阻力为 f N/m ,管道过渡段轴向应变为 q，
管罐结构设计
过渡段长度为 L m 、热伸长量为 s m*/ { double R; float q; R=tan((45+r/2)*3.14159/180); Q=H*(D*1E-3)*(D*1E-3)*R*R*c*g; printf("土壤反力：Q=%0.4f N\n",Q); y=-E*a*1E-2*(T-t)+v*p*D/(2*b); printf("管道嵌固段轴向应力：y=%0.4f MPa\n",y); A=3.14159*(D*1E-3)*(b*1E-3); f=3.14159*(D*1E-3)*c*g*H*u; L=(((x-y)*1E6)*A-Q)/f; printf("管道过渡段长度：L=%0.4f m\n",L); q=e/2-Q/(2*A*E*1E9); printf("管道过渡段轴向应变：q=%f\n",q); { s=q*L; printf("管道过渡段热伸长量：s=%f m\n",s); } } /*在地下管道的土壤嵌固段，由于温度变化直线管段所受的最大
P=(-v*(p*1E6)*D/(2*b)+(a*1E-2)*(E*1E9)*(T-t))*A; W=2*sqrt(K*(D*1E-3)*(E*1E9)*I); printf("I=%f m^4\n", I); printf("P=%0.4f N\n", P); printf("W=%0.4f N\n", W); } { if(P<=n*W) printf("经校核管道轴向是稳定的。\n"); else printf("管道轴向不稳定\n"); } } }
膨胀系数 a cm/(m·℃),重力加速为 g N/kg ,管道自由段轴向应力为 x
MPa ,管道自由段最大轴向应变为 e */
{
float D=720,b=10,p=3,t=13,T=45,E=210, c=1760,r=30;
float x, e , Q, y, A, f, L, s ;
float H=1.56,u=0.52,a=1.5E-3,v=0.3,g=9.8;
3/4
管罐结构设计
4/4