Numerical method of multi-layer elastic system by using stiffness matrix method
-
摘要: 从弹性力学基础理论出发, 采用刚度矩阵法, 推导了应用于直角坐标系下的三维多层弹性层状体系静力学数值解法。引入二维傅里叶变换及高斯积分求解法, 基于MATLAB数学软件平台编制计算程序, 实现三维多层弹性层状体系理论计算方法的数值求解。针对典型有砟轨道轨下基础结构, 采用提出的计算方法和编制的相应计算程序对其进行静力学分析, 并将所获得的计算结果与采用通用有限元程序ABAQUS的计算结果进行对比。分析结果表明: 采用提出的计算方法和通用有限元计算方法获得有砟轨道轨下基础最大竖向位移分别为1.50、1.95 mm, 最大竖向应力分别为0.34、0.21 MPa, 计算结果较为接近, 计算反映出来的各状态分量变换规律基本一致, 提出的计算方法及其相应计算程序可应用于多层弹性层状体系的静力学计算。Abstract: By using the theory of elastic mechanics and stiffness matrix method, a numerical method of static mechanics for calculating 3D multi-layer elastic system under rectangular coordinate system was developed.2D Fourier transform and Gauss integral method were used, a calculation program was given based on MATLAB mathematical software platform, and the numerical solution of 3D multi-layer elastic system was got.The static mechanics analysis of ballast track foundation structure was carried out by using the numerical method and program, and the calculation results were compared with the results obtained by using general finite element program ABAQUS.Analysis result shows that the maximum vertical displacements of ballast track foundation are 1.50, 1.95 mm respectively by using the numerical method and general finite element calculation method, and the maximum vertical stresses of the foundation are 0.34, 0.21 MPa respectively by using the two methods.The calculation results are close to each other, and the conversion rules of state components are basically the same according to the calculation, so the proposed numerical method and program can be applied to the static mechanics calculation of multi-layer elastic system.
-
车辆是道路交通事故的主要(必然) 当事方, 对道路交通事故进行模拟再现的直接目的是在逆向计算出事故车辆制动前或碰撞前瞬间行驶速度的前提下, 在计算机屏幕上形象再现事故的全过程, 其最终目的是为了帮助交通管理执法人员查证事故原因, 鉴定事故责任。在分析道路交通事故时, 一般都将整个事故过程分为3个阶段[1]
(1) 碰撞(瞬间) 前的运动过程, 当事方都进入事故不可逆时段但尚未碰撞接触之前的运动阶段;
(2) 碰撞接触过程, 当事方发生接触并且有明显力作用的阶段;
(3) 碰撞后运动过程, 当事方脱离碰撞接触之后的自由运动阶段。
目前, 模拟再现道路交通事故的主流技术路线如下:
(1) 首先求解当事车辆在碰撞后运动阶段的瞬态运动状态, 应用经典运动学和动力学理论, 建立车辆运动力学模型(包括车轮-地面力学模型) 和轨迹叠代计算模型等, 根据来自事故现场的碰撞位置和停止位置参数, 逆向计算出碰撞后脱离接触瞬间的车辆速度和角速度。
(2) 根据碰撞模型逆向求解碰撞前瞬间的车辆速度和角速度。
(3) 根据地面痕迹确定事故车辆在发现事故倾向前正常行驶时的速度和路线。
(4) 根据车辆运动力学模型(包括车轮-地面力学模型)、车辆行驶速度、初始位置与路线等, 按设定时间步长正向计算车辆的瞬时状态, 并适时输出到屏幕。
显然, 道路交通事故过程中的第(1) 和第(3) 阶段的定量计算都主要依托于车辆动力学模型, 其模型本身的正确性和精确度对道路交通事故的计算及模拟再现结果正确与否有着决定性的作用。
1. 车辆三维四轮动力学模型
本研究在改进以前的碰撞车辆二维四轮动力学模型[2-3]的基础上, 提出在进行道路交通事故过程中碰撞后运动阶段的动力学计算时, 可采用如图 1所示的碰撞车辆三维简化模型。在此阶段, 碰撞车辆在冲击或运动惯量的作用下, 从碰撞作用结束瞬间位置开始“自由运动”直至车辆最后停止, 此间车辆只受重力和车轮-地面作用力, 而诸如空气阻力等可忽略不计。若假设:
(1) 路面坡度为零, 无(不计) 坑洼不平;
(2) 忽略空气阻力, 忽略轮胎回正力矩的作用;
(3) 整车质心在车体上的位置不因悬架变形而改变。
则由图 1可得下列碰撞车辆动力学三维计算模型
˙Voz=(Fzfh+Fzfm+Fzrh+Fzrm-mg)/m˙Voσ=-[(Fυfh+Fυfm)cosδ+(Fσrh+Fσrm)+(Fζfh+Fζfm)sinδ]/m˙Voξ=-[(Fζfh+Fζfm)cosδ+(Fξrh+Fξrm)-(Fυfh+Fυfm)sinδ]/m
¨ψ=[Ζ0(Fξfh+Fξfm+Fξrh+Fξrm)-Lcf2(Fzfh+Fzfm)+Lcr2(Fzrh+Fzrm)]/(mρ2ψ)
¨φ=[Ζ0(Fσfh+Fσfm+Fσrh+Fσrm)+Bcm2(Fzfm+Fzrm)-Bch2(Fzfh+Fzrh)]/(mρ2ψ)
˙ω=Μ0/(mρ2θ)Μ0=[(Fζfh+Fζfm)sinδ+(Fυfh+Fυfm)cosδ]Lf-(Fσrh+Fσrm)(L-Lf)+B(Fξfm-Fξfh)/2+Μf
Μf=B[(Fζfm-Fζfh)cosδ-(Fυfm-Fυfh)sinδ]/2Fξfh=Fζfhcosδ-FυfhsinδFzfh=[mfg-Κfh(Lsfh-Lsp)+Dfh(Lcf˙ψ+Bch˙φ)]/cos(φ+ψ)Fzfm=[mfg-Κfm(Lsfm-Lsp)+Drh(Lcf˙ψ-Bch˙φ)]/cos(φ+ψ)Fξfm=Fζfmcosδ-FυfmsinδFσfh=Fυfhcosδ+FζfhsinδFσfm=Fυfmcosδ+FζfmsinδFzrh=[mrg-Κrh(Lsrh-Lsp)+Dfm(Lcf˙ψ-Bcm˙φ)]/cos(φ+ψ)Fzrm=[mrg-Κrm(Lsrm-Lsp)-Drm(Lcf˙ψ+Bcm˙φ)]/cos(φ+ψ)
图 2为碰撞车辆动力学简化模型的侧视分析示意图、正视分析示意图, 由图 2可得车辆各相关参数关系式
Lsfm=[ (Z0-rt) secψ-Lftg ψ]secφ+B2tgφ
Lsrm=[ (Z0-rt) secψ+Lrtg ψ]secφ+B2tgφ
Lsfh=[ (Z0-rt) secψ-Lftg ψ]secφ-B2tgφ
Lsrh=[ (Z0-rt) secψ-Lftg ψ]secφ-B2tgφ
Bcm=B2cosφ+(Ζ0+B2sinφ-re)tgφ
Bch=B2cosφ-(Ζ0-B2sinφ-re)tgφ
Lcf=Lfvcosψ- (Z0-Lfsinψ-re) tgψ
Lcr=Lrvcosψ+ (Z0+Lrsinψ-re) tgψ
模型中诸符号意义为: Z、X、Y为车辆质心在地面坐标系中的位置坐标值(m); ψ为过质心的车体横截面与Z轴向所成角度值(俯仰角, 逆时针为正) (rad); φ为车体纵向中心面与Z轴向所成角度值(侧倾角, 逆时针为正) (rad); θ为车体纵向中心面与X轴向所成角度值(横摆角, 逆时针为正) (rad); ω为车体横摆角速度(rad/s); δ为转向角度值(逆时针为正) (rad); V为速度(m/s); F为地面对轮胎的作用力(N); M为作用力矩(N·m); ρ为车体横摆转动惯性半径(m); m为汽车质量(kg); K为悬架等效刚度(N/m); D为悬架阻尼系数(Ns/m); L为轴距(m); B为轮距(m); Lf为车体质心距前轴距离(m); Lr为车体质心距后轴距离(m); Lsp为车辆静止状态时质心距车轴垂直平面的距离(m); re为车轮滚动半径(m)。
上述模型中各符号的下标意义: z、x、y为Z轴向、X轴向、Y轴向分量; f为前轴的值; r为后轴的值; o为质心处的值; σ、ξ为σ轴向、ξ轴向分量; υ、ζ为前轮有转向角δ时的车轮侧向及纵向; h为左侧车轮的值; m为右侧车轮的值; φ为绕ξ轴转动的量; ψ为绕σ轴转动的量; θ为绕Z轴转动的量。碰撞车辆运动力学计算模型中地面对轮胎的作用力可由简化的G.Gim and P.E Nikavesh轮胎理论模型[4-6]求得, 其纵向力特性(Fx-ss) 和横向力特性(Fy-α) 分为下列3种情况表述。
上述模型中诸符号意义: T为车轮的侧偏角(rad); V为车轮的侧顷角(rad); ss为纵向滑移率; ssc为纵向临界滑移率; sV为由侧顷角V引起的车轮滑移率(sV=sinV); sVc为sV的临界值; sT为由侧偏角引起的侧向滑移率(sT=tanT); sTc为侧向临界滑移率; l为轮胎与地面的接触面长度(m); la为纵向接触面长度中附着区长度(m); ln为附着接触面长度的无量纲值(ln=la/l); Cs为轮胎纵向刚度(N); CT为轮胎横向刚度(N); CV为轮胎侧顷刚度(N); _y为横向附着系数;
2. 三维模拟模型的计算误差界定
通常情况下, 实际道路交通事故现场实测数据, 尤其是由调查或简单计算而得的碰撞车速数据都是不准确的, 不宜用来精确界定模拟再现模型的计算误差。本研究引用参考文献[7]所示的由日本汽车研究所提供的16起(T1~T16) 在可精确控制和测定条件下进行的车对车实车碰撞试验的条件数据, 采用本文提出的车辆动力学三维模型进行逆向模拟计算。16起车对车碰撞试验的车辆碰撞车速计算值、实验实测值及其计算误差见表 1, 悬架刚度为13500 N/m, 悬架阻尼系数为1500 N·s/m。16例中各车辆计算碰撞车速相对试验测定值的最大绝对误差值为8.6 km/h (T7例); 计算碰撞车速的最大平均相对误差值为11.58% (T8例), 总体平均相对误差值为6.65%;甲车计算碰撞车速的最大相对误差值为14.15% (T3例); 乙车计算碰撞车速的最大相对误差值为14.29% (T7例)。
表 1 计算车速与实测车速比较Table 1. Comparison of computation vehicle speeds and measured vehicle speeds3. 三维模拟模型与二维四轮模拟模型使用效果比较
将采用车辆三维动力学模型进行模拟计算时的碰撞车速平均相对误差值与采用车辆二维动力学模型进行模拟计算时的碰撞车速平均相对误差值(根据参考文献[7]的模拟结果计算而得) 进行比较, 其结果如图 3所示, 16例中有10例(占62.5%) 是三维模型的平均相对误差值大于二维四轮模型的平均相对误差值; 三维模型的总体平均相对误差值(6.65%) 比二维四轮模型的总体平均相对误差值(5.22%) 大1.43%;三维模型的最大相对误差(T7例乙车14.29%) 比二维四轮模型的最大相对误差[6] (T9例乙车11.20%) 大3.09%。
上述16例针对实车碰撞实验数据的模拟计算结果比较表明, 在目前采用的假设条件及模型结构下, 三维模型的模拟计算精度会略差于二维四轮模型的模拟计算精度。
4. 采用三维计算模型对道路交通事故的模拟再现
在本验证比较例中, 采用三维模型比采用二维四轮模型时总体平均计算精度有所下降, 但对含有车辆滚翻及坠崖等空间运动的事故形态进行模拟再现时, 必须采用车辆三维模型。图 4为采用本文的车辆三维计算模型对一起含有车辆侧翻运动形式的实际道路交通事故的模拟再现画面帧。
作为模拟计算初始条件的是事故前车辆的行驶路线、位置和速度, 这些参数是在此前的事故逆向分析计算阶段求得: 甲车为以31.1 km/h的车速由西向北左转行驶的轻型货车, 乙车为以73.6 km/h的车速由东向西行驶的小轿车。
5. 结语
本文提出的用于道路交通事故模拟再现的车辆三维计算模型的总体(针对16例实车实验数据) 平均相对误差界定为6.65%, 相对本研究以前所用的二维四轮模型, 其模拟计算的总体平均精度下降了1.43%, 但考虑到采用车辆三维计算模型时在道路交通事故形态包容性方面的巨大优势, 计算精度的适度降低是可以接受的。
-
表 1 有砟轨道计算参数
Table 1. Calculating parameters of ballast track
表 2 最大竖向位移对比
Table 2. Comparison of maximum vertical displacements
-
[1] BURMISTER D M. The general theory of stresses and displacements in layered systems. Ⅰ[J]. Journal of Applied Physics, 1945, 16(2): 89-94. doi: 10.1063/1.1707558 [2] BURMISTER D M. The general theory of stresses and displacements in layered soil system. Ⅱ[J]. Journal of Applied Physics, 1945, 16(3): 126-127. doi: 10.1063/1.1707562 [3] BURMISTER D M. The general theory of stresses and displacements in layered soil system. Ⅲ[J]. Journal of Applied Physics, 1945, 16(5): 296-302. doi: 10.1063/1.1707590 [4] 王凯. N层弹性连续体系在圆形均布荷载作用下的应力与位移[J]. 土木工程学报, 1982, 15(2): 65-76.WANG Kai. Stress and displacements analysis of an N-layered elastic-continuous system under vertical load uniformly distributed on a circular area[J]. China Civil Engineering Journal, 1982, 15(2): 65-76. (in Chinese). [5] 王凯. N层弹性连续体系在双圆均布复合荷载作用下的力学计算[J]. 固体力学学报, 1983(1): 136-153. https://www.cnki.com.cn/Article/CJFDTOTAL-GTLX198301017.htmWANG Kai. Calculation of stress, strains and displacements in an N-layered elastic system under the multiple inward horizontal loads on circular areas[J]. Acta Mechanica Solida Sinica, 1983(1): 136-153. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GTLX198301017.htm [6] 王凯. N层弹性体系在多圆向心水平荷载作用下的力学计算[J]. 重庆交通学院学报, 1984(2): 50-64. https://www.cnki.com.cn/Article/CJFDTOTAL-CQJT198402005.htmWANG Kai. Calculation of stress, strains and displacements in an N-layered elastic system under the multiple inward horizontal loads on circular areas[J]. Journal of Chongqing Jiaotong University, 1984(2): 50-64. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-CQJT198402005.htm [7] 王凯, 姚炳卿. N层弹性体系在多圆旋转水平荷载作用下的力学计算[J]. 西安公路学院学报, 1986, 4(3): 15-30. https://www.cnki.com.cn/Article/CJFDTOTAL-XAGL198603001.htmWANG Kai, YAO Bing-qing. Calculation of an N-layered elastic system under the action of multiple rotating-horizontal loads distributed on the circular areas[J]. Journal of Xi'an Institute of Highway, 1986, 4(3): 15-30. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-XAGL198603001.htm [8] 张起森. 弹性层状体系理论的实验验证及应用[J]. 土木工程学报, 1985, 18(4);63-76. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC198504007.htmZHANG Qi-sen. Experimental verification of the elastic layer system theory and its application[J]. China Civil Engineering Journal, 1985, 18(4): 63-76. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC198504007.htm [9] 张起森, 郑健龙. 弹性层状体系考虑层间接触非线性的有限单元分析法[J]. 土木工程学报, 1989, 22(3): 63-75. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC198903009.htmZHANG Qi-sen, ZHENG Jian-long. Finite element analysis of elastic layer system considering nonlinear contact condition between layers[J]. China Civil Engineering Journal, 1989, 22(3): 63-75. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC198903009.htm [10] 黄卫. 弹性层状体系弯拉应力近似计算[J]. 岩土工程学报, 1995, 17(6): 52-54. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC506.005.htmHUANG Wei. Regression formulas of tensile stress in the elastic multilayer[J]. Chinese Journal of Geotechnical Engineering, 1995, 17(6): 52-54. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC506.005.htm [11] 任瑞波. 沥青路面结构计算方法与FWD应用技术的研究[D]. 哈尔滨: 啥尔滨建筑大学, 2000.REN Rui-bo. The study of asphalt pavement structure calculation and FWD application technique[D]. Harbin: Harbin University of Civil Engineering and Architecture, 2000. (in Chinese). [12] 任瑞渡, 钟岱辉, 孔军, 等. 沥青路面层状粘弹性半空间轴对称问题的求解[J]. 山东建筑工程学院学报, 2002, 17(4): 1-7. doi: 10.3969/j.issn.1673-7644.2002.04.001REN Rui-bo, ZHONG Dai-hui, KONG Jun, et al. The theoretical method for solving axisymmetrical problems in multilayered viscoelastic asphalt pavement half space[J]. Journal of Shandong Institute of Architecture and Engineering, 2002, 17(4): 1-7. (in Chinese). doi: 10.3969/j.issn.1673-7644.2002.04.001 [13] 钟阳, 王哲人, 郭大智. 求解多层弹性半空间轴对称问题的传递矩阵法[J]. 土木工程学报, 1992, 25(6): 37-43. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC199206004.htmZHONG Yang, WANG Zhe-ren, GUO Da-zhi. The transfer matrix method for solving axisymmetrical problems in multilayered elastic half space[J]. China Civil Engineering Journal, 1992, 25(6): 37-43. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC199206004.htm [14] 钟阳, 王哲人, 郭大智, 等. 求解多层弹性半空间非轴对称问题的传递矩阵法[J]. 土木工程学报, 1995, 28(1): 66-72. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC199501007.htmZHONG Yang, WANG Zhe-ren, GUO Da-zhi, et al. Transfer matrix method for solving non-axisymmetrical problems in multilayered elastic half space[J]. China Civil Engineering Journal, 1995, 28(1): 66-72. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC199501007.htm [15] 钟阳, 张永山. 求解多层弹性半空间轴对称动态问题的精确刚度矩阵法[J]. 力学季刊, 2003, 24(3): 395-400. doi: 10.3969/j.issn.0254-0053.2003.03.018ZHONG Yang, ZHANG Yong-shan. Explicit solution for axisymmetrical multilayered elastic half space problems by exact stiffness matrix method[J]. Chinese Quarterly of Mechanics, 2003, 24(3): 395-400. (in Chinese). doi: 10.3969/j.issn.0254-0053.2003.03.018 [16] 钟阳, 郭大智, 张肖宁. 轴对称弹性半空间问题一般解的新方法[J]. 哈尔滨建筑大学学报, 1995, 28(2): 23-27. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBJ502.003.htmZHONG Yang, GUO Da-zhi, ZHANG Xiao-ning. New way of the general solution of axial symmetry elastic half space problem[J]. Journal of Harbin University of Architecture and Engineering, 1995, 28(2): 23-27. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HEBJ502.003.htm [17] 钟阳, 陈静云, 王龙, 等. 求解动荷载作用下多层粘弹性半空间轴对称问题的精确刚度矩阵法[J]. 计算力学学报, 2003, 20(6): 749-755. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJG200306017.htmZHONG Yang, CHEN Jing-yun, WANG Long, et al. Explicit solution for dynamic response of axisymmetrical problems in multilayered viscoelastic half space by exact stiffness matrix method[J]. Chinese Journal of Computational Mechanics, 2003, 20(6): 749-755. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JSJG200306017.htm [18] HUANG Y H, LIN C, DENG X J, et al. Kentrack: a Computer Program for Hot-mix Asphalt and Conventional Ballast Railway Trackbeds[M]. Lexington and Lanham: Asphalt Institute and NAPA Publication, 1984. [19] ROSE J G, SU B, LONG W B. Kentrack: a railway trackbed structural design and analysis program[C]//AREMA. Proceedings of the AREMA 2003 Annual Conference and Exposition. Chicago: AREMA, 2003: 1-25. [20] ROSE J G, KONDURI K C. Kentrack—a railway trackbed structural design program[C]/7 AREMA. Proceedings of the AREMA 2006 Annual Conference and Exposition. Lanham: AREMA, 2006: 1-31. [21] ABRAMOWITZ M, STEGUN I A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables[M]. New Youk: Dover Publications, 1965. [22] 刘尧军, 郭增强, 赵玉成. 秦沈客运专线不同填料路基动力特性的试验研究[J]. 石家庄铁道学院学报, 2004. https://www.cnki.com.cn/Article/CJFDTOTAL-SJZT200402005.htm17(2): 19-22. LIU Yao-jun, GUO Zeng-qiang, ZHAO Yu-cheng. Research on dynamic characteristics of Qinshen high-speed railway roadbed with different filling materials[J]. Journal of Shijiazhuang Railway Institute, 2004, 17(2): 19-22. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-SJZT200402005.htm [23] WANG Kai-yun, HUANG Chao, ZHAI Wan-ming, et al. Progress on wheel-rail dynamic performance of railway curve negotiation[J]. Journal of Traffic and Transportation Engineering; English Edition, 2014, 1(3): 209-220. [24] 杨尧. 客运专线铁路基床填料动静三轴试验研究[D]. 成都: 西南交通大学, 2009.YANG Yao. Experimental study on typical filler for passenger dedicated railway subgrade roadbed by static and dynamic triaxial test[D]. Chengdu: Southwest Jiaotong University, 2009. (in Chinese). -