CN111783328B - Motor joint surface leakage risk determination method based on finite element analysis - Google Patents
Motor joint surface leakage risk determination method based on finite element analysis Download PDFInfo
- Publication number
- CN111783328B CN111783328B CN202010495158.6A CN202010495158A CN111783328B CN 111783328 B CN111783328 B CN 111783328B CN 202010495158 A CN202010495158 A CN 202010495158A CN 111783328 B CN111783328 B CN 111783328B
- Authority
- CN
- China
- Prior art keywords
- motor
- contact stress
- joint surface
- working conditions
- finite element
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention discloses a motor joint surface leakage risk determination method based on finite element analysis, which comprises the following steps: s1: establishing a motor model, and carrying out grid division and load and constraint input on the motor; s2: analyzing and solving the model; s3: searching a contact stress penetration area or a contact stress bandwidth insufficient area under each working condition; s4: and selecting proper working conditions by screening the working conditions. The invention provides a method for determining the leakage risk of a motor joint surface based on finite element analysis, which can find problems early and does not need to carry out a large number of working condition experiments on the motor. The contact stress of the joint surface can be rapidly analyzed, the leakage risk point of the motor joint surface is determined, the joint surface leakage problem is avoided, the initial design problem is eliminated, meanwhile, the test working conditions are screened according to the analyzed contact stress threshold value results under different working conditions, the necessary test working condition quantity is controlled, the research and development time period is controlled, and the research and development cost is saved.
Description
Technical Field
The invention relates to the technical field of motor manufacturing, in particular to a motor joint surface leakage risk determination method based on finite element analysis.
Background
The motor joint surface comprises a motor shell and an acceleration shell joint surface, and the motor shell and a motor end cover joint surface are insufficient in contact stress when the joint surface is insufficient, so that poor sealing can be caused, leakage of motor lubrication cooling liquid can occur, and damage to the motor can be caused. The traditional motor leakage problem is generally verified through tests, and the screening of test conditions is a problem due to the fact that the conditions are relatively large, the leakage problem generally occurs under fatigue endurance conditions, the test period is relatively long, moreover, the test is relatively large in artificial uncontrollable factors, true sources of the problems are difficult to eliminate, and therefore the method is completely dependent on the tests, and is neither economical nor time-efficient. By means of the finite element analysis method, the contact stress distribution of the motor joint surface can be rapidly analyzed, design risk points are found out according to the comparison of the analyzed contact stress result with a set contact stress threshold value and a set contact stress bandwidth, and therefore structural iteration improvement can be rapidly carried out. Meanwhile, analysis results under different working conditions are compared, test working conditions can be effectively selected, unnecessary test working conditions are eliminated, the research and development time period is further controlled, and research and development cost is saved.
The existing method is to make test prototypes according to 3d digital models, test the working conditions of the test papers after assembling, disassemble the test papers after finishing, analyze the distribution of the contact stress on the joint surface according to the distribution of the color concentration of the test papers after the test papers are extruded by the joint surface or show different color concentration changes by matching with a special scanner for data analysis.
The invention provides a mandrel strength optimization design method for a large direct-drive wind turbine generator, which is disclosed in China patent publication No. CN107665278A, publication No. 2018, 02-06, and comprises the following implementation steps: establishing a mandrel three-dimensional solid model for a conical mandrel of a wind turbine generator, importing the mandrel three-dimensional solid model into finite element analysis software for grid division and material attribute definition, establishing a mandrel finite element analysis model, calculating the conical mandrel strength under the condition of limiting load bearing according to the mandrel finite element analysis model, obtaining conical mandrel stress test data under different orthogonal data by changing the structural parameters of the conical mandrel, taking the conical mandrel strength as a constraint condition and the structural parameters of the conical mandrel as optimization variables, and optimizing an objective function C by adopting a genetic algorithm according to the conical mandrel stress test data to obtain the optimal structural parameters. This application can not screen experimental conditions, can not reach the effect of shortening experimental time.
Disclosure of Invention
The invention provides a method for determining the leakage risk of a motor joint surface based on finite element analysis, which can discover problems early and does not need to perform a large number of working condition experiments on a motor, in order to solve the problems of the prior art that the problems cannot be discovered early and the problems of a large number of working condition experiments on the motor need to be performed.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
the technical scheme adopted for solving the technical problems is as follows: a motor joint surface leakage risk determination method based on finite element analysis comprises the following steps:
s1: establishing a motor model, and carrying out grid division and load and constraint input on the motor;
s2: analyzing and solving the model;
s3: searching a contact stress penetration area or a contact stress bandwidth insufficient area under each working condition;
s4: and selecting proper working conditions by screening the working conditions. The contact stress of the joint surface can be rapidly analyzed, the leakage risk point of the motor joint surface is determined according to the contact stress threshold value, the joint surface leakage problem is avoided, the initial design problem is eliminated, the research and development period is shortened, and the research and development cost is reduced. Meanwhile, according to the contact stress threshold value results analyzed under different working conditions, test working conditions are screened, the necessary test working condition quantity is controlled, the research and development time period is further controlled, and the research and development cost is saved. .
Preferably, the step S1 includes the steps of:
s11: selecting a motor shell, a speed reducer shell and a motor end cover as a motor model;
s12: grid division is carried out on a motor shell, a speed reducer shell, a motor end cover and bolts in the motor model;
s13: and inputting the bolt pretightening force load and the bearing inner load into a motor model.
Preferably, the specific process of step S12 is as follows: the motor shell, the reducer shell and the motor end cover are meshed by adopting first-order tetrahedron C3D4 units.
Preferably, the bolts are meshed by adopting a first-order hexahedral C3D8 unit, and the joint part of the root parts of the bolts and the machine shell adopts a common node grid.
Preferably, the step S3 includes the steps of:
s31: setting a cpress value limit threshold and a contact stress bandwidth threshold;
s32: and adjusting the lower limit value of the CPress value until the contact stress penetration region happens to be in the region exceeding the threshold value of the CPress value limit on the joint surface or the contact stress bandwidth exceeds the threshold value of the contact stress bandwidth happens to be in the region exceeding the threshold value of the CPress value limit on the joint surface, and recording the CPress value.
Preferably, the contact stress bandwidth threshold is at least one unit thickness.
Preferably, the step S4 includes the steps of:
s41: ordering the working conditions according to the contact stress threshold value;
s42: and selecting m working conditions with the minimum contact stress threshold as screening test working conditions. .
Therefore, the invention has the following beneficial effects: (1) The contact stress of the joint surface can be rapidly analyzed, the leakage risk point of the motor joint surface is determined according to the contact stress threshold value, the joint surface leakage problem is avoided, the initial design problem is eliminated, the research and development period is shortened, and the research and development cost is reduced. Meanwhile, according to the contact stress threshold value results analyzed under different working conditions, test working conditions are screened, the necessary test working condition quantity is controlled, the research and development time period is further controlled, and the research and development cost is saved.
Drawings
FIG. 1 is a flow chart of the present invention
FIG. 2 is a schematic diagram of a motor according to the present invention
FIG. 3 is a schematic view of a contact stress test according to the present invention
In the figure: 1. the device comprises a suspension bracket, a motor shell, a motor end cover, a speed reducer shell and bolts.
Detailed Description
The invention is further described below with reference to the drawings and detailed description.
Examples: a motor joint leakage risk determination method based on finite element analysis, as shown in FIG. 1, comprises the following steps:
s1: establishing a motor model, and carrying out grid division and load and constraint input on the motor;
s11: selecting a motor shell, a speed reducer shell and a motor end cover as a motor model;
s12: dividing grids of a motor shell 2, a speed reducer shell 4, a motor end cover 3 and bolts 5 in the motor model; as shown in fig. 2, the motor housing, the reducer housing, and the motor end cover are meshed with first order tetrahedron C3D4 units, wherein the size of the units at the joint positions should be less than 2mm and at least 3 units must be distributed in the thickness direction. The screw thread part is removed and simplified when the connecting bolt is modeled, the screw diameter takes the screw pitch diameter value, the bolt adopts a first order hexahedron C3D8 unit, the joint grid is adopted at the joint part of the root of the bolt and the casing, and the size of the bolt grid is within 2mm. In order to enhance analysis convergence, the joint bolt and nut are in contact with at least 3 layers of grids on the joint surface of the shell, and the grid size is smaller than 2mm. The inner ring grid of each shell bearing chamber is established and coupled to the bearing center point in a distributed manner and used as a loading point of bearing load. The suspension bracket adopts a first order tetrahedron C3D4 unit to divide grids, suspension rubber in a connection mounting hole of the suspension bracket and a vehicle body is simplified into a three-way rigidity unit, a rigidity value is input according to a design value, and six-degree-of-freedom constraint grounding is carried out. The suspension bracket is connected with the motor by rigid coupling. The contact pairs are arranged on the matching surfaces, and corresponding friction coefficients are set according to the material properties of the different contact pairs.
S13: and inputting the bolt pretightening force load and the bearing inner load into a motor model. The load includes a bolt preload load and a bearing inner load.
The bolt torque is converted into a pretightening force load through an empirical formula, and the formula is as follows:
f=t/(kd), where F is the bolt preload, k is the equivalent coefficient of friction, typically taking the value 0.2, d is the bolt pitch diameter value. The preload is applied to the selected screw unit plane by abaqus software. The selected unit plane is positioned between the joint surface and the screw rod, and the unit plane at the middle part of the screw rod can be taken.
The motor outputs different torques under different working conditions, and the internal load transmitted to the bearing by the transmission system can be determined according to the motor output torques under different working conditions. The load in the bearing can be extracted through romax gear transmission system analysis software, load decomposition is carried out through romax analysis software according to different working conditions, and forces Fx, fy and Fz in three directions corresponding to the different working conditions are obtained through the load decomposition. The decomposed load force is applied to the bearing coupling point, and note that the direction of decomposed and extracted load is input to abaqus to be consistent in coordinates.
And setting a load linear loading curve and a constraint point position in abaqus.
Multiple analysis steps may be provided to improve convergence. The first step only loads small pretightening force such as 10N, the second step loads all pretightening force, and the third step loads the bearing inner load.
S2: analyzing and solving the model; computing solutions using abaqus nonlinear solver
S3: searching a contact stress penetration area or a contact stress bandwidth insufficient area under each working condition;
s31: setting a cpress value limit threshold and a contact stress bandwidth threshold;
s32: and adjusting the lower limit value of the CPress value until the contact stress penetration region happens to be in the region exceeding the threshold value of the CPress value limit on the joint surface or the contact stress bandwidth exceeds the threshold value of the contact stress bandwidth happens to be in the region exceeding the threshold value of the CPress value limit on the joint surface, and recording the CPress value. The contact stress bandwidth threshold is at least one cell thickness.
And checking a contact stress cpress cloud chart at the Abaqus post-processing module, adjusting the lower limit value of the cpress value from 0 to 0 until a contact stress through region or a contact stress bandwidth insufficient region happens to the joint surface, and considering that the joint surface has joint surface leakage risk under the working condition. The contact stress penetration and the contact stress bandwidth deficiency refer to:
setting two color segments in the abs post-processing cpress cloud image to display, such as red and blue, when the cpress value on the joint surface exceeds a threshold value, displaying red and blue when the cpress value is smaller than the threshold value, as shown in fig. 2, when the cpress cloud image threshold value is adjusted to a, the condition that the thickness of the joint surface is penetrated in the right blue region, which indicates that the lubricating fluid is at risk of leakage from the region, is called contact stress penetration; when the cpress cloud pattern limit value is adjusted to b, just the contact stress bandwidth deficiency occurs, and the contact stress bandwidth deficiency leakage is called. The contact stress threshold c=min (a, b), i.e. the smaller of a and b, is defined.
The contact stress bandwidth refers to the length of a red area along the thickness direction on the joint surface, a contact stress bandwidth threshold value is defined as t (at least one unit thickness), and when the contact stress bandwidth is regulated to a value b, the minimum contact stress bandwidth on the joint surface is just smaller than the contact stress bandwidth threshold value t, so that the joint surface is considered to have a leakage risk.
The optimization process comprises the following steps:
the method is characterized in that a target contact stress threshold value d under the current target working condition is initially defined, when the contact stress threshold value c < d analyzed in the step S32 is designed, the design is not satisfied, at the moment, the rigidity of the area is improved by reducing the bolt hole distances at the two ends of the contact stress bandwidth insufficient area and/or reducing the bolt hole distances at the two ends of the contact stress penetration area, after the bolt hole distances are adjusted, analysis is submitted again until the optimization is carried out until c is more than or equal to d, the effect that c is more than or equal to d can be achieved when the bolt hole distances are reduced to a certain hole distance, and the maximum value of the bolt hole distances at the two ends of the contact stress bandwidth insufficient area and the bolt hole distances at the two ends of the contact stress penetration area can be sought through limited adjustment, so that the number of bolts is saved, the tightness requirement is met under the working condition, and leakage cannot occur.
S4: and selecting proper working conditions by screening the working conditions.
S41: ordering the working conditions according to the contact stress threshold value;
s42: and selecting m working conditions with the minimum contact stress threshold as screening test working conditions.
Assuming that n working conditions are calculated and analyzed, the contact stress threshold value of the working condition 1 is c 1 The contact stress threshold value of the working condition 2 is c 2 The contact stress value of the working condition i is c i Up to a contact stress threshold c for operating mode N n And (3) according to the requirement, supposing that m tests are needed, sorting the stress threshold values from small to large, wherein the smaller the contact stress threshold value is, the worse the working condition is, and then the corresponding working condition of m before the ranking of the front stress threshold value can be screened out for testing, wherein the m working conditions are screening test working conditions.
The invention can rapidly analyze the contact stress of the joint surface, determine the leakage risk point of the motor joint surface according to the contact stress threshold value, avoid the problem of joint surface leakage, eliminate the problem of initial design, shorten the research and development period and reduce the research and development cost. Meanwhile, according to the contact stress threshold value results analyzed under different working conditions, test working conditions are screened, the necessary test working condition quantity is controlled, the research and development time period is further controlled, and the research and development cost is saved.
Claims (5)
1. A motor joint surface leakage risk determination method based on finite element analysis is characterized by comprising the following steps:
s1: establishing a motor model, and carrying out grid division and load and constraint input on the motor;
s2: analyzing and solving the model;
s3: searching a contact stress penetration area or a contact stress bandwidth insufficient area under each working condition;
s4: selecting proper working conditions by screening the working conditions;
step S1 comprises the steps of:
s11: selecting a motor shell, a speed reducer shell and a motor end cover as a motor model;
s12: grid division is carried out on a motor shell, a speed reducer shell, a motor end cover and bolts in the motor model;
s13: inputting the bolt pretightening force load and the bearing inner load into a motor model;
step S3 comprises the steps of:
s31: setting a cpress value limit threshold and a contact stress bandwidth threshold;
s32: and adjusting the lower limit value of the CPress value until the contact stress penetration region happens to be in the region exceeding the threshold value of the CPress value limit on the joint surface or the contact stress bandwidth exceeds the threshold value of the contact stress bandwidth happens to be in the region exceeding the threshold value of the CPress value limit on the joint surface, and recording the CPress value.
2. The method for determining the risk of leakage of the joint surface of the motor based on finite element analysis according to claim 1, wherein the specific process of step S12 is as follows: the motor shell, the reducer shell and the motor end cover are meshed by adopting first-order tetrahedron C3D4 units.
3. The method for determining the leakage risk of the joint surface of the motor based on finite element analysis according to claim 1, wherein the bolts are meshed by using a first-order hexahedral C3D8 unit, and the joint part of the root parts of the bolts and the machine shell adopts a common node mesh.
4. A method of determining risk of leakage at a motor interface based on finite element analysis according to claim 1, wherein the contact stress bandwidth threshold is at least one unit thickness.
5. The method for determining the risk of leakage of the joint surface of the motor based on finite element analysis according to claim 1, wherein the step S4 comprises the steps of:
s41: ordering the working conditions according to the contact stress threshold value;
s42: and selecting m working conditions with the minimum contact stress threshold as screening test working conditions.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010495158.6A CN111783328B (en) | 2020-06-03 | 2020-06-03 | Motor joint surface leakage risk determination method based on finite element analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010495158.6A CN111783328B (en) | 2020-06-03 | 2020-06-03 | Motor joint surface leakage risk determination method based on finite element analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111783328A CN111783328A (en) | 2020-10-16 |
| CN111783328B true CN111783328B (en) | 2023-10-03 |
Family
ID=72753405
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010495158.6A Active CN111783328B (en) | 2020-06-03 | 2020-06-03 | Motor joint surface leakage risk determination method based on finite element analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111783328B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113357349B (en) * | 2021-06-18 | 2022-09-30 | 中国第一汽车股份有限公司 | Prediction method for sealing pressure of joint surface of speed reducer shell |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427658A (en) * | 2019-07-12 | 2019-11-08 | 浙江零跑科技有限公司 | A kind of motor stator system fit tolerance based on finite element analysis determines method |
-
2020
- 2020-06-03 CN CN202010495158.6A patent/CN111783328B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427658A (en) * | 2019-07-12 | 2019-11-08 | 浙江零跑科技有限公司 | A kind of motor stator system fit tolerance based on finite element analysis determines method |
Non-Patent Citations (1)
| Title |
|---|
| 弯曲载荷下特殊螺纹接头密封性能的有限元分析;吴稀勇;闫龙;陈涛;史彬;;钢管(06);第76-79页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111783328A (en) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| GB2545334B (en) | Driveline modeller | |
| CN109977464B (en) | Prediction method of piston cutting deformation based on BP neural network | |
| CN104077428A (en) | Remote finite element analysis method serving for industry alliance | |
| CN111783328B (en) | Motor joint surface leakage risk determination method based on finite element analysis | |
| CN112733294B (en) | Design method of engine cylinder gasket | |
| CN114491863A (en) | Reliability simulation analysis method for threaded connection pair of main bearing seat of engine | |
| CN115422814A (en) | Closed-loop optimization design method for complex electromechanical product driven by digital twin | |
| CN109255141B (en) | Optimization method for cross section shape of forward conceptual design of automobile body | |
| CN116822353A (en) | Digital twin model rapid construction method in manufacturing process | |
| CN116341312A (en) | Split dynamic vibration response result evaluation method for engine | |
| CN112287484B (en) | Complex engineering system reliability design method based on vector agent model | |
| CN112597676B (en) | Automobile suspension vibration source sensitivity identification method and device | |
| CN116595823A (en) | Method, device, equipment and medium for detecting mechanical property of power battery | |
| CN115438543A (en) | System fatigue reliability prediction method for large-scale aviation planetary mechanism | |
| CN107491570B (en) | Cylinder cover gasket sealing detection method based on finite element analysis | |
| CN107665286A (en) | One kind is on automobile engine bearings dynamic analysis method | |
| CN113609742A (en) | Wind generating set main shaft optimization method for overall structure | |
| CN112231858B (en) | Failure prediction method and device for multi-row screw connection structure of composite material cylinder | |
| CN113505435A (en) | Carbon fiber shell analysis method based on response surface model | |
| CN113239476A (en) | Rubber bushing constitutive parameter reverse calibration method based on binomial Mooney-Rivlin model | |
| CN112784379B (en) | Optimal design method of impeller rotor system | |
| CN115048747A (en) | Method for optimizing engine crankshaft structure | |
| CN117217037A (en) | Method for analyzing strength and fatigue safety coefficient of engine connecting rod | |
| CN117852249A (en) | Crankshaft fatigue life prediction method considering rotational speed influence | |
| CN114528724A (en) | Complete engine dynamics split type modeling method |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 310051 1st and 6th floors, no.451 Internet of things street, Binjiang District, Hangzhou City, Zhejiang Province Applicant after: Zhejiang Zero run Technology Co.,Ltd. Address before: 310051 1st and 6th floors, no.451 Internet of things street, Binjiang District, Hangzhou City, Zhejiang Province Applicant before: ZHEJIANG LEAPMOTOR TECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |