用PDS(Powertrain Dynamic Simulator)快速建立V6发动机Crank Train动力学模型,然后导入Virtual.Lab Motion进行仿真计算,并和理论计算结果进行对比。
结论:
- 1、曲轴系统惯性不平衡力对发动机Z向产生6阶力冲击,X向和Y向无力冲击。
- 2、曲轴系统惯性不平衡力对发动机X向产生3阶力矩冲击(Roll),Y向(Pitch)和Z向(Yaw)产生一致的1阶和2阶力矩冲击。
- 3、仿真计算的结果和理论计算结果一致,但是理论计算一般只考虑2阶以内,而仿真计算可到几千HZ。
仿真计算过程:
1、PDS建模,设置发动机转速6000(1阶频率=100HZ),运行2个循环(曲轴转4圈)。
2、设置平衡块参数,将Reciprocating Mass Coefficient=0,这是为了后面和理论计算做比较。
The Reciprocating Mass Coefficient affects the magnitude of mass used in the counterweights that are attached to the crankshaft in the PDS model. The formulas used to define the counterweight (or "bob weight") mass are typically based on an idealization which breaks the conrod mass properties into a so-called "big end" & "little end", with the "big end" mass identified as a "rotating mass", and the "little end" identified as a "reciprocating mass". Nearly all hand-derived solutions will define the counterweight mass as a combination of these two masses. However many differ from each other regarding how much additional mass to lump into the reciprocating component in order to account for the additional mass of the piston, wrist pin, ring pack, etc.
It is possible to find analytically derived imbalance equations that are tabulated according to engine type in a wide variety of engine design and mechanical design reference texts. One such derivation can be found in Taylor1. Taylorís formulation operates under the assumption that the Counterweight masses have no contribution from the reciprocating mass.With respect to the equations below, the assumption is that q=0.
q = Reciprocating Mass Coefficient
Although it may be more typical to run models with q=1.0 or 0.5, we will reset the Reciprocating Mass Coefficient inside PDS to 0.0 for later correlation against theoretical results derived using Taylorís equation.
3、PDS生成建模文件.def,导入Motion
4、用后处理工具绘制发动机壳体(block)受力及FFT图,发动机受Z向(Vertical)6阶冲击力,X向(Fore-Aft)和Y向(Lateral)不受冲击力。
5、用后处理工具绘制发动机壳体(block)受力矩及FFT图,发动机受X向(Roll)3阶冲击力矩,Y向(Pitch)和Z向(Faw)受到一致的1阶和2阶冲击力矩。
6、与理论计算对比
Now let‘s compare our results against known theoretical solution found in Taylor, The Internal Combustion Engine in Theory and Practice, Vol. 2. The key parameters used for both PDS and the hand calculations are summarized in Table 1 below.
Parameter | Value | Description |
m | 5.1E-04 | Big End ConRod Mass (Mg) |
m1 | 1.0E-04 | Little End ConRod Mass (Mg) |
m2 | 4.56E-04 | Piston Mass (Mg) |
L | 151.75 | ConRod Length (mm) |
R | 43.0 | Crank Throw (mm) |
a | 98.0 | Bore Spacing (mm) |
b | 20.0 | Bore Offset (mm) |
Poff | 0.0 | Wrist Pin Offset (mm) |
N | 6000.0 | Crank RPM (real) |
Table 1: Key Input Parameters
60∞ V6 |
1st order Force |
1st order Couple |
2nd order Force |
2nd order Couple | |
Taylor Analytical Formula | 0 | 1.5 Q a Circular CW | 0 | 1.5 R/L Q a Circular CCW
| |
Analytical Reference Value | 0 | 1,387,459 N-mm | 0 | 393,151 N-mm | |
PDS Simulation | 0 | 1,387,153 N-mm Pitch 1,387,153 N-mm Yaw | 0 | 397,820 N-mm Pitch 397,820 N-mm Yaw |
Table 2: Theory vs. PDS
Regarding Table 2, M = m1+m2, and Q = N^2MR, where N = 6000RPM = 628.3185 rad/sec. Furthermore, a "circular" moment simply implies the pitch and yaw couples are meant to be equal in magnitude (as they are in this case), whereas an "elliptical" moment would mean they are unequal.
Thus, it can be seen that PDS generates accurate imbalance loads. It can be noted that most theory tables such as those presented in Table 2 present results only for 1st and 2nd order, assuming that higher order effects are negligible. PDS models do not make such assumptions, as the Virtual.Lab Motion solver works in the time domain in such a way as to track frequencies well into the kilohertz range.
Furthermore, PDS can extend much further than most hand-derived formulas, because one can examine the influence of many other parameters such as variable bore spacing, different bank-angles producing ësplit piní crank arrangements, piston pin offsets, combustion loads with misfire events, and even the inclusion of a Valve Train subsystem to incorporate the effects of valve line imbalances.
7、References