Flexible fiber fluid-solid coupling method

Abstract

The invention belongs to the technical field of flexible fiber washing and conveying analog simulation, and relates to a flexible fiber fluid-solid coupling method, which comprises the following steps of: firstly, respectively constructing a fluid model containing flexible fibers and a flexible fiber model; constructing a bidirectional fluid-solid coupling system based on the flexible fiber model and the fluid model; then carrying out analogue simulation and physical experiment on the motion evolution form of the single flexible fiber in the flow field by utilizing a bidirectional fluid-solid coupling system, and proposing evaluation indexes of the flexibility of the flow field of the flexible fiber and the vibration amplitude of the flow field of the flexible fiber; and finally, correcting the flexible fiber fluid-solid coupling model through evaluation indexes, and further exploring bending, winding, entanglement clustering and opening conditions of the flexible fibers. According to the flexible fiber fluid-solid coupling method, simulation of motion evolution behaviors of the flexible fibers in a flow field can be effectively achieved, and good expansibility and adaptability are achieved.

Description

Technical Field

[0001] The present invention belongs to the technical field of flexible fiber water washing transportation simulation, and relates to a flexible fiber fluid-solid coupling method. Background Art

[0002] The large-capacity flexible fiber water transportation technology is one of the key technologies for producing green fibers. It has the advantages of high efficiency, large capacity, strong adaptability, and small environmental impact, and is of great significance to the new energy strategy. However, the forces acting on flexible fibers in the flow field are complex, and their motion postures are unpredictable. During the transportation process, due to the interaction between fibers and fluids and the mutual contact and collision between fibers, the fiber morphology changes, such as bending, winding, clustering, etc. This directly affects the quality of fiber products and poses a risk of production line stoppage, seriously affecting production efficiency and cost.

[0003] Flexible fibers are prone to bending, winding, entanglement, and loosening during the water washing transportation process, which cannot be observed during the actual production preparation process. It is necessary to show through simulation results the tendency of flexible fibers to bend, wind, entangle, and loosen during the water washing transportation process. However, at present, there is a lack of theoretical guidance in actual production, and it is difficult to carry out research on the flow evolution behavior.

[0004] To solve the above problems, the literature (Influence of the flow field state in the wool conveying pipe on the opening treatment effect of fiber bundles [J]. Advanced Textile Technology, 2024, 32(3): 29-37. DOI: 10.19398 / j.att.202308011) used the fluid simulation software Cradle CFD to perform coupled calculations on the motion process of fiber bundles in the flow field of the wool conveying pipe, compared the distribution characteristics of the flow field at different positions in the wool conveying pipe and the fiber motion trajectories at different times, and explored the opening mechanism of fiber bundles during transportation in the wool conveying pipe and the influence of the flow field state in the wool conveying pipe on fiber motion. However, this simulation method does not consider the mechanical properties of flexible fibers and cannot accurately correct the coupled model, and the simulation results are not persuasive.

[0005] Therefore, it is of great significance to study a flexible fiber fluid-solid coupling method to solve the problems existing in the prior art. Summary of the Invention

[0006] The purpose of the present invention is to solve the problems existing in the prior art and provide a flexible fiber fluid-solid coupling method.

[0007] To achieve the above purpose, the technical solution adopted by the present invention is as follows:

[0008] A flexible fiber fluid-structure interaction method is applied to the dynamic simulation of the movement evolution characteristics of large aspect ratio flexible fiber water washing and conveying. The large aspect ratio means that the aspect ratio is 100 - 500. For example, for Lyocell short cut fibers, the length is commonly 6 - 50 mm (usually 38 mm in enterprise production and preparation), and the fiber diameter is commonly 60 - 100 μm. The flexible fiber fluid-structure interaction method includes the following steps:

[0009] Step 1: Construct a fluid model containing flexible fibers based on the Computational Fluid Dynamics (CFD) method;

[0010] Step 2: Describe the flexible fibers using the Discrete Element Method (DEM) combined with the sphere-column unit model to obtain multiple fiber sphere-column units of the same length. The multiple fiber sphere-column units of the same length are connected and combined to form a flexible fiber model;

[0011] Step 3: Introduce the porosity of the fluid phase and establish a mass control volume equation containing the fluid and flexible fibers; introduce the interaction force between the fluid and flexible fibers into the source term of the force in the momentum control volume equation to obtain a momentum control volume equation containing the source term of the flexible fiber acting force;

[0012] Step 4: The fluid model containing flexible fibers and the flexible fiber model are interpolated and coupled with each other through the trilinear method in the finite element volume method. The fluid model containing flexible fibers and the flexible fiber model together form a two-way fluid-structure interaction system; there is an interaction force between the flexible fiber model and the fluid model in the two-way fluid-structure interaction system; the coupling simulation process is as Figure 1 shown; through the divided finite element mesh, the flow field is spatially discretized, and trilinear interpolation is used to transfer data between the fluid grid nodes and the discrete particle positions. The steps are as follows:

[0013] Through the divided finite element mesh, the flow field is spatially discretized, and trilinear interpolation is used to transfer data between the fluid grid nodes and the discrete particle positions. The steps are as follows:

[0014] ① Locate the fluid grid unit where the fiber is located, and locate its grid unit according to the fiber coordinates The vertex coordinates of the unit are to , and find the eight vertex coordinates of the unit and their corresponding fluid variable velocities;

[0015] ② Calculate the distance between each finite element grid node and the particle center position. The relative distances (weight coefficients) in each direction are:

[0016] ;

[0017] Among them, , and calculate in the same way.

[0018] ③ Calculate the vertex weights. The weight of each vertex is the product of the weight coefficients in three directions. For vertex , its weight calculation formula is:

[0019] ;

[0020] ④ According to the trilinear interpolation formula, weighted sum the values of the eight vertices of the grid to obtain the interpolation result at the fiber position. Let the values of the eight vertices be , and the interpolation formula is:

[0021] ;

[0022] Expand it to:

[0023] ;

[0024] Among them, is the coordinate is the weight value of the vertex;

[0025] is the coordinate of the vertex.

[0026] In the two-way coupling, the fluid variables at time t are interpolated to the fiber position, and the force on the fiber also needs to be fed back to the fluid grid. At this time, the interpolation weight is also determined by the weight of the trilinear interpolation. The force of each particle will be distributed to the surrounding eight nodes, and the force received by each node is the fiber force multiplied by the corresponding weight to obtain the interpolation result at the fiber position. Subsequently, at time t+1, time t+2, time t+3... until the calculation time is reached.

[0027] Step 5: Use the two-way fluid-structure coupling system to simulate and physically experiment on the motion evolution form of a single flexible fiber in the flow field to obtain the minimum bending amount of the flexible fiber at any time (i.e., the bending distance between the two endpoints when the amplitude displacement of the flexible fiber is the largest) and the maximum displacement amount and the minimum displacement amount of the end point of the flexible fiber, and then calculate the indexes for evaluating the interaction force between the fluid and the flexible fiber, namely the flow field curvature of the flexible fiber and the flow field amplitude of the flexible fiber;

[0028] In the simulation process, the motion posture animation of the flexible fiber in the flow field is decomposed into frames at 0.001 s and saved. Through image processing technology, the value of the minimum bending amount is measured. Substitute it into the following formula (1) to obtain the value of the flow field curvature of the flexible fiber; measure each frame result through image processing technology to measure the maximum value and the minimum Value, substitute its two values into the following formula (2) to obtain the central axis of the amplitude of the flexible fiber flow field , on the central axis of the amplitude of the flexible fiber flow field , based on the central axis of the amplitude of the flexible fiber flow field, measure the displacement of the flexible fiber in the time series (at each continuous time interval) , substitute it into the following formula (3) to obtain the amplitude of the flexible fiber flow field.

[0029] Step 6: Introduce a drag force model correction factor into the source term of the force in the momentum control volume equation , correct the strength of the force between the fluid and the flexible fiber through the evaluation index in Step 5, and perform simulations according to the number of fibers in the actual working condition, so as to judge the bending and winding, entanglement and loosening of the flexible fibers

[0030] The calculation formula for the curvature of the flexible fiber flow field is as follows

[0031] (1);

[0032] Among them, S is the curvature of the flexible fiber flow field is the initial length of the flexible fiber, with the unit of mm is the minimum bending amount of the flexible fiber at any moment, with the unit of mm

[0033] The calculation formula for the amplitude of the flexible fiber flow field is as follows

[0034] (2);

[0035] (3);

[0036] Among them is the central axis of the amplitude of the flexible fiber flow field, with the unit of mm is the maximum displacement of the end point of the flexible fiber, with the unit of mm is the minimum displacement of the end point of the flexible fiber, with the unit of mm is the amplitude of the flexible fiber flow field in the time series (the time series is the simulation time interval, and the simulation time interval in the present invention is set to 0.001 s), with the unit of mm is the displacement of the flexible fiber in the time series, with the unit of mm

[0037] As a preferred technical solution

[0038] For a flexible fiber fluid-solid coupling method as described above, the specific steps for constructing the fluid model in Step 1 are as follows

[0039] Step 1.1: Based on the physical geometric model of the conveying pipeline, extract the fluid domain, divide it using the finite element method to obtain fluid meshes, and introduce the volume fraction occupied by the flexible fibers in the fluid meshes to obtain the fluid model;

[0040] Step 1.2: Discretize the fluid dynamics equations based on the finite element method to obtain a non-linear equation system with the grid node parameter values as the solution variables; solve the non-linear equation system to obtain the velocity and pressure changes of the flow field at the grid nodes;

[0041] Step 1.3: Based on the fluid meshes, establish a fluid model containing fibers.

[0042] In a flexible fiber fluid-structure interaction method as described above, in the flexible fiber model of Step 2, each fiber sphere-column unit is composed of a cylinder and two end spheres. As Figure 3 shown, adjacent fiber sphere-column units are connected by sphere elastic members. Forces and torques are transmitted through the spherical nodes in the fiber sphere-column units. The fiber sphere-column units can undergo morphological evolutions such as bending, torsion, and stretching under the external forces received, thereby changing the flexibility of the fibers.

[0043] In a flexible fiber fluid-structure interaction method as described above, the flexible fiber model of Step 2 includes the dynamic equations of the forces exerted by the fluid on the flexible fibers. The dynamic equations of the forces exerted by the fluid on the flexible fibers include the force translation equation and the torque rotation equation;

[0044] The force translation equation is described as:

[0045] (4);

[0046] The torque rotation equation is described as:

[0047] (5);

[0048] Where, is the mass of a single flexible fiber;

[0049] is the translational velocity of the th fiber node sphere;

[0050] is the rotational angular velocity vector of the node sphere;

[0051] is the gravitational acceleration;

[0052] is the time;

[0053] is the moment of inertia;

[0054] is the contact force of the flexible fiber and ;

[0055] is the total drag force of the fluid on the flexible fiber;

[0056] is the flexible fiber and tangential moment;

[0057] is the flexible fiber and normal moment.

[0058] A flexible fiber fluid-structure interaction method as described above, the calculation formula is:

[0059] (6);

[0060] wherein, is the drag force of the th flexible fiber;

[0061] is the fluid unit volume;

[0062] is the number of flexible fibers in the fluid unit;

[0063] is the interpolation weight function (distributing the contribution of fibers to surrounding grid nodes according to trilinear interpolation method);

[0064] is the position of the flexible fiber;

[0065] is the fluid grid center coordinate.

[0066] A flexible fiber fluid-structure interaction method as described above, based on the finite element volume cell method, solves equations (4) and (5), and the steps are as follows:

[0067] First, equations (4) and (5) are respectively simplified to:

[0068] (7);

[0069] (8);

[0070] In the formula, is the acceleration of fiber i; is the moment of inertia of fiber i, is the angular acceleration of fiber i; is the force on fiber i, is the resultant external moment;

[0071] Then, by integrating, the velocity of the fiber at time t is obtained:

[0072] (9);

[0073] In the formula, is the acceleration of fiber i at moment; is the angular acceleration of fiber i at moment;

[0074] Integrate again to update the motion attitude of the fiber:

[0075] (10);

[0076] In the formula, is the displacement of fiber i obtained by solving, is the displacement of fiber i at moment; By integrating with respect to time, the displacement of the fiber at the next moment is obtained, and so on, to obtain the change in fiber displacement within the simulation time series;

[0077] is the rotation angle of fiber i obtained, is the rotation angle of fiber i at moment; By integrating with respect to time, the rotation angle of the fiber at the next moment is obtained, and so on, to obtain the change in fiber rotation angle within the simulation time series.

[0078] As described above, in a flexible fiber fluid-structure interaction method, the mass control volume equation of the fluid and the flexible fiber and the momentum control volume equation including the source term of the flexible fiber acting force in step 3 are shown in equations (11) and (12) respectively;

[0079] (11);

[0080] The momentum conservation equation including the source term of the flexible fiber acting force is described as:

[0081] (12);

[0082] Among them, is the fluid porosity (the fluid porosity represents the volume fraction of the fluid in the control volume);

[0083] is the fluid density;

[0084] is the fluid velocity;

[0085] is the fluid pressure;

[0086] is the Hamiltonian operator;

[0087] is the fluid viscosity;

[0088] is the time;

[0089] is the acceleration due to gravity;

[0090] is the stress tensor;

[0091] is the reaction force of the flexible fiber on the fluid, where and satisfy Newton's third law, and there is a relationship between the action force and the reaction force, that is ;

[0092] The fluid porosity The calculation formula is:

[0093] (13);

[0094] The sum of the fiber volume fraction and the fluid volume fraction in the entire control volume is 1, that is, the fiber volume is distributed to the adjacent fluid grid nodes or the centers of the grid cells according to its position in the fluid grid;

[0095] Among them, is the volume of the k-th ball-column unit fiber in the control volume;

[0096] V is the volume of the control volume;

[0097] is the interpolation weight function (distributing the contribution of the fiber to the surrounding grid nodes according to the trilinear interpolation method);

[0098] is the position of the fiber;

[0099] is the coordinate of the fluid grid center.

[0100] A flexible fiber fluid-structure interaction method as described above. In step 2, the discrete element method is used in combination with the ball-column unit model to describe the flexible fiber, which means using physical property parameters and mechanical property parameters for description. The physical property parameters are fiber diameter, fiber length, fiber density, fiber number, fiber elastic modulus, and fiber Poisson ratio. The mechanical property parameters are: the normal force of the ball-column unit and the tangential force , the normal force at the node ball , the tangential force , the torsional moment , the bending moment , the resultant torsional moment and the resultant bending moment ;

[0101] The parameters of the flexible fiber are measured according to the physical fiber of the research object. The fiber diameter is 60 - 120 μm, the fiber length is 5 - 100 mm, the fiber density is 750 - 2000 kg / m 3 , the number of fibers is 1 - 9999, the elastic modulus of the fiber is 100 - 260 MPa, and the Poisson's ratio of the fiber is 0.2 - 0.35;

[0102] The calculation formulas for each mechanical property parameter are as follows:

[0103] (14);

[0104] (15);

[0105] (16);

[0106] (17);

[0107] (18);

[0108] (19);

[0109] (20);

[0110] (21);

[0111] Among them, is the elastic modulus of the flexible fiber;

[0112] is the node ball damping coefficient, and the value is 0.035;

[0113] is the shear modulus of the flexible fiber, , v is the Poisson's ratio;

[0114] is the cross-sectional area of the spherical column unit;

[0115] is the polar moment of inertia of the spherical column unit cross-section, is the radius of the spherical column unit;

[0116] and are the normal displacement and tangential displacement of the spherical-cylindrical element, respectively;

[0117] and are the normal stiffness and shear stiffness at the spherical node, respectively;

[0118] is the distance between the spherical nodes;

[0119] and are the torsional stiffness and bending stiffness at the spherical node, respectively;

[0120] is the mass and the moment of inertia of the relative normal velocity between two spherical nodes;

[0121] is the mass and the moment of inertia of the relative tangential velocity between two spherical nodes;

[0122] is the mass and the moment of inertia of the relative torsional angular velocity between two spherical nodes;

[0123] is the mass and the moment of inertia of the relative bending angular velocity between two spherical nodes.

[0124] For a flexible fiber fluid-structure interaction method as described above, the calculation formula for the drag force of the mutual force between the flexible fiber model and the fluid model in the two-way fluid-structure interaction system in step 4 is:

[0125] (22);

[0126] where the momentum exchange coefficient between the fluid and the fiber phase is expressed as:

[0127] (23);

[0128] The drag coefficient C D is expressed as:

[0129] (24);

[0130] The calculation formula for the Reynolds number Re is:

[0131] (25);

[0132] wherein, is the fluid velocity;

[0133] is the flexible fiber velocity;

[0134] is the fluid porosity;

[0135] is the fluid density;

[0136] is the flexible fiber diameter;

[0137] is the fluid viscosity;

[0138] is the acceleration of gravity.

[0139] A flexible fiber fluid-structure interaction method as described above introduces a drag force model correction factor in the source term of the force in the momentum control volume equation. The corrected drag force model expression is:

[0140] (26).

[0141] A flexible fiber fluid-structure interaction method as described above performs mutual interpolation coupling calculation by the trilinear method in the finite element volume method in step 4. Specifically, equations (7) and (8) are rewritten as a general formula:

[0142] (27);

[0143] wherein, the fluid density vector with porosity is expressed as:

[0144] (28);

[0145] The convective terms of the three components in the spatial coordinate system in the fluid model , and are respectively expressed as:

[0146] (29);

[0147] (30);

[0148] (31);

[0149] The viscous fluxes of the three components in the spatial coordinate system in the fluid model , and is expressed as:

[0150] (32);

[0151] (33);

[0152] (34);

[0153] is the force source term, including the interaction force between the fiber and the fluid:

[0154] (35);

[0155] wherein, is the component of the fluid velocity on the axis in the space coordinate system;

[0156] is the component of the fluid velocity on the axis in the space coordinate system;

[0157] is the component of the fluid velocity on the axis in the space coordinate system;

[0158] is the convective term on the axis in the fluid dynamics equations;

[0159] is the convective term on the axis in the fluid dynamics equations;

[0160] is the convective term on the axis in the fluid dynamics equations;

[0161] is the viscous flux on the axis in the fluid dynamics equations;

[0162] is the viscous flux on the axis in the fluid dynamics equations;

[0163] is the viscous flux on the axis in the fluid dynamics equations;

[0164] is the fluid porosity, is the fluid pressure, is the fluid viscosity tensor;

[0165] , , are respectively the force source terms in three spatial directions;

[0166] is the viscous force perpendicular to the axis and along the axis direction;

[0167] is the viscous force perpendicular to the axis and along the axis direction;

[0168] is the viscous force perpendicular to the axis and along the axis direction;

[0169] is the viscous force perpendicular to the axis and along the axis direction;

[0170] is the viscous force perpendicular to the axis and along the axis direction;

[0171] is the viscous force perpendicular to the axis and along the axis direction;

[0172] is the viscous force perpendicular to the axis and along the axis direction;

[0173] is the viscous force perpendicular to the axis and along the axis direction;

[0174] is the viscous force perpendicular to the axis and along the axis direction.

[0175] Advantageous effects:

[0176] (1) A flexible fiber fluid-structure interaction method of the present invention, based on the spherical column unit model modeling method, establishes a flexible fiber fluid-structure interaction system, which can effectively simulate the motion evolution behavior of flexible fibers in a fluid field, and has good scalability and adaptability, and is applicable to the simulation of the motion evolution characteristics of large aspect ratio flexible bodies in a fluid field.

[0177] (2) A flexible fiber fluid-structure interaction method of the present invention. The discrete element method (DEM) is used to study particle dynamics, which can effectively simulate the dynamic characteristics of particles with different shapes. The computational fluid dynamics (CFD) method can accurately capture fluid field information. By self-defining a coupling program, the CFD-DEM model can simulate the solid-liquid two-phase flow system, capture particle collision and position attitude information, and achieve numerical simulation of fluid-solid interaction. This method provides a new theoretical basis for the entanglement behavior and agglomeration characteristics of flexible fibers with a large aspect ratio in the flow field, and provides guidance for future fiber applications in the new energy industry. BRIEF DESCRIPTION OF THE DRAWINGS

[0178] Figure 1 It is the coupling calculation flow chart of the flexible fiber fluid-structure interaction method of the present invention;

[0179] Figure 2 It is the schematic diagram of the bending degree and amplitude of a single fiber in the flow field of the present invention; wherein, Figure 2 (a) in it is the schematic diagram for solving the bending degree, Figure 2 (b) in it is the maximum displacement in the y-axis direction of the amplitude, Figure 2 (c) in it is the minimum displacement in the y-axis direction of the amplitude; the dotted line represents the amplitude central axis T, and the red line represents the fiber;

[0180] Figure 3 It is the schematic diagram of the flexible fiber model based on the sphere-column unit model of the present invention;

[0181] Figure 4 It is the verification experimental scheme diagram of the flexible fiber fluid-structure interaction method of the present invention;

[0182] Figure 5 It is the comparison diagram of the simulation and experiment of the motion form evolution of a single fiber in the flow field of the present invention;

[0183] Figure 6 It is the change diagram of the amplitude index of the fiber flow field under different correction factors of the present invention;

[0184] Figure 7 It is the change diagram of the bending degree index of the fiber flow field under different correction factors of the present invention;

[0185] Figure 8 It is the evolution behavior and velocity distribution diagram of the fiber in the flow field at different times. Among them, Figure 8 (a) in it is the evolution behavior and velocity distribution of the fiber at 0.025 s, Figure 8 (c) in it is the evolution behavior and velocity distribution of the fiber at 0.038 s, Figure 8 (e) in it is the evolution behavior and velocity distribution of the fiber at 0.089 s, Figure 8 (g) in it is the evolution behavior and velocity distribution of the fiber at 0.165 s,Figure 8 in (b), Figure 8 in (d), Figure 8 in (f) and Figure 8 in (h) are respectively Figure 8 in (a), Figure 8 in (c), Figure 8 in (e) and Figure 8 the partial enlarged views at corresponding moments of (g);

[0186] Wherein, 1 - node ball, 2 - ball - column unit, 3 - water source, 4 - aluminum profile frame, 5 - CCD high - definition camera, 6 - background board, 7 - lighting lamp, 8 - water flow velocity meter, 9 - universal rotary pipeline fixing clamp, 10 - hair - conveying pipe, 11 - recovery tank. Specific embodiments

[0187] The present invention will be further described below in conjunction with specific embodiments. It should be understood that these embodiments are only used to illustrate the present invention and not to limit the scope of the present invention. In addition, it should be understood that after reading the content taught by the present invention, those skilled in the art can make various changes or modifications to the present invention, and these equivalent forms also fall within the scope defined by the appended claims of this application.

[0188] A flexible fiber fluid - solid coupling method is applied to the dynamic simulation of the movement evolution characteristics of large - aspect - ratio flexible fiber water - washing transportation. The large - aspect - ratio means that the aspect - ratio is 100 - 500. The flexible fiber fluid - solid coupling method includes the following steps:

[0189] Step 1: Construct a fluid model containing flexible fibers based on the computational fluid dynamics method (CFD), specifically as follows:

[0190] Step 1.1: Based on the physical geometric model of the conveying pipeline, extract the fluid domain, divide it by the finite - element method to obtain fluid meshes, and introduce the volume fraction occupied by the flexible fibers in the fluid meshes to obtain the fluid model;

[0191] Step 1.2: Discretize the computational fluid dynamics equations based on the finite - element method to obtain a non - linear equation system with the grid - node parameter values as the solution variables; solve the non - linear equation system to obtain the velocity and pressure changes of the flow field at the grid nodes;

[0192] Step 1.3: Based on the fluid meshes, establish a fluid model containing fibers.

[0193] Step 2: Use the discrete element method (DEM) combined with the ball - column unit model to describe the flexible fibers, and obtain a plurality of fiber ball - column units with the same length. The plurality of fiber ball - column units with the same length are connected and combined to form a flexible fiber model;

[0194] The discrete element method is combined with the spherical column unit model to describe flexible fibers, which means using physical property parameters and mechanical property parameters for description. The physical property parameters are fiber diameter, fiber length, fiber density, fiber quantity, fiber elastic modulus, and fiber Poisson's ratio. The mechanical property parameters are: the normal force of the spherical column unit and the tangential force , the normal force at the node ball , the tangential force , the torsional moment , the bending moment , the combined torsional moment and the combined bending moment ;

[0195] The parameters of the flexible fiber are measured according to the physical fiber of the research object. The fiber diameter is 60 - 120 μm, the fiber length is 5 - 100 mm, the fiber density is 750 - 2000 kg / m 3 , the fiber quantity is 1 - 9999, the fiber elastic modulus is 100 - 260 MPa, and the fiber Poisson's ratio is 0.2 - 0.35;

[0196] The calculation formulas for each mechanical property parameter are as follows:

[0197] (14);

[0198] (15);

[0199] (16);

[0200] (17);

[0201] (18);

[0202] (19);

[0203] (20);

[0204] (21);

[0205] Among them, is the elastic modulus of the flexible fiber;

[0206] is the damping coefficient of the node ball, and the value is 0.035;

[0207] is the shear modulus of the flexible fiber, , v is Poisson's ratio;

[0208] is the cross-sectional area of the spherical column unit;

[0209] is the polar moment of inertia of the spherical column unit cross-section, is the radius of the spherical column unit;

[0210] and are the normal displacement of the spherical column unit and the tangential displacement of the spherical column unit respectively;

[0211] and are the normal stiffness and shear stiffness at the node ball respectively;

[0212] is the distance between the node balls;

[0213] and are the torsional stiffness and bending stiffness at the node ball respectively;

[0214] is the mass and the moment of inertia of the relative normal velocity between two node balls;

[0215] is the mass and the moment of inertia of the relative tangential velocity between two node balls;

[0216] is the mass and the moment of inertia of the relative torsional angular velocity between two node balls;

[0217] is the mass and the moment of inertia of the relative bending angular velocity between two node balls;

[0218] Each fiber spherical column unit is composed of a cylinder and spheres at both ends, as Figure 3 shown, adjacent fiber spherical column units 2 are connected by node balls 1;

[0219] The flexible fiber model includes the dynamic equations of the fluid acting on the flexible fiber. The dynamic equations of the fluid acting on the flexible fiber include the force translation equation and the moment rotation equation;

[0220] The force translation equation is described as:

[0221] (4);

[0222] The torque rotation equation is described as:

[0223] (5);

[0224] Wherein, is the mass of a single flexible fiber;

[0225] is the translational velocity of the th fiber node sphere;

[0226] is the rotational angular velocity vector of the node sphere;

[0227] is the gravitational acceleration;

[0228] is the time;

[0229] is the moment of inertia;

[0230] is the contact force between the flexible fiber and ;

[0231] is the total drag force of the fluid on the flexible fiber;

[0232] is the tangential torque between the flexible fiber and ;

[0233] is the normal torque between the flexible fiber and ;

[0234] The calculation formula of

[0235] is: (6);

[0236] Wherein, is the drag force of the th flexible fiber;

[0237] is the fluid unit volume;

[0238] is the number of flexible fibers in the fluid unit;

[0239] is the interpolation weight function;

[0240] is the position of the flexible fiber;

[0241] are the central coordinates of the fluid grid;

[0242] Based on the finite element volume element method, equations (4) and (5) are solved as follows:

[0243] First, equations (4) and (5) are respectively simplified to:

[0244] (7);

[0245] (8);

[0246] wherein is the acceleration of fiber i; is the moment of inertia of fiber i, is the angular acceleration of fiber i; is the force acting on fiber i, is the resultant external torque;

[0247] Then, integration is performed to obtain the motion velocity of the fiber at time t:

[0248] (9);

[0249] wherein, is the acceleration of fiber i at time; is the angular acceleration of fiber i at time;

[0250] Integrate again to update the motion attitude of the fiber:

[0251] (10);

[0252] wherein, is the displacement of fiber i obtained by solution, is the displacement of fiber i at time; is the rotation angle of fiber i obtained, is the rotation angle of fiber i at time.

[0253] Step 3: Introduce the porosity of the fluid phase and establish the mass control volume equation including the fluid and the flexible fiber; introduce the interaction force between the fluid and the flexible fiber into the source term of the force in the momentum control volume equation to obtain the momentum control volume equation including the source term of the flexible fiber acting force;

[0254] The mass control volume equation containing the fluid and flexible fibers and the momentum control volume equation containing the source term of the acting force of the flexible fibers are shown in equations (11) and (12) respectively;

[0255] (11);

[0256] The momentum conservation equation containing the source term of the acting force of the flexible fibers is described as:

[0257] (12);

[0258] Among them, is the fluid porosity (the fluid porosity represents the volume fraction of the fluid in the control volume);

[0259] is the fluid density;

[0260] is the fluid velocity;

[0261] is the fluid pressure;

[0262] is the Hamiltonian operator;

[0263] is the fluid viscosity;

[0264] is the time;

[0265] is the gravitational acceleration;

[0266] is the stress tensor;

[0267] is the reaction force of the flexible fibers on the fluid, among which, and satisfy Newton's third law, and there is a relationship of action and reaction, that is ;

[0268] The fluid porosity The calculation formula of is:

[0269] (13);

[0270] The sum of the fiber volume fraction and the fluid volume fraction in the entire control volume is 1, that is, the fiber volume is distributed to the adjacent fluid grid nodes or the centers of the grid elements according to its position in the fluid grid;

[0271] Among them, is the volume of the fiber in the k-th sphere-column unit in the control volume;

[0272] V is the volume of the control volume;

[0273] is the interpolation weight function (allocating the contribution of fibers to the surrounding grid nodes according to the trilinear interpolation method);

[0274] is the position of the fiber;

[0275] is the central coordinate of the fluid grid;

[0276] Step 4: The fluid model containing flexible fibers and the flexible fiber model are coupled by mutual interpolation through the trilinear method in the finite element volume method. The fluid model containing flexible fibers and the flexible fiber model together constitute a two-way fluid-structure interaction system;

[0277] The mutual interpolation coupling calculation is carried out through the trilinear method in the finite element volume method. Specifically, equations (7) and (8) are rewritten as a general formula:

[0278] (27);

[0279] Among them, the fluid density vector with porosity is expressed as:

[0280] (28);

[0281] The convective terms of the three components in the spatial coordinate system of the fluid model , and are respectively expressed as:

[0282] (29);

[0283] (30);

[0284] (31);

[0285] The viscous fluxes of the three components in the spatial coordinate system of the fluid model , and are expressed as:

[0286] (32);

[0287] (33);

[0288] (34);

[0289] is the force source term, including the interaction force between the fiber and the fluid:

[0290] (35);

[0291] Among them, is the component of the fluid velocity on the axis in the space coordinate system;

[0292] is the component of the fluid velocity on the axis in the space coordinate system;

[0293] is the component of the fluid velocity on the axis in the space coordinate system;

[0294] is the convective term on the axis in the fluid dynamics equations;

[0295] is the convective term on the axis in the fluid dynamics equations;

[0296] is the convective term on the axis in the fluid dynamics equations;

[0297] is the viscous flux on the axis in the fluid dynamics equations;

[0298] is the viscous flux on the axis in the fluid dynamics equations;

[0299] is the viscous flux on the axis in the fluid dynamics equations;

[0300] is the fluid porosity, is the fluid pressure, is the fluid viscosity tensor;

[0301] 、 、 are respectively the force source terms in the three directions of space;

[0302] is that the viscous force is perpendicular to the axis and along the axis direction;

[0303] The viscous force is perpendicular to the axis and along the axis direction;

[0304] The viscous force is perpendicular to the axis and along the axis direction;

[0305] The viscous force is perpendicular to the axis and along the axis direction;

[0306] The viscous force is perpendicular to the axis and along the axis direction;

[0307] The viscous force is perpendicular to the axis and along the axis direction;

[0308] The viscous force is perpendicular to the axis and along the axis direction;

[0309] The viscous force is perpendicular to the axis and along the axis direction;

[0310] The viscous force is perpendicular to the axis and along the axis direction;

[0311] There is an interaction force between the flexible fiber model and the fluid model in the two-way fluid-structure interaction system; the calculation formula for the drag force of the interaction force between the flexible fiber model and the fluid model is:

[0312] (22);

[0313] Among them, the momentum exchange coefficient between the fluid and the fiber phase is expressed as:

[0314] (23);

[0315] The drag coefficient C D is expressed as:

[0316] (24);

[0317] The calculation formula for the Reynolds number Re is:

[0318] (25);

[0319] Wherein, is the fluid velocity;

[0320] is the flexible fiber velocity;

[0321] is the fluid porosity;

[0322] is the fluid density;

[0323] is the flexible fiber diameter;

[0324] is the fluid viscosity;

[0325] is the acceleration of gravity.

[0326] The coupling simulation process is as Figure 1 shown. A 90° straight pipe model is established in 3D modeling software, and the stp file is saved. The geometric model file is imported into Spaceclaim for fluid domain extraction, and the fluid domain file is saved. Then it is imported into fluent for mesh generation to obtain the fluid model. Set the CFD fluid model parameters, define the fluid domain as water, the inlet condition as velocity inlet, and the outlet as default pressure outlet. Select the turbulence model, and use the SIMPLE method to perform numerical simulation calculations within the given time step, with iterative loops. Define the number of calculation time steps as 1000, the time step as 0.001, and the maximum number of iterations as 20 to perform convergent calculation of fluid forces to obtain the fluid phase velocity field and pressure field information; in the discrete element DEM environment, establish a flexible fiber model, define the flexible fiber model parameters, including physical properties and mechanical properties, solve the fiber dynamics equation, calculate the contact force between fibers, and update the fiber velocity and position.

[0327] Through the divided finite element mesh, the flow field is spatially discretized, and trilinear interpolation is used to transfer data between fluid mesh nodes and discrete particle positions. The steps are as follows:

[0328] ① Locate the fluid mesh element where the fiber is located. According to the fiber coordinates locate the mesh element it is in. The vertex coordinates of the element are to , and find the eight vertex coordinates of the element and their corresponding fluid variable velocities;

[0329] ② Calculate the distance between each finite element mesh node and the particle center position. The relative distances (weighting coefficients) in each direction are:

[0330] ;

[0331] Among them, , calculate in the same way .

[0332] ③ Vertex weight calculation. The weight of each vertex is the product of the weight coefficients in three directions. For vertex , its weight calculation formula is:

[0333] ;

[0334] ④ According to the trilinear interpolation formula, sum the values of the eight vertices of the grid with weights to obtain the interpolation result at the fiber position. Let the values of the eight vertices be , the interpolation formula is:

[0335] ;

[0336] Expand it to:

[0337] ;

[0338] Among them, is the coordinate is the weight value of the vertex;

[0339] is the coordinate of the vertex.

[0340] In the two-way coupling, the fluid variables at time t are interpolated to the fiber position. It is also necessary to feedback the force on the fiber to the fluid grid. At this time, the interpolation weight is also determined by the weight of the trilinear interpolation. The force of each particle will be distributed to the surrounding eight nodes, and the force received by each node is the fiber force multiplied by the corresponding weight to obtain the interpolation result at the fiber position. In turn, at time t+1, time t+2, time t+3... until the calculation time is reached, the calculation ends.

[0341] Step 5: Use the two-way fluid-structure interaction system to simulate and physically experiment on the motion evolution form of a single flexible fiber in the flow field to obtain the minimum bending amount of the flexible fiber at any time (i.e., the bending distance between the two endpoints when the amplitude displacement of the flexible fiber is the largest) and the maximum displacement amount and the minimum displacement amount of the end point of the flexible fiber, and then calculate the indexes for evaluating the interaction force between the fluid and the flexible fiber, namely the bending degree of the flexible fiber in the flow field and the vibration amplitude of the flexible fiber in the flow field;

[0342] Step 6: Introduce the drag model correction factor , the strength of the force between the fluid and the flexible fiber is corrected by the evaluation index in step 5, and the simulation is carried out according to the number of fibers in the actual working condition, so as to judge the bending, winding, entanglement and loosening of the flexible fiber;

[0343] Introduce the drag force model correction factor into the source term of the force in the momentum control volume equation , the corrected drag force model expression is:

[0344] (26);

[0345] Among them, Momentum exchange coefficient between fluid and fiber phase;

[0346] The calculation formula of the curvature of the flexible fiber flow field is as follows:

[0347] (1);

[0348] Among them, S is the curvature of the flexible fiber flow field, is the initial length of the flexible fiber, in mm, is the minimum bending amount of the flexible fiber at any time, in mm;

[0349] The calculation formula of the amplitude of the flexible fiber flow field is as follows:

[0350] (2);

[0351] (3);

[0352] Among them, is the central axis of the amplitude of the flexible fiber flow field, in mm, is the maximum displacement of the end point of the flexible fiber, in mm, is the minimum displacement of the end point of the flexible fiber, in mm, is the amplitude of the flexible fiber flow field in the time series (the time series is the simulation time interval, and the simulation time interval in the present invention is set to 0.001 s), in mm, is the displacement of the flexible fiber in the time series, in mm.

[0353] Among them, in the simulation process, the motion posture animation of the flexible fiber in the flow field is decomposed and saved in frames of 0.001 s, and through the image processing technology, the minimum bending amount value in the whole time is measured, and it is substituted into formula (1) to obtain the value of the curvature of the flexible fiber flow field; through the image processing technology, each frame result is measured, and the maximum value and the minimum Value, substitute its two values into Equation (2) to obtain the central axis of the amplitude of the flexible fiber flow field , on the central axis of the amplitude of the flexible fiber flow field , measure the displacement of the flexible fiber at each continuous time interval in the time series , substitute it into Equation (3) to obtain the amplitude of the flexible fiber flow field. The schematic diagram is as Figure 2 shown.

[0354] The following uses specific embodiments to illustrate a flexible fiber fluid-structure interaction method of the present invention, specifically as follows:

[0355] To further verify the correctness of the model, the motion evolution form of a single flexible fiber flow field is simulated and the experiment is corrected and verified. In the CFD environment, the fluid domain is defined as water, the inlet condition is a velocity inlet, and the flow velocity is 4 m / s. The outlet is defaulted to a pressure outlet. In the DEM environment, the fiber diameter is defined as 60 μm, the fiber length is 38 mm, and the fiber density is 1750 kg / m 3 , the fiber generation speed is 0.05 kg / s, the fiber elastic modulus is 122 MPa, and the fiber Poisson's ratio is 0.3.

[0356] The experimental scheme is as Figure 4 shown. A flexible fiber fluid-structure interaction physical experimental platform is built. The experimental platform consists of a water source 3, an aluminum profile frame 4, a CCD high-definition camera 5, a background board 6, a lighting lamp 7, a water flow meter 8, a universal rotating pipe fixing clip 9, a highly transparent acrylic hair conveying pipe 10, and a recovery tank 11. Among them, the size of the aluminum profile frame 4 is 400 mm × 400 mm × 400 mm, the specification of the water flow meter 8 is a 4-inch thread G1 / 2 water flow sensor speedometer, the water pressure resistance is less than 1.75 MPa, the accuracy is ±5%, and the measuring temperature thermistor is 50 kΩ. The background board 7 is matte black to avoid the interference of wall reflection on the fiber flow field shooting effect. The designed experimental conditions are consistent with the numerical simulation calculation conditions. In order to prevent air from entering the pipeline during the experiment and affecting the fiber flow change, the pipeline interface is sealed with waterproof raw tape.

[0357] When performing numerical simulation calculations on the flexible fiber fluid-structure interaction, by changing the drag model correction factor of the interaction source in the momentum equation over a large range, a large number of numerical simulation calculations are carried out to make the values of the evaluation indexes of the mutual force of the fiber flow field, such as the curvature and amplitude, as close as possible to the values obtained from the experiment. Through continuous simulation calculations, when the drag model correction factor , the experimental and numerical simulation results are relatively close, and the error is within 5%. Moreover, in the experiment and numerical simulation calculations, the motion characteristics and evolution behaviors of the fiber in the flow field are relatively consistent. The comparison between the physical experiment and the numerical simulation process is as Figure 5As shown; under different drag force model correction factors, the comparison results of the vibration amplitude of the evaluation indexes between the physical experiment and the numerical simulation are as Figure 6 shown. Under different drag force correction factors, the comparison results of the curvature of the evaluation indexes between the physical experiment and the numerical simulation are as Figure 7 shown.

[0358] Therefore, when the drag force model correction factor is adopted in the present invention, the numerical simulation of the evolution characteristics of the flexible fiber flow field during fiber transportation is carried out. Under the working conditions of a fiber generation speed of 0.05 kg / s and a water flow speed of 2 m / s, during the transportation process, the fiber evolution behaviors at different moments are as Figure 8 shown. It can be seen from the figure that the fiber initially flows vertically along with the water. When passing through the 90° elbow group, due to the change in the flow field pattern, the fiber undergoes irregular bending and deformation evolution, laminar flow and turbulent flow phenomena occur in the flow field, and the fiber is involved in the turbulent flow and undergoes winding and knotting behaviors.

[0359] The present invention realizes the simulation visualization of the production process of flexible fiber transportation. Given a fiber generation speed of 0.05 kg / s and a water flow speed of 4 m / s under the initial conditions, the simulation is carried out, and model correction and experimental verification are carried out. The actual production process of the fiber washing and transportation process is simulated and reproduced, providing theoretical guidance for studying the phenomena of fiber bending, winding, entanglement and opening in the flow field. It shows that the present invention can be well applied to the research of engineering problems in actual production transportation, has certain academic significance and engineering value, and is expected to promote wider applications in the industry and promote the development of the industry towards science and efficiency.

Claims

1. A flexible fiber fluid-solid coupling method, characterized in that: It is applied to the dynamic simulation of the motion evolution characteristics of flexible fibers with large aspect ratio during water washing and conveying. Large aspect ratio refers to the aspect ratio of 100~500. The flexible fiber fluid-solid coupling method includes the following steps: Step 1: Construct a fluid model containing flexible fibers based on the fluid dynamics method; Step 2: The flexible fiber is described by using a discrete element method combined with a spherical column unit model to obtain a plurality of fiber spherical column units with the same length, and the plurality of fiber spherical column units with the same length are connected to each other to form a flexible fiber model; Step 3: Introduce the porosity of the fluid phase and establish a mass control body equation containing the fluid and the flexible fiber; introduce the interaction force between the fluid and the flexible fiber into the source term of the force in the momentum control body equation to obtain the momentum control body equation containing the source term of the flexible fiber force; Step 4: The fluid model including the flexible fiber and the flexible fiber model are interpolated and coupled by a trilinear method in a finite element volume method, and the fluid model including the flexible fiber and the flexible fiber model together constitute a bidirectional fluid-solid coupling system; there is an interaction force between the flexible fiber model and the fluid model in the bidirectional fluid-solid coupling system; Step 5: Use the bidirectional fluid-solid coupling system to simulate and conduct physical experiments on the motion evolution of a single flexible fiber in the flow field to obtain the minimum bending amount of the flexible fiber at any time. and the maximum displacement of the flexible fiber end point and minimum displacement , and then the indices for evaluating the interaction between the fluid and the flexible fiber are calculated, namely, the curvature of the flexible fiber flow field and the amplitude of the flexible fiber flow field; Step 6: Introduce the drag model correction factor into the source term of the force in the momentum control body equation , the strength of the force between the fluid and the flexible fiber is corrected by the evaluation index in step 5, and the simulation is performed according to the number of fibers in the actual working condition, so as to judge the bending, winding, entanglement and loosening of the flexible fiber; The calculation formula of the curvature of the flexible fiber flow field is as follows: （1）； Where S is the curvature of the flexible fiber flow field, is the initial length of the flexible fiber, is the minimum bending amount of the flexible fiber at any time; The calculation formula of the flexible fiber flow field amplitude is as follows: （2）； （3）； in, is the center axis of the flexible fiber flow field amplitude, is the maximum displacement of the end point of the flexible fiber, is the minimum displacement of the end point of the flexible fiber, is the amplitude of the flexible fiber flow field in time series, is the time series flexible fiber displacement.

2. A flexible fiber fluid-solid coupling method according to claim 1, characterized in that: The steps for constructing the fluid model in step 1 are as follows: Step 1.1: Based on the physical geometric model of the transmission pipeline, the fluid domain is extracted and the fluid grid is obtained by dividing it using the finite element method. The volume fraction occupied by the flexible fiber in the fluid grid is introduced to obtain the fluid model. Step 1.2: Discretize the fluid dynamics equations based on the finite element method to obtain a nonlinear equation with grid node parameter values ​​as solution quantities; solve the nonlinear equations to obtain the velocity and pressure changes of the flow field on the grid nodes; Step 1.3: Based on the fluid grid, a fluid model including fibers is established.

3. A flexible fiber fluid-solid coupling method according to claim 2, characterized in that: In the flexible fiber model of step 2, each fiber ball-column unit is composed of a cylinder and spheres at both ends, and adjacent fiber ball-column units are connected by spherical elastic parts.

4. A flexible fiber fluid-solid coupling method according to claim 3, characterized in that: The flexible fiber model in step 2 includes a group of dynamic equations of the force exerted by the fluid on the flexible fiber, and the group of dynamic equations of the force exerted by the fluid on the flexible fiber includes a force translation equation and a torque rotation equation; The force translation equation is described as: （4）； The torque rotation equation is described as: （5）； in, is the mass of a single flexible fiber; For the The translation speed of the fiber node ball; is the angular velocity vector of the node ball; is the acceleration due to gravity; For time; is the moment of inertia; Flexible fiber and contact force; is the total drag of the fluid on the flexible fiber; Flexible fiber and The tangential moment of Flexible fiber and The normal moment of .

5. A flexible fiber fluid-solid coupling method according to claim 4, characterized in that: The calculation formula is: （6）； in, For the The drag force of a flexible fiber; is the fluid unit volume; is the number of flexible fibers in the fluid unit; is the interpolation weight function; is the location of the flexible fibers; is the center coordinate of the fluid grid.

6. A flexible fiber fluid-solid coupling method according to claim 5, characterized in that: Based on the finite element volume unit method, equations (4) and (5) are solved as follows: First, equations (4) and (5) are simplified to: （7）； （8）； In the formula, is the acceleration of fiber i; is the moment of inertia of fiber i, is the angular acceleration of fiber i; is the force on fiber i, is the resultant external torque; Then integrate to get the velocity of the fiber at time t: （9）； In the formula, For fiber i The acceleration of the moment; For fiber i Angular acceleration at time Perform integration again to update the fiber motion posture: （10）； In the formula, To solve the displacement of fiber i, For fiber i displacement of time; To obtain the rotation angle of fiber i, For fiber i The rotation angle at the moment.

7. A flexible fiber fluid-solid coupling method according to claim 6, characterized in that: The mass control body equation containing the fluid and the flexible fiber in step 3 and the momentum control body equation containing the flexible fiber force source term are shown in equations (11) and (12), respectively; （11）； The momentum conservation equation including the flexible fiber force source term is described as: （12）； in, is the fluid porosity; is the fluid density; is the fluid velocity; is the fluid pressure; is the Hamiltonian operator; is the fluid viscosity; is the acceleration due to gravity; is the stress tensor; is the reaction force of the flexible fiber on the fluid; Fluid Porosity The calculation formula is: （13）； in, is the volume of the kth spherical-column unit fiber in the control volume; V is the volume of the control volume; is the interpolation weight function; is the position of the fiber; is the center coordinate of the fluid grid.

8. A flexible fiber fluid-solid coupling method according to claim 7, characterized in that: In step 2, the flexible fiber is described by using the discrete element method combined with the spherical column unit model, which means that the physical property parameters and mechanical property parameters are used for description. The physical property parameters are fiber diameter, fiber length, fiber density, fiber number, fiber elastic modulus and fiber Poisson's ratio. The mechanical property parameters are: spherical column unit normal force and tangential force , the normal force at the node sphere , tangential force , Torsional torque , Bending moment , Total torsional moment Combined bending moment ; The fiber diameter is 60~120μm, the fiber length is 5~100mm, and the fiber density is 750~2000kg / m 3 , the number of fibers is 1~9999, the fiber elastic modulus is 100~260MPa, and the fiber Poisson's ratio is 0.2~0.35; The calculation formulas for each mechanical property parameter are as follows: （14）； （15）； （16）； （17）； （18）； （19）； （20）； （21）； in, is the elastic modulus of the flexible fiber; is the node ball damping coefficient, the value is 0.035; is the shear modulus of the flexible fiber, , v is Poisson’s ratio; is the cross-sectional area of ​​the spherical column unit; is the polar moment of inertia of the spherical cylinder unit section, is the radius of the spherical cylinder unit; and are the normal displacement of the spherical-cylindrical unit and the tangential displacement of the spherical-cylindrical unit respectively; and are the normal stiffness and shear stiffness at the node ball, respectively; is the distance between the node balls; and are the torsional stiffness and bending stiffness at the node ball, respectively; For quality and moment of inertia The relative normal velocity between the two node balls; For quality and moment of inertia The relative tangential velocity between the two node balls; For quality and moment of inertia The relative torsional angular velocity between the two node balls; For quality and moment of inertia The relative bending angular velocity between the two node balls.

9. A flexible fiber fluid-solid coupling method according to claim 8, characterized in that: The calculation formula for the interaction drag force between the flexible fiber model and the fluid model in the bidirectional fluid-solid coupling system described in step 4 is: （22）； Among them, the momentum exchange coefficient between fluid and fiber is The expression is: （23）； Drag coefficient C D The expression is: （24）； The Reynolds number Re calculation formula is: （25）； in, is the fluid velocity; is the flexible fiber speed; is the fluid porosity; is the fluid density; is the flexible fiber diameter; is the fluid viscosity.

10. A flexible fiber fluid-solid coupling method according to claim 9, characterized in that: Introducing the drag model correction factor into the source term of the force in the momentum control body equation , the modified drag model expression is: （26）。 11. A flexible fiber fluid-solid coupling method according to claim 10, characterized in that: In step 4, the interpolation coupling calculation is performed by the trilinear method in the finite element volume method. Specifically, equations (7) and (8) are rewritten as general formulas: （27）； Among them, the fluid density vector including porosity It is expressed as: （28）； Three component convection terms in the spatial coordinate system in the fluid model , and Respectively expressed as: （29）； （30）； （31）； Three component viscous fluxes in the spatial coordinate system of the fluid model , and It is expressed as: （32）； （33）； （34）； is the force source term, including the interaction force between the fiber and the fluid: （35）； in, is the fluid velocity in the space coordinate system Components on the axis; is the fluid velocity in the space coordinate system Components on the axis; is the fluid velocity in the space coordinate system Components on the axis; For the fluid dynamics equations Convection term on axis; For the fluid dynamics equations Convection term on axis; For the fluid dynamics equations Convection term on axis; For the fluid dynamics equations Viscous flux on the axis; For the fluid dynamics equations Viscous flux on the axis; For the fluid dynamics equations Viscous flux on the axis; is the fluid porosity, is the fluid pressure, is the fluid viscosity tensor; , , They are Force source terms in three directions of space; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction; is the viscous force perpendicular to axis, and along Axis direction.

Images