文章来源：锻压装备与制造技术.2023,58(05):7-10
文/赵文存，倪欢欢.合肥合锻智能（603011）制造股份有限公司
摘要：介绍多连杆机械压力机中八连杆机械压力机动力学特点。通过Vb.net软件编制专用程序，对八杆机构进行力学分析，绘制受力曲线图，对杆系轴瓦Pv值进行分析，对压力机设计人员具有实际应用价值。
关键词：多连杆机械压力机；八杆机构；动力学分析；Vb.net。
1 引言
多连杆机械压力机结构可靠、性能优良，广泛应用于薄板件的拉延成形工艺。
多连杆机械压力机中，主传动采用八连杆机构的设计其结构较为复杂，设计变量较多，计算难度大。八连杆机构杆系的动力学分析，可先建立数学模型，在计算机中完成计算分析。作者采用Vb.net软件编制八连杆机构动力学分析专用程序，绘制压力机全行程受力曲线，判断杆系受力极值、各轴瓦处PV极值，可为八连杆机械压力机设计提供直接计算依据。
2 动力学分析
八连杆机械压力机，其主传动结构简图如图1所示。偏心体1在偏心齿轮带动下转动，驱动整个传动系统。然后通过上拉杆2、上摇杆3、下摇杆5(上摇杆与下摇杆在摇杆轴4处焊接联接)、下拉杆6、角架7(角架为一个三角形焊接构件)、连杆8带动滑块做往复直线运动，完成冲压过程。
图1.八连杆机构主传动结构简图
将主传动中的结构简化为八杆杆组机构图，如图2所示。可采用解析法对八杆机构进行分析，用复数矢量法列出机构封闭矢量方程。
图2.八杆杆组机构图
将八杆机构进行拆解成3个模块(二级机构)，分别建立方程组，求解方程组得出滑块行程S(mm)
S=L1+L7+L8+L1*Sin(φ1)+L7*Sin(φ7)+L8*Sin(φ8).........................(1)
对滑块行程S求导可得滑块速度V(mm/s)
V=L1*ω1*Cos(φ1)+L7*ω7*Cos(φ7)+L8*ω8*Cos(φ8).......................(2)
对滑块速度V求导可得滑块加速度A(mm/s2)
A=-L1*ω1^2*Sin(φ1)+L7*(a7*Cos(φ7)-ω7^2*Sin(φ7))+L8*(a8*Cos(φ8)-ω8^2*Sin(φ8))..............................................................(3)
φ1—L1杆相对于X坐标轴转角(rad)；φ7—L7杆相对于X坐标轴转角(rad)；
φ8—L8杆相对于X坐标轴转角(rad)；
ω7—L7杆角速度(rad/s)；ω8—L8杆角速度(rad/s)；
a7—L7杆角加速度(rad/s2)；a8—L8杆角加速度(rad/s2)；
偏心体角速度为ω1(rad/s)
ω1=2*π*N/60..........................................................(4)
N—滑块连续行程次数；
要求出杆系受力，必须先求出偏心体1的驱动扭矩。
由于功率守恒公式可得F*V=T*ω1...........................................(5)
根据公式(2)、(4)、(5)即可求出偏心体1所需要的最大驱动扭矩Tg(N.m)
Tg=Fg*Vgω1............................................................(6)
Fg—公称力(kN)；Vg—滑块达到公称力时对应的速度(mm/s)；
滑块在行程任意位置所允许的载荷是不一样的，允许载荷为F(kN)。
F=Tg*ω1/V.............................................................(7)
V—滑块在行程任意位置时对应的速度(mm/s)；
由于多连杆压力机杆系运动速度较慢，偏心体转速一般为18spm左右，杆系惯性力对分析结果影响很小，杆系摩擦力与惯性力在程序分析中大约合计为滑块公称力的5%；滑块惯性力可由滑块重量和滑块加速度求出。将杆系分解成多个模块，列出平面任意力系的平衡方程：
杆系可分解为以下6个模块，分别进行受力分析。
图3.杆系模块分解
平面任意力系的平衡方程
ΣFxi=0................................................................(8)
ΣFyi=0................................................................(9)
ΣMo(Fi)=0............................................................(10)
根据公式(8)、(9)、(10)求解出各杆受合力............................(11)
对滑块所处的位置进行分析，滑块工作可分为三种区域：(1)拉伸开始至下死点前；(2)下死点前至180°；(3)情况(1)、(2)以外的区域。滑块所处的位置区域不同，受力情况也有所不同，在程序中判断出滑块所处区域后即可利用以下平衡方程式编制程序完成计算以及进行曲线绘制。
F8x-F8y/Tan(φ)=0
F7x+F8x=0
F7y+F8y=0
F7x1+F7x=0
F7y1+F7y=0
F6y1=F6x1*Tan(φ)
F6x=-F6x1
F6y=-F6y1
F5（FFIV）x=-F6x F5（FFIV）y=-F6y F5（FFIV）x1=-F5（FFIV）x F5（FFIV）y1=-F5（FFIV）y
各个连杆之间铰接处采用轴瓦连接，稀油润滑。设计时杆系PV值最大不允许超过许用值，PV值的最大值在程序中利用算法来判定。以一个20000kN多连杆机械压力机为例，程序界面如附图(4)。
图4.程序界面
3 部分程序代码
受力分析
Public Function F_L(ByVal K As Integer,ByValφt As Single,ByValα1As Single,ByValα2As Single)As Single
定义杆系受力函数,α1为拉伸开始处(附图3中St的赋值)对应的角度，α2为下死点前(附图3中Sg的赋值)处对应的角度
Dim L1,L2,L3,L4,L5,L6,L7,L8,α,β,X,Y As Single
……
V=L1*ω1*Cos(φ1)+L7*ω7*Cos(φ7)+L8*ω8*Cos(φ8)
F=Form1.TextBox59.Text*ω1/Abs(V)
A=-L1*ω1^2*Sin(φ1)+L7*(a7*Cos(φ7)-ω7^2*Sin(φ7))+L8*(a8*Cos(φ8)-ω8^2*Sin(φ8))
……
偏心轴处受力值相关分析
If(φt>α1Andφt<=α2)Then F8y=F ElseIf(φt>α2Andφt<180)Then F8y=Form1.TextBox13.Text/Form1.TextBox16.Text Else F8y=-Form1.TextBox13.Text*0.05/Form1.TextBox16.Text
-(Form1.TextBox18.Text/Form1.TextBox16.Text*A)
