CN109063248B - Gas-liquid-solid coupling calculation method for laser shock peening - Google Patents
Gas-liquid-solid coupling calculation method for laser shock peening Download PDFInfo
- Publication number
- CN109063248B CN109063248B CN201810672217.5A CN201810672217A CN109063248B CN 109063248 B CN109063248 B CN 109063248B CN 201810672217 A CN201810672217 A CN 201810672217A CN 109063248 B CN109063248 B CN 109063248B
- Authority
- CN
- China
- Prior art keywords
- fluid
- calculation
- solving
- model
- transient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Laser Beam Processing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a gas-liquid-solid coupling calculation method for laser shock peening, which mainly comprises the following steps: designing a transient gas-liquid-solid coupling system; geometric modeling and grid division; performing simulation pretreatment by FLUENT; solving the settings with a transient structure; and a gas-liquid-solid system is arranged in a coupling way. Aiming at various defects in the prior art, the invention provides a simulation method based on laser shock peening, and provides a gas-liquid-solid coupling simulation model, wherein the model comprises a two-phase fluid and a solid phase, and a simulation system capable of simultaneously performing transient gas-liquid-solid coupling simulation is required, so that ANSYS is used for performing the gas-liquid-solid coupling simulation. The method can effectively solve the problem of difficult optimization of laser shock peening parameters.
Description
Technical Field
The invention relates to the technical field of laser shock peening, in particular to a gas-liquid-solid coupling calculation method for laser shock peening.
Background
The laser shock strengthening technology is a new material surface strengthening means, it acts on the absorption coating coated on the surface of the metal target material through the laser beam with high power density (GW/cm 2 magnitude) and pulse width (ns magnitude) through the transparent constraint layer, the coating material absorbs the laser energy and gasifies rapidly to form high-temperature and high-pressure plasma, the plasma generates high-intensity shock wave to the metal surface under the action of the constraint layer, the shock wave acts on the material surface, and the material surface is strengthened. When the peak pressure of the shock wave exceeds the dynamic yield strength of the material, plastic strain is generated on the surface layer of the material, residual compressive stress with a certain depth is generated in the material due to the reaction of the material in the impact region when the laser action is finished, and the existence of the residual compressive stress causes the closing effect of cracks, so that the driving force for fatigue crack propagation is effectively reduced, and the service life of the part is prolonged.
The three-dimensional flat-topped Gaussian beam is a laser with a uniform flat-topped area in light energy spatial distribution, and the impact strengthening process is influenced by various variable factors due to a complex mechanism, so that great difficulty is brought to the optimization of laser impact strengthening process parameters. The method of multiple attempts, relying solely on experimental data and operating experience, requires a significant amount of time and money. In turn, finite element modeling methods are used to assist in the selection of laser shock peening process parameters, while accounting for the peening mechanism by analyzing changes in stress strain and displacement. In the aspect of finite element simulation, brasted and Brockman firstly adopt an ABAQUS/Explicit + ABAQUS/Implicit method to simulate laser shock strengthening, later people mostly adopt the finite element simulation method, but for laser shock strengthening, the method is time-consuming and needs to establish a plurality of analysis models aiming at different process parameters (spot radius, lap joint rate, shock route and the like), so that a gas-liquid-solid coupling simulation method for laser shock strengthening is urgently needed to simulate and analyze the laser shock strengthening.
At present, the prior art has the following defects: the research of the prior experimental technology on laser-induced forward transfer needs a lot of time and cost, when the finite element simulation method is used for laser shock strengthening, the time is consumed, the result of explicit analysis of each light spot needs to be continuously brought into implicit analysis, and multiple analysis models need to be established according to different process parameters (such as the radius of the light spot, the lap joint rate and the shock route).
Accordingly, further improvements and improvements are needed in the art.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a gas-liquid-solid coupling calculation method for laser shock peening.
The purpose of the invention is realized by the following technical scheme:
a gas-liquid-solid coupling calculation method for laser shock peening mainly comprises the following steps:
step S1: designing a transient gas-liquid-solid coupling system; the method mainly comprises the steps of firstly establishing a structural geometric model and a fluid geometric model in geometric modeling, then carrying out grid division, carrying out simulation pretreatment on a fluid module, carrying out solving setting on the structural module, then carrying out coupling mode setting on a system coupling module, and finally solving and post-processing a simulation result.
The modeling step of the transient gas-liquid-solid coupling system in the step S1 is as follows:
step S11: establishing a geometric model in a Geometry modeling module according to the distribution of fluid and a Structural region in the laser shock strengthening process, defining a boundary name, allocating the geometric model of the fluid region to a FLUENT fluid calculation simulation module for Transient fluid modeling, and allocating the geometric model of the Structural region to a Transient Structural Transient structure calculation simulation module for Transient structure modeling;
step S12: respectively dividing a fluid grid and a Structural grid in a FLUENT simulation module and a Transient Structural simulation module;
step S13: performing simulation pre-processing setting in a FLUENT fluid calculation simulation module;
step S14: carrying out solving setting in a Transient Structural Transient structure calculation simulation module;
step S15: setting a Coupling mode in a System Coupling module;
step S16: and storing and performing solution and simulation result post-processing.
Step S2: geometric modeling and grid division; the method mainly comprises the steps of firstly establishing the thickness of a three-dimensional model into a thin sheet with the size of a grid, setting symmetrical boundary conditions to reduce calculated amount, secondly naming boundaries and setting the boundary conditions, and finally adjusting a structure module to enable the grid of the structure module to be denser along with the gradual distance from a fluid-solid interface.
The geometric model establishing and mesh dividing in the step S2 specifically comprises the following steps:
step S21: establishing the thickness of the geometric model as a sheet with the thickness of one grid size, setting the front surface and the rear surface of the fluid area as symmetrical boundary conditions, dividing the upper half part of the fluid area into a FLUENT fluid area, and dividing the lower half part of the fluid area into a Transient Structural structure area;
step S22: after the geometric model is established, respectively allocating the geometric model of the fluid region and the geometric model of the Structural region to a FLUENT fluid calculation simulation module and a Transient Structural Transient structure calculation simulation module, and then carrying out mesh division in the respective modules;
step S23: in the fluid simulation module, a fluid grid is divided into a cubic structured grid, the grid size ensures smaller calculation amount and higher solving precision, and the grid size is set to be 10 mu m;
step S24: in the structural simulation module, the area closer to the fluid-solid interface needs higher solving precision, so the grid needs to be divided into more fine areas, and the area far away from the fluid-solid interface does not need too high grid density to save calculation amount, so the grid size is gradually increased.
And step S3: performing simulation pretreatment by FLUENT; the method mainly comprises the steps of defining a multiphase flow model, an energy model, a turbulence model, material attributes, setting boundary conditions, solving methods and control, and initially calculating and setting.
The simulation preprocessing in the step S3 includes the following steps:
step S31: defining a multi-phase flow model, wherein the fluid model comprises two-phase fluid of plasma bubbles and a constrained layer;
step S32: an energy model, which needs to be set on due to the compressible flow of the high pressure plasma bubble gas phase involved in the fluid model;
step S33: a turbulence model: in a flow field for generating plasma by laser induction, the Reynolds number of the flow field is calculated by the following formula:
wherein, the first and the second end of the pipe are connected with each other,
P 0 ＝P R -P ∞ in the formula L 0 Represents a length; u. u 0 Represents the speed; p is ∞ Representing the pressure of the environment; p is R Represents the internal pressure of the bubble; ρ represents the density of the fluid; v represents the kinematic viscosity of the fluid;
step S34: material attribute, setting the two fluid phases as a first phase which is a plasma bubble phase water-vapor, and the material attribute is ideal gas; the second phase is water-liquid or air, and the material properties are water and air respectively;
step S35: setting the boundary condition of a flow field according to the actual situation of the laser shock peening process, and setting the pressure outlet as the pressure outlet boundary condition by combining the named boundary name in geometric modeling;
step S36: solving method and control, calculating and simulating the compressible multiphase flow laser-induced plasma bubble transient flow field;
step S37: initializing, wherein the global initialization of a flow field is adopted during initialization, and the Region where the plasma bubble is located is marked by Region attachment of FLUENT self-contained;
step S38: calculating settings: setting the calculation time step to 1 × 10 -8 And s, the number of calculation steps is 100, data is stored once in each step, and the calculation time step length of each step can be increased and the number of calculation steps can be reduced appropriately to reduce the calculation amount when water constraint exists.
And step S4: solving the settings with the transient structure; the method mainly comprises the steps of defining materials, solving time, applying load, applying fixation and adding solving quantity.
The concrete steps of solving the transient structure in the step S4 are as follows:
step S41: defining materials, selecting steel in a material library as the material of a workpiece in the model, and keeping the material attribute at a default value;
step S42: solving time, namely matching the solving step length and the solving time of the fluid-solid module according to the time step length and the calculation time set defined in the fluid model;
step S43: applying a load, wherein a surface named wall-solid in geometric modeling is defined as an application surface of the load, and an application mode is defined as a fluid-solid interface so that the load of a flow field can be applied to the surface;
step S44: applying fixation, and defining a surface named fixed-support in geometric modeling as a surface of a fixed support to constrain the geometric body of the workpiece;
step S45: adding solving quantity, inserting response of deformation, strain and stress needing to be solved.
Step S5: a gas-liquid-solid system coupling device; the method mainly comprises the steps of setting solving time, data transmission and transmission sequence and calculating and solving.
The gas-liquid-solid system coupling setting in the step S5 includes the following specific steps:
step S51: solving time, matching the time step length and the calculation time defined in the fluid model and the structural model, and enabling the solving step length and the solving time of the system coupling module to be consistent with those of the other two modules;
step S52: data transmission, namely adding two surfaces named wall-fluid and wall-solid in geometric modeling as a data transmission surface, wherein the transmitted physical attribute is force;
step S53: the data transmission sequence is that the simulation result data of the transient fluid is used as the boundary condition of the transient structure calculation simulation module to be loaded, and then the calculation result of the transient structure is used as the boundary condition of the next calculation time step fluid simulation module to be loaded;
step S54: and (4) calculating and solving, namely setting data output as all time steps, and then saving the model to carry out calculation and solving.
Compared with the prior art, the invention also has the following advantages:
(1) The gas-liquid-solid coupling calculation method for laser shock peening, provided by the invention, has the advantages of simplicity in operation and easiness in implementation, and can effectively solve the problem of difficulty in optimization of laser shock peening parameters.
(2) According to the gas-liquid-solid coupling calculation method for laser shock peening, the thickness of the geometric model is established into a sheet with the thickness of one grid size, and the front surface and the rear surface of the fluid area are set as symmetrical boundary conditions, so that a three-dimensional model with small calculation can be obtained, and the calculation requirements of a small calculated amount and a simulation module are met.
Drawings
Fig. 1 is a flow chart of modeling a transient coupling system according to the present invention.
Fig. 2 is a schematic diagram of a fluid-solid coupling geometric model provided by the present invention.
FIG. 3 is a diagram illustrating simulation results with a water confining layer according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer and clearer, the present invention is further described below with reference to the accompanying drawings and examples.
Example 1:
as shown in fig. 1, 2 and 3, the present embodiment discloses a gas-liquid-solid coupling calculation method for laser shock peening. Firstly, designing a transient gas-liquid-solid coupling system; secondly, geometric modeling and grid division; then performing simulation pretreatment by FLUENT; setting a transient structure solution; finally, a gas-liquid-solid system is coupled. The simulation method mainly comprises 5 modules, and the specific modules are as follows:
1. transient gas-liquid-solid coupling system
Fig. 1 is a schematic diagram showing a modeling process of an ANSYS transient gas-liquid-solid coupling system and a data exchange mode thereof. The modeling process of the transient gas-liquid-solid coupling system comprises the following steps:
1) Establishing a geometric model in a Geometry modeling module according to the distribution of fluid and a Structural region in the laser shock strengthening process, defining a boundary name so as to play a role in marking when setting fluid and Structural boundary conditions and setting system coupling data transmission, allocating the geometric model of the fluid region to a FLUENT fluid calculation simulation module for Transient fluid modeling, and allocating the geometric model of the Structural region to a Transient Structural Transient structure calculation simulation module for Transient structure modeling;
2) Respectively dividing a fluid grid and a Structural grid in a FLUENT simulation module and a Transient Structural simulation module;
3) Carrying out simulation pretreatment setting in a FLUENT fluid calculation simulation module, sequentially defining a multiphase flow model, an energy model and a turbulence model, defining attributes of plasma bubbles and a constrained layer substance generated by laser induction, defining phase and boundary conditions, defining a solving method, solving control and initialization, and defining a calculation time step length, a solving step number and a result storage step length;
4) The method comprises the steps of carrying out solving setting in a Transient Structural Transient structure calculation simulation module, sequentially setting Structural materials and material attributes thereof, defining solving time step length and solving time, defining load boundary conditions and support fixed boundary conditions, and defining deformation, strain, stress and the like to be solved.
5) Setting a Coupling mode in a System Coupling module, sequentially defining solving time step length and solving time, defining a Coupling data transmission and loading boundary surface between a transient fluid calculation result and a transient structure calculation result, defining a Coupling data transmission and loading sequence, and defining a result storage step length.
6) And storing and performing solution and simulation result post-processing.
2. Geometric modeling and meshing
The Transient Structural Transient structure calculation simulation module in ANSYS can only solve three-dimensional problems, and correspondingly, the FLUENT module is also set to be three-dimensional solving, so that the calculated amount of the model is relatively large, and in order to reduce the calculated amount as much as possible, the thought of geometric modeling is as follows: the thickness of the geometric model is established into a sheet with the thickness of one grid size, and the front surface and the rear surface of the fluid area are set as symmetrical boundary conditions, so that a three-dimensional model with small calculation can be obtained, the calculation result is equal to a two-dimensional model, and the calculation requirements of a small calculation amount and a simulation module are met.
The geometric model of the laser shock peening fluid-solid coupling is shown in fig. 2, the upper half part is a FLUENT fluid region, and the lower half part is a Transient Structural region. In order to set boundary conditions and data transmission interfaces in the following modeling setting, three surfaces of the upper surface, the left surface and the right surface of the fluid region are named as pressure output for setting the fluid region as pressure outlet boundary conditions, the front surface and the rear surface are respectively named as symmetry-forward and symmetry-backward for setting the fluid region as symmetrical boundary conditions, and the surface adjacent to the structural region is named as wall-fluid for setting the fluid region as a wall surface and a coupling data exchange surface; the face above the structural area and adjacent to the fluid area is named wall-solid for setting it as a load-bearing face and a coupled data exchange face, and below named fixed-support for setting it as a structural fixed support boundary condition.
After the geometric model is built, the geometric model of the fluid region and the geometric model of the Structural region are respectively distributed to a FLUENT fluid calculation simulation module and a Transient Structural Transient structure calculation simulation module, and then grid division is carried out in the respective modules.
In the fluid simulation module, FLUENT supports structured and unstructured grids at the same time, in order to simplify the grid division process, the fluid grid is divided into a cubic structured grid, the grid size ensures less calculation amount and higher solving precision, and the grid size is set to be 10 μm.
In the structural simulation module, a region closer to a fluid-solid interface needs higher solving precision, so that a grid needs to be divided into thinner regions, and a region far away from the fluid-solid interface does not need too high grid density to save calculation amount, so that the grid size is gradually increased, specifically, the grid size of a wall-solid surface is defined to be 10 μm corresponding to the size of a fluid structured grid, the grid division mode of a geometric model of the whole structural region is defined to be a sweep, and the type of the grid is a whole square. In order to make the mesh grow larger as it becomes farther from the fluid-solid interface, the growth rate of the mesh is defined to be 1.20, the minimum mesh size is 10 μm, and the maximum mesh size is 100 μm.
3. FLUENT simulation preprocessing
1) Defining a multiphase flow model: the Fluid model contains two-phase Fluid of plasma bubbles and a constrained layer, so the model adopts a Volume of Fluid (VOF) multi-phase flow model. The VOF is a typical Euler multiphase flow model, the tracking of a moving multiphase flow interface is realized by solving the volume fraction of a target fluid in a unit grid, and when the volume of the target fluid in the unit grid is more than zero and less than one, the unit grid is the position of the fluid interface;
2) An energy model: due to the compressible flow of the high pressure plasma bubble gas phase involved in the flow model, its energy model needs to be set on;
3) A turbulence model: in a flow field for generating plasma by laser induction, the Reynolds number of the flow field is calculated by the formula:
P 0 ＝P R -P ∞ (3)
in the formula L 0 A length (m); u. u 0 Velocity (m/s); p is ∞ The pressure (Pa) of the environment; p is R The internal pressure (Pa) of the bubbles; density of rho fluid (kg/m) 3 ) (ii) a V kinematic viscosity of fluid (m) 2 In s). The Reynolds number obtained by calculation by using the formulas (1), (2) and (3) is smaller than the critical Reynolds number, so that the flow field is laminar, and in addition, the model only takes the flow field condition within 1 mu s into consideration, and sets the turbulent flow model as laminar flow;
4) The material property is as follows: in order to compare laser shock waves with a water constraint layer and a non-water constraint layer and the action of the laser shock waves on a workpiece, the two fluid phases are set as follows: the first phase is a plasma bubble phase water-vapor, with the material properties of an ideal gas. The second phase is water-liquid (in the case of a water-binding layer) or air (in the case of no water-binding layer), and the material properties are water and air, respectively;
5) Boundary conditions: setting boundary conditions of a flow field of the laser shock peening process according to the actual situation of the laser shock peening process, setting pressure outlets as boundary conditions of pressure outlets by combining with named boundary names during geometric modeling, setting the value of the pressure outlets to be gage pressure equal to 0 as the operation condition of a model is atmospheric pressure, setting wall-fluid as a wall surface, and setting symmetry-forward and symmetry-backward as symmetrical boundary conditions;
6) Solving method and control: for the calculation simulation of the compressible multiphase flow laser-induced plasma bubble transient flow field, the pressure-speed coupling system to be adopted is PISO, and the pressure difference format is PRESTO! The format, density and momentum equation discrete format all use Second Order Upwind to achieve higher accuracy. During solving control, the sub-relaxation factor setting values of all the physical quantities are obtained by practice of predecessors, the solving precision and the calculation convergence are comprehensively considered, so that the default values are taken as the criteria, the sub-relaxation factor of the pressure item is 0.3, the sub-relaxation factor of the momentum item is 0.7, and the sub-relaxation factor of the density, the volume force and the energy item is 1;
7) Initialization: during initialization, global initialization of a flow field is adopted, the Region where the plasma bubble is located is marked by Region attachment of FLUENT, the Region is a semicircle with the diameter of 100 mu m, the circle center is at the center of the bottom of the flow field, the Region pressure is set as the pressure of the plasma bubble, the phase is set as a first phase, namely a plasma bubble phase, and the rest regions are set as a second phase, namely a water phase or an air phase. Since the yield strength of a metal material is generally several hundred MPa or more, the pressure is set to 0.1X 10 in the case of a water-restraining layer and a non-water-restraining layer, respectively 9 Pa for comparing the two cases of shock waves, corresponding to the material, and then setting its pressure 2X 10 in the presence of a water-binding layer 9 Pa for mouldsThe impact of a laser shock wave on the order of GPa in the flow field and on the workpiece;
8) Calculating settings: setting the calculation time step to 1 × 10 -8 And s, the number of calculation steps is 100, data is stored once in each step, and the calculation time step length of each step can be increased and the number of calculation steps can be reduced appropriately to reduce the calculation amount when water constraint exists.
4. Transent Structural solution settings
The Structural modeling of the Transient Structural module is simpler than the modeling of the fluid, and the process of Structural modeling is also explained here:
1) Definition of materials: selecting a steel in a material library as a material of a workpiece in the model, and keeping the material attribute as default;
2) Solving time: matching the solving step length and the solving time of the fluid-solid module according to the time step length and the calculation time set defined in the fluid model;
3) Applying a load: defining a surface named wall-solid in geometric modeling as a load application surface, and defining an application mode as a fluid-solid interface so that the load of a flow field can be applied to the surface;
4) Applying and fixing: defining a face named fixed-support in geometric modeling as a face of a fixed support to constrain the geometry of the workpiece;
5) Adding a solution amount: insertion deformation, strain and stress, etc. are required to solve the response.
5. System Coupling settings
The system coupling setting is the last step of fluid-solid coupling modeling setting, and specifically comprises the following steps:
1) Solving time: matching time step lengths and calculation time defined in the fluid model and the structural model to enable the solving step lengths and solving time of the system coupling module and the other two modules to be consistent;
2) Data transmission: two surfaces named wall-fluid and wall-solid in geometric modeling are added as data transmission surfaces, and the transmitted physical attribute is force;
3) Data transfer sequence: loading simulation result data of the transient fluid as a boundary condition of a transient structure calculation simulation module, and loading a calculation result of the transient structure as a boundary condition of a next calculation time step fluid simulation module;
4) And (3) calculating and solving: and setting the data output as all time steps, and then saving the model to calculate and solve.
Example 2:
in this example, a shock wave with a smaller size is used as a wave source to investigate the influence of high-pressure LIP bubbles on the structure in the presence of a water confinement layer. Using an initial diameter of 100 μm and a pressure of 0.1X 10 9 A semicircle ideal gas bubble near the wall surface of Pa is used as a medium for generating pressure, and simulation is carried out under the condition that a water constraint layer exists. In order to more intuitively display the relationship between the bubble pressure, the bubble volume and the stress response of the structure, the left half of the fluid region in the result diagram at each moment is displayed as a pressure cloud diagram, the right half is displayed as a bubble phase diagram (the darkest color in the middle is an LIP bubble phase, the darker color in the periphery is an environment phase), and the structure region is displayed as a stress cloud diagram. Fig. 3 shows simulation results from 40ns to 400ns with a water confining layer.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents thereof, and all such changes, modifications, substitutions, combinations, and simplifications are intended to be included in the scope of the present invention.
Claims (6)
1. A gas-liquid-solid coupling calculation method for laser shock peening is characterized by comprising the following steps:
step S1: designing a transient gas-liquid-solid coupling system; firstly, establishing a structural geometric model and a fluid geometric model in geometric modeling, then carrying out grid division, carrying out simulation pretreatment on a fluid module, carrying out solving setting on the structural module, carrying out coupling mode setting on a system coupling module, and finally solving and post-processing a simulation result;
step S2: geometric modeling and grid division; firstly, establishing the thickness of a three-dimensional model as a sheet with the size of a grid, setting symmetrical boundary conditions to reduce the calculated amount, secondly naming the boundary and setting the boundary conditions, and finally adjusting a structural module to ensure that the grid is denser along with the gradual distance from a fluid-solid interface;
and step S3: performing simulation pretreatment by FLUENT; the method comprises the steps of defining a multi-phase flow model, an energy model, a turbulence model, material attributes, setting boundary conditions, solving methods and control, and setting initial calculation;
and step S4: solving the settings with a transient structure; defining materials, solving time, applying load, applying fixation and adding solving quantity;
step S5: a gas-liquid-solid system coupling device; the method comprises the steps of setting solving time, data transmission and transmission sequence and calculating and solving.
2. The gas-liquid-solid coupling calculation method for laser shock peening according to claim 1, wherein the modeling step of the transient gas-liquid-solid coupling system in step S1 is as follows:
step S11: establishing a geometric model in a Geometry modeling module according to the distribution of fluid and a Structural region in the laser shock strengthening process, defining a boundary name, allocating the geometric model of the fluid region to a FLUENT fluid calculation simulation module for Transient fluid modeling, and allocating the geometric model of the Structural region to a Transient Structural Transient structure calculation simulation module for Transient structure modeling;
step S12: respectively dividing a fluid grid and a Structural grid in a FLUENT simulation module and a Transient Structural simulation module;
step S13: performing simulation pre-processing setting in a FLUENT fluid calculation simulation module;
step S14: carrying out solving setting in a Transient Structural Transient structure calculation simulation module;
step S15: setting a Coupling mode in a System Coupling module;
step S16: and storing and performing solution and simulation result post-processing.
3. The gas-liquid-solid coupling calculation method for laser shock peening according to claim 1, wherein the geometrical model building and meshing in step S2 are specifically performed by the following steps:
step S21: establishing the thickness of the geometric model as a sheet with the thickness of one grid size, setting the front surface and the rear surface of the fluid area as symmetrical boundary conditions, dividing the upper half part of the fluid area into a FLUENT fluid area, and dividing the lower half part of the fluid area into a Transient Structural structure area;
step S22: after the geometric model is established, respectively allocating the geometric model of the fluid region and the geometric model of the Structural region to a FLUENT fluid calculation simulation module and a Transient Structural Transient structure calculation simulation module, and then carrying out grid division in the respective modules;
step S23: in the fluid simulation module, a fluid grid is divided into a cubic structured grid, the grid size ensures smaller calculation amount and higher solving precision, and the grid size is set to be 10 mu m;
step S24: in the structural simulation module, the area closer to the fluid-solid interface needs higher solving precision, so the grid needs to be divided into more fine areas, and the area far away from the fluid-solid interface does not need too high grid density to save calculation amount, so the grid size is gradually increased.
4. The gas-liquid-solid coupling calculation method for laser shock peening according to claim 1, wherein the pre-simulation processing in step S3 includes the steps of:
step S31: defining a multi-phase flow model, wherein the fluid model comprises two-phase fluid of plasma bubbles and a constrained layer;
step S32: an energy model, which needs to be set on due to the compressible flow of the high pressure plasma bubble gas phase involved in the fluid model;
step S33: a turbulence model: in laser-induced plasma generationIn the daughter flow field, the reynolds number of the flow field is calculated by the formula:
wherein the content of the first and second substances,
P 0 ＝P R -P ∞ in the formula L 0 Represents a length; u. u 0 Represents a speed; p ∞ Representing the pressure of the environment; p R Represents the internal pressure of the bubble; ρ represents the density of the fluid; v represents kinematic viscosity of the fluid;
step S34: the material attribute, setting the two fluid phases as a first phase which is a plasma bubble phase water-vapor, and setting the material attribute as ideal gas; the second phase is water-liquid or air, and the material properties are water and air respectively;
step S35: setting the boundary condition of a flow field according to the actual situation of the laser shock peening process, and setting the pressure outlet as the pressure outlet boundary condition by combining the named boundary name in geometric modeling;
step S36: solving method and control, calculating and simulating the compressible multiphase flow laser-induced plasma bubble transient flow field;
step S37: initializing, wherein global initialization of a flow field is adopted during initialization, and a Region where a plasma bubble is located is marked by Region addition of FLUENT;
step S38: calculating settings: setting the calculation time step to 1 × 10 -8 s, the number of calculation steps is 100, data is stored once in each step, and the calculation time step length of each step can be increased and the number of calculation steps can be reduced appropriately when water constraint exists so as to reduce the calculation amount.
5. The gas-liquid-solid coupling calculation method for laser shock peening according to claim 1, wherein the concrete steps of the transient structure solution in step S4 are as follows:
step S41: defining materials, selecting steel in a material library as the material of a workpiece in the model, and keeping the material attribute at a default value;
step S42: solving time, namely matching the solving step length and the solving time of the fluid-solid module according to the time step length and the calculation time set defined in the fluid model;
step S43: applying a load, wherein a surface named wall-solid in geometric modeling is defined as an application surface of the load, and an application mode is defined as a fluid-solid interface so that the load of a flow field can be applied to the surface;
step S44: applying fixation, and defining a surface named fixed-support in geometric modeling as a surface of a fixed support to constrain the geometric body of the workpiece;
step S45: adding solving quantity, inserting response of deformation, strain and stress needing to be solved.
6. The gas-liquid-solid coupling calculation method for laser shock peening according to claim 1, wherein the gas-liquid-solid system coupling setting in step S5 comprises the following specific steps:
step S51: solving time, matching the time step length and the calculation time defined in the fluid model and the structural model, and enabling the solving step length and the solving time of the system coupling module to be consistent with those of the other two modules;
step S52: data transmission, namely adding two surfaces named wall-fluid and wall-solid in geometric modeling as data transmission surfaces, wherein the transmitted physical attribute is force;
step S53: the data transmission sequence is that the simulation result data of the transient fluid is used as the boundary condition of the transient structure calculation simulation module to be loaded, and then the calculation result of the transient structure is used as the boundary condition of the next calculation time step fluid simulation module to be loaded;
step S54: and (4) calculating and solving, namely setting data output as all time steps, and then saving the model to carry out calculation and solving.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810672217.5A CN109063248B (en) | 2018-06-26 | 2018-06-26 | Gas-liquid-solid coupling calculation method for laser shock peening |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810672217.5A CN109063248B (en) | 2018-06-26 | 2018-06-26 | Gas-liquid-solid coupling calculation method for laser shock peening |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109063248A CN109063248A (en) | 2018-12-21 |
| CN109063248B true CN109063248B (en) | 2022-11-22 |
Family
ID=64821099
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810672217.5A Active CN109063248B (en) | 2018-06-26 | 2018-06-26 | Gas-liquid-solid coupling calculation method for laser shock peening |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109063248B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109885942B (en) * | 2019-02-25 | 2023-09-01 | 北京易盛泰和科技有限公司 | Open-top type climate chamber and gas uniformity simulation verification method thereof |
| CN111914383A (en) * | 2019-05-10 | 2020-11-10 | 南京航空航天大学 | T-shaped structure laser welding simulation method based on shell-entity coupling unit |
| CN112364362B (en) * | 2020-11-16 | 2023-12-29 | 宁波九寰适创科技有限公司 | Parallel multi-layer self-adaptive local encryption method oriented to fluid simulation direction |
| CN113935257B (en) * | 2021-09-30 | 2023-05-05 | 燕山大学 | Gas-liquid-solid three-phase flow simulation method |
| CN113935169B (en) * | 2021-10-14 | 2022-09-23 | 深圳泽森软件技术有限责任公司 | Physical simulation method, physical simulation device, computer equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104951629A (en) * | 2015-07-14 | 2015-09-30 | 无锡沃尔福汽车技术有限公司 | Direct fluid-solid coupled heat transfer analysis method for whole engine water jacket |
| CN105512433A (en) * | 2016-01-12 | 2016-04-20 | 北京航空航天大学 | Fluid-solid node two-phase flow modeling method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007122269A (en) * | 2005-10-26 | 2007-05-17 | Sony Corp | Method of coupled numerical simulation of fluid-structure, and storage device program for coupled numerical simulation of fluid-structure |
-
2018
- 2018-06-26 CN CN201810672217.5A patent/CN109063248B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104951629A (en) * | 2015-07-14 | 2015-09-30 | 无锡沃尔福汽车技术有限公司 | Direct fluid-solid coupled heat transfer analysis method for whole engine water jacket |
| CN105512433A (en) * | 2016-01-12 | 2016-04-20 | 北京航空航天大学 | Fluid-solid node two-phase flow modeling method |
Non-Patent Citations (2)
| Title |
|---|
| 圆柱壳容器内部可燃气体爆炸冲击作用的双向流固耦合模拟;胡可等;《建筑结构学报》;20160525;全文 * |
| 激光支持爆轰波的仿真研究;胡雅骥等;《宇航材料工艺》;20090415(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109063248A (en) | 2018-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109063248B (en) | Gas-liquid-solid coupling calculation method for laser shock peening | |
| Li et al. | 3D full coupling model for strong interaction between a pulsating bubble and a movable sphere | |
| Ingram et al. | Developments in Cartesian cut cell methods | |
| CN110222467B (en) | Analysis method for swing high-speed jet characteristic of underwater engine spray pipe | |
| Chan-Braun et al. | Force and torque acting on particles in a transitionally rough open-channel flow | |
| CN107217133B (en) | Laser impact intensified Finite Element Method | |
| Fourtakas et al. | On the coupling of incompressible SPH with a finite element potential flow solver for nonlinear free-surface flows | |
| Banari et al. | Efficient GPGPU implementation of a lattice Boltzmann model for multiphase flows with high density ratios | |
| Kawamura et al. | Large eddy simulation of a flow past a free surface piercing circular cylinder | |
| Shopov et al. | Interaction of a deformable bubble with a rigid wall at moderate Reynolds numbers | |
| Spandan et al. | Drag reduction in numerical two-phase Taylor–Couette turbulence using an Euler–Lagrange approach | |
| Zabala et al. | Wave-induced real-fluid effects in marine energy converters: Review and application to OWC devices | |
| Zhu et al. | Numerical study on vortex-induced vibration responses of a circular cylinder attached by a free-to-rotate dartlike overlay | |
| Schwartzentruber et al. | Multiscale particle-continuum simulations of hypersonic flow over a planetary probe | |
| Kim et al. | Shape optimization of a hyperloop pod’s head and tail using a multi-resolution morphing method | |
| Zhang et al. | Prediction of the hydrodynamic disturbance characteristics for two projectiles launched successively underwater based on the radial basis function neural network | |
| CN116499698A (en) | Pneumatic and kinematic mechanical coupling analysis method for magnetic levitation flight wind tunnel magnetic levitation platform | |
| CN108763692B (en) | Efficient wave making method for ship numerical pool | |
| Markus et al. | Complementary numerical–experimental benchmarking for shape optimization and validation of structures subjected to wave and current forces | |
| Hunt et al. | Viscous coupling of shear-free turbulence across nearly flat fluid interfaces | |
| Hojeij et al. | Experimental study and numerical modeling of mixing by air injection in yield stress fluids using the OpenFOAM software | |
| Chen et al. | Internal structure of a jet nozzle for coalbed methane mining based on airfoil curves | |
| Miao et al. | Computation of three-dimensional multiphase flow dynamics by fully-coupled immersed flow (fcif) solver | |
| CN110414141B (en) | Three-dimensional numerical simulation method for liquid drop atomization in process of transonic flow of compressible fluid | |
| Liu et al. | A volume of fluid based method for consistent flux computation in large-density ratio two-phase flows and its application in investigating droplet bag breakup behavior |
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 |