专利名称：一种机床固定结合部动力学参数的识别方法
技术领域：
本发明涉及机床常见固定结合部(螺栓联接)动力学参数的识别方法， 特别是涉及对结合部刚度和阻尼参数的识别方法。
背景技术：
自从20世纪60年代，人们认识到机床结构中结合部对整机结构的重 要影响以来，研究者就对结合部的动力学特性和参数识别方法进行了广泛 而深入的研究，取得了令人瞩目的成果。但是，结合部的属性参数不能精 确控制，使得仍有许多问题有待更深入的研究。
无论是何种类型的结合部，其结合均属于"柔性结合"，这已经是一个
不争的事实。所谓"柔性结合"是指当结合部受到外加复杂动载荷作用时，
结合面间会产生多自由度、有阻尼的微幅震动(即变化的微小相对位移或 转动)，从而使结合部表现出既有弹性又有阻尼，既储存能量又耗散能量 的"柔性结合"的本质及其特性。因此，结合部的动力学特性中，刚度特性和 阻尼特性对机床结构的动态性能影响是显著的，而其质量特性对机床结构 的动态特性影响甚微。研究结合部的动力学特性就是对结合部刚度和阻尼 进行辨识。
90年代以前，建立结合部动力学模型最常用的方法是将结合部用一系 列孤立的粘弹单元来表示，如图1所示，亦即用这一系列粘弹单元来模拟
整个结合面动态特性。但该方法存在以下两个缺陷
1) 用该方法建模识别的结合部参数只适用于特定结构，当结合部尺寸 形状发生改变时，结合面参数移植困难，不具有通用性。
2) 这种模型中，仅考虑各自由度法向的刚度特性和阻尼特性，而忽略
了各自由度之间的耦合特性，更没有考虑结合部切方向的力学特性。
4青岛大学的张杰、童忠钫基于有限元思想，又提出了一种新的"理想结 合部"的"理想结合面元"模型，并进一步提出了一种基于"结合面元"的结合
部动力学模型，如图2所示(O-XYZ为整体坐标系，o'-x'r'为单元局部
坐标系)。这种模型假定各种条件下的结合部(包括不同的材料，预紧力， 表面粗糙度，加工方法，介质等)动态特性，均可用理想结合部的动态特 性与之等效。所谓"理想结合部"是指结合部之间的比压在整个结合部上为常 数，结合部上各点均匀接触，并在所有的接触点上具有相同力学性质的结 合部。将理想结合部进行有限元划分而得到理想结合部有限单元。但是该 模型的条件过于苛刻，对于接触面尺寸较大的结合部有一定的局限性。在 模型参数确定时比压需要准确确定，这是比较困难的。
确定结合部等效动力学模型之后，还需要确定其模型中各个自由度的
弹簧v阻尼器的等效动力学参数。
早期研究者通过直接测试的方法识别结合部参数，即直接测试结合部
处的力和位移(响应)之间的关系来研究结合部的刚度特性；通过反复加 载卸载，测试结合部的力和位移(响应)的迟滞回线，而用迟滞回线所包 含的面积大小来度量结合部处的阻尼特性。这种识别方法虽然很直接，物 理意义明显，但是通常结合部的刚度特性和阻尼特性都是非线性的，这就 使得所识别的参数与实际相差很大，不能很好的模拟结合部复杂的动态特 性。
对于平面联结结合部的刚度和阻尼特性，基尔萨诺娃、科巴赫、索恩 利等人进行了深入的研究，提出了一种"螺栓结合部刚度计算法"。假设接触
面上的法向位移义与平均接触压力&的指数函数关系式为；1 = "/^("/^
是由结合面的加工方法和材料等决定的常数，可从经验数据表中查取)，从 而，结合面的法向刚度为^ =丄尸 1,; 二平面接触的结合面在平均剪切力
rf义 "附
K作用下的切向位移5可表示为5二~&， ^为结合面剪切柔度， 为剪切动刚度。
5日本京都大学的吉村允孝对机床结合部的研究表明，尽管结合部接触 面积不同，但是只要结合部平均接触压力相同，单位面积结合部的动态性 能数据是相同的，并总结得到了结合部单位面积动态性能数据表。 一旦结 构确定后，相应的结合部也就确定了，这时只需根据结合面的条件从数据 表中查出相应的单位面积的特性数据，结合面的刚度和阻尼就可以通过对 结合面积积分获得，这就是"吉村允孝积分法"。
后期研究者又提出了一种基于机械阻抗的频响函数试验识别法。其基 本思想是利用子结构和整体结构的频响函数(FRF)提取结构结合部的动 力学参数，并考虑了测试频响函数时测试误差的影响。但是，在计算结合 部动态参数时需要对子结构的频响函数矩阵求逆，而通过矩阵求逆运算后， 因噪声干扰引入的微小的测量误差必然导致很大的识别误差。虽然之后的 研究者对这个问题提出了解决方法，在一定程度上提高识别精度。然而， 无论那种方法，在识别结合部参数时，均只是测试各子结构的频率响应函 数，并没有任何有关结合部的特性参数(如法向预载荷、表面粗糙度、金 属材料、表面的加工方式等)。因此，所识别的参数再准确，也只能适用 于正在识别的动力系统，没有普适性。
发明内容
本发明的目的在于提供一种机床固定结合部动力学参数的识别方法， 充分考虑了多个关键自由度方向的阻尼和刚度及其相互耦合关系，避免了 矩阵求逆所带入的二次误差，精度高。
一种机床固定结合部动力学参数的识别方法，具体为以结合部单元 结构阻尼矩阵q和结合部单元刚度矩阵《为参数变量，求取使得 (- 2M + ^C趋近于单位矩阵的最优解C，和《，其中质量矩阵M为 子结构质量矩阵M。，阻尼矩阵C由结合部单元结构阻尼矩阵C，和子结构阻
6尼矩阵C。组装而成，刚度矩阵《由结合部单元刚度矩阵《与子结构刚度矩 阵K。组装而成，//(W)为角频率w的频率响应函数，/为虚数单位；
所述结合部单元结构阻尼矩阵q和结合部单元刚度矩阵&的设定形式
分别考虑了结合部单元的各节点在三个平动自由度方向的阻尼和刚度，且 満足C,/g&， g为结构损耗因子。
一种机床固定结合部动力学参数的识别方法，具体为以结合部单元 结构阻尼矩阵C,和结合部单元刚度矩阵《为参数变量，求取使得 (- 2M + C + iO//(w)趋近于单位矩阵的最优解C,和《，其中质量矩阵M为子
结构质量矩阵M。，阻尼矩阵c为结合部单元结构阻尼矩阵c;的扩展矩阵，
刚度矩阵《由结合部单元刚度矩阵《与子结构刚度矩阵K。组装而成，//(w) 为角频率^的频率响应函数，/为虚数单位；
所述结合部单元结构阻尼矩阵C,和结合部单元刚度矩阵&的设定形式
分别考虑了结合部单元的各节点在三个平动自由度方向的阻尼和刚度，且 满足C,dg&， g为结构损耗因子。
所述结合部单元结构阻尼矩阵q和结合部单元刚度矩阵《均为 Px")x(h")维矩阵，n为结合部单元的节点个数。结合部单元的节点个数 可为8或16或18，选择的节点个数越多，精度越高，但是计算复杂程度也 提高，因此优选8节点。
本发明以有限元理论和振动力学理论为基础，建立了一种结合部的"参 数化"动力学模型，该模型充分考虑到了结合部的各种影响因素以及动力学 参数间的耦合关系；然后，以此为基础，利用多自由度动态系统阻抗矩阵与频响矩阵的互逆性这一基本性质，采用理论模态计算与实验模态测试相 结合的方法对结合部刚度和阻尼进行识别，识别精度高，可以更加准确表 征结合部更丰富的动力学特性。
本发明首先获得有限元分析的理论值(实体子结构的刚度、阻尼矩阵) 与实验模态分析得到的实验值(整体结构的频响函数)，然后利用优化设 计方法，以要求的结合部参数(在研究中指结合部刚度矩阵、阻尼矩阵) 为设计变量，以阻抗矩阵和频响函数之积与单位矩阵之差最小为优化设计 目标，通过多次迭代得到的设计变量值即为识别出的结合部参数。
本发明所建立的新的结合部动力学模型的理论计算结果与实验比较误 差不超过7%，与现有技术相比，精度得到了极大的提高，说明了我们所建 立的模型与求解方法的先进性。
图1为结合部弹簧阻尼器模型示意图； 图2为矩形理想结合部单元示意图3为螺栓联结形式示意图，图(3a)为"线形"式，图(3b)"阵列"
式；
图4为结合部螺栓联结单元模型示意图,图(4a)"线形式"联结单元模 型,图(4b)"阵列式"联结单元模型；
图5为结合部有限单元模型示意图6为结合部单元受力模型示意图； 图7为结合部动力学模型示意图； 图8为实验装置示意图9为实验试件示意图。
具体实施例方式
本发明的具体实施步骤如下
81) 对常见机床固定结合部的各种形式进行调研，给出固定结合部 有限单元模型的划分准则，提取其本质属性，建立"参数化"结合部单元动力 学模型；
2) 对子结构进行有限元分析，提取子结构质量矩阵M。和刚度矩阵
&，子结构一般都是由钢铁加工而成，这种材料本身的材料阻尼是很小的，
在这种情况下，机床结构的内阻尼的90%来自于机床结合部，故子结构的
阻尼矩阵C。与结合部阻尼矩阵q相比可忽略；
3) 根据"参数化"结合部单元动力学模型的节点对应关系，利用有
限元矩阵组装理论，将子结构刚度矩阵《。和结合部单元的刚度矩阵《进行
组装得到总刚K;结合部质量特性对机床结构的动态特性影响甚微可以忽 略，因此总质^ = ^。，阻尼矩阵c是由结合部结构阻尼矩阵q扩展得到， 即C中结合部8节点位置所对应的分块矩阵为C，，其余元素为零；^和C,即 为待求解动力学参数；
4) 通过对真实结构提取出来的包含结合部单元的縮小模型进行模 态实验，测定螺栓联结的整个结构的频响函数，根据有限元理论组装得到 的含有待求解参数的整个结构的阻抗矩阵与实测频响函数的互逆性，利用 最小二乘法，通过多次迭代求解设计变量值。
5) 结合部不可能脱离机械结构系统而存在，验证求解得到的结合 部动力学参数的有效性，必须把相应参数带入整个结构系统中，对该整体 模型进行理论模态分析,并与实验模态分析结果进行比较验证。
1.调研，分析，建模
目前研究者们所建立的结合部动力学模型都没有充分体现结合部的具 体属性(结合部的材料、表面形貌、预紧力、是否存在介质(润滑油等)、 结合部尺寸与几何形状等)，将这种模型命名为"非参数化"结合部模型，这 也是限制目前结合部研究成果通用性的主要原因。
有鉴于此，我们提出建立一种能反映结合部具体属性的"参数化"结合 部动力学模型。这种模型所建立的结合部刚度K、阻尼C的表达式包含了 结合部具体属性，可以想象，只要结合部的具体参数一致，就可以使用相
9同的结合部动力学参数模型。所以结合部"参数化"模型具有通用性，可以为 机床结构的整体动力学模型的建立提供统-的理论和数据支持。
通过调研发现，机床固定结合部多为螺钉联接，且多为"线形"式和"阵 列"式两种形式，称螺钉联接位置为联结点，如图3所示。根据弹性力学
St.Venant局部影响原理和固定结合部应力仿真计算结果，对结合部的动力 学特性做如下假设。对于"线形"连接形式相邻两个螺钉之间的结合部动力 学特性，仅受这两个螺钉力学状态的控制，而与这两个螺钉之外的其他螺 钉的力学状态无关(如图4 (a)所示)；对于"阵列式"连接相邻四个螺 钉之间的结合部的动力学特性，仅受这四个螺钉力学状态的控制，而与这 四个螺钉之外的其他螺钉的力学状态无关(如图4 (b)所示)。
可见，不论是"线形式"连接还是"阵列式"连接，都可以将结合部单元 模型表示成图5所示的形式，其中I、 II分别是螺栓联结的上下子结构单 元，m是假想的结合部单元，1-5、 2-6、 3-7、 4-8是结合部单元对应的8个 节点。每」个结合部单元有8个节点，每个节点具有3个自由度(三个平 动自由度)，每个单元总共有24个自由度。结合部单元的运动将通过节点 1与节点5、节点2与节点6、节点3与节点7、节点4与节点8之间的相 对运动表现出来。只要能准确建立这些节点位移与节点受力之间的关系， 就等于建立了其动力学模型。
对于一种结合部，存在六种不同形式的动态力，即六个自由度的广义 力(如图6所示)。也可产生六个自由度上的阻尼力，并且具有不同的阻 尼损耗因子。这些动态力分别是Z方向的正向力F:， X、 Y方向上的剪切 力《,&，绕X、 Y轴的弯矩M^A^，以及绕Z轴的剪切扭矩M^。
对于结合部来说，它的质量是可以忽略的，在建立其动力学模型时， ^R考虑其弹性和阻尼特性。
设各节点位移为^，各节点所受的力为^， Hl，2,…,7，8; j=l,2,3。首 先推导结合部有限单元的刚度矩阵。如前所述，结合部的运动特性是通过 节点i与节点5、节点2与节点6、节点3,与节点7、节点4与节点8之间 的相对运动表现出来，而这些节点之间的相对运动可表示为(d!j- d5j)、
10(d2j-d6j)、 (d3j-d7j)、 (d4j-d8j)， j=l，2,3。在不考虑阻尼的情况下,根据刚 度影响系数法，则
|j《(dln- d5n) + |j《(d2n- d6n) + |;《(d3n- d6n) + |j《(d4n- d8n) = ^ (1)
其中，尺i为刚度影响系数，i,m=l、 2、 3、 4，表示节点；n,j=l、 2、 3， 表示自由度方向；《i具体物理意义为仅在节点m与节点(m+4)的n方向
产生单位相对位移，而相应在i节点j方向所需施加的力。
在平衡条件下，有/17 =-/5/ 、 /2y. =-/6/ 、 /3/ =-/7y. 、 /4广-/8J. (j=l,2,3)
令结合部单元节点位移向量为—(1 ， 12， 3 ，1,2,3 ", 1 ， 82， 3) 结合部单元节点弹性力向量为 = "1，,12，,13，721，"^22'^23".'^81'^82,^83) 因此，(l)式写成矩阵形式为 [A][D，]呵Fw]
由上述推导可知，结合部单元刚度矩阵[《]具有对称性和可分块性；其
/ ft r,、
分块形式为 ;一，，
I-A X 」24x24
其中K'为结合部单元刚度矩阵[&]进行"四等分块"后的左上角矩阵，利 用该性质，可以简化[《]的计算量，只要求得[^']12>;12，利用该性质即可得到 [&]24x24 。矩阵[]12><12中元素的具体排列形式如下所示。
ii结合部单元的阻尼通常为结构阻尼，结构阻尼的阻尼力^与振动位移 X成正比，相位比位移超前90。，结构阻尼系数77，则力^/77X。又^与刚度 A成正比，77 = gA， g为结构损耗因子(也称结构阻尼比，为一经验常数)。
从而力=勿^ = ^^。仿照结合部单元刚度矩阵[《]的建立过程，可以得到结 合部单元的结构阻尼力矩阵[F^]与振动位移[DJ之间的平衡方程为[D,H^]，其中[C,],[《]
所以，结合部单元的结构阻尼矩阵矩阵形式为[q]=
姊， -,.沐' —'.沐'姊'
乂24x24
根据上述结合部单元的受力状况分析，可以将其形象的表示为图7所
示的模型，本发明将其称之为结合部"参数化"动力学模型，该模型单元为八 节点六面体型单元，即该单元是具有六面体空间拓扑形状的假想单元，八
节点即为图示8个角点。图7中明确绘出了各个方向的刚度及其耦合关系， 阻尼只是在图中用C表示，其阻尼耦合关系同理于刚度。 2.识别方法原理说明
对于多自由度振动系统,其动力学微分方程为M[i3]+C[3]+《[D] = [F]。
其中M[i3]为惯性力[&]， c[z3]为阻尼力[i^]，〖[D]为弹性力[ig， [F]为外
加激励力，[D]为系统位移矩阵。
以该方程表示的系统阻抗矩阵(也称动刚度矩阵)为- 2M + / C + iT， 其逆矩阵为系统的频率响应函数矩阵H(w)。
所以可以得到 (-w2M + / C += ￡ ， E为单位矩阵
子结构一般都是由钢铁加工而成，这种材料本身的材料阻尼是很小的，
V V V V V w w <
& XI Kl & w &4 4 4
rl f— rl w. ^,. ，^
V V 1<2 & & &
12在这种情况下，机床结构的内阻尼的90%来自于机床结合部，故子结构的 阻尼矩阵Q与结合部阻尼矩阵G相比可忽略。系统阻尼力矩阵[d]，其中[c]为系统阻尼姬阵，由于系统阻尼仅考虑结合部的结
构阻尼，故[c]只与结合部单元结构阻尼矩阵c,相关。
所以动力学微分方程可简化为MC[Z)]+= [F]，其阻抗矩阵为 + C + K ，使其与试验采集到的频响函数矩阵H(w)的乘积趋近于单位矩 阵的最优解即可求得C,和《。式-w2M + C + K中的系统阻尼矩阵C是由结合 部单元结构阻尼矩阵C,扩展得到，若假设结合部8个节点编号在结构整体
有限元模型中的编号仍为1 8，则C矩阵形式可以表示为C-f^1 "]，即C中
、o o乂
结合部8节点位置所对应的分块矩阵为C,，其余元素为零；刚度矩阵《由 子结构刚度矩阵《。和结合部单元的刚度矩阵《组装得到；质量矩阵M = M。。 对于各子结构，质量矩阵M。'，阻尼矩阵(比例阻尼，可忽略)，刚度矩阵《。可 通过有限元理论获得；对于结合部，忽略其质量矩阵，其刚度矩阵&和阻 尼矩阵C,为本发明所求动力学参数。 3.实验*求解，验证
本实验系统由固定结合部试件(子结构i、 n、结合部m)、加速度 传感器w、连接螺钉、激励力锤vi (含力传感器v)、电荷放大器vn、
LMS-CADA-X信号采集及处理系统W1、 PC机IX、柔性悬挂绳、支撑架、
电缆屏蔽线等实验设备组成。将结合部试件用柔性绳自由悬挂，采用单点
激励多点测试(SIMO)的方法，对结合部模型进行频响测试，实验装置示 意图如8所示。
实验时，将结合部实验试件，如图9所示(子结构i、 n、结合部m)
悬垂，在某点用力锤VI进行激励，采用多个加速度传感器IV进行响应信号 采集，经LMS信号采集及处理系统VDI处理后，从PC机IX中提取整体结构的频响函数，阻抗矩阵可以利用有限元理论组装得到，利用多自由度动态 系统阻抗矩阵与频响矩阵的互逆性这一基本性质，即可求得结合部未知刚 度、阻尼参数。
图9所示的试件的基本参数如下
实验材料45#钢；
结合部尺寸60mmx31.5mmx29mm; 结构件尺寸190mmxl90mmxl00mm; 预紧力矩为45N-m;
实验划分结点数为72个； 两个螺栓连接，螺栓M10
将采用本发明所识别的结合部刚度矩阵& 、阻尼矩阵c,代入建立的结
合部动力学模型之中，并将该模型与构成结合部的各子结构的有限元动力 学模型组装，可得到含结合部的结构整体动力学模型，对该整体模型进行 理论模态分析,并与实验结果进行比较。比较结果如表1所示。
各阶模态振型实验结果 CHz)理论结果(Hz)误差
1偏航4094090%
2测翻4714514.2%
3俯仰8718532.1%
4左右平动167317655.5%
5前后平动213022786.95%
6上下平动256826473,%
表1各阶固有频率
从上表可以看出,理论与试验计算结果误差最大不超过7%。可见，采 用这种方法识别的结合部参数是具有很高精度的，能够很好的模拟结合部 的真实特性。
14将求解结果与"螺栓结合部刚度计算法"和"吉村允孝积分法"进行比 较，结果如表2所示.
各阶模态实验结果/Hz理论结果/Hz吉村允孝螺栓刚度计算
1偏航40940911087
2测翻471451385456
3俯仰871853755886
4左右平16731765391309
5前后平21302278394310
6上下平2568264718522430
表2各阶模态固有频率比较
表2表明,现有技术所述结合部动力学模型的理论计算结果与实验结果 之间的误差是比较大的,而我们所建立的新的结合部动力学模型的理论计算 结果与实验比较误差不超过7%，说明了我们所建立的模型以及求解方法的 先进性。
同样是基于频率响应函数的结合部参数识别技术研究，本方法却能够 很好的模拟结合部的真实结构。之所以会出现这种结果，是因为我们所建 立的结合部动力学模型，综合考虑到了结合部各个方向的刚度和阻尼特性， 并没有人为指定某一特定类型的结构来模拟结合部模型，充分考虑到了结 合部的各个参数及其耦合关系对结合部特性的影响。
如果把结合部的特性和参数作为一个"未知箱子"，那么，我们建立结 合部动力学模型时则是采用"黑箱建模"方法，即仅仅根据系统的输入输出数 据建模，即使对系统的结构和参数一无所知，也能够对求解模型有一个很 好的模拟。而相对于前人的"灰箱建模"来说，由于他们在建立模型时加入了 一定的经验假设，势必带来一定的偏差。
本项发明中的创新点如下 充分考虑了结合部各个自由度之间的耦合关系。
给出了固定结合部有限单元模型的划分准则"线形"式和"阵列"式。
15 建模精度与求解结果精度得到极大提高。
提出了结合部"参数化"动力学模型的概念，结合部"参数化"动力学模 型的通用性，为机床结构的动力学研究提出了一种新的思路，也为机床结 构的动态优化设计提供了一种全新的计算方法。
权利要求
1、一种机床固定结合部动力学参数的识别方法，具体为以结合部单元结构阻尼矩阵C1和结合部单元刚度矩阵K1为参数变量，求取使得(-ω2M+iωC+K)H(ω)趋近于单位矩阵的最优解C1和K1，其中质量矩阵M为子结构质量矩阵M0，阻尼矩阵C由结合部单元结构阻尼矩阵C1和子结构阻尼矩阵C0组装而成，刚度矩阵K由结合部单元刚度矩阵K1与子结构刚度矩阵K0组装而成，H(ω)为角频率ω的频率响应函数，i为虚数单位；所述结合部单元结构阻尼矩阵C1和结合部单元刚度矩阵K1的设定形式分别考虑了结合部单元的各节点在三个平动自由度方向的阻尼和刚度，且满足C1＝igK1，g为结构损耗因子。
2、 一种机床固定结合部动力学参数的识别方法，具体为以结合部单 元结构阻尼矩阵q和结合部单元刚度矩阵X,为参数变量，求取使得 (- 2M + C + iQ//(w)趋近于单位矩阵的最优解q和《，其中质量矩阵M为子结构质量矩阵M。，阻尼矩阵c为结合部单元结构阻尼矩阵q的扩展矩阵，刚度矩阵K由结合部单元刚度矩阵《与子结构刚度矩阵《。组装而成，i/(w) 为角频率^的频率响应函数，/为虚数单位；所述结合部单元结构阻尼矩阵C,和结合部单元刚度矩阵&的设定形式 分别考虑了结合部单元的各节点在三个平动自由度方向的阻尼和刚度，且 满足C,/g&， g为结构损耗因子。
3、 根据权利要求1或2所述的机床固定结合部动力学参数的识别方 法，其特征在于，所述结合部单元结构阻尼矩阵q和结合部单元刚度矩阵《 均为(3xw)x(3xw)维矩阵，n为结合部单元的节点个数。
4、 根据权利要求3所述的机床固定结合部动力学参数的识别方法，其特征在于，所述结合部单元为八节点六面体型单元。
5、 根据权利要求4所述的机床固定结合部动力学模型参数的识别方 法，其特征在于，所述结合部单元结构阻尼矩阵q和结合部单元刚度矩阵《 均为24x24维矩阵。
6、 根据权利要求5所述的机床固定结合部动力学参数的识别方法，其 特征在于，所述八节点六面体型单元内包括两个相邻联结点。
7、 根据权利要求5所述的机床固定结合部动力学参数的识别方法，其 特征在于，所述八节点六面体型单元内包括四个相邻联结点。
全文摘要
本发明提供了一种机床固定结合部动力学参数的识别方法，将结合部单元结构阻尼矩阵和结合部单元刚度矩阵作为参数变量，以位移阻抗矩阵和位移频响函数之积与单位矩阵之差最小为优化设计目标，通过多次迭代得到的参数变量值即为识别出的结合部参数。本发明考虑了多个关键自由度方向的阻尼和刚度及其相互耦合关系，避免了矩阵求逆所带入的二次误差，精度高，更加准确表征结合部更丰富的动力学特性。
文档编号G01N19/00GK101458205SQ20081019682
公开日2009年6月17日 申请日期2008年8月29日 优先权日2008年8月29日
发明者强 徐, 斌 李, 毛宽民 申请人:华中科技大学