刚性球与平面弹性接触的临界参数计算接触问题作为应用力学的一个分支在工程中会经常遇到实际上，在有机械部 分的工业设备中，几乎无一例外地存在着接触现象典型的例子有:齿轮间的接触， 轴承中滚子与坐圈的接触，凸轮机构中凸轮与传动件的接触，火车的车轮与铁轨的接 触等［1-2］许多工程表面的接触问题，在宏观上一般可以简化为回转体接触，即便 在微观上,实际表面的接触也是椭球状的微凸体接触对于类似这样的接触问题，都 可以简化成球体与平面的接触模型，如果接触过程为弹性变形，则可采用经典的Hertz模型来进行计算处理图1反映了刚性球与 平面接触过程的变形演化趋势当压入深度较小时,材料处于完全弹性接触状态，随 着压入深度的增加，材料内部发生屈服，开始出现塑性变形，当压入深度达到一定值 后,接触区域呈现完全塑性变形在真实的接触过程中，总是希望两体之间的接触处 于弹性状态，此时工件的变形较小，使用寿命也会很高但是当压入深度超过某一值 后,材料就会发生屈服，出现塑性变形，如果工件长期在塑性状态下工作，将会对其使 用寿命产生很大的影响由于材料的弹塑性变形的非线性使得接触问题复杂化，因 此，获得从弹性接触进入弹塑性接触的临界点，即Hertz接触的临界参数显得尤为重要。
图1刚性球与平面接触的变形演化示意图目前，在工程上采用有限元分析方法来仿真接触体的变形、应力分布、接触面 积等得到了广泛应用，但其计算时间长，软硬件成本较高由于弹性接触问题在工程 实际中普遍存在，如何采用一种行之有效的方法进行工件的接触强度分析和校核，建 立符合工程实际的设计和校核公式，一直是工程技术人员和广大科研工作者的一个研究方向因此，我们以弹性接触理论和弹塑性力学为基础，建立刚性球 与平面弹性接触的临界接触参数计算模型,为构件的接触变形分析提供参考1 Hertz接触理论Hertz在研究半径为R的弹性球与弹性半空间变形体的接触时，为了简化分析， 做了如下假设:在压入深度为h时，其接触处于弹性变形阶段，两者之间的接触是无 摩擦的局部变形，其接触区域相对于球体而言很小，则接触面的投影形状为圆形图 2反映了球与半空间之间的弹性接触示意,a表示接触半径因此，两个接触体可以 认为是弹性半空间,在接触区域上受到了相同的接触压力Hertz假设接触区域的压力分布为抛物线形状,其表达式为：//= pa J I - I r /n J-(1)式中,p0为接触中心处的最大接触压力;「为接触点距接触中心的径向距离。
对式(1 )在整个接触区域0 -a区间进行积分，即可得到其接触的合力F为：2 ,(2)接触半径a的表达式为：(3) 式中,E*为等效弹性模量,其表达式为：1 I 一铲 1 -房(4)式中,v、vi分别为两接触体的泊松比,E、Ei为其弹性模量根据Hertz接触理论,半径为R的刚性球与弹性平面接触，其接触力F,最大接触压力p0,接触半径a与压入深度h之间的关系分别是：〜号矿R'W国=2JW/jg J R!i(5)根据Hertz接触理论，对于不同半径的球之间的接触，同样也可以等效为球与平 面的接触，此时半径R为两球的等效半径，其计算表达式为：(6)式中,R1、R2分别为两球的半径,+表示两球为外接触，-表示两球为内接 触2临界参数计算根据图2所示的坐标，采用最大接触压力p0对z轴上各点的应力场进行无量纲化,正应力和剪应力的表达式分别为： q 二 n t 1 ~—=-f !+ / I - — awl.曲f — J "7「I + % J/J (l U = - tr~㈣ 在一孔=耳广上7）式中，Rr、RH和Rz分别是径向、周向和z向的正应力;Srz、SrH和RzH是 相应的剪应力根据式（9广式（11）可以发现:Rz与泊松比v无关，但是Rr和RH与泊松比v相关。
z轴上的无量纲应力Rr/p0、RH/p0和Rz/p0与无量纲位置z/a和泊松比v 之间的变化关系曲线如图3所示图3无量纲应力R/p0、Y/p0与无量纲深度z/a和泊松比v的关系从图3中可以发现，当z/a小于1时,Rr、RH和Rz随深度的增加而增加，但是Rz增加的速度明显没有另外两个快当z/a大于1后,Rr和RH变化很小，且趋近于0,而Rz随z/a增加而继续增加泊松比v越小,Rr和RH越大，这是因为材料的可 压缩性增加（泊松比小），应力也会随之增加在弹性变形区域，应力随着F或h增加而增加，最终导致材料屈服在z轴方向,因此,Rr、RH和Rz为主应力[6]根据Von Mises屈服由于Srz、SrH和SzH为 0,准则，在z轴上的点的屈服应力可以表示为：『二■实三㈣rdf 1+ -l-{ I + ci /1 - —iiJthmf— ) Hft 2(8)对于某一个v,式(13)所表示的屈服应力值随着z/a变化而变化,其变化关系曲 线如图3所示从图3中可以看出:无量纲屈服应力Y/p0随着泊松比v的增加而减 小，而且Y/p0的最大值发生在接触中心的正下方，即最初的屈服发生在材料内部 当该点的屈服应力达到屈服强度Ry时，此时的压入深度hy即为初始屈服时的临界 压入深度。
为获得临界屈服位置，将式(13)相对于z/a进行微分，并令该值为0,即 可得到最大Von Mises应力下的无量纲位置F0=z0/a] A 3 %|I 1+ —..「= 0% S+ f 1+ 0 (9)将式(9)改写为：踞】r I 门 i [f I + J 1 + %(10)对于大部分金属材料，泊松比v在0到0.5之间变化由于式(5)的表达式比较复杂，因此采用数值方法求解该超越方程图4所示为该等式的曲线表达式及 其拟合曲线，从图4中可见，在初始屈服时，其拟合函数关系式为：r UI5 监~ 1一 151 5(11)式(11)的拟合误差小于0.31%因此，无量纲的初始屈服位置F0和泊松比v之间可以用近似线性函数表示：%=壮3幻7+ r丹.016 9(12)从等式(12)中可以看出,F0随着v增加而增加因此，大的v导致初始屈服发生在接触中心正下方更深的位置025 心亍 0.55 065(1« 【搭 0.4 0.3 、Q.2G-02 -44 0.6eh'xt图4泊松比v与无量纲屈服点z0/a的关系当初始屈服发生时q二 P(/£「14 瑞厂'■ r i.+ r i/i ■%1皿侦由｛ r 〃勺J(13) 式中,Fy为弹塑性材料的屈服强度。
布I十 r//l-Fj'b'1(14) 贝p片r.：,耿(15)同样，式(14)也是一个超越方程Cv与v的关系曲线及其拟合曲线如图5所示,拟合函数的表达式为：C「= 0. S42 4尸+ Q S79 &r + 1-3(1) 5(16)对于大部分金属材料,v=0.3则根据式(16)可以计算得到(17) 因此，对于v=0.3的材料，在初始屈服时，临界压入深度hy为：(18) 将式(18)带入到式(5)中，得到其临界压入载荷Fy为:叭=门险点3 h I u ill AHiM ibi^^ iw to ri h h A(19) 同样地，根据式⑺得到其临界接触半径ay为皿二 2. 306 4 Bd I ri Hi(20) 从式(18广式(20)可看出，只要知道材料的弹性模量、屈服强度和球半径,即可计算出刚性球与平面弹性接触时的临界压入深度、临界压入载荷和临界接触半径3 2 i uT id―iri G^- ri's图5 Cv与泊松比v的关系从上面的分析还可以看出，对于泊松比v=0.3的材料，当&Rz-Rr&=0.62p0时，在接触区域中心正下方深度为0.481 1a处,材料开始屈服，此时接触中心处的 接触压力p0=1.613Ry。
如果压入深度继续增加，则接触中心下方的材料塑性变形区 域逐渐扩大，进入弹塑性变形阶段在此变形阶段，由于存在塑性变形，此时，两体 的接触不再满足Hertz弹性接触的使用条件，因此，不能再用Hertz接触理论来分析 两体的接触过程需要指出的是，由于上面分析的是球与平面的弹性接触，因此，弹性模量和泊松 比对接触响应起主要作用当材料发生屈服后，屈服强度就会对接触响应产生较大 的影响在实际工程应用中，很多工件的表面在加工过程中或多或少会产生加工硬 化，在此情况下，可以认为表层材料的屈服强度得到提高，因此，采用上面的计算公式 时，屈服强度需要进行必要的修正，即根据硬化程度而采用合适的等效屈服强度，同 样可以采用式（18广式（20）计算获得弹性接触时的临界参数3结论以弹性接触理论和弹塑性力学为基础，通过分析刚性球与平面弹性接触时的应 力变化关系，采用数值方法获得了初始屈服发生时的位置位于接触中心的正下方 0.481 1a处，此时接触中心处的接触压力为1.613Ry根据Hertz接触理论，建立了 两体弹性接触的临界接触参数计算公式，如果已知材料的弹性模量、屈服强度（或等 效屈服强度）和球半径，即可获得发生初始屈服时的临界参数,为构件的弹塑性接触 变形分析提供了参考。
电厂分散控制系统故障分析与处理作者：单位：摘要:归纳、分析了电厂DCS系统出现的故障原因，对故障处理的过程及注意 事项进行了说明为提高分散控制系统可靠性，从管理角度提出了一些预防措施建 议，供参考关键词:DCS故障统计分析预防措施随着机组增多、容量增加和老机组自动化化改造的完成，分散控制系统以其系 统和网络结构的先进性、控制软件功能的灵活性、人机接口系统的直观性、工程设 计和维护的方便性以及通讯系统的开放性等特点，在电力生产过程中得到了广泛应 用，其功能在DAS、MCS、BMS、SCS、DEH系统成功应用的基础上，正逐步向MEH、 BPC、ETS和ECS方向扩展但与此同时，分散控制系统对机组安全经济运行的影 响也在逐渐增加；因此如何提高分散控制系统的可靠性和故障后迅速判断原因的能 力，对机组的安全经济运行至关重要本文通过对浙江电网机组分散控制系统运行 中发生的几个比较典型故障案例的分析处理，归纳出提高分散系统的可靠性的几点 建议，供同行参考1考核故障统计浙江省电力行业所属机组，目前运行的分散控制系统，有TELEPERM-ME、 MOD300，INFI-90，NETWORK-6000，MACS?和 MACS-?，XDPS-400，A/I。
DEH 有 TOSAMAP-GS/C800， DEH-IIIA等系统笔者根据各电厂安全简报记载，将近几年 因分散控制系统异常而引起的机组故障次数及定性统计于表1书状1*遍铮炉域4二胃伸弊4*n井11m i1£an2WI平312a：w表1热工考核故障定性统计2热工考核故障原因分析与处理根据表1统计，结合笔者参加现场事故原因分析查找过程了解到的情况，下面 将分散控制系统异常(浙江省电力行业范围内)而引起上述机组设备二类及以上故障 中的典型案例分类浅析如下：2.1测量模件故障典型案例分析测量模件“异常”引起的机组跳炉、跳机故障占故障比例较高，但相对来讲故 障原因的分析查找和处理比较容易，根据故障现象、故障首出信号和SOE记录，通 过分析判断和试验，通常能较快的查出“异常”模件这种“异常”模件有硬性故 障和软性故障二种，硬性故障只能通过更换有问题模件，才能恢复该系统正常运行; 而软性故障通过对模件复位或初。