抽油机方案计算.ppt

图为悬点pr随theta变化 x=acos((L.^2-1.2043)./(5.64*L)); fi=x+beta; pr=3.935.*(0.372*pi-fi); （3）扭矩因数 图为扭矩因数tf随theta变化 tf=-1.064*3.935.*sin(theta3-theta2)./(2.82*sin(theta3-theta4)); 3.3等效驱动力矩 等效驱动力矩 图为等效驱动力矩Med随n的变化 Med=151.3508.*Md; 3.4悬点载荷分析 1,2曲线是抽油机上冲程 3,4曲线是抽油机下冲程 图中 —抽油杆在液柱中的重量36000N —抽油杆液面以上整个柱塞面积上的液柱重量30000N —悬点冲程长度m —抽油杆在液体载荷 作用下的静变形， 0.2m 图中四点左边如下 由曲线形状得1,2,3,4段曲线的方程为： 特s， ， ， 带入公式①得如下结果： 其中以下死点为位移零点，向上为位移正方向，任意时间悬点的位移 为： 已知如下参量： 将 带入③得 图为悬点载荷F随theta的分段变化 F=(180000*pr+30000).*(theta=0theta=0.1378*pi)+(66000).*(theta=0.9631*pitheta=0.1378*pi)+(180000*pr-489025.17).*(theta=0.9631*pitheta=1.2136*pi)+(30000).*(theta=1.2136*pitheta=2*pi); 图为L随theta的变化曲线 theta=0:0.005:2*pi L=sqrt(19.3223-9.0759*cos(2.219*pi-theta)) 3.6等效力矩、等效阻力矩 因为构件重量对结果影响不大，故为了简化方便计算可以忽略全部构件的重量: 又因为 ①AB 段时 ②BC段时 ③CD段时 ④DA段时 图为等效阻力矩Mef与theta的分段关系 Mef=((180000.*pr+30000)*3.935*1.064/2.82.*sin(theta3-theta2)./sin(theta3-theta4)).*(theta=0theta=0.1378*pi)+(66000.*sin(theta3-theta2)./sin(theta3-theta4)*3.935).*(theta=0.9631*pitheta=0.1378*pi)+((180000*pr-489025.17)*3.935.*sin(theta3-theta2)./sin(theta3-theta4)).*(theta=0.9631*pitheta=1.2136*pi)+(30000*3.935.*sin(theta3-theta2)./sin(theta3-theta4)).*(theta=1.2136*pitheta=2*pi); 1）简化后模型 曲柄摇杆机构简图如右图所示 2）等效转动惯量Je 将抽油机中所有运动件质量（从电动机转子到抽油机悬点）向曲柄处简化，用作用于曲柄轴上的转动惯量代替系统中的所有质量的转动惯量，这样简化要满足动能相等条件，即： : : 图为Je随theta的变化 Je=107549.7594+102195.885.*(1.064/2.82.*sin(theta3-theta2)./sin(theta3-theta4)).^2; 3.8系统运动微分方程 非保守系统 lagrange方程 注意：角速度只是时间的函数。 整理的等效构件运动微分方程为 将 : ②ＢＣ段时: ③ＣＤ段时： ④ＤＡ段时: 图为 关于 t的运动曲线 图为 关于 的运动曲线 p=(6*2*pi/10+3.55)/2 h=0.0001 n=0 p_rec(1)=p; x_rec(1