|
我们来讨论一下bb枪的外弹道需要研究一些什么东西,我总结如下一些物理量:' s: J$ {! f, U5 I) u5 f' U
4 G; N$ N% o* m) J, |1. bb初速度(或者bb初始动能E)---在外弹道中,bb初速v0变为常数,方向当然沿枪管向前& N9 n5 R% D" B5 @1 y
2. hop自转角速度w---这个转速直接影响bb获得的升力(这个升力当然还和初速有关),
7 m$ A/ E" S5 Z: `9 _! Q3 V3. bb速度v随时间t的变化规律,函数v=T(t)( B8 y5 Q- w" @/ @' g
4. bb位移s随v的变化规律,s=S(v),(这个是中间函数,不需要求得也可以。)
/ L1 P ]; U$ y1 L M" {& ~5. bb位移随时间t的变化规律,s=T’(t)5 R& G* O2 m4 ?0 l3 [" @5 p* o$ N; }
6. bb获得的hop升力F的变化规律,这个比较复杂(涉及很多空气动力学的内容),不过可以肯定的是F=F(w,v0,t),即这个升力至少是bb自转转速w,bb初速v0和时间t的函数,并且单调下降。
9 r' ?+ s) _ G( [+ w/ G7. 空气阻力f的变化规律及影响(这里的空气阻力指压差阻力,下文会对空气阻力作详细论述)。这里f=-kv^2 (负号指为阻力,与bb运动方向相反;k为6mm直径bb的空气阻力系数,外弹道学中有更一般的阻力公式f=-0.5pSCMav^2 (p,空气密度;S迎风面积,CMa外弹道阻力常数,Ma,马赫数速度与当地音速的比值,v,bb速度。以后我将用这个公式估算出常数k的大小,以简化我们bb场合的运算。)4 o- V+ h: b: L* k P9 y
& N6 }8 m; p7 r0 T1 D6 V# n: K4 Y) b: i e" m
如果我们求出这些未知的函数,那么我认为bb的外弹道问题就基本解决了。我们先从哪里入手呢?* x, Z* k D& N: w I: p* B
: ]& ^+ M, z) [/ l6 t同时考虑所有这些因素是复杂而且没有办法研究的,我们先从特殊的,省略次要因素的情况入手得出结论后再叠加这些次要因素修正。
; a$ q7 { {4 C' J- y" d6 e$ ~
9 a5 y3 t4 B) _' _
7 ]! o) I8 X$ E1. 水平射击无hop的情况(研究水平方向bb运动变化规律):5 k! N' |+ G( ?$ r7 B9 B; x
5 C- b5 W c& N先说结论,
# t6 C$ l0 | Rbb弹水平方向速度变化方程:
2 E! _. x. W6 R) H2 _1 ^0 m# }位移方程:s=T’(t)= (m/k)[ln(tkv0+m)-ln(m)] (1)
; P2 b- Q, i3 }4 x% Z- z/ ]8 T速度方程:v=T(t)=mv0/(tkv0+m) (2)
; P5 l5 ]% s3 P2 @ k; {9 J$ n5 _5 X常数k=4.37x10^(-6) (3)" k& {" E' t6 n! E2 c+ }
/ N, t: E, U4 u- F& X+ b
推导过程如下:. g9 Z1 B8 G1 H) r7 N; W8 V/ m* _. A8 S
首先我们先忽略hop自转的影响,因为这个hop自转的本质是提供升力减慢bb的下落速度以提高射程,对前进方向的速度影响可以忽略,况且我们感兴趣的是bb前进方向的速度。其次我们先研究水平射击的情况。在这个情况下,bb的重力完全抵消hop升力,在bb前进方向,bb仅受到空气阻力,这样便于研究。8 J: E' b. g! H* {
我们建立方程:) M s4 T6 c. k$ W
由于水平射击,根据牛顿运动定律我们得出:
2 `6 A/ i' V7 Um(dv/dt)=-kv^2
, M+ @: q; R7 A由这个公式直接积分得出v=T(t) 比较困难;不过我们可以由上式得出% Z# k% E: ~9 G, x( Y+ B6 W3 u
dt/dv=m/(-kv^2)0 w! h. n9 n! b ~7 N- `0 N& D
dt= m/(-kv^2)dv
2 r) J; W& a& H( N4 F" [) H我们可以直接积分得出t于t(v)的函数关系(结合积分上下限),* L% _. F- w: Z6 S8 u- m
∫dt=t=∫[m/(-kv^2)]dv
( F9 A* z5 X6 ]) }6 \7 _$ S=-(m/k)∫(1/v^2)dv. G% C+ k4 E- X# f, X
=-(m/k)[-1/v-(-1/v0)]; w( U& L) A) _. j. B
=m/k(1/v-1/v0),% n+ I) m. \0 k8 J6 ^
即t= m/k(1/v-1/v0)。* d: l* |/ ?+ ?7 B, P
1 I* e2 ]% ~ s/ i, Y" f, Q$ Q7 q
变换上式,我们就导出bb速度v与时间t的函数
/ ?2 {0 f! L0 \* {/ ev=T(t)=mv0/(ktv0+m),' D) X( P9 V/ m$ _
. z" R. k* m2 t) z; `$ [1 J9 `' W
为了便于后面积分方便我们再将v变换为
; p+ ]$ Y: H! O: [' s5 l% u5 Fv=T(t)=mv0/(ktv0+m)
. _ x3 u% ^% s) G; f$ u=(m/k){1/[t+m/(kv0)]}
6 U ], `' r" Q0 y; j
* u6 s$ U* d+ o& ~3 h- B" g$ H6 Y3 e9 Z* G/ Q& M; s
将v对dt再积分一次得到位移(结合积分上下限,外弹道开始时t=0)
3 f; q8 p% O" i+ Rs=∫vdt=∫T(t)dt=∫(m/k){1/[t+m/(kv0)]}dt
% w k! f* ^( o8 G( P=(m/k)∫{1/[t+m/(kv0)]}dt' ? r- q0 r! _, W5 E: r, i# r
= (m/k)∫{1/[t+m/(kv0)]}d[t+ m/(kv0)]4 W1 g( X4 f# w. H+ R% u
=(m/k){ln[t+m/(kv0)]-ln[m/(kv0)]}
' O' z4 W- w6 m' w=(m/k)[ln(tkv0+m)-ln(kv0)-ln(m)+ln(kv0)]& V! c" T' r' A: L
=(m/k)[ln(tkv0+m)-ln(m)]
$ C% `1 X. ~' d; e9 _0 Q6 f
- m# q b) @7 a. b: m0 v至此我们导出bb外弹道位移公式和bb外弹道速度公式(水平方向):
; N6 k2 _' B% p8 j: ms=T’(t)= (m/k)[ln(tkv0+m)-ln(m)] (1)* P7 T4 v( j% p- O9 D" M
v=T(t)=mv0/(tkv0+m) (2)
" ]; E) c* @- n6 ~# T
" T3 Z4 Z7 j4 y. T0 a: w由公式可以看出将以上两个公式用于实际计算的最大困难是如何确定这个阻力系数k,我们来估算这个k值。
- h5 C9 z% D9 N* T由f=-0.5pSCMav^2, 我们得到k=0.5pSCMa,我查了弹道表(外弹道学1943年阻力定律),不过弹道表中最小从Ma=0.7开始,即初速238m/s。没有我们想要的数值我们bb一般最多170左右,即Ma=0.5,郁闷!我估算一下取值Ma=0.5时,CMa=0.15# m: r- |6 c3 p9 l
Ma CMa
J! K) h# }; r, w0.7 0.1574 a9 A! c% @; V/ ~0 S+ L
0.8 0.159' O, V9 } n U: {) U3 d
0.9 0.184" S0 O* E5 b+ U7 F
1.0 0.325
. _" d6 l: Y% y% W1.1 0.378
- `$ C: {. {2 ]7 r! @+ X3 M继续网上找资料,找到一个球体空气阻力C的系数取值为0.24。
) D6 E1 o/ ?& E
/ o* R+ c4 {4 O0 t9 e3 i7 J% Q' z考虑到弹道学研究的子弹前部是尖头的锥形的,其形状和我们的bb弹前部半圆形有很大区别。可以想象,圆形的bb在空气中的阻力肯定大于尖头的锥形子弹在空气中的阻力,所以我认为CMa不能取弹道表中的0.15而应取网上的值0.24。
2 a, j- }$ n: O
( |7 k% F) h2 r- N( G取空气密度1.29kg/m^3 S=3.14(0.006)^2/4= 0.00002826
; g3 z+ E, V& W4 xk=0.5x1.29x0.00002826x0.24=0.0000026=4.37x10^(-6)。
- z! l5 L. C' ^( k% c* m
$ v! s4 ~( L( m: A至此我们求的外弹道的两个重要方程及常数k3 B y/ {. b0 C# e' E
s=T’(t)= (m/k)[ln(tkv0+m)-ln(m)] (1)
7 Z$ u- t( \+ a6 O4 vv=T(t)=mv0/(tkv0+m) (2)
' `8 F( m, {: k8 @5 fk=4.37x10^(-6) (3)1 i0 _5 `, p4 B+ {8 _8 \
4 F) W: W. v% ^8 v0 T$ I7 P% k6 |8 o1 e3 z
下面实际应用计算,以我的bar10用0.2gbb大概初速145左右,那么假如灌底的力是120,问在多大射程内还能灌底?
0 r0 O O) ]: S; a% [- h0 Q$ T5 `# K/ \$ v( z3 z
求解:4 z& W% @0 N: a7 ?7 r! X2 {) K$ g
由(2)我们得到 v=T(t)=mv0/(tkv0+m)=0.0002x145/(tx0.00000437x145+0.0002)=120: G, N1 T4 l, {' ~" M' E& X- a
解得t=0.07 代入(1)* O% R' A# {' @+ V4 f- l
得到2 | C6 s% l3 u& l
s=T’(t)= (m/k)[ln(tkv0+m)-ln(m)]
- p, G1 o& Y$ \: t5 R6 Z=(0.0002/0.00000437)[ln(0.07x0.00000437x145+0.0002)-ln(0.0002)]
8 c8 o4 w1 L& O5 b& ]4 j# k=9.16m
; @+ f7 a+ O" D# Z5 ?0 v即在9米左右还能贯底。) {1 R, m! |* Q! C
& v9 h, _0 }, k/ x" b
: k/ Q$ y Q' d( }. j V g2. 水平射击有hop的情况(重点研究竖直方向bb运动变化规律,-马格努斯效应):6 T$ _+ x5 _! R- Y/ V5 V
5 c6 R+ H1 z2 l( z5 B) D
7 s5 w4 g/ i6 b
推导过程如下:- m9 m: p' x( J' n
当水平射击考虑hop影响时,主要研究的是hop升力对bb竖直方向的影响。由hop产生升力的原理我们可以推论出hop自转对bb水平方向的运动影响不大,hop只是增大bb的水平射程。bb水平方向的运动规律还是可以用情况1的方程描述。. g# g3 F4 _0 d0 A: L
- s( y- Y* m, ^3 d/ p( ?: t下面我们来重点研究一下bb竖直方向的运动规律。
5 W7 ], [5 Z; N# g
$ ~/ d9 u- ^7 D( i$ p9 }4 _0 \7 Z
Y. f) @- N# R# f, `+ ]bb在空气中的受力模型(收到三个力,粘滞阻力、压差阻力和上旋升力):, ?" j$ @* G9 a$ q
bb在空中飞行时,bb在向前运动的同时也在旋转。在讨论旋转球体问题时,我们以bb为参照物建立质心参考系。在这个质心参考系中,空气在运动。如果空气是理想流体,那么与bb之间就没有相互作用,bb的运动与它是否旋转无关。但实际空气是粘性流体,它对旋转bb有阻力、同时也受到bb的作用,在bb上下面的空气对bb相对速度不同,产生了压力差,形成了一个“上旋升力”,使bb改变运动方向。而且bb在前进中受到了粘滞阻力和压差阻力。2 }+ d$ c8 W/ s5 P- r0 _
+ E- u) X* h Y& f9 I7 j
a粘滞阻力:6 T$ c" f9 Q+ m, _- Y2 c4 K0 ]) h8 J
由于空气粘滞阻力的存在,使bb旋转角速度减慢。不过因为空气的粘性系数非常小,又因为bb表面很光滑,所以这个粘滞阻力对bb产生的影响非常小,可以忽略。
5 Q! t3 C$ `" T) S# K
' R: l2 G3 Y. f- ^1 T$ e9 ~% vb压差阻力:6 G7 C8 O$ M% ~! J% Z. a
当空气运动遇到bb时,空气会被bb分开,从bb的不同侧面流过。如果空气具有一定的粘性,靠近bb的那部分空气的速度将减慢,在bb的后面一侧形成“真空”地带,离bb较远处的空气将向这个“真空”地带补充,bb后方形成了“尾流”。此时,bb前后两部分空气内单位体积分子数不同,前后侧面受到空气的压力不同,使得bb受到空气的阻力,这种阻力称为压差阻力。在空气动力学中所说的“物体运动时受到空气与速度二次方成正比的阻力”,指的就是空气对物体的压差阻力,压差阻力的大小与物体运动速度的平方成正比,即f∝v^2,可以写为f = kv^2压差阻力的大小与流体的密度、物体的速度有关。这个压差阻力对bb影响巨大,是bb动能衰减的主要原因,情况1已经对这个阻力的影响研究过了,在此就不重复了。
: s7 }: X( E9 s: c9 M' O7 T/ `; H1 d1 L5 ~& Z1 V5 n6 ]8 v
c 附加作用力:
- M5 L: J4 q8 H空气在物体光滑表面的流动为层流,贴近物体的薄层由于空气粘度和bb共同转动,运动规律遵从伯努利方程。伯努利方程是物质能量守恒的体现,说明在流体稳定流动的每一个时刻,空间每一个区域内的能量不变,其数学表达式为:) R8 \4 \2 Q" i
P1+pgh1+pv1^2/2= P2+pgh2+pv2^2/2 3 k) O+ Y3 _9 l# m
式中 P1、P2、h1、h2、v1、v2、分别为空气在1、2两点处的压强、高度和速度,p空气密度;g重力加速度。+ d- \, ^0 S7 C8 W, O
根据伯努利方程,即使1、2两点的高度相同,当它们之间的空气速度有差别时,两点间也产生压强差。空气从bb表面流过。因为从bb上方流过的空气速度v1比从bb下方流过的空气速度v2大,所以在bb上方的气体的压强p1比bb下方气体的压强p2 小。产生了一个向上的力。结果bb一面向前走,一面承受一个把它推向上的力,当这个力和重力相当时,bb就能沿水平方向前运动而不是抛物线。
. X. v' ~. Q+ a, X! D
. f& r1 f- z! R3 Y8 ?% J+ a下面我们来建立方程:( J3 Q. @6 P6 C* V W) N% c
* F ?; ?: p: f+ `3 p( B
显然bb在内弹道中冲开hop胶后在内管中完成加速后出枪口瞬间是“上旋”的。
" D2 o' w, E* C& j/ F取一般情况下9 Z# S! ~( w8 j4 u# x
设bb半径为r,bb的质心以初速度v向前运动,bb自旋转的角速度为w。 5 x, \9 l) W9 |& p! a
选用bb质心参考系,由相对运动可知,此时空气整体上以速度v向后运动;与bb接触的空气在bb表面的作用下绕bb转动。由空气粘度得,在bb表面处的空气与bb面相对静止,以角速度ω绕bb中心旋转。对于质心,bb上侧空气向后方运动,相对bb的速度加大;bb下侧空气向前方运动,相对bb的速度减小。bb上侧空气的速度为v1=v+rw、下侧空气的速度为v2=v-rw。上侧空气速度v1大于下侧空气速度v2/ z$ S5 p2 D4 _) D5 {0 a+ \( _) \
根据伯努利方程可知:
! _0 ~/ D, q1 [% E, X& FP1+pgh1+pv1^2/2= P2+pgh2+pv2^2/2 (5)* x( _, U! e% R# B, y$ S
因为bb上下面高度h1、h2(仅差6mm)基本相同,所以方程可写为:/ \! x1 ^% i6 z/ ]8 ?
P1+pv1^2/2= P2+ pv2^2/22 w0 {- K" G* i. V7 S, v4 A8 G
P2-P1= pv1^2/2- pv2^2/2
6 w# D% ^# M- G+ }& n( n# @=p[(v+rw)^2/2-(v-rw)^2/2] y1 k% ^0 V! M' W* A
=2prwv* ^; G4 F7 p0 b: j( j5 \
又由bb的截面积得到" R- A$ S" y# ~3 d% z! P& T
升力函数) y- }+ J. ^( I
F(t)= (P2-P1)S=2prwv*3.14r*r=6.28pwvr^3
/ o; o: p8 `4 i6 M8 ]4 ?由v=T(t)=mv0/(ktv0+m)$ e" r; a0 o& b2 u5 l
bb竖直方向加速度(向上为正),5 |+ ^) N/ \2 Q, p2 u; `
a=F(t)/m-g
: B5 G* `" a4 I: e* f! A \=6.28pwvr^3/m-g
) l5 O& s$ u0 B% E=6.28pwr^3 [mv0/(ktv0+m)]/m-g
( V% M- D9 Y+ e=6.28pwr^3 [v0/(ktv0+m)]-g
' r! D7 a: L. I# c0 w7 Z2 ~/ N! z2 a4 }8 W# c
说明:由这个加速度a对时间t积分两次即可获得bb竖直偏移的位移S1以时间t的函数,( I( f! A! J- E9 H
又由水平方向bb的位移公式:/ [" t" C# F3 P7 _# r/ x6 [9 [
s=T’(t)= (m/k)[ln(tkv0+m)-ln(m)] 我们可以求出t以水平位移s的函数,将这个t代回到偏移函数S1中,我们就获得了
9 B4 z5 f& i6 M2 s1 F& @bb竖直偏移S1以bb初速v0,自转角速度w和水平位移s的函数。
+ F2 Z K. J q8 Q4 q" T/ Q+ q8 G7 a4 c# B: W
由于没有获得最后的数学表达式,努力补习数学中...........争取早日导出方程......0 T) b. r, U, E+ C) D/ Z+ t
* e. M0 k3 b2 E- a
& d; I: a) q! A6 A; S( Q7 r
" B5 ]0 ^1 v. ^3 m8 o
1 q* F2 d, {6 Z6 B, b! n2 P- Y/ b9 E
3. 枪管有仰角θ(可定义正为向上,负为向下)射击的情况:, I8 |) F. o4 b% `( n
9 X& b% w8 _2 C) k先说结论,; a6 E1 t4 `3 T4 q
枪口向上时,bb速度比水平时慢,hop上旋幅度可能变大也可能变小,在此不能确认,与仰角θ,初速都有关系;枪口向下时,bb速度比水平时快,hop上旋幅度肯定比水平是大。$ c2 |4 ~$ P3 P7 d9 C5 X
" b. T; }1 }8 ]# R' n8 K! [/ I+ e( K- Y/ ~ z+ R
简要推导过程如下(无耐心看的兄弟请略过):
1 i* S7 T3 D; f, z+ X# r/ f! x6 ]当枪管有仰角时,原方程组
3 J% R% {# {& v% x* t- t! v( W' x
水平方向:m(dv/dt)=-kv^2. ?9 e+ B3 P3 e4 X
垂直方向:m(dv1/dt)=mg-F(w,v0,t)2 T' k8 o2 g& x. K' v5 P5 H
5 ^: c) Y( S; w2 g& D
就变为如下的一般形式方程组:: F) {, i+ ]; M7 z
& g4 k$ ?+ a$ G% @% k
bb弹初速方向:m(dv/dt)=-kv^2-mgsinθ
/ i, y2 F! S4 gbb弹初速法线方向:m(dv1/dt)=mgcosθ-F(w,v0,t)0 D; |+ j! }" w! z: P S9 ]
1 i$ H5 B+ W' I4 N
速度,位移的推导过程和上面是一样的,在此就不重复了。只是定性的讲一下,$ L: L q* E. D4 A; Y
枪口向上时,重力的一部分转换成阻力,而抵消hop的重力则变小了;枪口向下时,重力的一部分转换成动力,而抵消hop的重力也变小了。
+ e+ l4 W7 O" t( F" w; L, {4 @( W4 \% ~( D; H, p
1 w6 V1 L1 [6 ]" |
4. 风偏对精度的影响计算" y2 |: Y1 R6 _' G
3 |* a: K; a0 j
; N" R. O& i5 m9 {风力等级表7 @0 P2 f: K$ R6 R
风级和符号 名称 风速(米) 陆地物象 海面波浪 浪高(米)# F: k5 a: [% v
0 无风 0.0-0.2 烟直上 平静 0.06 @* Q( z3 G3 Y8 B2 { `0 S9 N
1 软风 0.3-1.5 烟示风向 微波峰无飞沫 0.1
0 M5 C4 G8 W" j" H4 r 2 轻风 1.6-3.3 感觉有风 小波峰未破碎 0.2
5 C8 p N5 H+ v4 D 3 微风 3.4-5.4 旌旗展开 小波峰顶破裂 0.63 k. ]4 r0 B" h; j
4 和风 5.5-7.9 吹起尘土 小浪白沫波峰 1.06 z' l7 B* u8 N3 |7 F7 [: i0 B
5 劲风 8.0-10.7 小树摇摆 中浪折沫峰群 2.0
/ Q6 X [; K8 ^1 ]" R. L* D0 n( L 6 强风 10.8-13.8 电线有声 大浪到个飞沫 3.0
8 k- o' y; E: E2 `/ v 7 疾风 13.9-17.1 步行困难 破峰白沫成条 4.0
# u0 l& ?4 N" K9 N 8 大风 17.2-20.7 折毁树枝 浪长高有浪花 5.5' O+ e. M( H! b
9 烈风 20.8-24.4 小损房屋 浪峰倒卷 7.0# w7 B+ O9 @9 g+ b1 _$ f5 ~
还是以我的bar10用0.2gbb大概初速145左右,射击20米距离,外弹道时间t=0.173 (求解过程我是用excel编好公式后算的,很方便,而且可以求任意时间、任意位移的值)1 ?- [1 I% r6 v% s( w! A* r9 h
我们假设用微风即,风速6m/s讨论,假设风的影响是垂直弹道的。5 V: _( [1 P0 E
风压的计算问题。0 f0 B( B I0 d3 m; G5 @5 c! r
我们知道,风压就是垂直于气流方向的平面所受到的风的压力。根据伯努利方程得出的风-压关系,风的动压为
, h2 J- _* m* q% C3 P wp = 0.5·ro·v^2 (1)
3 t/ ~+ D# x3 d0 R5 S% h其中wp为风压[kN/m2],ro为空气密度[kg/m3],v为风速[m/s]。8 b# u/ T6 o1 n
由于空气密度(ro)和重度(r)的关系为 r=ro·g, 因此有 ro=r/g。在(1)中使用这一关系,得到
3 j0 y8 {( J) R6 z0 x; z1 N5 U" o6 d. t9 B wp = 0.5·r·v^2/g (2)+ o1 x' o* M; t8 o
此式为标准风压公式。在标准状态下(气压为1013 hPa, 温度为15°C), 空气重度 r=0.01225 [kN/m3]。纬度为45°处的重力加速度g=9.8[m/s2], 我们得到; U9 {1 V) Y9 X1 n+ i5 E1 D
wp = v^2/1600 (3)
* c) l. x" ^& c% A* R得到标准单位的wp=v^2/1.6, bb受风压(4级风)f=6^2/1.6x3.14x(0.006)^2/4=0.00064N; W7 M% K4 m+ _) i( g) ]8 x, g
结合前面得出的t=0.173s3 A4 v; v+ R; k- T1 q
得到20米风偏距离=0.5ft^2/m=0.5x0.00064x0.173^2/0.0002=0.047m
/ k& E* f1 K4 t9 w0 V结论,4级风时bar10用0.2gbb大概初速145左右,射击20米距离会产生最大5cm左右风偏。
" V. U7 y5 |5 L6 ^. c4 J) a4 o4 a7 C. A% a
同样情况下如果换0.32gbb的话,初速114.6, 外弹道时间t=0.2007 + A# E) J; X% f, v* s x7 }# }
得到20米风偏距离=0.5ft^2/m=0.5x0.00064x0.2007^2/0.00032=0.04m. 4cm左右风偏。, ^- q' H8 w7 t/ g4 M. R
% [* c' G3 V( e0 l4 \! M7 x, U[ 本帖最后由 小职员 于 2008-7-15 21:51 编辑 ] |
|