CN107423479A - The finite element modeling method of needling preform unit cell - Google Patents
The finite element modeling method of needling preform unit cell Download PDFInfo
- Publication number
- CN107423479A CN107423479A CN201710352165.9A CN201710352165A CN107423479A CN 107423479 A CN107423479 A CN 107423479A CN 201710352165 A CN201710352165 A CN 201710352165A CN 107423479 A CN107423479 A CN 107423479A
- Authority
- CN
- China
- Prior art keywords
- mrow
- acupuncture
- laying
- unit cell
- pinprick
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Nonwoven Fabrics (AREA)
Abstract
The finite element modeling method of needling preform unit cell, it is related to the modeling technique of needling preform, in order to meet the mechanical property of needle-punched composite, damage and the analysis of failure demand.The invention is determination acupuncture zone radius R, needle plate plants pin parameter, the thickness parameter of fiber composite material laying order and each laying, Needle-Punched Process Parameters, calculate the size of precast body unit cell, the coordinate position of the pinprick of each laying, establish the two-dimensional geometry model of precast body unit cell, two-dimensional grid unit is generated by two-dimensional geometry model, 3D solid unit is generated by two-dimensional grid unit, identify the overlapping cases in acupuncture region in identical laying and different layings, assign corresponding performance parameter, and assign corresponding performance parameter for the non-acupuncture region of unit cell, obtain the FEM model of needling preform unit cell.The present invention is applied to establish the FEM model of needling preform unit cell.
Description
Technical field
The present invention relates to the modeling technique of needling preform, and in particular to the finite element modeling side based on Needle-Punched Process Parameters
Method.
Background technology
Needling technique is a kind of pre- fiber preparation process, can be by pricker by the face in fiber composite material laying
Fiber is incorporated into thickness direction, forms the needling fiber beam of vertical distribution, each fiber laying is combined closely, strengthens precast body
Interlayer mechanical property.Fiber composite material mainly includes carbon cloth and net tire, and wherein carbon cloth ensure that precast body in laying
There is higher fiber volume fraction, net tire is the main source of interlayer fiber in plane.Needling preform is as a kind of unique
Tridimensional network, have the characteristics that fiber be evenly distributed, be easy to densification shaping.Needling technique overcomes two-dimensional fiber knot
The shortcomings that structure interlayer performance is weak, with braiding, woven compared with traditional textile technologies such as knittings with technique is simple, cost is low and property
Valency than it is high the advantages of, occupy staple market status in numerous prefabricated preparation process at present.
Needle-Punched Process Parameters are numerous, Parameter adjustable scope is big, and not similar shape can be prepared by adjusting needling process scheme
Shape, various sizes of carbon fiber needle punched fabric, including bulk, tabular, tubular, taper and complex-curved shape fiber preform etc..
Different needling densities, depth of needling, or even the fibre of variable cross-section, variable density can also easily and flexibly be obtained by changing acupuncture parameter
Dimensional fabric, embody the good designability of needling process.It can be prepared using acupuncture carbon fiber precast body as enhancing skeleton various
High performance acupuncture C/C and C/C-SiC composites, needle-punched composite are wide with its excellent mechanical property and thermophysical property
It is general to be applied to aerospace field, as aeroplane brake discses, solid propellant rocket outlet spouts, larynx lining and variously-shaped heat are prevented
Protect part etc..Also there is wide application space in many civil area needle-punched composites, such as various anti-heat insulation structurals, industry
High temperature stove and heat-resisting crucible etc..
The fibre structure of precast body has decisive influence to the mechanical property of composite, from meso-scale to composite wood
Expect modeling analysis, be the effective means of research material effective mechanics properties and damage process.This process is needed to composite
Fibre structure carry out careful observation and analysis, established according to fibre structure and represent elementary volume, volume element, then utilize component material
Performance obtain the overall mechanical properties of composite.At present, built for traditional textile composites such as braiding, woven and knittings
Comparative maturity, people can build vertical thin sight analysis model according to microscopic structural parameters such as fibre bundle size, angle of weave
It is vertical to represent elementary volume, volume element, study mechanical property and the derogatory behaviour of composite.But the characteristics of needling process itself, causes pin
It is sufficiently complex to pierce precast body fibre structure, and the work such as the fibre structure of needling preform and needling density, depth and arranging needle form
Skill parameter is directly related.On the one hand, by repeated puncture, the pinprick in precast body is intensive, unordered, and many positions are weighed
Multiple acupuncture；On the other hand, the carbon fiber in acupuncture region is in pinprick it is observed that obvious by a certain degree of damage
Fibrous fracture, deflection and migration, the geometric shape for being needled rear fiber are difficult to describe.Under macro-scale, needling preform is
A kind of network structure of uniform, the unordered arrangement of pinprick, but under micro- meso-scale, the fibre structure at local acupuncture position is very
Complexity, there is certain uncertainty.Therefore, the representative elementary volume, volume element for needle-punched composite being established based on morphology observations result is
Extremely difficult, the fibre structure of needle-punched composite, forecast material also up to the present can be described without carefully sight analysis model
The mechanical property of material, and the material damage based on needling preform microscopical structure and failure analysis work also have no report.
The content of the invention
The invention aims to meet mechanical property, damage and the analysis of the failure demand of needle-punched composite, from
And provide the finite element modeling method of needling preform unit cell.
The finite element modeling method of needling preform unit cell of the present invention, comprises the following steps：
Step 1: the fibre structure of scanning needling preform, obtains micro-image, to each acupuncture region in micro-image
Measure and be averaging, obtain the mean radius R in acupuncture region；
Step 2: the needle plate according to corresponding to needling preform, which obtains needle plate, plants pin parameter；
Step 3: according to micro-image, the laying order of the fiber composite material laying of needling preform and each paving are determined
The thickness parameter of layer；
Step 4: according to micro-image, Needle-Punched Process Parameters are determined；
Step 5: the parameter obtained according to step 2 and step 3 calculates the size of precast body unit cell, according to step 2 and
The parameter that step 4 obtains calculates the coordinate position of the pinprick of each laying；
Step 6: the coordinate position of the pinprick of each laying in step 5 is projected in z=0 plane, with reference to pre-
The size of body unit cell processed, the acupuncture region two-dimensional geometry model of mean radius size in establishment step one at each pinprick,
And then obtain precast body unit cell two-dimensional geometry model；Z-axis is perpendicular to laying；
Step 7: two-dimensional grid unit is generated by the precast body unit cell two-dimensional geometry model generated in step 6, with reference to
Depth of needling in the thickness and Needle-Punched Process Parameters of each laying, 3D solid unit, identification are generated by two-dimensional grid unit
The overlapping cases in acupuncture region in identical laying and different layings, corresponding performance parameter is assigned, and be the non-acupuncture area of unit cell
Domain assigns corresponding performance parameter, obtains the FEM model of needling preform unit cell.
Preferably, the needle plate that the step 2 obtains, which plants pin parameter, includes the spacing L of adjacent two rows pricker1With same row
The spacing L of interior neighbouring lancets2；
In step 5, the method that the parameter obtained according to step 2 and step 3 calculates the size of precast body unit cell is：
The L obtained according to step 21And L2Obtain the minimum period property arrangement range L of precast body surface pin holeU×LU, its
Middle LUFor L1And L2Least common multiple；The thickness parameter of each laying obtained according to step 3, obtain precast body unit cell total thickness
Spend LT；The size of precast body unit cell is LU×LU×LT。
Preferably, the Needle-Punched Process Parameters that the step 4 obtains are the relative position (x of needle plate and precast bodyij,
yij)；
In step 5, the parameter obtained according to step 2 and step 4 calculates the coordinate position of the pinprick of each laying
Method is：
The L obtained according to step 21And L2, and the relative position (x of the needle plate of step 4 and precast bodyij,yij) calculate often
Coordinate (the x of the pinprick of one layingijk, yijk, zijk)；xij,yijRespectively the x-axis of the i-th laying, jth time acupuncture and y-axis are sat
Mark；xijk, yijk, zijkRespectively time acupuncture of the i-th laying, jth, the x-axis of k-th pinprick, y-axis, z-axis coordinate；
Relative position (the x of needle plate and precast bodyij,yij) when starting acupuncture, sat with the origin of coordinates from nearest pinprick
Cursor position.
Preferably, in the step 7, the mean radius R in the acupuncture region obtained according to step 1, step 4 obtain
Needle-Punched Process Parameters in depth of needling and the obtained coordinate of the pinprick of each laying of step 5, judge the acupuncture region
It is whether overlapping with other acupuncture regions.
Preferably, it is described to judge whether the method overlapping with other acupuncture regions is in the acupuncture region：
For two acupuncture regions in same fiber composite material laying, the coordinate of pinprick is respectively (xijk, yijk,
zijk) and (xilm, yilm, zilm), two pinprick distances are for d：
If d<2R, then it is overlapping acupuncture region to judge the two acupuncture regions, is otherwise independent acupuncture region；
For two acupuncture regions in different fiber composite material layings, the coordinate of pinprick is respectively (xijk, yijk,
zijk) and (xlmn, ylmn, zlmn), the i-th laying is located on l layings, and two pinprick plane projection distances are that d ' is：
Two acupuncture region z to distance be d "
D "=| zijk-zlmn|-di (5)
If d '<2R, and d " ＜ 0, then it is overlapping acupuncture region to judge the two acupuncture regions, is otherwise independent pin
Spine area domain, diFor the depth of needling of the i-th laying.
The present invention can obtain corresponding precast body FEM model according to any needling process scheme, establish needling process
With the incidence relation between precast body fibre structure.Precast body single cell model, Ke Yifang are established by inputting various process parameters
Just quantity, acupuncture region relative distance, independent acupuncture region and the Repeated Acupuncture Domain Volume of pinprick in each laying of output
The quantitative informations such as content, intuitively reflect pinprick arrangement density, the uniformity of diverse location situations such as.The present invention can make
For the preferred of needling process scheme and the important supplementary means of design, adjacent acupuncture position can be further studied based on this
Mutually influence situation.According to the distributing position in acupuncture region, FEM model is divided into independent acupuncture region, Repeated Acupuncture
Region and non-acupuncture region, corresponding material properties are assigned to different region, and (including elastic performance and thermophysical property are joined
Number), so as to analyze the effective mechanics properties of needle-punched composite and failure damage process.
Brief description of the drawings
Fig. 1 is needling preform forming process schematic diagram；1 is needle plate, and 2 be pinprick, and 3 be i-th layer of fiber composite material, 4
For acupuncture region；
Fig. 2 is the structural representation of needle plate；5 be pricker, and 6 pierce to swing to hook, and dotted line frame is a cloth pin unit cell in figure；
Fig. 3 is the pinprick arrangement schematic diagram on needling preform surface；7 be needling preform；
Fig. 4 be prefabricated acupuncture once, rotated, then the pinprick arrangement schematic diagram after acupuncture once；
Fig. 5 be needling preform each laying in pinprick arrangement schematic diagram；
Fig. 6 be Fig. 5 each laying in pinprick project pinprick arrangement schematic diagram in z=0 plane；
Fig. 7 is the schematic diagram in cylindrical needle spine area domain；
Fig. 8 is the schematic diagram in two acupuncture regions in same compound material laying；
Fig. 9 is the schematic diagram in two acupuncture regions in different composite material laying；
Figure 10 is the FEM model of needling preform unit cell；
Figure 11 is the acupuncture area distribution position view of needling preform unit cell.
Embodiment
Embodiment：Present embodiment is illustrated with reference to Fig. 1 to Figure 11, the acupuncture described in present embodiment is prefabricated
The finite element modeling method of body unit cell, comprises the following steps：
Step 1: the fibre structure of scanning needling preform, obtains micro-image, to each acupuncture region in micro-image
Measure and be averaging, obtain the mean radius R in acupuncture region；
Step 2: the needle plate according to corresponding to needling preform, which obtains needle plate, plants pin parameter；
Step 3: according to micro-image, the laying order of the fiber composite material laying of needling preform and each paving are determined
The thickness parameter of layer；
Step 4: according to micro-image, Needle-Punched Process Parameters are determined；
Step 5: the parameter obtained according to step 2 and step 3 calculates the size of precast body unit cell, according to step 2 and
The parameter that step 4 obtains calculates the coordinate position of the pinprick of each laying；
Step 6: the coordinate position of the pinprick of each laying in step 5 is projected in z=0 plane, with reference to pre-
The size of body unit cell processed, the acupuncture region two-dimensional geometry model of mean radius size in establishment step one at each pinprick,
And then obtain precast body unit cell two-dimensional geometry model；Z-axis is perpendicular to laying；
Step 7: two-dimensional grid unit is generated by the precast body unit cell two-dimensional geometry model generated in step 6, with reference to
Depth of needling in the thickness and Needle-Punched Process Parameters of each laying, 3D solid unit, identification are generated by two-dimensional grid unit
The overlapping cases in acupuncture region in identical laying and different layings, corresponding performance parameter is assigned, and be the non-acupuncture area of unit cell
Domain assigns corresponding performance parameter, obtains the FEM model of needling preform unit cell.
In present embodiment, the needle plate plant pin parameter that step 2 obtains includes the spacing L of adjacent two rows pricker1With same row
The spacing L of interior neighbouring lancets2；
In step 5, the method that the parameter obtained according to step 2 and step 3 calculates the size of precast body unit cell is：
The L obtained according to step 21And L2Obtain the minimum period property arrangement range L of precast body surface pin holeU×LU, its
Middle LUFor L1And L2Least common multiple；The thickness parameter of each laying obtained according to step 3, obtain precast body unit cell total thickness
Spend LT；The size of precast body unit cell is LU×LU×LT。
In present embodiment, the Needle-Punched Process Parameters that step 4 obtains are needle plate and the relative position (x of precast bodyij,
yij)；
In step 5, the parameter obtained according to step 2 and step 4 calculates the coordinate position of the pinprick of each laying
Method is：
The L obtained according to step 21And L2, and the relative position (x of the needle plate of step 4 and precast bodyij,yij) calculate often
Coordinate (the x of the pinprick of one layingijk, yijk, zijk)；xij,yijRespectively the x-axis of the i-th laying, jth time acupuncture and y-axis are sat
Mark；xijk, yijk, zijkRespectively time acupuncture of the i-th laying, jth, the x-axis of k-th pinprick, y-axis, z-axis coordinate；
Relative position (the x of needle plate and precast bodyij,yij) when starting acupuncture, sat with the origin of coordinates from nearest pinprick
Cursor position.
In present embodiment, in step 7, the mean radius R in the acupuncture region obtained according to step 1, step 4 obtain
Needle-Punched Process Parameters in depth of needling and the obtained coordinate of the pinprick of each laying of step 5, judge the acupuncture region
It is whether overlapping with other acupuncture regions.
In present embodiment, judge whether the method overlapping with other acupuncture regions is in the acupuncture region：
For two acupuncture regions in same fiber composite material laying, the coordinate of pinprick is respectively (xijk, yijk,
zijk) and (xilm, yilm, zilm), two pinprick distances are for d：
If d<2R, then it is overlapping acupuncture region to judge the two acupuncture regions, is otherwise independent acupuncture region；
For two acupuncture regions in different fiber composite material layings, the coordinate of pinprick is respectively (xijk, yijk,
zijk) and (xlmn, ylmn, zlmn), the i-th laying is located on l layings, and two pinprick plane projection distances are that d ' is：
Two acupuncture region z to distance be d "
D "=| zijk-zlmn|-di (5)
If d '<2R, and d " ＜ 0, then it is overlapping acupuncture region to judge the two acupuncture regions, is otherwise independent pin
Spine area domain, diFor the depth of needling of the i-th laying.
Specific embodiment：
Step 1: being scanned observation to needling preform by surface sweeping Electronic Speculum SEM, pinprick is measured by displaing micro picture
Size.Pinprick in precast body is formed in acupuncture course.It is dispersed with pricker and swings to hook thorn, pricker is pierced into laying
Fiber is even fractureed to thickness direction deflection, migration in many faces afterwards, and after pricker is extracted, acupuncture is left on precast body
Hole.Due to chopped strand random arrangement in net tire compound material, and there are uncertain, different pinpricks in fibrous fracture elongation
Size and fibre structure it is different.It is pin to define the pinprick region that nearby fiber deflects, migrates or fractureed in face
Spine area domain.Needling preform surface is observed, it is found that acupuncture region is generally circular；Needle-punched composite is cut into slices
Observation, it is found that fiber is in " people " font in the inner face of acupuncture region.By the way that substantial amounts of observed image is measured and counted
It is R (mm) to acupuncture zone leveling radius.In the present embodiment, acupuncture region is reduced to the cylindrical region that radius is R；
Step 2: according to needle plate, the plant pin parameter of pricker on needle plate is determined.Pricker is according to certain rows of installation of spacing
On needle plate, as shown in Figure 2.Needle plate width is Wp(cm), the spacing of adjacent two rows pricker is L1(cm), adjacent thorn in same row
The spacing of pin is L2(cm) the minimum period property arrangement scope of pricker on needle plate, can be obtained according to above-mentioned size；
Step 3: according to displaing micro picture, the laying parameter of needling preform is determined.It is first in the forming process of needling preform
First by compound material according to 0 ° without the laminated structures of dimension cloth, net tire and 90 ° without dimension cloth together.Often lay certain thickness compound
Material, in precast body surface acupuncture several times, then proceed to lay compound material, continue acupuncture, move in circles progress, until precast body
Reach predetermined thickness and needling density.O-xyz coordinate systems, 0 ° and 90 ° of machine directions point are established along compound material long fibre direction
X directions and y directions are not corresponded to, and thickness direction corresponds to z directions, and it is T to define i-th layer of fiber composite material overlay thicknessi(cm), as schemed
Shown in 1.
Step 4: according to displaing micro picture, Needle-Punched Process Parameters are determined, needling preform is obtained with reference to Needle-Punched Process Parameters
Unit cell size and the spatial distribution position for calculating unit cell internal needle spine area domain.In modeling process it needs to be determined that Needle-Punched Process Parameters
It is as follows：
(1) stepping-in amount s (cm), represent per acupuncture once, the distance that fiber composite material advances：
Wherein v (cm/min) is cloth outputting speed, and f (pin/min) is punch frequency；
(2) acupuncture frequency ni, for the acupuncture number after the i-th laying fiber composite material of laying, (compound material leads to below needle plate
Cross and once represent acupuncture once)；
(3) needling density Dn(pin/cm2), represent unit area acupuncture hole number in precast body：
Dn=Mp∑ni (2)
Wherein Mp(pin/cm2) it is bestock, represent the quantity of pricker installed on needle plate in unit area；
The depth of needling d of (4) i-th layers of i-th layingi(cm), as shown in Figure 1；
(5) when starting acupuncture, the relative position (x of needle plate and precast bodyij,yij), that is, when starting acupuncture with the origin of coordinates from
Nearest pinprick coordinate position, wherein i represent the i-th laying, and j represents jth time acupuncture, as shown in Figure 3.
Step 5: according to Step 2: precast body minimum period property is calculated in the technological parameter in step 3 and step 4
Unit cell size.In acupuncture course, in order to ensure the uniformity in precast body x directions and the arrangement of y directions pin hole, after every acupuncture once all
By precast body in 90 ° of x-y faces internal rotation, acupuncture is then proceeded to.Because pricker according to certain spacing is rows of is arranged on pin
On plate, therefore the pin hole row of precast body surface is in periodic arrangement after reciprocal acupuncture.According to the plant pin spacing L on needle plate1
And L2The minimum period property arrangement scope that precast body surface pin hole is calculated is LU×LU(cm2), wherein LUFor L1And L2Most
Small common multiple.Further, according to precast body gross thickness LT(LT=Σ Ti) to determine precast body unit cell size be LU×LU×LT
(cm3), TiFor the thickness of the i-th laying.According to plant pin spacing L1、L2With needle plate and the relative position (x of precast bodyij,yij) calculate
Coordinate (the x of pinprick into each layingijk, yijk, zijk), wherein i represents the i-th laying, and j represents jth time acupuncture, and k represents the
K pinprick, (it is illustrated as shown in Figure 5 in figure by taking the needling preform of three layers of compound material composition as an example)；
Step 6: pinprick coordinate position in each laying is projected in z=0 plane, existed by Python
Precast body unit cell and the two-dimensional geometry model in acupuncture region are established in finite element software Abaqus, as shown in Figure 6.Precast body surface
The minimum period property arrangement range L of pin holeU×LU, acupuncture region geometry model radius is R；
Step 7: mesh generation generation two-dimensional grid unit, grid cell are carried out to the two-dimensional geometry model in step 6
Type uses 3 node shell units (i.e. two-dimentional S3 units).Ordered using " generateMeshByOffset " in Abaqus, with
Based on two-dimensional cell, stretched generation 3D solid unit (i.e. C3D6 units), stretching gross thickness is the total thickness of precast body
Spend LT.According to the coordinate (x of the pinprick of each layingijk, yijk, zijk), acupuncture zone radius R and depth of needling di, utilize
" getByBoundingCylinder " order in Abaqus chooses the unit in each acupuncture region in different layings and establishes unit
Set.Illustrated by taking k-th of pinprick that the i-th laying, jth time acupuncture are formed as an example, the acupuncture region includes circle shown in Fig. 7
All units in the range of cylinder, cylinder radius R, it is highly di, upper and lower central coordinate of circle is respectively (xijk, yijk, zijk) and
(xijk, yijk, zijk-di).Position coordinates according to acupuncture region is that the coordinate of pinprick judges whether adjacent needles spine area domain weighs
Close, if misaligned, it is independent acupuncture region to judge the acupuncture region；If overlapped, it is Repeated Acupuncture to judge the acupuncture region
Region；Judge whether any two acupuncture region repeats specifically to carry out in two kinds of situation as follows：
(1) for two acupuncture regions in same compound material laying, pin hole coordinate is respectively (xijk, yijk, zijk)
(xilm, yilm, zlmn), two pinprick distances are d
If d<2R, then the two acupuncture regions are judged for Repeated Acupuncture region, as shown in figure 8, being otherwise independent pin
Spine area domain；
(2) for two acupuncture regions in different composite material laying, pin hole coordinate is respectively (xijk, yijk, zijk) and
(xlmn, ylmn, zlmn), two pinprick plane projection distances are d '
Acupuncture region z to distance be d "
D "=| zijk-zlmn|-di (5)
If d '<2R, and d "<0, then the two acupuncture regions are judged for Repeated Acupuncture region, as shown in figure 9, no
It is then independent acupuncture region；
The needling preform unit cell FEM model of any technological parameter can be established according to the present invention.With three layers of compound material
Illustrated exemplified by the needling preform being formed by stacking, precast body molding technique parameter as shown in Table 1 and Table 2, the list of the coordinate of table 2
Position is mm：
The needling preform molding technique parameter of table 1
| Parameter | L1 | L2 | S | WP | R | T1 |
| Numerical value | 0.5cm | 0.6cm | 1.5cm | 1.5cm | 0.6mm | 0.3cm |
| Parameter | T2 | T3 | d1 | d2 | d3 | |
| Numerical value | 0.3cm | 0.3cm | 0.3cm | 0.6cm | 0.6cm |
The relative position coordinates of needle plate and precast body in the acupuncture course of table 2
| First layer | (X11, Y11)=(00) | (X12, Y12)=(3,0) |
| The second layer | (X21,Y21)=(0,4) | (X22,Y22)=(3,4) |
| Third layer | (X31,Y31)=(0,8) | (X32,Y32)=(3,8) |
Finally, needling preform unit cell FEM model is divided into three parts：Non- acupuncture region, independent acupuncture region and
Repeated Acupuncture region, as shown in FIG. 10 and 11.Calculate the volume integral in each region in FEM model automatically by Python
Number, corresponding performance parameter is assigned to different zones, Equivalent Mechanical performance and failure procedure of unit cell etc. can be calculated.
Present embodiment carries out substantial amounts of microscopic observation using SEM scan images to needling preform, passes through the side of statistics
Method obtains the mean radius in acupuncture region, and acupuncture position is reduced into cylindrical region in modeling process.In FEM model
Acupuncture area size can represent the average value of acupuncture position actual size in material.It is pre- that present embodiment is based on actual acupuncture
Make body formed technological parameter to be modeled, the forming process of precast body, resulting unit cell are accurately simulated by numerical method
The spatial arrangement position in acupuncture region is consistent with the regularity of distribution at acupuncture position in real material in model.
Present embodiment can obtain the finite element unit cell of different technical parameters needling preform by adjusting modeling parameters
Model：Acupuncture number by adjusting each laying compound material can realize needling density different in unit cell, not even with paving
The pinprick arrangement of layer variable density；Relative position when adjusting each acupuncture between needle plate and compound material laying can simulate pin hole
The different arrangement form in precast body, i.e., different " needle trackings ", simulate uniform, the scattered arrangement in acupuncture region in unit cell or concentrate
The different situations such as arrangement；The multiple depth of needling of adjustment can realize acupuncture region in the different arrangement shape of single cell model thickness direction
Formula, realize that pricker pierces through the simulation with the different situations such as puncture completely to compound material laying part.
Claims (5)
1. the finite element modeling method of needling preform unit cell, it is characterised in that comprise the following steps：
Step 1: the fibre structure of scanning needling preform, obtains micro-image, each acupuncture region in micro-image is carried out
Measurement and averaging, obtain the mean radius R in acupuncture region；
Step 2: the needle plate according to corresponding to needling preform, which obtains needle plate, plants pin parameter；
Step 3: according to micro-image, the laying order of the fiber composite material laying of needling preform and each laying are determined
Thickness parameter；
Step 4: according to micro-image, Needle-Punched Process Parameters are determined；
Step 5: the parameter obtained according to step 2 and step 3 calculates the size of precast body unit cell, according to step 2 and step
Four obtained parameters calculate the coordinate position of the pinprick of each laying；
Step 6: the coordinate position of the pinprick of each laying in step 5 is projected in z=0 plane, with reference to precast body
The size of unit cell, the acupuncture region two-dimensional geometry model of mean radius size in establishment step one at each pinprick, and then
Obtain precast body unit cell two-dimensional geometry model；Z-axis is perpendicular to laying；
Step 7: two-dimensional grid unit is generated by the precast body unit cell two-dimensional geometry model generated in step 6, with reference to each
Depth of needling in the thickness and Needle-Punched Process Parameters of laying, 3D solid unit is generated by two-dimensional grid unit, identification is identical
The overlapping cases in acupuncture region in laying and different layings, corresponding performance parameter is assigned, and assigned for the non-acupuncture region of unit cell
Corresponding performance parameter is given, obtains the FEM model of needling preform unit cell.
2. the finite element modeling method of needling preform unit cell according to claim 1, it is characterised in that
The needle plate that the step 2 obtains, which plants pin parameter, includes the spacing L of adjacent two rows pricker1And in same row between neighbouring lancets
Away from L2；
In step 5, the method that the parameter obtained according to step 2 and step 3 calculates the size of precast body unit cell is：
The L obtained according to step 21And L2Obtain the minimum period property arrangement range L of precast body surface pin holeU×LU, wherein LUFor
L1And L2Least common multiple；The thickness parameter of each laying obtained according to step 3, obtain precast body unit cell gross thickness LT；
The size of precast body unit cell is LU×LU×LT。
3. the finite element modeling method of needling preform unit cell according to claim 2, it is characterised in that
The Needle-Punched Process Parameters that the step 4 obtains are the relative position (x of needle plate and precast bodyij,yij)；
In step 5, the method for the coordinate position of the pinprick of each laying of parameter calculating obtained according to step 2 and step 4
For：
The L obtained according to step 21And L2, and the relative position (x of the needle plate of step 4 and precast bodyij,yij) calculate each paving
Coordinate (the x of the pinprick of layerijk, yijk, zijk)；xij,yijThe respectively x-axis and y-axis coordinate of the i-th laying, jth time acupuncture；
xijk, yijk, zijkRespectively time acupuncture of the i-th laying, jth, the x-axis of k-th pinprick, y-axis, z-axis coordinate；
Relative position (the x of needle plate and precast bodyij,yij) i.e. start acupuncture when, with the origin of coordinates from nearest pinprick coordinate bit
Put.
4. the finite element modeling method of needling preform unit cell according to claim 1, it is characterised in that the step 7
In, depth of needling in the Needle-Punched Process Parameters that the mean radius R in the acupuncture region obtained according to step 1, step 4 obtain and
The coordinate of the pinprick for each laying that step 5 obtains, judge whether the acupuncture region is overlapping with other acupuncture regions.
5. the finite element modeling method of needling preform unit cell according to claim 4, it is characterised in that the judgement should
Whether the method overlapping with other acupuncture regions is in acupuncture region：
For two acupuncture regions in same fiber composite material laying, the coordinate of pinprick is respectively (xijk, yijk, zijk)
(xilm, yilm, zilm), two pinprick distances are for d：
<mrow>
<mi>d</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
If d<2R, then it is overlapping acupuncture region to judge the two acupuncture regions, is otherwise independent acupuncture region；
For two acupuncture regions in different fiber composite material layings, the coordinate of pinprick is respectively (xijk, yijk, zijk)
(xlmn, ylmn, zlmn), the i-th laying is located on l layings, and two pinprick plane projection distances are that d ' is：
<mrow>
<msup>
<mi>d</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>l</mi>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>l</mi>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Two acupuncture region z to distance be d "
D "=| zijk-zlmn|-di (5)
If d '<2R, and d " ＜ 0, then it is overlapping acupuncture region to judge the two acupuncture regions, is otherwise independent acupuncture area
Domain, diFor the depth of needling of the i-th laying.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710352165.9A CN107423479A (en) | 2017-05-18 | 2017-05-18 | The finite element modeling method of needling preform unit cell |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710352165.9A CN107423479A (en) | 2017-05-18 | 2017-05-18 | The finite element modeling method of needling preform unit cell |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107423479A true CN107423479A (en) | 2017-12-01 |
Family
ID=60425908
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710352165.9A Pending CN107423479A (en) | 2017-05-18 | 2017-05-18 | The finite element modeling method of needling preform unit cell |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107423479A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108416139A (en) * | 2018-03-06 | 2018-08-17 | 武汉凌云建筑装饰工程有限公司 | Big plate width Pierced Aluminum Plank Equivalent calculation method |
| CN112080387A (en) * | 2020-08-25 | 2020-12-15 | 南开大学 | Cell nucleus operation method based on cell nucleus position dynamic drift modeling |
| WO2021004513A1 (en) * | 2019-07-11 | 2021-01-14 | 南京航空航天大学 | Ceramic matrix composite bolt prefabricated body-structure integrated design method |
| CN114492009A (en) * | 2022-01-17 | 2022-05-13 | 北京理工大学 | Fiber path determination and performance calculation method of needle-punched twill cloth composite material |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106066913A (en) * | 2016-05-31 | 2016-11-02 | 西北工业大学 | Complex composite material structure equivalent material performance multi-dimension computational methods |
-
2017
- 2017-05-18 CN CN201710352165.9A patent/CN107423479A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106066913A (en) * | 2016-05-31 | 2016-11-02 | 西北工业大学 | Complex composite material structure equivalent material performance multi-dimension computational methods |
Non-Patent Citations (6)
| Title |
|---|
| JUNBO XIE ET.AL: "Effect of needling parameters on the effective properties of 3D needled C/C-SiC composites", 《COMPOSITES SCIENCE AND TECHNOLOGY》 * |
| 张泽书: "汽车内饰用GMT复合材料的制备与研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑（月刊）》 * |
| 张盼: "三维针刺C_C复合材料的微结构建模及力学性能预测", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑（月刊）》 * |
| 徐健等: "非织造针刺机针板的计算机设计方法及分析", 《产业用纺织品》 * |
| 胡吉永等: "《纺织结构成型学 2 多维成形》", 31 March 2016 * |
| 赵宁涛等: "计算机高级语言在针刺机针板杂乱布针中的应用", 《非织造布》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108416139A (en) * | 2018-03-06 | 2018-08-17 | 武汉凌云建筑装饰工程有限公司 | Big plate width Pierced Aluminum Plank Equivalent calculation method |
| WO2021004513A1 (en) * | 2019-07-11 | 2021-01-14 | 南京航空航天大学 | Ceramic matrix composite bolt prefabricated body-structure integrated design method |
| CN112080387A (en) * | 2020-08-25 | 2020-12-15 | 南开大学 | Cell nucleus operation method based on cell nucleus position dynamic drift modeling |
| CN112080387B (en) * | 2020-08-25 | 2021-10-01 | 南开大学 | Cell nucleus operation method based on cell nucleus position dynamic drift modeling |
| CN114492009A (en) * | 2022-01-17 | 2022-05-13 | 北京理工大学 | Fiber path determination and performance calculation method of needle-punched twill cloth composite material |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107423479A (en) | The finite element modeling method of needling preform unit cell | |
| Gereke et al. | A review of numerical models for 3D woven composite reinforcements | |
| Isart et al. | Internal geometric modelling of 3D woven composites: A comparison between different approaches | |
| US10303157B2 (en) | Additive manufacturing method for three-dimensional object | |
| DE112013001946B4 (en) | Emission signal visualization device | |
| CN108859128B (en) | The path planning fill method of complicated more hole regions | |
| CN106127856B (en) | The method that the strata model of column containing crack is made based on CT scan and 3D printing | |
| Xia et al. | Stress-based tool-path planning methodology for fused filament fabrication | |
| CN103752823B (en) | Triangular mesh type laser scanning method for selective laser sintering | |
| Volpato et al. | The influence of support base on FDM accuracy in Z | |
| Soltani et al. | Prediction of permeability of realistic and virtual layered nonwovens using combined application of X-ray μCT and computer simulation | |
| Abdiwi et al. | Characterising and modelling variability of tow orientation in engineering fabrics and textile composites | |
| CA2815968A1 (en) | Novel composite parts, methods and apparatus for manufacturing the same | |
| CN106553345B (en) | A kind of Method of printing and print control unit of more material 3D objects | |
| CN113192576B (en) | Modeling method of 2.5D woven composite material single cell model | |
| CN105774093B (en) | A kind of gradient-structure punctures fabric and preparation method thereof | |
| CN110348429B (en) | Method for identifying warps and wefts of plain weave structure | |
| Soltani et al. | 3D fiber orientation characterization of nonwoven fabrics using X-ray micro-computed tomography | |
| CN109190167A (en) | The lateral microstructure generation method of unidirectional long fiber reinforcement composite material | |
| CN107464285A (en) | The Meshing Method and device of a kind of threedimensional model | |
| CN110660129A (en) | Three-dimensional orthogonal fabric micro-geometric structure modeling method based on digital unit method | |
| Shan et al. | Additive manufacturing of non-planar layers using isothermal surface slicing | |
| Sajjad et al. | Impact of multiple infill strategy on the structural strength of single build FDM printed parts | |
| Eskandarnia et al. | Effect of fabric structure on in-plane and through-plane hydraulic properties of nonwoven geotextiles | |
| CN113239433A (en) | Steel bridge deck pavement layer rolling construction simulation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171201 |
|
| RJ01 | Rejection of invention patent application after publication |