Distributed numerical simulation method for heat-force-vibration interaction in grinding process

Abstract

The invention relates to a heat-force-vibration interaction distributed numerical simulation method in the grinding process, and belongs to the field of numerical simulation. According to the method, a numerical simulation system server client communication architecture is established through TCP / IP, simulation starting is controlled by a server, a mode is selected firstly during starting, the number of clients required is determined according to the mode, mechanical simulation is an excitation signal, the mechanical simulation is placed on the server, and a temperature field simulation model and a vibration model are placed on the clients. And a multi-thread method is utilized to circularly wait for client link, when a link signal is received, a client address port number is recorded, a first client is set to execute temperature field simulation, and a second client is set to execute vibration simulation. According to the method, distributed simulation of force-heat-vibration interaction is realized, computing resources are increased, the computing time is shortened, and the computing efficiency is improved.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of numerical simulation, and particularly relates to a distributed numerical simulation method for thermal-force-vibration interaction in a grinding process. BACKGROUND

[0002] Grinding is an important finishing process for mechanical products, and is often applied in the fields of aerospace, automobiles and navigation, such as gears, worms, blades and crankshafts. The physical field in the grinding process is an important factor affecting the manufacturing error and surface morphology of a workpiece, and therefore, the simulation of the physical field in combination with material properties, geometric parameters, processing environment factors and machine tool dynamic characteristics is the basis for predicting the processing effect of the workpiece and realizing process adjustment and processing control.

[0003] Common existing processing process simulation software includes Abaqus, Deform, Ansys and COMSOL. The above software first establishes the energy conservation condition formula of the thermal force and vibration in the processing process, and derives the energy conservation model to form a differential equation. The constitutive equation is used to represent the nonlinear and linear stress-strain mapping relationship of the material, the three-dimensional grid is divided for the simulated object, the differential equation is combined to form an element grid dynamic characteristic representation model, and the whole dynamic characteristic representation model is assembled. The initial conditions, constraint conditions, tool-work contact conditions and load conditions are set on the grid nodes of the object, the grid reconstruction method is combined with the above model to realize the modeling of the time-varying dynamic model of the tool-work cutting domain at different times in the grinding process. Then, the optimization direction of the numerical solution of the whole field at each iteration step is established by using the residual amount method, and the numerical solution is simulated by using the least square method, Newton iteration method and Galerkin algorithm.

[0004] Compared with the patent CN104657564A (numerical simulation method for abrasive flow machining based on molecular dynamics), the grinding process simulation focusing on the microscale (atomic level) is carried out, the molecular dynamics method is used to study the change of the crystal structure of the abrasive, and the micro parameters such as atomic displacement and bond angle change are concerned. The application is aimed at the simulation of the interaction of macroscopic physical fields (thermal-force-vibration), and realizes distributed numerical calculation based on the CS architecture, and is concerned about the overall dynamic characteristics of the process system and the simulation stability under complex geometric conditions. The patent solves the macroscopic field coupling problem through the distributed architecture, avoids the grid distortion problem of the traditional finite element software under complex geometry, and CN104657564A is only applicable to atomic level modeling of microscale, and the application scene and methodology of the two are completely different.

[0005] Compared with patent CN111015370B (grinding monitoring method based on thermal coupling), real-time temperature and force signals are collected by sensors, and the processing state is judged by threshold value, which belongs to real-time monitoring technology. This patent focuses on the system architecture design of numerical simulation, supports distributed simulation of thermal-force-vibration coupling, improves computing efficiency through server-client separation model, and does not rely on real-time data acquisition. This patent proposes an open simulation architecture that allows custom mechanical and thermal models (such as replacing ANSYS modules), while CN111015370B is only a fixed threshold monitoring method that cannot be flexibly adjusted.

[0006] Compared with patent CN114239188A (cut-in grinding theoretical model simulation prediction method), the scratch morphology is predicted based on the statistical analysis of the wheel topography and the rebound effect, and the dynamic model of the interaction between single abrasive particles and workpieces is emphasized. This patent uses a distributed architecture to solve complex geometric simulation problems, avoids grid distortion, supports multi-physical field (thermal-force-vibration) coupling calculation, and is not limited to a specific grinding type. The system architecture of this patent can adapt to various machining conditions (such as shafts and irregular shapes), while CN114239188A only predicts the scratch rebound for cut-in grinding, with a relatively narrow application range.

[0007] Compared with patent CN118395714A (grinding workpiece surface micro-topography modeling method based on simulation), the surface micro-topography model is generated through coordinate transformation and silica gel replication of the wheel topography, relying on experimental data and geometric mapping. This patent realizes distributed simulation based on the CS architecture, supports thermal-force-vibration interaction calculation, and does not involve specific topography modeling algorithms, but provides flexible model replacement capabilities. This patent uses a server-client separation design to distribute computing tasks and improve efficiency, while CN118395714A focuses on the topography generation method and does not solve the grid distortion problem in complex working conditions.

[0008] Compared with patent CN118821545A (multi-pass macro-micro multi-field coupling simulation method), the grain size and dislocation density are continuously predicted for multi-pass plastic forming, integrating the JMAK model and the constitutive equation. This patent realizes multi-physical field coupling simulation through a distributed architecture, supports the embedding of self-programmed models or commercial software (such as ANSYS), and is not limited to material microstructure evolution. This patent proposes an open system architecture that allows free replacement of mechanical, thermal, and vibration models, while CN118821545A has a fixed model and focuses on material microstructure evolution, with lower flexibility.

[0009] In summary, the existing cutting process numerical simulation software is closed in terms of mechanics, heat and vibration, and researchers can only change the working conditions and develop them by changing the contact conditions, boundary conditions, load conditions, motion trajectories and material constitutive equations. However, the above finite element method is not suitable for all working conditions, and not all working conditions need to simulate all physical fields of the grinding process. For example, the machining process of blade and shaft parts has complex geometric characteristics of workpieces, so the grid used to characterize the numerical model often has quality problems. When simulating the cutting process by using the above software, negative volume grid distortion often occurs during grid reconstruction, which leads to simulation failure or large simulation error due to grid quality. Therefore, the above commercial software cannot be used to simulate all working conditions. Moreover, the cutting stability problem and the workpiece surface quality problem do not need to obtain the accurate physical quantity fluctuation of the whole cutting process, and only part of the parameters are needed. At this time, the existing commercial software is not suitable, so the present application combines geometry, modal test method, part of the simulation software function, mechanical process control theory and socket programming method to realize a set of numerical simulation method and system for complex working condition grinding with strong openness. SUMMARY

[0010] (I) Technical problems to be solved

[0011] The technical problem to be solved by the present application is to provide a grinding process heat-force-vibration interaction distributed numerical simulation method to solve the problem that the existing cutting process numerical simulation software is closed in terms of mechanics, heat and vibration, and the existing finite element method is not suitable for all working conditions and not all working conditions need to simulate all physical fields of the grinding process.

[0012] (II) Technical scheme

[0013] In order to solve the above technical problems, the present application provides a grinding process heat-force-vibration interaction distributed numerical simulation method, which comprises the following steps:

[0014] S1, a system server and a client communication system are established by a parallel communication socket method, the communication between the server and the client is realized through TCP / IP protocol, the server runs a grinding force model, the client 1 runs a temperature field simulation model, and the client 2 runs a vibration model;

[0015] S2, after the server starts simulation, the grinding force model is run to obtain the grinding force result, and the simulation time step, the grinding force and the contact position information are sent to the client 1 and the client 2 for setting to realize time synchronization;

[0016] S3, the client 1 carries a temperature field simulation model, is realized through ANSYS APDL, the heat source method model is directly integrated into the client 1 with finite element model, the temperature field result is obtained by combining the grinding force result simulation, the temperature field file is returned to the server;

[0017] S4, the client 2 carries a vibration model, is realized through programming or through ANSYS, the vibration model is directly integrated into the client 2 with finite element model, the vibration result is obtained by combining the grinding force result simulation, the vibration result file is returned to the server;

[0018] S5, after the server receives the file returned by the client 1 and the client 2, judges the running time step, if the running time step reaches the upper limit of time step, then end simulation, output grinding force result, temperature field result and vibration result, if not end then continue S2, S3, S4 step, until simulation ends.

[0019] (Three) beneficial effects

[0020] The present application proposes a kind of grinding process heat-force-vibration interaction distributed numerical simulation method, the present application has following technical effects:

[0021] 1, realize the distributed simulation of force-heat-vibration interaction, increase computing resource, shorten computing time and improve computing efficiency.

[0022] 2, using the framework, mechanics, thermal, vibration model can be replaced by self-development model, broaden the possibility of special working condition solution calculation model for researchers.

[0023] 3, using the framework, mechanics model can be replaced by based on geometry model and developed, can be used for the simulation calculation of complex curved surface machining condition, avoid the grid reconstruction attachment and problem in finite element simulation. DETAILED DESCRIPTION

[0024] Figure 1 for the numerical simulation method of heat-force-vibration interaction of the present application;

[0025] Figure 2 for the grinding force file format chart;

[0026] Figure 3 for the temperature field file format chart;

[0027] Figure 4 for the vibration quantity file format chart. DETAILED DESCRIPTION

[0028] To make the purpose, content and advantages of the present application clearer, the specific embodiments of the present application are further described in detail below in conjunction with the drawings and examples.

[0029] The application provides a grinding process thermal-force-vibration interaction distributed numerical simulation method based on a CS architecture, by which distributed numerical simulation of thermal-force coupling simulation, cutting vibration simulation and workpiece surface topography simulation of shaft and special-shaped parts can be realized.

[0030] The application provides a grinding process thermal-force-vibration interaction distributed numerical simulation method based on a CS architecture, by which distributed numerical simulation of thermal-force coupling simulation, cutting vibration simulation and workpiece surface topography simulation of shaft and special-shaped parts can be realized. Figure 1 The application provides a grinding process thermal-force-vibration interaction distributed numerical simulation method based on a CS architecture, by which distributed numerical simulation of thermal-force coupling simulation, cutting vibration simulation and workpiece surface topography simulation of shaft and special-shaped parts can be realized.

[0031] S1, a system service end and a client communication system are established through a parallel communication socket method, communication between the service end and the client is realized through a TCP / IP protocol, the service end runs a grinding force model, the client 1 runs a temperature field simulation model, and the client 2 runs a vibration model.

[0032] Since the mechanical effect in the grinding process is an excitation source generated by process system vibration and grinding temperature, and the mechanical effect is directly related to material and geometric model, the grinding process geometric model and material characteristic parameter program are written on the service end. The grinding force model can refer to formula 1 and formula 2, the model mainly considers the shear effect in the grinding process, the grinding force model can be replaced according to the research working condition requirement, and the overall grinding force can be calculated according to the given process parameters, combined with the time step and the model shown in formula 1 and formula 2, wherein,

[0033] Formula 1 is a single abrasive particle mechanical model, wherein, F n and F t are normal grinding force and tangential grinding force of abrasive particles, k c is a shear coefficient, A is a grinding cross-sectional area, p is unit area normal pressure, s is a blunt area, is a ratio of normal grinding force to tangential grinding force, mu is a friction factor, and k tc is a tangential coefficient.

[0034] Formula 2 is a complex curved surface Taylor series expansion model. Wherein, G is a tool path surface model, theta and alpha are tool path surface parameters, n x , n y , nz is the component of the three-dimensional normal vector of the workpiece surface point on the coordinate system, x, y, z are the workpiece surface points, u is the tangent vector, and the grinding force is calculated according to the above parameters.

[0035] According to the given simulation form, i.e. thermal coupling simulation or thermal-mechanical coupling simulation, the number of clients to be monitored is determined. If it is thermal coupling simulation, only one client needs to be started for temperature field simulation. If it is thermal-mechanical coupling simulation of the application, two clients need to be started for simulation. The multi-thread monitoring method is used to monitor the client address and port number accessing the server. When the number of access clients reaches the set requirement, the simulation can be started.

[0036]

[0037] S2, after the server starts the simulation, the grinding force model is run to obtain the grinding force result, and the simulation time step, grinding force and contact position information are sent to the client 1 and the client 2 for setting to realize time synchronization.

[0038] The grinding force result and the grinding position and area are formed in the format as shown in Figure 2 , wherein circle 1 is the time step, circle 2 is the grinding force, and circle 3 is the contact position. Then the information is transmitted to the client 1 and the client 2 through a binary data stream, and then the loop is started and suspended. The fixed port number information of the client is scanned and monitored. Only when the reply information corresponding to the given number of clients is received at the current time step, the next time step calculation is performed.

[0039] S3, the client 1 carries the temperature field simulation model, which is realized by ANSYS APDL, and the heat source method model and the finite element model are directly integrated into the client 1. The temperature field result is obtained by combining the grinding force result and returned to the server.

[0040] The heat flow density Q model calculated by the heat source method is shown in formula 3, wherein Q is the heat flow density, ε is the heat exchange coefficient, F t is the tangential grinding force, v is the linear velocity, and S is the contact area between the abrasive grain and the workpiece. The simulation is performed through the ANSYS batch processing mode. The client 1 first generates the APDL txt file (which stores the numerical information used for simulation and calculation), reads the finite element three-dimensional grid model of the workpiece built in advance in the APDL program, sets the transient thermal simulation, sets the heat exchange coefficient, sets the boundary temperature condition, calculates the heat flow density according to the heat source method, sets the heat flow density according to the contact position simulated in step S2, and sets the Newton iteration solver. The heat source energy on the workpiece is solved by iteration. The above is all the work involved in the program in APDL, which is automatically formed when the temperature field simulation is formed.

[0041] Then save the APDL txt file and the finite element model to a fixed path, use the DOS system to open the ANSYS batch mode to read the APDL txt file for simulation, use APDL to apply the heat source and combine the ANSYS transient thermal simulation function to simulate the temperature field. The partial differential equation of heat conduction is shown in Equation 4, where k is the thermal diffusion coefficient, U is the temperature field distribution, t is time, and V is the unit volume.

[0042] Before returning the file to the server, use the windows manager to scan whether the batch processing is completed. If the execution is completed, read the temperature field file and return it to the server to prevent returning the result file of the previous time step to the server. The file structure returned to the server is as follows: Figure 3 As shown, circle 1 is the time step, circle 2 is the temperature, circle 3 is the expansion, and circle 4 is the position.

[0043]

[0044] S4, client 2 carries the vibration model, which is implemented through programming or ANSYS, and the vibration model and the finite element model are directly integrated into client 2, and the vibration results are obtained by simulation in combination with the grinding force results, and the vibration result file is returned to the server;

[0045] If the program is implemented, the finite element unit stiffness and mass matrices K and M, the linear damping matrix C, and the assembly method program are integrated into Client 2. The spindle model is characterized by diameter and length, and the K and M matrices are solved. The K and M matrices can be derived from the Timoshenko beam. The assembly method program merges the K and M matrices into the K and M matrices of the overall model based on the nodes. The NewMark algorithm and modal shape method are implemented in Client 2 to solve the vibration curve during the output process. Then, the next step is to suspend before receiving the mechanical simulation signal from the server and continuously monitor the server information. The simulation starts only when the server receives the grinding force information, solves the grinding force application position on the spindle, reads the finite element model, applies the grinding force to solve the vibration, generates the vibration result file, and returns it to the server. The simulation is implemented in Client 2. When the simulation is complete, it resumes suspending and monitoring the server information.

[0046] If the simulation is carried out through the ANSYS batch mode, the client 2 firstly generates the APDL txt file, reads the finite element three-dimensional model of the grinding wheel-spindle system realized modeling by using the APDL txt file, and sets to the transient dynamics simulation mode, then according to the reasonable grinding wheel position simulated in step S2, loads the load to the grinding force applying position by using the APDL program, sets the NewMark method solver, solves the grinding force applying position on the spindle, all of the above are the work involved in APDL, and are automatically formed by the program when the APDL file is generated. Then the APDL txt file and the grinding process simulation model are saved to the fixed path, the ANSYS batch mode is opened by using the DOS system to read the APDL txt file for simulation, the grinding force is applied by using the APDL to solve the vibration, and the vibration curve simulation is realized. Before returning the file to the server, it is scanned by using the windows manager whether the batch processing is executed to end, the vibration result file is read if the execution ends, and the server is returned, so as to prevent the result file of the previous time step from being returned to the server. The file structure returned to the server is as shown in Figure 4

[0047] S5, after the server receives the files returned by the client 1 and the client 2, the running time step is judged, if the running time step reaches the upper limit of the time step, the simulation is ended, and the grinding force result, the temperature field result and the vibration result are output, if it is not ended, steps S2, S3 and S4 are continuously carried out until the simulation is ended.

[0048] The key points of the application

[0049] 1, the application proposes a thermal-force-vibration interaction numerical simulation method based on the CS architecture, and is protected for the thermal-force-vibration interaction as shown in the framework. Figure 1 The numerical simulation system service client communication architecture is established by using TCP / IP, the service server controls the simulation to start, when starting, the mode is selected first, the number of required clients is determined according to the mode, the mechanical simulation is placed in the service server, and the temperature field simulation model and the vibration model are placed in the client. The multi-thread method is used to circularly wait for the client link, when the link signal is received, the client address port number is recorded, the first client is set to execute the temperature field simulation, and the second client is set to execute the vibration simulation.

[0050] ​2、When using finite element software to simulate the client, the order of the process is controlled by the windows manager scanning method. Since the simulation results of the time step and the simulation results of the previous time step are controlled by the APDL file, the file name and storage path are usually fixed, and it is difficult to distinguish when the simulation ends. In order to realize the synchronization of time step, the windows manager is used to monitor whether the ANSYS batch processor is closed, and the batch processor is closed as the end signal to confirm the end of simulation, the simulation results are updated, and the next simulation time step is triggered.

[0051] Effects of the present application

[0052] 1、Realize the distributed simulation of force-heat-vibration interaction, increase the computing resources, shorten the computing time and improve the computing efficiency.

[0053] 2、Using the framework, the mechanical, thermal and vibration models can be replaced with self-developed models, which can broaden the possibility of researchers to study special working conditions and solve calculation models.

[0054] 3、Using the framework, the mechanical model can be replaced with the development based on the geometric model, which can be used for simulation calculation of complex curved surface machining conditions, avoiding the grid reconstruction attachment and problems in finite element simulation.

[0055] The above only describes the preferred embodiments of the present application, and it should be pointed out that for ordinary skilled in the art, without departing from the technical principles of the present application, a number of improvements and modifications can be made, and these improvements and modifications should be considered as the protection scope of the present application.

Claims

1. A distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process, characterized in that: The method comprises the following steps: S1. Establish a system server and client communication system through a parallel communication socket method, and realize communication between the server and the client through the TCP / IP protocol. The server runs the grinding force model, the client 1 runs the temperature field simulation model, and the client 2 runs the vibration model; S2. After the server starts the simulation, it runs the grinding force model to obtain the grinding force results, and sends the simulation time step, grinding force, and contact position information to Client 1 and Client 2 for setting and achieving time synchronization; S3, client 1 carries the temperature field simulation model, which is implemented through ANSYS APDL. The heat source method model and the finite element model are directly integrated into client 1. The temperature field results are obtained by combining the grinding force results simulation, and the temperature field file is returned to the server; S4, client 2 carries the vibration model, which is implemented through programming or ANSYS, and the vibration model and the finite element model are directly integrated into client 2, and the vibration results are obtained by simulation in combination with the grinding force results, and the vibration result file is returned to the server; S5. After the server receives the files returned by client 1 and client 2, it determines the running time step. If the running time step reaches the time step upper limit, the simulation ends and the grinding force results, temperature field results and vibration results are output. If it is not ended, it continues with steps S2, S3, and S4 until the simulation ends.

2. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 1, characterized in that: In the S1, The grinding force model is Equation 1 and Equation 2. This model considers the shear effect during the grinding process and calculates the overall grinding force based on the given process parameters, the time step, and the model of Equation 1 and Equation 2. Formula 1 is the mechanical model of a single abrasive particle, where F n and F t is the normal grinding force and tangential grinding force of the abrasive, k c is the shear coefficient, A is the grinding cross-sectional area, p is the positive pressure per unit area, s is the grinding blunt area, is the ratio of normal grinding force to tangential grinding force, μ is the friction factor, k tc Tangential coefficient; Formula 2 is the Taylor series expansion model of complex surfaces, where G is the tool path surface model, θ and α are the tool path surface parameters, and n x 、n y 、n z is the component of the three-dimensional normal vector of a point on the workpiece surface in the coordinate system, x, y, z are the workpiece surface points, and u is the tangent vector. The grinding force is calculated based on the above parameters; 3. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 1, characterized in that: The multi-threaded monitoring method is used to monitor the client addresses and port numbers accessing the server, and the simulation is started when the number of accessing clients reaches the set requirements.

4. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 1, characterized in that: In S2, the simulation time step, grinding force, and contact position information are transmitted to client 1 and client 2 respectively through binary data streams, and then a loop is started to suspend and scan the fixed port number information of the listening client. The next time step calculation is performed only when and only when the reply information of the given number of clients is received in the current time step.

5. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 1, characterized in that: In S3, the heat source method is used to calculate the heat flux density Q model as shown in Formula 3: Where Q is the heat flux, ε is the heat exchange coefficient, F t is the tangential grinding force, v is the linear velocity, S is the contact area between the abrasive and the workpiece, and the simulation is performed through the ANSYS batch mode. The client 1 first generates an APDL txt file to store numerical information used for simulation and calculation, reads the pre-modeled finite element three-dimensional mesh model of the workpiece in the APDL program, sets the transient thermal simulation, sets the heat exchange coefficient, sets the boundary temperature conditions, calculates the heat flux density according to the heat source method, realizes the setting of the heat flux density according to the contact position simulated in step S2, and sets it to a Newton iterative solver to solve the heat source energy on the workpiece through iteration. The above are all the work involved in the program in APDL, which is automatically formed when the temperature field simulation is formed.

6. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 5, characterized in that: In S3, the APDL txt file and the finite element model are saved to a fixed path, and the APDL txt file is read into the ANSYS batch mode using the DOS system for simulation. The APDL is used to apply a heat source and the temperature field is simulated in combination with the ANSYS transient thermal simulation function. The partial differential equation of heat conduction is shown in Equation 4. Where k is the thermal diffusivity, U is the temperature field distribution, t is the time, and V is the unit volume; Before returning the file to the server, use the windows manager to scan whether the batch processing has been completed. If the execution is completed, read the temperature field file and return it to the server to prevent returning the result file of the previous time step to the server.

7. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 6, characterized in that: In S3, the temperature field file information returned by the client 1 to the server includes: time step, temperature, expansion amount and position.

8. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 1, characterized in that: In the S4, if the programming is implemented, the finite element unit stiffness mass matrix K, M, linear damping matrix C and assembly method program are integrated into client 2, the spindle model is characterized by diameter length, and the K matrix and M matrix are solved. The K matrix and M matrix are derived according to the Timoshenko beam. The assembly method program merges the K and M matrices into the K and M matrices of the overall model according to the nodes, and implements the NewMark algorithm\modal vibration method in client 2 to solve the vibration curve in the output processing, and then proceeds to the next step, suspends before receiving the mechanical simulation signal from the server, and continuously monitors the server information in a loop. The simulation is started when and only when the grinding force information from the server is received, the grinding force application position on the spindle is solved, the finite element model is read, the grinding force is applied to solve the vibration, the vibration result file is generated and returned to the server, and the simulation work is implemented in client 2; when the simulation is completed, it continues to suspend and monitor the server information.

9. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 1, characterized in that: In said S4, if the simulation is performed through the ANSYS batch mode, the client 2 first generates an APDL txt file, uses the APDL txt file to read the modeled grinding wheel-spindle system finite element three-dimensional model, and sets it as a transient dynamic simulation model, and then uses the APDL program to load the load to the grinding force application position based on the reasonable grinding wheel position simulated in step S2, sets it to the NewMark method solver, and solves the grinding force application position on the spindle. The above are all APDL-related work, which are automatically generated by the program when the APDL file is generated; the APDL txt file and the grinding process simulation model are saved to a fixed path, and the DOS system is used to open the ANSYS batch mode to read the APDL txt file for simulation, and the grinding force is applied using APDL to solve the vibration and realize vibration curve simulation; before returning the file to the server, the windows manager is used to scan whether the batch execution is completed. If the execution is completed, the vibration result file is read and returned to the server to prevent the result file of the previous time step from being returned to the server.

10. The distributed numerical simulation method for thermal-mechanical-vibration interaction in a grinding process according to claim 8 or 9, characterized in that: In S4, the information of the vibration result file returned by the client 2 to the server includes: time step and vibration amount.

Images