资源目录
压缩包内文档预览：
编号：66729908
类型：共享资源
大小：81.74MB
格式：ZIP
上传时间：2020-04-04
上传人：机****料
认证信息
个人认证
高**（实名认证）
河南
IP属地：河南
50
积分
- 关 键 词：
-
温室
卷帘
设计
- 资源描述：
-
压缩包已打包上传。下载文件后为完整一套设计。【清晰，无水印，可编辑】dwg后缀为cad图纸，doc后缀为word格式，所见即所得。有疑问可以咨询QQ 197216396 或 11970985
- 内容简介：
-
温室窗帘和窗户通风设备的物理模型 我们采用实践与理论相结合的方法研究了温室的通风设备，主要是研究温室卷帘与开窗设备在低气压下的使用寿命及其在不同风速和气压下的波动情况。这些波动与平均风速有关，而平均风速是通过能量谱分析得来。由此推知，气压与平均风速有关，而气体紊流对平均风速的估算有很大影响。温室通风设备与启动电位之间的关系是基于对流体机械的突破性研究。关于气流与风压和气温的变化相对应这个预测符合实验得出的数据，大体上看，他们之间差异小于 20%.符号标志符号标志 A 面积 z 高度，mCc 对流影响系数 z0 表面粗糙度，mCkl kolmogorov 常量（0.5） Cwf 摩擦系数 希腊字母符号希腊字母符号Cu 紊乱动能常量 开窗下垂面与框架间的角度，度F(n) 能量谱密度，m2/s 热能膨胀系数，K-1g 重力加速度，m/s*s 孔积率 h 高度，m 冯.卡曼常量（0.4）H 溶度参量，m 动力黏度，NsKp 渗透率，m2 密度，千克/立方米L 长度参量，m 气压系数 n 频率，HzP 大气压，PaP0 绝对气压，PaPst 相对气压，Pa 下标符下标符Pw 风压，Pa fr 循环运动Q 流速，立方米、秒 I 内侧面ru 动能消耗比率，平方米、秒 l 超长度T 绝对温度，K rmsw 平均风速的平方根u 流体速度，立方米/平方毫秒 s 次长度u* 摩擦速度，m/s st 静态V 容积，立方米 w 风Y 惯性系数1.绪论绪论温室的自然通风设备是一个极其复杂的过程，这一过程取决于温室的所有参数（额定值，位置和窗户的几何形状以及漏孔面积等等）和外部环境条件。这种现象持续不断的利害关系会导致自然通风设备在质量和能量上的不平衡，并进一步严重影响室内环境。自从 1954 年 Morris 和 Neale1 第一次对温室开窗设备中的自然通风设备进行实验性研究以来，人们越来越多地考虑到潜在的物理现象对通风设备影响。几十年后，Baiely2和 Cotton3以及 Miguel 以及其他 4 的实验论证了在植物和温室顶棚之间使用遮阳网的潜在好处。近期以来，在开窗设备5处采用幕帘下是为了防止昆虫进入，这样即可减少用化学杀虫剂。所以，幕帘表面的气流研究同样是一门重要的课题。现有的对这个课题的分析方法大致可以分为以下两种：（1）：经验主义和半经验主义的研究方法研究穿过6.13和幕帘2.5.14的气流。（2）：以流体机械公式为指导数字化的解决问题。15.17 属于第一种类型的所有模式基于的理论思想是在气流与驱动气流循环运动的潜在物之间有一种简单的非线性关系。这种分析认为温室室内的空气是不可压缩的，这一分析纯碎根据经验或贝努利方程所得。尽管简单的假设是与具体实验下得到的改正系数结合在一起的，还是大大限制了这项研究成果的有效范围。 第二种模式中，采用一种数字化模拟程序来解决惯性方程，动能和能量方程，决定速度，温度及压力和温室内最终气流循环模式。 本文介绍的这项研究有以下两个方面：（1）加深了解气流通过带孔卷帘（保温，遮光，防虫）和天窗的交换过程中的物理特性。更重要的是论述带孔卷帘和天窗使用寿命的影响因素以及启动差位的各项参数。（2）推出一种简单而精算的计算风与温室框架对卷帘和天窗的影响的公式。2.理论理论 无约束的和约束的流体运动经过有孔设备和天窗的物理现象表明这种现象与天窗或带孔设备的特征参数和驱动循环运动势差的特征参数有关。公式精确描述这个现象必须考虑到这些因素。附录 A 中附有该公式。根据附录所示，穿过任何一天窗或带孔卷帘的气流与驱动势差有关，如下：（/）Q/ t + Kp-1Q +YA-1Kp-1/2QQ+0.5(AHCc2)-1QQ= -A(pw/H+pst/H) (1)和 Q=A其中，Q 是流量， 是气体密度，A 是工作面积， 是气体动力黏度， 是材料的孔积率（单位体积容纳流体的体积） ，Kp 是材料的渗透度（流体通过媒介体的能力） ，Cc表示对流影响系数，H 是渗透参数，pw 是风压差，pst 是与温度相关的压力差，Y 在附录 A 中有定义。为了便于计算气流通过多层次密封箱或天窗的空间时的运动，必须建立计算机联网服务。薄膜温室是一个分为两层的密封区域（如图 1） ，图 2 说明了气流的联网计算原理。 自然通风设备物理模型自然通风设备物理模型 在只有一个天窗和卷帘的温室内，网络对循环运动的控制由三个接点组成，这三个接点由两个电阻连接在一起（图 2 a） 。当卷帘有空隙时，网络控制也是由三个接点组成，但这三个接点由三个电阻串接在一起。 2.1 天窗与带孔卷帘设备的气流特征参数天窗与带孔卷帘设备的气流特征参数 可渗透材料的气流特征参分为孔积率和渗透度。18. 19 孔积率表示单位材料容纳气体的体积与总材料容纳气体体积之比，大小在 0 至 1 之间（01） 。对一个孔隙来说，孔积率规定为 1（=1） ，因为整个孔隙都充满气体。 任何一种可渗透材料都能透过气体，这种能力的度量衡称之为渗透度且符合气体运动原理渗透度不仅与流体运动黏度有关还与粒子扩散时与阻碍物的碰撞频率有关。18.20 就孔隙材料而言，气体分子的碰撞频率大于 103Hz，因为气体运动黏度的数值为 10-5，渗透度的数值即小于 10-7.而对于开口设备（窗，门）来说，碰撞频率接近于零所以渗透度 Kp- 。与以上一致，连接两个接端（毛细孔与大开口）的数量即为渗透度 Kp。对于空隙材料来说气体可认为是不可压缩的（1/Cc20） 。由公式（1）推出：（/）Q/ t + Kp-1Q +YA-1Kp-1/2QQ=-A(pw/H+pst/H) （2）对于开口（=1）材料的渗透度 Kp- ,公式（1）的左边第二和第三项可忽略，导出：Q/ t+0.5(AHCc2)-1QQ=-A(pw/H+pst/H) （3）公式（2）即为著名的 Forchheimer 公式，适用于多孔隙材料。对于微弱气体可舍去二次项，得出 Darcy 定律。当气流静止时公示(3)简化为贝努利方程。2.2 驱动势差驱动势差 引起气流的驱动势差是因为温差（静因素）或风速（导致风压变化）或两者同时作用时产生的。任何一种原因都可以产生一种梯度势差引起稳定循环。风速的波动会产生一种额外的循环运动。通过温室天窗的气流波动可以分为震动循环运动和旋窝型穿透运动。震动循环运动是由于风的波动和室内气体的可压缩性。旋窝型穿透运动是由于热气流的涡流对密封室内的气流造成转动的影响。2.2.1 温室对气流运动的影响温室对气流运动的影响 当多孔薄膜或天窗的内外有温差时，就会产生一种静态压力，这种压力可以导致气流运动。设想两边温度不同，根据以下公式，两边压力也不同：st=gh （4）和 = T其中， T 表示多次测量的绝对温度差，g 表示重力加速度，h 表示垂直高度差， 表示热膨胀系数。2.2.2 风速对气流运动的影响（风压）风速对气流运动的影响（风压）测量时风速是一个不定量。在 t 时刻的瞬时值可以人为是组成部分的平均值和波动值的总和，那就是 w和uw。一段时间间隔内的平均风压与风速有关（附录 B）： Pw=0.5( 2w +uwuw) (5)其中， w表示平均风速，uw表示风速的波动值。有一些风速计可以直接读出平均风速和平均风速的平方根。假设uw有一个 Gaussian 概率分布，可以简便用下面的公式代替公式（5）：Pw=0.5 2w+-1urmsw-1 (6)其中 urmsw 表示平均风速的平方根。普通情况下的风速计仅仅提供平均风速的读数，因此，不仅能显示平均风速，而且不忽略风的波动产生的影响的风速测量计是非常有用的。因此要将动态风速的平均值与静态风速的平均值联系在一起。动能消耗比率() 根据 Kolmogorov 定律21： r2/3=0.75ckl-1F(n)( w/2)2/3n5/3 (7)其中，F（n）为能量谱密度，n 表示频率，ckl为 kolmogorov 常量（0.5） 。再利用离表面 z 距离的动能湍流与动能消耗比率21之间的关系0.5 uwuw=c-0.5rk(z+z0)2/3 (8)风速波动的平方根值和平均值的关系如下：uwuw=3 F(n)( w/2)2/3n5/3 c-0.5rk(z+z0)2/3 （9）其中，k 是冯.卡曼常量（0.4） ，c是常量（0.99） ，z 为测量风速的高度，而z0为表面粗糙度（z0的值在参考文献 22 处可寻） 。为了获得 ，弄清楚 F(n)和 n 是必要的，而且需要实施一种能量谱分析的方法。这种能量谱分析法23是用来测量频率连续不断变化的摆动过程中变量方差的。这种分析不仅可以确定 值，而且可以获得一些阐明和刻划湍流波动结构的额外信息，并且可以确定风域内主要涡流的频率。风速通常在一个参考高度测量。应该导出一个面积系数，该系数与相对基本风域的参考高度的风域有关。这个系数是由大气风界层的纵断面上的风速决定的并遵循风向定位封闭室（温室） 。这些数值都由风道分界层的多次测量所得。3.实验研究实验研究为了证实波动存在的潜能，同时也检验模型的可适用性，我们进行了如下实验。实验是在两个向东（E-W）方位建造的温室中进行的。两个温室尺寸相同（如图 3）：屋檐高 4.5米，屋顶坡角 22 。 ，宽 4.1 米，长 6.6 米。每个温室内在离地面 2.9 米处的水平面装有一个LS11 型号的隔热板（Miguel 以及其他等等18验证渗透度 7*10-10m2，孔积率 0.99） 。在实验的第一部分，每个温室内的薄膜都是合格的。实验的第二部分，在卷帘薄膜的中间打开了 0.20 米*3.80 米细小的口子，看上去像一条水平裂缝。 图 3.两个尺寸相同的玻璃温室图解每个温室内部都被牢牢密封住，墙上装有绝热板条，地面盖有聚苯乙烯泡沫绝热层，除了屋顶的窗户，每两个窗户间装有 2.05 米*0.90 米大小的活动板条。开窗斜度可升高 30度。装有铝制表面的水平圆柱形电加热器（8 个 3.80 米，直径为 0.05 米的柱面，成对排列，每两个间隔 0.35 米，每队间隔 1.15 米） ，如图 4 所示。在实验中，隔热板上下方的温度会不同，加热器被用来抑制这种温度的不断变化。图 4. 加热系统排列的俯视图在每个温室中安装 25 个铜制镍铜合金温差电偶，用来测量隔热板上下及户外的温度。它们被均匀分布：10 温差电偶被分布在隔热板下方，10 个分布在上方，还有 5 个曝露在室外。为了保证获得一个快速响应8（响应频率大约为 12Hz） ，温差电偶都是用非常细的电线（直径约为 2.5*10-5米）做成。用带有薄膜层的压力转换器来测量压力，室内测量板的上部和下部（每个部分测量三次） 。室外瞬时风压的测量位置在天窗上 0.20 米处。风速的测量用一个快速响应的风速计（响应频率为 9.5Hz）,把它安置在距天窗 0.20 米处，同时也可测方向。气流的测定用到一种自动扫描气体的仪器。在实验中，用到了恒定流速和衰减率9.25.穿过隔板狭缝的气流量用总波动量减去通过隔板的波动量。每个温室中，在地面上用两个小鼓风机和打满孔的桶把扫描气体（N2O）分离出来。在实验中，每个温室内部空气取样检验都要在 18 个不同的位置获得（9 个在隔板上层，9个在隔板下层） ，并在红外线下进行分析。 两个取样检验在 1996 年 2 月至 3 月之间的 32 天进行：（1） 气流会穿过温室天窗和卷帘是由风速产生的。实验必须在室内室外温差少于2.00.5oC，风速高于 1.5 米每秒下进行。（2） 气流会穿过多孔卷帘和卷帘上的矩形狭缝（0.02 米*3.80 米）仅仅是由于卷帘上部和下部的温差产生的（稳定条件下） 。为了减小风压的影响，测定时要关闭所有窗户（背风面的窗户可开 2o） 。实验的第一部分（气流的产生是风的作用） ，在风速和室内外温差存在下收集的数据频率 10 分钟内都为 8HZ 。实验的第二部分（气流的产生是缓慢聚集起来的） ，在卷帘上部和下部温差稳定（T1 oC）下采集的数据频率为 1.66*10-2HZ（每 60 秒内） ，稳定情况下获得。4.结果和论述结果和论述4.1 风速和风压 谱分析是在风速的采集在 10 分钟内的频率都为 8HZ的条件下进行的。取样频率是通过对可获得的最高频率分析测定的，即取样频率的一半（Nyquist 频率）在本实验中为4HZ。由于风力的特征在背风面和迎风面会有不同，9所以谱分析要分别在两处进行。关闭温室内窗户所得的结果绘制成图，如图 5 和 6 所示。图 5. 三种不同风速下的风速能量谱密度：1.27m/s(*)，3.49m/s(o)，5.50m/s(+)，在迎风面距屋顶 0.20 米处测得（-为-5/3 的斜度）图 5 所示为三种不同风速下的风速能量谱密度：1.27m/s，3.49m/s，5.50m/s，在迎风面距屋顶 0.20 米处测得。图 6. 风速为 0.52m/s(*)，2.24m/s(O)，3.37m/s(+)时，在背风面距屋顶 0.20 米处测得能量谱密度（-倾斜度为-5/3）图 6 所示为背风处测得的相应风速能量谱密度（风速为 0.52m/s，2.24m/s，3.37m/s，在背风面距屋顶 0.20 米处测得） 。如图 5 和图 6 中所描绘的，能量谱均衡的分布，频率/能量为-5/3 的范围，符合Kolmogorov 定律。类似的可获得迎风面和背风面的风速能量谱。两处能量谱的主要波动最高点在频率低于 0.1，0.2HZ 处，频率高于 1HZ就没什么特别作用了，即低频在风域中占主要地位。事实上，风速中主要的能量大的涡流就发生在这个低频范围内。这证实了 Kaimal等26和 Bot.8的研究学说。参数 可以从 图 5 图 6 中能量谱密度和频率值计算出。计算中的表面粗糙度为 0.04米22，结果如表 1.4.2 穿过卷帘和天窗的气流穿过卷帘和天窗的气流对于气流仅仅因为温差所引起（即对于稳定温度条件下 Q/t0）的如图 7 和图 8所描绘的（P w0） 。图 7. 气流穿过卷帘对应于用公式（2）和（4）预算的卷帘温度差值（测量数据（*） ）图 8. 气流穿过卷帘中心的矩形狭缝对应于用公式（3）和（4）预算的卷帘温度差（测量数据（*） ）图 7 所示为实验中气流穿过卷帘对应于用公式（2）和（4）预算的卷帘温度差值。图8 所示为实验中气流穿过卷帘中心的矩形狭缝对应于用公式（3）和（4）预算的卷帘温度差。对于气流仅仅因为风速引起的（0）,如图 911 所示。在这些图中，气流穿过天窗或卷帘与多孔隔板由公式（2）得出的值，天窗由公式（3）得出的值相对应。图 9. 预测气流穿过迎风面（*）和背风面（O）时与迎风面的相对压力（开窗度为4o）图 10. 预测气流穿过迎风面（*）和背风面（O）时与背风面的相对压力（开窗度为4o）图 11. 预测气流穿过隔板与迎风面的相对压力（任意面有一扇窗打开，开窗度为20o）迎风面（*）,背风面（O）图 12. 预测气流穿过隔板中心狭缝与迎风面的相对压力（任意面有一扇窗打开，开窗度为20o） 迎风面（*）,背风面（O）图 9 所示为实验中预测气流穿过迎风面（*）和背风面（O）时与迎风面的相对压力（开窗度为 4o）图 10 所示为实验中 预测气流穿过迎风面（*）和背风面（O）时与背风面的相对压力（开窗度为 4o）图 11 所示为实验中预测气流穿过隔板与迎风面的相对压力。其中，温室中迎风或背风面的一扇窗户要打开（开窗度为 20o） 。图 12 所示为实验中预测气流穿过隔板中心狭缝与迎风面的相对压力。其中，温室中迎风或背风面的一扇窗户要打开（开窗度为 20o） 。图 7 和图 11 描绘了由于存在温差，各自的风域推动气流穿过有孔的隔板，它们显示出了两种不同的波动规律：Darcy 波动规律（如图 7）和 Forchheimer 波动规律（如图 11） 。这些发现归功于 Bear,Bachmat19和 Bailey2。对于 Reynolds 数据（Re=uKp1/2/）小于整体值时，气流与推动势成比例，即温室与周围空气间正常温差范围（T25K）内产生的情形。对于 Reynolds 数据大于整体值时，必须加上气流的二次项（Forchheimer 波动规律） ，即波动是由风速高于 0.25m/s 时产生的。 实验气流值与模型预测等效的气流值相比较，所得的误差总是小于 20%，出了图 12 之外。图中显示的大部分分散值是由于气流穿过裂缝所致，这些气流是从总的波动量中减去穿过隔板的波动量获得。由于差别很小，所以误差相应就很大。5.结论结论1. 对于换气设备起重要作用的、能量大的涡流是由频率低于 0.10.2HZ范围内湍急的风速下引起的。2. 据估计，这种湍急的风速在平均风速中占 13%至 52%，也就是说在总风压中占了重要地位。3. Forchheimer 方程式方程（2）描述了气流穿过有孔的隔板（孔集中在隔板中间） 。可是，对于温室与周围空气间温差在正常范围内所引起的气流波动，气流的二次项可忽略，所以方程就缩写成了 Darcy 法则（Kp-1Q=-Ap/HP） 。4. 通过理论预测的穿过隔板和通风口的气流与实验值很接近，大体上它们之间的差别小于20%。参考文献参考文献1 Morris L G; Neal F E the infra-red carbon dioxide gas analyzer and its use in greenhouse research. National Institute of Agricultural Engineering Report，Silsoe，UnitedKingdom，19542 Bailey B J Glasshouse thermal screens: air flow through permeable materials. Departmental Note no. DN/G/859/04013. National Institute of Agricultural Engineering, Silsoe, United Kingdom, 19783 Bailey B J; Cotton R F Greenhouse thermal screens: influence of single and double screens on heat loss and crop environment. Departmental Note no. DN/G/982/04013. National Institute of Agricultural Engineering, Silsoe, United Kingdom, 19804 Miguel A F; Silva A M; Rosa R Solar irradiation inside a single span greenhouse with shading screens . Journal of Agricultural Engineering Research, 1994, 59, 61-725 Kosmas S R; Riskowski G L; Christianson L Force and static pressure resulting from airflow through screens. Transactions of ASAE, 1993, 36, 1467-14726 Okada M; Takakura T Guide and data for greenhouse air conditioning. 3: heat loss due to air infiltration of heated greenhouses. Journal of Agricultural Meteorology (Tokyo), 1973, 28,223-2307 Kozai T; Sass S A simulation natural ventilation for a multi-span greenhouse. Acta Horticulturae, 1978, 87, 39-498 Bot G P A Greenhouse climate: from physical processes to a dynamic model. PhD dissertation, Agricultural University of Wageningen, The Netherlands, 19839 De Jong T Natural ventilation of large multi-span greenhouses. PhD dissertation, Agricultural University of Wageningen, The Netherlands, 199010 Montero J I; Anton A; Biel C Natural ventilation in polyethylene greenhouses with and without shading screens. Proceedings of the International Seminar and British-Israel Workshop on Greenhouse Technology, Bet Dagon, Israel, 1990, 65-7111 Fernandez J E; Bailey B J Measurements and prediction of greenhouse ventilation rates .Agricultural and Forest Meteorology, 1992, 58, 229-24512 vans Ooster A Using natural ventilation theory and dynamic heat balance modeling for real time prediction of ventilation rates in naturally ventilated livestock houses. Ag. Eng. 94 Milano, Report N. 94-C-012, 1994 13 Bullard T; Menses J F; Merrier M; Papadakos G The mechanisms involved in the natural ventilation of greenhouses. Agricultural and Forest Meteorology, 1996, 79, 61-7714 Batemans L Assessment of criteria for energetic electiveness of greenhouse screens. PhD dissertation. University of Ghent, Belgium, 198915 Okushima L; Sass S; Nara M A support system for natural ventilation design of computational aerodynamics. Alta Horticultural, 1989, 248, 129-134 16 Outworker E N J; Voskamp J P; Alaskan Y Climate simulation and validation for an aviary system for laying hens. Ag. Eng. 94 Milano, Report N.94-C-063,199417 Mastitis A; Bot G P A; Pacino P; Scarascia-Mugnozza G Analysis of the decency of greenhouse ventilation using computational fluid dynamics. Agricultural and Forest Meteorology, 1997, 85, 217-22818 Miguel A F; Van de Barak N J; Silva A; Bot G P A Analysis of the airflow characteristics of greenhouse screening materials. Journal of Agriculture Engineering Research, 1997, 67, 105-11219 Bear J; Bathmat Y Theory and Applications of Transport Phenomena in Porous Media. Dordrecht: Kluwer Academic20 LeBron G; Clout A thermo dynamical modeling of fluid flows through porous media: application to natural convection. International Journal of Heat and Mass Transfer, 1986, 29, 381-39021 Tenneco H; Lumley J L A First Course in Turbulence. Cambridge, MA: MIT Press, 197222 Euro code ENV 1991-2-4. Basis of design and action on structures: wind action. NNT Delft, 199523 Turkey J W The sampling theory of power spectrum estimates. Symposium on Application of Autocorrelation Analysis to Physical Problems. Washington, Office Naval Research 1950, 47-6724 Clarke R M Observational studies in the atmospheric boundary layer. Quarterly Journal Royal Meteorology Society, 1970, 96, 91-11425 Sherman M H Tracer-gas techniques for measuring ventilation in a single zone. Building and Environment, 1990, 2526 Kemal J C; Wynyard J C; Izumi Y; Cote O R Spectral Characteristics of surface-layer turbulence. Quarterly Journal Royal Meteorology Society, 1972, 98, 563-58927 Miguel A F; Van de Barak N J; Silva A; Bot G P A Forced fluid motion through openings and pores. Building and Environment 1998(in press)28 Walker I S; Wilson D J Evaluating models for superposition of wind and stack elect in air infiltration. Building and Environment, 1993, 28, 201-21029 Peixoto J P; Fort A Physics of Climate. New York, USA: American Institute of Physics, 1992附录附录 A: 穿过通风口和气孔的流速运动方程 对于呈线性的气流穿过一个可渗透的材料，其运动方程可描述为18.27(/)u/t+Kp-1u+YKp-1/2uu+(/2)uu/j（/） （2u/j2）=-pt0/j (A1)其中 Y=4.3610-2-2.12总压力28 Pt0=Pw+Pst 公式中，u 是变速度， 是密度，Pw 压力是由风或力学原因造成，Pst 压力是累积作用的产生的（热气压） ，Pt0 是总压力，pt0/j 是在 j 方向
- 温馨提示:
1: 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
人人文库网所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
川公网安备: 51019002004831号