Relationship Analysis of Port Traffic Accidents and Port Environment Factors
-
摘要: 为了评价港口交通环境的安全性, 合理确定各环境要素的权重, 在应用因子分析方法分析全国九个港口环境要素数据的基础上, 对各港口的综合因子得分与港口事故数进行线性回归, 得到回归方程, 通过回归方程中各要素系数确定各环境要素对港口交通事故的影响程度。分析结果表明: 能见度、助航标志、VTS管理程度和分道通航程度对交通安全影响较大, 航道长度和交叉点数影响较小, 基本上反映了各环境要素对港口交通事故的影响程度, 可见因子分析方法可定量确定各环境要素对交通安全的影响, 可用于港口交通事故的原因分析与评价。Abstract: In order to evaluate the safety of port traffic environment and determine the weighting of each environment factor, the environment factors data of nine ports were analyzed by factor analysis method, the relationship of port integrated factor scores and port traffic accident numbers was studied, their regression equation was deduced, each factor effect degree on port traffic accidents were determined by each factor coefficient in regression equation.Analysis result indicates that visibility, navigational aids, VTS management degree and traffic separation degree have greater effect on traffic safety, the length of channel and the number of turning points have less effect on traffic safety, which basically reflects the effect degrees of environment factors on port traffic accidents.It is shown that factor analysis method can quantitatively determine the effect degree of each environment factor on traffic accidents, it can be applied to the relationship analysis and assessment of port traffic accidents.
-
0. 引言
为了修筑高质量的沥青路面, 沥青路面施工材料越来越多地采用新工艺制备的改性沥青[1], 但目前的传统指标难于表征沥青的路用性能, 不能有效区别道路沥青性能的优劣。在“八五”科技攻关项目中开展的道路沥青及沥青混合料路用性能的专题研究中, 推荐的针入度指数、当量软化点、当量脆点等评价指标对于重交沥青有较好的区分, 但对于改性沥青、天然沥青及高模量沥青不能给予合理的评价。与此同时, 美国也在其战略公路研究计划(SHRP)中提出一整套基于性能分级的沥青结合料路用性能规范, 其中采用动态剪切试验, 提出车辙因子指标来表征沥青的高温性能, 但研究证明, 由于加载方式的缺陷, 车辙因子并不适用于评价改性沥青[2-4]。目前在欧洲, 大多数国家以在60 ℃时非常低的频率或剪切速率下获得的粘度作为沥青的性能评价指标[5-7]。本文通过动态频率扫描及长时间单循环蠕变恢复试验, 对沥青高温评价指标进行研究, 分析不同流变模型对沥青拟合的影响, 得出了针对沥青高温性能指标的零剪切粘度适用范围。
1. 试验材料及仪器
试验材料分别选取韩国SK重交沥青、埃索改性沥青和高模量改性沥青。
试验仪器采用的是应力控制的高级流变仪Advanced Rheometer-2000(AR2000), 在控制应力中, 可以任意施加或释放应力, 并能直接测量样品的实际特性。
试验通过动态频率扫描及静态蠕变恢复来测试不同沥青在不同加载模式下的应变响应, 并且采用不同的流变模型对试验结果进行拟合对比。
2. 试验方法及数据分析
2.1 动态频率扫描
零剪切粘度(ZSV)是剪切速率接近于0时的粘度极限值, 可通过动态模式加载试验来获得。本文采用25 mm平行板, 对基质沥青在0.1~10.0 Hz, 埃索改性沥青和高模量沥青在0.1~100.0 Hz的频率范围内对沥青结合料进行动态频率扫描, 温度采用结合路面实际高温的65 ℃, 应变控制统一采用5%, 保证3种沥青均在线粘弹性范围内。试验得到的3种沥青的频率扫描结果见图 1~3。
动粘度η'是复数粘度η*的实轴分量, 其表征因沥青的粘性抵抗造成的能量损失, 即
η'=G″/ω'
式中: G″为损失模量; ω'为加载频率。
图 1中, 基质沥青的粘度随着频率的减小很快达到稳定, 在1 Hz内就表现出牛顿流体特性, 说明在高温条件下, 基质沥青的变形基本上都是不可恢复的粘性部分变形, 对此进行图形外延很容易得到沥青的零剪切粘度。图 2中埃索改性沥青的粘度并没有随着频率的减小而趋于稳定, 并且动粘度η'始终低于复数粘度η*, 表明沥青的弹性成分始终占有一部分的比例; 相位角δ表征沥青的粘弹比例, 其随着频率的减小也有降低的趋势, 表明随着频率的减小, 改性沥青的弹性成分所占比例逐渐增大, 说明在试验的低剪切速率下改性沥青无法进入粘性流动状态, 不能通过试验的图形外延直接得到零剪切粘度。图 3中高模量沥青的粘弹特性与频率关系基本上等同于SBS改性沥青, 但相位角δ与频率的关系表现的更为复杂, 随着频率的减小, 相位角δ呈现先增大后降低的凸性曲线, 在1~10 Hz之间处于稳定阶段, 在两端的δ都逐渐降低, 高模量沥青表现出的这种性质区别于改性沥青。加载频率的高低模拟的是路面行车速度的快慢, 10 Hz加载频率等同于路面行车速度为60 km·h-1, 在高速行驶的条件下, 高模量沥青的弹性成分所占比例增加, 有利于抵抗车辙的形成, 其复数粘度与加载频率的关系等同于埃索改性沥青, 也不能通过试验的图形外延直接得到零剪切粘度。
基质沥青在加载时变形呈线性增长, 卸载后变形基本是不可恢复的, 表明其变形主要是流动变形, 而改性沥青在荷载作用下的非线性粘弹性和材料本身延迟弹性的影响, 与车辙相关的变形不只是粘性流动变形, 还包括未恢复的延迟弹性变形。对基质沥青适用的图形外延方法, 对于改性沥青并不适用。由于在试验的范围内随着频率的降低很难达到粘度稳定的状态, 而这个纯粹的粘性常数一定存在, 因此, 需要采用流变学模型拟合ZSV。比较常用的用于拟合ZSV的模型有Cross和Carreau模型。
Cross模型描述的是伪塑性非牛顿流体的应力应变关系, 即
η-η∞η'0-η∞=11+(kω)m (1)
式中: η为粘度; η'0为零剪切粘度; η∞为剪切速率很大时粘度达到的最小值, 称为第二牛顿流区域粘度; ω为稳定状态的剪切速率; k为常数, 具有时间量纲的材料参数; m为常数, 无量纲材料参数。
对于一般试验条件下(如频率扫描的范围在0.1~10.0 Hz之间)得到的粘度η, 可以假设η'0≫η≫η∞, 于是式(1)可以简化为
η=η'01+(kω)m (2)
SYBILSKI D的研究表明, 对同种沥青来说, k值随试验温度的升高而显著降低, m值受温度的影响较小, 主要体现的是材料的剪切敏感性, 也称为速率指数, 当m趋近于0时, 沥青结合料的流变特性趋近于牛顿流体。
Carreau模型的表达式为
η-η∞η'0-η∞=1[1+(kω)2]m/2 (3)
类似于Cross模型的简化, 式(3)简化为
η=η'0[1+(kω)2]m/2 (4)
通常采用式(2)、(4)进行模型拟合, 以获得ZSV, 拟合的结果见图 4、5。
对比图 4、5可知, Cross模型拟合的ZSV值明显大于Carreau模型拟合的ZSV值, 从模型的本构方程上分析, k与剪切速率ω的乘积在剪切速率趋近于0时, 其值很小, 对2个模型的参数拟合影响不显著。当m趋近于0时, 沥青结合料的流变特性趋近于牛顿流体。而当0 < m < 1时, 2种模型拟合的粘度都是随着剪切速率的减小而增大, 称之为假塑性流体, 此时Cross模型本构方程右侧的分母值要大于Carreau模型的, 在对相同的η进行拟合时, Cross模型拟合的η'0也相应地要大于采用Carreau模型拟合的η'0。美国通过加速加载试验比较高温评价指标的研究表明[8]: 采用Carreau模型获得的ZSV值与加速加载试验的车辙相关系数达到0.54, 要优于采用Cross模型的ZSV值。由于改性沥青和高模量沥青的弹性比例较大, 通过试验不能直接获得ZSV值, 推荐通过动态加载Carreau流变模型进行参数拟合ZSV来评价这种沥青的高温性能。
2.2 静态单循环蠕变恢复
图 6为Burgers的流变模型, 它由一组Maxwell模型G0、η0和一组Kelvin模型G1、η1所组成。
Burgers模型的本构方程为
γ=τ0G0+τ0G1(1-e-tG1/η1)+tτ0η0=γe+γde+γv (5)
式中: γ为剪应变; τ0为恒定的剪应力(Pa); G0为Maxwell模型中的弹性模量(Pa); G1为Kelvin模型中的弹性模量(Pa); η1为Kelvin模型中粘性系数(Pa·s); η0为Maxwell模型中的粘性系数(Pa·s); t为蠕变时间, 即加载时间; γe、γde、γv分别为瞬时弹性应变响应、延迟弹性应变响应和粘性应变响应。
为了求解η0, 可将式(5)两边同除以常量τ0, 方程变形为
γτ0=1G0+1G1(1-e-tG1/η1)+tη0 (6)
若以柔量的方式表达, 可得模型总粘度为
J(t)=Je+Jde(1-e-t/Jdeη1)+Jv (7)
式中: Je为弹性部分, 称为瞬时或玻璃态剪切柔量; Jv为剪切模量粘性部分; Jde为延迟剪切柔量, 其用离散延迟谱拟合时, 通常记为φ(t), 有
φ(t)=n∑i=1Jde(1-e-t/λi) (8)
式中: λi为第i次延迟时间。
离散延迟谱可以考虑更复杂的延迟弹性表现, 常用来描述线粘弹性材料的剪切柔量。相对于Burgers模型, 较为简单的离散延迟谱模型包含了2个Kelvin和1个Maxwell串联, 即延迟弹性部分采用了2个Kelvin元件。
通过变换后的Burgers模型本构方程就是应变与时间t的关系, 对试验数据进行拟合, 就可以获得零剪切粘度。对Burgers流变模型进行分析时, 若通过无限长的时间作用后, Burgers模型中的Maxwell的粘性响应部分会无限接近稳定的粘性流动状态, 即变形速率基本上不随时间变化, 此时得到的粘度参数也可作为零剪切粘度, 通过长时间单循环的蠕变恢复试验可以实现。相关研究表明只有剪应力不超过牛顿流体范围, 才能达到稳定的粘流状态, 而只有在蠕变应力不高于30 Pa时, 沥青才能够获得足够低的剪切速率进入牛顿区域, 因此, 试验采用的应力水平为30 Pa, 温度定为60 ℃, 加载时间分别定为5 h蠕变-5 h恢复和15 h蠕变-15 h恢复2种方式, 对改性和高模量沥青进行研究时, 2种沥青长时间蠕变恢复的应变响应见图 7。
由蠕变恢复试验计算ZSV, 即Burgers模型参数中Maxwell元件中的η0时, 可以从加载和卸载2个阶段分别获得。对加载阶段的蠕变响应进行Burgers流变模型拟合, 由Burgers模型的本构方程可知, 参数的最后一项可求得η0; 从卸载阶段的应变恢复计算时, 恢复阶段可以视作在加载阶段的基础上又施加了一个反向的应力。由Boltzman叠加原理, 结合式(7)的柔量表达式, 可以得到卸载阶段的柔量曲线表达式为
J(t)=t0/η0+Jde[e-(t-t0)/Jdeη1-e-t/Jdeη1] (9)
长时间恢复, 即t-t0趋于无穷大时, 式(7)可整理为
J(t)=t0/η0 (10)
式(10)两边同乘以τ0得
runr=τ0t0/η0 (11)
式中: runr为一次蠕变恢复作用下的不可恢复变形, 即粘性流动变形; t0为初始时刻。
表 1列出了采用Burgers模型和离散延迟谱模型拟合得到的ZSV。
表 1 长时间蠕变恢复试验拟合的ZSVTable 1. Fitting ZSV from long-time creep-recovery test沥青 蠕变-恢复时间/h 蠕变阶段η0/(Pa·s) 恢复阶段η0/(Pa·s) 离散延迟谱η0/(Pa·s) 埃索改性 5-5 6.35×104 6.10×104 7.32×104 15-15 7.76×104 7.67×104 8.98×104 高模量 5-5 3.49×105 3.04×105 4.43×105 15-15 4.52×105 4.17×105 5.06×105 从表 1可以看出, 无论是5 h加载、5 h卸载, 还是15 h加载、15 h卸载的试验时间, 采用离散延迟谱拟合的η0均大于Burgers模型的, 而从蠕变阶段得到的η0又大于从恢复阶段得到的。随着加载时间的增加, 拟合所得到的η0也在增大, 埃索改性沥青15 h蠕变得到的η0比5 h得到的η0增加了22.2%, 高模量沥青的增加了29.5%。可知随着时间的增加, 沥青的粘度变化也很大, 在5-5 h的加、卸载时间内沥青尚未达到稳定的粘性流动状态, 而15-15 h的加、卸载时间是否使沥青达到了稳流状态也不能判定, 从这种静态加载的拟合原理来看, 由于延迟弹性的影响, 如果要把η0从沥青的粘弹性质中完全分离出来, 则需要无限长的作用时间。对比上面提到的采用动态加载模式拟合得到的η0可知, 采用静态加载拟合得到的ZSV均小于动态加载的ZSV, 但采用Carreau模型拟合的ZSV与采用离散延迟谱的拟合15-15 h的ZSV比较接近, 也就是说, 采用离散延迟谱对静载模式下进行ZSV拟合与Carreau模型对动载模式下进行拟合的ZSV相关性最好, 并且除了动态加载的Cross拟合模型外, 其余4个模型的拟合ZSV基本在一个数量级上。由上述对动态加载的分析可知, 离散延迟谱和用Carreau模型拟合的ZSV可以相互验证, 并且与混合料的相关性较好, 可以共同作为评价沥青的高温指标。
3. 结语
本文对不同类型的沥青进行了动、静态加载模式试验, 基于流变模型分析, 得到以下结论。
(1) 采用动态加载和静态加载模式拟合ZSV作为评价沥青高温性能的指标, 由于其加载的模式和拟合的原理不同, 因此, 获得的ZSV也有差异。
(2) 动态加载模式下, 基质沥青可以通过图形外延直接获得ZSV, 改性沥青由于延迟弹性的影响, 需要借助流变模型进行拟合。
(3) 动态加载的Carreau模型拟合结果与静态加载的15-15 h离散延迟谱模型拟合的ZSV相近, 二者具有较好的相关性, 可以共同作为评价改性沥青的高温路用性能指标。
-
表 1 原始数据
Table 1. Primal Data
港口序号 x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 Y 1 27.00 40.40 3.00 9.50 5.00 0.10 1.50 57.86 5.00 5.00 99.98 29.20 2 11.30 8.60 3.00 14.00 2.00 0.25 1.20 30.74 4.70 4.18 99.98 3.20 3 9.80 80.20 1.50 20.20 1.00 0.14 1.60 41.20 5.00 4.55 99.98 4.80 4 24.20 77.40 0.50 4.20 2.00 0.24 2.10 13.98 4.50 4.06 99.98 6.40 5 54.00 94.80 3.00 7.40 8.00 0.07 8.10 45.70 4.90 4.44 99.92 39.00 6 16.10 139.00 2.90 3.60 7.00 0.36 0.68 14.00 4.80 4.50 99.92 23.00 7 28.00 21.40 5.00 13.60 6.00 0.20 17.60 21.80 5.00 4.21 99.92 3.00 8 25.00 73.30 3.50 63.70 1.30 0.20 2.90 42.00 5.00 5.00 99.90 37.00 9 13.30 11.30 3.50 40.80 18.00 0.10 0.50 32.80 4.60 4.11 99.90 40.00 表 2 旋转后因子贡献及贡献率
Table 2. Contributivenesses and Contributiveness Ratios of Factors
1 2 3 4 5 特征值 2.434 3 2.151 0 2.054 2 2.027 3 1.483 6 贡献率 22.129 7 19.554 3 18.674 6 18.430 0 13.487 7 累计贡献率 22.129 7 41.684 0 60.358 6 78.788 6 92.276 3 表 3 载荷矩阵
Table 3. Load Matrix
x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 1 0.032 2 0.181 9 0.196 8 0.399 2 -0.541 3 -0.072 5 -0.013 0 0.591 3 0.847 8 0.913 9 0.008 3 2 0.031 0 -0.030 9 0.583 5 0.680 8 0.666 1 -0.079 5 0.040 6 -0.012 1 0.023 6 0.089 8 -0.940 6 3 0.528 2 -0.274 8 0.009 5 0.035 4 0.349 4 -0.955 6 0.011 5 0.776 5 0.151 7 0.187 4 0.042 1 4 0.522 5 -0.176 3 0.652 3 -0.265 7 0.008 9 -0.071 6 0.971 0 -0.103 1 0.421 5 -0.153 4 -0.259 4 5 0.576 7 0.909 2 -0.318 5 -0.336 3 -0.058 6 0.136 0 -0.026 1 -0.122 9 0.034 6 0.180 5 -0.195 9 表 4 因子得分矩阵
Table 4. Factor Score Matrix
x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 1 -0.095 2 0.069 6 0.096 8 0.190 1 -0.294 2 0.093 5 -0.012 3 0.166 0 0.355 8 0.375 2 0.018 1 2 -0.023 0 0.102 6 0.188 9 0.346 8 0.322 0 0.031 5 -0.100 9 -0.050 7 -0.051 0 0.064 5 -0.472 6 3 0.289 7 -0.099 9 -0.122 9 -0.075 9 0.229 4 -0.502 2 -0.066 8 0.349 1 -0.063 1 -0.003 7 0.077 2 4 0.213 0 -0.112 3 0.291 7 -0.216 5 -0.096 4 0.035 5 0.518 0 -0.109 0 0.206 1 -0.120 8 -0.018 9 5 0.421 7 0.616 7 -0.210 2 -0.182 9 0.075 0 0.029 0 -0.056 1 -0.062 6 -0.029 7 0.103 5 -0.211 0 表 5 因子得分
Table 5. Factor Scores
港口序号 1 2 3 4 5 6 7 8 9 F1 0.984 1 -0.312 9 0.757 3 -1.270 2 -0.292 8 -0.104 4 -0.114 5 1.689 1 -1.335 6 F2 -0.745 4 -0.817 4 -0.914 9 -1.151 2 0.090 0 0.496 6 0.133 6 1.218 9 1.689 9 F3 1.206 2 -0.538 7 0.076 5 -0.422 9 1.535 0 -1.612 1 -0.590 0 -0.421 6 0.767 5 F4 -0.212 3 -0.186 0 -0.675 5 -0.598 7 0.843 7 -0.302 0 2.343 0 -0.407 8 -0.804 3 F5 -0.290 8 -1.342 2 -0.219 8 0.591 0 1.554 2 1.365 5 -0.746 2 0.025 7 -0.937 4 表 6 综合因子得分
Table 6. Integrated Factor Score
港口序号 1 2 3 4 5 6 7 8 9 综合因子得分F 0.237 3 -0.590 6 -0.163 8 -0.667 3 0.655 2 -0.106 8 0.240 3 0.500 4 -0.104 6 事故数Y 29.20 3.20 4.80 6.40 39.00 23.00 3.00 37.00 40.00 -
[1] 何晖光, 吴兆麟, 方祥麟. 海上交通环境的综合评价方法[J]. 大连海事大学学报, 1997, 23 (3): 36-41. https://www.cnki.com.cn/Article/CJFDTOTAL-DLHS703.008.htmHE Hui-guang, WUZhao-lin, FANG Xiang-lin. Integrated Assessment Method of Marine Traffic Environment[J]. Journal of Dalian Maritime University, 1997, 23 (3): 36-41. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DLHS703.008.htm [2] 王芳. 主成分分析与因子分析的异同比较及应用[J]. 统计教育, 2003, 5 (5): 14-17. https://www.cnki.com.cn/Article/CJFDTOTAL-TJJY200305005.htmWANG Fang. Comparison and Application of Principal Component Analysis and Factor Analysis[J]. Statistical Education, 2003, 5 (5): 14-17. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TJJY200305005.htm [3] 于建英, 何旭宏. 数据统计分析与Spss应用[M]. 北京: 人民邮政出版社, 2004. [4] 任中奇. 多元统计分析在多指标综合评价问题中的应用[J]. 辽宁工程技术大学学报: 自然科学版, 1994, 13 (1): 97-100. https://www.cnki.com.cn/Article/CJFDTOTAL-FXKY401.026.htmREN Zhong-qi. The Application of Multiple Statistics Analysisin Multiple Indexes Comprehensive Evaluation[J]. Journal of Liaoning Technology University: Natural Science Edition, 1994, 13 (1): 97-100. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-FXKY401.026.htm [5] 岑况. 因子分析中变量和样本空间的统一[J]. 地质科技情报, 1994, 13 (4): 93-97.CEN Kuang. Combination of Variable and Object Spacein Factor Analysis[J]. Geological Science and Technology Information, 1994, 13 (4): 93-97. (in Chinese) [6] 郑士源. 国际干散货航运市场的评价[J]. 交通运输工程学报, 2004, 4 (4): 89-92. http://transport.chd.edu.cn/article/id/200404022ZHENG Shi-yuan. Evaluation of International Dry Bulk Shipping Market[J]. Journal of Traffic and Transportation Engineering, 2004, 4 (4): 89-92. (in Chinese) http://transport.chd.edu.cn/article/id/200404022 [7] 藤井弥平. 海上交通工学[M]. 东京: 海文堂, 1981. [8] 郑中义, 吴兆麟, 杨丹. 港口船舶事故致因的灰色关联分析模型[J]. 大连海事大学学报, 1997, 23 (2): 61-64. https://www.cnki.com.cn/Article/CJFDTOTAL-DLHS702.014.htmZHENG Zhong-yi, WU Zhao-lin, YANG Dan. Analysis Model of Accident's Main Causes on Port Vessels Incidence by Grey System Theory[J]. Journal of Dalian Maritime University, 1997, 23 (2): 61-64. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DLHS702.014.htm -

计量
- 文章访问数: 351
- HTML全文浏览量: 133
- PDF下载量: 330
- 被引次数: 0