路面设计原理与方法 7.积分计算子程序(同前) 8.余项计 在某一区段内无穷积分积的执行程序可见多层体系计算程序。在程序中余项值的计算由控 制变量INTT控制执行。当INTT值为零时,表明计算点在层状体系顶面(Z=0),则要求计算 余顶值,当INTT为1时,表明计算点在层状体系顶面以下点(Z≠0),则不要计算余项。余 项计算值代表符VS11,RS11,TS11及W11分别表示σ2,σ,τm及ω在有限积分段的余项 值 9.多层体系应力、应变及位移计算程序说明 本程序适用于多层体系结构,对双圆或多圆均布荷载可采用应力迭加原理得到。执行程序如 下 程序中主要变量说明 INTT一计算点位置符(INTT=0,表示要求余项) IC—一积分次数; VS(…),RS(…),TS(…),w(·)分别存放每次积分后o,0,,τ2,w的对应值 DEL—一积分的相对精度 WT(·),D(·)分别为高斯积分宽度和高斯积分点。 其它变量名同子程序SOL 六、双圆或多圆荷载应力的计算 现行的路面设计规范多用双圆均匀荷载,利用本程序进行应力计算时,必须对源程序进行修 正。修正的基本方法是根据应力迭加原理。要求荷轮隙、轮印中心及轮印中心两侧的应力或 位移,可利用计算两点位移或应力迭加的方法 七.弹性多层体系应力、位移分析程序(AP01) 1.程序功能 弹性多层体系应力,位移分析程序适用于N层组成的多层结构体系,具有如下功能 1)适用于多层弹性体系,层数不限,在此最大值定为L=6 2)每个层次的弹性模量和泊桑比不受限制; 3)适用于计算各层体系任意一点的应力,位移计算,可同时算出多个点的应力及位移,计算 点最大值为24点 4)荷载为单园垂直均布荷载,作用于上层顶面 5)对双圆荷载,则利用单圆荷载进行应力迭加。 2.程序的输入输出说明 输入变量 NL一一层状体系的层数(NL≤6) NS一一应力、应变及位移计算点数(NS≤24) NIC一一积分最多次数(NIC≤40) INTT一计算点状态参数(INTT=0,表面点,INTT≠0,内部点) DEL—一近似积分的精度,常用0.0001 CR—荷载圆的半径(cm): CP一荷载的单位接触压力Kg/cm^2: R(NS)—一每个计算点离荷载中心的径向坐标值(cm) Z(NS)—一每个计算点离表面的垂直坐标值(cm); E(NL)—一每一层的弹性模量,kg/cm^2 PR(N)——每一层的泊桑比 HA(NH)一一每一结构层的厚度(NH=NL-1)cm 输出变量 STRESS Z——Z方向正应力kg/cm^2: STRESS R—R方向正应力Kg/cm^2: STRESS T—一T方向正应力Kg/cm^2 STRESS ZR—ZR方向剪应力Kg/cm^2 STRAIN Z——Z方向应变 STRAIN R——R方向应变 第44页路面设计原理与方法 第 44页 7．积分计算子程序(同前) 8．余项计算 在某一区段内无穷积分积的执行程序可见多层体系计算程序。在程序中余项值的计算由控 制变量 INTT 控制执行。当 INTT 值为零时，表明计算点在层状体系顶面(Z=0)，则要求计算 余顶值，当 INTT 为 1 时，表明计算点在层状体系顶面以下点(Z≠0)，则不要计算余项。余 项计算值代表符 VS11，RS11，TS11 及 W11 分别表示σz，σr，τzr及ω在有限积分段的余项 值。 9．多层体系应力、应变及位移计算程序说明 本程序适用于多层体系结构，对双圆或多圆均布荷载可采用应力迭加原理得到。执行程序如 下： 程序中主要变量说明 INTT——计算点位置符(INTT=0，表示要求余项)； IC——积分次数； VS(…)，RS(…)，TS(…)，ｗ(·)分别存放每次积分后σz，σｒ，ｚｒ,w 的对应值 DEL——积分的相对精度 WT(·)，D(·)分别为高斯积分宽度和高斯积分点。 其它变量名同子程序 SOL 六、双圆或多圆荷载应力的计算 现行的路面设计规范多用双圆均匀荷载，利用本程序进行应力计算时，必须对源程序进行修 正。修正的基本方法是根据应力迭加原理。要求荷轮隙、轮印中心及轮印中心两侧的应力或 位移，可利用计算两点位移或应力迭加的方法。 七．弹性多层体系应力、位移分析程序(AP01) １．程序功能 弹性多层体系应力，位移分析程序适用于 N 层组成的多层结构体系，具有如下功能： 1)适用于多层弹性体系，层数不限，在此最大值定为 L=6； 2)每个层次的弹性模量和泊桑比不受限制； 3)适用于计算各层体系任意一点的应力，位移计算，可同时算出多个点的应力及位移，计算 点最大值为 24 点； 4)荷载为单园垂直均布荷载，作用于上层顶面； 5)对双圆荷载，则利用单圆荷载进行应力迭加。 ２．程序的输入输出说明 输入变量 NL——层状体系的层数(NL≤6)； NS——应力、应变及 位移计算点数(NS≤24)； NIC——积分最多次数(NIC≤40)； INTT——计算点状态参数(INTT=0，表面点，INTT≠0，内部点)； DEL——近似积分的精度，常用 0.0001； CR——荷载圆的半径(cm)； CP——荷载的单位接触压力 Kg/cm^２； R(NS)——每个计算点离荷载中心的径向坐标值(cm)； Z(NS)——每个计算点离表面的垂直坐标值(cm)； E(NL)——每一层的弹性模量，kg/cm^２； PR(NL)——每一层的泊桑比； HA(NH)——每一结构层的厚度(NH=NL-1)cm； 输出变量 STRESS Z——Z 方向正应力 kg/cm^２； STRESS R——R 方向正应力 Kg/cm^２； STRESS T——T 方向正应力 Kg/cm^２； STRESS ZR——ZR 方向剪应力 Kg/cm^２； STRAIN Z——Z 方向应变； STRAIN R——R 方向应变；