找回密码
 注册会员

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 321|回复: 0

计算机视觉二维刀具预调仪的研究

[复制链接]
发表于 2010-9-12 10:59:10 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转磨削论坛

您需要 登录 才可以下载或查看,没有账号?注册会员

×

1 概述

5 F D7 K( S5 u5 l
对刀仪是数控机床及加工中心的必要设备,它的作用是迅速而准确地确定数控机床及加工中心刀库中每把刀具的几何参量,以便在加工过程中用程序控制刀具的进给量。现有的对刀仪产品大部分是机械光学式的,依靠人眼通过影屏对刀具进行瞄准测量,从而带来主观误差,使对刀精度受到影响。目前国产的一些对刀仪只能测出刀形两个坐标方向上轮廓顶点的坐标值,仅能用于切削圆柱面和端面的程序编制。在某些曲面的加工中,刀具与工件的接触点并非上述轮廓顶点(如锥面加工),且接触点可能随时变化(如球面加工)。刀尖轮廓具有一定的尺寸,只有已知刀形的中心坐标及随切点角度变化的刀具圆弧半径,才能由程序对刀具半径进行补偿,从而给出理想的切点运动轨迹,这样对刀具刀尖部分轮廓的测量就是十分必要的。 # Q* y, o9 v7 I
国外的某些对刀仪产品已经能够进行准确的刀形测量。根据对刀具半径的补偿,可以精确加工出任意形状的回转面。它采用计算机视觉检测系统,以面阵CCD摄像机摄取刀形,通过通用图像板将刀形图像转换成每个像元有256个灰度级的灰度矩阵送入计算机,利用软件进行刀形轮廓提取。这种方法运算量大,成本相应较高,测量速度相对较慢。 % Y6 ^, R% y j; F. Q+ m
本文所研究的计算机视觉二维刀形检测系统与一般的计算机视觉检测系统不同。它没有采用通用的图像板,而是针对刀形轮廓成像时具有很强的黑白对比特性,由自行研制的硬件解调板采用硬件方法直接提取刀形轮廓。该系统同时实现了刀形的自动测量及刀具预调工作。在刀形测量数据处理方面,本系统采用分段圆弧拟合的方法给出刀形曲线的圆心及分段圆弧半径,便于程控加工中调用。本系统具有测量速度快,精度高,重量轻,结构简单的特点,刀形测试及刀具调整的精度为3µm(3s)。利用该方法设计的刀形检测系统已成功用于天津机床光学仪器厂的对刀仪产品。
% z9 P+ C2 K, E# [! Z & S' m* k! q& |. a9 `5 ?' O: @* e) T h$ {6 w+ Q7 h" \; n) [7 r- c

图1 视频信号处理电路框图
; Y, j6 J0 A2 M/ ^. N( Q2 u" b

2 边缘提取硬件电路

0 O0 @) B; E" j% R& J7 p
为了减小像敏元间距给测量带来的影响,提高系统分辨率,本系统根据Purll D.C.提出的解调测量思想,以硬件的方法从视频信号中解调出图像轮廓边缘特征点,从而将空间尺寸量转换成时间方波。利用插入高频稳定的晶振时钟脉冲,方便地实现了尺寸细分。硬件视频信号处理电路框图如图1所示。 $ O! u! f- D! N; G, C2 c6 R8 f, j! k+ v
电路应该解决三个主要问题:(1)将CCD摄像机输出的视频采样信号恢复成连续的视频信号:(2)在连续的视频信号中确定拐点的位置作为边缘特征点:(3)与PC机的时序相匹配。
/ K$ ^! J8 q1 D2 _1 R- i2 {* \% _2 y6 {' w- N& Q& Q" Y$ U O8 \5 p6 ^) K: Z8 q+ {4 S$ y4 F! j& n

图2 处理电路各部的输出波形
$ _; P/ p: U3 Z1 m7 w8 o
根据刀具边缘成像具有很强的黑白对比度的特性,我们选用具有线性相移的低通滤波器作为对视频信号进行解调的网络。其输出波形是已恢复的连续信号,如图2a所示。为确定该连续信号的拐点(该点对应被测件的轮廓点),对连续信号进行二次微分。一次微分电路、绝对值电路、二次微分电路的输出波形如图2b、c、d所示。二次微分电路输出波形中前4个波形分别对应视频信号的行消隐前沿、行同步前沿、行同步后沿和行消隐后沿。后面波形的过零点对应图2a中的拐点。通过过零比较器和脉冲选择逻辑电路,得到与轮廓点对应的脉冲,如图2e所示。选择行消隐信号的后沿作为每一扫描行的起始点,如图2f所示。起始点信号与轮廓点信号通过方波形成电路得到系列方波,如图2g所示。在方波中插入高频稳定的脉冲并由计数器记录该脉冲数,则每一行的边缘轮廓点的X坐标值得以确定。Y 坐标值由扫描行序数确定。
. v0 g) O9 `: g7 L3 b4 {4 {" ~) q& N7 X" J6 T; i) d, c( L( q& T$ k# T% `& ?; h6 M6 L0 s; {; f

图3 与PC机的时序匹配
/ h$ U. U+ f; n8 _ \" C7 T
图3表明了扫描行时序与PC机的匹配情况。计数器的清零应在测量起始点之前,因此选择行同步信号后沿作为计数器清零信号。行消隐信号的前沿作为向PC机申请中断的信号。在中断服务程序中完成数据采集工作。在中断申请信号与计数器清零信号之间要有足够的时间,以保证中断服务程序执行过程中数据采集可靠。场信号作为向PC机申请中断的另一信号,用以区分奇偶场并将奇偶两场合并为一幅完整的图像。
8 {9 s4 p) i2 m2 \7 r; }0 c6 A ; I' g9 m% R/ G4 o9 m. D1 t2 V# I' A0 J; [" v, F' R8 j* ]: k1 G5 A ^. N% g6 v. }% |2 c9 k$ K/ W+ @% [3 M

图4 对刀仪零点棒

图5 各个坐标系之间的关系
$ i5 I( F3 P" I1 g! a2 ?' D

3 对刀仪定标及坐标系统变换

1 W3 L% A6 c Y: [
在测量正式进行之前,需要对测量系统定标,即确定在X方向上每个脉冲数所代表的量值及在Y 方向上相邻两行扫描线之间的距离。 4 ^7 N. c2 H+ f+ [
本测量系统借用对刀仪中的零点棒作为定标中的标准件。与普通对刀仪不同之处在于采用端点中心处有标准圆球的阶梯轴作为零点棒,如图4所示。将视觉系统的行扫描线水平安装。通过视觉系统读出阶梯轴上阶梯所对应的脉冲数,用阶梯的标准值除以该脉冲数,得到脉冲当量即X标值。将视觉系统精确旋转45°(此时为视觉系统的工作位置),对标准棒母线进行测量,然后对测点进行直线拟合,选相距较远的两拟合点,其坐标值为(X1,Y1)、(X2,Y2)。由:(Y2-Y1)/(X2-X1)=∆Y/∆X=cos45°=1,Y2-Y1=n·∆Y,得出Y方向的标值,即行距∆Y。其中:n为二采样点的行间隔数,∆Y为行间距。 2 y- [1 ^* ~7 c ]+ t! ?* q
在刀形检测系统中有4个坐标系:机床坐标系(X,Y),由安装在机床X,Y方向导轨上的光栅尺输出读数:视觉检测系统中与(X,Y)平行的坐标系(x,y):视觉检测系统坐标系(I,J),其横轴与CCD摄像机的行扫描方向平行:被测刀具的极坐标系(R,q),极点选择在刀尖圆弧的圆心处,q=0°时极径的方向与视觉检测系统的轴重合。为使测量能够正确进行,需将各个坐标系统一。视觉检测系统的测量值最终要统一到机床坐标系中。图5给出各个坐标系之间的关系。 7 A4 x/ H D% r; [1 s
面阵CCD相对机床坐标系呈45°安装,以避免刀形受检部位与CCD扫描线平行而无清晰的交点,因此视觉坐标系首先要旋转45°,使之与机床坐标系平行。x=I·cosq-J·sinq:y=I·sinq-J·cosq。旋转后视觉系统的坐标值转换到参考坐标系中,若(x,y)坐标系的坐标原点在(X,Y)坐标系中的坐标为(x0,y0),则(x,y)在(X,Y)中的坐标为:X=x+x0:Y=y+y0。视觉系统坐标值通过参考坐标系转换到机床坐标系中可用下式表示: / a: H7 z. k' D: E $ o* U0 V: s2 n( g7 O) P% }' P2 g N! J9 d7 G! {5 @0 Q6 ~) @- t! A6 B" p! V% ~& o0 \; {0 ]3 d/ a8 _9 `. Y! Y j1 n" P' O+ |3 h' Z& o3 S4 n9 z1 X' B% J, Q% |7 s3 k! r/ B: `% q! A6 h0 R" p0 |/ ~" o/ f: ]1 i6 @ U( K, p4 `, P. x0 e/ U; \! M) A1 {. K+ k2 O+ i. Y- S: ^0 K7 s4 b0 G' l; ^! ]. b" F- x7 h) w0 e. |- T; }8 e" [3 ?5 s4 s/ C7 G; k; }& s7 X: F* u) e- v7 w5 g$ } Y' O0 |, }* H/ A/ k5 K( q* o- D- P* h- ?" ~( Y( D2 R3 ^& n/ D! f2 N! l3 l3 U* U/ e1 Z! f8 M) {& o4 Q' R! E4 f% m1 W' d R6 O( z- C- B* H* ]0 \ H
[ X ] = [ cosq -sinq ][ I ] + [ x0 ]
Y sinq cosq J y0
式中,q=45°。这里并没有考虑参考坐标系与机床坐标系的小夹角a,在测量进行之前应预先测量出夹角a的大小,而上式中的q应为45°-a
' K& |- ]; ?/ `" `% ?& m 6 g# @" N" B6 L7 @& A: t/ j# E/ W6 @. B1 M- e5 Z9 C6 u$ i9 k6 r1 g' a, q; Y" O9 m5 a! i; }/ H. D7 R e. w3 V- `! d) Q) y2 w# a' _* {0 I( C8 R! Q/ J0 a. S, _+ T, ]5 Q: ^# W' C: ]8 k! T3 `1 ?1 V2 ?: E0 |0 o4 q6 V+ e. p- O' i+ d3 p$ F, }* @. h+ v/ w( f8 p* \; p. X1 V8 D9 W+ C! L* c4 Y& [# s- `" i8 |: T% x' m A# g& k9 g' j6 U7 m2 i8 J- L9 f; ^" Y6 b6 c$ \. `3 p$ J. u( Q2 W6 y) O' m5 G% L: a" ^5 j8 h9 s1 O! ?! q2 D8 m/ N1 R# q" {9 c* E3 T( q- d5 g; x" q( Z7 V% w: Y1 B9 L- Y J: u+ @7 a& @0 Z3 O; A7 l! j1 j/ @/ F( t6 K! [- x3 b5 _. a( g) y" }9 c1 g: x2 `) X2 ?+ e" d. V" ?3 f: j! D/ S1 ]* i9 `* }( i7 `$ r8 t: R* Q+ u2 R: u4 u2 K# _) A3 a& c& m1 q1 [2 @1 I( x8 T# ^& y; e5 E- x# T$ A! A; Z% d8 B/ i' o1 M1 |0 L# E& ^$ l0 Y4 u' W) T" @, o" j0 L. Q/ c! S: r- T% |) m# ?' C9 N n+ ^& S* [, M5 V1 l: i) P, n* |& N( j) C$ Q; A, E) n7 M2 E. X4 ^7 O( a' {0 y8 h4 X8 j6 z( l$ W, H) r' O* Z4 s# o; f- Z( p" ^6 w7 Z" z9 |1 P$ z. J3 E3 }- M" O6 v* c$ o9 z+ ~4 G7 \ a# F( Z9 m" r, M1 H9 K3 I% w- T( e2 ^# A& p3 O, f+ W5 F' {4 P6 B, u( P; ^6 U5 s+ Z+ v; G# a5 H4 z! m W% w7 r: O% s; V4 a: ?& z4 _/ J7 y2 B) `7 r2 G, v0 Z( S4 I0 z. p3 e0 l. {. r' W* X: e* c+ Z) |* h9 G; q+ U0 y6 D$ D$ T# i" v- g' p8 j S; o) C" V: W/ i5 v" p! a7 k P, O- ]$ j3 U( d) L8 f! a2 I/ l( ]5 y3 z- z! L- M8 V; s2 ]$ F8 U' ]6 K& @6 X' r, t4 p- M* J+ G/ r7 R+ P z& r& g. r; }; y' q( g5 v7 K1 e! ~) Y, A+ @5 j/ N$ A9 ^; j4 r' C I7 [* Y9 G k# f8 H" }1 f' K& N2 k: Y& }7 W3 u; |: S6 e# H! K( l# z: w" N2 @* o$ x' p4 g- p0 t3 ~8 J; D" Q: d( [0 @8 U1 w- k! m! b6 H0 P% c1 J1 R5 k3 F3 \% ^' T' [! y: \6 f0 }, ~. u9 M3 k% `* t! \9 b% R( p1 ?0 K6 f, b4 I) m/ {' x9 n) D0 t6 P1 a- |% |5 Z, c) Q# g+ c4 n8 |2 H! r* _) A% ]5 m/ ^' V4 G G# F4 u: L1 P7 t: _9 @- u1 T$ }6 W
表1 (单位:µm)
测序 圆心坐标(X,Y) 圆弧半径 Xtop Ytop
1 1110.1,770.2 502.3 249 860
2 1109.2,768.6 503.4 251 859
3 1110.3,769.9 502.8 248 859
4 1110.1,769.7 503.0 251 859
5 1109.1,769.7 502.5 251 860
6 1110.6,770.5 503.7 248 859
7 1109.7,768.4 502.4 250 861
8 1110.0,770.3 502.4 249 860
9 1110.5,769.4 501.8 248 859
10 1110.0,770.4 502.3 251 860
3s X:1.5 Y:2.2 1.7 3 2
/ l4 }# c- _$ U

4 测量软件

% q- y8 P) W+ ?) n- n% g& L
本系统涉及的测量软件包括数据采集程序与数据处理程序两部分。数据采集程序完成将被测刀具的边缘轮廓点坐标采入计算机的任务。数据处理程序完成被测刀具的各分段圆弧半径的计算及X,Y方向的顶点和加工初始点的计算,依据采入计算机的刀具边缘轮廓点坐标,用最小二乘法进行分段圆弧拟合及计算X,Y方向各轮廓点坐标的极值,具体算法从略。 ; N$ w6 S; k* Q0 R- ^

5 实验结果

; C) x4 |# |* v
采用公称直径为1mm的标准钢球作为测量对象,来模拟对被测刀具的刀尖圆弧测量。取中心角为120°圆弧上的轮廓点坐标进行最小二乘拟合计算,每隔一段时间测量一次。表1中第2、3列为该圆弧圆心坐标和半径的10次测量数据,第4、5列为被测球在X方向与Y方向上两个顶点坐标的测量数据,其精度可代表对刀仪的对刀精度。
$ M9 R' O8 q3 h3 l, H : ?* G7 ?* o5 M% T8 L8 L( A- K- b: r$ l5 [! S: a/ l+ Q0 F _# g7 H Y+ g5 n* k- K) @7 m' @2 |/ D0 u0 ]% A/ @8 ~2 G5 C3 E! p6 h3 a* I4 I( x4 f" L' G- r# r' x; O9 D( e' x0 A, {9 k: G3 u4 f+ ?& a' \% P* R. e+ R n5 S* ]4 C4 f) z8 N# d! ~1 q% }, b/ d+ }. W* S$ z) T; N6 k& h; v) ?/ \: z/ o, v& r0 c* p7 o" r2 ^; a+ f8 {3 U; x i( z7 ^0 E! }0 u, `0 E. _ U' B- K# j* f! ~$ ^, e) V% N0 o3 L1 H& h% w$ y6 |& i3 a w+ U& T8 I3 v- p, ^/ d6 q$ \( `' @- d+ E1 b" Y+ e# _: v5 S$ E& @* e; f9 S" R4 v2 G$ l1 K8 x! q8 b+ O+ i+ f% ]( \6 v$ J4 N2 k3 X% }( O: @9 G7 s5 j) B& ~/ d6 F; o1 E8 v3 w# P: [9 U, V, k8 _' |: M0 q4 E) j( F1 h/ G% t* r% ] S6 _4 H7 I# J7 j* i x, z' B+ I9 x2 R7 m6 l" w5 Y& D |, ?0 g9 c; J) H @0 F4 E5 E) `- T9 j& e* u7 L6 M5 C4 r' v$ C8 j {" K% r: s! l/ _ K* s' I; ]. }" |# O+ ^8 |% I* E" a* h0 a- d- t+ P+ ~5 E2 F3 b) J R! [/ }/ e" ^: C, y" o& Z/ t. g0 v; g9 J1 l. D! t1 ^3 v( M# w4 i. ?5 z# N- z6 x0 d- W, d, Q# [; z! G+ C- W1 ?4 B) ?; I5 v( s0 U# w1 ]1 Y+ ]) M( V9 K! E* G$ z$ L3 D2 p8 h; g+ W/ w( X4 Q# `" O+ |2 [! Y* b- O; _4 \! ?" B) g6 w1 u% ~5 e6 }. z6 H3 Y3 |# p9 q6 O, J! g0 t0 N1 Q7 V, O1 f+ T3 I: C( R" S! g/ @' ? K3 D' [1 X- R' B1 L6 J+ K( K" C8 M9 j: }/ i$ p8 F' S1 x' x# Y0 O+ H, j& y* z- }* v! Q# q+ f9 s3 o$ Z3 ]) Y7 W% D' R7 r, a7 ?6 l. m, z; G! I6 m, _0 O# A& r+ I5 c. ^0 O5 u. @/ c8 k. E" h" j- r, r A r0 n% W/ l: e' q1 i# G1 s% M# ~; v2 h& o3 V6 A; Y* U! w1 ?: d3 X, y" @& ^: b0 F! S7 _5 v/ R' C' x. {( {' S: [9 d, p: f0 o4 K ` q% d; \9 w# n* H$ H% ~% G: p5 f, s8 W. B4 B8 f3 d, W0 J3 j) E: ~& ]8 c1 c" h2 a6 }% x5 q; U: i" @1 S# }/ |3 g y: k( H/ l3 A* W M' d* E( E5 k8 C0 N: g& x! H i1 F/ M: M8 U w8 H+ u) h$ z, o5 x: }6 i- o) q5 D3 m5 Q9 q9 [: [8 U5 u3 L: @( w4 N1 N9 Q$ @5 g4 j+ Q g1 P. X: w; l. F; Q2 I& y/ t% ~6 p, P6 S5 l% V2 D$ K0 u; z9 @7 ^* |8 G$ |/ Q7 g- u- ]. [6 m2 _! I/ [9 P) z& |' O# J5 W- I# @$ f- a; g8 n/ ?2 ~& y+ B' f$ v* N7 h: t, E" ?2 F+ l( [) F1 D3 J) }2 M9 q6 ^/ K7 n8 ]; d8 _8 t1 ]; R& T' e! Y1 X) E+ p% X; O0 r# ]) K+ m& H6 b1 A- r/ h6 R7 H" x4 x5 Z4 u1 W6 J5 M$ A4 K& G/ M3 t- k+ r/ B/ m1 n/ Y1 K& B7 ~# i) B# C) l0 A- a8 l4 \& i" c* _+ d. J5 b: J: _0 c" @* ]) g' r: N$ C0 {; Y* @# V5 N" W+ O t- `* z& K0 U3 `1 p* ^4 N5 `6 e3 r) J X) ]# I, F" W5 D) [! [$ r. M. }( O! V% ~& C9 g/ x _0 ~% s9 [9 ~* }. H7 [- U% a# d3 y+ `# @6 \, G0 i$ ^6 {7 T9 R9 W5 O2 t. A$ f9 o# |# t5 B% r4 ?5 V3 [' X9 \' S* A+ ~3 E) u7 y, v7 j9 `0 R6 h* U0 q' W; w1 z) G" ?" w2 r$ J) d! k1 l' }) } _+ ]6 T* L8 b& {& |1 v5 }( C6 ~! b6 K4 B9 k" }! a( s8 T3 i0 G, e# W$ Y' A/ O! m* j- }- e# z, |1 N; U/ E6 ]$ ^3 N' b; t s5 u' h9 m0 Q: V$ B7 E! N- x3 {! Z& Q/ l; w/ N9 g9 k+ z2 C+ ~! j e( B+ {6 }3 Q& K/ U0 ?5 j5 q5 \8 I( m- c& n/ y& m* T6 ]8 w% x; d c3 R, g! n1 m. W! q) v" Z& @; ~1 n% e2 `9 {6 y+ }/ m9 ^7 z5 f U5 |. w+ O9 Q7 @7 ~6 i. j( ~9 N* D7 R- r! |8 g5 u9 J' S; ?, K9 }2 s" M# {4 `, ?0 j" s5 c4 g- w6 ~3 @' P2 J' y( Z& l1 T- }+ q* L2 P; p. d+ ^9 E7 d3 u* k, N0 |2 n) F* X! Y$ j; K/ V- O y3 z* {8 o/ d6 W2 d; b( K, e/ i! O G- D/ t. O- O' t- X. u0 k* p R- ^7 d# C* u1 S( w# t x& \4 [$ l. L/ m7 Z6 U% p6 G0 a$ s$ \9 U# F2 c* m1 i* m9 N+ b1 W/ v2 ]* p2 u/ {" v: |0 ~' X, y6 g5 o1 U8 L, N/ S0 z2 o; m# d, f$ d" T, T' t, L0 ~* Z( H, Y) c( Q R) D" N. E, [! a6 [7 G2 m& d4 k% ~& K
表2 (单位:µm)
测序 10° 20° 30° 40° 50° 60° 70° 80° 90° 100° 110° 120°
1 504.3 504.3 503.4 503.3 503.4 503.7 504.0 504.2 504.6 504.5 504.0 503.9
2 504.7 504.6 503.8 503.8 504.0 504.2 504.5 504.8 505.1 504.9 504.2 504.1
3 505.1 505.2 504.3 504.4 504.5 504.7 505.0 505.3 505.6 505.4 504.7 504.5
4 505.0 505.1 504.4 504.4 504.5 504.7 505.0 505.3 505.6 505.4 504.5 504.3
5 504.7 504.6 503.9 503.9 504.0 504.1 504.5 504.8 505.1 504.9 504.3 504.1
6 504.1 504.0 503.1 503.1 503.2 503.5 503.8 504.0 504.4 504.3 503.7 503.6
7 504.8 504.9 504.0 504.0 504.1 504.3 504.6 504.9 505.2 505.0 504.3 504.1
8 505.0 505.0 504.2 504.1 504.3 504.5 504.8 505.1 505.4 505.2 504.4 504.2
3s 1.1 1.2 1.3 1.3 1.2 1.2 1.2 1.3 1.2 1.1 1.0 0.9
( R) m2 X' X/ @. v- f
对分段圆弧半径的测量,将120°中心角的范围平均分成12段,每段取8次测量结果做统计分析。计算仍采用最小二乘法拟合,但计算时应加上圆心固定的约束条件。表2给出其测量数据,取各段中最大误差作为刀形测量误差,经计算±3smax=2.6µm。 0 w$ B+ m) @( Z& ?1 P* M$ q

5 结论

$ l8 ~' e2 X7 [- A- t* e
该视觉检测系统应用解调测量方法对视频信号进行处理,用硬件方法直接提取轮廓特征点,测量精度高,速度快。应用在国产对刀仪上可方便地进行刀具预调工作和对刀具的刀形参数进行高精度测量。通过对被测刀具刀尖半径的实时补偿,能够提高数控加工的精度。该系统结构简单,重量轻,可方便地进行机上对刀。
& b* F0 m% B, U/ c8 v+ X. H
您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

中国磨床技术论坛
论 坛 声 明 郑重声明:本论坛属技术交流,非盈利性论坛。本论坛言论纯属发表者个人意见,与“中国磨削技术论坛”立场无关。 涉及政治言论一律删除,请所有会员注意.论坛资源由会员从网上收集整理所得,版权属于原作者. 论坛所有资源是进行学习和科研测试之用,请在下载后24小时删除, 本站出于学习和科研的目的进行交流和讨论,如有侵犯原作者的版权, 请来信告知,我们将立即做出整改,并给予相应的答复,谢谢合作!

中国磨削网

QQ|Archiver|手机版|小黑屋|磨削技术网 ( 苏ICP备12056899号-1 )

GMT+8, 2025-1-10 00:32 , Processed in 0.146035 second(s), 23 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表