CN112700824B - Method for establishing fiber concrete mesoscopic model - Google Patents

Method for establishing fiber concrete mesoscopic model Download PDF

Info

Publication number
CN112700824B
CN112700824B CN202110040802.5A CN202110040802A CN112700824B CN 112700824 B CN112700824 B CN 112700824B CN 202110040802 A CN202110040802 A CN 202110040802A CN 112700824 B CN112700824 B CN 112700824B
Authority
CN
China
Prior art keywords
fiber
fibers
aggregate
model
distance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110040802.5A
Other languages
Chinese (zh)
Other versions
CN112700824A (en
Inventor
李厚民
柯俊宏
汪洋
吴克洋
黄笑宇
邓维超
李子毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hubei University of Technology
Original Assignee
Hubei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hubei University of Technology filed Critical Hubei University of Technology
Priority to CN202110040802.5A priority Critical patent/CN112700824B/en
Publication of CN112700824A publication Critical patent/CN112700824A/en
Application granted granted Critical
Publication of CN112700824B publication Critical patent/CN112700824B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like

Abstract

The invention provides a method for establishing a fiber concrete mesoscopic model, which comprises the following steps: establishing a random aggregate model, randomly generating fibers in the random aggregate model, acquiring related parameters of the randomly generated fibers and related parameters of all aggregates, calculating the related parameters of the fibers and the related parameters of each aggregate, judging whether intrusion phenomena do not occur between the fibers and each aggregate, if so, the fibers meet the requirements, and otherwise, rejecting and regenerating the fibers; and calculating whether the intrusion phenomenon is not generated between the randomly generated fibers meeting the requirements and all the fibers meeting the requirements in the model, if so, the fibers meet the requirements of the model, otherwise, removing and regenerating the fibers and continuing to perform the calculation of the steps, thereby establishing the fiber concrete mesoscopic model without the intrusion phenomenon. The calculation method is simple and consumes less time, and can ensure that no invasion phenomenon is generated between aggregates and fibers and between fibers in the established model.

Description

Method for establishing fiber concrete mesoscopic model
Technical Field
The invention belongs to the technical field of fiber concrete simulation, and particularly relates to a method for establishing a fiber concrete mesoscopic model.
Background
On a microscopical level, fiber concrete is considered to be a multiphase composite material consisting of coarse and fine aggregates, an aggregate interface layer, cement mortar, fibers, and an interface layer between fibers and mortar. In order to establish a more realistic numerical model which can reflect the distribution of the fibers and the aggregates in the fiber concrete, it is necessary to ensure that no invasion phenomenon exists between the aggregates and the fibers and between the fibers during modeling. When aggregate and fiber judgment is carried out, the Schopper schroer proposes to calculate the distance from the fiber to the centroid of the aggregate in a vector projection mode, but the method is not easy to expand to three dimensions and has a complex calculation mode. When the position relation between the fibers is carried out, the method proposes that the position relation of a circle with the fiber length as the diameter is compared to carry out early-stage rapid screening, and whether the fibers are intersected or not is judged by judging whether the fibers are spanned or not through a vector method. The method is complicated, has more calculation processes and takes long time for screening.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides the method for establishing the fiber concrete mesoscopic model, the method is simple in calculation method and less in time consumption, and the phenomenon that the aggregate and the fiber are not invaded in the established model can be ensured.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method for establishing a fiber concrete mesoscopic model comprises the following steps:
s1: establishing a random aggregate model, randomly generating fibers in the random aggregate model, acquiring related parameters of the randomly generated fibers and related parameters of all aggregates, and calculating the related parameters of the fibers and the related parameters of each aggregate to judge whether intrusion phenomenon does not occur between the randomly generated fibers and each aggregate, wherein if yes, the fibers meet the requirements, and if not, the fibers are removed and the related parameters of the fibers are regenerated;
s2: and calculating whether the intrusion phenomenon is not generated between the randomly generated fibers meeting the requirements in the step S1 and all the fibers meeting the requirements in the model, if so, the fibers meet the requirements of the model, otherwise, relevant parameters of the fibers are removed and regenerated, and the calculation of the steps is continued, so that the fiber concrete mesoscopic model without the intrusion phenomenon is established.
Further, on a microscopic layer, establishing an xoy plane coordinate system, and randomly generating related parameters of the fibers and related parameters of the aggregates;
wherein the relevant parameter of the fiber comprises the fiber length LlengthFiber start point coordinate (X)start,Ystart) The rotation angle theta of the fiber is an included angle between the fiber and the x axis;
the related parameters of the aggregate comprise the centroid coordinate (X) of the aggregateAggregate material,YAggregate material) Radius R of aggregateAggregate materialAnd an aggregate interface thickness T.
Further, the method for judging whether the invasion phenomenon is generated between the fiber and the aggregate is as follows:
a. calculating the coordinates (X) of the termination point of the randomly generated fibers according to the related parameters of the fibersend,Yend) The calculation formula is as follows: xend＝Xstart+Llength×cosθ；Yend＝Ystart+Llength×sinθ；
b. By randomly generating the coordinates (X) of the starting point of the fibresstart,Ystart) Coordinate of end point (X)end,Yend) And the coordinates of the centroid of the aggregate (X)Aggregate material,YAggregate material) Calculating the area S of a triangle formed by randomly generated fiber starting points, end points and an aggregate centroid;
c. a triangular area S formed by randomly generating a start point and an end point of the fiber and an aggregate centroid and a fiber length LlengthCalculating the distance L from the centroid of the aggregate to the fiberdistance；
d. If the distance L from the centroid of the aggregate to the fiberdistanceGreater than or equal to the corresponding aggregate radius RAggregate materialAnd the sum of the thickness T of the interface layer, judging that the fiber and the aggregate do not have invasion, continuously judging whether the fiber and all the aggregates do not have invasion through the steps, and if so, ensuring that the fiber meets the requirement; otherwise, rejecting the fiber and regenerating the related parameters of the fiber and returning to the step a for re-calculation.
Further, the calculation formula of the area S of the triangle is:
Figure BDA0002895764570000021
Ldistancethe calculation formula of (2) is as follows: l isdistance＝2S/Llength。
Further, in step S2, the acquired relevant parameters of the fiber include:
on a microscopic level, an xoy plane coordinate system is established, and the related parameters of a certain fiber meeting the requirements in the model comprise the length L of the fiber1lengthThe coordinates of the starting point of the fiber (X1)start,Y1start) The rotation angle theta of the fiber1Angle of rotation theta1Is the angle between the fiber and the x-axis;
randomly generating relevant parameters of the ith fiber in the model, wherein the relevant parameters respectively comprise the length Li of the ith fiberlengthCoordinates of the starting point of the ith fiber (Xi)start,Yitart) Rotation angle θ of ith fiberiRotation angle θ of ith fiberiIs the angle between the ith fiber and the x-axis.
Further, after the ith fiber and all the aggregates are obtained through the calculation in the step S1 and do not have invasion phenomenon, it is determined whether the ith fiber and all the fibers which have been in the model and meet the requirements do not have invasion phenomenon, and the determination method is as follows:
calculating the coordinates of the termination point of the ith fibre (Xi)end,Yiend) The calculation formula is as follows:
Xiend＝Xistart+Lilength×cosθi；
Yiend＝Yistart+Lilength×sinθi；
ii, calculating the area S of a triangle formed by the starting point of the ith fiber and the starting point and the ending point of a certain fiber in the modelstart；
Calculating the center point of the ith fiber, and the starting point and the ending point of a certain fiber existing in the modelArea S of trianglecenter；
Iv, calculating the termination point of the ith fiber and the area S of a triangle formed by the starting point and the termination point of a certain fiber in the modelend；
V, respectively calculating the distance Ls from the starting point of the ith fiber to a certain fiber in the modeldistanceDistance Lc from the midpoint of the ith fiber to a fiber in the modeldistanceAnd the distance Le from the termination point of the ith fiber to a certain fiber existing in the modeldistance；
Vi, if 2Lcdistance＝Lsdistance+LedistanceThen calculate the line segment X1start XiendSlope of (1) and line segment XiendX1endAnd judging whether the two slopes are equal to each other, and then calculating the line segment X1start XistartSlope of (1) and line segment Xistart X1endJudging whether the two slopes are equal, and if the two slopes are not equal, judging that the ith randomly generated fiber and a certain fiber existing in the model do not generate an invasion phenomenon; according to the method, whether the ith fiber generated randomly does not produce invasion phenomenon with all fibers existing in the model or not is judged, if yes, the ith fiber meets the requirement of the model; if the condition is not met, the relevant parameters of the ith fiber are regenerated, and the step is returned to calculate again.
Further, the calculation formula of the area is as follows:
Figure BDA0002895764570000031
Figure BDA0002895764570000032
Figure BDA0002895764570000041
according to area SstartArea ScenterSum area SendThe following can be determined:
Lsdistance＝2Sstart/L1length，Lcdistance＝2/ScenterL1length，Ledistance＝2Send/L1length。
compared with the prior art, the invention has the beneficial effects that: the calculation method is simple and consumes less time, and can ensure that no invasion phenomenon is generated between aggregates and fibers and between fibers in the established model, so that the established model is more in line with the actual situation, the simulation degree is high, and a foundation is provided for the subsequent simulation test and research of the concrete.
Drawings
FIG. 1 is a flow chart of a method for building a fiber concrete mesoscopic model according to an embodiment of the invention
FIG. 2 is a schematic view showing the non-invasion of fibers and aggregates in the embodiment of the present invention;
FIG. 3 is a schematic view showing invasion of fibers and aggregates in the embodiment of the present invention;
FIG. 4 is a schematic view of the non-intrusion of fibers into one embodiment of the invention;
FIG. 5 is a schematic view showing the invasion phenomenon of the fiber and the fiber in the embodiment of the present invention
Fig. 6 is a schematic structural diagram of a model established according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive efforts based on the embodiments of the present invention, shall fall within the scope of protection of the present invention.
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
The present invention is further illustrated by the following examples, which are not to be construed as limiting the invention.
The invention provides a method for establishing a fiber concrete mesoscopic model, which comprises the following steps:
s1: establishing a random aggregate model, randomly generating fibers in the random aggregate model, acquiring related parameters of the randomly generated fibers and related parameters of all aggregates, and calculating the related parameters of the fibers and the related parameters of each aggregate to judge whether intrusion phenomena do not occur between the randomly generated fibers and each aggregate, wherein if yes, the fibers meet the requirements, and if not, the fibers are removed and the related parameters of the fibers are regenerated;
specifically, as shown in fig. 1, the random fiber feeding method is based on a traditional random aggregate model. Firstly, storing relevant parameters of each aggregate put in randomly, establishing an xoy plane coordinate system on a microscopic level, and randomly generating relevant parameters of fibers and the aggregates; wherein the relevant parameter of the fiber comprises the fiber length LlengthFiber start point coordinate (X)start,Ystart) The rotation angle theta of the fiber is an included angle between the fiber and the x axis; the related parameters of the aggregate comprise the centroid coordinate (X) of the aggregateAggregate material,YAggregate material) Radius R of aggregateAggregate materialAnd an aggregate interface thickness T.
Then judging whether intrusion phenomenon is generated between the randomly thrown fibers and the aggregates, wherein the specific method comprises the following steps:
a. calculating the coordinates (X) of the termination point of the fiber according to the related parameters of the fiberend,Yend) The calculation formula is as follows: xend＝Xstart+Llength×cosθ；Yend＝Ystart+Llength×sinθ；
b. Through the coordinates (X) of the starting point of the fibrestart,Ystart) Coordinate of end point (X)end,Yend) And the coordinates of the centroid of the aggregate (X)Aggregate material,YAggregate material) Calculating the area S of a triangle formed by the fiber starting point, the fiber ending point and the aggregate centroid,
Figure BDA0002895764570000051
c. a triangular area S formed by the fiber starting point, the fiber ending point and the aggregate centroid and a fiber length LlengthCalculating the distance L from the centroid of the aggregate to the fiberdistance，LdistanceThe calculation formula of (2) is as follows: l isdistance＝2S/Llength；
d. If the distance L from the centroid of the aggregate to the fiberdistanceGreater than or equal to the corresponding aggregate radius RAggregate materialAnd the sum of the thickness T of the interface layer, judging that the fiber and the aggregate do not have invasion, continuously judging whether the fiber and all the aggregates do not have invasion through the steps, and if so, ensuring that the fiber meets the requirement; otherwise, rejecting the fiber and regenerating the related parameters of the fiber and returning to the step a for re-calculation.
For convenience of explanation, a specific example is taken for explanation, as shown in fig. 2, the starting point of the first fiber generated randomly is point B and its coordinates are (X1, Y1), the coordinates of the ending point C of the first fiber (X2, Y2) can be obtained by four basic parameters of the first fiber, and in addition, the coordinates of the centroid a of the aggregate according to the parameters of the aggregate are (X3, Y3), and then the formula is used for obtaining the coordinates of the centroid a of the aggregate
Figure BDA0002895764570000052
Calculate the area S of triangle ABCABC。
Furthermore, the area S due to triangle ABCABCCan also be passed through
Figure BDA0002895764570000053
Thus, the distance L from the centroid of the aggregate to the fiber can be determineddistanceI.e. length of AD, Ldistance＝2SABCand/BC, wherein the length of BC is the length of the first fiber. If the length of AD is greater than the radius R of the aggregateAggregate materialIf the sum of the thickness T of the interface layer with the aggregate is judged to be the same, the fiber and the aggregate are separated from each other, and the intrusion phenomenon between the fiber and the aggregate is not generated. Such as byAnd if the condition that the generated first fiber meets the condition that the first fiber is not invaded by all aggregates is obtained through calculation, generating the 2 nd fiber, and judging the invasion between the 2 nd fiber and the aggregates according to the method, and if the condition is met, judging the invasion between the fibers. When the first fiber and the aggregate are penetrated, as shown in FIG. 3, the distance L from the centroid A of the aggregate to the fiber BC in FIG. 3 is determined by the above methoddistanceAnd if the length of the AD is smaller than the sum of the radius of the corresponding aggregate and the thickness of the aggregate interface layer, determining that an invasion phenomenon occurs between the fibers and the aggregate, rejecting the fibers, randomly generating related parameters of the first fiber again, and continuously calculating according to the method.
S2: and calculating whether the intrusion phenomenon is not generated between the randomly generated fibers meeting the requirements in the step S1 and all the fibers meeting the requirements in the model, if so, the fibers meet the requirements of the model, otherwise, relevant parameters of the fibers are removed and regenerated, and the calculation of the steps is continued, so that the fiber concrete mesoscopic model without the intrusion phenomenon is established.
In this embodiment, if a first fiber is randomly generated and satisfies that the first fiber and all aggregates do not have an intrusion phenomenon, four random parameters of the first fiber are obtained, where the four random parameters respectively include a length L of the first fiber1lengthFirst fiber start point coordinate (X1)start,Y1start) The rotation angle theta of the first fiber1Angle of rotation theta1Is the angle between the fiber and the x-axis;
then acquiring four random parameters of the ith randomly generated fiber, wherein i is more than or equal to 2 and respectively comprises the fiber length LilengthFiber start point coordinate (Xi)start,Yitart) The rotation angle theta of the fiberiAngle of rotation thetaiIs the angle between the fiber and the x-axis;
after the ith fiber and all aggregates are obtained through calculation and do not generate the intrusion phenomenon, whether the intrusion phenomenon is generated between the ith fiber and the fiber meeting the requirements in the model is judged, in this embodiment, whether the intrusion phenomenon is generated between the ith fiber and the first fiber is judged, and the specific method is as follows:
calculating the coordinates of the termination point of the ith fibre (Xi)end,Yiend) The calculation formula is as follows:
Xiend＝Xistart+Lilength×cosθi；
Yiend＝Yistart+Lilength×sinθi；
ii, calculating the area S of a triangle formed by the starting point of the ith fiber, the starting point of the first fiber and the ending point of the first fiberstart，
Figure BDA0002895764570000061
Calculating the area S of a triangle formed by the midpoint of the ith fiber, the starting point of the first fiber and the ending point of the first fibercenter，
Figure BDA0002895764570000071
Iv, calculating the area S of a triangle formed by the termination point of the ith fiber, the starting point of the first fiber and the termination point of the first fiberend，
Figure BDA0002895764570000072
V, respectively calculating and calculating the distance Ls from the starting point of the ith fiber to the 1 st fiberdistanceDistance Lc from the midpoint of the ith fiber to the 1 st fiberdistanceAnd the distance Le from the starting point of the ith fiber to the 1 st fiberdistanceWherein, Lsdistance＝2Sstart/L1length，Lcdistance＝2/ScenterL1length，Ledistance＝2Send/L1length；
Vi if two fibres do not invadeThen a right-angle trapezoid can be formed by the distance from the starting point of the ith fiber to the 1 st fiber, and the length of the ith fiber. The distance from the starting point of the ith fiber to the 1 st fiber and the distance from the starting point of the ith fiber to the 1 st fiber are parallel edges, the length of the ith fiber is the oblique edge of a right trapezoid, the distance from the middle point of the ith fiber to the 1 st fiber is the median line of the trapezoid, and according to the theorem of the median line of the trapezoid, 2Lc existsdistance＝Lsdistance+LedistanceConsidering the special case that the first fiber is perpendicular to the ith fiber, the line segment X1 is calculatedstart XiendSlope of (1) and line segment Xiend X1endAnd judging whether the two slopes are equal or not, and then calculating a line segment X1startXistartSlope of (1) and line segment XistartX1endAnd determining whether the slopes are equal to each other, and if neither of the slopes is equal to each other, determining that the first fiber and the ith fiber are not invaded. According to the method, whether the ith fiber generated randomly and all fibers in the model do not produce invasion is judged, if yes, the ith fiber meets the requirements of the model; if the condition is not met, regenerating the relevant parameters of the ith fiber and returning to the step i for calculation again.
Specifically, as shown in fig. 4, AC is a newly generated fiber, point B is the midpoint of AC, and DE is an existing fiber. And (3) making a perpendicular line of DE (Dee) through the point A, taking a foot F, making a perpendicular line of DE through the point B, taking a foot G, making a perpendicular line of DE through the point C, and taking a foot H. Still according to the formula
Figure BDA0002895764570000073
The areas of the triangles ADE, BDE and CDE are obtained, the coordinates of three vertexes of any one of the three triangles are respectively represented by (X1, Y1), (X2, Y2) and (X3, Y3) in the calculation formula, and then the distances from the starting point, the middle point and the ending point of the newly generated fiber to the existing fiber DE, namely the lengths of AF, BG and CH, are obtained. If the quadrangle formed by the points AFHC is a right trapezoid and BG is the right trapezoidThe median line of the graph, 2 × BG ═ AF + CH, according to the median theorem. Therefore, it is considered that the above condition is satisfied if the two fibers do not have the intrusion phenomenon. However, consider the situation shown in FIG. 5, where the nascent fiber AC is perpendicular to the existing fiber DE and the point C is on the fibril DE. The conditions of the right trapezoid and the median line described above are also met for this particular case, but the two fibers still have an intersection and therefore need to be excluded. The exclusion method is as follows: whether the slopes from the starting point A of the newly-generated fiber AC to the starting point E and the ending point D of the existing fiber DE are equal or not is determined, whether the slopes from the ending point C of the newly-generated fiber AC to the starting point E and the ending point D of the existing fiber DE are equal or not is determined, and if the DC slope is equal to the CE slope in FIG. 5, it is determined that the newly-generated fiber and the existing fiber are invaded, and the newly-generated fiber and the existing fiber need to be discarded; if not, the newly generated fiber is judged to have no invasion phenomenon with the existing fiber.
In this embodiment, a conventional random aggregate model is constructed, the size of the random aggregate model is selected to be 150mm by 150mm, the aggregates are selected to be 5-20mm continuous gradation, the particle size range is divided into two types, namely 5-10mm and 10-20mm, the number of the aggregates with the particle size of 5-10mm is 40, and the number of the aggregates with the particle size of 10-20mm is 13 according to the wararan formula. And the fibers were randomly generated in the aggregate, and 30mm long fibers were used, and the number of the fibers was 25. The simulation modeling is carried out by adopting the method, the result after modeling is shown in figure 6, the fiber and the aggregate have no invasion, and the fiber also have no invasion, thereby meeting the mesoscopic modeling requirement of the fiber concrete.
While the invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention.

Claims (5)

1. A method for establishing a fiber concrete mesoscopic model is characterized by comprising the following steps:
s1: establishing a random aggregate model, randomly generating fibers in the random aggregate model, acquiring related parameters of the randomly generated fibers and related parameters of all aggregates, and calculating the related parameters of the fibers and the related parameters of each aggregate to judge whether intrusion phenomenon does not occur between the randomly generated fibers and each aggregate, wherein if yes, the fibers meet the requirements, and if not, the fibers are removed and the related parameters of the fibers are regenerated;
s2: calculating whether intrusion phenomena are not generated between the randomly generated fibers meeting the requirements in the step S1 and all the existing fibers meeting the requirements in the model, if so, the fibers meet the requirements of the model, otherwise, relevant parameters of the fibers are removed and regenerated, and the calculation of the steps is continued, so that a fiber concrete mesoscopic model without intrusion phenomena is established;
in step S1, establishing a xoy plane coordinate system on a microscopic layer, and randomly generating related parameters of the fibers and related parameters of the aggregates;
wherein the relevant parameter of the fiber comprises the fiber length LlengthFiber start point coordinate (X)start,Ystart) The rotation angle theta of the fiber is an included angle between the fiber and the x axis;
the related parameters of the aggregate comprise the centroid coordinate (X) of the aggregateAggregate material,YAggregate material) Radius R of aggregateAggregate materialAnd an aggregate interface thickness T;
the method for judging whether the invasion phenomenon is generated between the fiber and the aggregate is as follows:
a. calculating the coordinates (X) of the termination point of the randomly generated fibers according to the related parameters of the fibersend,Yend) The calculation formula is as follows: xend＝Xstart+Llength×cosθ；Yend＝Ystart+Llength×sinθ；
b. By randomly generating the coordinates (X) of the starting point of the fibresstart,Ystart) Coordinate of end point (X)end,Yend) And the coordinates of the centroid of the aggregate (X)Aggregate material,YAggregate material) Calculating the area S of a triangle formed by randomly generated fiber starting points, end points and an aggregate centroid;
c. by randomly generating the starting points, the end points and the aggregates of the fibresArea S of triangle formed by centroid and length L of fiberlengthCalculating the distance L from the centroid of the aggregate to the fiberdistance；
d. If the distance L from the centroid of the aggregate to the fiberdistanceGreater than or equal to the corresponding aggregate radius RAggregate materialAnd the sum of the thickness T of the interface layer, judging that the fiber and the aggregate do not have invasion, continuously judging whether the fiber and all the aggregates do not have invasion through the steps, and if so, ensuring that the fiber meets the requirement; otherwise, rejecting the fiber and regenerating the related parameters of the fiber and returning to the step a for re-calculation.
2. The method for building the fiber concrete mesoscopic model according to claim 1, wherein the calculation formula of the area S of the triangle is as follows:
Figure FDA0003615823630000021
Ldistancethe calculation formula of (c) is: l isdistance＝2S/Llength。
3. The method for building a fiber concrete mesoscopic model according to claim 1, wherein in step S2, the obtained relevant parameters of the fibers include:
on a microscopic level, an xoy plane coordinate system is established, and the related parameters of a certain fiber meeting the requirements in the model comprise the length L of the fiber1lengthThe coordinates of the starting point of the fiber (X1)start,Y1start) The rotation angle theta of the fiber1Angle of rotation theta1Is the angle between the fiber and the x-axis;
randomly generating relevant parameters of the ith fiber in the model, wherein the relevant parameters respectively comprise the length Li of the ith fiberlengthCoordinates of the starting point of the ith fiber (Xi)start,Yitart) The rotation angle theta of the ith fiberiRotation angle θ of ith fiberiIs the ith fiberAngle to the x-axis.
4. The method for building the fiber concrete mesoscopic model according to claim 3, wherein after the ith fiber is obtained through the calculation of step S1 and no invasion phenomenon is generated with all aggregates, whether the ith fiber and all existing fibers meeting the requirements in the model generate no invasion phenomenon is judged, and the judging method is as follows:
calculating the coordinates of the termination point of the ith fiber (Xi)end,Yiend) The calculation formula is as follows:
Xiend＝Xistart+Lilength×cosθi；
Yiend＝Yistart+Lilength×sinθi；
ii, calculating the area S of a triangle formed by the starting point of the ith fiber and the starting point and the ending point of a certain fiber in the modelstart；
Calculating the area S of a triangle formed by the midpoint of the ith fiber and the starting point and the ending point of a certain fiber existing in the modelcenter；
Iv, calculating the termination point of the ith fiber and the area S of a triangle formed by the starting point and the termination point of a certain fiber in the modelend；
V, respectively calculating the distance Ls from the starting point of the ith fiber to a certain fiber in the modeldistanceDistance Lc from the midpoint of the ith fiber to a fiber in the modeldistanceAnd the distance Le from the termination point of the ith fiber to a certain fiber existing in the modeldistance；
Vi, if 2Lcdistance＝Lsdistance+LedistanceThen calculate the line segment X1start XiendSlope of (1) and line segment Xiend X1endAnd judging whether the two slopes are equal or not, and then calculating a line segment X1start XistartSlope of (1) and line segment Xistart X1endAnd whether the two slopes are equal, as described above, if they are not equal, then it is determined to be randomThe generated ith fiber does not produce invasion phenomenon with a certain fiber existing in the model; according to the method, whether the ith fiber generated randomly does not produce invasion phenomenon with all fibers existing in the model or not is judged, if yes, the ith fiber meets the requirement of the model; if the condition is not met, the relevant parameters of the ith fiber are regenerated, and the step is returned to calculate again.
5. The method for establishing the fiber concrete mesoscopic model according to claim 4, wherein the area calculation formulas are respectively as follows:
Figure FDA0003615823630000031
Figure FDA0003615823630000032
Figure FDA0003615823630000033
according to area SstartArea ScenterSum area SendThe following can be determined:
Lsdistance＝2Sstart/L1length，Lcdistance＝2/ScenterL1length，Ledistance＝2Send/L1length。
CN202110040802.5A 2021-01-13 2021-01-13 Method for establishing fiber concrete mesoscopic model Active CN112700824B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110040802.5A CN112700824B (en) 2021-01-13 2021-01-13 Method for establishing fiber concrete mesoscopic model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110040802.5A CN112700824B (en) 2021-01-13 2021-01-13 Method for establishing fiber concrete mesoscopic model

Publications (2)

Publication Number Publication Date
CN112700824A CN112700824A (en) 2021-04-23
CN112700824B true CN112700824B (en) 2022-06-14

Family

ID=75514297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110040802.5A Active CN112700824B (en) 2021-01-13 2021-01-13 Method for establishing fiber concrete mesoscopic model

Country Status (1)

Country Link
CN (1) CN112700824B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006133150A1 (en) * 2005-06-08 2006-12-14 Mousseau Nicholas M Process of fabricating fiber reinforced concrete masonry units, and products made by said process
CN106874599A (en) * 2017-02-17 2017-06-20 武汉大学 The method for quickly generating the concrete random convex polyhedral aggregate model that cobble rubble is mingled with
CN108932393A (en) * 2018-07-13 2018-12-04 同济大学 A kind of concrete two dimension meso-mechanical model modeling method intrinsic based on aggregate geometry
CN110162849A (en) * 2019-05-07 2019-08-23 南京理工大学 A kind of modeling method of hybrid fiber concrete
CN110395963A (en) * 2019-08-16 2019-11-01 北京市高强混凝土有限责任公司 A kind of construction refuse regenerated concrete

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006133150A1 (en) * 2005-06-08 2006-12-14 Mousseau Nicholas M Process of fabricating fiber reinforced concrete masonry units, and products made by said process
CN106874599A (en) * 2017-02-17 2017-06-20 武汉大学 The method for quickly generating the concrete random convex polyhedral aggregate model that cobble rubble is mingled with
CN108932393A (en) * 2018-07-13 2018-12-04 同济大学 A kind of concrete two dimension meso-mechanical model modeling method intrinsic based on aggregate geometry
CN110162849A (en) * 2019-05-07 2019-08-23 南京理工大学 A kind of modeling method of hybrid fiber concrete
CN110395963A (en) * 2019-08-16 2019-11-01 北京市高强混凝土有限责任公司 A kind of construction refuse regenerated concrete

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Mesoscale modelling of dynamic tensile behaviour of fibre reinforced concrete with spiral fibres;Z Xu等;《Cement and Concrete Research》;20121231;第1457-1493页 *
基于细观尺度的钢纤维混凝土损伤破坏数值模拟研究;邵兵;《中国优秀博硕士寻味论文全文数据库（硕士）工程科技Ⅱ辑（月刊）》;20180315;全文 *

Also Published As

Publication number Publication date
CN112700824A (en) 2021-04-23

Similar Documents

Publication Publication Date Title
Turkan et al. Tracking of secondary and temporary objects in structural concrete work
JP2008186440A (en) Topology evolution optimization computing method for structural design
Soroushian et al. Morphological operations, planar mathematical formulations, and stereological interpretations for automated image analysis of concrete microstructure
CN111008423B (en) Optimal design method for double-curved-surface curtain wall
CN104820730B (en) A kind of method to be conflicted based on BIM solution regular reinforcements with prestress pipe
WO2021087842A1 (en) Method for measuring three-dimensional roughness on surface of concrete
CN112700824B (en) Method for establishing fiber concrete mesoscopic model
CN105809735A (en) Topology maintenance method based on three-dimensional geometric body combination
JP6201059B2 (en) Ground shape estimation program, ground shape estimation device, and ground shape estimation method
CN108829954A (en) A kind of reinforcing bar configuration method based on structural model profile features
CN112634455B (en) Method for repairing three-dimensional model ridge line by using cut triangular surface patches
CN105354878A (en) Electronic map data making method and apparatus
CN112712579B (en) Uniform level data parsing system
Zhao et al. A hybrid crack detection approach for scanning electron microscope image using deep learning method
Binjin et al. A novel construction quality control and management method based on BIM and 3D laser scanning technology
Dhondt Unstructured 20-node brick element meshing
Su et al. View frustum culling algorithm based on optimized scene management structure
CN116451305A (en) Coarse aggregate geometric modeling method based on real concrete bridge
CN103882867B (en) Against making the visual analysis method of permanent part column relative settlement under frame conditions
Chen et al. Quantitative surface quality evaluation for 3D-printed concrete with coarse aggregate through 3D scanning
Ren-xi et al. Generation method of normal vector from disordered point cloud
CN111400969B (en) Method for accelerating generation of unstructured right-angle grid
CN109544678B (en) Three-dimensional streamline parallel placement method with basically uniform distribution
Feng et al. Edge detection of concrete mesostructure based on DIS operator
CN115908542A (en) BIM-based method for rapidly calculating distance between polyhedrons

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant