糟糕 忘了贴代码 抱歉
n = 0.00067; c = 29.099; S2 = 9.425*10^(-5); h =
8000 T2 = 356; P1 = 100000; V1 = 0.000015; T1 = 280; S1 =
1.70*10^(-6); rou = 1000; h1 = 0.1 g = 10;
NDSolve[{T'[t]*c*n - S2*h*(T2 - T[t]) -
P1*V1*S2/T1*T[t]/(V1 + S2*x[t])*x'[t] == 0,
P1*V1/T1*T[t]/(V1 + S2*x[t]) - P1 +
1/2*rou*(1 - S2^2/S1^2)*x'[t]^2 + rou*(x[t] - h1)*(x''[t] + g) ==
0, x[0] == x'[0] == 0, T[0] == T1}, {T, x}, {t, 0, 3}]
Plot[x, {t, 0, 0.5}]
n = 0.00067; c = 29.099; S2 = 9.425*10^(-5); h =
8000 T2 = 356; P1 = 100000; V1 = 0.000015; T1 = 280; S1 =
1.70*10^(-6); rou = 1000; h1 = 0.1 g = 10;
NDSolve[{T'[t]*c*n - S2*h*(T2 - T[t]) -
P1*V1*S2/T1*T[t]/(V1 + S2*x[t])*x'[t] == 0,
P1*V1/T1*T[t]/(V1 + S2*x[t]) - P1 +
1/2*rou*(1 - S2^2/S1^2)*x'[t]^2 + rou*(x[t] - h1)*(x''[t] + g) ==
0, x[0] == x'[0] == 0, T[0] == T1}, {T, x}, {t, 0, 3}]
Plot[x, {t, 0, 0.5}]