CN113076671B - Damper simplification method in finite element simulation analysis and computer equipment - Google Patents
Damper simplification method in finite element simulation analysis and computer equipment Download PDFInfo
- Publication number
- CN113076671B CN113076671B CN202110348791.7A CN202110348791A CN113076671B CN 113076671 B CN113076671 B CN 113076671B CN 202110348791 A CN202110348791 A CN 202110348791A CN 113076671 B CN113076671 B CN 113076671B
- Authority
- CN
- China
- Prior art keywords
- oil film
- finite element
- damper
- guide rail
- spring
- 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
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Vibration Prevention Devices (AREA)
Abstract
The invention relates to the technical field of finite element simulation analysis, and discloses a damper simplification method and computer equipment in finite element simulation analysis. Establishing a three-dimensional solid model of a guide rail and a damper; carrying out mesh division on the three-dimensional entity model to obtain a finite element model; calculating oil film dynamic parameters on a contact surface between the guide rail and the damper; adding a spring damping unit in the finite element model according to the oil film dynamics parameter; and carrying out mechanical simulation analysis on the finite element model added with the spring damping unit. The spring damping unit is adopted to replace an oil film, only structural force simulation can be carried out on the guide rail and the damper, complex calculation combining fluid simulation and structural force simulation is not needed, and therefore the finite element model establishing process of the damper for the guide rail pair in finite element simulation analysis is simplified. Meanwhile, the calculation amount of finite element analysis is simplified, the analysis efficiency is increased, and the calculation precision is ensured not to be lost.
Description
Technical Field
The invention relates to the technical field of finite element simulation analysis, in particular to a damper simplifying method and computer equipment in finite element simulation analysis.
Background
The precise sliding table transmission device formed by the rolling linear guide rail pair is a key transmission part of a high-end numerical control machine tool, has the characteristics of high rigidity, high bearing capacity, high speed, high-precision positioning and the like, but the unique defect of poor vibration resistance compared with the surface-to-surface contact of the traditional sliding guide rail is that the vibration resistance is poor because the joint part formed by the sliding block rolling body and the guide rail is in point or line contact, and the dynamic characteristic of the numerical control machine tool is poor directly caused. The contact rigidity of the slider-guide rail joint part is usually a weak link in the rigidity constitution of the whole machine, thereby becoming a great factor influencing the precision of the whole machine. Meanwhile, the contact damping characteristic of the joint part of the sliding block and the guide rail is poor, and the dynamic characteristic of the whole machine is severely restricted. The static and dynamic performance of the guide rail system will greatly affect the performance of the complete machine of a high-end numerical control machine tool, and therefore, the contact damping characteristic of the slider-guide rail joint part needs to be simulated and analyzed. When the traditional simulation is carried out on the damping vibration attenuation platform of the precision rolling linear guide rail pair system, the finite element simulation of the structural mechanics can be realized only by carrying out complex calculations such as fluid simulation and the like. How to simplify the damper for the precise guide rail pair in finite element simulation analysis and not influence the calculation precision is very important.
Disclosure of Invention
Based on this, it is necessary to provide a method and a computer device for simplifying the damper for the precision guide rail pair in the finite element simulation analysis, aiming at the problem of how to simplify the damper for the precision guide rail pair in the finite element simulation analysis without affecting the calculation accuracy.
A damper simplifying method in finite element simulation analysis comprises the steps of establishing a three-dimensional solid model of a guide rail and a damper; carrying out mesh division on the three-dimensional entity model to obtain a finite element model; calculating oil film dynamic parameters on a contact surface between the guide rail and the damper; adding a spring damping unit in the finite element model according to the oil film dynamics parameter; and carrying out mechanical simulation analysis on the finite element model added with the spring damping unit.
In the damper simplification method in the finite element simulation analysis, the three-dimensional entity models of the guide rail and the damper are drawn in the three-dimensional drawing software, and the model is introduced into the finite element pretreatment software and is subjected to meshing so as to obtain the finite element model. The outer surface of the guide rail and the inner surface of the damper are connected by adopting a spring damping unit. And the parameters of the spring damping unit are determined according to the dynamic parameters of an oil film on the contact surface between the guide rail and the damper. And importing the finite element model added with the spring damping unit into engineering simulation software for mechanical calculation. Since the oil film is a fluid, complicated calculations such as fluid simulation are required when the oil film is simulated. The spring damping unit is adopted to replace an oil film, only structural force simulation can be carried out on the guide rail and the damper, and fluid-solid simulation is not needed, so that simplification in finite element analysis of the damper is realized. Meanwhile, the calculation amount of finite element analysis is simplified, the analysis efficiency is increased, and the calculation precision is ensured not to be lost.
In one embodiment, the oil film dynamics parameters include oil film normal stiffness and oil film normal damping.
In one embodiment, the calculating of the oil film dynamic parameter on the contact surface between the guide rail and the damper includes obtaining a reynolds equation of the contact surface according to the length and the width of the contact surface between the guide rail and the damper; setting the boundary condition of the Reynolds equation and calculating the fluid pressure on the contact surface; integrating the area of the fluid pressure on the contact surface to obtain the oil film bearing capacity; and calculating the oil film normal stiffness and the oil film normal damping according to the oil film bearing capacity.
In one embodiment, a plurality of contact surfaces with different length and width dimensions are arranged between the guide rail and the damper, and the oil film dynamic parameter is calculated according to the length and width dimensions of each contact surface.
In one embodiment, the spring damping unit comprises a spring, and the elastic parameter of the spring is determined according to the oil film dynamic parameter.
In one embodiment, the adding a spring damping unit in the finite element model according to the oil film dynamics parameter includes determining a radial stiffness of the spring according to the oil film normal stiffness, determining a radial damping of the spring according to the oil film normal damping, and determining a length of the spring according to a gap between the guide rail and the damper; and connecting the guide rail and each contact surface of the damper by using the spring damping unit.
In one embodiment, the three-dimensional solid model is divided by using hexahedral mesh.
In one embodiment, the three-dimensional solid model is subjected to meshing, and when the finite element model is obtained, the Jacobian of all the elements in the finite element model is ensured to be not less than 0.7.
In one embodiment, the finite element model is imported into CAE finite element simulation software for statics and dynamics analysis.
A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor when executing the computer program performs the steps of the damper reduction method in finite element simulation analysis as described in any one of the above embodiments.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the specification, and other drawings can be obtained by those skilled in the art without inventive labor.
FIG. 1 is a schematic method flow diagram of a damper simplification method in a finite element simulation analysis according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a simulation of a finite element model according to an embodiment of the present invention;
FIG. 3 is a schematic cross-sectional view of a three-dimensional solid model according to one embodiment of the present invention;
FIG. 4 is a schematic flow chart of a method for calculating oil film dynamics parameters according to an embodiment of the present invention;
FIG. 5 is a schematic view showing the length and width dimensions of the contact surface between the guide rail and the damper according to one embodiment of the present invention;
FIG. 6 is a schematic flow chart of a method for adding a spring damping unit to a three-dimensional solid model according to an embodiment of the present invention;
FIG. 7 is a cross-sectional view of a finite element model with spring damping elements added according to one embodiment of the present invention.
Detailed Description
To facilitate an understanding of the invention, the invention will now be described more fully with reference to the accompanying drawings. Preferred embodiments of the present invention are shown in the drawings. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
The high-end equipment manufacturing industry in China faces the embarrassing situation that key equipment and parts are seriously imported, the static and dynamic characteristics of the high-end numerical control machine tools made in China are poor, the core technology of basic parts is limited by people for a long time, and long-term and healthy development of the high-end numerical control machine tools in China is always restricted. Because the matching capability of the basic components is relatively insufficient, the performance level of the core components can not meet the requirements, and the research, development and manufacturing level and capability of the whole machine are further influenced. The contact rigidity of the slider-guide rail joint part is usually a weak link in the whole machine rigidity constitution, thereby becoming a big factor influencing the whole machine precision, and meanwhile, the contact damping characteristic of the slider-guide rail joint part is poor, thereby seriously restricting the dynamic characteristic of the whole machine. Therefore, the static and dynamic performance of the guide rail system will greatly affect the performance of the complete machine of the high-end numerical control machine.
A method of dissipating the energy of the object vibration as much as possible in the damping layer when the object vibrates is called damping vibration. Damping refers to an action that opposes the relative motion of an object and converts the energy of the motion into heat or other energy that can be dissipated. The damping reduces the resonance amplitude of the mechanical structure, thereby avoiding the structural damage caused by the dynamic stress reaching the limit. The damping can reduce the resonance amplitude of the mechanical structure, thereby avoiding the structural damage caused by the dynamic stress reaching the limit. Meanwhile, the processing precision, the measuring precision and the working precision of various machine tools, instruments and the like can be improved. Various machines, especially precision machine tools, need higher shock resistance and dynamic stability when working in a dynamic environment, and the dynamic performance of the machines can be greatly improved through various damping treatments. Therefore, the vibration resistance of the system can be effectively improved by increasing the viscous damping of the guide rail system by using the damping element, so that the static and dynamic properties of the guide rail system are improved.
Finite element simulation refers to the simulation of a real physical system (geometric and load working conditions) by using a mathematical approximation method. With simple and interacting elements (i.e., cells), a finite number of unknowns can be used to approximate a real system of infinite unknowns. The damping performance of the damper in precision equipment such as a high-grade main control machine tool can be verified by performing simulation on the damping vibration attenuation platform of the precision rolling linear guide rail pair system by a finite element simulation method. The method has great theoretical and practical significance for vibration control, dynamic design and optimal design of high-end numerical control machines and improvement of the overall technical level of the high-end intelligent equipment industry in China.
Fig. 1 is a flow chart of a method for simplifying a damper in a finite element simulation analysis according to an embodiment of the present invention, wherein the method includes the following steps S100 to S500.
Step S100: and establishing a three-dimensional solid model of the guide rail and the damper.
Step S200: and carrying out mesh division on the three-dimensional solid model to obtain a finite element model.
Step S300: and calculating the dynamic parameters of an oil film on the contact surface between the guide rail and the damper.
Step S400: and adding a spring damping unit in the finite element model according to the oil film dynamics parameter.
Step S500: and carrying out mechanical simulation analysis on the finite element model added with the spring damping unit.
In step S100, a three-dimensional CAD solid model of the guide rail and the damper is created using three-dimensional CAD software. The three-dimensional CAD software can be Solidworks/Proe/Creo/NX UG/CATIA/Inventor software and the like. After the three-dimensional solid models of the damper and the guide rail are built, the exported file format is an intermediate format such as xt/step/stp/iges, and the file format is convenient to call the three-dimensional solid model file in finite element preprocessing software.
When performing the finite element analysis on the three-dimensional solid model in step S200, the model needs to be first gridded in order to make the model become a finite element. The mesh division is to divide the model into a plurality of small units, after the mesh division, the displacement increment of the unit node is the basic unknown quantity in the finite element iteration process, and the mesh quality after the structure dispersion directly influences the solving time and the correctness of the solving result. The three-dimensional solid model file is imported into finite element preprocessing software, and high-precision meshing is performed on the three-dimensional solid model file, so that the finite element model shown in fig. 2 can be obtained. FIG. 2 is a schematic diagram of a finite element model according to an embodiment of the present invention.
The oil film damping vibration attenuation mechanism is that vibration energy is consumed by extruding an oil film, the vibration amplitude is reduced, and lubricating oil in the oil film is extruded by surface force generated by vibration when the oil film is extruded. An oil film is formed in a gap between the damper and the guide rail, and the damper extrudes the oil film to absorb vibration energy of the rotor, so that system vibration is reduced, and system running stability is improved. In step 300, by reasonably selecting relevant parameters of the oil film damping structure according to the surface condition of the contact surface between the guide rail and the damper, different oil film damping values can be obtained to adapt to the optimal damping requirements of different guide rail systems, and the guide rail system has remarkable vibration absorption and vibration resistance effects on different guide rail systems.
However, oil films are fluid rather than solid, and therefore fluid simulation is required to simulate oil films. The guide rail and the damper need to be subjected to solid simulation, so that complex simulation calculation combining fluid simulation and structural force simulation is needed when a system is simulated. The optimal damping requirement in the guide rail system is determined by calculating oil film dynamic parameters on a contact surface between the guide rail and the damper, and a spring damping unit which can be used for simulating the damping action of an oil film between the guide rail and the damper is selected according to the oil film dynamic parameters. Therefore, in step S400 of the present embodiment, a spring damping unit is added to the finite element model according to the oil film dynamic parameter, so as to use the spring damping unit instead of the oil film. The damping effect of the oil film in the system is simulated by utilizing the spring damping unit, the structural force simulation can be carried out on the system only, and the complex calculation combining the fluid simulation and the structural force simulation is not needed, so that the model establishment process of the guide rail auxiliary damper in finite element simulation analysis is simplified.
In step 500, a finite element model with spring damping elements added is derived from the finite element preprocessing software, and then the finite element model is imported into the engineering simulation software for mechanical calculation. The spring damping unit is adopted to replace an oil film, only structural force simulation can be carried out on the guide rail and the damper, and fluid-solid simulation is not needed, so that simplification in finite element analysis of the damper is realized. Meanwhile, the calculation amount of finite element analysis is simplified, the analysis efficiency is increased, and the calculation precision is ensured not to be lost.
In one embodiment, the oil film dynamics parameter includes oil film normal stiffness K n And oil film normal damping C n . Oil film kinetic parameters are main data reflecting damping performance of the oil film kinetic parameters and are closely related to various parameters of the oil film. In the embodiment, the oil film normal stiffness and the oil film normal damping are mainly used for representing the oil film damping performance. FIG. 3 is a cross-sectional view of a three-dimensional solid model according to an embodiment of the invention. Calculating oil film according to contact parameters of contact surfaces of damper and guide railNormal stiffness and oil film normal damping. The contact parameters comprise the fluid property of an oil film and the geometric dimension of a contact surface between the guide rail and the damper.
Fig. 4 is a flowchart illustrating a method of calculating an oil film dynamic parameter according to an embodiment of the present invention, where in one embodiment, calculating the oil film dynamic parameter on the contact surface between the guide rail and the damper includes the following steps S310 to S340.
Step S310: and acquiring a Reynolds equation of the contact surface according to the length and the width of the contact surface between the guide rail and the damper and the fluid property of the oil film.
Step S320: setting the boundary condition of Reynolds equation and calculating the fluid pressure on the contact surface.
Step S330: the fluid pressure is integrated over the contact surface area to obtain the oil film loading.
Step S340: and calculating the normal stiffness and the normal damping of the oil film according to the oil film bearing capacity.
FIG. 5 is a schematic view showing the length and width of the contact surface between the guide rail and the damper according to one embodiment of the present invention, wherein the contact surface between the guide rail and the damper has a length L x Width L of y A rectangular flat plate of (1). In step S310, the length L of the rectangular flat plate is acquired x Width L of y And obtaining the Reynolds equation of the oil film in contact with the rectangular flat plate on the contact surface. The expression of the Reynolds equation is:
wherein h is the film thickness of the oil film; eta is the dynamic viscosity of the oil film fluid; p is the squeeze film pressure of the oil film fluid; t is time.
In step S320, a boundary condition is set for the reynolds equation for squeezing the film by the rectangular plate, and the oil film pressure distribution on the rectangular plate is calculated. In this example, the boundary conditions set by the Reynolds equation for a rectangular plate extruded film are:
by solving the above equation by using a separation variable method, the following expression of the fluid pressure distribution of the oil film can be obtained:
assuming that the damper moves relative to the guide rail in a simple harmonic mode, the vibration velocity upsilon satisfies the following relation:
wherein upsilon is 0 For velocity vibration amplitude, ω is the vibration frequency.
In step S330, the fluid pressure is integrated over the entire rectangular plate to obtain the following load capacity generated by the extrusion fluid pressure:
in step S340, according to the definition of the stiffness, the oil film normal stiffness is calculated, and the fluid stiffness of the normal squeeze film is obtained as follows:
according to the definition of the damping, the oil film normal damping is calculated, and the normal squeeze oil film fluid damping is obtained as follows:
according to the method for simplifying the damper in finite element simulation analysis, the spring is used for simulating surface contact when an oil film is contained between the damper and the guide rail, the contact proportion is represented by the Reynolds equation and the contact factor, software simulation can be performed on the vibration state of the system through a simulation means, fluid simulation is not needed, and the simulation efficiency is greatly improved.
In one embodiment, a plurality of contact surfaces with different length and width dimensions are arranged between the guide rail and the damper, and oil film dynamic parameters are respectively calculated according to the length and width dimensions of each contact surface. As can be seen from the cross section shown in fig. 3, the damping slider inner surface and the guide rail outer surface together include 9 contact surfaces, and the lengths and widths of the rectangular flat plates on different contact surfaces are different, so that different oil film dynamic parameters need to be calculated for the contact surfaces with different lengths and widths respectively to meet the optimal damping requirements required by different contact surfaces between the guide rail and the damper, and the damping slider has significant vibration absorption and vibration resistance effects for different contact surfaces.
In one embodiment, the oil film between the damper and the guide rail is replaced by a Spring damping unit comprising a Spring (Spring), wherein the elastic parameter of the Spring is determined according to the oil film dynamic parameter. In the three-dimensional solid model, the damping action of a spring on a contact surface between a guide rail and a damper is utilized to simulate an oil film damping vibration attenuation mechanism between the damper and the guide rail. The spring with the same damping performance as the oil film is adopted to replace the oil film, only structural force simulation can be carried out on the guide rail and the damper, and fluid-solid simulation is not needed, so that simplification in finite element analysis of the damper is realized. Meanwhile, the calculation amount of finite element analysis is simplified, the analysis efficiency is increased, and the calculation precision is ensured not to be lost.
Fig. 6 is a flowchart illustrating a method for adding a spring damping unit in a three-dimensional solid model according to an embodiment of the present invention, wherein in an embodiment, adding a spring damping unit in a finite element model according to the oil film dynamic parameter includes the following steps S410 to S420.
Step S410: and determining the radial stiffness of the spring according to the normal stiffness of the oil film, determining the radial damping of the spring according to the normal damping of the oil film, and determining the length of the spring according to the gap between the guide rail and the damper.
Step S420: the spring damping unit is used to connect the guide rail with each contact surface of the damper.
In this embodiment, the value of the radial stiffness of the spring is the same as the value of the normal stiffness of the oil film, the value of the radial damping of the spring is the same as the value of the normal damping of the oil film, and the length of the spring is the same as the gap between the guide rail and the damper. As shown in fig. 3, there are 9 contact surfaces between the inner surface of the damping slider and the outer surface of the guide rail, and the oil film normal stiffness K of each of the 9 contact surfaces needs to be calculated n And oil film normal damping C n Value, oil film normal stiffness K calculated on different contact surfaces in step S410 n And oil film normal damping C n A spring on the contact surface is selected. For example, when selecting the spring parameter of the spring added between the contact surface of the upper surface of the damper and the guide rail, the radial stiffness of the spring is selected as the oil film normal stiffness K calculated according to the length and width of the contact surface of the upper surface of the damper and the guide rail n The radial damping of the spring is selected as the oil film normal damping C calculated according to the length and width of the contact surface of the upper surface of the damper and the guide rail n While making the initial length of the spring and the gap between the upper surface of the damper and the guide rail uniform.
In step 420, a spring is added to each contact surface between the inner surface of the damping slider and the outer surface of the guide rail, and 9 springs are used to connect each contact surface of the guide rail and the damper, respectively, so as to obtain a cross section of the finite element model shown in fig. 7. FIG. 7 is a cross-sectional view of a finite element model with spring damping elements added according to one embodiment of the present invention. The spring damping unit is adopted to replace an oil film, so that the simulation efficiency is improved, and meanwhile, the high-efficiency simulation of the damping vibration attenuation platform of the precision rolling linear guide rail pair system can be realized.
It should be understood that although the steps in the flowcharts of fig. 1, 4, and 6 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 1, 4, and 6 may include multiple steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed alternately or alternately with other steps or at least some of the other steps or stages.
In one embodiment, the three-dimensional solid model is divided by using hexahedral mesh. The finite element meshing is a crucial step for carrying out finite element numerical simulation analysis, and directly influences the accuracy of a subsequent numerical calculation analysis result. Meshing involves the shape of the cell and its topology type, the cell type, the choice of mesh generator, the density of the mesh, the number of cells, and the geometric voxels, etc. Because the stiffness matrixes of different units are different and the solving modes of numerical integration are different, in practical application, reasonable units are needed to be adopted for simulation solving.
Meanwhile, due to limited computing resources, the number of grid grids needs to be reasonably controlled, the number is excessive, the computing result is very large, and more time is consumed, so that when static and dynamic computing is performed, grids above 4 layers are often divided in the thickness direction, and the grids around holes are densely divided, so that the influence of stress concentration can be eliminated (the stress concentration is caused by a plurality of conditions, the constraint is simplified during simulation, and the stress concentration in a local area is caused), the grids above 4 layers are also favorable for checking whether the stress penetrates or not, and the follow-up further analysis is facilitated.
In this embodiment, hypermesh is selected as the finite element preprocessing software to perform meshing on the three-dimensional solid model, and the meshing operation performed by using hypermesh is simple and convenient. Meanwhile, the hexahedral mesh is used for dividing the three-dimensional solid model, and the method has the advantages of high precision, few meshes, high calculation efficiency and the like.
In one embodiment, the three-dimensional solid model is subjected to meshing, and when the finite element model is obtained, the Jacobian of all the elements in the finite element model is ensured to be not less than 0.7. The Jacobian of the grid refers to the ratio of the smallest Jacobian to the largest Jacobian in the Jacobian matrix values at the integration points within the cell. For example, there is only one integration point in the triangle unit, with a Jacobian of 1. When the grid division is checked, the Jacobian is commonly used for checking the degree of the shape of the unit, 1 is an ideal shape, the Jacobian is required to be more than 0.7 in general structure calculation, and the calculation result is credible. If the value of the Jacobian is too low, the calculation result is influenced, and the calculation result is not converged. Therefore, in this embodiment, it is necessary to ensure that the Jacobian of each element in the finite element model is not less than 0.7, and prevent the Jacobian of the element from being too low to affect the calculation result.
In one embodiment, the finite element model is imported into CAE finite element simulation software for statics and dynamics analysis. In this embodiment, the CAE finite element simulation software may be ANSYS or other simulation software.
A computer device comprising a memory and a processor, said memory storing a computer program, wherein said processor when executing said computer program performs the steps of the damper reduction method in a finite element simulation analysis as described in any of the preceding embodiments.
A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the damper simplification method in a finite element simulation analysis as set forth in any one of the above embodiments.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database or other medium used in the embodiments provided herein can include at least one of non-volatile and volatile memory. Non-volatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical storage, or the like. Volatile Memory can include Random Access Memory (RAM) or external cache Memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.
In the description herein, references to the description of "some embodiments," "other embodiments," "desired embodiments," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, schematic depictions of the above terms do not necessarily refer to the same embodiment or example.
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (9)
1. A damper simplification method in finite element simulation analysis is characterized by comprising the following steps:
establishing a three-dimensional solid model of the guide rail and the damper;
carrying out mesh division on the three-dimensional entity model to obtain a finite element model;
calculating oil film dynamic parameters on a contact surface between the guide rail and the damper; the oil film dynamic parameters comprise oil film normal stiffness and oil film normal damping;
adding a spring damping unit in the finite element model according to the oil film dynamics parameter;
performing mechanical simulation analysis on the finite element model added with the spring damping unit;
wherein the calculating of the oil film dynamics parameter on the contact surface between the guide rail and the damper comprises the following steps:
acquiring a Reynolds equation of a contact surface between the guide rail and the damper according to the length and the width of the contact surface and the fluid property of an oil film;
setting the boundary condition of the Reynolds equation and calculating the fluid pressure on the contact surface;
integrating the area of the fluid pressure on the contact surface to obtain the oil film bearing capacity;
and calculating the oil film normal stiffness and the oil film normal damping according to the oil film bearing capacity.
2. The method of claim 1, wherein said meshing said three-dimensional solid model to obtain a finite element model comprises:
making the thickness of the three-dimensional solid model at least 4 layers of grids;
and carrying out secondary meshing on the peripheral region of the hole in the three-dimensional solid model to obtain the finite element model.
3. The method of claim 1, wherein the guide rail and the damper include a plurality of contact surfaces with different length and width dimensions, and the oil film dynamics parameter is calculated according to the length and width dimensions of each contact surface.
4. The method according to claim 1, characterized in that the spring damping unit comprises a spring, the spring parameters of which are determined from the oil film dynamics parameters.
5. The method of claim 4, wherein the adding a spring damping unit in the finite element model according to the oil film dynamics parameter comprises:
determining the radial stiffness of the spring according to the oil film normal stiffness, determining the radial damping of the spring according to the oil film normal damping, and determining the length of the spring according to the gap between the guide rail and the damper;
and connecting the guide rail and each contact surface of the damper by using the spring damping unit.
6. The method of claim 1, wherein the three-dimensional solid model is partitioned using a hexahedral mesh.
7. The method of claim 1, wherein the three-dimensional solid model is gridded to ensure that the Jacobian of all elements in the finite element model is not less than 0.7 when the finite element model is obtained.
8. The method of claim 1, wherein the finite element model is imported into CAE finite element simulation software for statics and dynamics analysis.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor realizes the steps of the method of any one of claims 1 to 8 when executing the computer program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110348791.7A CN113076671B (en) | 2021-03-31 | 2021-03-31 | Damper simplification method in finite element simulation analysis and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110348791.7A CN113076671B (en) | 2021-03-31 | 2021-03-31 | Damper simplification method in finite element simulation analysis and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113076671A CN113076671A (en) | 2021-07-06 |
| CN113076671B true CN113076671B (en) | 2023-04-07 |
Family
ID=76614133
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110348791.7A Active CN113076671B (en) | 2021-03-31 | 2021-03-31 | Damper simplification method in finite element simulation analysis and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113076671B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109027017A (en) * | 2018-08-15 | 2018-12-18 | 重庆交通大学 | A kind of Space Rolling Bearing state of wear appraisal procedure |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5421655A (en) * | 1987-05-29 | 1995-06-06 | Kmc, Inc. | Fluid dampened support having variable stiffness and damping |
| CN104111051B (en) * | 2014-07-08 | 2016-08-17 | 华中科技大学 | A kind of measure the device of oil film thickness between rolling linear guide and antivibrator |
| CN104143026B (en) * | 2014-07-31 | 2017-07-14 | 北京工业大学 | It is a kind of to consider to tilt the static pressure oil pad Calculation Methods for Performance with heat |
-
2021
- 2021-03-31 CN CN202110348791.7A patent/CN113076671B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109027017A (en) * | 2018-08-15 | 2018-12-18 | 重庆交通大学 | A kind of Space Rolling Bearing state of wear appraisal procedure |
Non-Patent Citations (1)
| Title |
|---|
| 基于动网格模型的液体动静压轴承刚度阻尼计算方法;熊万里等;《机械工程学报》;20121205(第23期);第118-126页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113076671A (en) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108984829B (en) | Calculation method and system for stacking process of rock-fill concrete rock-fill body | |
| CN110210079B (en) | Machine tool supporting part quality matching method oriented to complete machine dynamic characteristics | |
| CN110941923A (en) | Method for determining sensitive parameters of air spring structure | |
| CN102902859A (en) | Hybrid sliding bearing fluid-solid coupling design method based on phase change model | |
| CN116757146B (en) | Distributed random walk parasitic capacitance extraction method, device, equipment and medium | |
| CN113177290A (en) | Satellite component temperature field prediction method based on depth proxy model normalization | |
| CN113076671B (en) | Damper simplification method in finite element simulation analysis and computer equipment | |
| CN101241520A (en) | Model state creation method based on characteristic suppression in finite element modeling | |
| CN111881578B (en) | Graph database based mechanical product digital twin model layered modeling method | |
| CN117034851A (en) | Region decomposition method based on pruning and graph segmentation optimization | |
| CN107480364B (en) | BIM (building information modeling) model dynamic organization method for engineering construction information management | |
| CN112989524B (en) | Ball bearing rolling body modeling method and system and storage medium | |
| CN112115616B (en) | Dynamic characteristic analysis method and device of power transmission tower | |
| CN113806951A (en) | Elastic simulation method for natural adjacent point search based on half-edge data structure | |
| CN103324803A (en) | Modeling method and device of spinning machine | |
| CN116227295B (en) | Battery pack finite element modeling method, device, computer equipment and storage medium | |
| US8271518B2 (en) | Nearest neighbor search method | |
| CN113027921B (en) | Method and device for acquiring pressure distribution of oil film of static and dynamic pressure bearing | |
| CN114048570A (en) | Vibration isolation efficiency analysis method and device of vibration isolator, terminal and storage medium | |
| CN113282976B (en) | Powder bed construction method based on COMSOL | |
| CN110750940B (en) | Fuzzy evaluation method for comprehensive performance of rolling linear guide rail | |
| CN117454520A (en) | Thickness distortion scaling model construction method and system for train thin-wall energy-absorbing pipe | |
| CN116663366A (en) | Sliding valve core radial force calculation method, system, equipment and medium | |
| Li et al. | Researches on the Fast Determination of Cell Type of Cartesian Grid | |
| CN117421963A (en) | Combined optimization method and device for rotary machine |
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 |