A hybrid fiber modeling method considering different distribution densities

Abstract

The present application belongs to the technical field of simulation modeling of concrete, and specifically provides a hybrid fiber modeling method that takes into account different distribution densities. By randomly generating fibers, combining rough judgment and precise judgment of the spatial intersection of fibers, supplemented by the translation and rotation of fibers, and using spatial grids to calculate the fiber density of different specimen heights, three-dimensional spatial hybrid fibers are quickly and randomly generated; this method considers the mixing of straight fibers and end hook fibers in the process of establishing a hybrid fiber geometric model, and improves the efficiency of the hybrid fiber model through the translation and rotation of fibers. At the same time, affected by the actual construction process of the project, the unevenness of fiber density caused by the change of fiber density with the specimen height is taken into account, thereby improving the accuracy of the model simulation results.

Description

Technical Field

[0001] The present application belongs to the technical field of simulation modeling of concrete, and specifically provides a hybrid fiber modeling method considering different distribution densities. Background Art

[0002] At present, the study of the mechanical bridging effect of meso-scale fibers in concrete is mainly carried out through experiments and numerical simulation methods, and the spatial distribution of fibers is one of the main factors affecting the mechanical properties of fiber concrete.

[0003] In establishing the hybrid fiber geometric model, the spatial modeling is difficult due to the different geometric shapes of the hybrid fibers. Especially for the mixing of straight fibers and end hook fibers, it involves the determination and modeling of lines and surfaces, which makes the spatial stacking of the two very complicated, and it is difficult to establish a hybrid effect that conforms to the actual situation. In addition, in actual engineering, affected by the construction process, the fiber density may change with the height of the specimen, and in the traditional modeling process, the unevenness of the fiber density is often ignored, and the fiber modeling is mainly based on straight fibers. The traditional modeling method is slow and has poor results, resulting in low accuracy of the simulation results. Summary of the invention

[0004] The present application provides a hybrid fiber modeling method that takes into account different distribution densities to solve the above-mentioned problems existing in the prior art, namely, the hybrid fiber model in the prior art is inefficient and cannot consider the mixing of straight fibers and end hook fibers. In actual engineering, affected by the construction process, the height of the concrete specimen may also affect the density of the hybrid fibers, thereby resulting in the problem of low accuracy of the simulation results.

[0005] The present application provides a hybrid fiber modeling method considering different distribution densities, the method comprising:

[0006] Determine the geometric dimensions of the hybrid fiber concrete, construct the hybrid fiber concrete model, mesh the hybrid fiber concrete model, and number the meshes;

[0007] Obtain the fiber volume content, fiber geometric properties and volume ratio of straight fibers to end hook fibers of hybrid fiber concrete, and determine the required number of straight fibers and end hook fibers;

[0008] According to the meshed hybrid fiber concrete model, the geometric position information of the existing end hook fibers and straight fibers is determined;

[0009] According to the geometric position information of the end hook fiber and the straight fiber, the fiber is generated at the random position in the hybrid fiber concrete model, and the grid number of the new fiber is recorded;

[0010] Determine the generation state of new fibers in the hybrid fiber concrete model, judge the stacking state of the generated fibers and the existing fibers, perform spatial rotation and translation processing on the spatially overlapping fibers to generate new fibers, and determine the new fiber positions;

[0011] According to the new fiber position, the generation state of the new fiber is determined again, and the stacking state of the new fiber and the existing fiber is judged again; according to the meshed hybrid fiber concrete, the fiber density of the fibers at different heights in the non-overlapping state is obtained, and when it exceeds the actual density, the fiber is generated again at a random position in the hybrid fiber concrete model;

[0012] Repeat the process of generating non-overlapping fibers. When the number of straight fibers and end hook fibers reaches the required number, the modeling is completed.

[0013] Optionally, the obtaining of the fiber volume content, fiber geometric properties, and volume ratio of straight fibers to end hook fibers of the hybrid fiber concrete to determine the required number of straight fibers and end hook fibers specifically includes:

[0014] According to the geometric properties of the fiber, the fiber volume is obtained:

[0015] According to the fiber volume, volume content and the volume ratio of straight fibers to end hook fibers, the number of straight fibers and end hook fibers required in the concrete matrix is ​​obtained.

[0016] Optionally, the determining of the geometric position information of the existing end hook fibers and straight fibers according to the meshed hybrid fiber concrete model specifically includes:

[0017] According to the gridded hybrid fiber concrete, the midpoint coordinates of the end hook fiber and the straight fiber are determined;

[0018] The geometric position information of the end hook fiber and the straight fiber is determined according to the midpoint coordinates, and the geometric position information includes the end point coordinates of the straight fiber and the corner point and end point coordinates of the end hook fiber.

[0019] Optionally, determining the generation state of the new fiber in the hybrid fiber concrete model specifically includes:

[0020] Determine whether the fiber coordinates extend beyond the edge of the concrete specimen;

[0021] If so, it is determined that the generated state in the hybrid fiber concrete is an ungenerated state;

[0022] If not, the generating state is determined to be a generating new fiber state.

[0023] Optionally, the determining of the stacking state of the fibers in the generated state and the existing fibers specifically includes:

[0024] Get the mesh number of the generated fiber;

[0025] Determine whether the grid where the new fiber is located has the same grid number as the existing fiber;

[0026] If so, the repeated fiber number is recorded, and a determination is made as to whether the new fiber spatially overlaps with the repeated fibers of the grid;

[0027] If not, the new fiber and the generated fiber are not overlapped;

[0028] If the new fiber overlaps with the geometric space of the mesh repeating fibers, the new fiber overlaps with the already generated fibers;

[0029] Otherwise, the new fibers do not overlap with the already generated fibers.

[0030] Optionally, performing spatial rotation and translation processing on the spatially overlapping fibers to generate new fibers specifically includes:

[0031] Rotate and translate the fibers in space;

[0032] When the fiber coordinates fall within the concrete specimen, new fibers are obtained.

[0033] Optionally, the method of obtaining fiber densities of fibers at different heights in a non-overlapping state according to the meshed hybrid fiber concrete, and when the fiber densities exceed the actual density, generating fibers again at random positions in the hybrid fiber concrete model, specifically includes:

[0034] According to the gridded hybrid fiber concrete, determine the grid height where the center point of each equally divided fiber is located;

[0035] According to the grid height, the volume sum of all fiber equal segments in the grid of the same height is obtained;

[0036] According to the volume sum of all fiber equal segments, the fiber density at different heights is obtained;

[0037] Determine whether the fiber density at the height where the fiber is located exceeds the actual density, wherein the actual density is 10% of the allowable density;

[0038] If so, the fiber density state is exceeded, and fibers are generated again at random positions within the hybrid fiber concrete model;

[0039] If not, the fiber density state is an unsaturated state, and the fiber information of the unsaturated density state is recorded.

[0040] Compared with the prior art, the beneficial effects of the present application are as follows: the present application provides a hybrid fiber modeling method that takes into account different distribution densities. The method grids the hybrid fiber concrete, combines rough and precise judgments to determine the stacking state of the hybrid fibers, and quickly adjusts the fiber position in space through translation and rotation processing, thereby improving the modeling efficiency and modeling effect of the hybrid fibers. In addition, the present application also takes into account the unevenness of the fiber density caused by the change of fiber density with the height of the specimen, and obtains the fiber density of different grid heights through spatial grids, thereby greatly improving the accuracy of the simulation results. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0042] Figure 1 A flow chart of a hybrid fiber modeling method considering different distribution densities provided in an embodiment of the present application Figure 1 ;

[0043] Figure 2 A flow chart of a hybrid fiber modeling method considering different distribution densities provided in an embodiment of the present application Figure 2 ;

[0044] Figure 3 It is a schematic diagram of the modeling results of the fiber density varying with height in the case of a mixture of end hooks and straight fibers provided in an embodiment of the present application;

[0045] The above drawings have shown clear embodiments of the present application, which will be described in more detail later. These drawings and text descriptions are not intended to limit the scope of the present application in any way, but to illustrate the concept of the present application to those skilled in the art by referring to specific embodiments. DETAILED DESCRIPTION

[0046] In order to make the purpose, technical solutions and advantages of this application clearer, the technical solutions in this application will be clearly and completely described below in conjunction with the drawings in this application. Obviously, the described embodiments are part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0047] The terms "first", "second", "third", "fourth", etc. (if any) in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the numbers used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in sequences other than those illustrated or described herein, for example.

[0048] In the embodiments of the present application, the words "exemplary" or "for example" are used to indicate examples, illustrations or descriptions. Any embodiment or design described as "exemplary" or "for example" in the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "exemplary" or "for example" is intended to present related concepts in a specific way.

[0049] At present, due to the complexity of hybrid fibers, in establishing the hybrid fiber geometric model, especially for the mixing of straight fibers and end hook fibers, the spatial stacking is very complicated, which will lead to low efficiency and poor effect of the mixing effect. The traditional modeling process often ignores the unevenness of fiber density, resulting in low accuracy of simulation results.

[0050] In response to the above problems, the present application provides a hybrid fiber modeling method that takes into account different distribution densities. This method grids the hybrid fiber concrete, combines rough and precise judgments to determine the stacking state of the hybrid fibers, and quickly adjusts the fiber position in space through translation and rotation processing, thereby improving the modeling efficiency of the hybrid fibers. In view of the unevenness of fiber density caused by the change of fiber density with the height of the specimen, the fiber density at different grid heights is obtained through spatial grids, thereby greatly improving the accuracy of the simulation results.

[0051] The technical solution of the present application and how the technical solution of the present application solves the above-mentioned technical problems are described in detail below with specific embodiments. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments. The embodiments of the present application will be described below in conjunction with the accompanying drawings.

[0052] Figure 1 This is a schematic diagram of a modeling method for mixed fibers with different distribution densities provided in an embodiment of the present application. Figure 1 .like Figure 1 As shown, this embodiment shows a hybrid fiber modeling method considering different distribution densities, including:

[0053] S101: Determine the geometric dimensions of hybrid fiber concrete.

[0054] Among them, the lengths of the three sides are x, y, and z respectively, and the volume of the concrete matrix is:

[0055] V=x·y·z.

[0056] The present application does not impose any special restrictions on the specific implementation method of obtaining the geometric dimensions of the hybrid fiber concrete, as long as the geometric dimensions of the hybrid fiber concrete can be obtained.

[0057] S102: Meshing the hybrid fiber concrete according to the geometric dimensions and numbering the meshes.

[0058] The hybrid fiber concrete is meshed to facilitate the subsequent generation of new fibers, and the meshes are numbered to determine the stacking state of the generated new fibers.

[0059] In this step, the concrete specimen space can be divided into small grids, divided into Nivx, Nivy, and Nivz parts along the X, Y, and Z directions respectively, and the grids are numbered in sequence along the X, Y, and Z directions.

[0060] For example, (1, 1, 1) is numbered 1,…, (Nivx, 1, 1) is numbered Nivx,…, (Nivx, Nivy, 1) is numbered Nivx×Nivy,…, (Nivx, Nivy, Nivxz) is numbered Nivx×Nivy×Nivz.

[0061] S103: According to the fiber characteristics, the required number of straight fibers and end hook fibers is obtained.

[0062] The fiber characteristics may include, for example, fiber volume content, fiber geometric properties, and the volume ratio of straight fibers to end hook fibers.

[0063] S104: Determine geometric information of the end hook fiber and the straight fiber.

[0064] The geometric information may include, for example, the midpoint coordinates of the end hook fiber and the straight fiber, the coordinates of the straight fiber, the plane where the end hook fiber is located, the end point coordinates of the straight fiber, and the corner point and end point coordinates of the end hook fiber.

[0065] S105: Generate fibers at random positions in space.

[0066] After obtaining the geometric information of the end hook fiber and the straight fiber, a spatial rotation and translation process is performed according to the geometric information of the fiber to generate fibers at random positions in space.

[0067] You can choose random angles (γ, θ, φ) around the XYZ axes and random coordinates in space (x1, y1, z1) to rotate and translate the fiber space according to the rotation matrix.

[0068] S106: Determine the stacking state of the newly generated fibers and the fibers already generated in the concrete.

[0069] After obtaining new fibers, the stacking state of the new fibers and the generated fibers can be determined according to the position information of the fibers, which can improve the modeling effect of the mixed fibers.

[0070] In this step, a rough judgment is first performed. After finding the position containing the fiber in the grid, the current fiber number is recorded in the grid containing the fiber, and the stacking state of the new fiber and the generated fiber is judged according to the grid number; then a precise judgment is performed. If there are other fibers in the grid where the new fiber is located, the distance between the fibers in the grid is accurately calculated to determine whether they overlap.

[0071] S107: performing translation and rotation processing on the fiber overlap.

[0072] If the new fiber overlaps with other fibers in the concrete space, the fiber is translated and rotated in the opposite direction of the overlap. After the translation and rotation processing, step S106 is performed to determine the stacking state of the position-adjusted fiber and other fibers.

[0073] In this step, the fibers are slightly translated and rotated to reduce the number of times the fibers are regenerated, thereby reducing the number of cycles and increasing the speed of fiber generation in concrete.

[0074] S108: Determine fiber density at different fiber heights.

[0075] Among them, after generating new fibers, the grid height where the center point of each equally segmented fiber is located is judged to determine the fiber density at different fiber heights, which can ensure the accuracy of the fiber modeling simulation effect.

[0076] S109: When the number of straight fibers and the number of end hook fibers reach the required number, the modeling is completed.

[0077] The present embodiment provides a modeling method for hybrid fibers with different distribution densities. The hybrid fiber concrete is gridded, the fiber position is quickly adjusted in space through translation and rotation of the hybrid fibers, and the stacking state of the hybrid fibers is determined, thereby improving the modeling efficiency and modeling effect of the hybrid fibers. In addition, the present application also takes into account the unevenness of the fiber density caused by the change of the fiber density with the height of the specimen, and obtains the fiber density at different grid heights through spatial grids, thereby improving the accuracy of the simulation results.

[0078] Figure 2 A flow chart of a hybrid fiber modeling method considering different distribution densities provided in an embodiment of the present application Figure 2 .like Figure 2 As shown, this embodiment is Figure 1Based on the embodiment, a modeling method of hybrid fibers considering different distribution densities is described in detail. This embodiment shows a modeling method of hybrid fibers considering different distribution densities, including:

[0079] Step 1: Determine the geometric dimensions of the hybrid fiber concrete.

[0080] Step 2: Divide the concrete specimen space into small grids and number the grids.

[0081] Steps 1 and 2 are similar to the above-mentioned steps S101 and S102 and are not described in detail here.

[0082] Step 3, according to the fiber volume content, fiber geometric characteristics and the volume ratio of straight fibers to end hook fibers, the number of straight fibers and end hook fibers is calculated as follows:

[0083] Step 3.1, first, calculate the volume of a single straight fiber based on the fiber geometric properties:

[0084] V stright =π·d s 2 ·L s 4

[0085] Among them, d s is the diameter of the straight fiber, L s is the length of the straight fiber,

[0086] Step 3.2, calculate the volume of a single end hook fiber:

[0087] V hooked =π·d h 2 ·L h 4

[0088] Among them, d h is the diameter of the end hook fiber, L h is the total length of the end hook fiber,

[0089] Step 3.3, calculate the required fiber quantity in the concrete matrix based on the fiber volume content:

[0090]

[0091] Where ratio is the volume ratio of the fiber, N s is the number of straight fibers, and n is the ratio of hooked fibers to straight fibers.

[0092] Step 3.4, calculate the number of end hook fibers according to the ratio of straight fibers to end hook fibers:

[0093] N h =n·N s

[0094] Among them, N h is the number of end hook fibers.

[0095] Step 4, determine the position information of the end hook fiber and the straight fiber.

[0096] Among them, the midpoint of the end hook fiber and the straight fiber is selected at the coordinate origin (0, 0, 0), the straight fiber is located on the X-axis, the straight segment of the end hook fiber is located on the X-axis, and the plane of the end hook fiber is selected as the XOY plane. Then the coordinates of the end point of the straight fiber, the corner point and the end point of the end hook fiber can be expressed as follows:

[0097] Straight fiber endpoint coordinates: (-L s / 2,0,0),(L s / 2,0,0)

[0098] Coordinates of end hook fiber corner point 1: (-L hs / 2,0,0),(L hs / 2,0,0)

[0099] Coordinates of the 2nd corner point of the end hook fiber: (-L hs / 2-L i ˙cosθ,L i ˙sinθ,0)、(L hs / 2+L i ˙cosθ,L i ˙sinθ,0)

[0100] End hook endpoint coordinates: (-L hs / 2-L i ˙cosθ-L is , L i ˙sinθ,0)、(L hs / 2+L i ˙cosθ+L is , L i ˙sinθ,0)

[0101] Among them, L i L is the inclined length of the bent part of the end hook fiber, is It is the horizontal length of the bent part of the end hook fiber.

[0102] Step 5, using the spatial translation formula and the rotation matrix, generate fibers at random positions in space.

[0103] In this step, random angles (γ, θ, φ) around the XYZ axes and random spatial coordinates (x1, y1, z1) can be selected to rotate and translate the fiber space according to the spatial rotation matrix;

[0104] (x f ,yf ,z f )=(x1,y1,z1)+Dr˙(x',y',z')

[0105] Among them, (x', y', z') are the fiber node coordinates before coordinate transformation, (x f ,y f ,z f ) is the fiber node coordinate after coordinate transformation, and Dr is the rotation matrix:

[0106]

[0107] In this step, the fiber node coordinates before coordinate transformation are spatially rotated by using a rotation matrix to obtain the fiber node coordinates after coordinate transformation.

[0108] like Figure 3 As shown, by processing the translation and rotation of the hybrid fibers, the hybrid fibers can be quickly generated in the spatial hybrid fiber concrete model, which improves the modeling efficiency of the hybrid fibers and enables the establishment of the hybrid fiber concrete spatial model in more than 10 seconds.

[0109] Step 6: Check whether the coordinates of the two end points of the fiber exceed the edge of the concrete specimen.

[0110] If the coordinates are outside the concrete specimen, return to step 5 until the coordinates of the fiber are contained inside the specimen.

[0111] Step 7, find the position containing the fiber in the grid, and record the current fiber number in the grid containing the fiber.

[0112] In this step, a matrix from 1 to Nivx×Nivy×Nivz is created. If the fiber intersects a grid in space, the fiber number is recorded in the corresponding number.

[0113] For straight fibers, find the coordinate matrix (X cubic , Y cubic , Z cubic ), and then calculate the grid numbers in sequence:

[0114] V num =X cubic +Y cubic Nivx+Z cubic ·Nivx·Nivy

[0115] Among them, V num is the number of the grid where the fiber is located.

[0116] For the end hook fiber, divide it into a bent parallel segment, a bent oblique segment, and a middle segment, and find the grid position and its number in space in turn.

[0117] Step 8, roughly determine whether the new fiber overlaps with the existing fiber, and check whether the grid where the new fiber is located has the same grid number as the existing fiber. If it does, continue with the following steps; if not, jump to step 10.

[0118] This step is a rough judgment, which checks whether other fiber numbers already exist in the grid number where the new fiber is located. If not, no further judgment is required.

[0119] Step 9, accurately determine whether the new fiber overlaps with the existing fiber, record the repeated fiber numbers, and calculate in sequence whether the new fiber overlaps with the repeated fibers in the grid. Specifically:

[0120] If it is a straight fiber, calculate the fiber vector and two mutually perpendicular vectors perpendicular to the fiber, divide the fiber into 10 equal parts, and divide it into 12 equal parts along the vertical direction of the fiber according to the angle, and calculate the coordinates of each equal-division point along the length and each equal-division angle along the vertical direction, with the length being the fiber radius. A total of 120 points are used to determine whether the fibers are crossed or overlapped.

[0121] The coordinates of the 120 points are sequentially compared with the 120 points of the fiber containing the common grid to determine whether the distance is less than the arc length of the 12-square angle. If the distance between the dividing point of the new fiber and the dividing point of the existing fiber is less than the arc length of the 12-square angle, then go to step 10, otherwise go to step 11.

[0122] If it is an end hook fiber, the fiber is divided into a bent parallel section, a bent oblique line section, and a middle section. The middle section is treated as a straight fiber and divided into 10 sections along the length direction and the fiber surface is divided into 12 equal parts along the fiber cross section direction. At the same time, according to the length of the middle section, the bent oblique section and the bent parallel section are divided according to the length along the length direction, and the fiber surface is divided into 12 equal parts along the cross section direction.

[0123] Step 10, calculate the direction vector of the closest new fiber dividing point and the existing fiber dividing point, take the existing fiber dividing point as the starting point, translate the coordinates of the new fiber according to this direction vector, and each translation length is the fiber radius. After each translation, check whether the coordinates of the two end points of the fiber exceed the edge of the concrete specimen. If they do not exceed the edge of the concrete specimen, proceed to step 7-step 9, otherwise continue this step and continue to translate the fiber. If the translation is 10 times and the following steps are not entered, proceed to step 11.

[0124] Step 11, increase the initial rotation angle of the fiber in three directions by a random value within 10°, and rotate the fiber slightly. After the rotation is completed, check whether the coordinates of the end points on both sides of the fiber exceed the edge of the concrete specimen. If they do not exceed the edge of the concrete specimen, proceed to step 7-step 10, otherwise continue with this step. If the following step is not entered after 10 rotations, proceed to step 5-step 10.

[0125] Step 12, determine the grid height where the center point of each equally divided fiber segment is located, calculate the volume sum of all equally divided fiber segments in the grid at the same height, and determine whether the fiber volume at this height exceeds 10% of the allowable density. If the density of the fiber in the fiber-containing grid layer does not exceed 10% of the allowable density, continue with the following steps, otherwise repeat steps 5 to 10.

[0126] Figure 3 This is a schematic diagram of the modeling results of the fiber density changing with height in the case of a mixture of end hooks and straight fibers provided in this embodiment. Figure 3 As shown, since the density of the mixed fibers is affected by the fiber height, the density of the mixed fibers of the end hook fibers and the straight fibers gradually decreases along the direction of the arrow.

[0127] Step 13, record the fiber number and fiber type, the end point coordinates and corner point coordinates of the fiber space, and the grid number where the fiber is located, for subsequent fiber judgment.

[0128] Step 14: Repeat steps 4 to 13 until the number of straight fibers reaches N. s , the number of hook fibers reaches N h .

[0129] When the number of recorded straight fibers and end hook fibers reaches the number required by the calculation, the fiber random space modeling ends.

[0130] The present embodiment provides a hybrid fiber modeling method that takes into account different distribution densities. The hybrid fiber concrete is gridded, and the grid position of the newly generated fibers and the spatial distance between each point of the new fiber at the same grid fiber dividing points are determined, so as to roughly and accurately judge the stacking state of the hybrid fibers. The spatially overlapping fibers are fine-tuned by spatial rotation and translation, thereby improving the modeling efficiency and modeling effect of the hybrid fibers. At the same time, the non-uniformity of the fiber density caused by the change of the fiber density with the height of the specimen is also taken into account. The fiber density changes at different grid heights are obtained by spatial gridding, thereby improving the accuracy of the simulation results.

[0131] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0132] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.

Claims

1. A hybrid fiber modeling method considering different distribution densities, characterized in that: The method comprises: Determine the geometric dimensions of the hybrid fiber concrete, construct the hybrid fiber concrete model, mesh the hybrid fiber concrete model, and number the meshes; Obtain the fiber volume content, fiber geometric properties and volume ratio of straight fibers to end hook fibers of hybrid fiber concrete, and determine the required number of straight fibers and end hook fibers; According to the meshed hybrid fiber concrete model, the geometric position information of the existing end hook fibers and straight fibers is determined; According to the geometric position information of the end hook fiber and the straight fiber, the fiber is generated at the random position in the hybrid fiber concrete model, and the grid number of the new fiber is recorded; Determine the generation state of new fibers in the hybrid fiber concrete model, judge the stacking state of the generated fibers and the existing fibers, perform spatial rotation and translation processing on the spatially overlapping fibers to generate new fibers, and determine the new fiber positions; According to the new fiber position, the generation state of the new fiber is determined again, and the stacking state of the new fiber and the existing fiber is judged again; according to the meshed hybrid fiber concrete, the fiber density of the fibers at different heights in the non-overlapping state is obtained, and when it exceeds the actual density, the fiber is generated again at a random position in the hybrid fiber concrete model; The method of obtaining fiber densities of fibers at different heights in a non-overlapping state according to the meshed hybrid fiber concrete, and when the fiber densities exceed the actual density, generating fibers at random positions in the hybrid fiber concrete model again, specifically includes: According to the gridded hybrid fiber concrete, determine the grid height where the center point of each equally divided fiber is located; According to the grid height, the volume sum of all fiber equal segments in the grid of the same height is obtained; According to the volume sum of all fiber equal segments, the fiber density at different heights is obtained; Determine whether the fiber density at the height where the fiber is located exceeds the actual density, where the actual density is 10% of the allowable density; If so, the fiber density state is exceeded, and fibers are generated again at random positions within the hybrid fiber concrete model; If not, the fiber density state is unsaturated, and the fiber information of the unsaturated density state is recorded; Repeat the process of generating non-overlapping fibers. When the number of straight fibers and end hook fibers reaches the required number, the modeling is completed.

2. The method according to claim 1, characterized in that The method of obtaining the fiber volume content, fiber geometric characteristics, and volume ratio of straight fibers to end hook fibers of the hybrid fiber concrete and determining the required number of straight fibers and end hook fibers specifically includes: According to the geometric properties of the fiber, the fiber volume is obtained: According to the fiber volume, volume content and the volume ratio of straight fibers to end hook fibers, the number of straight fibers and end hook fibers required in the concrete matrix is ​​obtained.

3. The method according to claim 1, characterized in that The method of determining the geometric position information of the existing end hook fibers and straight fibers according to the meshed hybrid fiber concrete model specifically includes: According to the gridded hybrid fiber concrete, the midpoint coordinates of the end hook fiber and the straight fiber are determined; The geometric position information of the end hook fiber and the straight fiber is determined according to the midpoint coordinates, and the geometric position information includes the end point coordinates of the straight fiber and the corner point and end point coordinates of the end hook fiber.

4. The method according to claim 1, characterized in that The determining of the generation state of the new fiber in the hybrid fiber concrete model specifically includes: Determine whether the fiber coordinates extend beyond the edge of the concrete specimen; If so, it is determined that the generated state in the hybrid fiber concrete is an ungenerated state; If not, it is determined that the generating state is a generating new fiber state.

5. The method according to claim 1, characterized in that The determining of the stacking state of the fibers in the generated state and the existing fibers specifically includes: Get the mesh number of the generated fiber; Determine whether the grid where the new fiber is located has the same grid number as the existing fiber; If so, the repeated fiber number is recorded, and a determination is made as to whether the new fiber spatially overlaps with the repeated fibers of the grid; If not, the new fiber and the generated fiber are not overlapped; If the new fiber overlaps with the geometric space of the mesh repeating fibers, the new fiber overlaps with the already generated fibers; Otherwise, the new fibers do not overlap with the already generated fibers.

6. The method according to claim 1, characterized in that The performing spatial rotation and translation processing on the spatially overlapping fibers to generate new fibers specifically includes: Rotate and translate the fibers in space; When the fiber coordinates fall within the concrete specimen, new fibers are obtained.

Images