大庆石油学院机械系机械设计教研室12轴设计计算程序的编制1.程序的使用范围a)本程序可以对转轴进行受力分析、弯扭合成校核和疲劳强度精确校核。b)可进行多个危险截面的校核,若校核不合格、可根据加大直径或改进材料重新计算,直至合格为止。2.数学模型的建立、数表和线图的公式化a)绝对尺寸系数的拟合公式31793.0558869364.2d(1)式中d为轴的截面直径,mm.b)圆截面钢材的扭转剪切尺寸系数de310867.19513.0(2)式中d为轴的截面直径,mm.c)钢材的敏性系数323232323232016709.0139645.039766.033835.0)()(0139.011988.035439.041319.0)()(0179332.0142175.0377568.046663.0)()(0183968.0143857.0371525.051812.0)()(01257435.01.02711768.06547196.0)()(00886175.0069183.0176398.0814.0)()(rrrqqfrrrqqerrrqqdrrrqqcrrrqqbrrrqqa(3)式中r为圆角半径,mm;序号(a)~(f)代表的材料强度极限,B见表1表1轴的材料强度极限B、B(MPa)序号(a)(b)(c)(d)(e)(f)备注B1400980700560420350计算qB1250840560420--计算q若轴的材料强度极限不等于表1中的B值时,可按插入法计算。d)弯曲疲劳的表面质量系数的拟合公式)(4647815.38)()(000001825987.0)()(0002575974.0)()(000075.0963.0)()(1)(6442821499.0未加工粗车精车磨削抛光BBBBedcba(4)3扭转剪切疲劳的表面质量系数如无试验资料时,可取≈。e)轴上键槽处的有效应力集中系数)(kk的拟合公式0.0019.0k0.0011kBB(5)f)轴肩圆角处的理论应力集中系数)(的数表程序化见子程序DKA和子程序DKB。g)零件与轴过盈配合处的)(kk值的数表程序化见子程序DKA和子程序DKB。h)支反力的计算公式轴上所受外力情况常见的可分为三类(见图1、图2、图3):HHtHHtiHVVriVVriVQRFRLLLLLQLFRRQFRLLMLLLQLFRka11232323311232132331)()()1(1).(见图(1)4HHtHHtiHVVriVVriVQRFRLLQLLFRRQFRLLQMLLFRkb1122321112231211)())()2(2).(见图(2)52112321323211223212132321)())()3(3).(tHtHttiHVrriVVriVFRFRLLLLFLLFRRFFRLLLMMLFLLFRkc见图(3)i)轴的弯矩计算公式轴的弯矩计算公式简图如图4所示。现给出轴上各段的弯矩计算公式。(a)第一段(1Lx)xRMxRMMHHVV11101(4)式中x为所计算的截面到轴左端点的距离,mm,0M为轴左端上的集中弯矩。(b)第二段(211LLxL))()(121112112LxRMMLxRMMMHHHVVV(5)6(c)第三段(32121LLLxLL))()(21323212223LLxRMMLLxRMMMHHHVVV(6)说明:在求K=1,2,3中的某一种情况的弯矩时,应根据轴的受力简图对照轴的弯矩计算简图(图4)做相应的外力代换。如K=3时,式(3)、(4)、(5)中的外力根据轴的受力简图(3)代换为,00M,23VVFR,23tHFR,12tHFR,12VVFR11HHRR。d)总弯矩M的计算公式22ViHiMMMi=1,2,3e)计算弯矩caM的计算公式22)(TMMca式中α为折算系数,见教材第375页;T为扭矩,,/1055.95nPT单位为N·Mpa。3确定轴的强度校核计算程序变量名轴的强度校核计算程序变量名列于表1中。名称程序变量名公式符号单位备注轴传递的功率轴的转速齿轮齿数齿轮模数法向压力角齿轮螺旋角圆锥齿轮分度圆锥角圆锥齿轮齿宽系数带(链)轮的压轴力PNZ1,Z2CM1,CM2ALPHANBETA1BETA2DELTADRQPnz1,z2m1,m2αnβ1,β2δφRQkWr/minmm(°)(°)(°)N实型量法向模数φR=0.25~0.357压轴力与水平面的夹角轮齿径向力轮齿圆周力轴上集中弯矩轴的计算长度轴上扭矩垂直面支反力水平面支反力总弯矩弯矩计算弯矩最大计算弯矩垂直面弯矩折算系数轴上所受外力情况轴上外力方向HVQQ轴上外力Ft1˙Ft2方向轴上外力Fr1˙Fr2方向轴上集中弯矩M1,M2方向`轴径危险截面位置材料强度极限材料许用弯曲疲劳极限材料弯曲疲劳极限材料剪切疲劳极限ALPHAQFR1,FR2FT1,FT2M1,M2L1,L2,L3TRV1,RV2RH1,RH2MXL,MXRMH,MVMC1,MC2MCAMVAALPHAKQ1,Q2E1,E2C1,C2S1,S2DXSIGMBSIGMABSIGM1TAU1Fr1,Fr2Ft1,Ft2M1,M2L1,L2,L3TRV1,RV2RH1,RH2MMH,MVMcaMcaMVαdxσB[σ-1]bσ-1τ-1(°)NNN.mmmmN.mmNNN.mmN.mmN.mmN.mmN.mmmmmmMPaMpaMpaMpa实型量实型量L、R表示有集中弯矩的左、右截面MH,MV分别为水平和垂直面弯矩x处截面右侧弯矩K=1,2,3与轴受力简图同向,赋正1与轴受力简图同向,赋正1与轴受力简图同向,赋正1与轴受力简图同向,赋正1距轴左端点的距离8轴截面上的弯曲应力轴截面上的扭转切应力轴肩圆角半径轴肩直径比弯扭合成应力应力幅平均应力切应力幅平均切应力理论应力集中系数理论应力集中系数敏性系数敏性系数有效应力集中系数有效应力集中系数绝对尺寸系数绝对尺寸系数表面质量系数表面质量系数综合系数强化系数计算公式选择系数数据文件选择系数综合系数材料特性系数材料特性系数法向应力计算安全系数扭转切应力计算安全系数计算安全系数设计安全系数SIGMBBTAUTRDDSIGMCASIGMASIGMMTAUATAUMASIGMATAUQSQTKSKTESETBSBTKSIGMBEATQKBKKKTAUPSPTSSSTSCASσbτTrD/dσσaσaσmτaτmασατqσqτkσkτεσετβσβτKσβqKτφσφτSσSτScaSMpaMpaMmMpaMpaMpaMpaMpa对正应力对切应力对正应力对切应力实型量,对正应力实型量,对切应力对正应力对切应力对正应力对切应力对正应力,实型量见附表3-9KB=1、2、3、4、5的含义见式(4)KK=1、2、3的含义见图1对切应力,实型量对正应力对切应力94.轴的强度校核计算程序框图图1。轴的强度校核计算程序框图输入:Q,P,N,CM1,CM2,Z1,Z2,K,L1L2,L3,ALPHN,BETA,SIGMAB输入:K2,XSIGMCA=10*MCA/D**3计算MC1,MC2,MCASIGMCASIGMAB打印:弯扭合成校核结果BA计算RV1,RV2,RH1,RH2计算FR1,FR2,FT1,FT2,M1,M2,T输入:D开始NY10A输入新的X?(X/Y)是Y?输入:K2K2=?输入:SIGMB,SIGM1,TAU1,BEATQ,SK3=?计算SS,ST,SCA计算ES,ET,KS/ES,KT/ET计算键槽的KS,KTSCA》S打印:计算结果BENDDCEF计算KS,KT计算BS,BT输入:KBNYYN01321图1。轴的强度计算程序设计框图(续1)11C计算qσ,qτ打开ασ数据文件输入:R,DD用子程序DK读ασ用子程序DKB查ασ打开ατ数据文件用子程序DK读ατ用子程序DKB查ατ计算Kτ,KτEDKK=?打开直径50的k数据文件打开直径30的k数据文件打开直径大于100的k数据文件用子程序DK读k(k)用子程序DKB查k(k)F输入:KK,C132图1。轴的强度计算程序设计框图(续2)125.较典型的程序段在数组DK(8,11)中用插值法查找已存入的系数ασ(ατ)或κσ/εσ(κτ/ετ)的子程序DKB。输入:DX为已知的圆角半径与轴径的比值r/d,或轴上零件与轴配合种类;DY为轴肩直径比D/d,或材料的强度极限σB;KD为在数组DK(8,11)中所存元素的最大列数;DK为已读入系数的数组。输出DA、DB为查得的系数,若KD=11时,DA为所求值;否则DB为所求值。6.上机计算1)题目:见教材第238页例题。2)输入数据:启动程序后,由键盘输入以下数据:Q=0.0,ALPHAQ=0.0,P=9.4,N=93.6,CM1=4.0,CM2=0.0,Z1=95.0,Z2=0.0,BETA1=9.11,BETA2=0.0,DELTA=0.0,DR=0.0,L1=106.0,L2=79.0,L3=149.0,K=1,ALPHA=0.59,SIGMAB=650.03)输出计算打印结果如下:FIRSTHANDDATAQ=?P=?N=?ALPHAQ=?CM1=?CM2=?BETA1=?BETA2=?DELTA=?L1=?L2=?L3=?K=?ALPHA=?SIGMAB=?T=?DESIGNRESULDFR1=?FT1=?FR2=?FT2=?M1=?M2=?FV1=?FV2=?FH1=?FH2=?X=?MH=?MV=?MC1=?MC2=?MCA=?D=?SIGMCA=?SIGMAB=?SS=?ST=?SCA=?S=?D=?SIGMCA=?SIGMAB=?SS=?ST=?SCA=?S=?13