Removing nonlinear equations from CVVD mechanism

In the previous post An Ideal Continuously Variable Valve Duration, an ideal continuously variable valve duration model was created. This model was useful for understanding how a CVVD model works however it was slow because it introduced nonlinear systems of equations. This post looks at modifying this model to remove these nonlinear systems of equations.

What is a nonlinear system equation and why do they typically make a model slow?

In this case the nonlinear systems being referred to are algebraic loops or coupled equations that utilize an iterative loop to solve. For further details please refer to Introduction to model translation and symbolic processing. To solve these equations the inverse of the Jacobian of the nonlinear system of equations has to be calculated and an iterative process is used, both of which can be computationally demanding.

Also having nonlinear systems of equations in a model can reduce model robustness.

Why does the CVVD model contain nonlinear systems?

To work out what variables are present in the nonlinear systems in the Simulation->Setup on the Translation tab select List non-linear iteration variables as in Figure 1.

Figure 1.  'List non-linear iteration variables' is selected in Simulation setup to gain further information about the nonlinear equations
Figure 1. ‘List non-linear iteration variables’ is selected in Simulation setup to gain further information about the nonlinear equations

When the model is translated with List non-linear iteration variables selected then the translation log displays which variables are in the nonlinear systems of equations as in Figure 2.

Figure 2.  CVVD model with translation log displaying the nonlinear iteration variables
Figure 2. CVVD model with translation log displaying the nonlinear iteration variables

There are two nonlinear systems in the CVVD model in Figure 2. One of the nonlinear systems contains the cam shaft slider position and CVVD disc angle and the other contains the cam lobe pivot angle and the cam lobe slider position. The iteration variables in each nonlinear system are dependent on each other and so are solved in an iterative process.

How can the nonlinear systems be removed?

In this case, it seems likely that there is a function for the cam lobe angle that is dependent on the angle of the timing axle and the relative location of the CVVD disc. Finding this relationship would simplify the system and hopefully will remove all nonlinear systems.

Determining the cam lobe angle function

By using geometry it is possible to determine the cam lobe function based on the timing axle and the relative location of the CVVD disc. In Figure 3 the CVVD disc angle is calculated.

Figure 3.  The calculation of the CVVD disc angle
Figure 3. The calculation of the CVVD disc angle

It is assumed that when the timing axle angle is 0 that this aligns with the direction of the displacement s as in Figure 3. An angular offset could be included if this were not the case.

Figure 4.  The calculation of the cam lobe angle
Figure 4. The calculation of the cam lobe angle

The angle of the cam lobe pivot is calculated from the equation in Figure 4 by using the atan2 function in Modelica.

The equations to solve for the cam lobe angle unfortunately make use of atan2. This will return an angle that is wrapped between -pi and pi; however, we require the cam lobe pivot angle to be continuous. Fortunately, in the VeSyMA library there is an external function called VeSyMA.Roads.Functions.External.continuousS which will return the continuous solution from the wrapped cam lobe angle. The equations used to calculate the cam lobe angle are below:

Figure 5.  Modelica equations to calculate the cam lobe angle
Figure 5. Modelica equations to calculate the cam lobe angle

Calculation of the cam lobe torque

In the previous section the angle for the cam lobe pivot is calculated; however, the torque is also required. This CVVD mechanism can be seen as an ideal variable gear, so the equations from the Modelica.Mechanics.Rotational.Components.IdealGear were used to calculate the torque. The gearing ratio is calculated from the angular speed of the timing axle divided by the angular speed of the cam lobe. The equations for calculating the torque are in Figure 6.

Figure 6.  Equations to calculate the cam lobe torque
Figure 6. Equations to calculate the cam lobe torque

Comparison between previous and current model

A test rig was created to run the lobes with the CVVD mechanism in place as in Figure 7.

Figure 7.  Test rig used to test CVVD mechanism
Figure 7. Test rig used to test CVVD mechanism

The new CVVD mechanism ran about 30% faster and has no nonlinear systems. The model is about 5% slower than the fixed cam lobe model.

Written by: Garron Fish – Chief Engineer

Please get in touch if you have any questions or have got a topic in mind that you would like us to write about. You can submit your questions / topics via: Tech Blog Questions / Topic Suggestion

CONTACT US

Got a question? Just fill in this form and send it to us and we'll get back to you shortly.

Sending

© Copyright 2010-2023 Claytex Services Ltd All Rights Reserved

Log in with your credentials

Forgot your details?