本帖最后由 伟光 于 2016-12-3 22:43 编辑
5 M' U( s9 h5 _! C+ D3 ^& }8 m# Y4 G# j2 k, }; h# l6 q7 W
三代核电压水堆采用非能动安全策略，
0 @4 L) ?, G3 _: W# r3 r6 R. Y采用无轴封的屏蔽式核主泵，& X- _- S+ G M- K) B' h2 A
核主泵为一回路压力边界，' S7 f' P3 M! n5 O% ~! U
推力轴承为水润滑，
$ H* e- v6 j/ ?, S; ~' ~- e8 W惰转飞轮与推力盘复合结构，8 E% S% M# f% h0 l$ Y; P1 u
推力瓦往往为石墨或碳碳复合材料，不耐冲击( D7 Y- O4 n1 d* @; M4 i
在设计阶段如何估算安装时安全落入高度？ K; W5 L8 J5 C) S2 A1 P6 e
+ e, Q6 A+ X& \0 q3 B e% Z' b
上图的物理模型可以进一步简化为下图碰撞数学模型
6 x. K/ I9 i/ @: ?
4 Y7 y# R7 A8 [0 c# {2 [$ T3 l推力瓦的数量为n，转子的质量抽象为一个质点m，转子下降速度为v，
4 L4 {% \4 v+ G5 g1 n# o. o& u+ y推力盘距离推力瓦的高度为h，推力瓦简化一个两点简支梁，
9 F3 z7 I9 L1 C; Z b5 Z# W并且做以下假设：' T( g- f3 ~/ O4 x! U+ ]
a、不考虑推力盘在轴向的微小偏转
- \% `: b6 p# R- G, ?( ~: j6 z5 |b、不考虑推力瓦之间载荷不均匀性
" t- U# c3 Q6 f0 o9 Z; ^4 L$ o5 d/ h9 Ac、不考虑推力瓦具体限位结构带来的周向和径向微小偏转
& @; V* X j7 {9 @d、假设推力轴承为刚性支撑，冲击过程推力盘变形折算至推力瓦处，
$ t; X$ V/ a6 q) fe、假设被冲击的推力瓦重量与转子重量相比很小而且可忽略不计，5 M' `. A* S% Z2 l$ P! |
这个时候转子与推力瓦接触后，可认为附着与推力瓦上而成为一个系统，+ Y- K: J/ p& g3 }& ^) r
在随后的过程中转子沿着冲击方向的运动即可用一个单自由度系统的弹性扰动来模拟：' m- c+ l* T! m1 s& o
+ R9 m U1 X& E' Z) H1 t' @3 I. i. T（dangdangdang，敲下黑板）
# d U# `$ h) A0 s0 _" i2 J; Q e# H2 Q7 Q) U. B$ A
+ q) r7 f$ J3 x力的本质是什么？相互作用，惯性系中的相互作用，
{- |* W! p# |* W# ^4 e5 C在这个碰撞工况下，梁的变形应力F1，就是相互作用，
# f9 l2 }( ~& M3 m碰撞物所在的惯性体系中的等效力F2=Ma., {: k& k8 |7 x; O8 f6 w& q
' C; Z6 Z! f+ C# o: m: Y1 P7 o" q. c' w0 G( I
好啦，冲击力F=F1+F2，8 m9 R" z; M, d( R& e
假设梁的静态变形值为Δ，即位移，' k6 d; c6 B) B& N3 v
根据胡克定律F1=k*Δ，k为梁弹性系数，% B) O, @" C$ V2 [0 c7 q$ C% e
既然是位移，一阶导数为速度，二阶导数为加速度a，
* y& ~- D* n" R所以F2=（Δ的二阶导）*m/n
' e- l, a! L# Q+ |1 d; d
8 u$ s2 x% B I9 `! c
- e" R. P; E, @: Q% q: v7 O7 d整理一下：. p% E# S; ^3 I
/ y+ W" D& B: @8 D4 \5 L1 Y, A' ^+ Z7 [0 L& v: |
哦原来是一个二阶微分方程，对他求解一下：
, K6 j+ L0 i4 G, ]: x
+ Q# O2 J8 }9 w7 {
- M& A5 P: ]/ O9 A8 D根据前面假设，Δ一阶导为V0，Δ在0秒时值为Δ，继续简化为：
- I9 V9 R" P% o U2 c5 L
( h' H$ S G2 x3 K, n- w
' l9 B3 }2 z; S3 l8 N+ L在转子插入屏蔽电机的过程，转子是匀速点动下降，3 I- m# d; r; X+ O% ~5 ~
转子的推力盘在与推力瓦接触前，初速度为0，8 R$ Q1 n2 S( x+ m( J; p, w* o$ ?# ^
在最后一次点动下降之前，推力盘距离推力瓦高度为h；
. O, M9 D6 t0 Y; v在最后一次下降过程中，3 ~+ V) c* o: s4 u9 R
推力盘与推力瓦即将接触的一个微小时间前数值为v0，
: o! S/ x% F1 x, Z; c; D转子的重力势能G转化为转子动能E，根据能量守恒原则，重力势能的变化等于下落动能的增加。
0 H7 L; `3 E% W1 t/ ?7 e
4 ~0 w& r3 P3 D" k
4 L, S x+ b. M! N, n
9 D% U5 F( a8 \- O可见，安装过程中的动载荷冲击是不能简单忽略的，, ~$ X7 }5 G6 d5 _. Q0 f
即便推力盘在与推力瓦之间无限接近处初速度为0的状态下发生碰撞，, A2 Y' K; d* x6 a5 @
碰撞产生的冲击载荷至少为静载荷的两倍。5 t; P* b0 ?& S4 _
: W$ D9 |; G' A0 q& T; r) f根据某型核主泵水润滑推力轴承静态载荷变形计算，累计等效变形为0.07mm，
7 m. b6 i; @4 G! \. F* c9 m+ v5 E! b0 @4 a g9 |
: O& g* A; U- ~: X
某型核主泵推力轴承装配对应的动载荷系数8 `$ l: R! B! V- \& ~: r4 K' E
0 |$ ^' Q3 `% \# I0 y
对于某型核主泵，转子最后一次下落时的高度h不同时，
7 F4 S5 O& e# T# U6 D其动态载荷系数K变动较为明显，如下图所示，- H$ Q5 P$ {6 L8 S3 k
随着高度h从0mm增加至100mm，动态载荷系数从2增加至54.46。5 w; A# R8 s+ R! r0 j% p: U- o1 F
1 S3 Q/ D& T) E+ s
) ^1 f8 P- e7 r+ f% U6 I
" R5 C1 x* P- a0 D6 C核主泵水润滑推力轴承推力瓦材料一般选用石墨或碳碳复合材料，
: c7 @. ~2 D5 K- t2 x6 r* Y三代核电瓦面材料的抗折强度一般应大于90MPa，/ Q' w+ Y7 U: {" S0 u/ L. }) L4 |
可以根据最小抗压强度推算装配过程发生冲击载荷时转子的最大安装落入距离hmax；
) p" U; c# z. ^5 S2 c5 K, {$ O' U瓦面抗折强度为P，最大取值90MPa，此推力瓦面瓦型为阶梯型，最大点静载荷为P0，取值约1.5MPa，
& @$ F, T1 a7 Z6 i8 g) y: p* E& v5 C
7 l) D9 K: Q- j v) v7 T5 a; ?
( T, I8 T) p9 Y/ g2 Y5 w' z% F9 }% R+ r
从上可知，在核主泵推力轴承支撑结构固定时，
/ }8 y7 J T% y/ y+ P瓦面抗折强度增加时，最大安装落入距离增加，
2 v1 ]8 ^$ X% H9 K瓦面最大点静载荷增加时，最大安装落入距离减小，其分布规律如下图所示。% l% U7 D# F4 L- |, H% T
# U! C- \3 o2 l* h' D
6 I7 y" n" X3 S# ^
. h4 W- P, D6 u/ e* p将P=90，P0=1.5代入，可得hmax=121.8mm。
3 l y4 n A# Y6 j2 a) O- N根据以上分析，此核主泵插转子过程中，一旦落入距离大于122mm时，* Q# k4 X7 U& J( n) F
瓦面的局部区域将发生部分浅表层脱落，在安装过程中，推力瓦表面将产生的破碎颗粒。) M- J- B% `; F% F
+ O9 g. o' Y8 p1 [. e" N8 o+ b0 ]! w可见当初速度为0时，冲击载荷主要和两个因素有关，静载荷变形和碰撞距离。
0 f4 w( k) Y8 y0 c. |0 d E4 X
$ m" a8 m' N+ l2 ^3 dPS：江湖传言天下武功唯快不破乎？为何武功高手劈树不砍树梢，请诸君自解。
7 O+ d3 O ` I
. [. k# ?% D" v& N# Z$ p. Q/ v/ o2 h
) T6 @8 P, _4 Z3 r" s$ w9 ^7 O1 d# x3 a3 G
6 Q Y# j" {1 y2 _7 j9 O |
发表于 2016-12-3 22:39:31
