+高级检索
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于高分辨率模态柔度矩阵桥梁损伤识别的少量传感器方法  PDF

  • 马宏伟 1,2
  • 郑晓杰 3,4
  • 夏子立 5
  • 闫禹 5
  • 聂振华 3,4
1. 东莞理工学院 生态环境与建筑工程学院,广东 东莞 523808; 2. 广东省城市生命线工程智慧防灾与应急技术重点实验室,广东 东莞 523808; 3. 暨南大学 力学与建筑工程学院,广东 广州 510632; 4. “重大工程灾害与控制”教育部重点实验室,广东 广州 510632; 5. 港珠澳大桥管理局,广东 珠海 519060

中图分类号: TN911.6

最近更新:2024-07-29

DOI: 10.16339/j.cnki.hdxbzkb.2024074

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

为了降低监测系统成本和提高桥梁损伤识别精度,本文利用少量传感器获得高分辨率模态振型. 在考虑路面粗糙度的情况下,基于柔度变化率进行损伤识别. 在不同路面粗糙度水平下,对移动车辆荷载作用下桥梁上所测位移响应进行特征正交分解,再将得到的主成分矩阵中高频动力分量信息过滤,获得高分辨率的模态振型. 利用少量传感器方法得到的高分辨率模态振型计算柔度变化率,将其作为桥梁损伤定位的指标. 为了验证该方法的有效性,对一简支梁进行模拟和试验. 结果显示,该方法在考虑不同粗糙度水平情况下,利用少量传感器得到高分辨率模态振型的柔度变化率能准确定位损伤.

温度变化、环境腐蚀和材料性能退化等因素,导致桥梁结构性能降低甚至出现损伤,严重影响交通正常运营和人民生命财产安全. 因此,桥梁结构的安全监测和损伤识别尤为重要,已成为当前研究的热

1.

由于结构动力参数变化要比静力参数包含更多结构状态信息,所以多数研究基于结构的动力参数开展损伤识别. 而基于模态参数的方法原理简单、计算方便,在损伤识别研究中得到广泛关注,其中以“动力指纹”类为例,模态振

2、固有频3、曲率模4、柔度矩5以及协调模态置信因6等都具有良好的损伤识别效果. 聂振华7利用可视化的图像形式直观地展现结构基于模态的动力指纹指标,验证了该方法的可靠性.

结构在发生损伤后,会导致刚度降低,柔度增加. 柔度具有对损伤更加敏感的特点,Raghavendrachar

8利用固有频率、模态振型以及柔度三个损伤指标,对一座三跨混凝土桥进行了损伤识别,通过对比各损伤指标的识别效果,得到基于柔度的损伤识别指标的灵敏性明显高于基于固有频率以及模态振型的损伤识别指标. 而且对于实际工程中一些大型复杂结构,仅能较精确地测得其低阶的模态参数,而柔度矩阵与固有频率的平方成反比,因此只需测得结构的前几阶模态参数即可得到精度较好的柔度矩阵.基于柔度矩阵的以上优点,国内外学者以柔度作为损伤指标对结构的损伤识别问题进行了研究.

无论是曲率模态还是模态柔度等动力指纹方法,均需要布置大量传感器才能获得较高分辨率的指标效果,高昂的成本导致其在实际工程中难以被推广应用. 而现有方法只能提供稀疏的、低分辨率的模态振型,难以满足基于模态振型的损伤定位. 因此利用少量传感器获取高分辨率模态振型具有重要意义. 一些学者对高分辨模态振型做了研究. Yang

9提出了利用移动车辆的动力响应识别桥梁高分辨率振型的间接方法. Eshkevari10将二阶盲源识别与频响函数和集成经验模态分解相结合,提出了基于车载传感器数据的桥梁高分辨率模态识别方法. 然而,这些方法仍有缺点. 在间接方法中,车辆的响应不仅包括桥梁的动力信息,还包括车辆本身的动力信息以及路面粗糙度增强的噪声,这给模态参数的识别带来了困难.

针对上述问题,在考虑不同路面粗糙度的前提下,本文提出在移动车辆荷载作用下,通过对桥梁上少量传感器测得的位移信号进行特征正交分解(Proper Orthogonal Decomposition,POD)的方法来提取高分辨率的模态振型,然后利用高分辨率柔度变化率作为损伤指标进行损伤定位.

1 基本理论

1.1 桥梁高分辨率模态振型识别的少量传感器方法

一个移动荷载匀速通过长度为l的均质简支梁,如图1所示,其动力微分方程为:

fig

图1  移动荷载作用下的简支梁模型

Fig.1  Simple bridge model under a moving load

ρ2ut2+cut+EI4ux4=fδx-vt (1)

式中:uρcEI分别表示简支梁的竖向位移、密度、阻尼系数、弹性模量、惯性矩;fv分别表示竖向移动荷载和移动速度.

采用振型叠加法,振型叠加的公式表示为:

u(x,t)=n=1qn(t)φn(x) (2)

式中:qn表示第n阶振型坐标;φn表示第n阶振型.

简支梁在移动荷载作用下的竖向位移响

11表示为:

ux,t=2fρln=11ωn21-Sn22+2ξnSn21-Sn2sinΩnt-2ξnSncosΩnt+e-ξnωnt2ξnSncosωdnt+Sn1-ξn22ξn2+Sn2-1sinωdntsinnπxl (3)

式中:ωn为第n阶结构固有频率;Sn=πv/lωnSn为无量纲参数;ζn为第n阶阻尼比;Ωn为移动荷载的第n阶激励频率,Ωn=nπv/lωdn为阻尼振动频率.

对竖向位移响应的解析解进行POD分解,由于大多数桥梁由钢、混凝土组成,阻尼系数较小且材料密度较大,式(3)中阻尼比ζn是一个无穷小量,故可以忽略阻尼影响,得到无阻尼主成分公

11为:

C=UV=2fρl1ω12sinπvtl-πvlω1sinω1t1ω22sin2πvtl-πvlω2sinω2t1ωm2sinmπvtl-πvlωmsinωmtT
(t=t1,,tk) (4)

式中:U为位移信号响应矩阵;V为位移响应矩阵的协方差矩阵的特征向量.

式(4)可以看出,各阶主成分包含两部分参量:sinmπvtl为第m阶模态振型,πvlωmsinωmt为第m阶动力分量信息. 为了获得模态振型,需将动力分量信息过滤.

Nie

12提出了移动平均滤波器,用于直接滤除特定频率的信号. 滤波器公式表达为:

C(i,n)¯=fnfsj=i-K/2j=i+K/2C(j,n),n=1,,m (5)

式中:fn为桥梁的第n个频率;fs为采样频率;K是由fs/fn计算的滤波长度.

利用低通滤波器把高频动力分量信息过滤可得到第m阶模态振型φm.

该方法通过少量传感器即可获得高分辨率模态振型,如布置三个传感器即可获得前3阶高分辨率模态振型.

1.2 柔度变化率

Pandey

13提出了柔度变化率的结构损伤识别方法,该方法后期被证明是可靠的. 但由于传感器用量的限制,传统模态识别方法只能得到与传感器数量一致的稀疏振型向量,当结构损伤位于两个传感器之间时,柔度矩阵方法无法精准定位结构损伤,从而制约了该方法在工程中的应用. 而本文方法得到的高分辨率模态振型,为柔度矩阵方法在精准定位损伤上提供了新的解决途径. 结构的柔度矩阵F近似表示为:

F=i=1n1ωi2φiφiT (6)

式中:ωi为第i阶固有频率;φi为第i阶振型.

Ko

14提出柔度变化率,表示为:

D=diag(Fu)-diag(Fd)diag(Fu) (7)

式中:FuFd分别为结构损伤前后柔度矩阵,diag()为取矩阵对角值.

2 数值模拟

2.1 仿真模型及损伤工况

有限元模型中,使用ANSYS中plane42单元建立简支梁模型. 一个质量为1 000 kg的移动质量块从简支梁的一端移动到另外一端,速度设置为1 m/s. 采用矩形截面钢材,矩形截面尺寸为0.1 m×0.2 m,密度为7 850 kg/m3,弹性模量为210 GPa,泊松比为0.3,全梁长l=20 m,将简支梁模型沿长度方向划分400个单元,高度方向划分10个单元. 在实际情况中路面的粗糙度会显著影响损伤指标的识别效果,为了更加符合实际应用情况,分别考虑了两个不同粗糙度水平A、B的影响. 在简支梁上均匀布置7个传感器,采集位移信号. 采样频率设置为200 Hz. 简支梁模型如图2所示.

fig

图2  简支梁模型

Fig.2  Simple beam model

路面粗糙度是影响系统分析的关键因素,较差的路面会显著增加车辆荷载和桥梁响应的振动幅度. 基于国际标

15提供的路面功率谱函数,利用谐波叠加16模拟路面粗糙度分布,公式如下:

G(ni)=G(n0)(nin0)-2 (8)
r(x)=i=1N2G(ni)cos(2πnix+θi) (9)

式中:ni为每单位长度的空间频率;n0为0.1 cycle/m;G(n0)为位移功能密度函数值;θi为0~2π之间的随机相位角.

在梁的0.4l 和0.7l 位置处设置两个单损,通过沿梁高方向从桥底到中性面依次去除相应的单元来设置不同程度的损伤,截面损失率分别为10%、20%和30%,以及在梁的0.4l和0.7l位置处设立两个损伤,工况设置如表1所示. 为验证本方法的有效性,考虑了两个不同粗糙度水平的影响. 以A粗糙度水平下工况1为例,图3为简支梁的位移响应,其中的相对位置为无量纲单位.

表1  工况设置
Tab.1  Scenarios setting
工况损伤位置损伤程度/%
1 0 0
2 0.4l 10、20、30
3 0.7l 10、20、30
4 0.4l、0.7l 20
5 0.4l、0.7l 40
fig

图3  A粗糙度水平下工况1的位移响应

Fig.3  Displacement response at Scenario 1 of roughness A

2.2 数值模拟结果

在路面粗糙度的影响下,通过对桥梁上布置的7个传感器获取的位移响应进行POD处理,再通过低通滤波器将动力分量信息过滤,即可获得简支梁的前7阶模态振型.

特征值表示其对应分量在振动响应中所占据的能量,所以利用主成分分析中第n个分量的贡献率来估计能量意义上对应的贡献率(Contribution Ratio,CR

17-18. 因此,第n阶模态贡献率表示为:

CR=λnj=17λj,n=1,,7 (10)

式中:λn对应协方差矩阵C的第n个分量的特征值.

利用式(10)计算前7阶模态贡献率,从表2可看出,前3阶的模态参与率达到了99.99%. 因此,从节约成本和低阶就能获取较完备的柔度矩阵的角度来考虑,只在简支梁上均匀布置3个传感器,分别对应图2的1号、4号和7号传感器的位置. 前3阶的固有频率如表3所示.

表2  模态贡献率
Tab.2  Contribution ratio of each mode ( % )
阶数粗糙度1st2nd3rd4th5th6th7th
A 99.593 03 0.389 62 0.015 47 0.001 53 0.000 25 0.000 06 0.000 01
B 99.576 88 0.405 87 0.015 34 0.001 56 0.000 25 0.000 06 0.000 01
表3  简支梁的固有频率
Tab.3  Natural frequencies of the simple beam model
模态阶数1st2nd3rd
频率/Hz 1.14 4.57 10.26

传统方法利用3个传感器只能获得3个振型数据,基于3个数据对工况3的损伤因子柔度变化率进行计算,如图4所示. 可以看到,由于振型数据点有限,3个不同损伤程度的工况均不能精确对损伤位置进行定位. 因此采用传统模态识别得到的柔度变化率方法在传感器不足的情况下失效,无法推广到实际工程中.

fig

图4  模态分析损伤识别结果

Fig.4  Damage identification results of modal analysis

以工况1为例,在车辆速度为1 m/s时,采样点个数为4 000,对桥梁位移信号进行POD处理,获得模态振型个数也为4 000,因此得到的桥梁模态振型具有较高的空间分辨率. 在A、B两个粗糙度水平下,3个传感器获得的前3阶高分辨率模态振型归一化后如图5所示.

fig

(a)  A粗糙度前3阶振型

fig

(b)  B粗糙度前3阶振型

图5  A、B粗糙度识别的高分辨率模态振型

Fig.5  The identified high-resolution modal shapes of roughness A and B

利用A粗糙度水平的高分辨率模态振型计算柔度变化率,图中红色三角形表示损伤位置. 损伤识别结果如图6所示,可以看到0.4l与0.7l损伤位置有峰值,可以对损伤进行精确定位.

fig

(a)  A粗糙度水平0.4l处10%、20%和30%损伤

fig

(b)  A粗糙度水平0.7l处10%、20%和30%损伤

fig

(c)  A粗糙度水平0.4l和0.7l处20%和40%损伤

图6  A粗糙度水平下损伤识别结果

Fig.6  Damage identification results of roughness A

增大路面粗糙度水平到B后,分别计算各个工况的柔度变化率损伤指标值. 如图7(a)、(b)所示,可以看到在B粗糙度水平下,在0.4l与0.7l两个单损工况下,在损伤位置均有峰值,可以对损伤进行精确定位;如图7(c)所示,在0.4l与0.7l两个双损工况下,损伤位置均有峰值,可以对损伤进行精确定位.

fig

(a)  B粗糙度水平0.4l处10%、20%和30%损伤

fig

(b)  B粗糙度水平0.7l处10%、20%和30%损伤

fig

(c)  B粗糙度水平0.4l和0.7l处20%和40%损伤

图7  B粗糙度水平下损伤识别结果

Fig.7  Damage identification results of roughness B

以上结果均表明,在分别考虑了两个不同路面粗糙度水平影响的情况下,基于高分辨率模态振型计算得到的柔度变化率损伤因子仍然能够对损伤进行精确定位,而且随着损伤程度的增大,各个工况的柔度变化率也在增大,可以定性比较各个工况损伤程度的大小. 同时也可以看到,由于边界效应的影响,各个工况在梁两端的曲线会有比较大的波动.

2.3 噪声的鲁棒性

在土木工程结构健康监测中,不可避免地受到环境中各种噪声的影响. 为了验证该方法的噪声鲁棒性,把数值模拟中得到位移响应添加70 dB和 90 dB的高斯白噪音来模拟环境因素的影响. 通过信噪比衡量噪声的水平,信噪比公式如下:

SNR=10lgx2(n)y2(n)=20lgASAN (11)

式中:x(n)表示原始信号;y(n)表示噪声;AS是信号x(n)的均方根;AN是信号y(n)的均方根;信噪比的单位为分贝(dB).

以A粗糙度水平单损和双损工况为例,损伤识别结果分别如图8图9所示. 在添加70 dB、90 dB噪声后,除了图8(a)中10%损伤程度的峰值不太明显之外,其余工况在损伤位置处均出现峰值,同时可以看到由于边界效应和噪声的影响,梁两端曲线也存在比较大的波动,但是在损伤位置是存在最大峰值的,仍然能精确识别出损伤位置,因此该方法对噪声具有一定的鲁棒性.

fig

(a)  90 dB噪声下A粗糙度水平0.4l处10%、20%和30%损伤

fig

(b)  70 dB噪声下A粗糙度水平0.4l处10%、20%和30%损伤

图8  噪声影响下A粗糙度水平下单损识别结果

Fig.8  Single damage identification results of roughness A with noise

fig

(a)  90 dB噪声下A粗糙度水平双损

fig

(b)  70 dB噪声下A粗糙度水平双损

图9  噪声影响下A粗糙度水平下双损识别结果

Fig.9  Multiple damage identification results of roughness A with noise

3 试验研究

为了进一步验证本文提出的损伤识别方法的可行性,利用简支梁试验模型对其进行验证. 简支梁为一空心矩形钢梁,钢梁长度L=6 m,弹性模量E= 210 GPa,截面参数为200 mm×100 mm×3 mm(长×宽×厚),材料密度为7 850 kg/m3,小车模型质量为10.5 kg,采用直流电动机来牵引小车在轨道上匀速移动,通过改变电动机的转速来控制小车的移动速度. 沿梁长度方向均匀布置5个测点采集位移信号,从左端到右端依次标记传感器序号为1~5,本试验利用DH5922N动态数据采集系统进行信号采集,采样频率设置为500 Hz. 为了检验该方法对不同小车速度激励的敏感性,试验中分别设置小车速度为 0.5 m/s和0.75 m/s. 试验装置简图如图10所示.

fig

图10  试验装置简图

Fig.10  Simple model of experimental device

由于试验梁表面并不是完全理想化的光滑,会存在加工以及环境因素造成梁表面的不平度,所以以此来等同于粗糙度,研究其对该方法的影响. 试验装置图如图11所示.

fig

图11  试验装置图

Fig.11  Experimental device

在梁0.68l位置处设立一处损伤,在梁0.38l、0.68l位置处设置两处损伤. 利用切割机对简支梁底部进行切割来制造损伤,在0.68l处设置一处损伤,损伤深度为7 mm,在0.38l处设置一处损伤,损伤深度分别为4 mm和7 mm,试验损伤工况如表4所示.

表4  试验工况
Tab.4  Experimental scenarios
工况损伤位置损伤深度/mm小车速度/(m·s-1
1 0 0 0.5、0.75
2 0.68l 7 0.5、0.75
3 0.38l、0.68l 4、7 0.5、0.75
4 0.38l、0.68l 7、7 0.5、0.75

以工况1中0.75 m/s小车速度为例,位移传感器采集到的位移响应如图12(a)所示. 傅里叶频谱图如图12(b)所示.

fig

(a)  位移响应

fig

(b)  傅里叶频谱

图12  0.75 m/s位移响应及傅里叶频谱

Fig.12  Displacement responses and Fourier

spectrums of 0.75 m/s

以工况1中0.75 m/s小车速度为例,5个传感器采集位移信号并进行POD分析,可以得到前5阶高分辨率模态振型,如图13所示,可以看到识别的第4和第5阶振型效果较差,这是由于在测量响应中与模态相关的信息太弱,因此不能很好地识别出高阶模态振型.

fig

图13  0.75 m/s工况1前5阶振型

Fig.13  First 5 mode shapes of scenarios 1 with 0.75 m/s

计算各个工况下的损伤因子柔度变化率,红色三角形标记点表示损伤位置,在小车速度为0.5 m/s时,从图14(a)可以看出,在0.68l单损工况下,在损伤位置有峰值,精确定位了损伤;如图14(b)所示,在0.38l与0.68l两个双损工况下,在损伤位置均有峰值,说明该方法对多损伤工况同样有效,而且随着损伤程度的增大,柔度变化率也在增大.

fig

(a)  单损工况

fig

(b)  双损工况

图14  0.5 m/s的损伤识别结果

Fig.14  Damage identification results of 0.5 m/s

为了检验小车速度对该方法的影响,将小车速度增加到0.75 m/s时,损伤识别结果如图15所示. 从图15(a)可以看出,在0.68l单损工况下,在损伤位置有峰值,仍然可以对损伤进行精确的定位;如图15(b)所示,在0.38l与0.68l两个损伤位置的双损工况下,在损伤位置均有峰值,均可对损伤进行精确的定位,而且随着损伤程度的增大,柔度变化率也在增大.

fig

(a)  单损工况

fig

(b)  双损工况

图15  0.75 m/s的损伤识别结果

Fig.15  Damage identification results of 0.75 m/s

为了进一步检验传感器数量减少时该方法的有效性,将上述试验中5个传感器数量减少到3个(分别取图10中的1号、3号和5号传感器),利用POD处理位移响应得到高分辨率模态振型值,计算柔度变化率损伤指标.在0.5 m/s与0.75 m/s的速度下,损伤识别结果分别如图16~图17所示,无论是在单损还是双损工况下,在损伤位置处均存在峰值,成功地对损伤位置进行精确定位.

fig

(a)  单损工况

fig

(b)  双损工况

图16  0.5 m/s时3个传感器的损伤识别结果

Fig.16  Damage identification results of 0.5 m/s with three sensors

fig

(a)  单损工况

fig

(b)  双损工况

图17  0.75 m/s时3个传感器的损伤识别结果

Fig.17  Damage identification results of 0.75 m/s with three sensors

4 结 论

本文在考虑了桥面粗糙度影响的情况下,获得移动荷载作用下桥梁少量传感器的位移响应,通过对位移响应进行特征正交分解(POD),得到高分辨率模态振型,并以此计算损伤因子柔度变化率,对桥梁进行损伤定位. 本文主要结论如下:

1)基于少量传感器的柔度变化率损伤指标能够精准识别损伤位置,而且在减少传感器的数量时,仍然具有很好的损伤识别效果.

2)该方法克服了传统模态方法只能获得稀疏振型的缺点,并且考虑了不同路面粗糙度水平和噪声对该方法的影响,更加适合实际工程.

3)利用少量传感器就能实现损伤识别,大大降低了成本.

参考文献

1

马宏伟聂振华桥梁安全监测最新研究进展与思考[J].力学与实践2015372):161-170 [百度学术] 

MA H WNIE Z HRecent advances and review of bridge safety monitoring[J].Mechanics in Engineering2015372):161-170(in Chinese) [百度学术] 

2

SHIH H WTHAMBIRATNAM D PCHAN T H TDamage detection in slab-on-girder bridges using vibration characteristics[J].Structural Control and Health Monitoring20132010):1271-1290 [百度学术] 

3

WICKRAMASINGHE W RTHAMBIRATNAM D PCHAN T H Tet alVibration characteristics and damage detection in a suspension bridge[J].Journal of Sound and Vibration2016375254-274 [百度学术] 

4

PANDEY A KBISWAS MSAMMAN M MDamage detection from changes in curvature mode shapes[J].Journal of Sound and Vibration19911452):321-332 [百度学术] 

5

PANDEY A KBISWAS MExperimental verification of flexibility difference method for locating damage in structures[J].Journal of Sound Vibration19951842):311-328 [百度学术] 

6

HUNT D LApplicrion of an enhanced coordinate modal assurance criteria[C]//10th International Modal Analysis ConferenceSan Diego,Califomia. CASociety for Experiment Mechanics, Inc1992 [百度学术] 

7

聂振华程良彦马宏伟基于结构动力特性的损伤检测可视化方法[J].振动与冲击20113012):7-13 [百度学术] 

NIE Z HCHENG L YMA H WVisualization method for structural damage detection based on its dynamic characteristics[J].Journal of Vibration and Shock20113012):7-13(in Chinese) [百度学术] 

8

RAGHAVENDRACHAR MAKTAN A EFlexibility by multireference impact testing for bridge diagnostics[J].Journal of Structural Engineering19921188):2186-2203 [百度学术] 

9

YANG Y BLI Y CCHANG K CConstructing the mode shapes of a bridge from a passing vehicle: a theoretical study[J]. Smart Structures and Systems2014135):797-819 [百度学术] 

10

ESHKEVARI S SMATARAZZO T JPAKZAD S NBridge modal identification using acceleration measurements within moving vehicles[J].Mechanical Systems and Signal Processing2020141106733 [百度学术] 

11

NIE Z HSHEN Z FLI Jet alOutput-only complete mode shape identification of bridges using a limited number of sensors[J].Mechanical Systems and Signal Processing2022178109246 [百度学术] 

12

NIE Z HNGO TMA H WReconstructed phase space-based damage detection using a single sensor for beam-like structure subjected to a moving mass[J].Shock and Vibration201720171-20 [百度学术] 

13

PANDEY A KBISWAS MDamage detection in structures using changes in flexibility[J].Journal of Sound and Vibration19941691):3-17 [百度学术] 

14

KO J MSUN Z GNi Y QMulti-stage identification scheme for detecting damage in cable-stayed Kap Shui Mun Bridge[J].Engineering Structures2002247):857-868 [百度学术] 

15

Mechanical vibration-Road surface profiles-Reporting of measured dataISO 8608:2016[S]. LondonBSI Standards Limited20161-48 [百度学术] 

16

王红岩王钦龙芮强车辆行驶路面的数字化建模方法研究[J].兵工学报2016377):1153-1160 [百度学术] 

WANG H YWANG Q LRUI Qet alResearch on digitized modeling method of riding road of vehicle[J]. Acta Armamentarii2016377):1153-1160(in Chinese) [百度学术] 

17

LIANG Y CLEE H PLIM S Pet alProper orthogonal decomposition and its applications:part Ⅰ:theory[J].Journal of Sound and Vibration20022523): 527-544 [百度学术] 

18

WU C GLIANG Y CLIN W Zet alA note on equivalence of proper orthogonal decomposition methods[J].Journal of Sound and Vibration20032655):1103-1110 [百度学术] 

作者稿件一经被我刊录用,如无特别声明,即视作同意授予我刊论文整体的全部复制传播的权利,包括但不限于复制权、发行权、信息网络传播权、广播权、表演权、翻译权、汇编权、改编权等著作使用权转让给我刊,我刊有权根据工作需要,允许合作的数据库、新媒体平台及其他数字平台进行数字传播和国际传播等。特此声明。
关闭