A lubricating oil temperature simulation method, simulation system, simulation device and computer readable storage medium

Abstract

The application relates to a lubricating oil temperature simulation method, a simulation system, simulation equipment and a computer readable storage medium. The simulation method comprises the following steps: S1, calculating a first influence value of an engine rotor rotating speed on a lubricating oil system temperature and a second influence value of an operating state of an engine fuel system on the lubricating oil system temperature; S2, respectively calculating a steady-state temperature of first lubricating oil and a steady-state temperature of second lubricating oil; S3, calculating an initial value temperature of the first lubricating oil and an initial value temperature of the second lubricating oil; and S4, calculating a temperature of the first lubricating oil and a temperature of the second lubricating oil. The application provides a lubricating oil temperature simulation method, a simulation system, simulation equipment and a computer readable storage medium, and can accurately simulate an aero-engine lubricating oil temperature signal.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of aero-engine design, and particularly relates to a lubricating oil temperature simulation method, a simulation system, a simulation device and a computer readable storage medium. BACKGROUND

[0002] After the controller of an aero-engine is designed, the controller needs to be verified through semi-physical simulation, and whether the control logic and parameters of the engine controller are reasonable is determined according to the semi-physical simulation result.

[0003] During the semi-physical simulation test process, the lubricating oil temperature signal of the aero-engine needs to be simulated and simulated signals are output in real time to the aero-engine controller for collection.

[0004] The lubricating oil temperature signal needs to run in a real-time operating system, and needs to respond to various parameters of the engine in real time to ensure that the lubricating oil temperature can output appropriate lubricating oil temperature signals according to the current running state of the engine.

[0005] In the prior art, for example, Chinese patent application CN 106815389 A provides a temperature simulation device. However, the device is a chip containing multiple blocks suitable for transaction level design, and the temperature is simulated in the form of a hardware circuit. The device cannot adjust various parameters in real time and cannot track the engine state in real time.

[0006] Chinese patent application CN 113591333 A relates to a construction method of a GIS device temperature simulation model based on digital twinning. The patent application establishes a digital twinning model of the GIS device to simulate the temperature of the device. The patent application cannot simulate the engine state in real time and cannot call the board card to output the thermal resistance to the controller.

[0007] Chinese patent application CN 111157120 B relates to a land surface temperature simulation method with spatial continuity. The patent application uses time and space data to calculate the land surface temperature. Although the patent application has a correction algorithm, the patent application is only suitable for thermal infrared remote sensing. SUMMARY

[0008] In view of the above problems of the prior art, the present application provides a lubricating oil temperature simulation method, a simulation system, a simulation device and a computer readable storage medium, which can accurately simulate the lubricating oil temperature signal of the aero-engine.

[0009] Specifically, the present application provides a lubricating oil temperature simulation method suitable for an aero-engine, which comprises the following steps:

[0010] S1, calculating a first influence value of the engine rotor speed on the lubricating oil system temperature and a second influence value of the running state of the engine fuel system on the lubricating oil system temperature;

[0011] S2, the aero-engine comprises a first route lubricating oil and a second route lubricating oil, and the steady state temperature of the first route lubricating oil and the steady state temperature of the second route lubricating oil are calculated based on the first influence value and the second influence value respectively;

[0012] S3, the initial value temperature of the first route lubricating oil is calculated based on the steady state temperature of the first route lubricating oil and the fan inlet temperature, and the initial value temperature of the second route lubricating oil is calculated based on the steady state temperature of the second route lubricating oil and the fan inlet temperature;

[0013] S4, the temperature of the first route lubricating oil is calculated based on the initial value temperature of the first route lubricating oil, and the temperature of the second route lubricating oil is calculated based on the initial value temperature of the second route lubricating oil.

[0014] According to one embodiment of the present application, in step S1, the engine rotor speed comprises a low pressure rotor speed and a high pressure rotor speed, the first contribution rate of the low pressure rotor real-time speed to the lubricating oil system temperature is calculated according to the ratio of the low pressure rotor real-time speed and the low pressure rotor design speed, and the second contribution rate of the high pressure rotor real-time speed to the lubricating oil system temperature is calculated according to the ratio of the high pressure rotor real-time speed and the high pressure rotor design speed;

[0015] The first influence value is calculated based on the first contribution rate and the second contribution rate.

[0016] According to one embodiment of the present application, in step S1, the second influence value of the operating state of the engine fuel system to the lubricating oil system temperature is calculated, comprising steps of:

[0017] The fuel flow, the current fuel temperature and the current lubricating oil temperature are obtained;

[0018] The last period lubricating oil temperature is calculated based on the current lubricating oil temperature;

[0019] The temperature difference value of the current fuel temperature and the last period lubricating oil temperature is calculated, and the product of the temperature difference value and the fuel flow is calculated;

[0020] The second influence value is obtained by looking up the table based on the product of the temperature difference value and the fuel flow.

[0021] According to one embodiment of the present application, in step S2, the steady state temperature of the first route lubricating oil is calculated, comprising steps of:

[0022] Exponential operation is performed based on the first influence value and the temperature rise influence coefficient of the first route lubricating oil;

[0023] The product of the exponential operation result and the steady state value of the first route lubricating oil is obtained;

[0024] The difference value of the product result and the second influence value is calculated as the steady state temperature of the first route lubricating oil;

[0025] calculating the steady state temperature of the second oil circuit, comprising the steps of:

[0026] performing an exponential operation based on the first influence value and the temperature rise influence coefficient of the second oil circuit;

[0027] obtaining the product of the exponential operation result and the steady state value of the second oil circuit;

[0028] calculating the difference between the product result and the second influence value as the steady state temperature of the second oil circuit.

[0029] According to one embodiment of the present application, in step S3, the initial value temperature of the first oil circuit is calculated, comprising the steps of:

[0030] performing a switching value selection based on the steady state temperature of the first oil circuit, the fan inlet temperature and the delay parameter, and the selection result is the initial value temperature of the first oil circuit;

[0031] calculating the initial value temperature of the second oil circuit, comprising the steps of:

[0032] performing a switching value selection based on the steady state temperature of the second oil circuit, the fan inlet temperature and the delay parameter, and the selection result is the initial value temperature of the second oil circuit.

[0033] According to one embodiment of the present application, in step S4, the temperature of the first oil circuit is calculated, comprising the steps of:

[0034] obtaining the sampling time of the first oil circuit and the time constant of the first oil circuit;

[0035] dividing the sampling time of the first oil circuit by the sum of the sampling time of the first oil circuit and the time constant of the first oil circuit to generate a first result, and dividing the time constant of the first oil circuit by the sum of the sampling time of the first oil circuit and the time constant of the first oil circuit to generate a second result;

[0036] calculating the product of the initial value temperature of the first oil circuit and the first result, obtaining the temperature of the first oil circuit in the last period, calculating the product of the temperature of the first oil circuit in the last period and the second result, and adding the two product results to obtain the temperature of the first oil circuit;

[0037] calculating the temperature of the second oil circuit, comprising the steps of:

[0038] obtaining the sampling time of the second oil circuit and the time constant of the second oil circuit;

[0039] divide the sampling time of the second lubricating oil by the sum of the sampling time of the second lubricating oil and the time constant of the second lubricating oil to generate a third result; divide the time constant of the second lubricating oil by the sum of the sampling time of the second lubricating oil and the time constant of the second lubricating oil to generate a fourth result;

[0040] multiply the initial value temperature of the second lubricating oil by the third result to obtain the temperature of the second lubricating oil in the last period, multiply the temperature of the second lubricating oil in the last period by the fourth result, and add the two multiplication results to obtain the temperature of the second lubricating oil.

[0041] The application further provides a lubricating oil system temperature simulation system suitable for the lubricating oil system temperature simulation method, comprising:

[0042] a first calculation unit configured to calculate a first influence value of the engine rotor speed on the lubricating oil system temperature;

[0043] a second calculation unit configured to calculate a second influence value of the operating state of the engine fuel system on the lubricating oil system temperature;

[0044] a third calculation unit configured to calculate the steady-state temperature of the first lubricating oil based on the calculation results of the first calculation unit and the second calculation unit

[0045] a fourth calculation unit configured to calculate the steady-state temperature of the second lubricating oil based on the calculation results of the first calculation unit and the second calculation unit;

[0046] a fifth calculation unit configured to calculate the initial value temperature of the first lubricating oil based on the calculation result of the third calculation unit and the fan inlet temperature;

[0047] a sixth calculation unit configured to calculate the initial value temperature of the second lubricating oil based on the calculation result of the fourth calculation unit and the fan inlet temperature;

[0048] a seventh calculation unit configured to calculate the temperature of the first lubricating oil based on the initial value temperature of the first lubricating oil;

[0049] an eighth calculation unit configured to calculate the temperature of the second lubricating oil based on the initial value temperature of the second lubricating oil.

[0050] According to an embodiment of the application, the first calculation unit comprises:

[0051] a first division module configured to calculate the ratio of the real-time speed of the low-pressure rotor to the design speed of the low-pressure rotor;

[0052] a second division module configured to calculate the ratio of the real-time speed of the high-pressure rotor to the design speed of the high-pressure rotor;

[0053] a first multiplication module configured to calculate a product of an output result of the first division module and a contribution ratio of the low-pressure rotor in the temperature rise of the lubricating oil system;

[0054] a second multiplication module configured to calculate a product of an output result of the second division module and a contribution ratio of the high-pressure rotor in the temperature rise of the lubricating oil system;

[0055] a first addition module configured to calculate a sum of the output result of the first multiplication module and the output result of the second multiplication module.

[0056] According to an embodiment of the present application, the second calculation unit comprises:

[0057] a first delay module configured to calculate a lubricating oil temperature of a previous period based on a current lubricating oil temperature;

[0058] a first subtraction module configured to calculate a temperature difference between the current lubricating oil temperature and an output result of the first delay module;

[0059] a third multiplication module configured to calculate a product of the fuel flow and the output result of the first subtraction module;

[0060] a table lookup module having a data table stored thereon, the data table being configured to represent an influence value of the product of the fuel flow and the temperature difference on the lubricating oil temperature, and the table lookup module being configured to look up the data table based on the output result of the third multiplication module to obtain a corresponding influence value.

[0061] According to an embodiment of the present application, the third calculation unit comprises:

[0062] a first exponential operation module configured to perform an exponential operation on the output result of the first calculation unit and a temperature influence coefficient of the first lubricating oil;

[0063] a fourth multiplication module configured to calculate a product of the output result of the first exponential operation module and a steady-state value of the first lubricating oil;

[0064] a second subtraction module configured to calculate a difference between the output result of the fourth multiplication module and the output result of the second calculation unit.

[0065] According to an embodiment of the present application, the fourth calculation unit comprises:

[0066] a second exponential operation module configured to perform an exponential operation on the output result of the first calculation unit and a temperature influence coefficient of the second lubricating oil;

[0067] a fifth multiplication module configured to calculate a product of the output result of the second exponential operation module and a steady-state value of the second lubricating oil;

[0068] a third subtraction module configured to calculate a difference between an output result of the fifth multiplication module and an output result of the second calculation unit.

[0069] According to an embodiment of the present application, the fifth calculation unit comprises:

[0070] a first delay parameter generation module configured to generate a delay parameter;

[0071] a first switch quantity selection module configured to select a switch quantity based on an output result of the third calculation unit, the fan inlet temperature and an output result of the first delay parameter generation module, and the selected result is the initial value temperature of the first lubricating oil.

[0072] According to an embodiment of the present application, the sixth calculation unit comprises:

[0073] a second delay parameter generation module configured to generate a delay parameter;

[0074] a second switch quantity selection module configured to select a switch quantity based on an output result of the fourth calculation unit, the fan inlet temperature and an output result of the second delay parameter generation module, and the selected result is the initial value temperature of the second lubricating oil.

[0075] According to an embodiment of the present application, the seventh calculation unit comprises:

[0076] a second addition module configured to calculate a sum of the sampling time of the first lubricating oil and the time constant of the first lubricating oil;

[0077] a third division module configured to calculate a division of the sampling time of the first lubricating oil by an output result of the second addition module;

[0078] a fourth division module configured to calculate a division of the time constant of the first lubricating oil by an output result of the second addition module;

[0079] a sixth multiplication module configured to calculate a product of the initial value temperature of the first lubricating oil and an output result of the third division module;

[0080] a second delay module configured to calculate a temperature of the first lubricating oil in a previous period based on the current temperature of the first lubricating oil;

[0081] a seventh multiplication module configured to calculate a product of an output result of the fourth division module and an output result of the second delay module;

[0082] a third addition module configured to calculate a sum of an output result of the sixth multiplication module and an output result of the seventh multiplication module.

[0083] According to an embodiment of the present application, the eighth calculation unit comprises:

[0084] a fourth adding module configured to calculate a sum of the sampling time of the second lubricating oil and the time constant of the second lubricating oil;

[0085] a fifth dividing module configured to calculate a result of the sampling time of the second lubricating oil divided by the output result of the fourth adding module;

[0086] a sixth dividing module configured to calculate a result of the time constant of the second lubricating oil divided by the output result of the fourth adding module;

[0087] an eighth multiplying module configured to calculate a product of the initial value temperature of the second lubricating oil and the output result of the fifth dividing module;

[0088] a third delaying module configured to calculate the temperature of the second lubricating oil in the last period based on the current temperature of the second lubricating oil;

[0089] a ninth multiplying module configured to calculate a product of the output result of the sixth dividing module and the output result of the third delaying module;

[0090] a fifth adding module configured to calculate a sum of the output result of the eighth multiplying module and the output result of the ninth multiplying module.

[0091] The application further provides a lubricating oil temperature simulation device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the lubricating oil temperature simulation method according to any one of the preceding embodiments when executing the computer program.

[0092] The application further provides a computer readable storage medium, wherein a computer program is stored on the computer readable storage medium, and the computer program implements the steps of the lubricating oil temperature simulation method according to any one of the preceding embodiments when executed by a processor.

[0093] The lubricating oil temperature simulation method, simulation system, simulation device and computer readable storage medium provided by the application can accurately simulate the lubricating oil temperature signal of an aero-engine by fully considering the rotor characteristics, fuel characteristics and hysteresis characteristics of the temperature of the aero-engine.

[0094] It should be understood that the above general description and the following detailed description of the application are exemplary and illustrative, and are intended to provide further explanation of the application as claimed. BRIEF DESCRIPTION OF DRAWINGS

[0095] The accompanying drawings are included to provide further explanation of the application, which are incorporated into and constitute a part of this application, and show embodiments of the application, and together with the specification serve to explain the principles of the application.

[0096] In the drawings:

[0097] Figure 1A flow chart of a lubricating oil temperature simulation method according to an embodiment of the present application is shown.

[0098] Figure 2 A structural diagram of a lubricating oil temperature simulation system according to an embodiment of the present application is shown.

[0099] Figure 3 A schematic diagram of a first calculation unit according to an embodiment of the present application is shown.

[0100] Figure 4 A structural diagram of a first calculation unit according to an embodiment of the present application is shown.

[0101] Figure 5 A schematic diagram of a second calculation unit according to an embodiment of the present application is shown.

[0102] Figure 6 A structural diagram of a second calculation unit according to an embodiment of the present application is shown.

[0103] Figure 7 A schematic diagram of a third calculation unit and a fourth calculation unit according to an embodiment of the present application is shown.

[0104] Figure 8 A structural diagram of a third calculation unit according to an embodiment of the present application is shown.

[0105] Figure 9 A structural diagram of a fourth calculation unit according to an embodiment of the present application is shown.

[0106] Figure 10 A schematic diagram of a fifth calculation unit according to an embodiment of the present application is shown.

[0107] Figure 11 A structural diagram of a fifth calculation unit according to an embodiment of the present application is shown.

[0108] Figure 12 A schematic diagram of a sixth calculation unit according to an embodiment of the present application is shown.

[0109] Figure 13 A structural diagram of a sixth calculation unit according to an embodiment of the present application is shown.

[0110] Figure 14 A schematic diagram of a seventh calculation unit according to an embodiment of the present application is shown.

[0111] Figure 15 A structural diagram of a seventh calculation unit according to an embodiment of the present application is shown.

[0112] Figure 16 A schematic diagram of an eighth calculation unit according to an embodiment of the present application is shown.

[0113] Figure 17 FIG. 8 shows a schematic diagram of an eighth computing unit of one embodiment of the present application. DETAILED DESCRIPTION

[0114] It should be noted that the embodiments and features of the embodiments in the present application can be combined with each other without conflict.

[0115] The technical solutions in the embodiments of the present application will be clearly and completely described with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all the embodiments. The following description of at least one exemplary embodiment is actually only illustrative, but not as any limitation on the present application and its application or use. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative labor are within the scope of protection of the present application.

[0116] It should be noted that the terms used herein are only intended to describe specific embodiments, and are not intended to limit the exemplary embodiments according to the present application. As used herein, the singular form is intended to include the plural form unless the context clearly indicates otherwise, and it should also be understood that when the terms "comprise" and / or "include" are used in the specification, there is a feature, step, operation, device, component and / or combination thereof.

[0117] Unless specifically stated otherwise, the relative arrangement of the components and steps illustrated in these embodiments and the numerical expressions and values set forth herein are not limiting of the scope of the present application. It should be understood that the various parts of the drawings are not necessarily drawn to scale, and that, for the purpose of convenience and clarity, not all components and steps can be shown in the drawings. Techniques, methods, and devices known to those of ordinary skill in the relevant art can not be discussed in detail, but can be assumed to be part of the present application. In all examples shown and discussed herein, any specific value should be interpreted as merely an example, and not as a limitation. Thus, other examples of the exemplary embodiments can have different values. It should be noted that like reference numerals and letters refer to like items in the following drawings, and thus, once an item is defined in one drawing, it need not be discussed further in subsequent drawings.

[0118] In the description of the present application, it should be understood that the orientation words such as "front, back, up, down, left, right", "transverse, vertical, perpendicular, horizontal" and "top, bottom" and the like indicated orientation or position relationship are generally based on the orientation or position relationship shown in the drawings, only for the convenience of describing the present application and simplifying the description, without making the opposite statement, these orientation words do not indicate and imply that the indicated device or element must have a particular orientation or be constructed and operated in a particular orientation, therefore, it cannot be understood as a limitation on the scope of protection of the present application; The orientation words "inner, outer" refer to the inner and outer of the contour of each component itself.

[0119] In addition, it should be noted that the use of "first", "second" and the like to limit parts is only for the convenience of distinguishing the corresponding parts, and the above words have no special meaning unless otherwise stated, therefore, it cannot be understood as a limitation on the scope of protection of the present application. In addition, although the terms used in the present application are selected from the commonly known and used terms, some terms mentioned in the specification of the present application may be selected by the applicant according to his or her judgment, and the detailed meaning of each term is described in the relevant part of the description. In addition, the present application is required to be understood not only by the actual terms used, but also by the meaning implied by each term.

[0120] Figure 1 The flow chart of the oil temperature simulation method of one embodiment of the present application is shown. As shown in the figure, an oil temperature simulation method suitable for an aero-engine, comprising the steps of:

[0121] S1, calculating the first influence value of the engine rotor speed on the oil system temperature, calculating the second influence value of the running state of the engine fuel system on the oil system temperature;

[0122] S2, the aero-engine contains first oil and second oil, based on the first influence value and the second influence value, the steady state temperature of the first oil and the steady state temperature of the second oil are calculated respectively;

[0123] S3, based on the steady state temperature of the first oil and the fan inlet temperature, the initial value temperature of the first oil is calculated, based on the steady state temperature of the second oil and the fan inlet temperature, the initial value temperature of the second oil is calculated;

[0124] S4, based on the initial value temperature of the first oil, the temperature of the first oil is calculated, based on the initial value temperature of the second oil, the temperature of the second oil is calculated.

[0125] Preferably, in step S1, the engine rotor speed includes a low-pressure rotor speed and a high-pressure rotor speed, a first contribution rate of the low-pressure rotor real-time speed to the oil system temperature is calculated according to a ratio of the low-pressure rotor real-time speed to the low-pressure rotor design speed, and a second contribution rate of the high-pressure rotor real-time speed to the oil system temperature is calculated according to a ratio of the high-pressure rotor real-time speed to the high-pressure rotor design speed.

[0126] A first influence value is calculated based on the first contribution rate and the second contribution rate.

[0127] Preferably, in step S1, a second influence value of an operating state of an engine fuel system on the oil system temperature is calculated, including steps of:

[0128] An oil flow, a current fuel temperature and a current oil temperature are obtained. The oil flow is a total flow of fuel combusted by the aero-engine, which is calculated by inverse calculation of the fuel metering valve opening in an actual engine or by an instruction value of an engine fuel system model in a real-time simulation system. The fuel temperature can be measured by a fuel temperature sensor in an actual aero-engine or calculated by a fuel temperature model in a real-time simulation system.

[0129] A last-cycle oil temperature is calculated based on the current oil temperature.

[0130] A temperature difference between the current fuel temperature and the last-cycle oil temperature is calculated, and a product of the temperature difference and the oil flow is calculated.

[0131] The second influence value is obtained by looking up a table based on the product of the temperature difference and the oil flow.

[0132] Preferably, in step S2, a steady-state temperature of first oil is calculated, including steps of:

[0133] Exponential operation is performed based on the first influence value and a temperature rise influence coefficient of the first oil.

[0134] A product of an exponential operation result and a steady-state value of the first oil is obtained.

[0135] A difference between the product result and the second influence value is calculated as the steady-state temperature of the first oil.

[0136] A steady-state temperature of second oil is calculated, including steps of:

[0137] Exponential operation is performed based on the first influence value and a temperature rise influence coefficient of the second oil.

[0138] A product of an exponential operation result and a steady-state value of the second oil is obtained.

[0139] The difference between the product result and the second influence value is calculated as the steady state temperature of the second oil path.

[0140] Preferably, in step S3, the initial value temperature of the first oil path is calculated, including steps of:

[0141] The switch value is selected based on the steady state temperature of the first oil path, the fan inlet temperature and a delay parameter, and the selection result is the initial value temperature of the first oil path. The initial value of the delay parameter is 0, and then is always 1.

[0142] The initial value temperature of the second oil path is calculated, including steps of:

[0143] The switch value is selected based on the steady state temperature of the second oil path, the fan inlet temperature and a delay parameter, and the selection result is the initial value temperature of the second oil path.

[0144] Preferably, in step S4, the temperature of the first oil path is calculated, including steps of:

[0145] The sampling time of the first oil path and the time constant of the first oil path are obtained;

[0146] The sampling time of the first oil path is divided by the sum of the sampling time of the first oil path and the time constant of the first oil path to generate a first result, and the time constant of the first oil path is divided by the sum of the sampling time of the first oil path and the time constant of the first oil path to generate a second result;

[0147] The product of the initial value temperature of the first oil path and the first result is calculated, the temperature of the first oil path in the last cycle is obtained, the product of the temperature of the first oil path in the last cycle and the second result is calculated, and the two product results are added to obtain the temperature of the first oil path;

[0148] The temperature of the second oil path is calculated, including steps of:

[0149] The sampling time of the second oil path and the time constant of the second oil path are obtained;

[0150] The sampling time of the second oil path is divided by the sum of the sampling time of the second oil path and the time constant of the second oil path to generate a third result, and the time constant of the second oil path is divided by the sum of the sampling time of the second oil path and the time constant of the second oil path to generate a fourth result;

[0151] The product of the initial value temperature of the second oil path and the third result is calculated, the temperature of the second oil path in the last cycle is obtained, the product of the temperature of the second oil path in the last cycle and the fourth result is calculated, and the two product results are added to obtain the temperature of the second oil path.

[0152] Figure 2A structural schematic diagram of the oil temperature simulation system of one embodiment of the present application is shown. As shown in the figure, the present application also provides an oil system temperature simulation system suitable for the oil system temperature simulation method described above. The oil system temperature simulation system comprises:

[0153] a first calculation unit for calculating a first influence value of the engine rotor speed on the oil system temperature;

[0154] a second calculation unit for calculating a second influence value of the operating state of the engine fuel system on the oil system temperature;

[0155] a third calculation unit for calculating the steady-state temperature of the first route oil based on the calculation results of the first calculation unit and the second calculation unit

[0156] a fourth calculation unit for calculating the steady-state temperature of the second route oil based on the calculation results of the first calculation unit and the second calculation unit;

[0157] a fifth calculation unit for calculating the initial value temperature of the first route oil based on the calculation result of the third calculation unit and the fan inlet temperature;

[0158] a sixth calculation unit for calculating the initial value temperature of the second route oil based on the calculation result of the fourth calculation unit and the fan inlet temperature;

[0159] a seventh calculation unit for calculating the temperature of the first route oil based on the initial value temperature of the first route oil;

[0160] an eighth calculation unit for calculating the temperature of the second route oil based on the initial value temperature of the second route oil.

[0161] Figure 3 A schematic diagram of the first calculation unit of one embodiment of the present application is shown. As shown in the figure, the first calculation unit I-1 has two input signals and one output signal, wherein the input signals are the low-pressure rotor real-time speed I-2 and the high-pressure rotor real-time speed I-3, and the output signal is the first influence value I-4. Figure 4 A structural schematic diagram of the first calculation unit of one embodiment of the present application is shown. Preferably, the first calculation unit comprises:

[0162] a first division module I-1-3 for calculating the ratio of the low-pressure rotor real-time speed I-1-1 and the low-pressure rotor design speed I-1-2.

[0163] A second division module I-1-12 is configured to calculate the ratio of the high-pressure rotor real-time rotating speed I-1-8 and the high-pressure rotor design rotating speed I-1-9. The low-pressure rotor real-time rotating speed I-1-1 and the high-pressure rotor real-time rotating speed I-1-8 are obtained by the high-pressure rotor rotating speed sensor and the low-pressure rotor sensor, and in the real-time simulation system, the engine model is used to calculate the real-time rotating speed.

[0164] A first multiplication module I-1-5 is configured to calculate the product of the output result of the first division module I-1-3 and the contribution ratio I-1-4 of the low-pressure rotor in the oil temperature rise.

[0165] A second multiplication module I-1-11 is configured to calculate the product of the output result of the second division module I-1-12 and the contribution ratio I-1-10 of the high-pressure rotor in the oil temperature rise.

[0166] A first addition module I-1-6 is configured to calculate the sum of the output result of the first multiplication module I-1-5 and the output result of the second multiplication module I-1-11, which corresponds to the first influence value I-1-7, and is equivalent to the output signal I-4 of the first calculation unit I-1 in the formula (1). Figure 3

[0167] Figure 5 A schematic diagram of the second calculation unit in one embodiment of the present application is shown. As shown in the figure, the second calculation unit II-1 has two input signals and one output signal, wherein the input signals are the fuel flow II-2 and the fuel temperature II-3, and the output signal is the second influence value II-4. The fuel flow II-2 is the total fuel flow of the engine participating in combustion, which is calculated by the fuel metering valve opening degree in the actual engine, and is calculated by the instruction value of the engine fuel system model in the real-time simulation system. The fuel temperature II-3 is measured by the fuel temperature sensor in the actual engine, and is calculated by the fuel temperature model in the real-time simulation system. The second influence value II-4 is used to represent the influence of the fuel system on the oil temperature rise.

[0168] Figure 6 A structural schematic diagram of the second calculation unit in one embodiment of the present application is shown. As shown in the figure, preferably, the second calculation unit II-1 includes:

[0169] A first delay module II-1-4 is configured to calculate the last period oil temperature based on the current oil temperature II-1-3.

[0170] A first subtraction module II-1-5 is configured to calculate the temperature difference between the current fuel temperature II-1-2 and the output result of the first delay module II-1-4.

[0171] ​a third multiplication module II-1-6 for calculating the product of the fuel flow II-1-1 and the output result of the first subtraction module II-1-5.

[0172] a lookup table module II-1-7 having a data table stored thereon. The data table is a one-dimensional table, the table header of which is the product of the fuel and the oil temperature difference and the fuel flow, and the data in the table is the influence (value) of the flow product on the oil temperature. The lookup table module II-1-7 looks up the data table based on the output result of the third multiplication module II-1-6 to obtain a corresponding second influence value II-1-8. The second influence value II-1-8 is equivalent to the output signal of the second calculation unit II-1 in Figure 5

[0173] Figure 7 A schematic diagram of the third calculation unit and the fourth calculation unit of one embodiment of the present application is shown. As shown, the input signals of the third calculation unit III-1 are the steady-state input signal III-3 and the second influence value III-4, and the output signal is the first-path oil steady-state temperature III-6. The input signals of the fourth calculation unit III-2 are the steady-state input signal III-3 and the second influence value III-5, and the output signal is the second-path oil steady-state temperature III-7. Figure 8 A structural schematic diagram of the third calculation unit of one embodiment of the present application is shown. Preferably, as shown, the third calculation unit III-1 comprises:

[0174] a first exponential operation module III-1-3 for performing exponential operation on the output result III-1-1 of the first calculation unit I-1 and the temperature influence coefficient III-1-2 of the first-path oil.

[0175] a fourth multiplication module III-1-5 for calculating the product of the output result of the first exponential operation module III-1-3 and the steady-state value III-1-4 of the first-path oil.

[0176] a second subtraction module III-1-7 for calculating the difference between the output result of the fourth multiplication module III-1-5 and the output result III-1-6 of the second calculation unit II-1, and outputting the difference as the output result III-1-8, which is the steady-state temperature of the first-path oil.

[0177] Figure 9 A structural schematic diagram of the fourth calculation unit of one embodiment of the present application is shown. As shown, preferably, the fourth calculation unit comprises:

[0178] a second exponential operation module III-2-3 for performing exponential operation on the output result III-2-1 of the first calculation unit I-1 and the temperature influence coefficient III-2-2 of the second-path oil;

[0179] ​A fifth multiplication module III-2-5 is used to calculate the product of the output result of the second exponential operation module III-2-3 and the steady state value III-2-4 of the second lubricating oil.

[0180] A third subtraction module III-2-7 is used to calculate the difference between the output result of the fifth multiplication module III-2-5 and the output result of the second calculation unit II-1, and output the result III-2-6, which is the steady state temperature of the second lubricating oil.

[0181] Figure 10 A schematic diagram of the fifth calculation unit in one embodiment of the present application is shown. As shown in the diagram, the input signals of the fifth calculation unit IV-1 are the steady state temperature IV-3 of the first lubricating oil and the fan inlet temperature IV-4, and the output signal is the initial value temperature IV-5 of the first lubricating oil. Figure 11 A structural schematic diagram of the fifth calculation unit in one embodiment of the present application is shown. As shown in the diagram, preferably, the fifth calculation unit IV-1 comprises:

[0182] A first delay parameter generation module IV-1-3 is used to generate a delay parameter. The initial value of the first delay parameter generation module IV-1-3 is 0, and then it is always a constant 1.

[0183] A first switch quantity selection module IV-1-4 is used to select a switch quantity based on the output result IV-1-1 of the third calculation unit, the fan inlet temperature IV-1-2 and the output result of the first delay parameter generation module IV-1-3, and the selection result is the initial value temperature IV-1-5 of the first lubricating oil.

[0184] Figure 12 A schematic diagram of the sixth calculation unit in one embodiment of the present application is shown. As shown in the diagram, the input signals of the sixth calculation unit IV-2 are the steady state temperature IV-6 of the second lubricating oil and the fan inlet temperature IV-7, and the output signal is the initial value temperature IV-8 of the second lubricating oil. Figure 13 A structural schematic diagram of the sixth calculation unit in one embodiment of the present application is shown. As shown in the diagram, preferably, the sixth calculation unit IV-2 comprises:

[0185] A second delay parameter generation module IV-2-3 is used to generate a delay parameter. The initial value of the second delay parameter generation module IV-2-3 is 0, and then it is always a constant 1.

[0186] A second switch quantity selection module IV-2-4 is used to select a switch quantity based on the output result IV-2-1 of the fourth calculation unit, the fan inlet temperature IV-2-2 and the output result of the second delay parameter generation module IV-2-3, and the selection result is the initial value temperature IV-2-5 of the second lubricating oil.

[0187] Figure 14 The schematic diagram of the seventh calculation unit of one embodiment of the present application is shown. As shown in the figure, the input signals of the seventh calculation unit V-1 are the initial value temperature V-3 of the first lubricating oil, the sampling time V-4 of the first lubricating oil and the time constant V-5 of the first lubricating oil, and the output signal is the first lubricating oil temperature V-6. Figure 15 The structural schematic diagram of the seventh calculation unit of one embodiment of the present application is shown. As shown in the figure, preferably, the seventh calculation unit V-1 comprises:

[0188] The second addition module V-1-11 is used for calculating the sum of the sampling time V-1-9 of the first lubricating oil and the time constant V-1-10 of the first lubricating oil;

[0189] The third division module V-1-12 is used for calculating the sampling time V-1-9 of the first lubricating oil divided by the output result of the second addition module V-1-11;

[0190] The fourth division module V-1-13 is used for calculating the time constant V-1-10 of the first lubricating oil divided by the output result of the second addition module V-1-11;

[0191] The sixth multiplication module V-1-4 is used for calculating the product of the initial value temperature V-1-1 of the first lubricating oil and the output result V-1-2 of the third division module V-1-12;

[0192] The second delay module V-1-7 is used for calculating the temperature of the first lubricating oil in the last period based on the current temperature of the first lubricating oil;

[0193] The seventh multiplication module V-1-5 is used for calculating the product of the output result V-1-3 of the fourth division module V-1-13 and the output result of the second delay module V-1-7;

[0194] The third addition module V-1-6 is used for calculating the sum of the output result of the sixth multiplication module V-1-4 and the output result of the seventh multiplication module V-1-5, i.e. the first lubricating oil temperature V-1-8.

[0195] Figure 16 The schematic diagram of the eighth calculation unit of one embodiment of the present application is shown. As shown in the figure, the input signals of the eighth calculation unit V-2 are the initial value temperature V-7 of the second lubricating oil, the sampling time V-8 of the second lubricating oil and the time constant V-9 of the second lubricating oil, and the output signal is the second lubricating oil temperature V-10. Figure 17 The structural schematic diagram of the eighth calculation unit of one embodiment of the present application is shown. As shown in the figure, preferably, the eighth calculation unit V-2 comprises:

[0196] a fourth addition module V-2-11 for calculating a sum of the sampling time V-2-9 of the second lubricating oil and the time constant V-2-10 of the second lubricating oil;

[0197] a fifth division module V-2-12 for calculating a result of the sampling time V-2-9 of the second lubricating oil divided by an output result of the fourth addition module V-2-11;

[0198] a sixth division module V-2-13 for calculating a result of the time constant V-2-10 of the second lubricating oil divided by the output result of the fourth addition module V-2-11;

[0199] an eighth multiplication module V-2-4 for calculating a product of the initial value temperature V-2-1 of the second lubricating oil and the output result V-2-2 of the fifth division module V-2-12;

[0200] a third delay module V-2-7 for calculating a temperature of the second lubricating oil in a previous period based on a current temperature of the second lubricating oil;

[0201] a ninth multiplication module V-2-5 for calculating a product of the output result V-2-3 of the sixth division module V-2-13 and an output result of the third delay module V-2-7;

[0202] a fifth addition module V-2-6 for calculating a sum of the output result of the eighth multiplication module V-2-4 and the output result of the ninth multiplication module V-2-5, i.e., the second lubricating oil temperature V-2-8.

[0203] The application further provides a lubricating oil temperature simulation device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of any of the lubricating oil temperature simulation methods when executing the computer program.

[0204] The application further provides a computer readable storage medium, which stores a computer program, wherein the computer program is executable on a processor to implement the steps of any of the lubricating oil temperature simulation methods.

[0205] The specific implementation manners and technical effects of the lubricating oil temperature simulation system, the simulation device and the computer readable storage medium can be referred to the lubricating oil temperature simulation method provided by the application, and will not be described herein.

[0206] The oil temperature simulation method, simulation system, simulation device and computer readable storage medium provided by the application can calculate the real-time value of the oil temperature through the real-time running parameters of the aero-engine, and comprehensively consider the influence of various factors in the real-time simulation parameters, including the rotor characteristics of the aero-engine, the fuel characteristics and the hysteresis characteristics of the temperature. Meanwhile, the simulation method in the application considers the calibration influence of the board card hardware, contains calibration correction parameters, and the result can be directly output to the board card to accurately simulate the output of the oil temperature sensor.

[0207] Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the embodiments disclosed herein can be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans can implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.

[0208] The various illustrative logical blocks, and circuits described in connection with the embodiments disclosed herein can be implemented or performed with a general purpose processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or any combination thereof designed to perform the functions described herein. The general purpose processor can be a microprocessor, but in the alternative, the processor can be any conventional processor, controller, microcontroller, or state machine. A processor can also be implemented as a combination of computing devices, such as a combination of a DSP and a microprocessor, a plurality of microprocessors, one or more microprocessors in conjunction with a DSP core, or any other such configuration.

[0209] The steps of a method or algorithm described in connection with the embodiments disclosed herein can be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such that the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium can be integral to the processor. The processor and the storage medium can reside in an ASIC. The ASIC can reside in a user terminal. In the alternative, the processor and the storage medium can reside as discrete components in a user terminal.

[0210] In one or more exemplary embodiments, the functions described can be implemented in hardware, software, firmware, or any combination thereof. If implemented in software as a computer program product, the functions can be stored on or transmitted over as one or more instructions or code on a computer-readable medium. Computer-readable media includes both computer storage media and communication media including any medium that facilitates transfer of a computer program from one place to another. A storage media can be any available media that can be accessed by a computer. By way of example, and not limitation, such computer-readable media can comprise RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and that can be accessed by a computer. Also, any connection is properly termed a computer-readable medium. For example, if the software is transmitted from a website, server, or other remote source using a coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DSL), or wireless technologies such as infrared, radio, and microwave, then the coaxial cable, fiber optic cable, twisted pair, DSL, or wireless technologies such as infrared, radio, and microwave are included in the definition of medium. Disk and disc, as used herein, includes compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk and blu-ray disc where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. Combinations of the above should also be included within the scope of computer-readable media.

[0211] Various modifications and variations can be made to the above exemplary implementations without departing from the spirit and scope of the application. Therefore, it is intended that the application cover modifications and variations of the application provided they come within the scope of the appended claims and their equivalents.

Claims

1. A lubricating oil temperature simulation method for an aero-engine, comprising the steps of: S1, calculating a first influence value of engine rotor speed on lubricating oil temperature, and calculating a second influence value of engine fuel system operating state on lubricating oil temperature; S2, the aero-engine comprises first lubricating oil and second lubricating oil, and based on the first influence value and the second influence value, the steady-state temperature of the first lubricating oil and the steady-state temperature of the second lubricating oil are calculated respectively; S3, based on the steady-state temperature of the first lubricating oil and the fan inlet temperature, the initial value temperature of the first lubricating oil is calculated, and based on the steady-state temperature of the second lubricating oil and the fan inlet temperature, the initial value temperature of the second lubricating oil is calculated; S4, based on the initial value temperature of the first lubricating oil, the temperature of the first lubricating oil is calculated, and based on the initial value temperature of the second lubricating oil, the temperature of the second lubricating oil is calculated.

2. The lubricating oil temperature emulation method according to claim 1, characterized by, In step S1, the engine rotor speed includes the low-pressure rotor speed and the high-pressure rotor speed, the first contribution rate of the low-pressure rotor real-time speed to the lubricating oil temperature is calculated according to the ratio of the low-pressure rotor real-time speed to the low-pressure rotor design speed, and the second contribution rate of the high-pressure rotor real-time speed to the lubricating oil temperature is calculated according to the ratio of the high-pressure rotor real-time speed to the high-pressure rotor design speed; The first influence value is calculated based on the first contribution rate and the second contribution rate.

3. The lubricating oil temperature emulation method of claim 1, wherein, In step S1, the second influence value of the engine fuel system operating state on the lubricating oil temperature is calculated, comprising the steps of: obtaining the fuel flow, the current fuel temperature and the current lubricating oil temperature; based on the current lubricating oil temperature, the last period lubricating oil temperature is calculated; the temperature difference value of the current fuel temperature and the last period lubricating oil temperature is calculated, and the product of the temperature difference value and the fuel flow is calculated; based on the product of the temperature difference value and the fuel flow, the second influence value is obtained by looking up the table.

4. The lubricating oil temperature simulation method according to claim 3, characterized by, In step S2, the steady-state temperature of the first lubricating oil is calculated, comprising the steps of: exponential operation is performed based on the first influence value and the temperature rise influence coefficient of the first lubricating oil; the product of the exponential operation result and the steady-state value of the first lubricating oil is obtained; the difference between the product result and the second influence value is calculated as the steady-state temperature of the first lubricating oil; the steady-state temperature of the second lubricating oil is calculated, comprising the steps of: exponential operation is performed based on the first influence value and the temperature rise influence coefficient of the second lubricating oil; the product of the exponential operation result and the steady-state value of the second lubricating oil is obtained; the difference between the product result and the second influence value is calculated as the steady-state temperature of the second lubricating oil.

5. The lubricating oil temperature emulation method of claim 1, wherein, In step S3, the initial value temperature of the first lubricating oil is calculated, comprising the steps of: switching value selection is performed based on the steady-state temperature of the first lubricating oil, the fan inlet temperature and the delay parameter, and the selection result is the initial value temperature of the first lubricating oil; the initial value temperature of the second lubricating oil is calculated, comprising the steps of: switching value selection is performed based on the steady-state temperature of the second lubricating oil, the fan inlet temperature and the delay parameter, and the selection result is the initial value temperature of the second lubricating oil.

6. The lubricating oil temperature emulation method of claim 1, wherein, In step S4, the temperature of the first lubricating oil is calculated, comprising the steps of: the sampling time of the first lubricating oil and the time constant of the first lubricating oil are obtained; Divide the sampling time of the first lubricating oil by the sum of the sampling time of the first lubricating oil and the time constant of the first lubricating oil to generate a first result; Divide the time constant of the first lubricating oil by the sum of the sampling time of the first lubricating oil and the time constant of the first lubricating oil to generate a second result; Calculate the product of the initial value temperature of the first lubricating oil and the first result to obtain the temperature of the first lubricating oil in the last cycle, calculate the product of the temperature of the first lubricating oil in the last cycle and the second result, and add the two product results to obtain the temperature of the first lubricating oil. The temperature of the second lubricating oil is calculated, including the steps of: Obtaining the sampling time of the second lubricating oil and the time constant of the second lubricating oil; Divide the sampling time of the second lubricating oil by the sum of the sampling time of the second lubricating oil and the time constant of the second lubricating oil to generate a third result; Divide the time constant of the second lubricating oil by the sum of the sampling time of the second lubricating oil and the time constant of the second lubricating oil to generate a fourth result; Calculate the product of the initial value temperature of the second lubricating oil and the third result to obtain the temperature of the second lubricating oil in the last cycle, calculate the product of the temperature of the second lubricating oil in the last cycle and the fourth result, and add the two product results to obtain the temperature of the second lubricating oil.

7. An oil temperature simulation system adapted to the oil temperature simulation method according to claim 1, characterized in that It includes: A first calculation unit for calculating a first influence value of engine rotor speed on lubricating oil temperature; A second calculation unit for calculating a second influence value of the operating state of the engine fuel system on the lubricating oil temperature; A third calculation unit for calculating the steady-state temperature of the first lubricating oil based on the calculation results of the first calculation unit and the second calculation unit A fourth calculation unit for calculating the steady-state temperature of the second lubricating oil based on the calculation results of the first calculation unit and the second calculation unit; A fifth calculation unit for calculating the initial value temperature of the first lubricating oil based on the calculation result of the third calculation unit and the fan inlet temperature; A sixth calculation unit for calculating the initial value temperature of the second lubricating oil based on the calculation result of the fourth calculation unit and the fan inlet temperature; A seventh calculation unit for calculating the temperature of the first lubricating oil based on the initial value temperature of the first lubricating oil; An eighth calculation unit for calculating the temperature of the second lubricating oil based on the initial value temperature of the second lubricating oil.

8. The oil temperature emulation system of claim 7, wherein, The first calculation unit includes: A first division module for calculating the ratio of the real-time speed of the low-pressure rotor and the design speed of the low-pressure rotor; A second division module for calculating the ratio of the real-time speed of the high-pressure rotor and the design speed of the high-pressure rotor; A first multiplication module for calculating the product of the output result of the first division module and the contribution proportion of the low-pressure rotor in the lubricating oil system temperature rise; A second multiplication module for calculating the product of the output result of the second division module and the contribution proportion of the high-pressure rotor in the lubricating oil system temperature rise; A first addition module for calculating the sum of the output results of the first multiplication module and the second multiplication module.

9. The oil temperature emulation system of claim 8, wherein, The second calculation unit includes: A first delay module for calculating the lubricating oil temperature in the last cycle based on the current lubricating oil temperature; a first subtraction module configured to calculate a temperature difference between the current fuel temperature and an output result of the first delay module; a third multiplication module configured to calculate a product of the fuel flow and the output result of the first subtraction module; a lookup table module having a data table stored thereon, the data table being configured to represent an influence value of the product of the fuel flow and the temperature difference on the oil temperature, and the lookup table module being configured to look up the data table based on the output result of the third multiplication module to obtain a corresponding influence value.

10. The oil temperature emulation system of claim 9, wherein, The third calculation unit comprises: a first exponential operation module configured to perform an exponential operation on the output result of the first calculation unit and a temperature influence coefficient of the first oil; a fourth multiplication module configured to calculate a product of the output result of the first exponential operation module and a steady-state value of the first oil; a second subtraction module configured to calculate a difference between the output result of the fourth multiplication module and the output result of the second calculation unit.

11. The oil temperature emulation system of claim 10, wherein, The fourth calculation unit comprises: a second exponential operation module configured to perform an exponential operation on the output result of the first calculation unit and a temperature influence coefficient of the second oil; a fifth multiplication module configured to calculate a product of the output result of the second exponential operation module and a steady-state value of the second oil; a third subtraction module configured to calculate a difference between the output result of the fifth multiplication module and the output result of the second calculation unit.

12. The oil temperature emulation system of claim 11, wherein, The fifth calculation unit comprises: a first delay parameter generation module configured to generate a delay parameter; a first switch quantity selection module configured to perform switch quantity selection based on the output result of the third calculation unit, the fan inlet temperature, and the output result of the first delay parameter generation module, and the selection result being an initial value temperature of the first oil.

13. The oil temperature emulation system of claim 9, wherein, The sixth calculation unit comprises: a second delay parameter generation module configured to generate a delay parameter; a second switch quantity selection module configured to perform switch quantity selection based on the output result of the fourth calculation unit, the fan inlet temperature, and the output result of the second delay parameter generation module, and the selection result being an initial value temperature of the second oil.

14. The oil temperature emulation system of claim 9, wherein, The seventh calculation unit comprises: a second addition module configured to calculate a sum of a sampling time of the first oil and a time constant of the first oil; a third division module configured to calculate a division of the sampling time of the first oil by the output result of the second addition module; a fourth division module configured to calculate a division of the time constant of the first oil by the output result of the second addition module; a sixth multiplication module configured to calculate a product of the initial value temperature of the first oil and the output result of the third division module; a second delay module configured to calculate a temperature of the first oil of a previous period based on a current temperature of the first oil; a seventh multiplication module configured to calculate a product of the output result of the fourth division module and the output result of the second delay module; a third addition module configured to calculate a sum of the output result of the sixth multiplication module and the output result of the seventh multiplication module.

15. The oil temperature emulation system of claim 9, wherein, The eighth calculation unit comprises: a fourth addition module configured to calculate a sum of a sampling time of the second oil and a time constant of the second oil; a fifth division module configured to calculate a division of the sampling time of the second oil by the output result of the fourth addition module; a sixth division module for calculating a time constant of the second lubricant divided by an output of the fourth addition module; an eighth multiplication module for calculating a product of an initial value temperature of the second lubricant and an output of the fifth division module; a third delay module for calculating a temperature of the second lubricant of a previous period based on a current temperature of the second lubricant; a ninth multiplication module for calculating a product of an output of the sixth division module and an output of the third delay module; a fifth addition module for calculating a sum of an output of the eighth multiplication module and an output of the ninth multiplication module.

16. An engine oil temperature simulation device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the steps of the lubricant temperature simulation method of any one of claims 1-6 when executing the computer program.

17. A computer readable storage medium having stored thereon a computer program, characterized in that, The computer program, when executed by the processor, implements the steps of the lubricant temperature simulation method of any one of claims 1-6.

Images