Help needed for implementing clutch

Code:

float averageWheelsAngularVelo = 0;

foreach(Wheel w in poweredWheels)

  averageWheelsAngularVelo += w.angularVelocity * drivetrainFraction;

        

ClutchSpeed=averageWheelsAngularVelo* ratio; // wheels average speed *ratio

        

ClutchTorque = clutch.GetTorque(engineAngularVelo, ClutchSpeed,engineTotalTorqueMax);



engineTorque = CalcEngineTorque(rpm) * throttle; 

engineFrictionTorque=CalcEngineFrictionTorque(rpm);



if (clutch.IsLocked()) 

    TransferedTorque=engineTorque;

else

    TransferedTorque=ClutchTorque;



engineTotalTorque = engineTorque - engineFrictionTorque - ClutchTorque;



// update engine angular velo

engineAngularVelo  += (engineTotalTorque) / engineInertia * Time.deltaTime;



if (ratio != 0 && clutch.GetClutchPosition()!=0)   //clutch engaged or engaging

{ 



    // Apply torque to wheels

    foreach(Wheel w in poweredWheels)

    {

        lockingTorque = (averageWheelsAngularVelo - w.angularVelocity) * differentialLockCoefficient;

        w.drivetrainInertia = inertia * drivetrainFraction;

        w.driveFrictionTorque = engineFrictionTorque * Mathf.Abs(ratio) * drivetrainFraction;

        w.driveTorque =  TransferedTorque * ratio * drivetrainFraction + lockingTorque;

        slipRatio += w.slipRatio * drivetrainFraction;

    }

}

Code:

// clutch is modeled as limited highly viscous coupling

public float GetTorque ( float n_engine_speed, float n_clutch_speed, float engine_max_torque)

{

    torque_capacity=1.25f*engine_max_torque;

    

    nengine_max_torque=engine_max_torque;

    engine_speed = n_engine_speed;

    clutch_speed = n_clutch_speed;

    new_speed_diff = engine_speed  - clutch_speed;

    locked = true;

            

    if (clutch_position ==0) {

        locked = false; 

        return 0;

    }    

    else{

        max_torque = clutch_position * torque_capacity;

        friction_torque = max_torque* (new_speed_diff/factor);    // viscous coupling (locked clutch)



        if (friction_torque > max_torque) // slipping clutch

        {

            friction_torque  = max_torque;

            locked = false; 

        }

        else if (friction_torque < -max_torque) // slipping clutch

        {

            friction_torque  = -max_torque;

            locked = false;

        }



        //friction_torque =Mathf.Min(engine_max_torque,friction_torque);

        last_torque = friction_torque;

        return friction_torque;

    }

}