This application is based upon and claims priority to Chinese patent application No. 201810690706.3 filed Jun. 28, 2018, which is incorporated herein by reference in its entirety.
The present disclosure relates to the technical field of fiber orientation prediction, and more particularly to a method and a device for predicting fiber orientation.
Fiber composite material is an important engineering material, which is widely used in engineering and has broad development prospect. The performance of fiber composite material is closely related to a distribution state of interior fiber orientation. Therefore, prediction of fiber orientation in the formation of fiber composite material is of great significance.
At present, finite element software is usually used in predicting fiber orientation, for determining directions of the strongest and weakest fiber orientations according to main strain directions in individual elements of a finite element model; and fiber orientation is predicted according to the magnitudes of a first, second and third main strain of each element. However, this method is directed to short fibers. Since the fibers in short fiber composite material are completely randomly blended or approximately randomly blended on a microstructure scale, a three-dimensionally woven short fiber reinforced composite material is an isotropic material.
However, a three-dimensionally woven long fiber composite material is a typical anisotropic material, and therefore, the above method is not suitable for predicting fiber orientation of the three-dimensionally woven long fiber composite material.
In order to overcome the problem in the related art, embodiments of the present disclosure provide a method and a device for predicting fiber orientation. The technical solutions are as follows:
According to a first aspect of one embodiment of the present disclosure, there is provided a method for predicting fiber orientation. The method includes:
obtaining a mechanical property parameter value of a fiber composite material;
obtaining a change value of a fiber property of the fiber composite material after deformation of the fiber composite material, according to the mechanical property parameter value; and
predicting fiber orientation in a forming process of the fiber composite material according to the change value of the fiber property and a preset fiber orientation model, the preset fiber orientation model indicating a mapping relationship between the change value of the fiber property and the fiber orientation after the deformation of the fiber composite material.
The technical solution provided in embodiment of the present disclosure may have the following beneficial effects: obtaining a mechanical property parameter value of a fiber composite material; obtaining a change value of a fiber property of the fiber composite material after deformation of the fiber composite material, according to the mechanical property parameter value; and predicting fiber orientation in a forming process of the fiber composite material according to the change value of the fiber property and a preset fiber orientation model, the preset fiber orientation model indicating a mapping relationship between the change value of the fiber property and the fiber orientation after the deformation of the fiber composite material. The present disclosure introduces a preset fiber orientation model for indicating a mapping relationship between the change value of the fiber property and the fiber orientation after deformation of the fiber composite material, so that fiber orientation in a forming process of the fiber composite material can be predicted from the fiber change value as obtained.
According to a second aspect of one embodiment of the present disclosure, there is provided a device for predicting fiber orientation. The device includes:
a processor, and
a memory configured to store instructions executable by the processor,
wherein the processor is configured to:
obtain a mechanical property parameter value of a fiber composite material;
obtain a change value of a fiber property of the fiber composite material after deformation of the fiber composite material, according to the mechanical property parameter value; and
predict fiber orientation in a forming process of the fiber composite material according to the change value of the fiber property and a preset fiber orientation model, the preset fiber orientation model indicating a mapping relationship between the change value of the fiber property and the fiber orientation after the deformation of the fiber composite material.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and, together with the description, serve to explain the principles of the present disclosure.
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. The following description refers to the accompanying drawings in which the same numbers in different drawings represent the same or similar elements unless otherwise represented. The implementations set forth in the following description of exemplary embodiments do not represent all implementations consistent with the disclosure. Instead, they are merely examples of apparatuses and methods consistent with aspects related to the disclosure as recited in the appended claims.
In order to solve the above technical problem, the present disclosure provides a method for predicting fiber orientation, by which fiber orientation of a three-dimensionally woven long fiber composite material can be predicted.
In step S101, a mechanical property parameter value of a fiber composite material is obtained.
The mechanical property parameter value of the fiber composite material can be obtained by using preset simulation software.
For example, the mechanical property parameter value of the composite material can be obtained by using the finite element software Digimat.
Of course, before obtaining the mechanical property parameter value of the fiber composite material by using the preset simulation software, it is also possible to determine whether the mechanical property parameter value of the fiber composite material can be obtained through experiments. If experiment conditions are not met, the mechanical property parameter value of the fiber composite material is obtained by using preset simulation software.
In step S102, a change value of a fiber property of the fiber composite material after deformation of the fiber composite material is obtained, according to the mechanical property parameter value.
In step S103, fiber orientation in a forming process of the fiber composite material is predicted according to the change value of the fiber property and a preset fiber orientation model. The preset fiber orientation model indicates a mapping relationship between the change value of fiber property and the fiber orientation after the deformation of the fiber composite material.
The fiber orientation refers to directions of fibers in a woven structure, which are generally perpendicular to each other after completion of weaving (as shown in
Referring to
It should be noted that the fiber composite material in the present disclosure includes a long fiber composite material; and the long fiber composite material includes a three-dimensionally woven long fiber composite material.
The technical solution provided in embodiment of the present disclosure may have the following beneficial effects: obtaining a mechanical property parameter value of a fiber composite material; obtaining a change value of a fiber property of the fiber composite material after deformation of the fiber composite material, according to the mechanical property parameter value; and predicting fiber orientation in a forming process of the fiber composite material according to the change value of the fiber property and a preset fiber orientation model, the preset fiber orientation model indicating a mapping relationship between the change value of the fiber property and the fiber orientation after the deformation of the fiber composite material. The present disclosure introduces a preset fiber orientation model for indicating a mapping relationship between the change value of the fiber property and the fiber orientation after deformation of the fiber composite material, so that fiber orientation in a forming process of the fiber composite material can be predicted from the fiber change value as obtained.
For example, the change value of the fiber property includes a change value of a fiber in a preset Cartesian coordinate system and a change value of a fiber in a relative coordinate system, the relative coordinate system being obtained according to a warp direction and a weft direction of the fiber composite material;
the change value of the fiber in the preset Cartesian coordinate system includes: a length change value of the fiber along respective coordinate axes in the preset Cartesian coordinate system, and a change value of an angle between a warp fiber and a weft fiber in the composite material in the preset Cartesian coordinate system; and
the change value of the fiber in the relative coordinate system includes: a length change value of the fiber along respective coordinate axes in the relative coordinate system, and a change value of an angle between a warp fiber and a weft fiber in the composite material in the relative coordinate system.
As shown in
In an embodiment, obtaining a change value of a fiber property of the fiber composite material after deformation of the fiber composite material includes the following substep:
obtaining a length change value of a fiber in the fiber composite material along respective coordinate axes in a preset Cartesian coordinate system and a change value of an angle between a warp fiber and a weft fiber in the fiber composite material in the preset Cartesian coordinate system, after the deformation of the fiber composite material.
The above substep may include the following substeps A1-A3:
A1: obtaining a first length value of the fiber along respective coordinate axes in the preset Cartesian coordinate system and a first angle value of an angle between the warp fiber and the weft fiber in the composite material in the preset Cartesian coordinate system, before deformation of the fiber,
A2: obtaining a second length value of the fiber along respective coordinate axes in the preset Cartesian coordinate system and a second angle value of an angle between the warp fiber and the weft fiber in the composite material in the preset Cartesian coordinate system, after the deformation of the fiber; and
A3: obtaining a length change value of the fiber along respective coordinate axes in the preset Cartesian coordinate system according to the second length value and the first length value, and obtaining a change value of the angle between the warp fiber and the weft fiber in the fiber composite material in the preset Cartesian coordinate system according to the second angle value and the first angle value.
After obtaining the first length value of the fiber along respective coordinate axes in the preset Cartesian coordinate system before deformation of the fiber and the second length value of the fiber along respective coordinate axes in the preset Cartesian coordinate system after the deformation of the fiber, the length change value of the fiber along respective coordinate axes in the preset Cartesian coordinate system can be obtained by subtracting the first length value from the second length value.
Similarly, after obtaining the first angle value of an angle between the warp fiber and the weft fiber in the composite material in the preset Cartesian coordinate system before deformation of the fiber and the second angle value of an angle between the warp fiber and the weft fiber in the composite material in the preset Cartesian coordinate system after the deformation of the fiber, the change value of the angle between the warp fiber and the weft fiber in the fiber composite material in the preset Cartesian coordinate system can be obtained by subtracting the first angle value from the second angle value.
The above steps can be implemented by simulation software, and the disclosure does not limit the type of the simulation software.
In an embodiment, obtaining a change value of a fiber property of the fiber composite material after the deformation of the fiber composite material includes the following substep:
obtaining a length change value of a fiber in the fiber composite material along respective coordinate axes in the relative coordinate system and a change value of an angle between a warp fiber and a weft fiber in the fiber composite material in the relative coordinate system, after the deformation of the fiber composite material.
The above substep may include the following steps B1-B3:
B1: obtaining a third length value of the fiber along respective coordinate axes in the relative coordinate system and a third angle value of an angle between the warp fiber and the weft fiber in the composite material in the relative coordinate system, before deformation of the fiber;
B2: obtaining a fourth length value of the fiber along respective coordinate axes in the relative coordinate system and a fourth angle value of an angle between the warp fiber and the weft fiber in the composite material in the relative coordinate system, after the deformation of the fiber, and
B3: obtaining a length change value of the fiber along respective coordinate axes in the relative coordinate system according to the fourth length value and the third length value, and obtaining a change value of the angle between the warp fiber and the weft fiber in the composite material in the relative coordinate system according to the fourth angle value and the third angle value.
After obtaining the third length value of the fiber along respective coordinate axes in the relative coordinate system before deformation of the fiber and the fourth length value of the fiber along respective coordinate axes in the relative coordinate system after the deformation of the fiber, the length change value of the fiber along respective coordinate axes in the relative coordinate system can be obtained by subtracting the third length value from the fourth length value.
Similarly, after obtaining the third angle value of an angle between the warp fiber and the weft fiber in the composite material in the relative coordinate system before deformation of the fiber and the fourth angle value of an angle between the warp fiber and the weft fiber in the composite material in the relative coordinate system after the deformation of the fiber, the change value of the angle between the warp fiber and the weft fiber of the fiber composite material in the relative coordinate system can be obtained by subtracting the third angle value from the fourth angle value.
The above steps can be implemented by simulation software, and the disclosure does not limit the type of the simulation software.
In an embodiment, prior to said predicting fiber orientation in the fiber composite material according to the change value of the fiber property and a preset fiber orientation model, the method further includes the following substeps C1-C3:
C1: obtaining a Cartesian component of stress-strain tensor according to a change value of the fiber in a preset Cartesian coordinate system. The change value of the fiber in the preset Cartesian coordinate system includes: a length change value of the fiber along respective coordinate axes in the predetermined Cartesian coordinate system, and a change value of an angle between a warp fiber and a weft fiber in the composite material in the predetermined Cartesian coordinate system.
For example, the Cartesian component of the stress-strain tensor may be expressed as:
where dε11 represents a length change value of the fiber along a first coordinate axis in the preset Cartesian coordinate system; dε22 represents a length change value of the fiber along a second coordinate axis in the preset Cartesian coordinate system; and dγ12 represents a change value of an angle between a warp fiber and a weft fiber in the composite material in the preset Cartesian coordinate system.
Taking
C2: obtaining a relative coordinate component of the stress-strain tensor according to a change value of the fiber in a relative coordinate system. The change value of the fiber in the relative coordinate system including: a length change value of the fiber along respective coordinate axes in the relative coordinate system, and a change value of an angle between a warp fiber and a weft fiber in the composite material in the relative coordinate system.
For example, the relative coordinate component of the stress-strain tensor can be expressed as:
where dε11 represents a length change value of the fiber along a first coordinate axis in the relative coordinate system; dε12 represents a length change value of the fiber along a second coordinate axis in the relative coordinate system; and dε
Taking
C3: obtaining the preset fiber orientation model according to a relationship among a Cartesian component of the stress-strain tensor, a relative coordinate component of the stress-strain tensor, and fiber orientation.
Since the three-dimensional woven long fiber composite material is a typical anisotropic material, non-uniformity and anisotropy of the material are the main difficulties. In this embodiment, in order to select a suitable preset fiber orientation model, a plane coordinate system (i.e., the relative coordinate system described above) that is consistent with the weft and warp of the woven fabric can be embedded in a shell element model. The relative coordinate components (i.e., an inverse stress component and a covariation strain component) of the stress-strain tensor in the constitutive relationship are introduced into the relative coordinate system described above. The preset fiber orientation model of the fiber composite material is derived by the transformation between the relative coordinate component of the stress-strain tensor and the Cartesian component of the stress-strain tensor.
Further, since stretching and shearing can be decoupled in the above constitutive relationship, the preset fiber orientation model can be obtained based on the above-described constitutive relationship.
In an embodiment, the relationship among the Cartesian component of the stress-strain tensor, the relative coordinate component of the stress-strain tensor and the fiber orientation includes:
the relative coordinate component of the stress-strain tensor is a product of the Cartesian component of the stress-strain tensor, and a preset trigonometric function value corresponding to an angular value of the fiber orientation.
For example, the relative coordinate component of the stress-strain tensor being a product of the Cartesian component of the stress-strain tensor, and a preset trigonometric function value corresponding to an angular value of the fiber orientation can be expressed as:
where the α represents an angle change value of the angle between the fiber and the first coordinate axis in the preset Cartesian coordinate system; and the θ represents an angle value of the fiber orientation.
Taking
In practical application, the finite element software Abaqus can be used to predict fiber orientation in a forming process of fiber composite material.
Specifically, using Abaqus' VUMAT user-defined subroutine function, a user can input the above fiber orientation model via VUMAT.
Since the Abaqus software is able to calculate the change value of the fiber property of the fiber composite material (as shown in
In this case, the distribution result of the fiber orientation is output as shown in
The following is a device embodiment of the present disclosure, which may be used to implement the method embodiments of the present disclosure.
a first obtaining module 11, configured to a mechanical property parameter value of a fiber composite material;
a second obtaining module 12, configured to obtain a change value of a fiber property of the fiber composite material after deformation of the fiber composite material, according to the mechanical property parameter value obtained by the first obtaining module 11; and
a predicting module 13, configured to predict fiber orientation in a forming process of the fiber composite material according to the change value of the fiber property obtained by the second obtaining module 12 and a preset fiber orientation model, the preset fiber orientation model indicating a mapping relationship between the change value of the fiber property and the fiber orientation after the deformation of the fiber composite material.
In an embodiment, as shown in
the first obtaining sub-module 121 is configured to obtain a length change value of a fiber in the fiber composite material along respective coordinate axes in a predetermined Cartesian coordinate system and a change value of an angle between a warp fiber and a weft fiber in the fiber composite material in the predetermined Cartesian coordinate system, after the deformation of the fiber composite material; and
the second obtaining sub-module 122 is configured to obtain a length change value of a fiber in the fiber composite material along respective coordinate axes in a relative coordinate system and a change value of an angle between a warp fiber and a weft fiber in the fiber composite material in the relative coordinate system after the deformation of the fiber composite material.
In an embodiment, as shown in
the third obtaining module 14 is configured to obtain a Cartesian component of stress-strain tensor according to a change value of the fiber in a preset Cartesian coordinate system, the change value of the fiber in the preset Cartesian coordinate system including: a length change value of the fiber along respective coordinate axes in the predetermined Cartesian coordinate system, and a change value of an angle between a warp fiber and a weft fiber in the composite material in the predetermined Cartesian coordinate system;
the fourth obtaining module 15 is configured to obtain a relative coordinate component of the stress-strain tensor according to a change value of the fiber in a relative coordinate system, the change value of the fiber in the relative coordinate system including: a length change value of the fiber along respective coordinate axes in the relative coordinate system, and a change value of an angle between a warp fiber and a weft fiber in the composite material in the relative coordinate system; and
the fifth obtaining module 16 is configured to obtain the preset fiber orientation model according to a relationship among a Cartesian component of the stress-strain tensor obtained by the third obtaining module 14, a relative coordinate component of the stress-strain tensor obtained by the fourth obtaining module 15, and fiber orientation.
In an embodiment, the Cartesian component of the stress-strain tensor is expressed as:
where dε11 represents a length change value of the fiber along a first coordinate axis in the preset Cartesian coordinate system; dε22 represents a length change value of the fiber along a second coordinate axis in the preset Cartesian coordinate system; and dγ12 represents a change value of an angle between a warp fiber and a weft fiber in the composite material in the preset Cartesian coordinate system.
In an embodiment, the relative coordinate component of the stress-strain tensor is expressed as:
where dε11 represents a length change value of the fiber along a first coordinate axis in the relative coordinate system; dε22 represents a length change value of the fiber along a second coordinate axis in the relative coordinate system; and dγ12 represents a change value of an angle between a warp fiber and a weft fiber in the relative coordinate system.
In an embodiment, the relationship among the Cartesian component of the stress-strain tensor, the relative coordinate component of the stress-strain tensor and the fiber orientation includes:
the relative coordinate component of the stress-strain tensor is a product of the Cartesian component of the stress-strain tensor, and a preset trigonometric function value corresponding to an angular value of the fiber orientation.
In an embodiment, said the relative coordinate component of the stress-strain tensor being a product of the Cartesian component of the stress-strain tensor, and a preset trigonometric function value corresponding to an angular value of the fiber orientation includes:
where the α represents an angle change value of the angle between the fiber and the first coordinate axis in the preset Cartesian coordinate system; the θ represents an angle value of the fiber orientation.
In an embodiment, as shown in
the third obtaining sub-module 111 is configured to obtain the mechanical property parameter value of the fiber composite material by using preset simulation software.
For the device in the above embodiment, the specific manners in which various modules operate have been described in detail in the embodiment relating to the above method, and therefore will not be repeated herein.
Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed here. This application is intended to cover any variations, uses, or adaptations of the disclosure following the general principles thereof and including such departures from the present disclosure as come within known or customary practice in the art. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the disclosure being indicated by the following claims.
It will be appreciated that the present disclosure is not limited to the exact construction that has been described above and illustrated in the accompanying drawings, and that various modifications and changes may be made without departing from the scope thereof. It is intended that the scope of the disclosure only be limited by the appended claims.
Number | Date | Country | Kind |
---|---|---|---|
201810690706.3 | Jun 2018 | CN | national |