CN114352404B - Calculation method for parameter design of reverse-triangular rotor engine - Google Patents
Calculation method for parameter design of reverse-triangular rotor engine Download PDFInfo
- Publication number
- CN114352404B CN114352404B CN202210018257.4A CN202210018257A CN114352404B CN 114352404 B CN114352404 B CN 114352404B CN 202210018257 A CN202210018257 A CN 202210018257A CN 114352404 B CN114352404 B CN 114352404B
- Authority
- CN
- China
- Prior art keywords
- rotor
- cylinder
- face
- area
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Supercharger (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
The invention discloses a method for calculating parameter design of an inverted triangle rotor engine, which comprises the steps of calculating an actual molded line of a rotor and an actual molded line of a cylinder body; moving the actual molded line of the rotor according to the rotation angle of the crankshaft as i/2, solving the length of a circumferential working surface line of the rotor and the length of a circumferential contact surface line of the cylinder body, and solving the area of a working volume plane; if i is less than or equal to 1440, i = i +1, and the process is returned; if i is greater than 1440, the single-cylinder working volume, the displacement and the combustion chamber pit volume are obtained according to the cylinder thickness, the maximum swing angle and the eccentricity; and calculating the length of the end face sealing piece, the total length of the sealing piece, the allowable end face area of air inlet and exhaust, the maximum sealing angle, the average sealing piece speed at different rotating speeds and the combustion top dead center face-to-face ratio. The invention provides a theoretical basis for the design and research and development of the reverse triangle rotor engine and shortens the research and development period. The method has the advantages of high efficiency and accuracy in calculation, parallel calculation support, convenience in operation, clean interface and the like.
Description
Technical Field
The invention belongs to the field of simulation design of internal combustion engines, and relates to a method for calculating parameter design of an anti-triangular rotor engine.
Background
With the increasing demand of unmanned aerial vehicles on aviation power, miniaturized aviation power devices with high energy density, strong cruising ability and low cost become one of the important research targets in the field of unmanned aerial vehicles. Unlike the conventional delta Rotor Engine, the reverse delta Rotor Engine is also called reverse delta Rotor Engine (Inverted-Wankel Engine), or cycloidal Rotor Engine (cycloidal Rotor Engine), which is a new Rotor Engine structure type newly proposed in recent years. Compared with the traditional triangle rotor engine, the reverse triangle rotor engine has the characteristics of high compression ratio and small air leakage, can be suitable for spark ignition and compression ignition working modes, and has the advantages of high thermal efficiency, compact structure, high power-weight ratio and high reliability.
The traditional method for designing the triangular rotor engine is mature, relevant theories exist in parameter design, and a triangular rotor engine module is also arranged in numerical simulation software in the working process of the internal combustion engine, so that the application is simple and convenient. However, there is no reverse-triangular rotor engine module in the present numerical simulation software for the working process of all internal combustion engines, and at the same time, the rotor and cylinder profile of the reverse-triangular rotor engine are theoretically complex, the structural parameters of the main parts of the engine are numerous and involve a large amount of theoretical formula calculation, and there are some parameters such as instantaneous arc length, area, working volume, etc. there is no corresponding calculation formula at present, the work of determining these parameters is very complicated, and the working process simulation of the rotor engine, the strength calculation of the engine parts, etc. need to be performed by these parameters, so it is necessary to use computer programming to determine many structural parameters and kinematic parameters of the rotor engine, so as to implement the parameter design of the reverse-triangular rotor engine, complete the performance prediction and optimization of the reverse-triangular rotor engine quickly, and shorten the research and development cycle effectively.
Disclosure of Invention
The invention aims to provide an MATLAB-based design parameter calculation method for an anti-triangle rotor engine, which is efficient, accurate and simple and convenient to operate.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a method for calculating parameter design of an inverse triangle rotor engine comprises the following steps:
1) Calculating the actual molded line of the rotor and the actual molded line of the cylinder body according to the eccentricity, the created radius, the translation distance and the distance of the rotor cylinder body of the reverse triangular rotor engine;
2) Starting from i =0, increasing 1 step length in each cycle, moving the actual molded line of the rotor according to the crank angle of i/2, and intercepting a circumferential working surface line L1 of the single-cylinder rotor, a circumferential contact surface line L2 of the cylinder body and a working volume plane Fa in a two-dimensional plane; calculating the length L1 of a rotor circumferential working surface line and the length L2 of a cylinder circumferential contact surface line by adopting arc length integration, and calculating the area of a working volume plane by adopting a function poly area;
3) If i is not more than 1440, i = i +1, return to step 2); if i is larger than 1440, the single-cylinder working volume, the displacement and the combustion chamber pit volume are obtained according to the cylinder thickness, the maximum swing angle and the eccentricity;
4) According to the rotor circumferential working surface line, the thickness of the cylinder body and the circumferential contact surface line of the cylinder body, the area of the working surface of the single-cylinder rotor and the area of the contact surface of the cylinder body are obtained; solving the area of the front end face and the rear end face by adopting a function area in MATLAB;
wherein L1 is a rotor circumferential working surface line, B is a cylinder thickness, and L2 is a cylinder circumferential contact surface line;
6) Determining the change rule of the rotor actual molded line, the cylinder body actual molded line, the areas of the front end surface and the rear end surface, the area of the cylinder body contact surface, the area of the single-cylinder rotor working surface and the single-cylinder working volume V along with the crank angle;
7) And calculating the length of the end face sealing piece, the total length of the sealing piece, the allowable end face area of air inlet and exhaust, the maximum sealing angle, the average sealing piece speed at different rotating speeds and the combustion top dead center face-to-face ratio.
Further, the actual profile of the rotor is calculated by:
in the formula, v represents X r -Or-Y r Characteristic rotation angles under a coordinate system; x r Is X-axis, or is origin, Y r Is a Y axis; x is a radical of a fluorine atom r As abscissa, y, of a point on the actual profile of the rotor r Is the ordinate of a point on the actual profile of the rotor.
Further, the cylinder actual profile is calculated by the following formula:
in the formula, u represents a characteristic corner under an X-O-Y coordinate system; v represents a characteristic corner under an Xr-Or-Yr coordinate system; x is the abscissa of the point on the cylinder actual-contour line, and y is the ordinate of the point on the cylinder actual-contour line.
Further, the characteristic rotation angle u in the X-O-Y coordinate system is calculated by the following formula (1-3):
further, the single cylinder working volume V and the displacement V h And the volume V of the combustion chamber pit r Calculated by the following formula:
in the formula, phi max Is the maximum swing angle, K is the shape parameter, B is the cylinder thickness,
is a swing angle, V max Is the theoretical maximum working volume, V, of the engine min Is the theoretical minimum working volume of the engine.
Further, the maximum rocking angle is calculated by:
φ max ＝sin -1 (3e/R)。
further, the end-face-seal piece length L3 is calculated by the following procedure: starting from p =1, increasing 1 step per cycle, calculating the end face sealing piece length L (p) in each cycle by using an equation (1-6), wherein the threshold value of p is 1080, and accumulating the end face sealing piece lengths calculated in each cycle after the cycle is finished to obtain an end face sealing piece length L3;
further, the total length L of the sealing piece all Calculated by the following formula (1-7):
L all ＝2L3+3B#(1-7)。
further, the intake-exhaust end face allowable area F1 is calculated by the following formula (1-8):
wherein the maximum sealing angle phi max ＝sin -1 (3e/R)；
Average velocity v of sealing pieces with different rotating speeds m Calculated by the following formula (1-9):
further, combustion top dead center face to face ratio (S/V) max Calculated by the following formula (1-10):
where sr is the combustion chamber area increment, F min The minimum value of the area of the non-dimensional end face ventilation window is obtained.
Compared with the prior art, the invention has the following beneficial effects:
(1) According to the invention, the rotor cylinder body curve is drawn by calculating the horizontal and vertical coordinates of the rotor cylinder body at each rotation angle, and the curves of the front end surface area, the rear end surface area, the cylinder body contact surface area, the single-cylinder rotor working surface area, the single-cylinder working volume and the like are solved by utilizing the poly area function, so that the mathematical description and the equation of each curve are simplified and are realistic. When the length of the end face sealing piece, the total length of the sealing piece, the area of the allowable end face of air intake and exhaust, the maximum sealing angle, the average speed of the sealing pieces with different rotating speeds and the face volume ratio of the combustion top dead center are calculated, the codes are simpler and clearer by setting parameters to replace the expression of a more complex formula, complicated mathematical expressions are simplified, and a series of design parameters are obtained more efficiently. According to the invention, the actual molded lines of the rotor and the cylinder body can be rapidly calculated, the calculation time is greatly shortened, and the design time is saved.
(2) The method can arbitrarily select reasonable basic parameters of the rotary engine to calculate, thereby obtaining actual molded line mathematical description of rotors and cylinders of the rotary engine with different basic parameters and a series of important parameters of the structural design of the rotary engine, and can perform efficient performance simulation calculation only by selecting different calculation modules according to different requirements, thereby having strong functions.
(3) The calculated curve graph and the like can be stored by self, the operation is convenient, and reliable reference and theoretical basis can be provided for the research and development and design of the rotor engine and the performance prediction of the engine.
Drawings
FIG. 1 is a schematic flow chart of a method for calculating parameter design of an inverted triangle rotor engine according to the present invention;
FIG. 2 is an interface schematic diagram of a calculation method for parameter design of an inverted triangle rotor engine according to the present invention.
FIG. 3 is a schematic diagram of an interface after inputting parameters in embodiment 1
FIG. 4 shows the actual profile, front and rear end surface areas F of the rotor cylinder body obtained by calculation h Cylinder contact surface area S rh Single cylinder rotor working face area S r The change rule curve of the single-cylinder working volume V along with the crank shaft angle is shown in the specification, wherein (a) is the actual molded line of the rotor cylinder body, and (b) is the area F of the front end face and the rear end face h And (c) is the cylinder contact surface area S rh (d) is the working surface area S of the single-cylinder rotor r And (e) is a single cylinder working volume V.
FIG. 5 is a graph showing data for rotor machine displacement volume, combustion chamber pocket volume, and end seal length.
FIG. 6 is a schematic diagram of calculating intake and exhaust phase shift distance input parameters.
Fig. 7 is a diagram showing the intake and exhaust phases of the rotor machine.
Detailed Description
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
The design and research and development of the reverse triangular rotor engine provide theoretical basis, and the research and development period is shortened. By applying the method, the parameters of the rotary engine can be calculated through MATLAB, the basic parameters of the rotary engine are input into software, and a calculation program simulates and calculates the parameter curve of the rotary engine without a real object of the rotary engine. The method has the advantages of high efficiency and accuracy in calculation, parallel calculation support, convenience in operation, clean interface and the like.
Referring to fig. 1, a method for calculating parameter design of an inverse triangle rotor engine includes the following steps:
1) According to basic geometrical parameters of the inverted triangle rotor engine: the eccentricity e, the created radius R, the cylinder thickness B, the translation distance a, the rotor cylinder distance delta and the actual compression ratio epsilon are calculated by utilizing a mathematical equation (1-1). The X-O-Y coordinate of the actual profile of the rotor is a gear ring/cylinder component which is fixed and X r -Or-Y r The coordinate is a gear/rotor component, the gear/rotor component rotates clockwise according to the motion characteristic of the planetary gear mechanism, and the origin O is simultaneously r Revolved counterclockwise around the origin O.
In the formula, v represents X r -Or-Y r Characteristic rotation angle/rad under a coordinate system, the value range is [0,6 pi]；X r Is X-axis, or is origin, Y r Is a Y axis; x is a radical of a fluorine atom r As abscissa, y, of a point on the actual profile of the rotor r The ordinate of a point on the actual profile of the rotor.
Calculating the actual molded line of the cylinder body according to the formula (1-2);
in the formula, u represents a characteristic corner/rad under an X-O-Y coordinate system; v represents a characteristic corner/rad under an Xr-Or-Yr coordinate system; x is the abscissa of the point on the cylinder actual-line, Y is the ordinate of the point on the cylinder actual-line, wherein the characteristic rotation angle u in the X-O-Y coordinate system is calculated by the following formula (1-3):
2) Starting from i =0, increasing 1 step length per cycle, moving the actual molded line of the rotor according to the crank angle of i/2, and intercepting a circumferential working surface line L1 of the single-cylinder rotor, a circumferential contact surface line L2 of the cylinder body and a working volume plane Fa in a two-dimensional plane;
3) Calculating the lengths of a rotor circumferential working surface line L1 and a cylinder circumferential contact surface line L2 by adopting arc length integration, and calculating the area of a working volume plane Fa by adopting a function poly area;
4) If i is not more than 1440, i = i +1, return to step 2); if i is more than 1440, the single-cylinder working volume V and the displacement V are obtained according to the formula (1-4) h And the volume V of the combustion chamber pit r ；
In the formula, phi max Is the maximum swing angle/rad, phi max ＝sin -1 (3 e/R). K is a shape parameter, namely the ratio of the radius to the eccentricity K = R/e; v max Is the theoretical maximum working volume, V, of the engine min Is the theoretical minimum working volume of the engine.
5) The working surface area S of the single-cylinder rotor is obtained by the formula (1-5) r Area S of contact surface with cylinder body rh (ii) a Calculating the area F of the front end face and the rear end face by adopting a function poly area in MATLAB h 。
Wherein L1 is a rotor circumferential working surface line, B is a cylinder thickness, and L2 is a cylinder circumferential contact surface line.
6) Obtaining structural parameters and kinematic parameters of the engine, including actual molded lines of the rotor, actual molded lines of the cylinder body, and areas F of the front end face and the rear end face h Area S of cylinder contact surface rh Single cylinder rotor working face area S r The single-cylinder working volume V and other parameters, and the change rule of the structural parameters and the kinematic parameters of the engine along with the crank angle is a curve.
7) Calculating the length L3 of the end face sealing piece and the total length L of the sealing piece all The area F1 of the allowable end face for air intake and exhaust, and the maximum sealing angle phi max Different rotation speeds are sealedAverage speed v of the seals m And combustion top dead center face to face ratio (S/V) max 。
The specific process is as follows: starting from p =1, increasing 1 step per cycle, calculating the end face sealing piece length L (p) in each cycle by using an equation (1-6), wherein the threshold value of p is 1080, and accumulating the end face sealing piece lengths calculated in each cycle after the cycle is finished to obtain an end face sealing piece length L3;
total length L of sealing piece all Obtained from the formula (1-7):
L all ＝2L3+3B#(1-7)
the intake and exhaust end face allowable area F1 is obtained by the formula (1-8):
wherein the maximum sealing angle phi max ＝sin -1 (3e/R)；
Average velocity v of sealing pieces at different rotating speeds m Obtained by the formula (1-9):
combustion top dead center face to face ratio (S/V) max Obtained by the formula (1-10):
where sr is the combustion chamber area increment, F min The minimum value of the area of the air exchange window of the dimensionless end face is obtained.
The method adopts a calculation method for parameter design of the reverse-triangular rotor engine, and basic geometric parameters of the reverse-triangular rotor engine, such as eccentricity e, creation radius R, cylinder thickness B and translation distance are inputa. Selecting different calculation modules according to the requirement for the distance delta of the rotor cylinder body and the actual compression ratio, quickly calculating parameters according to the selected basic geometric parameters and the calculation modules and displaying the parameters in the form of a curve graph, a data text or a dialog box, directly selecting different calculation modules for calculation, and storing the data so as to preliminarily know the structural parameters and the kinematic parameters of the rotor cylinder body, the molded lines of the front end surface and the rear end surface of the rotor cylinder body and the area F of the front end surface and the rear end surface of the rotor cylinder body h Area S of cylinder contact surface rh Rotor working face area S r The change rule of the single-cylinder instantaneous working volume V and the like along with the rotation angle of the crankshaft, the single-cylinder working volume V and the pit volume V of the combustion chamber r End face sealing piece length L3 and sealing piece total length L all The area F1 of the allowable end face for air intake and exhaust, and the maximum sealing angle phi max Average velocity v of sealing pieces at different rotating speeds m And combustion top dead center face to face ratio (S/V) max And adding the calculation of the related parameters on the basis;
the calculation module in the MATLAB APP DESIGNER specifically comprises the following modules, namely a calculation module 1, a curve display module 2, a data display module 3 and a data storage module 4, and the specific description is as follows:
the calculation module 1:
the calculation module 1 is used for calculating performance data of the reverse triangular rotor engine and providing a preliminary theoretical basis for design and research and development of the reverse triangular rotor engine, and the calculation module 1 is used for simulating and calculating the rotor engine according to the creation principle of a rotor profile and a cylinder profile, the planetary motion rule of a rotor and the mathematical relationship among all basic parameters, and calculating performance parameters of the rotor engine, such as the rotor cylinder profile, the front end face area F and the rear end face area F h Cylinder contact surface area S rh Rotor working face area S r The change rule of the single-cylinder instantaneous working volume V and the like along with the rotating angle of the crankshaft, the single-cylinder working volume V and the pit volume V of the combustion chamber r End face sealing piece length L3 and sealing piece total length L all The area F1 of the allowable end face for air intake and exhaust and the maximum sealing angle phi max Average velocity v of sealing pieces at different rotating speeds m And combustion top dead center face to face ratio (S/V) max Etc.;
the curve display module 2:
the curve display of the module 2 is to draw the data calculated by the calculation module 1 in a curve graph form, so that the molded lines of the cylinder body and the rotor and the front and rear end surface areas F of the rotor engine are displayed in a more intuitive mode h Area S of cylinder contact surface rh Rotor working face area S r The change rule of the single-cylinder instantaneous working volume V and the like along with the rotation angle of the crankshaft is a relation curve;
the data display module 3:
the data of the module 3 is displayed by using the data calculated by the calculation module 1, such as a single cylinder working volume V, a combustion chamber pit volume V r End face sealing piece length L3 and sealing piece total length L all The area F1 of the allowable end face for air intake and exhaust, and the maximum sealing angle phi max Average velocity v of sealing pieces at different rotating speeds m And combustion top dead center face to face ratio (S/V) max Etc. are displayed in the information dialog box and the text box for the user to observe and analyze; the data display module 3 displays data as main parameter display, and displays the data by using numerical values or character strings, and the data display module 3 displays the calculated and stored data in the form of texts and numerical values, and can modify parameters for recalculation;
the data saving module 4:
the data storage of the data storage module 4 is to store the data calculated by the calculation module 1 and the curve graph displayed by the curve display module 2 so as to analyze and arrange the post-design data.
Example 1
As shown in FIG. 2, after a calculation module is loaded, the calculation module can be switched in a tab at the top of a window, wherein the tab comprises geometrical data of a rotor engine, solution of a cylinder translation equation, intake and exhaust phase translation distances and program description so as to meet different calculation requirements. After the basic parameters are input, the curve graph is drawn in a newly opened window, and the result parameters are displayed in a dialog box mode.
As shown in fig. 3, basic parameters of the rotary engine, namely eccentricity e =12mm, creation radius R =76mm, cylinder thickness B =50mm, translation distance a =1.45mm, rotor-cylinder distance Δ =0.3mm, actual compression ratio epsilon =11 and rated engine speed n =7000rpm are sequentially input into a numerical box of a geometrical data tab of the rotary engine, and then calculation and drawing are performed. Plotted graphs such as cylinder and rotor profiles, rotor engine front and rear end face areas F h Area S of cylinder contact surface rh Rotor working face area S r The change rule of the single-cylinder instantaneous working volume V and the like along with the crank angle is shown in figure 4.
After obtaining the volume V shown in figure 4, the single-cylinder working volume V and the combustion chamber pit volume V can be carried out r End face sealing piece length L3 and sealing piece total length L all The area F1 of the allowable end face for air intake and exhaust and the maximum sealing angle phi max Average velocity v of sealing pieces at different rotating speeds m And combustion top dead center face to face ratio (S/V) max Etc., as shown in fig. 5.
As shown in fig. 6, the eccentricity e =12mm, the formation radius R =76mm, and the rotor shift distance a are sequentially input in the numerical box in the intake/exhaust phase shift distance tab 0 =1.45mm, IVC (intake air closing phase) =550 ° CA, IVO (intake air opening phase) =300 ° CA, EVO (exhaust air opening phase) =150 ° CA, EVC (exhaust air closing phase) =400 ° CA. Clicking the parameters after the calculation to obtain the result shown in fig. 7.
Other parameters required by strength and reliability calculation can be automatically increased on the basis of the design research later.
(1) According to the method, based on MATLAB APP DESIGNER, the abscissa and ordinate of the rotor cylinder at each rotation angle are calculated through a circulation structure, a curve of the rotor cylinder is drawn, and the front end face area F and the rear end face area F are solved by using for circulation and a poly area function h Area S of cylinder contact surface rh Single cylinder rotor working face area S r And the single-cylinder working volume V and other curves, thereby simplifying and realizing the mathematical description and the equation of each curve. Calculating the length L3 of the end face sealing piece and the total length L of the sealing piece all The maximum allowable intake/exhaust end surface area F1Sealing angle phi max Average velocity v of sealing pieces at different rotating speeds m And combustion top dead center face to face ratio (S/V) max In the process, the code is simpler and clearer by setting parameters to replace the expression of a more complex formula, the complicated mathematical expression is simplified, and a series of design parameters are obtained more efficiently. The user interactivity is good, and the interactive interface is simple and attractive. The software reduces the time of developing and debugging programs for users by using a natural language close to a mathematical expression, a rapid mathematical processing process, a powerful function module and a complete graphic processing function. The parameter calculation of the invention adopts MATLAB code combining parallel calculation and for-loop structure, can quickly calculate the actual molded line of the rotor and the cylinder body, greatly shortens the calculation time and saves the design time.
(2) The method can randomly select reasonable basic parameters of the rotary engine to calculate, thereby obtaining actual molded line mathematical description of rotors and cylinders of the rotary engine with different basic parameters and a series of important parameters of the structural design of the rotary engine, and can perform efficient performance simulation calculation only by selecting different calculation modules according to different requirements, thereby having strong functions.
(3) The calculated curve graph and the like can be stored by self, the operation is convenient, and reliable reference and theoretical basis can be provided for the research and development and design of the rotor engine and the performance prediction of the engine.
The above embodiments are merely illustrative of the technical ideas and features of the present invention, and the purpose thereof is to enable those skilled in the art to understand the contents of the present invention and implement the present invention, and not to limit the protection scope of the present invention. All equivalent changes and modifications made according to the spirit of the present invention should be covered in the protection scope of the present invention.
Claims (8)
1. A method for calculating parameter design of an inverse triangle rotor engine is characterized by comprising the following steps:
1) Calculating an actual rotor profile and an actual cylinder profile according to the eccentricity, the created radius, the translation distance and the distance of the rotor cylinder of the inverted triangle rotor engine;
2) Starting from i =0, increasing 1 step length per cycle, moving the actual molded line of the rotor according to the crank angle of i/2, and intercepting a circumferential working surface line L1 of the single-cylinder rotor, a circumferential contact surface line L2 of the cylinder body and a working volume plane Fa in a two-dimensional plane; calculating the length L1 of a rotor circumferential working surface line and the length L2 of a cylinder circumferential contact surface line by adopting arc length integration, and calculating the area of a working volume plane by adopting a function poly area;
3) If i is not more than 1440, i = i +1, return to step 2); if i is larger than 1440, the single-cylinder working volume, the displacement and the combustion chamber pit volume are obtained according to the cylinder thickness, the maximum swing angle and the eccentricity;
4) According to the rotor circumferential working surface line, the thickness of the cylinder body and the circumferential contact surface line of the cylinder body, the area of the working surface of the single-cylinder rotor and the area of the contact surface of the cylinder body are obtained; calculating the area of the front end face and the rear end face by adopting a function poly area in MATLAB;
in the formula, L1 is a rotor circumferential working surface line, B is a cylinder thickness, and L2 is a cylinder circumferential contact surface line;
6) Determining the change rule of the rotor actual molded line, the cylinder body actual molded line, the areas of the front end surface and the rear end surface, the area of the cylinder body contact surface, the area of the single-cylinder rotor working surface and the single-cylinder working volume V along with the crank angle;
7) Calculating the length of the end face sealing piece, the total length of the sealing piece, the allowable end face area of air inlet and exhaust, the maximum sealing angle, the average sealing piece speed at different rotating speeds and the combustion top dead center face volume ratio;
wherein the actual profile of the rotor is calculated by:
in the formula, v represents X r -Or-Y r Characteristic rotation angles under a coordinate system; x r Is X-axis, or is origin, Y r Is a Y axis; x is a radical of a fluorine atom r To turn toAbscissa, y, of point on sub-actual line r Is the ordinate of the point on the actual profile of the rotor;
the cylinder actual profile is calculated by:
in the formula, u represents a characteristic corner under an X-O-Y coordinate system; v represents a characteristic corner under an Xr-Or-Yr coordinate system; x is the abscissa of the point on the cylinder actual-contour line, and y is the ordinate of the point on the cylinder actual-contour line.
2. The method of claim 1, wherein the characteristic rotation angle u in the XO-Y coordinate system is calculated by the following equation (1-3):
3. the method for calculating the parameter design of the inverse-triangular rotor engine according to claim 1, wherein the single-cylinder working volume V and the displacement V h And the volume V of the combustion chamber pit r Calculated by the following formula:
in the formula, phi max Is the maximum swing angle, K is the shape parameter, B is the cylinder thickness,
is a swing angle, V max Is the theoretical maximum working volume of the engine, V min Is the theoretical minimum working volume of the engine.
4. The method of claim 3, wherein the maximum rocking angle is calculated by the following formula:
φ max ＝sin -1 (3e/R)。
5. the method for calculating the parameter design of the inverted triangle rotor engine according to claim 1, characterized in that the end face seal piece length L3 is calculated by the following process: starting from p =1, increasing 1 step per cycle, calculating the end face sealing piece length L (p) in each cycle by using an equation (1-6), wherein the threshold value of p is 1080, and accumulating the end face sealing piece lengths calculated in each cycle after the cycle is finished to obtain an end face sealing piece length L3;
6. the method of claim 1, wherein the total length L of the sealing plate is greater than the total length L of the sealing plate all Calculated by the following formula (1-7):
L all ＝2L3+3B (1-7)。
7. the method of claim 1, wherein the allowable intake and exhaust end surface area F1 is calculated by the following equation (1-8):
wherein the maximum sealing angle phi max ＝sin -1 (3e/R)；
Average velocity v of sealing pieces with different rotating speeds m Calculated by the following formula (1-9):
8. the method of claim 1, wherein the combustion top dead center face-to-face ratio (S/V) max Calculated by the following formula (1-10):
where sr is the combustion chamber area increment, F min The minimum value of the area of the air exchange window of the dimensionless end face is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210018257.4A CN114352404B (en) | 2022-01-07 | 2022-01-07 | Calculation method for parameter design of reverse-triangular rotor engine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210018257.4A CN114352404B (en) | 2022-01-07 | 2022-01-07 | Calculation method for parameter design of reverse-triangular rotor engine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114352404A CN114352404A (en) | 2022-04-15 |
| CN114352404B true CN114352404B (en) | 2022-12-09 |
Family
ID=81107166
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210018257.4A Active CN114352404B (en) | 2022-01-07 | 2022-01-07 | Calculation method for parameter design of reverse-triangular rotor engine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114352404B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0485465A1 (en) * | 1989-08-11 | 1992-05-20 | Mechanology | Toroidal hyper-expansion rotary engine, compressor, expander, pump and method. |
| CN1073744A (en) * | 1991-12-28 | 1993-06-30 | 张寿龄 | Two-stroke rotary polygonal piston engine |
| CN101158313A (en) * | 2007-11-13 | 2008-04-09 | 上海师范大学 | Circular arc tooth cycloid rotor engine |
| CN101852123A (en) * | 2009-07-24 | 2010-10-06 | 尚世群 | Cycloid rotor engine |
| CN111460677A (en) * | 2020-04-14 | 2020-07-28 | 大连理工大学 | Method for establishing rotor stacking precision prediction model based on geometric algebraic theory |
| CN113836687A (en) * | 2021-05-20 | 2021-12-24 | 北京工业大学 | Design method for combined cylinder profile of triangular rotor engine |
-
2022
- 2022-01-07 CN CN202210018257.4A patent/CN114352404B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0485465A1 (en) * | 1989-08-11 | 1992-05-20 | Mechanology | Toroidal hyper-expansion rotary engine, compressor, expander, pump and method. |
| CN1073744A (en) * | 1991-12-28 | 1993-06-30 | 张寿龄 | Two-stroke rotary polygonal piston engine |
| CN101158313A (en) * | 2007-11-13 | 2008-04-09 | 上海师范大学 | Circular arc tooth cycloid rotor engine |
| CN101852123A (en) * | 2009-07-24 | 2010-10-06 | 尚世群 | Cycloid rotor engine |
| CN111460677A (en) * | 2020-04-14 | 2020-07-28 | 大连理工大学 | Method for establishing rotor stacking precision prediction model based on geometric algebraic theory |
| CN113836687A (en) * | 2021-05-20 | 2021-12-24 | 北京工业大学 | Design method for combined cylinder profile of triangular rotor engine |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114352404A (en) | 2022-04-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fan et al. | Experimental and numerical investigation of the fluid flow in a side-ported rotary engine | |
| Taskiran et al. | Comparison of flow field and combustion in single and double side ported rotary engine | |
| Sadiq et al. | Development of rotary Wankel devices for hybrid automotive applications | |
| Wang et al. | Towards a comprehensive optimization of the intake characteristics for side ported Wankel rotary engines by coupling machine learning with genetic algorithm | |
| Fan et al. | Effects of different parameters on the flow field of peripheral ported rotary engines | |
| Zou et al. | Mathematical modeling and analysis of thermodynamic processes in a twin-rotor piston engine | |
| CN114352404B (en) | Calculation method for parameter design of reverse-triangular rotor engine | |
| Jasak et al. | Rapid CFD simulation of internal combustion engines | |
| Wang et al. | Implementation of a novel dual-layer machine learning structure for predicting the intake characteristics of a side-ported Wankel rotary engine | |
| Hanipah et al. | The real-time interaction model for transient mode investigations of a dual-piston free-piston engine generator | |
| Chang et al. | Numerical investigation of the combined effect of injection angle and injection pressure in a gasoline direct injection rotary engine | |
| Chen et al. | Analysis on structural characteristics of rotors in twin-rotor cylinder-embedded piston engine | |
| Pan et al. | Effects of intake and exhaust valve timing on the performance of an air‐powered rotary engine | |
| Yang et al. | Effect of ignition timing on the combustion process of a port injection free piston linear generator: A system level multi-physics coupling method | |
| Fan et al. | Numerical study on flow field in a peripheral ported rotary engine under the action of apex seal leakage | |
| CN114961985A (en) | Intelligent prediction method and system for performance of hydrogen fuel aviation rotor engine | |
| Fan et al. | The formation mechanisms of flow field in a compound intake ported rotary engine considering apex seal leakage | |
| Zhang et al. | Combustion simulation and key parameter optimization for opposite axial piston engine in small-scale | |
| King | The Development of the Szorenyi Four-Chamber Rotary Engine | |
| Shen | Enhancing the Efficiency of a Rotary Thermal Propulsion System through Advanced Simulation Techniques | |
| Xu et al. | Investigation on theoretic profile curve of the spherical cam transmission mechanism in a Twin-rotor piston engine | |
| CN110134991A (en) | Based on Fluent software to the method for layered-scavenging engine structure optimization | |
| Zhu et al. | Numerical investigation of the effects of intake pipe deflection angles on the in-cylinder flow and vortex structures of a cycloidal rotary engine | |
| Wang et al. | Parametric modeling and analysis of intake phases for side-ported Wankel rotary engines | |
| King | Analysis of the Design Features of the Wankel and Szorenyi Rotary Engines |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |