rusefi/unit_tests/tests/lua/test_lua_Leiderman_Khlystov...

118 lines
2.5 KiB
C++

#include "pch.h"
#include "rusefi_lua.h"
#include "lua_lib.h"
// Leiderman-Khlystov Coefficients for Estimating Engine Full Load Characteristics and Performance
TEST(LuaVag, LeidermaKhlystov) {
const char* magic = LUA_POW R"(
currentRpm = 2000
maxPowerHp = 148
maxPowerRpm = 6000
maxTorqueNm = 147
maxTorqueRpm = 3500
maxPowerKw = maxPowerHp * 0.7355
torqAtMaxPower = maxPowerKw * 9549 / maxPowerRpm
print('torqAtMaxPower ' ..torqAtMaxPower)
rpmCoef = maxPowerRpm / maxTorqueRpm
print('rpmCoef ' ..rpmCoef)
torqCoef = maxTorqueNm / torqAtMaxPower
print('torqCoef ' ..torqCoef)
torquePotential =(torqCoef -1) * 100
zz =(100 *(rpmCoef -1) *(rpmCoef -1))
print('torquePotential ' ..torquePotential)
print('zz ' ..zz)
ax = 1 -((torquePotential * rpmCoef *(2 - rpmCoef)) / zz)
bx = 2 *((torquePotential * rpmCoef) / zz)
cx =(torquePotential * rpmCoef * rpmCoef) / zz
print('ax ' ..ax .. ', bx ' ..bx .. ', cx ' ..cx)
rpmRatio = currentRpm / maxPowerRpm
abcMult = ax * rpmRatio + bx * pow(rpmRatio, 2) - cx * pow(rpmRatio, 3)
print ('abcMult ' ..abcMult)
print('rpmRatio ' .. rpmRatio)
LeidermanPower = maxPowerKw * abcMult
print('LeidermanPower ' .. LeidermanPower)
LeidermanTorque = (9550 * LeidermanPower) / currentRpm
print('LeidermanTorque ' .. LeidermanTorque)
function testFunc()
return LeidermanTorque
end
)";
EXPECT_NEAR_M3(testLuaReturnsNumberOrNil(magic).value_or(0), 156.463);
}
TEST(LuaVag, LeidermaKhlystov2) {
const char* magic = LUA_POW R"(
currentRpm = 2000
maxPowerHp = 148
maxPowerRpm = 6000
maxTorqueNm = 147
maxTorqueRpm = 3500
maxPowerKw = maxPowerHp * 0.7355
torqAtMaxPower = maxPowerKw * 9549 / maxPowerRpm
print('torqAtMaxPower ' ..torqAtMaxPower)
rpmCoef = maxPowerRpm / maxTorqueRpm
print('rpmCoef ' ..rpmCoef)
torqCoef = maxTorqueNm / torqAtMaxPower
print('torqCoef ' ..torqCoef)
torquePotential =(torqCoef -1) * 100
print('torquePotential ' ..torquePotential)
ax=2-(25/torquePotential);
bx=(50/torquePotential) - 1;
cx=25/torquePotential;
print('ax ' ..ax .. ', bx ' ..bx .. ', cx ' ..cx)
rpmRatio = currentRpm / maxPowerRpm
abcMult = ax * rpmRatio + bx * pow(rpmRatio, 2) - cx * pow(rpmRatio, 3)
print ('abcMult ' ..abcMult)
print('rpmRatio ' .. rpmRatio)
LeidermanPower = maxPowerKw * abcMult
print('LeidermanPower ' .. LeidermanPower)
LeidermanTorque = (9550 * LeidermanPower) / currentRpm
print('LeidermanTorque ' .. LeidermanTorque)
function testFunc()
return LeidermanTorque
end
)";
EXPECT_NEAR_M3(testLuaReturnsNumberOrNil(magic).value_or(0), 415.858);
}