Composites & Materials

Woven composite simulation from yarn unit cell and draping to part stiffness

The preform wrinkles on the diagonal, the stress model assumes 0/90 everywhere, and the 3D weave coupons do not say what the part will deliver. Unit cells, draping shear maps and a tested hand-over show which weave or angle is responsible.

Left, a 3D woven fabric block with a stepped edge showing layers of woven yarns. Right, a schematic of warp and weft yarns with a dashed unit cell, a draped-preform shear-angle map and a curve of shear force against shear angle with a locking point.
A 3D woven fabric block beside schematics of a woven unit cell, a draped shear-angle map and a locking curve.

Woven composite simulation predicts what a weave delivers, stiffness, damage onset and interlaminar behaviour, from a yarn-level unit cell, and what draping does to fibre angle in the real part. We build both as one chain around your existing solver, with a validation record and a surrogate of stated range, so more weave, layup and tooling questions are answered before the next coupon campaign or tool trial.

The preform drapes cleanly at the edges and wrinkles on the 45-degree diagonal, while the stress model has warp at 0 and weft at 90 degrees everywhere. The coupons came from a flat panel. Qian and colleagues report that woven prepreg formability is dominated by intraply shear and interply friction, which change fibre orientation. For a 3D weave, Friemann and colleagues note that test campaigns can be prohibitively expensive. Waiting on it: the allowables plan, preform trials, scrap, the weight target and approval dates.

What does woven composite simulation have to predict?

Three results at three scales. The unit cell, a yarn-level model of one repeat of the weave, gives stiffness and, with the right physics, damage onset. The draping model gives each point of the preform its warp and weft directions, the angle between them and the thickness. The structural model receives both; see also our article on choosing the modelling scale.

How does a yarn unit cell give stiffness, and where does it go wrong?

Mesh one repeat of the weave, tie opposite faces with periodic boundary conditions, and run six unit-strain cases, three tensions and three shears; averaged stress over strain fills the stiffness matrix, as Friemann and colleagues do. Each yarn is a transversely isotropic continuum along its path, with properties from fibre and matrix micromechanics, as in Feng and colleagues. Four choices decide the answer.

  • Crimp. A loaded wavy yarn tries to straighten, pressing on the yarns it wraps and lowering in-plane stiffness. Feng and colleagues attribute plain weave's inferior moduli to its maximum number of undulations; of 9,000 random single-material 6-by-6 patterns, 28.7% beat it in both tensile moduli.
  • Geometry. Friemann and colleagues call purely geometric generators far too idealised for as-manufactured material and weaving-process models costly per cell; a cell from an X-ray scan carries the real yarn paths. We use idealised cells for sweeps, scanned ones for validation and scatter.
  • Nesting. A periodic cell knows nothing of neighbouring layers. Jackson and Portanova describe 2D layers as expected to nest snugly; how they nest sets thickness, fibre content and resin-rich pockets, so we model the stack.
  • Mesh. Voxel meshes from scans automate the work, but Friemann and colleagues note artificial stress concentrations that can erroneously initiate damage; a conformal mesh needs every yarn separated first.
A woven unit cell: crimped yarn path, interlacing pattern, periodic faces, layer nesting, mesh and homogenised stiffness Top left, a section along one warp yarn in a plain weave layer: the warp yarn, drawn as a thick cyan band, passes over one weft yarn and under the next, and the weft yarns are cut across as ellipses inside the matrix. Under tension along the yarn, amber arrows show the wavy yarn pressing down on the weft yarn it crosses over and up on the one it crosses under, because it tries to straighten toward the dashed straight line. A dashed rectangle marks one repeat of the weave, the unit cell, with two pulsing rings on its left and right faces joined by an arrow to show that opposite faces are tied by periodic boundary conditions. Top right, two small weave diagrams, a plain weave with interlacing at every crossing and a satin with interlacing once in five, and below each the yarn profile, many undulations against few. Bottom left, two layers seen in section that alternate between sitting in phase and nesting into each other, which changes the stack thickness. Bottom middle, one yarn cross-section meshed two ways: curved conformal elements on the left half and a stair-step voxel mesh on the right half, with amber spots at the re-entrant corners where artificial stress concentrations arise. Bottom right, three horizontal bands for the warp, weft and in-plane shear stiffness, each with cyan dots for individual unit cells scattered inside the band and an amber ringed dot for the coupon test the band is checked against. The drawing teaches that a woven composite stiffness comes from the real yarn path and pattern, from how the cell is closed and meshed, and from a spread that is checked against test. No values are shown. WARP YARN IN SECTION INTERLACING STACKED LAYERS YARN MESH STIFFNESS WEFT CRIMP 1 PERIODIC FACES UNIT CELL 3 PLAIN SATIN MORE WAVES FEWER 2 IN PHASE NESTED 4 CONFORMAL VOXEL 5 E1 E2 G12 CELLS COUPON TEST 6
  1. Crimp is what makes a weave a weave. The warp yarn passes over one weft yarn and under the next, so its path undulates. Loaded along its length it tries to straighten toward the dashed line. That presses it onto the weft yarns it wraps, the amber arrows, loads the matrix through the thickness and lowers the in-plane stiffness. The cell has to carry that path, not a straight fibre.
  2. The pattern sets how often a yarn undulates. A plain weave interlaces at every crossing, a satin once in several. The yarn profile under each grid shows many waves against few, which is why two weaves of the same yarn and resin differ in stiffness. The pattern is a design variable, so it is a parameter of the cell.
  3. A cell is one repeat with matched faces. The dashed box is one repeat of the weave. Periodic boundary conditions tie the pulsing points on opposite faces so the cell behaves like a piece of infinite fabric. Three tensions and three shears applied in turn give the full stiffness matrix, and the faces must match in geometry as well as in displacement.
  4. Layers nest, and nesting changes the stack. A cell with periodic faces knows nothing of the layers around it. Neighbouring layers can sit in phase, or nest into each other, and the picture alternates between the two. Thickness, fibre volume fraction and the resin-rich pockets between layers all change, so the stack is modelled, not one layer.
  5. The mesh decides what the yarn surface does to stress. A cell built from a scan can be meshed as voxels, the right half, which automates the work but leaves a staircase. The amber spots are re-entrant corners where artificial stress concentrations can start damage the real yarn does not have. A conformal mesh, the left half, follows the surface but needs each yarn separated first.
  6. Stiffness arrives with a spread, and a test to meet. Each unit cell gives one point per stiffness component, the cyan dots, and cells differ because real yarn paths differ. The amber dot is the coupon measurement the band has to be checked against. The structural model receives the band and its range of validity, not a single number.
From yarn path to stiffness: crimp, pattern, periodic faces, nesting, mesh and spread. Illustrative schematic.

On a 3D-woven layer-to-layer angle interlock carbon composite, six cells of 6.5 by 21.5 by 4 mm from one scan gave elastic stiffnesses within 6.5% of experiments, with very little scatter between cells; damage remains future work for the authors. That is evidence for stiffness, none for strength.

Do 3D woven composites need a different interlaminar model?

In a 3D weave the warp yarns cross several weft layers, so reinforcement sits where a laminate has only resin. In a NASA Langley curved-beam test, Jackson and Portanova found 2D materials delaminated between layers while 3D weaves failed earlier by radial cracks between surface plies, which they trace to crack formation, not low through-thickness strength. In indentation, an orthogonal interlock was among the most damage-resistant materials.

A cohesive interface tuned to a laminate's delamination toughness knows nothing of a yarn crossing the crack. We model the through-thickness yarns explicitly, or fit a bridging law to your delamination tests, and say which; see composite modelling approaches. Friemann and colleagues note that heterogeneity of tens of millimetres forces coupon methods to be adjusted so the gauge holds a representative volume; a model-to-test comparison must state how many cells it holds. Fatigue is checked separately: composites fatigue simulation.

How does draping simulation give the local shear angle, and what is locking?

Yarns are quasi-inextensible, so fabric covers a doubly curved surface by trellising: the angle between warp and weft changes by the shear angle γ and cell area falls with cos γ, raising local thickness or fibre content. Boisse and colleagues review the two tests. The picture frame is highly constrained and can pick up parasitic fibre tension; the bias-extension test avoids that but assumes no warp-weft slip. A benchmark among seven groups was needed to agree good practice, and prepreg shear is measured at forming temperature.

Locking is where yarns touch and shear stiffness climbs; wrinkles tend to start. The review warns it is not a clean limit: wrinkling involves all the stresses, and angles above 60 degrees have been measured at high binder pressure without wrinkles. Fabric bends far more easily than shell theory predicts, so draping shells decouple bending from membrane stiffness, and friction between plies of different orientation dominates: with zero interply friction, simulations showed no wrinkles. Per Qian and colleagues, kinematic models that map a flat mesh onto the surface omit material and process.

Qian and colleagues validated a finite-element model on a four-ply twill carbon prepreg double dome. In the 0/90° ply the largest shear angle, 43.7°, sat in the dome-to-flat transition. Predicted angles differed from CT-tracer measurements by less than 14.1% at four positions there and less than 19.7% at three in the ±45° plies, and a wrinkle appeared in both simulation and scan. We carry that spread forward as input uncertainty. In explicit runs we keep kinetic energy small against internal energy, so mass scaling cannot create wrinkles.

Draping a woven fabric: trellis shear, the shear map on a dome, locking, wrinkling and the hand-over to the structural mesh Top left, a flat square patch of woven yarns above the same patch sheared into a parallelogram: the yarns keep their length, the angle between warp and weft changes by the shear angle, the cell area falls with the cosine of that angle and the yarns bunch closer together, drawn with a denser fill. Top middle, a circular top view of a draped dome with a distorted grid; four amber lobes of high shear sit on the diagonals near the rim and the centre and the axes stay low, and short blinking amber ridges mark wrinkles at one lobe. Top right, a curve of shear force against shear angle that rises gently, then steeply at a dashed locking line, with the region beyond it shaded amber for wrinkles, and a dot travelling along the curve. Bottom left, a side section of two plies lying on a tool: holder arrows press the ends down, arrows in the two plies point in opposite directions to show slip, a zig-zag marks friction at the interface, and the upper ply buckles into a wrinkle. Bottom middle, a fine draping mesh with a coarser structural mesh drawn over it; in each structural element a pair of cross marks shows the mapped warp and weft directions, tilting further from a right angle as the shear grows. Bottom right, a sheared unit cell with an arrow to a plate of four elements, one of them flagged in amber as a wrinkle. The drawing teaches that draping is a trellis deformation whose shear map, locking and wrinkles depend on holder force and ply friction, and that the fibre directions it produces must be mapped onto the structural mesh with a stated rule. No values are shown. FABRIC SHEARS DOME SHEAR MAP SHEAR RESPONSE PLIES AND HOLDER DRAPE TO STRUCTURE PER-ELEMENT INPUT γ FLAT SHEARED cos γ CELL AREA 1 LOW HIGH 2 LOCKING WRINKLES SHEAR ANGLE FORCE 3 0/90 ±45 4 DRAPE STRUCTURAL 5 ! SHEARED CELL ELEMENTS WRINKLE FLAG 6
  1. Yarns do not stretch, so the fabric shears. To cover a doubly curved surface the fabric trellises: warp and weft stay straight and the angle between them changes by the shear angle, the amber arc. The area of each cell falls with the cosine of that angle, so the same fibre packs into less area and the local thickness or fibre content rises.
  2. Shear is not uniform. On a dome draped with 0/90 fabric the large angles sit on the diagonals, the amber lobes, where the surface demands the most trellising, and the centre and the axes stay low. This map is a result of the draping model, and it is the field that has to be handed on.
  3. Locking: stiffness climbs and wrinkles follow. As the angle grows the yarns come into contact, the in-plane shear stiffness rises steeply and wrinkles tend to start, the shaded zone. The locking angle is not a clean limit: wrinkling involves all the stresses in the fabric, so holder tension and friction can move it.
  4. Holder force and ply friction decide the wrinkle. Plies of different orientation slip against each other, the opposed arrows. Friction turns the slip into compression and the upper ply buckles, while holder force adds tension that can suppress it. A draping model without both inputs predicts the wrong wrinkles, so they are calibrated.
  5. The mapping rule is part of the model. The draping mesh is fine and the structural mesh is coarser. Each structural element receives a warp and a weft direction, a pair no longer at a right angle, drawn as the crossed marks, and how angles are averaged across the finer mesh is a modelling choice with its own error near gradients.
  6. Per-element input, and wrinkles flagged. Each element carries its fibre angle and local fibre content into a stiffness taken from a sheared unit cell or from tests. Where the draping model predicts a wrinkle the element is flagged, because a smooth structural mesh does not contain one.
Draping: trellis shear, the shear map, locking, wrinkles and the hand-over to the structural mesh. Illustrative schematic.

How do draped fibre angles reach the structural model?

Per element the draping model gives two direction vectors and a thickness, not one ply angle: after shear, warp and weft are not orthogonal. The structural mesh is coarser, so the mapping rule, nearest element, area-weighted, angles or vectors interpolated, is a modelling choice with an error near gradients. Properties then come from sheared unit cells or tests as a function of shear angle and fibre content, and wrinkles are flagged, not modelled. We script it, test it on patches with known angle fields and report the stiffness difference from the nominal 0/90 model.

Which woven composite question needs which model?

Woven composite simulation models and the decision each one supports
QuestionModelReportsLimit to state
What stiffness does this weave give?Unit cell, periodic faces, six strain casesStiffness by componentYarn data, geometry, nesting
How scattered is the as-made material?Several scanned or randomised cellsSpread per componentElastic only; strength needs damage
Where does the preform shear and wrinkle?Draping shell: shear, bending, frictionShear map, wrinkle flagsCharacterisation scatter
What does the sheared fabric do in the part?Mapped warp and weft directionsLocal stiffness by elementMapping rule; wrinkles not modelled
Does a 3D weave resist delamination?Cell with through-thickness yarns, bridging lawCrack growth resistanceLaw fitted to your tests
Sweeps over weave, yarn size, layupSurrogate on designed cell runsModuli inside the sampled rangeRange only; gate to full cell

A validated unit cell is not a design allowable

Simulation does not replace the building-block test programme, coupon to component, including impact and compression-after-impact tests, or the certifier's, customer's or authority's own process and sign-off. FAA Advisory Circular 20-107B describes building-block tests and analyses at coupon, element, detail and subcomponent level, and sets design values for non-laminated forms at the scale that represents the material in the part, or by tests of material substructure with a test-validated analytical method. Predicted stiffness, strength, shear angles and wrinkles are estimates that depend on yarn and matrix data and scatter, geometry, mesh, draping characterisation and calibration. How far virtual evidence is accepted is agreed with the certifier, authority or customer and follows the risk of the decision. Nothing we build certifies, qualifies or approves anything, generates certified allowables or guarantees strength, damage tolerance, crash performance or life.

Where does AI help woven composite design, and where does it not?

  • Surrogates over weave and yarn parameters. Feng and colleagues solved repeated unit cells of 9,000 random 6-by-6 patterns for in-plane E1, E2 and G12 and trained networks to predict them and propose patterns for a target modulus. Different patterns can share a modulus, so the inverse problem is many-to-one: we return candidate sets.
  • Scan segmentation. Friemann and colleagues train a 3D U-net on fully synthetic labelled data; once trained, a volume becomes a finite element input file in under a minute, with unit cells still picked by hand.
  • Physics-encoded surrogates. Ghane and colleagues note that data-driven surrogates for woven composites often need big data, lack interpretability and extrapolate poorly, and embed physically recurrent yarn networks in a meso-to-macro model for better generalisation under cyclic loading. We keep the constitutive law and learn only what the physics leaves open.
  • Calibration. Shear, bending and friction parameters are identified against your test curves, scatter included.

AI does not repair a wrong yarn law, and a surrogate fitted to some weaves says nothing about another. Ours state their range and send other queries to the full cell. See our surrogate modelling, hybrid physics-AI, uncertainty quantification and generative design pages and articles on surrogate models and validating them.

From one unit cell to a workflow your team keeps

The return comes when your engineers run the next weave, resin or tool through the same chain.

From one simulation study to a permanent workflow A loop of six steps: physics model, validation, automated runs, AI surrogate, decision tool and field data, which feeds back into the physics model. YOUR WORKFLOW Better with every run physics + data + AI, kept in-house 1Physics model 2Validation 3Automated runs 4AI surrogate 5Decision tool 6Field data CFD · FEA TESTS · TRIALS DOE · HPC CLOUD OPTIMISE DECIDE MONITOR RECALIBRATE
  1. Physics model. CFD and FEA set up once, as templates: geometry clean-up, mesh, solver settings and post-processing that every later run reuses.
  2. Validation. Mesh and time-step studies, and comparison with tests or full-scale data, so the model's uncertainty is known before it is trusted.
  3. Automated runs. Scripted runs over the parameters and operating conditions that matter, on your hardware or in the cloud.
  4. AI surrogate. A machine-learning model trained on those runs answers in seconds and says how sure it is; low-confidence questions go back to the solver.
  5. Decision tool. Optimisation, trade-off charts or operating guidance, in a form your engineers and operators use every day.
  6. Field data. Test, trial and monitoring data flow back to recalibrate the model, so it improves with every project instead of being rebuilt.
Physics, validation, automation and AI in one loop that improves with every run. Illustrative schematic.
  • A parametric unit-cell template: pattern, yarn and spacing as parameters, six strain cases scripted, checked on a case with a known answer.
  • A scan-to-cell hand-over if you hold CT data, with rejected cells recorded.
  • A draping template of shear, bending and friction data, each with its test, temperature and scatter.
  • A mapping tool with a contract: warp and weft vectors, shear and thickness per element in a versioned, tested file.
  • A validation record: coupons, a draped trial panel measured optically or by CT, tolerances set before the run.
  • A surrogate with a stated range and a gate to the full cell.
  • Handover and training, so your team runs and extends it: our training and software development services.

How credible is a woven composite simulation?

Verification is convergence: mesh or voxel size and load increment in the cell, a periodic-boundary check on a known case, the draping mesh and energy balance. Validation uses data you hold or can generate: coupons, a draped panel measured for shear angle, one case held back. Uncertainty is sampled: yarn and matrix data, cell-to-cell geometry, shear and friction data. Test scatter can be large: Boisse and colleagues report a 26-laboratory compressibility benchmark of woven and non-crimp glass fabrics with coefficients of variation of maximum pressure up to 50%. The answer is a range with ranked contributors.

Questions engineers ask about woven composite simulation

What is a unit cell in woven composite simulation?
One repeat of the weave, meshed and loaded with periodic boundary conditions so it stands for the whole fabric. Its averaged response gives the stiffness the structural model uses.
Can I use the shear locking angle as a wrinkling limit?
Only as a flag. Wrinkling depends on all the stresses, and angles above 60 degrees have been measured without wrinkles at high binder pressure.
Does a 3D weave remove delamination?
It changes it. In NASA's curved-beam test the 3D weaves failed by radial cracking instead, so architecture and load case decide.
Can simulation replace coupon or draping trials?
No. It narrows what trials must answer. Sign-off stays with you, the certifier and your customer.

If a weave choice, preform trial or allowables plan is waiting on this, send us the weave style, yarn and resin data, the coupon and picture-frame or bias-extension curves you hold, a draped trial's measurements and the part geometry with layup. We reply with the unit-cell and draping chain we would set up, the validation we would run and what a surrogate could cover. Start with a conversation, or see our woven fabric composites page.

OLTRION Engineering

CAE, CFD and simulation team

The OLTRION engineering team runs CFD, FEA and multiphysics simulation programmes for aerospace, automotive, energy and industrial clients, and builds the machine-learning tools that speed them up. These articles come from that project work.

Where this applies

Related capabilities

Have a model that will not converge?

Talk to the engineers who wrote this

Send the geometry, the physics and the decision it has to support. We answer with an approach, not a sales call.

Talk to an engineer