更全的杂志信息网

障碍物对丙烷-空气爆炸火焰加速的影响*

更新时间:2016-07-05

燃气泄漏导致的火灾和爆炸事故层出不穷,尤其当空间内存在障碍物时,爆炸产生的冲击波对周围环境会造成更严重的损坏。

有障碍物管道内火焰传播的研究有:Chapman等[1]首先对有障碍物存在的管道内火焰传播规律进行研究,发现当火焰传播至障碍物时,火焰速度明显增大;Dobashi[2]、Masri等[3]、Ibrahim 等[4]定性分析了不同障碍物结构对爆燃压力所产生的影响,发现平板形障碍物导致最高超压,圆柱体导致最低超压;Ciccarelli等[5]研究了不同障碍物间距对火焰传播速度的影响,障碍物间距与阻塞率比率为5时,火焰速度最大。

随着流体力学和数值方法的发展和结合,数值模拟方法开始应用于受限空间的气体爆炸现象[6]。Dunn-Rankin等[7]、Naamansen等[8]对障碍物加速火焰传播过程进行二维模拟,发现火焰与障碍物作用后产生明显的超压现象;Gamezo等[9]对置障管道内火焰传播过程进行二维及三维数值模拟,发现初始加速过程是由燃烧产物的扩散引起的,后期阶段的加速则是由激波与火焰相互作用引起的;Sarli等[10]采用大涡模拟与亚格子燃烧模型,很好地预测了障碍物附近非稳态火焰的传播过程;Johansen等[11]对三种阻塞率管道内火焰传播的三维大涡模拟进行研究,得知火焰可以加速到相对于层流火焰速度两个数量级以上的传播速度,而燃烧模式仍处于“层流预混燃烧区”和“薄反应区”。

对于湍流,通常认为大涡模拟比RANS(Reynolds averaged Navier-Stokes)模型计算结果更加精确,但对于高雷诺数流动、壁面湍流流动等问题,需要非常精细的网格,其巨大的计算耗费在工程上应用不切实际[12],RANS方法的易实现性依然是解决湍流燃烧工程中最实际的方法。

本文中,使用工业上最常见的两方程模型(RNG k-ε、可实现k-ε以及SST k-ω湍流模型),对无障碍物管道内丙烷-空气预混火焰传播进行模拟,发现可实现k-ε模型的计算结果与实验数据最接近,认为该模型在模拟湍流燃烧流场中较为精确。鉴于此,采用可实现k-ε模型对不同阻塞率、不同障碍物间距、不同障碍物空间位置对丙烷-空气爆炸火焰加速的效果进行定量分析,重点分析障碍物附近的湍流涡对火焰面作用的详细机理,提供障碍物引发火焰加速现象的理论参考。

1 数值模型及验证

1.1 数学模型

选用基于C方程的Zimont模型对爆燃过程进行模拟,C为一个表征反应进程的标量,表达式为:

(1)

式中:(Sc)t代表Schmidt数;Sc代表反应源项,其表达式为Sc=ρuUtc|,其中,ρu表示预混气体密度,Ut代表湍流火焰速度,湍流燃烧速度采用Zimont等[13]提出的经验表达式:

图9中,障碍物间距大约为1倍管径时,流线最密集,涡量集中区恰好“填满”障碍物之间的体积,因此火焰不会过早的进入涡流区并与之相互作用。阻塞率较高时,涡团的规模较大,湍流强度的进一步增加可能导致湍流涡向反应区输运未燃气速度太快,燃料氧化不充分使得火焰淬火,化学反应速率降低;当障碍物间距逐渐加大,流线随着间距的扩大越来越稀疏,涡量集中区在障碍物间距中所占的比例也越来越小,火焰锋面提前进入涡流区,从而加剧了湍流,使部分火焰来不及燃烧,产生火焰“跳跃”现象;当间距增大到一定程度,涡团的尺度远小于障碍物间距,火焰将有可能接触到管壁,与管壁接触的部分火焰熄灭使它表面积减小,不利于火焰加速;另外,障碍物间距过小,涡团虽“填满”了障碍物之间的空隙,但涡团尺度小,湍流强度低,并没有对火焰起到很好的激励作用。阻塞率较低如0.3时,图5中不同障碍物间距下的火焰加速度相差不大,可认为不同障碍物间距中形成湍流涡尺度均较小,涡流对火焰激励作用较弱,火焰阵面在传播过程中不断与涡流进行相互作用,湍动能增加,增大的火焰面与接触壁面熄灭的火焰面积相抵消。当阻塞率增大到一定程度如0.9时,涡量集中区域在不同障碍物间距内所占比例均较大,并且进入涡团的火焰因氧化不充分快速熄灭,因此火焰与管壁接触的概率很小。由此可见,只有当阻塞率在一定范围内,障碍物间距才影响火焰的加速效果。

不少在小学语文教师在实际教学过程中,经常会发现学生们在学习与生活相关的语文知识时,他们吸收知识时就比其他内容要快。基于这样的情况,小学语文教师在实践教学时,可将生活与语文知识巧妙地结合起来,这样不仅可以进一步加深学生对语文知识的理解,还能激发他们对语文学习的积极性,从而使学生在较短的课堂上最大程度的吸收语文知识,并提高小学语文教学有效性。

(2)

式中:n代表反应产物种类,Yi代表组分i的质量分数,Yi,e代表完全燃烧时组分i的质量分数。 C=0时,表示反应未开始;C=1时,表示完全反应。

(3)

式中:A为模型常数;u'为湍流脉动速度;α为反应物分子热传输系数;lt为湍流特征尺度;Ul为层流燃烧速度。

1.2 实验系统

图3为管道末端压力实验结果与三种湍流模型计算结果的对比。由于管壁边界条件设为绝热,燃烧产生的热量全部用来加热燃烧产物而不能及时扩散出去,因此可实现k-ε模拟结果的压力峰值略大于实验结果,压力达到峰值后下降较为缓慢。而另外两种湍流模型,RNG k-ε和SST k-ω的计算结果比较接近,但均与实验数据相差较大。

图1 实验系统 Fig.1 Schematic illustration of experimental setup

1.3 实验结果与计算结果对比

图6是阻塞率为0.5时3种障碍物空间位置下的火焰发展示意图,时间范围为39~47 ms。不难看出,障碍物位于管道单侧时对火焰的阻挡作用最弱。障碍物位于管道中间时,火焰传播的通道直径仅是障碍物在管道单侧时通道直径的一半,当火焰进入两侧狭小通道,与通道的接触面积增加,由于壁面与火焰的相互作用导致火焰部分熄灭[14-16],因此通道直径越小,火焰中自由基碰到管壁消亡的机会越多;障碍物位于管道两侧时,当通道截面面积突然扩大,混合气体的气流速度迅速减小,增加了反应时间,火焰开始膨胀,接触到障碍物的火焰熄灭,在经历了周期性的收缩流与膨胀流之后,火焰在传播过程中会损失一部分热量。

采用上述物理模型,从距离点火端2.7 m处开始布置障碍物,障碍物的数量为5。模拟方案有2种:(1) 在阻塞率不变的前提下,分别选取间距为0.075、0.15、0.3、0.45、0.6 m(即约为0.5倍管径、1倍管径、2.5倍管径、3.5倍管径、5倍管径)的障碍物,研究间距对火焰传播速度的影响;(2)固定障碍物间距为0.15 m,在阻塞率不变的前提下,分别将障碍物置于管道单侧、管道两侧、管道中间,研究障碍物的空间位置对爆炸火焰的加速效果,图4为阻塞率0.5时3种障碍物空间位置的截面图。

实验装置如图1所示,由管道系统、配气系统、测试系统(包括压力传感器、火焰传感器)、辅助设施组成,管道壁厚5 mm,外径135 mm,管道由316型无缝不锈钢管构成,总长12 m,内壁进行抛光处理,设计最大承受压力6 MPa。为安全起见,在封闭管道末端有自动泄压安全阀。为确保管道密封性,在管道法兰连接处采用聚四氟垫片进行密封。丙烷-空气混合气配置浓度为4.5%,管道初始温度约25 ℃,点火装置采用EPT-6点火能量试验台,点火能量可调,最小点火能量50 mJ,最大点火能量1 J,实验点火能量取1 J。

2.政策的适时调整。适合的政策是经济社会发展的驱动力。但是,多数政策具有时效性的特点,如果超越一定的时空范围就不能发挥作用,就应当调整和变革。抗日战争时期,国难当前,国共两党建立了抗日民族统一战线,过去的许多政策不适应新的形势。毛泽东指出,“过去十年土地革命时期的许多政策,现在不应当再简单地引用”[1]762,尤其是在土地革命后期的许多过“左”的政策,如在经济上消灭资产阶级(过“左”的劳动政策)和富农(分坏田),打击知识分子等政策,“不但在今天抗日时期,一概不能采用,就是在过去也是错误的”[1]762。

图2 火焰速度的实验值与计算值的对比 Fig.2 Comparison of experimental and calculated flame velocities

图3 爆炸压力的实验值与计算值的对比 Fig.3 Comparison of experimental and calculated pressures

3种模型所得的计算结果与实验结果详细对比,见表1。由此可见,可实现k-ε模型计算的爆燃压力、火焰传播速度与实验结果吻合较好,因此选用可实现k-ε模型。

表1 湍流模型的计算结果与实验结果的对比 Table 1 Comparison between calculated results of turbulence models and experimental results

方法pm/MPaδpm/%vm/(m·s-1)δvm/%实验0.508-194.8-可实现k-ε0.565 011.22182.756.59RNG k-ε0.412 723.09145.0034.3SST k-ω0.363 739.68150.4029.52

1.4 数值模拟方案

分别选用RNG k-ε、可实现k-ε以及SST k-ω湍流模型与湍流火焰封闭燃烧模型对无障碍物管道中丙烷-空气爆炸过程进行模拟,图2给出火焰速度实验结果与数值模拟结果的对比。实验结果的最大火焰速度为195 m/s,可实现k-ε计算结果与实验测试结果最接近,且管道轴线上火焰传播速度变化趋势与实验结果基本一致,与实验结果误差的绝对值约为7%。

图4 障碍物空间位置-管道截面图 Fig.4 Schematic diagram of cross section of pipeline and obstacles spatial positions

2 模拟结果与分析

2.1 阻塞率与障碍物间距对火焰传播特性的影响

图5 不同障碍物间距下的火焰平均加速度 Fig.5 Average flame acceleration under different obstacles spacing

图5给出7种阻塞率下不同障碍物间距对火焰平均加速度的影响,纵坐标的火焰平均加速度是指火焰速度为零至达到峰值这个过程中的平均加速度。显然,火焰平均加速度随着阻塞率的增加呈先增大后减小的趋势。阻塞率小于0.7时,火焰平均加速度随阻塞率的增大而增大,障碍物对火焰的激励作用显著;阻塞率大于0.7时,火焰平均加速度急剧下降,可认为阻塞率超过一定值后由于障碍物对火焰的阻挡作用较强,火焰难以穿越障碍物,甚至发生淬火,火焰从淬火到复燃所需时间过长,加速度减小。另外,不同阻塞率下的火焰平均加速度随着障碍物间距变化的程度各不相同:障碍物阻塞率较小(0.3)时,各间距下的火焰平均加速度均在5 000~6 000 m2/s之间,障碍物间距对火焰加速影响不大;而对于较高的阻塞率(0.5~0.7),间距对火焰加速的影响比较明显,其中阻塞率为0.6时,障碍物间距为5倍管径下的火焰平均加速度为9 800 m2/s,间距为一倍管径时的平均加速度达到约13 000 m2/s左右,差别最大;但当阻塞率大于0.7时,间距对火焰平均加速度的影响程度有明显减小的趋势。由图5中还可以观察到,当障碍物间距为一倍管径时所有阻塞率下的火焰平均加速度最大,意味着火焰加速效果最好。

不得不说,企业在经营管理发展的过程当中,积极地做好对全体员工的思想政治教育工作具有重要意义。但是,当前很多企业却并没有重视起这一问题,在思想政治教育的过程当中,所采用的方式方法非常单一、陈旧,大多采用的是讲授和灌输的方式,将其作为一项政治任务,并没有意识到思想政治教育的重要性,态度不积极。由于思想政治教育内容本身就是存在着一定的枯燥性,如果企业仍采用陈旧呆板的教育方式的话,最终所取得的效果也必定不尽人意,无法形成良好的企业凝聚力,这对于企业接下来的发展来说是极为不利的。

2.2 障碍物空间位置对火焰加速效果的影响

图6 阻塞率为0.5时障碍物空间位置对火焰传播的影响 Fig.6 Influence of obstacles spatial position on flame propagation when blocking rate is 0.5

建立与实验系统对应的物理模型,模型壁面为无滑移、绝热边界条件,丙烷-空气在管道中混合均匀,且处于常温、常压下的静止状态。压力速度耦合采用PISO算法,对流项离散化采用二阶迎风格式。

图7~8是阻塞率分别为0.3、0.5、0.7时3种障碍物空间位置下的火焰传播速度及火焰到达各测点的时间。由图7可以看出,障碍物位于管道单侧时火焰传播速度最大,且明显高于其他两种空间位置。而图8中,障碍物位于管道单侧时火焰传播至各测点的时间最短,特别是阻塞率为0.5时,火焰传播至管道末端用时约95 ms,是其他两种空间位置下火焰传播用时的一半。

2)过厚或过薄的空气层厚度会影响辐射板表面温度分布的均匀性,空气层厚度每增加2 mm,表面平均温度平均升高0.5 ℃,供冷能力平均减少7.18 W/m2.当供水温度为14 ℃时,空气层厚度为12 mm的辐射板的均匀性最好.

图7 管道轴线上火焰速度分布图 Fig.7 Distribution of flame velocity on axis of pipeline

图8 火焰到达各测点的时间 Fig.8 Times at which flame reaches monitor points

2.3 障碍物对管道内流场的影响分析

提高火焰传播速度最普遍的方法就是在空间中放置障碍物,火焰推动未燃气体流动时遇到障碍物阻挡使气流方向发生改变[17]。最终形成以障碍物为边界有规则的流体质点,即湍流涡结构,涡量集中的区域具有高强度湍动能[18];随着火焰传播以及压缩波的推动,高温气团使流场中气流黏性系数提高,涡团规模将持续扩大。涡团的形成使火焰锋面两侧形成速度方向相反的剪切层[17,19],涡团中的未燃气体与火焰锋面附近的未燃气体被分隔开,因此火焰传播时不会立即穿透涡团,使得火焰锋面上下两侧形成较大尺度的速度梯度和密度梯度,从而促使火焰在翻越障碍物时,Kelvin-Helmholtz和Rayleigh-Taylor不稳定条件被触发,导致火焰锋面失稳,产生新的涡流,进一步影响火焰传播。图9是阻塞率分别为0.3、0.5、0.7、0.9时不同障碍物间距下火焰翻越第1个障碍物时的流场。结果显示,阻塞率越大湍流涡尺度越大,涡团的扰动大大加剧流场中的湍流,火焰更加不稳定,并且由于障碍物的阻挡导致局部压力快速上升,这时Richtmyer-Meshkov不稳定性也会影响火焰结构,压力波与火焰的相互作用,使湍流动能进一步加剧,这就是前驱冲击波对未燃气体压缩和加热作用的正反馈机制,因此阻塞率的增大在一定程度上加速了火焰传播。但阻塞率大于一定程度时,火焰的传播也会受到障碍物的严重阻碍,穿过障碍物的少量热燃烧产物没有得到能量的补充,无法及时点燃未燃气体,可能导致燃烧过程中的局部淬火,有效能量释放率大大降低。因此存在一个合适的障碍物尺度,使火焰加速机制最强,能量释放率达到最高。在本研究工况下,阻塞率为0.7时,正反馈效应最大。

C方程滤波后表达式为:

图9 障碍物不同间距下火焰翻越障碍物时的流场 Fig.9 Flow field when flame crosses obstacles at different obstacle spacings

以障碍物间距为1倍管径为例,分析空间位置对火焰加速的影响。图10是阻塞率为0.5时障碍物不同空间位置下火焰翻越障碍物的流场。障碍物位于管道中间时,火焰在障碍物两侧形成对称的火焰面,火焰传播的通道直径小,火焰中自由基碰到管壁消亡的机会也越多;同时在障碍物后方形成两个旋转方向相反的反向涡对,由于初始涡对尺度小、强度低,火焰在与涡团相互作用中没有发生明显的不稳定现象,因此火焰传播速度较低。当障碍物位于管道两侧时,由前面分析可知,火焰在经历周期性的收缩流与膨胀流之后,会损失一部分热量;并且每侧障碍物空隙中形成对称涡旋,与障碍物位于管道中间时的情况类似。而障碍物位于管道单侧时,湍流涡规模最大,火焰在短时间内不会进入涡团内部,随着化学反应进行,火焰锋面前后的速度梯度和燃烧释放形成的压力波促进涡团规模持续扩大,以致Kelvin-Helmholtz不稳定和Rayleigh-Taylor不稳定的耦合作用较强,最终使火焰产生较大尺度的褶皱;失稳火焰与未燃气体的快速反应又产生新的涡流,进一步提高了预混气团的输运速率,热量从高温区快速转移到低温区,火焰不断加速。

图10 火焰翻越不同空间位置障碍物时的流场 Fig.10 Vector velocity field when flame crosses over obstacles at different spatial positions

3

(1)障碍物对火焰传播的激励作用显著,但随着阻塞率的增加,激励作用呈现先增大后减小的趋势。阻塞率为0.7时,火焰加速效果最佳,能量释放率最高;阻塞率大于0.7时,障碍物对火焰的阻挡作用较强,导致燃烧过程中局部淬火,有效能量释放率大大降低,不利于火焰加速。

环境、气氛、心理状态、身体状态等也影响与老人交流的质量。室温是否合适、环境是否安静和私密、交流气氛如何、座椅是否舒适、谈话的角度、交谈的角度、内容时间是否充裕等,要优先注意。

增强地方在支出层面的自主权。在下一步事权与支出责任划分中,要赋予地方更多的自主权,如支出标准,在中央总体决策的框架内,增加地方在支出方面自主决策的空间,增加地方自求平衡的空间,解决基层财政困难的问题;同时,要切实打破法定挂钩导致的固化和其他一些限制安排如预算稳定调节基金5%的限制等,体现因地制宜的因素,提升地方支出安排的自主权。

(2)阻塞率在0.5~0.7之间,障碍物间距才影响火焰的加速效果。阻塞率较低时,障碍物间距对火焰加速影响较小;阻塞率较高时,障碍物间距对火焰加速影响较大,若涡团的尺度小于障碍物间距,火焰因氧化不充分和与管壁接触而部分熄灭,不利于燃烧反应进行,当障碍物间距为一倍管径,涡团恰好“填满”障碍物的空隙,火焰不会过早进入涡流区,此时火焰加速最快;当阻塞率增大到一定程度,涡团在不同障碍物间距内所占比例均较大,火焰与管壁接触的概率很小。阻塞率为0.7时,障碍物间距为一倍管径下的火焰加速度峰值为14 000 m/s2

(3)障碍物的空间位置对火焰加速的影响非常显著。阻塞率相同时,当障碍物位于管道单侧,火焰失稳现象明显,火焰传播速度最高;当障碍物位于管道中间和管道两侧,对火焰的激励作用较弱,火焰面没有明显褶皱,火焰传播较慢。

要多渠道吸引优秀人才投身于贫困地区职业农民教育事业;积极加强“双师型”培训师资队伍建设,选择兼具农业专业知识与农业生产技术技能的“土专家”“土秀才”,成立“讲师团”,对贫困农民能够进行田间地头式的指导,面对面讲授诸如作物病虫害防治、田间管理、科学施肥等农业知识,传授农业生产技术,答疑解惑,交流经验。

参考文献

[1] CHAPMAN W R, WHEELER R V. VI.-The propagation of flame in mixtures of methane and air. Part V. The movement of the medium in which the flame travels[J]. Journal of the Chemical Society (Resumed), 1927:38-46.

[2] DOBASHI R. Experimental study on gas explosion behavior in enclosure[J]. Journal of Loss Prevention in the Progress Industries, 1997,10(2):83-89.

[3] MASRI A R, IBRAHIM S S, NEHZAT N, et al. Experimental study of premixed flame propagation over various solid obstructions[J]. Experimental Thermal and Fluid Science, 2000,21(1):109-116.

[4] IBRAHIM S S, MASRI A R. The effects of obstructions on overpressure resulting from premixed flame deflagration[J]. Journal of Loss Prevention in the Process Industries, 2001,14(3):213-221.

[5] CICCARELLI G, FOWLER C J, BARDON M. Effect of obstacle size and spacing on the initial stage of flame acceleration in a rough tube[J]. Shock Waves, 2005,14(3):161-166.

[6] CATLIN C A, FAIRWEATHER M, IBRAIIM S S. Predictions of turbulent, premixed flame propagation in explosion tubes[J]. Combustion and Flame, 1995,102(1):115-128.

[7] DUNN-RANKIN D, MCCANN M A. Overpressures from nondetonating, baffle-accelerated turbulent flames in tubes[J]. Combustion and Flame, 2000,120(4):504-514.

[8] NAAMANSEN P, BARALDI D, HJERTAGER B H, et al. Solution adaptive CFD simulation of premixed flame propagation over various solid obstructions[J]. Journal of Loss Prevention in the Process Industries, 2002,15(3):189-197.

[9] GAMEZO V N, OGAWA T, ORAN E S. Numerical simulations of flame propagation and DDT in obstructed channels filled with hydrogen-air mixture[J]. Proceedings of the Combustion Institute, 2007,31(2):2463-2471.

[10] SARLI V D, BENEDETTO A D, RUSSO G. Sub-grid scale combustion models for large eddy simulation of unsteady premixed flame propagation around obstacles[J]. Journal of Hazardous Materials, 2010,180(1):71-78.

[11] JOHANSEN C, CICCARELLI G. Modeling the initial flame acceleration in an obstructed channel using large eddy simulation[J]. Journal of Loss Prevention in the Process Industries, 2013,26(4):571-585.

[12] 孙明波,汪洪波,梁剑寒.复杂湍流流动的混合RANS/LES方法研究[J].航空计算技术,2011,41(1):24-28.

SUN Mingbo, WANG Hongbo, LIANG Jianhan. Evaluation of hybrid RANS/LES method for complex turbulent flow simulations[J]. Aeronautical Computing Technique, 2011,41(1):24-28.

[13] ZIMONT V L, BATTAGLIA V. Joint RANS/LES approach to premixed flame modeling in the context of the TFC combustion model[J]. Flow, Turbulence and Combustion, 2006,77(1):305-331.

[14] 陈鹏,李艳超,黄福军,等.方孔障碍物对瓦斯火焰传播影响的实验与大涡模拟[J].爆炸与冲击,2017,37(1):21-26.

CHEN Peng, LI Yanchao, HUANG Fujun, et al. LES approach to premixed methane/air flame propagation in the closed duct with a square-hole obstacle[J]. Explosion and Shock Waves, 2017,37(1):21-26.

[15] SHEN Xiaobo, HE Xuechao, SUN Jinhua. A comparative study on premixed hydrogen-air and propane-air flame propagations with tulip distortion in a closed duct[J]. Fuel, 2015,161:248-253.

[16] LV Xianshu, ZHENG Ligang, ZHANG Yugui. Combined effects of obstacle position and equivalence ratio on overpressure of premixed hydrogen-air explosion[J]. International Journal of Hydrogen Energy, 2016,41(39):17740-17749.

[17] 王公忠,张建华,李登科,等.障碍物对预混火焰特性影响的大涡数值模拟[J].爆炸与冲击,2017,37(1):68-76.

WANG Gongzhong, ZHANG Jianhua, LI Dengke, et al. Large eddy simulation of impacted obstacle’s effects on premixed flame’s characteristics[J]. Explosion and Shock Waves, 2017,37(1):68-76.

[18] 张欣,王涛,侯效森,等.低热值气体发动机涡团与火焰面相互作用[J].内燃机学报,2016,34(6):537-542.

ZHANG Xin, WANG Tao, HOU Xiaosen, et al. Flame-vortex interaction of lower heating value gases fueled engine[J]. Transaction of CSICE, 2016,34(6):537-542.

[19] JOHANSEN C, CICCARELLI G. Modeling the initial flame acceleration in an obstructed channel using large eddy simulation[J]. Journal of Loss Prevention in the Process Industries, 2013,26(4):571-585.

周宁,王文秀,张国文,宗永迪,赵会军,袁雄军
《爆炸与冲击》 2018年第5期
《爆炸与冲击》2018年第5期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号