CN116680842A - Rolling type confluence ring structural parameter simulation optimizing method, equipment and medium - Google Patents
Rolling type confluence ring structural parameter simulation optimizing method, equipment and medium Download PDFInfo
- Publication number
- CN116680842A CN116680842A CN202310976666.XA CN202310976666A CN116680842A CN 116680842 A CN116680842 A CN 116680842A CN 202310976666 A CN202310976666 A CN 202310976666A CN 116680842 A CN116680842 A CN 116680842A
- Authority
- CN
- China
- Prior art keywords
- ring
- simulation
- elastic ring
- rolling type
- elastic
- 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.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 89
- 238000005096 rolling process Methods 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 claims abstract description 58
- 230000006835 compression Effects 0.000 claims abstract description 10
- 238000007906 compression Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims description 24
- 238000003860 storage Methods 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000013515 script Methods 0.000 description 29
- 238000012360 testing method Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 101100243399 Caenorhabditis elegans pept-2 gene Proteins 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 101150000679 OPTN gene Proteins 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 101100243401 Caenorhabditis elegans pept-3 gene Proteins 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01R—ELECTRICALLY-CONDUCTIVE CONNECTIONS; STRUCTURAL ASSOCIATIONS OF A PLURALITY OF MUTUALLY-INSULATED ELECTRICAL CONNECTING ELEMENTS; COUPLING DEVICES; CURRENT COLLECTORS
- H01R39/00—Rotary current collectors, distributors or interrupters
- H01R39/64—Devices for uninterrupted current collection
-
- 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
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to the technical field of simulation, in particular to a method, equipment and medium for optimizing structural parameters of a rolling type confluence ring in a simulation manner, wherein the method comprises the following steps: constructing a simulation model of a target rolling type confluence ring; wherein the fitting offset of the elastic ring in the simulation model is set equal to the compression amount of the elastic ring; simulating the simulation model by adopting a plurality of analysis steps to obtain a corresponding simulation result; the multi-analysis step comprises an initial analysis step, a pressure loading analysis step, a rebound analysis step and a rotation analysis step; and acquiring the target structural parameters of the target rolling type confluence ring based on the simulation result. The invention can improve the simulation efficiency of the rolling type confluence ring.
Description
Technical Field
The invention relates to the technical field of simulation, in particular to a method, equipment and medium for optimizing structural parameters of a rolling type confluence ring in a simulation mode.
Background
The bus ring is a key component of radar equipment, is also an important device for transmitting electric signals of a fixed part and a rotating part, has the main functions of realizing transmission of electric power or electric signals between the fixed equipment and the rotating equipment, avoiding mutual winding of wires when the equipment rotates relatively, and can be widely applied to various fields such as aerospace, aviation, navigation, various precision instruments and the like. The bus ring may be divided into a sliding friction type bus ring and a rolling friction type bus ring according to the difference of contact friction modes. With the continuous development of the bus ring technology, the rolling bus ring is widely applied to civil products and munitions. The structural parameters of the rolling confluence ring directly determine the performance and the service life of the rolling confluence ring. Currently, structural parameter design of a rolling bus ring is generally obtained based on a large number of real tests, and the mode based on the large number of real tests can increase manufacturing cost. And the simulation test is also based on the simulation test, but the simulation test is basically performed on the basis of a single terminal, so that the simulation efficiency is low and the simulation effect is not ideal.
Disclosure of Invention
Aiming at the technical problems, the invention adopts the following technical scheme:
the embodiment of the invention provides a rolling type confluence ring structural parameter simulation optimizing method, wherein the rolling type confluence ring at least comprises an inner ring, an outer ring and an elastic ring arranged between the inner ring and the outer ring, and the method comprises the following steps:
s100, constructing a simulation model of a target rolling type confluence ring; wherein the assembly offset of the elastic ring in the simulation model is set to D er -（D 1 or -D 2 ir ) The method comprises the steps of carrying out a first treatment on the surface of the Wherein D is er The outer diameter D of the elastic ring of the target rolling type confluence ring 1 or An inner diameter D of an outer ring of the target rolling type confluence ring 2 ir An outer diameter of an inner ring of the target rolling type confluence ring;
s200, simulating the simulation model by adopting a simulation step of multiple analysis steps to obtain a corresponding simulation result; the multi-analysis step comprises an initial analysis step, a pressure loading analysis step, an elastic ring free rebound analysis step and an inner ring rotation analysis step;
and S300, acquiring the structural parameter optimal value of the target rolling type confluence ring based on the simulation result.
Embodiments of the present invention also provide a non-transitory computer readable storage medium having stored therein at least one instruction or at least one program loaded and executed by a processor to implement the foregoing method.
The embodiment of the invention also provides an electronic device comprising a processor and the non-transitory computer readable storage medium.
The invention has at least the following beneficial effects:
according to the rolling type converging ring structural parameter simulation optimizing method, the modeling is adopted to preset the assembly offset and the multi-analysis-step setting, so that the convergence of simulation solution can be ensured, the actual operating condition of the converging ring is restored, and the optimizing efficiency of the converging ring product structure can be improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for simulating and optimizing structural parameters of a rolling type confluence ring according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a simulation model of a rolling bus ring.
Reference numerals:
1: an outer ring; 2: an inner ring; 3: an elastic ring.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to fall within the scope of the invention.
The embodiment of the invention provides a rolling type confluence ring structural parameter simulation optimizing method which is used for acquiring optimal structural design parameters when the structural design of a rolling type confluence ring is performed. In the embodiment of the present invention, the rolling type bus ring may have a conventional structure, as shown in fig. 2, and may include at least an inner ring 2, an outer ring 1, and an elastic ring 3 disposed between the inner ring 2 and the outer ring 1. When the structure design is carried out, the diameters and the wall thicknesses of the inner ring and the outer ring are fixed, and the wall thickness of the elastic ring is also fixed. The outer diameter of the elastic ring is larger than the difference between the inner diameter of the outer ring and the outer diameter of the inner ring, so that after the assembly of the three parts is completed, the elastic ring has a certain compression amount, and the compression amount can be adjusted by changing the diameter of the elastic ring.
Further, the method for optimizing the structural parameters of the rolling type bus ring according to the embodiment of the present invention may include the steps shown in fig. 1:
s100, constructing a simulation model of a target rolling type confluence ring; wherein the assembly offset of the elastic ring in the simulation model is set to D er -（D 1 or -D 2 ir ) The method comprises the steps of carrying out a first treatment on the surface of the Wherein D is er The outer diameter D of the elastic ring of the target rolling type confluence ring 1 or An inner diameter D of an outer ring of the target rolling type confluence ring 2 ir Is the outer diameter of the inner ring of the target rolling type confluence ring.
Further, S100 may specifically include:
s110, constructing a basic three-dimensional model of the target rolling type confluence ring based on the inner diameter and the outer diameter of the outer ring, the inner diameter and the outer diameter of the inner ring, the initial outer diameter and the initial inner diameter of the elastic ring and the initial interference amount between the elastic ring and the inner ring of the target rolling type confluence ring.
In the embodiment of the invention, the simulation model can be realized based on a super computer, and specifically, set simulation software can be adopted to construct a basic three-dimensional model of the target rolling type bus ring by using a Linux command. In one exemplary embodiment, a front view of the constructed basic three-dimensional model may be as shown in FIG. 2.
In the embodiment of the invention, the center position G of the elastic ring er Can be set as (-D) RC 2+θ/2, 0), where D RC For the diameter of the reference circle, the reference circle is a circle located intermediate the inner and outer rings, as shown by the broken line in FIG. 2, it is apparent that D RC =（D 1 or +D 2 ir ) 2, wherein D 1 or An inner diameter D of an outer ring of the target rolling type confluence ring 2 ir For the outer diameter of the inner ring of the target rolling type confluence ring, θ is the compression amount of the elastic ring, θ=d er -（D 1 or -D 2 ir ）。
The inventor of the invention finds through experiments that when modeling, as the diameter of the elastic ring is larger than the difference between the inner diameter of the outer ring and the outer diameter of the inner ring, if the circle center of the elastic ring is positioned on a reference circle, the interference quantity theta/2 exists between the elastic ring and both the inner ring and the outer ring, the complexity of setting boundary conditions during simulation can be increased, and the convergence of simulation results is greatly reduced.
Therefore, when modeling the elastic ring, the inventor of the invention shifts the center of the elastic ring by θ/2 to the original point of sketch coordinate along the direction of connecting the original point of sketch coordinate and the center of the elastic ring, and the center coordinate of the elastic ring is (-D) RC 2+theta/2, 0), the interference amount of the elastic ring and the outer ring is 0, the interference amount of the elastic ring and the inner ring is theta, namely the assembly offset of the elastic ring is theta, by the method, the setting complexity of boundary conditions during simulation can be reduced, and smooth convergence of simulation calculation is ensured.
S120, performing simulation pretreatment on the basic three-dimensional model to obtain a corresponding basic simulation model; the pre-simulation processing at least comprises finite element mesh division of the basic three-dimensional model and setting corresponding simulation conditions.
In the embodiment of the invention, the inner ring, the outer ring and the elastic ring can be subjected to finite element meshing by using hexahedral mesh.
In an embodiment of the present invention, the simulation conditions at least include: the material properties of the inner ring, the elastic ring and the outer ring; a simulation step; a contact relationship; boundary conditions.
Wherein the material characteristics include a selected material and corresponding material parameters, which may include density, modulus of elasticity, poisson's ratio, stress-strain values, and the like. The materials of the inner ring, the outer ring and the elastic ring may be selected based on actual needs.
The simulation step comprises a plurality of analysis steps, and specifically comprises the following steps: an initial analysis step, a pressure loading analysis step, a rebound analysis step and a rotation analysis step. In an exemplary embodiment, the types of the analysis steps are all Dynamic, the Dynamic type of the Implicit type set is dynamically Implicit, the Time period is set to 1, the nlgeom is set to On, and the load change manner with Time is Ramp linearly over step.
Further, in an embodiment of the present invention, the contact relationship may include a contact relationship between the inner ring and the elastic ring and a contact relationship between the outer ring and the elastic ring.
Wherein, the setting of the contact relation includes:
the contact relation between the inner ring and the elastic ring is set as a first contact relation in the rebound analysis step, and the first contact relation is inherited in the rotation analysis step.
The contact relationship between the outer ring and the elastic ring is set as a second contact relationship in the initial analysis step, and the second contact relationship is inherited in the pressure loading analysis step, the rebound analysis step, and the rotation analysis step.
Wherein the contact characteristic in the first contact relationship is a first contact characteristic and the contact characteristic in the second contact relationship is a second contact characteristic, each contact characteristic including a friction contact algorithm, a friction factor, and a contact type. Wherein in one exemplary embodiment, the friction contact algorithm selects a penalty function, and the value of the friction factor is set to a preset value, determined based on the material of the contact part. The contact type is selected to be a hard contact.
Further, in an embodiment of the present invention, the boundary condition includes: rigid body constraint relation between the reference point and the inner ring; rigid body constraint relation between the reference point and the outer ring; a pressure load applying region on the elastic ring and a pressure load applied to the pressure load applying region; constraint conditions of the elastic ring; the reference point is a point on the Z axis in a coordinate system to which the simulation model belongs, and can be set based on actual requirements. The origin of the coordinate system to which the simulation model belongs is the center point of the target rolling type confluence ring, the X axis is the horizontal diameter direction of the inner ring, the Y axis is the vertical diameter direction of the inner ring, and the Z axis is the central axis of the target rolling type confluence ring, as shown in fig. 2. The magnitude of the pressure load is set so that the deformation amount of the elastic ring is equal to or larger than the compression amount of the elastic ring.
Further, the setting of the boundary condition may include:
(1) And setting a rigid body constraint relation between the reference point and the inner ring in the initial analysis step as a first rigid body constraint relation, and inheriting the first rigid body constraint relation in the pressure loading analysis step, the elastic ring free rebound analysis step and the inner ring rotation analysis step.
Specifically, a 6-degree-of-freedom full constraint boundary condition of the reference point is established in the initial analysis step, and this constraint condition is repeated in the pressure loading analysis step, the rebound analysis step, and the rotation analysis step. In the rebound analysis step, the angular velocity of the inner ring in the axial direction, which is the reference point, is set to a, and if the rotational speed of the inner ring is b (in r/min), a=2pi b/60.
(2) And setting a rigid body constraint relation between the reference point and the outer ring in the initial analysis step as a second rigid body constraint relation, and inheriting the second rigid body constraint relation in the pressure loading analysis step, the rebound analysis step and the rotation analysis step.
Specifically, a 6-degree-of-freedom full constraint boundary condition of a reference point is set in the initial analysis step as a second rigid body constraint relationship, and the second rigid body constraint relationship is inherited in the pressure loading analysis step, the rebound analysis step, and the rotation analysis step.
(3) Providing a pressure load application area on the elastic ring and a pressure load applied to the pressure load application area, and providing that the pressure load is activated only in the pressure analysis step. The elastic ring is elastically deformed by applying pressure load, and the deformation is larger than or equal to the interference quantity theta of the elastic ring and the inner ring, so that the compression deformation assembly of the elastic ring is realized.
In the embodiment of the present invention, the pressure load applying area is obtained by dividing the outer ring surface of the elastic ring by the first dividing surface and the second dividing surface, and as a specific area shown by a symbol a in fig. 2, the pressure load applying area is an outer ring surface area of the elastic ring located between the first dividing surface and the second dividing surface and close to the inner ring. The first dividing surface is obtained by rotating the reference surface around the central axis of the elastic ring clockwise by a set angle alpha, and the second dividing surface is obtained by rotating the reference surface around the central axis of the elastic ring anticlockwise by a set angle alpha, namely the pressure load applying area is close to the inner ring and the corresponding central angle is equal to 2 alpha. The reference plane is a plane passing through the central axis of the elastic ring and the coordinate origin of the corresponding simulation model.
In one exemplary embodiment, α=5°.
(4) In the initial analysis step, constraint boundary conditions of the upper end face and the lower end face of the elastic ring in the Z axis are set, so that the elastic ring can only move in the XY plane.
Those skilled in the art know that the constraint boundary conditions of the upper end face and the lower end face of the elastic ring in the Z axis are set in the initial analysis step, so that the elastic ring can only move in the XY plane.
S130, adjusting the outer diameter of the elastic ring in the basic simulation model for n-1 times according to the mode that the outer diameter of the elastic ring is sequentially increased, and adjusting the related variable parameters based on the outer diameter of the elastic ring after each adjustment to obtain n simulation models. n can be a custom value, which can be determined based on the structure of the actual design rolling manifold ring.
Further, S130 may specifically include:
s1301, set j=1;
s1302, if j is less than or equal to n-1, executing S1303.
S1303, setting the outer diameter D of the elastic ring of the jth adjustment simulation model 2j or =D 20 or +j*△d；D 20 or The elastic ring is used as the outer diameter of the elastic ring of the basic simulation model, deltad is a set step value, and the units are the same as the units of the outer diameter and can be set based on actual needs.
S1304, based on D 2j or Obtaining the inner diameter D of the corresponding elastic ring 2j ir Interference quantity theta j And a pressure load P applied to the pressure load application region j The method comprises the steps of carrying out a first treatment on the surface of the Wherein P is j =D 2j or ×P j-1 /D 2（j-1） or ，D 2（j-1） or Adjusting the outer diameter, P, of the elastic ring of the simulation model for the j-1 th j-1 The pressure load applied to the pressure load application area by the simulation model is adjusted for the j-1 th.
In the embodiment of the invention, the pressure load corresponding to the basic simulation model, namely the pressure load corresponding to the model with the minimum outer diameter of the elastic ring, can be determined based on a test.
S1305, based on the basic simulation model and the n-1 adjustment simulation models, n simulation models are obtained and used as n solving files.
Further, S130 may be generated in batch by existing industrial simulation software such as abaqus. Specifically, in the optimization stage of the converging ring structure, if the elastic ring D er N design values of diameter of (D) er1 、D er2 、D er3 、……、D ern That is, the interference amount θ of the elastic ring and the inner ring is n, respectively θ 1 、θ 2 、θ 3 、……、θ n . By rapidly modifying the values of the diameter of the elastic ring in the Python file and the values of the press load of the corresponding interference quantity in the Inp file, n solving files Opt1.Inp, opt2.Inp, opt3.Inp, … … and Optn. Inp are efficiently generated.
And S200, simulating the simulation model by adopting a simulation step of multiple analysis steps to obtain a corresponding simulation result.
In an embodiment of the present invention, S200 may be implemented based on a supercomputer. Further, S200 specifically includes:
and taking the simulation model as a solving file, submitting each solving file to a corresponding computing node resource, and solving to obtain a corresponding solving result which is taken as a simulation result.
In the embodiment of the invention, because the operating system of the supercomputer is a Linux system, shell scripts are adopted to execute the operating instructions of the user. According to the simulation demand, two scripts of task submission and task scheduling are required to be written, so that each solution file is submitted to the corresponding computing node resource through the task submission script and the task scheduling script. The task submitting script is used for calling the computing node resource, identifying the solving file and calling a solver in the industrial simulation software to complete the calculation of the solving file so as to generate a solving result. The task scheduling script is used for calling the task submitting script to realize the submitting of all solution files, and the task scheduling script specifically comprises:
s201, setting an r subfolder in the specified folder, and moving an r solving file into the r subfolder; r has a value of 1 to n;
s202, copying the task submission script into an r subfolder, and modifying the name of a solving file submitted by the task submission script in the r subfolder into a corresponding designated file name C1r;
s203, executing a task submitting script in the r subfolder to submit the corresponding solving file;
s204, setting r=r+1, executing S201 if r is less than or equal to n, otherwise, exiting the current control program.
Specifically, taking Abaqus as an example, the setting of the task submission script and the task scheduling script includes:
s1, naming the name of a task submitting script as a subset, wherein the storage position of the task submitting script is a folder Opt, and the realized functions are as follows: invoking a computing node resource of the supercomputer, identifying an Abaqus solving file (such as Opt1. Inp) in a folder, invoking an Abaqus solver installed in the supercomputer system to complete the calculation of the solving file, generating a solving result Optx.odb (x=1, 2,3, … …, n) file, and storing the solving result Opt in the Opt.
S2, naming the name of the task scheduling script as dispatch.sh, and realizing the functions in the following order:
(1) Creating a new folder Optx in the Opt folder; (2) Moving the Optx.inp solving file in the Opt folder to the Optx of the folder; (3) Copying the subset sh in the Opt folder to the folder Optx; (4) enter folder Optx; (5) Modifying the name of a solving file submitted by submit.sh in the file Optx as Optx.inp; (6) Adopting a task submission command to submit a task for solving the sub. (7) returning to the upper folder Opt; (8) The commands (1) to (7) are repeated in sequence until the n solving tasks (opt 1.Inp to opt. Inp) are submitted to end.
Further, in S200, the h solution file passes through the corresponding computing node resource R h Invoking the corresponding solver S h The method comprises the following steps of:
s210, computing node resource R h Based on the h-th calculationSolving the h solving file under the simulation condition corresponding to the solving file; if the elastic analysis step is not converged, executing S212, otherwise, obtaining a corresponding solving result and storing the solving result into a designated folder; h has a value of 1 to n.
S212, obtaining the displacement delta x of the node tangent to the inner ring of the elastic ring corresponding to the h solving file h I.e. the difference in horizontal coordinates between the initial state and the compressed state, if Deltax h ＜θ h If the current pressure load is insufficient, the current pressure load corresponding to the h-th solution file needs to be increased, and S210 is executed.
Those skilled in the art will appreciate that obtaining the displacement of the node at which the elastic ring is tangent to the inner ring, corresponding to the h solution file, may be a prior art.
In the embodiment of the invention, the increment value of the current pressure load is equal to the current pressure load multiplied by a preset coefficient b, and b is more than 0 and less than 1.b may be an empirical value. In one exemplary embodiment, b=0.05 to 0.1, preferably b=0.05.
Those skilled in the art will appreciate that the elastic analysis step non-convergence may be determined based on prior art techniques.
S300, acquiring target structural parameters of the target rolling type bus ring based on the simulation result.
In an embodiment of the present invention, S300 may be implemented by a supercomputer. Further, S300 may specifically include:
s310, processing the n simulation results to obtain n corresponding processing results, wherein each processing result comprises stress values of all nodes of the corresponding elastic ring, contact force values of all nodes of the outer ring surface of the elastic ring, and maximum stress values and maximum contact force values.
In the embodiment of the invention, the processing of the solving result is realized through the solving result extraction script and the result processing script. The solving result extraction script is used for extracting stress values of all nodes (namely the top points of each hexahedron) of the elastic ring, contact force values of all nodes of the outer ring surface of the elastic ring, and maximum stress values and maximum contact force values in the rebound analysis step; and the result processing script is used for calling the solving result extraction script to process the solving result so as to obtain the processing result.
Further, in an exemplary embodiment, the embodiments of the present invention perform efficient processing of rolling-type torus solution results based on Python scripts. Specifically, the method can comprise the following steps:
s10, writing a Python script file handle.py, wherein the storage position of the Python script file is a folder Opt, and the function of the Python script file is to extract stress values of all nodes of the Step-revving analysis Step elastic ring and contact force values of all nodes of the outer ring surface of the elastic ring in the Optx.odb result, obtain the maximum value of the stress values and the contact force values, and finally output a result file Maxx.txt.
S20, writing a result processing script handle. Sh, wherein the function of the result processing script handle. Sh is to call a Python script file and perform post-processing on a simulation result odb. The function is realized in the following order:
(1) Creating a new folder Max in the Opt folder; (2) Copy handle. Py in the Opt folder into folder Optx; (3) enter folder Optx; (4) Modifying the name of the solving result file processed by handle. Py in the folder Optx to be Optx.odb; (5) Executing a handle.py script to perform the result data processing of the Optx.odb, and generating a data processing result file Maxx.txt; (6) copying maxx.txt into folder Max; (7) returning to the upper folder Opt; (8) And (3) sequentially repeating the commands (2) to (7) until the processing of the Opt1.Odb to Optn. Odb n solving results is finished.
S320, based on the n processing results, obtaining a target processing result to obtain an optimal value of the structural parameter of the target rolling type confluence ring.
In the embodiment of the invention, a data processing result folder Max is read from a super computer system to a local computer, n.txt processing result files are finally processed, and an optimal result is obtained according to a structural parameter optimization standard, so that a target structural parameter of the rolling type confluence ring, namely an optimal structural parameter, particularly the outer diameter of an elastic ring is obtained. In the embodiment of the invention, according to the structural parameter optimization standard, the performance of the rolling type confluence ring which is actually designed can be determined, and the invention is not particularly limited.
In summary, the rolling type confluence ring structural parameter simulation optimizing method provided by the embodiment of the invention has at least the following advantages:
(1) The rolling type confluence ring flow based on the supercomputer can efficiently carry out batch modeling, boundary condition setting, grid division, calculation task solving and result post-processing;
(2) The preset of the elastic ring assembly offset can reduce the complexity of setting the simulation boundary conditions and ensure the convergence of the solving result;
(3) The compression assembly, free contact extrusion and friction rolling of the rolling type confluence ring elastic ring are realized by adopting multiple analysis steps, and the actual loading condition of the confluence ring can be restored;
(4) Under the condition of different compression amounts of the converging ring, the compression load can be determined by testing for a plurality of times, the preset proportion value is adopted, and then the load value is increased through the circulation command, so that the situation that the calculation task is not converged can be avoided.
Embodiments of the present invention also provide a non-transitory computer readable storage medium that may be disposed in an electronic device to store at least one instruction or at least one program for implementing one of the methods embodiments, the at least one instruction or the at least one program being loaded and executed by the processor to implement the methods provided by the embodiments described above.
Embodiments of the present invention also provide an electronic device comprising a processor and the aforementioned non-transitory computer-readable storage medium.
Embodiments of the present invention also provide a computer program product comprising program code for causing an electronic device to carry out the steps of the method according to the various exemplary embodiments of the invention as described in the specification, when said program product is run on the electronic device.
While certain specific embodiments of the invention have been described in detail by way of example, it will be appreciated by those skilled in the art that the above examples are for illustration only and are not intended to limit the scope of the invention. Those skilled in the art will also appreciate that many modifications may be made to the embodiments without departing from the scope and spirit of the invention. The scope of the present disclosure is defined by the appended claims.
Claims (10)
1. The rolling type confluence ring structure parameter simulation optimizing method is characterized by comprising the following steps of:
s100, constructing a simulation model of a target rolling type confluence ring; wherein the assembly offset of the elastic ring in the simulation model is set to D er -（D 1 or -D 2 ir ) The method comprises the steps of carrying out a first treatment on the surface of the Wherein D is er The outer diameter D of the elastic ring of the target rolling type confluence ring 1 or An inner diameter D of an outer ring of the target rolling type confluence ring 2 ir An outer diameter of an inner ring of the target rolling type confluence ring;
s200, simulating the simulation model by adopting a simulation step of multiple analysis steps to obtain a corresponding simulation result; the multi-analysis step comprises an initial analysis step, a pressure loading analysis step, an elastic ring free rebound analysis step and an inner ring rotation analysis step;
and S300, acquiring the structural parameter optimal value of the target rolling type confluence ring based on the simulation result.
2. The method according to claim 1, wherein S100 specifically comprises:
s110, constructing a basic three-dimensional model related to a target rolling type confluence ring based on the inner diameter and the outer diameter of an outer ring of the target rolling type confluence ring, the inner diameter and the outer diameter of an inner ring, the initial outer diameter of an elastic ring, the initial inner diameter and the initial interference quantity between the elastic ring and the inner ring;
s120, performing simulation pretreatment on the basic three-dimensional model to obtain a corresponding basic simulation model; the pre-simulation processing at least comprises finite element mesh division of the basic three-dimensional model and setting corresponding simulation conditions;
s130, adjusting the outer diameter of the elastic ring in the basic simulation model for n-1 times according to the mode that the outer diameter of the elastic ring is sequentially increased, and adjusting the related variable parameters based on the outer diameter of the elastic ring after each adjustment to obtain n simulation models.
3. The method according to claim 2, wherein S200 specifically comprises:
and taking the simulation model as a solving file, submitting each solving file to a corresponding computing node resource, and solving to obtain a corresponding solving result which is taken as a simulation result.
4. A method according to claim 3, characterized in that in S200 the h-th solution file is passed through the corresponding computing node resource R h Invoking the corresponding solver S h The method comprises the following steps of:
s210, computing node resource R h Solving the h solving file based on simulation conditions corresponding to the h solving file; if the elastic analysis step is not converged, executing S212, otherwise, obtaining a corresponding solving result and storing the solving result into a designated folder; h has a value of 1 to n;
s212, obtaining the displacement delta x of the node tangent to the inner ring of the elastic ring corresponding to the h solving file h If Deltax h ＜θ h The current pressure load corresponding to the h-th solution file is increased, and S210 is performed.
5. A method according to claim 3, wherein S300 comprises:
s310, processing n simulation results to obtain n corresponding processing results, wherein each processing result comprises stress values of all nodes of the corresponding elastic ring, contact force values of all nodes on the outer ring surface of the elastic ring, and maximum stress values and maximum contact force values;
s320, based on the n processing results, obtaining a target processing result to obtain an optimal value of the structural parameter of the target rolling type confluence ring.
6. The method according to claim 2, wherein the simulation conditions comprise at least: the material properties of the inner ring, the elastic ring and the outer ring; a simulation step of multiple analysis steps; a contact relationship; boundary conditions;
the contact relation comprises a contact relation between the inner ring and the elastic ring and a contact relation between the outer ring and the elastic ring;
the boundary conditions include: rigid body constraint relation between the reference point and the inner ring; rigid body constraint relation between the reference point and the outer ring; a pressure load applying region on the elastic ring and a pressure load applied to the pressure load applying region; constraint conditions of the elastic ring; the reference point is a point on a Z axis in a coordinate system to which the simulation model belongs, wherein an origin of the coordinate system to which the simulation model belongs is a center point of the target rolling type confluence ring, an X axis is a horizontal diameter direction of the inner ring, a Y axis is a vertical diameter direction of the inner ring, and the Z axis is a central axis of the target rolling type confluence ring; the magnitude of the pressure load is set so that the deformation amount of the elastic ring is equal to or larger than the compression amount of the elastic ring.
7. The method of claim 6, wherein the pressure load application region is divided by a first dividing surface and a second dividing surface facing an outer circumferential surface of the elastic ring;
the first dividing surface is obtained by clockwise rotation of the reference surface around the central axis of the elastic ring by a set angle alpha, and the second dividing surface is obtained by anticlockwise rotation of the reference surface around the central axis of the elastic ring by a set angle alpha; the reference surface is a plane passing through the central axis of the elastic ring and the coordinate origin of the corresponding simulation model; the pressure load application region is adjacent to the inner ring and the corresponding central angle is equal to 2α.
8. The method of claim 7, wherein α = 5 °.
9. A non-transitory computer readable storage medium having stored therein at least one instruction or at least one program, wherein the at least one instruction or the at least one program is loaded and executed by a processor to implement the method of any one of claims 1-8.
10. An electronic device comprising a processor and the non-transitory computer readable storage medium of claim 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310976666.XA CN116680842B (en) | 2023-08-04 | 2023-08-04 | Rolling type confluence ring structural parameter simulation optimizing method, equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310976666.XA CN116680842B (en) | 2023-08-04 | 2023-08-04 | Rolling type confluence ring structural parameter simulation optimizing method, equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116680842A true CN116680842A (en) | 2023-09-01 |
| CN116680842B CN116680842B (en) | 2023-10-20 |
Family
ID=87789530
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310976666.XA Active CN116680842B (en) | 2023-08-04 | 2023-08-04 | Rolling type confluence ring structural parameter simulation optimizing method, equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116680842B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107045571A (en) * | 2017-03-27 | 2017-08-15 | 杭州远算科技有限公司 | A kind of die casting simulation control method based on cloud computing |
| CN110532637A (en) * | 2019-08-05 | 2019-12-03 | 上海电气风电集团有限公司 | One kind being suitable for double feedback electric engine electric brush slip ring system temperature field emulated computation method |
| CN112966467A (en) * | 2021-03-17 | 2021-06-15 | 武汉大学 | Flexible PCB structure design method based on wet chemical etching combined simulation |
| US20210220884A1 (en) * | 2020-01-20 | 2021-07-22 | Yanshan University | Dynamic Contact Heat Transfer Simulation Device for Rolling Heavy-Load Deformation Zone |
| CN114818118A (en) * | 2022-04-01 | 2022-07-29 | 中国第一汽车股份有限公司 | Air spring system structural member design method, device, equipment and storage medium |
| CN115618698A (en) * | 2022-12-20 | 2023-01-17 | 天河超级计算淮海分中心 | Finite element simulation method based on falling scene, electronic equipment and storage medium |
-
2023
- 2023-08-04 CN CN202310976666.XA patent/CN116680842B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107045571A (en) * | 2017-03-27 | 2017-08-15 | 杭州远算科技有限公司 | A kind of die casting simulation control method based on cloud computing |
| CN110532637A (en) * | 2019-08-05 | 2019-12-03 | 上海电气风电集团有限公司 | One kind being suitable for double feedback electric engine electric brush slip ring system temperature field emulated computation method |
| US20210220884A1 (en) * | 2020-01-20 | 2021-07-22 | Yanshan University | Dynamic Contact Heat Transfer Simulation Device for Rolling Heavy-Load Deformation Zone |
| CN112966467A (en) * | 2021-03-17 | 2021-06-15 | 武汉大学 | Flexible PCB structure design method based on wet chemical etching combined simulation |
| CN114818118A (en) * | 2022-04-01 | 2022-07-29 | 中国第一汽车股份有限公司 | Air spring system structural member design method, device, equipment and storage medium |
| CN115618698A (en) * | 2022-12-20 | 2023-01-17 | 天河超级计算淮海分中心 | Finite element simulation method based on falling scene, electronic equipment and storage medium |
Non-Patent Citations (3)
| Title |
|---|
| YIZHOU YANG ET AL.: "Simulation and experimental analysis of rolling element bearing fault in rotor-bearing-casing system", ENGINEERING FAILURE ANALYSIS, vol. 92, pages 205 - 221 * |
| 杨金川 等: "滚动汇流环装配力学特性分析与优化设计", 工程设计学报, vol. 23, no. 4, pages 364 - 370 * |
| 胡扬: "一种新型滚动集电环的设计研究", 中国优秀硕士学位论文全文数据库工程科技Ⅱ辑, vol. 2020, no. 02, pages 042 - 246 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116680842B (en) | 2023-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110298443B (en) | Neural network operation device and method | |
| AU2019268193B2 (en) | Conversion method, apparatus, computer device, and storage medium | |
| US11281824B2 (en) | Authoring loading and boundary conditions for simulation scenarios | |
| US10331821B1 (en) | Automatic variable fidelity simulation | |
| EP2426564B1 (en) | System and method for solving quadratic programming problems with bound constraints utilizing a semi-explicit quadratic programming solver | |
| CN113221478B (en) | Numerical analysis method and system for actively controlling rotor wing | |
| WO2021057807A1 (en) | Deep learning model generation method and apparatus, device, and storage medium | |
| Zhang et al. | Adaptive ANCF method and its application in planar flexible cables | |
| CN107273569B (en) | Aerodynamic configuration drag reduction optimization method based on mesh deformation technique | |
| EP2887241A1 (en) | Cad-based initial surface geometry correction | |
| CN110955941B (en) | Vector field-based composite material structure optimization design method and device | |
| WO2018058427A1 (en) | Neural network computation apparatus and method | |
| CN116245049A (en) | Node type non-structural grid boundary correction method, device, equipment and medium | |
| CN116680842B (en) | Rolling type confluence ring structural parameter simulation optimizing method, equipment and medium | |
| Mastenko et al. | Generative design of a frame type construction | |
| JP2019046038A (en) | Evaluation function converter and program | |
| CN117010109B (en) | Rolling type confluence ring structural parameter optimization method based on super calculation | |
| CN112414668B (en) | Wind tunnel test data static bomb correction method, device, equipment and medium | |
| JP6901423B2 (en) | Information processing equipment, information processing terminals, and programs | |
| JP2021535522A (en) | Artificial intelligence execution model setting method and artificial intelligence execution acceleration system for accelerating the execution of artificial intelligence | |
| CN116933386A (en) | Aircraft pneumatic data fusion method based on MCOK proxy model | |
| CN113392436B (en) | CAE analysis software and office software integrated control method and platform | |
| CN110580528A (en) | network model generation method and device | |
| CN113076667A (en) | Equivalent simulation method, equivalent simulation system, computer device and computer-readable storage medium | |
| CN102063308A (en) | Method for controlling processing flow of seismic prospecting data |
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 |