CN111563894B - Method and system for measuring bending stiffness of continuous fiber reinforced material - Google Patents
Method and system for measuring bending stiffness of continuous fiber reinforced material Download PDFInfo
- Publication number
- CN111563894B CN111563894B CN202010396164.6A CN202010396164A CN111563894B CN 111563894 B CN111563894 B CN 111563894B CN 202010396164 A CN202010396164 A CN 202010396164A CN 111563894 B CN111563894 B CN 111563894B
- Authority
- CN
- China
- Prior art keywords
- bending
- curve
- image
- bending deformation
- shape
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/64—Analysis of geometric attributes of convexity or concavity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/68—Analysis of geometric attributes of symmetry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30172—Centreline of tubular or elongated structure
Abstract
The invention relates to a method and a system for measuring bending stiffness of a continuous fiber reinforced material. The method comprises the following steps: acquiring a bending deformation image of a continuous fiber reinforced material sample; extracting a central line of a bending deformation shape by adopting an image processing method according to the bending deformation image; fitting the central line of the bending deformation shape by adopting a uniform quartic B-spline curve to obtain a fairing fitting curve; determining curvature and bending moment according to the fairing fitting curve; and determining the bending rigidity of the continuous fiber reinforced material sample according to the curvature and the bending moment. The invention can accurately measure the bending rigidity of the continuous fiber reinforced material.
Description
Technical Field
The invention relates to the field of bending stiffness measurement, in particular to a method and a system for measuring bending stiffness of a continuous fiber reinforced material.
Background
Composite materials have found widespread use in the aerospace and automotive industries in recent years due to their superior properties of low density, high stiffness and high strength. The composite material consists of a matrix and a reinforcement, wherein the rigidity of the reinforcement is far greater than that of the matrix, and the rigidity and the strength are provided for the composite material.
In structural engineering applications, continuous fiber materials have excellent mechanical properties and good formability and are often used as reinforcement materials for composite materials. Resin Transfer Molding (RTM) is currently the predominant form-making method for continuous fiber resin-based composites. In RTM, the first stage is the formation of dry fiber reinforcement to make a preform, and after formation the orientation and distribution of the fibers in the preform has a significant effect on the resin infusion in the second stage (mainly affecting resin permeability), and the orientation of the reinforcing fibers plays a crucial role in the mechanical properties of the final material. Therefore, numerical simulation of the forming process is essential, which can accurately predict the deformation of the fiber reinforcement, optimize the forming process conditions, avoid defects, and improve efficiency. For accurate modeling, an accurate knowledge of the mechanical behavior of the fiber-reinforced material is required. The fiber reinforcement material undergoes primarily tensile, in-plane shear and bending deformation during forming. Currently, most of the traditional simulation models only consider the stretching and in-plane shear deformation of fiber reinforced body materials, and adopt the film deformation theory to neglect the bending rigidity of the materials. One of the main reasons for this problem is the lack of an accurate measurement of the flexural stiffness of the material that characterizes the fiber reinforcement.
Disclosure of Invention
The invention aims to provide a method and a system for measuring the bending stiffness of a continuous fiber reinforced material, which can accurately measure the bending stiffness of the continuous fiber reinforced material.
In order to achieve the purpose, the invention provides the following scheme:
a method of measuring bending stiffness of a continuous fiber reinforced material, comprising:
acquiring a bending deformation image of a continuous fiber reinforced material sample;
extracting a central line of a bending deformation shape by adopting an image processing method according to the bending deformation image;
fitting the central line of the bending deformation shape by adopting a uniform quadruplicate B-spline curve to obtain a fairing fitting curve;
determining curvature and bending moment according to the fairing fitting curve;
and determining the bending rigidity of the continuous fiber reinforced material sample according to the curvature and the bending moment.
Optionally, the acquiring a bending deformation image of the continuous fiber reinforced material sample specifically includes:
and acquiring a bending deformation image of the continuous fiber reinforced material sample by using a CCD (charge coupled device) camera.
Optionally, the extracting a central line of the bending deformation shape by using an image processing method according to the bending deformation image specifically includes:
smoothing, contrast enhancement, thresholding and shape smoothing are carried out on the curved image to obtain a processed image;
and performing skeletonization on the processed image by adopting a thinning method, and extracting a central line of the bending deformation shape.
Optionally, fitting the central line of the bending deformation shape by using a uniform quartic B-spline curve to obtain a fairing fitting curve, specifically including:
fitting the central line of the bending deformation shape by adopting a piecewise curve function to obtain a first fairing fitting curve;
and optimizing the first fairing fit curve by adopting an energy functional minimization and least square error combined method to obtain a second fairing fit curve.
A continuous fiber reinforced material bending stiffness measurement system comprising:
the image acquisition module is used for acquiring a bending deformation image of the continuous fiber reinforced material sample;
the image processing module is used for extracting a central line of the bending deformation shape according to the bending deformation image by adopting an image processing method;
the curve fitting module is used for fitting the central line of the bending deformation shape by adopting a uniform quartic B-spline curve to obtain a smooth fitting curve;
the curvature/bending moment determining module is used for determining curvature and bending moment according to the fairing fitting curve;
and the bending rigidity determining module is used for determining the bending rigidity of the continuous fiber reinforced material sample according to the curvature and the bending moment.
Optionally, the image obtaining module specifically includes:
and the image acquisition unit is used for acquiring the bending deformation image of the continuous fiber reinforced material sample by adopting a CCD (charge coupled device) camera.
Optionally, the image processing module specifically includes:
the image processing unit is used for carrying out smoothing, contrast enhancement, thresholding and shape smoothing on the curved image to obtain a processed image;
and the central line extraction unit is used for performing skeletonization on the processed image by adopting a thinning method and extracting the central line of the bending deformation shape.
Optionally, the curve fitting module specifically includes:
the fitting unit is used for fitting the central line of the bending deformation shape by adopting a piecewise curve function to obtain a first fairing fitting curve;
and the optimization unit is used for optimizing the first fairing fitting curve by adopting an energy functional minimization and least square error combined method to obtain a second fairing fitting curve.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
compared with the existing bending test method, the bending test method is simple to operate and set, allows the measurement of large enough bending deformation, obtains a wide bending moment-curvature relation, can accurately measure the bending rigidity of the continuous fiber reinforced material, and can be applied to the measurement of the bending rigidity of other materials. Based on the quartic uniform B-spline curve, selecting proper function segmentation quantity, overcoming the defect that a single polynomial curve function is limited by fitting precision and is easy to cause oscillation, and obtaining a smoother curve; in addition, the energy functional minimization and the least square error minimization are combined, the curve fitting error and the curve energy function are considered, the defect that unnecessary swinging is easily caused by a least square curve fitting method is overcome, a more reasonable curve is obtained, and the bending rigidity measurement accuracy is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a schematic diagram of a bending stiffness measurement method of the present invention;
FIG. 2 is a schematic view of a bending test apparatus according to the present invention;
FIG. 3 is a flow chart of a method for measuring bending stiffness of a continuous fiber reinforced material according to the present invention;
FIG. 4 is a schematic view of the curved shape smoothing process of the present invention;
FIG. 5 is a skeletonization and midline plot of the present invention;
FIG. 6 is a graph of the effect of the number of uniform quadruplicate B spline line segments on curvature map fitting error and smoothness;
FIG. 7 is an experimental plot of the present invention;
FIG. 8 is a structural diagram of a system for measuring bending stiffness of a continuous fiber reinforced material according to 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 drawings in 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 derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a method and a system for measuring the bending rigidity of a continuous fiber reinforced material, which can accurately measure the bending rigidity of the continuous fiber reinforced material.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Recent studies have shown that bending stiffness is closely related to the shape and morphology of wrinkles, and therefore, to accurately predict defects that may occur during the forming process of wrinkles and the like, bending deformation must be considered in numerical simulation. Since relative slippage between fibers is highly likely to occur, the bending stiffness of the continuous fiber reinforcement is not directly related to the in-plane tensile modulus, much less than that of a classical continuous material structure, and therefore the values need to be measured experimentally. Therefore, the invention provides a method for measuring the bending stiffness of a continuous fiber reinforced material, which adopts a cantilever beam clamping mode to carry out a bending test, utilizes a CCD (charge coupled device) camera to shoot a bending deformation shape image, calculates the bending deformation shape curvature of the fiber reinforced material based on an improved uniform B-spline curve theory, and further calculates the bending stiffness of the material; the method for measuring the bending stiffness of the continuous fiber reinforced material can also be applied to other material bending stiffness test methods.
FIG. 1 is a schematic diagram of a bending stiffness measurement method of the present invention. Firstly, shooting a sample bending deformation shape image by a high-resolution CCD camera (2048 multiplied by 2048 pixels), then extracting a central line of the bending deformation shape by adopting an image processing technology, finally fitting a bending curve function based on an optimized uniform B spline curve, calculating curvature and corresponding bending moment along the central line by using the fitting curve, wherein the slope of the bending moment-curvature curve is the bending rigidity of the sample.
FIG. 2 is a schematic view of a bending test apparatus according to the present invention. The cross section of the sample was left white with some dry white powder (to increase contrast with the background (black) and facilitate taking the midline in later image processing), the sample was clamped and deformed in cantilever fashion using a metal holder, a CCD camera (for taking a picture of the sample deformation) was connected to the acquisition software, and the image acquired by the acquisition software was transferred to the image processing software ImageJ for image processing. In image processing, the scale is placed on the same plane as the sample front edge and the scale is used to calculate the scaling factor (pixels/mm).
FIG. 3 is a flow chart of a method for measuring bending stiffness of a continuous fiber reinforced material according to the present invention. As shown in fig. 3, a method for measuring bending stiffness of a continuous fiber reinforced material includes:
step 101: the method for acquiring the bending deformation image of the continuous fiber reinforced material sample specifically comprises the following steps:
and acquiring a bending deformation image of the continuous fiber reinforced material sample by using a CCD (charge coupled device) camera.
A high-resolution CCD camera (2048 x 2048 pixels) is used for shooting a static image of the bending deflection shape of the sample, and the shot image is transmitted to acquisition software for processing.
Step 102: extracting a central line of the bending deformation shape by adopting an image processing method according to the bending deformation image, and specifically comprising the following steps of:
and carrying out smoothing, contrast enhancement, thresholding and shape smoothing on the bending deformation image to obtain a processed image. For the above four processes, the specific operations are as follows: 1) Selecting a Gaussian blur filter to carry out smoothing processing on the bending deformation image, and reducing the noise level in the image; 2) Contrast enhancement optimizes signal-to-noise ratio and normalizes the signal for processing and visualization; 3) Performing thresholding to separate pixels having luminance values within a desired range from pixels not within the desired range (also referred to as segmentation) and to separate the bending deflection shape of the sample from other parts; 4) The shape smoothing process is performed to make the curved shape smoother. FIG. 4 is a schematic diagram of the smoothing process for curved shapes according to the present invention, wherein a is the shape of noise and b is the shape of noise after smoothing.
And performing skeletonization on the processed image by adopting a thinning method, and extracting a central line of the bending deformation shape. Extracting a central line of the bending deflection shape of the test piece by using a skeletonization method defined by a thinning method in ImageJ, and specifically operating the following steps: using the scaling factor, physical cartesian coordinates of points on the skeleton are generated, and the curve represented by these discrete points is the centerline (deflection curve) of the curved deformed shape, with the results shown in fig. 5, which is the skeletonized and centerline plot of the present invention.
Step 103: adopting uniform four-time B-spline curve to fit the central line of the bending deformation shape to obtain a fairing fitting curve, which specifically comprises the following steps:
and 1031, fitting the central line of the bending deformation shape by adopting a piecewise curve function to obtain a first fairing fitting curve.
The monopolynomial curve function usually requires a high order to have a good fitting accuracy and is prone to give rise to oscillations, whereas through experimental observations the bending deflection curve of a continuous fiber reinforcement is usually composed of a plurality of curve shapes, so that a piecewise curve function is used. For a given curve, its curvature is plotted against the arc length, resulting in a curvature map of the curve. The curve has good smoothness if there are no or few slope discontinuities in the curvature map. Selection of the number of segments for a uniform quartic B-spline: too few segments are used to represent the characteristics of the curve represented by these sample points, and too many segments are used to make the fitted curve follow the noise and have many unwanted fluctuations. Increasing the number of segments can reduce the fitting error between the sample data point and the fitting curve, but the smoothness of the curvature graph on the fitting curve becomes poor, and both the fitting precision and the curvature smoothness should be considered, so as to achieve a better fitting effect solution. In the invention, curve fitting is respectively carried out for the number of segments of one, three and four, so as to respectively obtain the results shown in FIG. 6, and in the process, when the curvature diagram begins to fluctuate, although the tolerance requirement cannot be met, the increase of the number of segments is stopped, and finally, the segments are selected to be divided into three segments. FIG. 6 is a graph showing the effect of the number of uniform quadruplicate B spline segments on curvature map fitting error and smoothness, according to the present invention, wherein FIGS. 6 (a), (B) and (c) show the results obtained when one, three and four segments of fitting are used, respectively (left of the graph: sample data points and fitting curves; right of the graph: corresponding curvature distribution diagram): (a) a fitting error of 5.71 when a piece of fitting is adopted; (b) fitting error is 3.31 when three-stage fitting is adopted; and (c) adopting a four-segment fitting, wherein the fitting error is 3.1).
Step 1032: and optimizing the first fairing fitted curve by adopting a method combining energy functional minimization and least square error to obtain a second fairing fitted curve.
The energy functional minimization is combined with the least square error method, the defect of the least square curve fitting method is overcome, and the energy functional is used as the fairing term, so that the curve and the curved surface are smoother. And a fitting target function form is adopted, and curve fitting errors and a curve energy function are considered.
Step 104: and determining the curvature and the bending moment according to the fairing fitted curve.
The calculation process of the curvature and bending moment in the invention is as follows:
the uniform quartic B-spline curve is a piecewise polynomial curve composed of multiple segments, each segment is a quartic B-spline function, and m segments are respectively numbered as 1,2 i-1 And D i Are the start and end of the ith segment. The position vector of the point p on the ith segment is:
wherein X and Y are global Cartesian coordinates, u (i) Is the local parameter value of the point p along the i-th segment, B, varying from 0 to 1 i+k-1 Is a control point, M k,4 Is the basis function of a quartic B-spline curve:
from the local function of the p-point in equation (1), a global functional form is derived that represents the p-point position vector:
where m is the number of bus segments of the uniform quartic B-spline and t is the global parameter value for point p that varies from 0 to 1 along the uniform quartic B-spline.
Local parameter value u (i) The relationship to the global parameter value t is as follows:
u (i) ＝mt-[mt] (4)
[ mt ] represents the maximum integer not exceeding the product of m and t, i is the line segment number where the point p is located, and can be obtained according to the global parameter value:
i＝[mt]+1 (5)
N h,4 (t) global basis functions called uniform quartic B-spline curves:
combining the minimization of the energy functional with the minimization of the least square error, and considering the curve fitting error and the curve energy function to obtain the target function as follows:
where n is the total number of sample points, w j α, β, τ are weighting factorsA seed; c (t) is a position vector of the global parameter value t on the uniform quartic B-spline curve, C ' (t), C ' (t) and C ' (t) represent the first, second and third derivatives of C (t) with respect to the global parameter t, respectively;
is the sample point P j Wherein chord length methods are used to calculate their values;
is a vector that stores C (t) control points; p = [ P = 0 P 1 … P n ] T Is a vector that stores sample points; a is a scalar (n + 1) × (m + 4) matrix, and W is a memory weighting factor W j An (n + 1) × (n + 1) matrix of:
n, N 'and N' are the first, second and third derivatives of N, respectively, with respect to the global parameter t; n is a vector containing the global basis functions of the uniform quartic B-spline, N = [ N = 0，4 (t)···N m+3，4 (t)] T (ii) a Due to the similarity of the mechanical properties with the material, will K =: [ integral ] n t (αN′(N′) T +βN″(N″) T +τN″′(N″′) T ) dt is called the system stiffness matrix; item(s)
Referred to as the error term of the fit,
called the fairing term, which acts to smoothly fit the curve.
To implement the objective function solution, some boundary constraints are introduced. The fitted curve is required to pass through the first sampling point (the first sampling point is the pinch point of the bending deformation curve), the second derivative of the fitted curve of the last sampling point is zero (the last sampling point is the free end of the bending deformation curve), and the corresponding mathematical constraints are as follows:
the objective function under the boundary conditions is:
r matrix satisfies
D＝[C(0) C″(1)] T Ψ is a Lagrangian multiplier vector, = [ ψ = 1 ψ 2 ] T The other terms have the same meaning, and the minimum is obtained
The following equation is obtained:
solving equation (15) yields equation (16)
Calculating the curvature of any point C (t) on the uniform quartic B-spline curve by using a formula (16);
calculating the bending moment of any point A on the flexibility curve under the action of gravity through a formula (17);
wherein s is the abscissa of the curve of the point A; q is the weight per unit length of the sample (N/mm);
a unit vector that is the direction of gravity; l is the total length of the deflection curve; u is the Frey's inner coordinate of point B moving from A to F along the curve; du is the differential of u.
Step 105: and determining the bending rigidity of the continuous fiber reinforced material sample according to the curvature and the bending moment.
And calculating a bending moment-curvature curve according to the curvature and the bending moment, wherein the bending moment-curvature curve is a bending rigidity curve of the material. ( Note: the bending moment in the bending moment-curvature curve refers to the bending moment of the material per unit width, which is equal to M/d, and d is the width of the material sample in the bending test. )
Continuous fiber reinforced composite materials such as Hexcel G1151 and 2.5D interlayer angle interlocking Interlock are selected to verify the measurement method, and experimental results prove that the bending stiffness measurement method is simple and convenient to use and effective in calculation results, and the calculation results are shown in figure 7. Fig. 7 is a graph of experimental bending according to the present invention, where a is a curvature graph using Hexcel G1151, b is a curvature-bending moment graph using Hexcel G1151, c is a curvature graph using 2.5D inter-layer angle Interlock, and D is a curvature-bending moment graph using 2.5D inter-layer angle Interlock, and it should be noted that bending may occur during bending, but when bending is small, the influence on the test result may be ignored.
The invention provides a method for measuring bending stiffness of a continuous fiber reinforced material, which is used for the research of determining the bending stiffness of the continuous fiber reinforced material by an experimental method; the method has the main technical idea that the geometric characteristics of the bending deformation of the test piece are obtained through an image processing technology, a fairing bending deformation central line curve is obtained through fitting based on an improved uniform quartic B-spline curve fitting method, the bending deformation curvature and the bending moment of the continuous fiber reinforced material can be accurately measured and calculated through combination of experimental calculation, and the bending rigidity of the fiber reinforced material is further obtained.
FIG. 8 is a structural diagram of a system for measuring bending stiffness of a continuous fiber reinforced material according to the present invention. As shown in fig. 8, a continuous fiber reinforced material bending stiffness measurement system includes:
the image acquisition module 201 is used for acquiring a bending deformation image of the continuous fiber reinforced material sample.
And the image processing module 202 is configured to extract a central line of the curved deformation shape according to the curved deformation image by using an image processing method.
And the curve fitting module 203 is used for fitting the central line of the bending deformation shape by adopting a uniform quartic B-spline curve to obtain a smooth fitting curve.
And a curvature/bending moment determining module 204, configured to determine a curvature and a bending moment according to the fairing fit curve.
And the bending rigidity determining module 205 is used for determining the bending rigidity of the continuous fiber reinforced material sample according to the curvature and the bending moment.
The image obtaining module 201 specifically includes:
and the image acquisition unit is used for acquiring the bending deformation image of the continuous fiber reinforced material sample by adopting a CCD camera.
The image processing module 202 specifically includes:
and the image processing unit is used for carrying out smoothing, contrast enhancement, thresholding and shape smoothing on the curved image to obtain a processed image.
And the central line extraction unit is used for performing skeletonization on the processed image by adopting a thinning method and extracting the central line of the bending deformation shape.
The curve fitting module 203 specifically includes:
and the fitting unit is used for fitting the central line of the bending deformation shape by adopting a piecewise curve function to obtain a first fairing fitting curve.
And the optimization unit is used for optimizing the first fairing fit curve by adopting a method combining energy functional minimization and least square error to obtain a second fairing fit curve.
In the present specification, the embodiments are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principle and the embodiment of the present invention are explained by applying specific examples, and the above description of the embodiments is only used to help understanding the method and the core idea of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (4)
1. A method for measuring bending stiffness of a continuous fiber reinforced material is characterized by comprising the following steps:
acquiring a bending deformation image of a continuous fiber reinforced material sample;
extracting a central line of the bending deformation shape by adopting an image processing method according to the bending deformation image;
the extracting the central line of the bending deformation shape by adopting an image processing method according to the bending deformation image specifically comprises the following steps:
smoothing, contrast enhancement, thresholding and shape smoothing are carried out on the bending deformation image to obtain a processed image;
performing skeletonization on the processed image by adopting a thinning method, and extracting a central line of a bending deformation shape;
the smoothing, contrast enhancement, thresholding and shape smoothing processing are performed on the curved image to obtain a processed image, and specifically include:
selecting a Gaussian blur filter to carry out smoothing processing on the bending deformation image, and reducing the noise level in the image; contrast enhancement optimizes signal-to-noise ratio and normalizes the signal for processing and visualization; performing thresholding processing to separate pixels with brightness values within a desired range from pixels without brightness values within the desired range, so as to separate the bending deflection shape of the sample from other parts; performing shape smoothing processing to make the curved shape smoother;
the skeletonization is executed on the processed image by adopting a thinning method, and the central line of the bending deformation shape is extracted, which specifically comprises the following steps:
extracting a central line of the bending deflection shape of the test piece by using a skeletonization method defined by a thinning method in ImageJ, and specifically operating the following steps: generating physical Cartesian coordinates of points located on the skeleton using the scaling factor, the curve represented by the discrete points being a centerline of the curved deformed shape;
fitting the central line of the bending deformation shape by adopting a uniform quadruplicate B-spline curve to obtain a fairing fitting curve;
the step of fitting the central line of the bending deformation shape by adopting the uniform quartic B-spline curve to obtain the fairing fitting curve specifically comprises the following steps:
fitting the central line of the bending deformation shape by adopting a piecewise curve function to obtain a first fairing fitting curve;
optimizing the first fairing fitted curve by adopting a method combining energy functional minimization and least square error to obtain a second fairing fitted curve;
determining curvature and bending moment according to the fairing fitting curve;
and determining the bending rigidity of the continuous fiber reinforced material sample according to the curvature and the bending moment.
2. The method for measuring the bending stiffness of the continuous fiber reinforced plastic material according to claim 1, wherein the step of acquiring the bending deformation image of the continuous fiber reinforced plastic material sample specifically comprises the following steps:
and collecting the bending deformation image of the continuous fiber reinforced material sample by adopting a CCD camera.
3. A continuous fiber reinforced material bending stiffness measurement system, comprising:
the image acquisition module is used for acquiring a bending deformation image of the continuous fiber reinforced material sample;
the image processing module is used for extracting a central line of the bending deformation shape according to the bending deformation image by adopting an image processing method;
the image processing module specifically comprises:
the image processing unit is used for carrying out smoothing, contrast enhancement, thresholding and shape smoothing on the bending deformation image to obtain a processed image;
a central line extraction unit, configured to perform skeletonization on the processed image by using a thinning method, and extract a central line of a bending deformation shape;
the image processing unit specifically includes:
selecting a Gaussian blur filter to carry out smoothing processing on the bending deformation image, and reducing the noise level in the image; contrast enhancement optimizes signal-to-noise ratio and normalizes the signal for processing and visualization; performing thresholding processing to separate pixels with brightness values within a desired range from pixels without brightness values within the desired range, so as to separate the bending deflection shape of the sample from other parts; performing shape smoothing processing to make the curved shape smoother;
the centerline extraction unit specifically includes:
extracting a central line of the bending deflection shape of the test piece by using a skeletonization method defined by a thinning method in ImageJ, and specifically operating the following steps: generating physical Cartesian coordinates of points located on the skeleton using the scaling factor, the curve represented by the discrete points being a centerline of the curved deformed shape;
the curve fitting module is used for fitting the central line of the bending deformation shape by adopting a uniform quartic B-spline curve to obtain a fairing fitting curve;
the curve fitting module specifically comprises:
the fitting unit is used for fitting the central line of the bending deformation shape by adopting a piecewise curve function to obtain a first fairing fitting curve;
the optimization unit is used for optimizing the first fairing fit curve by adopting an energy functional minimization and least square error combination method to obtain a second fairing fit curve;
the curvature/bending moment determining module is used for determining curvature and bending moment according to the fairing fitting curve;
and the bending rigidity determining module is used for determining the bending rigidity of the continuous fiber reinforced material sample according to the curvature and the bending moment.
4. The system for measuring bending stiffness of continuous fiber reinforced material according to claim 3, wherein the image acquisition module specifically comprises:
and the image acquisition unit is used for acquiring the bending deformation image of the continuous fiber reinforced material sample by adopting a CCD (charge coupled device) camera.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010396164.6A CN111563894B (en) | 2020-05-12 | 2020-05-12 | Method and system for measuring bending stiffness of continuous fiber reinforced material |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010396164.6A CN111563894B (en) | 2020-05-12 | 2020-05-12 | Method and system for measuring bending stiffness of continuous fiber reinforced material |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111563894A CN111563894A (en) | 2020-08-21 |
| CN111563894B true CN111563894B (en) | 2022-11-29 |
Family
ID=72074629
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010396164.6A Active CN111563894B (en) | 2020-05-12 | 2020-05-12 | Method and system for measuring bending stiffness of continuous fiber reinforced material |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111563894B (en) |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102914473A (en) * | 2012-09-27 | 2013-02-06 | 大连民族学院 | Method for recognizing cross-section bending moment and curvature relation of reinforced concrete beam |
| CN103605171A (en) * | 2013-09-10 | 2014-02-26 | 国家电网公司 | All-sky imaging instrument and cloud layer characteristic analysis method based on all-sky imaging instrument |
| CN104597056B (en) * | 2015-02-06 | 2017-04-19 | 北京中科纳新印刷技术有限公司 | Method for detecting ink-jet printing ink dot positioning accuracy |
| US9852348B2 (en) * | 2015-04-17 | 2017-12-26 | Google Llc | Document scanner |
| CN108917591A (en) * | 2018-04-27 | 2018-11-30 | 中国铁道科学研究院集团有限公司基础设施检测研究所 | Rail profile autoegistration method and device under a kind of dynamic environment |
| CN109204885B (en) * | 2018-11-07 | 2024-01-30 | 哈尔滨工业大学 | Shape memory unfolding mechanism |
-
2020
- 2020-05-12 CN CN202010396164.6A patent/CN111563894B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111563894A (en) | 2020-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104697462B (en) | A kind of aerial blade surface-type feature parameter extracting method based on axis | |
| CN106407502B (en) | Blade section profile parameter evaluation method based on optimal matching | |
| CN106372283B (en) | A kind of thin wall obtained towards digital photography surveys three-dimensional appearance Processing Method of Point-clouds | |
| CN109118574A (en) | A kind of fast reverse modeling method extracted based on three-dimensional feature | |
| Wang et al. | Digital image correlation (DIC) based damage detection for CFRP laminates by using machine learning based image semantic segmentation | |
| Liang et al. | Curvature determination in the bending test of continuous fibre reinforcements | |
| Wang et al. | Improvement of springback prediction accuracy for Z-section profiles in four-roll bending process considering neutral layer shift | |
| CN102768181B (en) | Method for quantitatively evaluating orientation degree of short fiber reinforced composite fibers | |
| CN109583377B (en) | Control method and device for pipeline model reconstruction and upper computer | |
| CN108614917B (en) | Consider the bending pipes springback Prediction method of elastic moding and section elliptical distortion | |
| CN111563894B (en) | Method and system for measuring bending stiffness of continuous fiber reinforced material | |
| CN113516695A (en) | Point cloud registration strategy in laser profilometer flatness measurement | |
| Liu et al. | Optimization of post-stretching elongation in stretch bending of aluminum hollow profile | |
| CN112082496A (en) | Concrete internal deformation measurement method and system based on improved digital volume image correlation method | |
| CN106021779A (en) | Method for calculating bending rigidity of main box section of high-aspect-ratio wing | |
| Ruiz-Gironés et al. | Defining an 2-disparity Measure to Check and Improve the Geometric Accuracy of Non-interpolating Curved High-order Meshes | |
| CN108536932B (en) | Aviation blade stacking axis verticality calculation method based on mutual torsion constraint condition | |
| CN113468673B (en) | Section optimization method for shape of associated joint | |
| CN111444578B (en) | Automatic calibration method of variable modulus model parameters based on bending process | |
| Yao et al. | Robust locally weighted regression for profile measurement of magnesium alloy tube in hot bending process | |
| CN115186376A (en) | Small-aspect-ratio flying wing aircraft rolling torque model and parameter identification method thereof | |
| CN107102056A (en) | The flexible construction method of the unit of defect and magnetic leakage signal | |
| Tricoche et al. | A fast and robust method for visualizing separation line features | |
| Jiang et al. | Urine sediment image segmentation based on level set and Mumford-Shah model | |
| Tricoche et al. | Fast and robust extraction of separation line features |
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 |