11606 发表于 2006-11-30 08:37

请高手帮帮我,急

有谁做过这样的仿真吗?

我要根据上面的simulink图得到仿真曲线,我现在不会设置参数。

1.        曲柄滑块机构的机构组成
2.        曲柄滑块机构的向量图如图1所示。
曲柄l1 的长度为l1 (l1 = 30mm ),旋转角位移为θ1 ,旋转角速度为ω1 ,旋转角加速度为α1 ;连杆l2 的长度为l2 (l2= 120mm ),旋转角位移为θ2 , 旋转角速度为ω2 ; 滑块m 的位移量为0,速度为.l0 ,加速度为l0 。
根据曲柄滑块机构的向量图,把曲柄滑块机构的向量按x和y坐标轴方向分解可以得出:
l1 cosθ1 + l2cosθ2 = l0                                             (1)
l1 sinθ1 + l2sinθ2 = 0                                             (2)
式(1) 、(2)对时间求导,并注意到θ1 = ω1 和θ2 = ω2得:
-l1ω1 sinθ1 -l2ω2 sinθ2 = l0                                    (3)
l1ω1 cosθ1 + l2ω2 cosθ2 = 0                                       (4)
写成如下矩阵的形式得:      (5)
这就是曲柄滑块机构的速度方程。

3 曲柄滑块机构的加速度方程

式(3) 、(4)对时间再次求导,并注意到: θ.= 1 ω1 、θ= ω2 和ω1 =α1 、ω2 =α2得:
-l1α1 sinθ1 -l1ω21cosθ1 -l2α2 sinθ2 +-l2ω22cosθ2 =l0               (6)
l1α1 cosθ1 -l1ωl2 sinθ2 221sinθ1 +l2α2 cosθ2 +-l2ω22sinθ2            (7)
写成如下矩阵形式得:
这就是曲柄滑块机构的加速度方程。

2 曲柄滑块机构的匀角速度运动学仿真

211 匀角速度输入时曲柄滑块机构的Simulink 仿真模型

图2 曲柄匀角速运动时Simulink 仿真模型  
slipblock_speed.m文件如下:


曲柄连杆机构曲柄以ω1 = 188. 5 rad/s 作匀角速度运动(这是除草机单缸发动机的正常工作速度),主要研究在这一输入作用下的曲柄连杆机构连杆的运动角速度变化规律以及滑块的位移、速度的变化规律。应用MATLAB 算法建立计算曲柄滑
块机构的速度方程(5)的MATLAB 函数模块slip2 block_speed (该模块in: ω1 、θ1 、θ2 ; out: ω2 、ddl0t) 。应用该MATLAB Function 模块建立了Simulink 仿真模型如图2所示。仿真初始条件: θ1 = 0 rad 、θ2 = 0 rad 、l0 = 150 mm (仿真模型符号对应: alpha1 -α1 ; omega1 -ω1 ; omega_2 -ω2 ; theta1 -θ1 ; theta2 -θ2 ; 以下同) 。

212 匀角速度输入时曲柄滑块机构的Simulink 仿真曲线分析
曲柄连杆机构的曲柄以ω1 = 188. 5 rad/s 作匀角速度运动时(对应图2的仿真模型), 连杆的运动角速度ω2 (omega_2) 的仿真曲线如图3 ( a) 所示。滑块m 的位移仿真曲线如图3 ( b) 所示, 速度仿真曲线如图3 ( c) 所示。

从仿真曲线可以看出,当曲柄以ω1 = 188. 5 rad/s 作匀角速度运动时, 连杆的角速
ω2≈ 48 rad /s(图3 ( a) 所示);滑块m 的位移变化范围是90 mm 2150 mm (图3 ( b) 所示);滑块m 的速度变化范围大约是±0. 6 ×104 mm/s (图3 ( c) 所示) 。而且
它们在0. 07 s内变化了2个周

wf1254 发表于 2006-12-10 18:24

分别双击t积分模块heta1, theta2, l0, 在initial condition栏中输入初始角度和位移就是了,但要注意三者的初始值要满足几何约束关系,三者只有一个独立量,给定其中一个,其它两个就定了,计算出来输入就是了。

不知你问的是否这个意思?
页: [1]
查看完整版本: 请高手帮帮我,急