A design method for impact-resistant chiral vibration isolation structure combining Grasshopper and Abaqus

Abstract

The present invention discloses a method for designing an impact-resistant chiral vibration isolation structure combining Grasshopper and Abaqus, and belongs to the technical field of computer-aided structural design and anti-impact and vibration isolation. The steps for implementing the method include: S1: batch designing impact-resistant chiral vibration isolation structures based on Grasshopper, generating and storing models in sequence; S2: batch reading models based on the Abaqus platform and conducting mechanical simulation; S3: batch reading finite element simulation results, analyzing the relationship between structural mechanical properties and structural parameters; S4: determining and optimizing the design scheme of the impact-resistant chiral vibration isolation structure. The present invention uses Grasshopper to perform parametric modeling of chiral structures, and at the same time builds a computing platform for large-scale finite element simulation and post-processing analysis of solid structures, providing a fast and diversified optimization scheme for designing devices with both equipment anti-impact and vibration isolation performance, which helps to achieve good impact resistance and low-frequency vibration isolation performance under a wider spectrum for complex and changeable industrial application scenarios.

Description

Technical Field

[0001] The invention belongs to the technical field of computer-aided structural design and shock resistance and vibration isolation, and in particular relates to a design method for an shock-resistant chiral vibration isolation structure. Background Art

[0002] Equipment impact resistance and vibration reduction and isolation have always been the focus of attention in the industrial field. They are related to equipment safety, performance and life, and are widely used in the national defense and military fields. Strengthening the equipment's impact resistance and enhancing the equipment's vibration reduction and isolation performance in a wide frequency domain are of great significance for expanding the applicability of equipment, improving equipment performance, and extending equipment service life. The vibration isolation devices of existing equipment have the following technical bottlenecks that need to be solved: they are mostly designed for vibration isolation needs, and additional rigid limit protection is often used for impact, which makes the equipment very easy to be damaged by overload when impacted, and is prone to secondary impact damage; the traditional vibration isolation single frequency vibration isolation principle is used, and the vibration reduction and isolation effect in a wide low-frequency range is not good; the impact resistance and vibration reduction and isolation scenarios applicable to a single device are relatively single, and it is difficult to flexibly adapt to complex and changing industrial application scenarios.

[0003] By using metamaterial design technology, specific mechanical properties can be achieved through special material structure design, which has broad application prospects in the field of vibration and noise reduction. The chiral energy-locking structure, with its unique deformation mechanism, can convert axial vibration and impact into circumferential deformation through torsion, thereby enhancing the axial impact resistance and low-frequency vibration reduction effect of the device.

[0004] The main component of the chiral energy lock structure is a spatial cylindrical rod. In the design, not only the diameter and length of the rod itself must be considered, but also the angle between the rod and the horizontal plane, and the rotation angle of the rod end in the horizontal plane. Existing parametric design software such as CATIA, Rhino, PRO / E, etc. have the following three defects:

[0005] 1. Parametric design is not logical and cannot simultaneously associate the size of the rod, its angle with the horizontal plane, the rotation angle of the end of the rod within the horizontal plane and many other variables, and define the intrinsic functional relationship between the variables.

[0006] 2. The visualization level of parametric design is not high, the interactive interface is not friendly, and the design results cannot be seen immediately by directly dragging parameters.

[0007] 3. The parametric design space is small, and it is not possible to generate and verify a large number of chiral structure models, so it is not possible to provide more design options to meet actual application needs. Summary of the invention

[0008] In view of this, the present invention provides a design method for impact-resistant chiral vibration isolation structure combining Grasshopper and Abaqus, and uses Grasshopper to perform parametric modeling of chiral structure, realizes multi-parameter design and construction of structure, real-time dynamic adjustment and display of structural geometry, and batch generation of parametric models; at the same time, a computing platform for large-scale solid structure finite element simulation and post-processing analysis is built, which provides a fast and diversified optimization scheme for designing devices with both equipment impact resistance and vibration isolation performance, and helps to achieve good impact resistance and low-frequency vibration isolation performance in a wider spectrum for complex and changeable industrial application scenarios.

[0009] A method for designing a shock-resistant chiral vibration isolation structure combining Grasshopper and Abaqus, the method comprising the following steps:

[0010] S1: Batch design of impact-resistant chiral vibration isolation structures based on Grasshopper, generate and store models in sequence;

[0011] S2: Batch read models and carry out mechanical simulation based on the Abaqus platform;

[0012] S3: Batch read finite element simulation results and analyze the relationship between structural mechanical properties and structural parameters;

[0013] S4: Determine and optimize the design scheme of impact-resistant chiral vibration isolation structure based on actual application requirements.

[0014] Furthermore, the step S1 comprises:

[0015] S11: Determine the composition and material selection of the impact-resistant chiral vibration isolation structure;

[0016] The structure is composed of a chiral energy locking unit, a metal gasket, upper and lower adapter plates, bolts and nuts, wherein the chiral energy locking unit includes a cylindrical rod, an intermediate connecting piece, a unit upper end interface and a unit lower end interface, the intermediate connecting piece is an annular hollow shell, a horizontal metal gasket slot is processed on the outer circumferential surface of the shell, and two partitions are provided at positions corresponding to the metal gasket slots on the inner circumferential surface of the shell, and the two partitions extend from the inner circumferential surface of the shell and are in the same horizontal plane. After the metal gasket is inserted into the metal gasket slot, the metal gasket and the two partitions together form a planar structure, and after the structure is impacted, the structural deformation of the unit upper end interface and the unit lower end interface is limited;

[0017] The outer shapes of the upper interface of the unit and the lower interface of the unit are both cylindrical structures. The upper end of the upper interface of the unit has an external thread structure, and the lower end of the lower interface of the unit has an internal thread structure. The internal and external thread specifications of the upper interface and the lower interface are consistent;

[0018] The upper interface of the unit and the lower interface of the unit are respectively connected to the upper and lower ends of the middle connecting piece through cylindrical rods, the cylindrical rods are arranged periodically along the circumferential direction of the middle connecting piece, and the cylindrical rods at the upper and lower ends of the middle connecting piece are arranged in mirror symmetry;

[0019] Multiple chiral energy locking units are sequentially connected together in the vertical direction through the connection interfaces at the upper and lower ends to form a chiral energy locking unit group. Multiple chiral energy locking unit groups are distributed in the horizontal plane and fixedly connected through upper and lower adapter plates and bolts and nuts to form a whole. The metal gasket is inserted into the metal gasket slot of the chiral energy locking unit. The upper end interface and the lower end interface of the unit are both made of PLA or ABS or other polymer materials with good elasticity through 3D printing.

[0020] S12: Using the commercial 3D NURBS industrial modeling software Rhinoceros (Rhino) and its supporting parametric design plug-in Grasshopper, a parametric design platform for impact-resistant chiral vibration isolation structures was built, which includes four modules: parameter adjustment, unit cell modeling, array modeling, and model storage.

[0021] S13: parameter adjustment module. The function of the parameter adjustment module is to change the design parameters of the impact-resistant chiral vibration isolation structure, the angle θ between the cylindrical rod and the horizontal plane, the torsion angle β between the upper and lower end planes of the intermediate connecting member, and the ratio γ between the length of the cylindrical rod and the radius of the cylindrical rod. The parameter adjustment module uses the time counting timer command to sequentially read the parameter combination compiled by Matlab, so as to continuously adjust the structural design parameters. At the same time, the output of this module will also sequentially modify the storage file name of the corresponding model.

[0022] S14: Unit cell modeling module. The unit cell modeling module is used to design and generate chiral energy locking units. According to the given parameter combination, the circle drawing, offset, rotation and sine and cosine calculation commands are used to generate the trajectory lines of the center points of the upper and lower ends of the cylindrical rod. The divide command is used to divide the points into equal parts, and the cylindrical rods are arranged in an annular direction according to the equal points. The upper and lower end points are connected in sequence using the pipe command to generate a solid rod structure; the upper end of the cylindrical rod is rotated to a specific angle, and it is determined whether the annular arrangement interferes;

[0023] S15: Array modeling module, which is used to design and generate metal gaskets, unit upper end interfaces, and unit lower end interfaces; mirror the cylindrical rods arranged in the annular direction obtained above, and use Boolean operation commands to design intermediate connectors and slots; design the unit upper end interface and unit lower end interface, and perform Boolean operations on the combined unit upper end interface and the annular cylindrical rods; debug and generate the size coefficients and scaling ratios of the metal gaskets, unit upper end interfaces, and unit lower end interfaces to ensure reasonable coordination of the structure and meet the feasibility of actual processing and manufacturing;

[0024] S16; model storage module, the model storage module is used to sequentially generate model file names and store models, set the file storage path and file name according to the established parameter combination, number them in sequence according to the parameter sequence, and then bake them into the Rhino space and store them separately; the model storage module outputs four geometric models in total, namely a single cylindrical rod, a chiral energy locking unit, a longitudinal array of chiral energy locking units and a metal gasket, and the output geometric model is modified according to actual needs.

[0025] Furthermore, the process of batch reading models and performing mechanical simulation based on the Abaqus platform in step S2 includes:

[0026] S21: Using the commercial finite element software ABAQUS, a mechanical calculation finite element model for a cylindrical rod in a single chiral energy locking unit is established through the CAE interface; under axisymmetric coordinates, two constraints are set on the lower end of the cylindrical rod, one is completely fixed, and the other is only allowed to rotate in the circumferential direction, which is used to evaluate the influence of structural torsion on the static performance of the structure; the upper end of the cylindrical rod is rigidly coupled with the reference point, and is set to only have vertical displacement for displacement loading. The loading displacement is set according to the established model, and the value is twice the vertical height of the cylindrical rod; the material is selected as a polymer material, the modulus is 1MPa, and the density is 1.2e-9t / mm 3 ; The cylindrical rod adopts solid unit, and the loading adopts the quasi-static setting in implicit dynamics. The force, torque, displacement of the upper reference point and the strain energy of the whole model are set as historical output variables;

[0027] S22: Write a python script command to import a single cylindrical member model into ABAQUS and generate an inp calculation file. In the script, all objects in the Abaqus module should be imported to provide access to the default model database, mainly including importing the symbolic constant module for calling the predefined constants in Abaqus; importing all objects in the caeModules module; the script defines the model as a solid and assigns material properties, sets the minimum unit size and divides the grid, sets the vertical loading displacement, and finally generates an inp calculation file, and sets multi-core CPU calculation according to the actual situation of the computer; through this python script command, the single cylindrical member model is imported through the script in the CAE interface of S21;

[0028] S23: For a given parameter combination, Matlab is used to compile a python script file that batch imports cylindrical member models into ABAQUS and batch generates inp calculation files; based on the python script command written in S22, the array storing the structural parameter number and the array storing the structural parameter are used to batch modify the model name and parameters in the python script command according to the text ASCII code to generate a new python script command, and the python script command is run in the CAE interface of S21 to batch import cylindrical member models and generate inp calculation files;

[0029] S24: For a given parameter combination, Matlab is used to generate a script file that can batch submit multiple cylindrical bar model inp calculation files, and can batch complete the finite element calculation of cylindrical bar models. At the same time, multi-core CPU calculation settings can be performed according to actual conditions.

[0030] Furthermore, the process of batch reading finite element simulation results in step S3 and analyzing the relationship between structural mechanical properties and structural parameters includes:

[0031] S31: Write a python script command to extract the finite element calculation results of a single cylindrical member; import all objects in the Abaqus module in the script to provide access to the default model database, including importing the symbolic constant module for calling the predefined constants in Abaqus; import all objects in the caeModules module;

[0032] S32: For a given parameter combination, a python script file for batch extracting finite element calculation results is compiled using Matlab; based on the python script command written in S31, the model name and parameters in the python script command are replaced by batch modification according to the text ASCII code using the array storing the structural parameter number and the array storing the structural parameter;

[0033] S33, using Matlab to compile a batch finite element simulation result processing program, systematically analyze the static performance of cylindrical rods under different parameter combinations, including load, strain energy and torque.

[0034] Furthermore, the process of determining and optimizing the design scheme of the impact-resistant chiral vibration isolation structure in step S4 includes: clarifying the overall external dimensions and the number of array arrangements of the impact-resistant chiral vibration isolation structure according to the size restrictions in the actual application scenario, and then setting the parameters of the required chiral energy locking units; carrying out multi-parameter chiral energy locking element mechanics simulations in batches according to the above steps S2 and S3, and analyzing the relationship between the mechanical properties of the structure and the structural parameters; taking the maximum bearing capacity and ultimate deformation of the impact-resistant chiral vibration isolation structure as the optimization goals, and determining and optimizing the design scheme of the impact-resistant chiral vibration isolation structure.

[0035] Beneficial effects:

[0036] 1. The present invention is a design method for impact-resistant chiral vibration isolation structures combining Grasshopper and Abaqus. The chiral structure is parametrically modeled by utilizing the visualization feature of the Grasshopper operation process. Combined with the plug-in in Grasshopper, the multi-parameter design and construction of the structure, real-time dynamic adjustment and display of the structural geometry, batch generation of parametric models and other goals can be achieved, which greatly improves the speed of model generation and modification, thereby effectively improving work efficiency.

[0037] 2. The present invention can be extended to the design of most chiral surround unit cells and array structures. It has wide universality and only needs to adjust some parameters to make the design perfectly meet diverse needs. The parametric modeling process adopted by the invention can accurately control the size and positioning of the model, which is convenient for later processing, manufacturing and assembly. Based on the method proposed in the present invention, it can provide diversified optimization schemes for designing devices with both equipment impact resistance and vibration isolation performance, which helps to achieve good impact resistance and low-frequency vibration isolation performance in a wider spectrum for complex and changeable industrial application scenarios.

[0038] 3. The present invention makes parametric design more intrinsically logical, and can simultaneously associate many variables such as the rod size, the angle between the rod and the horizontal plane, the rotation angle of the rod end in the horizontal plane, and define the intrinsic functional relationship between the variables; at the same time, the parametric design has a high degree of visualization and a friendly interactive interface, and the design results can be seen instantly by directly dragging the parameters, which effectively improves the design effect; in addition, the parametric design space of the present invention is large and the variables are rich. Because it is based on programs and modular ideas, it provides the possibility for a large number of chiral structure models to be generated and verified, and provides more design options to meet actual application needs.

[0039] 3. In step S22 of the present invention, the python script command is used to implement the import of a single cylindrical rod model in the CAE interface of S21 through a script; compared with manual modeling in the CAE interface, work efficiency can be greatly improved.

[0040] 4. In step S23 of the present invention, for a given parameter combination, Matlab is used to compile a python script file that batch imports cylindrical rod models into ABAQUS and batch generates inp calculation files; based on the python script command written in S22, an array for storing structural parameter numbers and an array for storing structural parameters are used to batch modify the model names and parameters in the python script command according to the text ASCII code to generate a new python script command, and the python script command is run in the CAE interface of S21 to batch import cylindrical rod models and generate inp calculation files; this method avoids the operation of manually modifying and running the S22 python script command alone, and greatly improves the time for generating inp calculation files for multiple cylindrical rod models.

[0041] 5. In step S31 of the present invention, a python script command is written to extract the finite element calculation results of a single cylindrical rod; all objects in the Abaqus module are imported into the script to provide access to the default model database, including importing the symbolic constant module for calling the predefined constants in Abaqus; all objects in the caeModules module are imported; the python script command can open the finite element calculation result file of a single cylindrical rod, output the physical quantities such as the axial force of the reference point, the torque of the reference point, the structural kinetic energy, the structural strain energy and so on and store them as a document; compared with manually extracting the results in the CAE interface, the work efficiency can be greatly improved.

[0042] 6. In step S32 of the present invention, for a given parameter combination, Matlab is used to compile a python script file for batch extracting finite element calculation results; based on the python script command written in S31, an array for storing structural parameter numbers and an array for storing structural parameters are used to batch modify the model name and parameters in the python script command according to the text ASCII code to generate a new python script command; by running the python script command in the CAE interface of S31, the finite element calculation results of cylindrical rods can be batch extracted. This method avoids the operation of manually modifying and running the S31 python script command separately, and greatly improves the time required for batch extraction of finite element calculation results of multiple cylindrical rods.

[0043] 7. Step S33 of the present invention uses Matlab to compile a batch finite element simulation result processing program to systematically analyze the static properties of cylindrical rods under different parameter combinations, including load, strain energy and torque; the program distinguishes and analyzes the two situations of bottom sliding and no sliding, and can intuitively see the influence of the bottom constraint on the static properties of the chiral energy locking unit cell structure. At the same time, the program can analyze the influence of the angle θ between the cylindrical rod and the horizontal plane, the torsion angle β of the upper and lower end planes, and the ratio γ of the rod length to the rod radius on the static properties of the cylindrical rod, including the maximum force value, locked energy, released energy and torque. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 A flow chart of the steps in designing a method for the invention;

[0045] Figure 2 It is a schematic diagram of the impact-resistant chiral vibration isolation structure;

[0046] Figure 3 Schematic diagram of the structure of the chiral energy locking unit;

[0047] Figure 4 It is a diagram of the static performance of cylindrical rods under different parameter combinations (load, strain energy, torque, etc.);

[0048] Figure 5 It is a schematic diagram of the distribution of the maximum force, locked energy, released energy, torque, etc. of the cylindrical rod under different parameters;

[0049] Figure 6 Schematic diagram of the numerical ratio of maximum force, locked energy, released energy, torque, etc. for cylindrical rods with different parameters.

[0050] Among them, 1-chiral locking unit, 11-cylindrical rod, 12-intermediate connecting piece, 13-metal gasket slot, 14-unit upper end interface, 15-unit lower end interface, 2-metal gasket, 3-upper and lower adapter plates, 4-bolts and nuts. DETAILED DESCRIPTION

[0051] The present invention is described in detail below with reference to the accompanying drawings and embodiments.

[0052] As attached Figure 1 As shown, the present invention provides a design method for an impact-resistant chiral vibration isolation structure combining Grasshopper and Abaqus, and the implementation steps include:

[0053] S1: Batch design of impact-resistant chiral vibration isolation structures based on Grasshopper, generate and store models in sequence;

[0054] S2: Batch read models and carry out mechanical simulation based on the Abaqus platform;

[0055] S3: Batch read finite element simulation results and analyze the relationship between structural mechanical properties and structural parameters;

[0056] S4: Determine and optimize the design scheme of impact-resistant chiral vibration isolation structure based on actual application requirements.

[0057] In a preferred embodiment of the present invention, the above step S1 includes:

[0058] S11, determine the main components and material selection of the impact-resistant chiral vibration isolation structure.

[0059] As attached Figure 2 As shown, the impact-resistant chiral vibration isolation structure includes a chiral energy locking unit 1, a metal gasket 2, an upper and lower adapter plate 3, and upper and lower fixing bolts and nuts 4. A plurality of chiral energy locking units 1 are sequentially connected together in the vertical direction through the connection interfaces at the upper and lower ends to form a chiral energy locking unit group, and a plurality of chiral energy locking unit groups are distributed in the horizontal plane and fixedly connected through the upper and lower adapter plates 3 and the upper and lower fixing bolts and nuts 4 to form a whole, and the metal gasket 2 is inserted into the metal gasket slot of the chiral energy locking unit 1.

[0060] As attached Figure 3 As shown, the chiral energy lock unit 1 includes: a cylindrical rod 11, an intermediate connecting member 12, a unit upper end interface 14 and a unit lower end interface 15. Among them, the intermediate connecting member 12 is a ring-shaped hollow shell, and a horizontal metal gasket slot 13 is processed on the outer circumferential surface of the shell. There are two partitions at the position corresponding to the metal gasket slot 13 on the inner circumferential surface of the shell. The two partitions extend from the inner circumferential surface of the shell and are in the same horizontal plane. After the metal gasket 2 is inserted into the metal gasket slot 13, the metal gasket 2 and the two partitions together form a planar structure. After the structure is impacted, the structural deformation of the unit upper end interface 14 and the unit lower end interface 15 is limited.

[0061] The upper interface 14 of the unit is designed with a protruding thread, and the lower part is connected to the cylindrical rod 11; the lower interface 15 of the unit is designed with an inner concave thread, and the upper part is connected to the cylindrical rod 11. The internal and external thread specifications of the upper interface 14 and the lower interface 15 are consistent, so that multiple chiral energy locking units 1 can be combined by thread connection. The cylindrical rod 11, the intermediate connector 12, the upper interface 14 of the unit and the lower interface 15 of the unit are all made of PLA or ABS or other polymer materials with good elasticity and 3D printing molding.

[0062] S12, using the commercial 3D NURBS industrial modeling software Rhinoceros (Rhino) and its supporting parametric design plug-in Grasshopper, built a parametric design platform for impact-resistant chiral vibration isolation structures. The overall design process includes four modules: parameter adjustment, unit cell modeling, array modeling, and model storage.

[0063] S13, parameter adjustment module, is used to change the design parameters of the impact-resistant chiral vibration isolation structure. There are three key parameters, including the angle θ between the cylindrical rod and the horizontal plane, the torsion angle β of the upper and lower end planes, and the ratio γ between the rod length and the rod radius. This module uses the timer command to sequentially read the parameter combination compiled by Matlab, so as to continuously adjust the structural design parameters. At the same time, the output of this module will also sequentially modify the storage file name of the corresponding model.

[0064] S14, chiral energy lock unit cell modeling module. In the unit cell modeling module, it is used to design and generate impact-resistant chiral unit cell structures. According to the given parameter combination, the circle drawing, offset, rotation, and sine and cosine calculation commands are used to generate the trajectory lines of the center points of the upper and lower ends of the cylindrical rod respectively; the divide command is used to divide the points equally, and 20 is selected here, that is, 20 cylindrical rods are arranged in a ring, and the upper and lower endpoints are connected in sequence, and the pipe command is used to generate a solid rod structure; the upper end of the cylindrical rod is rotated to a specific angle, and it is determined whether the ring arrangement interferes.

[0065] S15, Chiral energy-locking structure array modeling module. The array modeling module is used to design and generate metal gaskets, upper and lower end interfaces of unit cells, and chiral structure arrays. The cylindrical rods arranged in an annular direction are mirrored, and the intermediate connectors and slots are designed using commands such as Boolean operations; the upper and lower end interfaces are designed, and the circumferential tubes are combined with Boolean operations. The size coefficients and scaling ratios of the metal gaskets and the upper and lower end interfaces of the unit cells are debugged to ensure that the structural coordination is reasonable and meets the feasibility of actual processing and manufacturing.

[0066] S16, model storage module. The model storage module is used to generate model file names and store models in sequence. According to the given parameter combination, set the file storage path and file name. The files can be numbered in the order of the parameters, and then baked into the Rhino space and stored separately. A variety of model storage formats can be set, such as stp, igs, etc. Four geometric models are output here, namely a single cylindrical rod, an impact-resistant chiral unit cell structure, an impact-resistant chiral unit cell structure longitudinal array, a metal gasket, etc. The output geometric model can be modified according to actual needs.

[0067] In a preferred embodiment of the present invention, the above step S2 includes:

[0068] S21, using the commercial finite element software ABAQUS, a finite element model for mechanical calculation of cylindrical rods in a single chiral energy lock unit cell is established through the CAE interface. Under axisymmetric coordinates, two constraints are set on the lower end of the cylindrical rod, one is completely fixed, and the other is only allowed to rotate in the circumferential direction, which is used to evaluate the influence of structural torsion on the static performance of the structure; the upper end of the cylindrical rod is rigidly coupled to the reference point, and is set to only have vertical displacement for displacement loading. The loading displacement is set according to the established model, and the value is twice the vertical height of the cylindrical rod. The material is selected as polymer material, the modulus is 1MPa, and the density is 1.2e-9t / mm 3 The cylindrical bar adopts solid elements, and the loading adopts the quasi-static setting in implicit dynamics. The force, torque, displacement of the upper reference point and the strain energy of the entire model are set as historical output variables.

[0069] S22, write a python script command to import a single cylindrical member model into ABAQUS and generate an inp calculation file. In the script, all objects in the Abaqus module should be imported to provide access to the default model database, mainly including importing the symbolic constant module to call the predefined constants in Abaqus; importing all objects in the caeModules module. The script defines the model as an entity and assigns material properties, sets the minimum unit size and divides the grid, sets the vertical loading displacement, and finally generates an inp calculation file. Multi-core CPU calculation can be set according to the actual computer situation. Through this python script command, the import of a single cylindrical member model can be realized through the script in the CAE interface of S21. Compared with manual modeling in the CAE interface, it can greatly improve work efficiency.

[0070] S23, for a given parameter combination, use Matlab to compile a python script file that can batch import cylindrical member models into ABAQUS and batch generate inp calculation files. Based on the python script command written in S22, using the array that stores the structural parameter numbers and the array that stores the structural parameters, the model name and parameters in the python script command are replaced by batch modification according to the text ASCII code, and a new python script command can be generated. By running the python script command in the CAE interface of S21, the cylindrical member models can be batch imported and inp calculation files can be generated. This method avoids the operation of manually modifying and running the S22 python script command separately, and greatly improves the time to generate inp calculation files for multiple cylindrical member models.

[0071] S24, for a given parameter combination, uses Matlab to generate a script file that can batch submit multiple cylindrical bar model inp calculation files, and can batch complete the finite element calculation of cylindrical bar models. At the same time, multi-core CPU calculation settings can be performed according to actual conditions.

[0072] In a preferred embodiment of the present invention, the above step S3 includes:

[0073] S31, write a python script command to extract the finite element calculation results of a single cylindrical member. In the script, all objects in the Abaqus module should be imported to provide access to the default model database, mainly including importing the symbolic constant module to call the predefined constants in Abaqus; importing all objects in the caeModules module. The python script command can open the finite element calculation result file of a single cylindrical member, output the axial force of the reference point, the torque of the reference point, the structural kinetic energy, the structural strain energy and other physical quantities and store them as a document. Compared with manually extracting results in the CAE interface, it can greatly improve work efficiency.

[0074] S32, for a given parameter combination, use Matlab to compile a python script file that can batch extract finite element calculation results. Based on the python script command written in S31, the array storing the structural parameter number and the array storing the structural parameters are used to batch modify the model name and parameters in the python script command according to the text ASCII code to generate a new python script command. By running the python script command in the CAE interface of S31, the finite element calculation results of cylindrical rods can be extracted in batches. This method avoids the operation of manually modifying and running the S31 python script command separately, greatly improving the time required for batch extraction of finite element calculation results of multiple cylindrical rods.

[0075] S33, using Matlab to compile batch finite element simulation result processing program, systematically analyze the static performance of cylindrical rods under different parameter combinations, such as Figure 4 As shown in the figure, the horizontal axis represents the calculation time step in the finite element analysis, and the vertical axis includes target physical quantities such as load, strain energy, torque, etc., which can intuitively analyze the influence of different parameter combinations on the above results; at the same time, the program distinguishes and analyzes the two situations of bottom sliding and no sliding, such as Figure 4 As shown in the comparison above and below, it can be intuitively seen that the bottom constraint has a significant impact on the static performance of the chiral energy-locking unit cell structure. At the same time, the program can analyze the influence of parameters such as the angle θ between the cylindrical rod and the horizontal plane, the torsion angle β between the upper and lower end planes, and the ratio of the rod length to the rod radius γ on the static performance of the cylindrical rod. Figure 5As shown in the figure, under a certain upper and lower end plane torsion angle β, the grayscale values ​​of each figure include physical quantities such as maximum force value, locked energy, released energy, torque, etc. The horizontal axis represents the ratio of rod length to rod radius γ, and the vertical axis represents the angle θ between the cylindrical rod and the horizontal plane. Figure 5 It is easy to see how the target physical quantity changes with the design parameters. Figure 6 The comparison results (ratios) of the maximum force, locked energy, released energy, torque and other physical quantities under two different upper and lower end plane torsion angles β with different parameter combinations can intuitively show the influence of different upper and lower end plane torsion angles β. The specific Matlab program will not be elaborated here.

[0076] In a preferred embodiment of the present invention, the above step S4 includes:

[0077] According to the size restrictions in the actual application scenarios, the overall dimensions of the impact-resistant chiral vibration isolation structure and the number of array arrangements and other parameters are clarified, and then the parameters of the required chiral energy-locking unit cell structure are set; according to the above steps S2 and S3, multi-parameter chiral energy-locking unit cell structure mechanical simulations are carried out in batches to analyze and obtain the relationship between the structural mechanical properties and the structural parameters; the maximum bearing capacity and ultimate deformation of the impact-resistant chiral vibration isolation structure are used as optimization targets, and the design scheme of the impact-resistant chiral vibration isolation structure is determined and optimized.

[0078] In summary, the above are only preferred embodiments of the present invention and are not intended to limit the protection scope of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A design method for impact-resistant chiral vibration isolation structure combining Grasshopper and Abaqus, characterized in that: The steps to implement this method include: S1: Batch design of impact-resistant chiral vibration isolation structures based on Grasshopper, generate and store models in sequence; S2: Batch read models and carry out mechanical simulation based on the Abaqus platform; S3: Batch read finite element simulation results and analyze the relationship between structural mechanical properties and structural parameters; S4: Determine and optimize the design scheme of impact-resistant chiral vibration isolation structure according to actual application requirements; The process of batch reading models and performing mechanical simulation based on the Abaqus platform in step S2 includes: S21: Using the commercial finite element software ABAQUS, a mechanical calculation finite element model for a cylindrical rod in a single chiral energy locking unit is established through the CAE interface; under axisymmetric coordinates, two constraints are set on the lower end of the cylindrical rod, one is completely fixed, and the other is only allowed to rotate in the circumferential direction, which is used to evaluate the influence of structural torsion on the static performance of the structure; the upper end of the cylindrical rod is rigidly coupled with the reference point, and is set to only have vertical displacement for displacement loading. The loading displacement is set according to the established model, and the value is twice the vertical height of the cylindrical rod; the material is selected as a polymer material, the modulus is 1MPa, and the density is 1.2e-9t / mm 3 ; The cylindrical rod adopts solid unit, and the loading adopts the quasi-static setting in implicit dynamics. The force, torque, displacement of the upper reference point and the strain energy of the whole model are set as historical output variables; S22: Write a python script command to import a single cylindrical member model into ABAQUS and generate an inp calculation file. In the script, all objects in the Abaqus module should be imported to provide access to the default model database, mainly including importing the symbolic constant module for calling the predefined constants in Abaqus; importing all objects in the caeModules module; the script defines the model as a solid and assigns material properties, sets the minimum unit size and divides the grid, sets the vertical loading displacement, and finally generates an inp calculation file, and sets multi-core CPU calculation according to the actual situation of the computer; through this python script command, the single cylindrical member model is imported through the script in the CAE interface of S21; S23: For a given parameter combination, Matlab is used to compile a python script file that batch imports cylindrical member models into ABAQUS and batch generates inp calculation files; based on the python script command written in S22, the array storing the structural parameter number and the array storing the structural parameter are used to batch modify the model name and parameters in the python script command according to the text ASCII code to generate a new python script command, and the python script command is run in the CAE interface of S21 to batch import cylindrical member models and generate inp calculation files; S24: For a given parameter combination, Matlab is used to generate a script file that can batch submit multiple cylindrical bar model inp calculation files, and can batch complete the finite element calculation of cylindrical bar models. At the same time, multi-core CPU calculation settings can be performed according to actual conditions.

2. The method for designing a shock-resistant chiral vibration isolation structure combining Grasshopper and Abaqus as claimed in claim 1, characterized in that: The step S1 comprises: S11: Determine the composition and material selection of the impact-resistant chiral vibration isolation structure; The structure is composed of a chiral energy locking unit, a metal gasket, upper and lower adapter plates, bolts and nuts, wherein the chiral energy locking unit includes a cylindrical rod, an intermediate connecting piece, a unit upper end interface and a unit lower end interface, the intermediate connecting piece is an annular hollow shell, a horizontal metal gasket slot is processed on the outer circumferential surface of the shell, and two partitions are provided at positions corresponding to the metal gasket slots on the inner circumferential surface of the shell, and the two partitions extend from the inner circumferential surface of the shell and are in the same horizontal plane. After the metal gasket is inserted into the metal gasket slot, the metal gasket and the two partitions together form a planar structure, and after the structure is impacted, the structural deformation of the unit upper end interface and the unit lower end interface is limited; The outer shapes of the upper interface of the unit and the lower interface of the unit are both cylindrical structures. The upper end of the upper interface of the unit has an external thread structure, and the lower end of the lower interface of the unit has an internal thread structure. The internal and external thread specifications of the upper interface and the lower interface are consistent; The upper interface of the unit and the lower interface of the unit are respectively connected to the upper and lower ends of the middle connecting piece through cylindrical rods, the cylindrical rods are arranged periodically along the circumferential direction of the middle connecting piece, and the cylindrical rods at the upper and lower ends of the middle connecting piece are arranged in mirror symmetry; A plurality of chiral energy locking units are sequentially connected together in the vertical direction through the connection interfaces at the upper and lower ends by threads to form a chiral energy locking unit group. The plurality of chiral energy locking unit groups are distributed in the horizontal plane and are fixedly connected through the upper and lower adapter plates and bolts and nuts to form a whole. The metal gasket is inserted into the metal gasket slot of the chiral energy locking unit. The upper end interface of the unit and the lower end interface of the unit are both made of polymer materials including PLA or ABS by 3D printing molding. S12: Using the commercial 3D NURBS industrial modeling software Rhinoceros and its supporting parametric design plug-in Grasshopper, a parametric design platform for impact-resistant chiral vibration isolation structures was built, which includes four modules: parameter adjustment, unit cell modeling, array modeling, and model storage. S13: parameter adjustment module. The function of the parameter adjustment module is to change the design parameters of the impact-resistant chiral vibration isolation structure, the angle θ between the cylindrical rod and the horizontal plane, the torsion angle β between the upper and lower end planes of the intermediate connecting member, and the ratio γ between the length of the cylindrical rod and the radius of the cylindrical rod. The parameter adjustment module uses the time counting timer command to sequentially read the parameter combination compiled by Matlab, so as to continuously adjust the structural design parameters. At the same time, the output of this module will also sequentially modify the storage file name of the corresponding model. S14: Unit cell modeling module. The unit cell modeling module is used to design and generate chiral energy locking units. According to the given parameter combination, the circle drawing, offset, rotation and sine and cosine calculation commands are used to generate the trajectory lines of the center points of the upper and lower ends of the cylindrical rod. The divide command is used to divide the points into equal parts, and the cylindrical rods are arranged in an annular direction according to the equal points. The upper and lower end points are connected in sequence using the pipe command to generate a solid rod structure; the upper end of the cylindrical rod is rotated to a specific angle, and it is determined whether the annular arrangement interferes; S15: Array modeling module, which is used to design and generate metal gaskets, unit upper end interfaces, and unit lower end interfaces; mirror the cylindrical rods arranged in the annular direction obtained above, and use Boolean operation commands to design intermediate connectors and slots; design the unit upper end interface and unit lower end interface, and perform Boolean operations on the combined unit upper end interface and the annular cylindrical rods; debug and generate the size coefficients and scaling ratios of the metal gaskets, unit upper end interfaces, and unit lower end interfaces to ensure reasonable coordination of the structure and meet the feasibility of actual processing and manufacturing; S16; model storage module, the model storage module is used to sequentially generate model file names and store models, set the file storage path and file name according to the established parameter combination, number them in sequence according to the parameter sequence, and then bake them into the Rhino space and store them separately; the model storage module outputs four geometric models in total, namely a single cylindrical rod, a chiral energy locking unit, a longitudinal array of chiral energy locking units and a metal gasket, and the output geometric model is modified according to actual needs.

3. The method for designing a shock-resistant chiral vibration isolation structure combining Grasshopper and Abaqus as claimed in claim 1, characterized in that: The process of batch reading finite element simulation results and analyzing the relationship between structural mechanical properties and structural parameters in step S3 includes: S31: Write a python script command to extract the finite element calculation results of a single cylindrical member; import all objects in the Abaqus module in the script to provide access to the default model database, including importing the symbolic constant module for calling the predefined constants in Abaqus; import all objects in the caeModules module; S32: For a given parameter combination, a python script file for batch extracting finite element calculation results is compiled using Matlab; based on the python script command written in S31, the model name and parameters in the python script command are replaced by batch modification according to the text ASCII code using the array storing the structural parameter number and the array storing the structural parameter; S33, using Matlab to compile a batch finite element simulation result processing program, systematically analyze the static performance of cylindrical rods under different parameter combinations, including load, strain energy and torque.

4. The method for designing a shock-resistant chiral vibration isolation structure combining Grasshopper and Abaqus as claimed in claim 3, characterized in that: The process of determining and optimizing the design scheme of the impact-resistant chiral vibration isolation structure in step S4 includes: clarifying the overall external dimensions and the number of array arrangements of the impact-resistant chiral vibration isolation structure according to the size restrictions in the actual application scenario, and then setting the parameters of the required chiral energy locking units; carrying out multi-parameter chiral energy locking element mechanics simulations in batches according to the above steps S2 and S3, and analyzing the relationship between the mechanical properties of the structure and the structural parameters; taking the maximum bearing capacity and ultimate deformation of the impact-resistant chiral vibration isolation structure as the optimization goals, and determining and optimizing the design scheme of the impact-resistant chiral vibration isolation structure.

Images