@theerapat-s28/ts-structural-eng-tools 2.3.0 → 3.0.0

package/dist/fea/frame/solver/frame-system.d.ts

@@ -1,14 +1,20 @@
 import type { SectionProperties } from "../../../core/types/section-properties.type";
 import type { Matrix } from "../../../utils/matrix";
 import { type SkylineMatrix } from "../../../utils/skyline";
+import { type FlatShellElement } from "../elements/shell/flat-shell";
 import { type CondensedMember } from "../elements/static-condensation";
-import type { FrameId, FrameModel, FrameNode } from "../types/model.type";
+import type { FrameId, FrameModel, FrameNode, ShellModifiers, ShellSection, ShellType } from "../types/model.type";
 import { type AutoRestraint } from "./auto-restraint";
 import { type BoundaryConditions } from "./boundary-conditions";
 import { type AssemblyMember } from "./global-assembly";
-/** A member with everything that does not depend on the loads. */
+/**
+ * A member, or a piece of one, with everything that does not depend on the loads. A member is
+ * cut into pieces where generated shell-mesh nodes lie inside it; otherwise it is one piece.
+ */
 export interface SystemMember extends AssemblyMember {
     id: FrameId;
+    /** Distance of the piece's start from the member's start (mm); 0 for a whole member. */
+    offset: number;
     section: SectionProperties;
     length: number;
     lambda: Matrix;
@@ -19,12 +25,43 @@ export interface SystemMember extends AssemblyMember {
     /** The member's own end displacements from its end nodes' (differ across releases). */
     memberEndDisplacements: CondensedMember["memberEndDisplacements"];
 }
+/** A model member and the pieces it is analysed as (one unless the shell mesh splits it). */
+export interface MemberGroup {
+    id: FrameId;
+    /** Indices into `FrameSystem.members`, from the member's start to its end. */
+    pieces: number[];
+    section: SectionProperties;
+    length: number;
+    lambda: Matrix;
+}
+/** A shell area with everything that does not depend on the loads. */
+export interface SystemShell {
+    id: FrameId;
+    /** The shell area the element belongs to (its own id when the area is one element). */
+    areaId: FrameId;
+    nodeIndices: number[];
+    type: ShellType;
+    section: ShellSection;
+    modifiers: ShellModifiers;
+    /** Axes, corner transformations and local and global stiffness. */
+    element: FlatShellElement;
+    /** ∫ Nᵢ dA of each corner on the mean plane (mm²): its share of a uniform surface load. */
+    cornerAreas: number[];
+}
 /** The load-independent part of a frame analysis: geometry, stiffness and its factorisation. */
 export interface FrameSystem {
     nodes: FrameNode[];
     nodeIndex: Map<FrameId, number>;
+    /** Index into `memberGroups` of each model member. */
     memberIndex: Map<FrameId, number>;
+    /** Model members, in model order. */
+    memberGroups: MemberGroup[];
+    /** The pieces that are assembled: one per member unless the shell mesh splits it. */
     members: SystemMember[];
+    shellIndex: Map<FrameId, number>;
+    /** Element indices of each shell area. */
+    shellsByArea: Map<FrameId, number[]>;
+    shells: SystemShell[];
     dofCount: number;
     boundary: BoundaryConditions;
     /** Equation number of each global DOF in the factorised system, −1 for a fixed DOF. */
@@ -39,12 +76,14 @@ export interface FrameSystem {
     profileSize: number;
 }
 /**
- * Builds the member matrices of a model that passed `validateFrameModel`, assembles the free
- * stiffness in skyline storage with nodes in reverse Cuthill–McKee order, and factorises it
- * once for all load cases. Joint rotations that no member resists (e.g. at truss joints) are
- * held by a spring along each null direction, which carries no force.
+ * Builds the member and shell matrices of a model that passed `validateFrameModel`, assembles
+ * the free stiffness in skyline storage with nodes in reverse Cuthill–McKee order, and
+ * factorises it once for all load cases. Joint rotations that nothing resists (e.g. at truss
+ * joints) are held by a spring along each null direction, which carries no force.
  *
- * @throws {FrameAnalysisError} 304 if the structure is a mechanism; `details` names a node and
- *   a DOF involved in it
+ * @throws {FrameAnalysisError} 328 if a node connected only to shells has a translation that
+ *   nothing resists (a membrane's out-of-plane or a plate's in-plane direction), with `nodeId`,
+ *   `shellId` and `dof`; 304 if the structure is a mechanism; `details` names a node and a DOF
+ *   involved in it
  */
 export declare const prepareFrameSystem: (model: FrameModel) => FrameSystem;

package/dist/fea/frame/solver/frame-system.js

@@ -6,6 +6,8 @@ const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.erro
 const skyline_1 = require("../../../utils/skyline");
 const beam_3d_stiffness_1 = require("../elements/beam-3d-stiffness");
 const end_releases_1 = require("../elements/end-releases");
+const flat_shell_1 = require("../elements/shell/flat-shell");
+const quad_geometry_1 = require("../elements/shell/quad-geometry");
 const static_condensation_1 = require("../elements/static-condensation");
 const coordinate_transform_1 = require("../transformation/coordinate-transform");
 const validate_model_1 = require("../validation/validate-model");
@@ -14,46 +16,152 @@ const boundary_conditions_1 = require("./boundary-conditions");
 const global_assembly_1 = require("./global-assembly");
 const DOF_NAMES = ["ux", "uy", "uz", "rx", "ry", "rz"];
 /**
- * Builds the member matrices of a model that passed `validateFrameModel`, assembles the free
- * stiffness in skyline storage with nodes in reverse Cuthill–McKee order, and factorises it
- * once for all load cases. Joint rotations that no member resists (e.g. at truss joints) are
- * held by a spring along each null direction, which carries no force.
+ * Builds the member and shell matrices of a model that passed `validateFrameModel`, assembles
+ * the free stiffness in skyline storage with nodes in reverse Cuthill–McKee order, and
+ * factorises it once for all load cases. Joint rotations that nothing resists (e.g. at truss
+ * joints) are held by a spring along each null direction, which carries no force.
  *
- * @throws {FrameAnalysisError} 304 if the structure is a mechanism; `details` names a node and
- *   a DOF involved in it
+ * @throws {FrameAnalysisError} 328 if a node connected only to shells has a translation that
+ *   nothing resists (a membrane's out-of-plane or a plate's in-plane direction), with `nodeId`,
+ *   `shellId` and `dof`; 304 if the structure is a mechanism; `details` names a node and a DOF
+ *   involved in it
  */
 const prepareFrameSystem = (model) => {
     const { nodes, members, restraints } = model;
+    // shell elements: after `expandModel` they carry their area's id
+    const shells = (model.shells ?? []);
     const nodeIndex = new Map(nodes.map((n, i) => [n.id, i]));
     const memberIndex = new Map(members.map((m, i) => [m.id, i]));
+    const shellIndex = new Map(shells.map((sh, i) => [sh.id, i]));
     const dofCount = 6 * nodes.length;
     const boundary = (0, boundary_conditions_1.buildBoundaryConditions)(restraints, nodeIndex, dofCount);
-    const systemMembers = members.map((member) => {
-        const start = nodes[nodeIndex.get(member.startNodeId)];
-        const end = nodes[nodeIndex.get(member.endNodeId)];
-        const { length, lambda } = (0, coordinate_transform_1.memberAxes)([start.x, start.y, start.z], [end.x, end.y, end.z], member.angle ?? 0);
+    const point = (id) => {
+        const n = nodes[nodeIndex.get(id)];
+        return [n.x, n.y, n.z];
+    };
+    const systemMembers = [];
+    const memberGroups = members.map((member) => {
+        // pieces between the member's ends and the mesh nodes inside it (see `expandModel`)
+        const chain = [member.startNodeId, ...(member.through ?? []), member.endNodeId];
+        const whole = (0, coordinate_transform_1.memberAxes)(point(member.startNodeId), point(member.endNodeId), member.angle ?? 0);
         const { E, G, A, I22, I33, J } = member.section;
-        const releases = (0, end_releases_1.memberEndReleases)(member.releases);
-        const condensed = (0, static_condensation_1.condenseMember)((0, beam_3d_stiffness_1.beam3dLocalStiffnessRightHanded)(E, G, A, I22, I33, J, length), releases);
+        const pieces = [];
+        let offset = 0;
+        for (let k = 0; k + 1 < chain.length; k++) {
+            const { length, lambda } = chain.length === 2
+                ? whole
+                : (0, coordinate_transform_1.memberAxes)(point(chain[k]), point(chain[k + 1]), member.angle ?? 0);
+            // releases only at the member's own ends
+            const releases = (0, end_releases_1.memberEndReleases)({
+                start: k === 0 ? member.releases?.start : undefined,
+                end: k === chain.length - 2 ? member.releases?.end : undefined,
+            });
+            const condensed = (0, static_condensation_1.condenseMember)((0, beam_3d_stiffness_1.beam3dLocalStiffnessRightHanded)(E, G, A, I22, I33, J, length), releases);
+            pieces.push(systemMembers.length);
+            systemMembers.push({
+                id: member.id,
+                offset,
+                section: member.section,
+                startIndex: nodeIndex.get(chain[k]),
+                endIndex: nodeIndex.get(chain[k + 1]),
+                T: (0, coordinate_transform_1.transformationMatrix)(lambda),
+                k: condensed.k,
+                hasReleases: releases.length > 0,
+                condenseFixedEndForces: condensed.condenseFixedEndForces,
+                memberEndDisplacements: condensed.memberEndDisplacements,
+                length,
+                lambda,
+            });
+            offset += length;
+        }
         return {
             id: member.id,
+            pieces,
             section: member.section,
-            startIndex: nodeIndex.get(member.startNodeId),
-            endIndex: nodeIndex.get(member.endNodeId),
-            T: (0, coordinate_transform_1.transformationMatrix)(lambda),
-            k: condensed.k,
-            hasReleases: releases.length > 0,
-            condenseFixedEndForces: condensed.condenseFixedEndForces,
-            memberEndDisplacements: condensed.memberEndDisplacements,
-            length,
-            lambda,
+            length: whole.length,
+            lambda: whole.lambda,
+        };
+    });
+    const systemShells = shells.map((shell) => {
+        const nodeIndices = shell.nodeIds.map((id) => nodeIndex.get(id));
+        const points = nodeIndices.map((n) => [nodes[n].x, nodes[n].y, nodes[n].z]);
+        const { thickness, E, nu } = shell.section;
+        const m = shell.modifiers ?? {};
+        const type = shell.type ?? "shell";
+        const element = (0, flat_shell_1.flatShellElement)(points, {
+            E,
+            nu,
+            t: thickness,
+            type,
+            membraneModifiers: [m.f11 ?? 1, m.f22 ?? 1, m.f12 ?? 1],
+            bendingModifiers: [m.m11 ?? 1, m.m22 ?? 1, m.m12 ?? 1],
+        }, shell.angle ?? 0);
+        return {
+            id: shell.id,
+            areaId: shell.areaId ?? shell.id,
+            nodeIndices,
+            type,
+            section: shell.section,
+            modifiers: m,
+            element,
+            cornerAreas: (0, quad_geometry_1.cornerAreas)(element.axes.corners),
         };
     });
+    const shellsByArea = new Map();
+    systemShells.forEach((sh, s) => {
+        const list = shellsByArea.get(sh.areaId);
+        if (list)
+            list.push(s);
+        else
+            shellsByArea.set(sh.areaId, [s]);
+    });
+    // Everything that adds stiffness, as elements over 6 DOFs per node; members have 2 nodes
+    const memberElements = systemMembers.map((m) => ({
+        nodes: [m.startIndex, m.endIndex],
+        kGlobal: (0, global_assembly_1.memberGlobalStiffness)(m),
+    }));
+    const shellElements = systemShells.map((sh) => ({
+        nodes: sh.nodeIndices,
+        kGlobal: sh.element.kGlobal,
+    }));
+    const elements = [...memberElements, ...shellElements];
+    // ---- Translations no shell resists at nodes without members (membrane / plate only) --
+    const memberNodes = new Set(systemMembers.flatMap((m) => [m.startIndex, m.endIndex]));
+    const shellOnlyNodes = new Set(systemShells.flatMap((sh) => sh.nodeIndices).filter((n) => !memberNodes.has(n)));
+    const unresisted = (0, auto_restraint_1.findUnresistedTranslations)([...shellOnlyNodes].sort((a, b) => a - b), // issues in model order
+    shellElements, boundary);
+    if (unresisted.length > 0) {
+        (0, validate_model_1.throwIssues)(unresisted.map(({ node, axis }) => {
+            const lead = axis.reduce((best, c, i) => (Math.abs(c) > Math.abs(axis[best]) ? i : best), 0);
+            return (0, validate_model_1.frameIssue)(frame_analysis_error_1.FrameErrors.SHELL_DOF_UNSUPPORTED, {
+                nodeId: nodes[node].id,
+                shellId: systemShells.find((sh) => sh.nodeIndices.includes(node))?.id,
+                dof: DOF_NAMES[lead],
+            });
+        }));
+    }
     // ---- Equation numbering: nodes in RCM order, fixed DOFs skipped ------------------
     const adjacency = nodes.map(() => []);
-    for (const m of systemMembers) {
-        adjacency[m.startIndex].push(m.endIndex);
-        adjacency[m.endIndex].push(m.startIndex);
+    for (const element of memberElements) {
+        for (const a of element.nodes) {
+            for (const b of element.nodes)
+                if (a !== b)
+                    adjacency[a].push(b);
+        }
+    }
+    // shells share edges: add each neighbour once (members keep their own order unchanged)
+    if (shellElements.length > 0) {
+        const known = adjacency.map((list) => new Set(list));
+        for (const element of shellElements) {
+            for (const a of element.nodes) {
+                for (const b of element.nodes) {
+                    if (a !== b && !known[a].has(b)) {
+                        known[a].add(b);
+                        adjacency[a].push(b);
+                    }
+                }
+            }
+        }
     }
     const equation = new Int32Array(dofCount).fill(-1);
     const free = [];
@@ -66,18 +174,13 @@ const prepareFrameSystem = (model) => {
             }
         }
     }
-    const kGlobal = systemMembers.map(global_assembly_1.memberGlobalStiffness);
-    const autoRestraints = (0, auto_restraint_1.findAutoRestraints)(nodes.length, systemMembers.map((m, i) => ({
-        startIndex: m.startIndex,
-        endIndex: m.endIndex,
-        kGlobal: kGlobal[i],
-    })), boundary);
+    const autoRestraints = (0, auto_restraint_1.findAutoRestraints)(nodes.length, elements, boundary);
     // ---- Profile, assembly and factorisation -----------------------------------------
     let factor = null;
     if (free.length > 0) {
         const first = free.map((_, e) => e);
-        const memberEquations = systemMembers.map((m) => (0, global_assembly_1.memberDofs)(m).map((d) => equation[d]));
-        for (const eqs of memberEquations) {
+        const elementEquations = elements.map((el) => (0, global_assembly_1.elementDofs)(el.nodes).map((d) => equation[d]));
+        for (const eqs of elementEquations) {
             const active = eqs.filter((e) => e >= 0);
             const low = Math.min(...active);
             for (const e of active)
@@ -85,14 +188,15 @@ const prepareFrameSystem = (model) => {
                     first[e] = low;
         }
         const K = (0, skyline_1.createSkyline)(first);
-        memberEquations.forEach((eqs, index) => {
-            // a member's 12 DOFs are distinct, so each unordered pair and each diagonal is added once
-            for (let a = 0; a < 12; a++) {
+        elementEquations.forEach((eqs, index) => {
+            // an element's nodes, hence its DOFs, are distinct: each pair and diagonal is added once
+            const k = elements[index].kGlobal;
+            for (let a = 0; a < eqs.length; a++) {
                 if (eqs[a] < 0)
                     continue;
-                for (let b = 0; b < 12; b++) {
+                for (let b = 0; b < eqs.length; b++) {
                     if (eqs[b] >= eqs[a])
-                        (0, skyline_1.skylineAdd)(K, eqs[a], eqs[b], kGlobal[index][a][b]);
+                        (0, skyline_1.skylineAdd)(K, eqs[a], eqs[b], k[a][b]);
                 }
             }
         });
@@ -100,7 +204,7 @@ const prepareFrameSystem = (model) => {
             if (boundary.spring[dof] > 0)
                 (0, skyline_1.skylineAdd)(K, e, e, boundary.spring[dof]);
         });
-        // k·n·nᵀ on the node's free rotations, inside the profile: its members span all 6 node DOFs
+        // k·n·nᵀ on the node's free rotations, inside the profile: its elements span all 6 node DOFs
         for (const { node, axis, stiffness } of autoRestraints) {
             for (let i = 0; i < 3; i++) {
                 const ei = equation[6 * node + 3 + i];
@@ -129,7 +233,11 @@ const prepareFrameSystem = (model) => {
         nodes,
         nodeIndex,
         memberIndex,
+        memberGroups,
         members: systemMembers,
+        shellIndex,
+        shellsByArea,
+        shells: systemShells,
         dofCount,
         boundary,
         equation,

package/dist/fea/frame/solver/frame-system.js.map

@@ -1 +1 @@
-{"version":3,"file":"frame-system.js","sourceRoot":"","sources":["../../../../src/fea/frame/solver/frame-system.ts"],"names":[],"mappings":";;;AAAA,uEAA2E;AAC3E,gFAAoE;AAGpE,gDAO4B;AAE5B,qEAAgF;AAChF,2DAA6D;AAC7D,yEAAuF;AACvF,iFAA0F;AAE1F,iEAAuE;AACvE,qDAA0E;AAC1E,+DAAyF;AACzF,uDAA2F;AAE3F,MAAM,SAAS,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AAoChE;;;;;;;;GAQG;AACI,MAAM,kBAAkB,GAAG,CAAC,KAAiB,EAAe,EAAE;IACnE,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,UAAU,EAAE,GAAG,KAAK,CAAC;IAC7C,MAAM,SAAS,GAAG,IAAI,GAAG,CAAkB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3E,MAAM,WAAW,GAAG,IAAI,GAAG,CAAkB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,QAAQ,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IAClC,MAAM,QAAQ,GAAG,IAAA,6CAAuB,EAAC,UAAU,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;IAE1E,MAAM,aAAa,GAAmB,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QAC3D,MAAM,KAAK,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,WAAW,CAAW,CAAC,CAAC;QACjE,MAAM,GAAG,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAW,CAAC,CAAC;QAC7D,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,IAAA,iCAAU,EACnC,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAC3B,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EACrB,MAAM,CAAC,KAAK,IAAI,CAAC,CAClB,CAAC;QACF,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC;QAChD,MAAM,QAAQ,GAAG,IAAA,gCAAiB,EAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;QACpD,MAAM,SAAS,GAAG,IAAA,oCAAc,EAC9B,IAAA,mDAA+B,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,CAAC,EAC7D,QAAQ,CACT,CAAC;QACF,OAAO;YACL,EAAE,EAAE,MAAM,CAAC,EAAE;YACb,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,UAAU,EAAE,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,WAAW,CAAW;YACvD,QAAQ,EAAE,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAW;YACnD,CAAC,EAAE,IAAA,2CAAoB,EAAC,MAAM,CAAC;YAC/B,CAAC,EAAE,SAAS,CAAC,CAAC;YACd,WAAW,EAAE,QAAQ,CAAC,MAAM,GAAG,CAAC;YAChC,sBAAsB,EAAE,SAAS,CAAC,sBAAsB;YACxD,sBAAsB,EAAE,SAAS,CAAC,sBAAsB;YACxD,MAAM;YACN,MAAM;SACP,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,qFAAqF;IACrF,MAAM,SAAS,GAAe,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAClD,KAAK,MAAM,CAAC,IAAI,aAAa,EAAE,CAAC;QAC9B,SAAS,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;QACzC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC;IAC3C,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,UAAU,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,IAAA,6BAAmB,EAAC,SAAS,CAAC,EAAE,CAAC;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YACzB,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;gBACzB,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;gBAC5B,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACjB,CAAC;QACH,CAAC;IACH,CAAC;IAED,MAAM,OAAO,GAAG,aAAa,CAAC,GAAG,CAAC,uCAAqB,CAAC,CAAC;IACzD,MAAM,cAAc,GAAG,IAAA,mCAAkB,EACvC,KAAK,CAAC,MAAM,EACZ,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC3B,UAAU,EAAE,CAAC,CAAC,UAAU;QACxB,QAAQ,EAAE,CAAC,CAAC,QAAQ;QACpB,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;KACpB,CAAC,CAAC,EACH,QAAQ,CACT,CAAC;IAEF,qFAAqF;IACrF,IAAI,MAAM,GAAyB,IAAI,CAAC;IACxC,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACpC,MAAM,eAAe,GAAG,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAA,4BAAU,EAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxF,KAAK,MAAM,GAAG,IAAI,eAAe,EAAE,CAAC;YAClC,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;YAChC,KAAK,MAAM,CAAC,IAAI,MAAM;gBAAE,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBAAE,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAC7D,CAAC;QAED,MAAM,CAAC,GAAG,IAAA,uBAAa,EAAC,KAAK,CAAC,CAAC;QAC/B,eAAe,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE;YACrC,0FAA0F;YAC1F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC5B,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,SAAS;gBACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC5B,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;wBAAE,IAAA,oBAAU,EAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5E,CAAC;YACH,CAAC;QACH,CAAC,CAAC,CAAC;QACH,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE;YACtB,IAAI,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC;gBAAE,IAAA,oBAAU,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1E,CAAC,CAAC,CAAC;QACH,4FAA4F;QAC5F,KAAK,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,cAAc,EAAE,CAAC;YACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBACtC,IAAI,EAAE,GAAG,CAAC;oBAAE,SAAS;gBACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;oBACtC,IAAI,EAAE,IAAI,EAAE;wBAAE,IAAA,oBAAU,EAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrE,CAAC;YACH,CAAC;QACH,CAAC;QAED,MAAM,MAAM,GAAG,IAAA,uBAAa,EAAC,CAAC,EAAE,sCAAqB,CAAC,CAAC;QACvD,IAAI,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC;YACf,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC;YACrC,IAAA,4BAAW,EAAC;gBACV,IAAA,2BAAU,EAAC,kCAAW,CAAC,kBAAkB,EAAE;oBACzC,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE;oBACrC,GAAG,EAAE,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC;iBACxB,CAAC;aACH,CAAC,CAAC;QACL,CAAC;QACD,MAAM,GAAG,CAAC,CAAC;IACb,CAAC;IAED,OAAO;QACL,KAAK;QACL,SAAS;QACT,WAAW;QACX,OAAO,EAAE,aAAa;QACtB,QAAQ;QACR,QAAQ;QACR,QAAQ;QACR,IAAI;QACJ,cAAc;QACd,MAAM;QACN,WAAW,EAAE,MAAM,CAAC,CAAC,CAAC,IAAA,qBAAW,EAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;KAC9C,CAAC;AACJ,CAAC,CAAC;AA/HW,QAAA,kBAAkB,sBA+H7B"}
+{"version":3,"file":"frame-system.js","sourceRoot":"","sources":["../../../../src/fea/frame/solver/frame-system.ts"],"names":[],"mappings":";;;AAAA,uEAA2E;AAC3E,gFAAoE;AAGpE,gDAO4B;AAE5B,qEAAgF;AAChF,2DAA6D;AAC7D,6DAAuF;AACvF,mEAA8D;AAE9D,yEAAuF;AACvF,iFAA0F;AAW1F,iEAAuE;AACvE,qDAI0B;AAC1B,+DAAyF;AACzF,uDAK2B;AAE3B,MAAM,SAAS,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AA0EhE;;;;;;;;;;GAUG;AACI,MAAM,kBAAkB,GAAG,CAAC,KAAiB,EAAe,EAAE;IACnE,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,UAAU,EAAE,GAAG,KAAK,CAAC;IAC7C,iEAAiE;IACjE,MAAM,MAAM,GAAG,CAAC,KAAK,CAAC,MAAM,IAAI,EAAE,CAAyC,CAAC;IAC5E,MAAM,SAAS,GAAG,IAAI,GAAG,CAAkB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3E,MAAM,WAAW,GAAG,IAAI,GAAG,CAAkB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,UAAU,GAAG,IAAI,GAAG,CAAkB,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,QAAQ,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IAClC,MAAM,QAAQ,GAAG,IAAA,6CAAuB,EAAC,UAAU,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;IAE1E,MAAM,KAAK,GAAG,CAAC,EAAW,EAA4B,EAAE;QACtD,MAAM,CAAC,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,EAAE,CAAW,CAAC,CAAC;QAC7C,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC;IACF,MAAM,aAAa,GAAmB,EAAE,CAAC;IACzC,MAAM,YAAY,GAAmB,OAAgD,CAAC,GAAG,CACvF,CAAC,MAAM,EAAE,EAAE;QACT,oFAAoF;QACpF,MAAM,KAAK,GAAG,CAAC,MAAM,CAAC,WAAW,EAAE,GAAG,CAAC,MAAM,CAAC,OAAO,IAAI,EAAE,CAAC,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC;QAChF,MAAM,KAAK,GAAG,IAAA,iCAAU,EACtB,KAAK,CAAC,MAAM,CAAC,WAAW,CAAC,EACzB,KAAK,CAAC,MAAM,CAAC,SAAS,CAAC,EACvB,MAAM,CAAC,KAAK,IAAI,CAAC,CAClB,CAAC;QACF,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC;QAChD,MAAM,MAAM,GAAa,EAAE,CAAC;QAC5B,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GACtB,KAAK,CAAC,MAAM,KAAK,CAAC;gBAChB,CAAC,CAAC,KAAK;gBACP,CAAC,CAAC,IAAA,iCAAU,EAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC;YAC1E,yCAAyC;YACzC,MAAM,QAAQ,GAAG,IAAA,gCAAiB,EAAC;gBACjC,KAAK,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC,SAAS;gBACnD,GAAG,EAAE,CAAC,KAAK,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC,CAAC,SAAS;aAC/D,CAAC,CAAC;YACH,MAAM,SAAS,GAAG,IAAA,oCAAc,EAC9B,IAAA,mDAA+B,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,CAAC,EAC7D,QAAQ,CACT,CAAC;YACF,MAAM,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,CAAC;YAClC,aAAa,CAAC,IAAI,CAAC;gBACjB,EAAE,EAAE,MAAM,CAAC,EAAE;gBACb,MAAM;gBACN,OAAO,EAAE,MAAM,CAAC,OAAO;gBACvB,UAAU,EAAE,SAAS,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAW;gBAC7C,QAAQ,EAAE,SAAS,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAW;gBAC/C,CAAC,EAAE,IAAA,2CAAoB,EAAC,MAAM,CAAC;gBAC/B,CAAC,EAAE,SAAS,CAAC,CAAC;gBACd,WAAW,EAAE,QAAQ,CAAC,MAAM,GAAG,CAAC;gBAChC,sBAAsB,EAAE,SAAS,CAAC,sBAAsB;gBACxD,sBAAsB,EAAE,SAAS,CAAC,sBAAsB;gBACxD,MAAM;gBACN,MAAM;aACP,CAAC,CAAC;YACH,MAAM,IAAI,MAAM,CAAC;QACnB,CAAC;QACD,OAAO;YACL,EAAE,EAAE,MAAM,CAAC,EAAE;YACb,MAAM;YACN,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,MAAM,EAAE,KAAK,CAAC,MAAM;YACpB,MAAM,EAAE,KAAK,CAAC,MAAM;SACrB,CAAC;IACJ,CAAC,CACF,CAAC;IAEF,MAAM,YAAY,GAAkB,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE;QACvD,MAAM,WAAW,GAAG,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,SAAS,CAAC,GAAG,CAAC,EAAE,CAAW,CAAC,CAAC;QAC3E,MAAM,MAAM,GAAG,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,EAAQ,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClF,MAAM,EAAE,SAAS,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,KAAK,CAAC,OAAO,CAAC;QAC3C,MAAM,CAAC,GAAG,KAAK,CAAC,SAAS,IAAI,EAAE,CAAC;QAChC,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,IAAI,OAAO,CAAC;QACnC,MAAM,OAAO,GAAG,IAAA,6BAAgB,EAC9B,MAAM,EACN;YACE,CAAC;YACD,EAAE;YACF,CAAC,EAAE,SAAS;YACZ,IAAI;YACJ,iBAAiB,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;YACvD,gBAAgB,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;SACvD,EACD,KAAK,CAAC,KAAK,IAAI,CAAC,CACjB,CAAC;QACF,OAAO;YACL,EAAE,EAAE,KAAK,CAAC,EAAE;YACZ,MAAM,EAAE,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,EAAE;YAChC,WAAW;YACX,IAAI;YACJ,OAAO,EAAE,KAAK,CAAC,OAAO;YACtB,SAAS,EAAE,CAAC;YACZ,OAAO;YACP,WAAW,EAAE,IAAA,2BAAW,EAAC,OAAO,CAAC,IAAI,CAAC,OAAO,CAAC;SAC/C,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,MAAM,YAAY,GAAG,IAAI,GAAG,EAAqB,CAAC;IAClD,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE;QAC7B,MAAM,IAAI,GAAG,YAAY,CAAC,GAAG,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;QACzC,IAAI,IAAI;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;;YAClB,YAAY,CAAC,GAAG,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACxC,CAAC,CAAC,CAAC;IAEH,yFAAyF;IACzF,MAAM,cAAc,GAAuB,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QACnE,KAAK,EAAE,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,CAAC,QAAQ,CAAC;QACjC,OAAO,EAAE,IAAA,uCAAqB,EAAC,CAAC,CAAC;KAClC,CAAC,CAAC,CAAC;IACJ,MAAM,aAAa,GAAuB,YAAY,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAClE,KAAK,EAAE,EAAE,CAAC,WAAW;QACrB,OAAO,EAAE,EAAE,CAAC,OAAO,CAAC,OAAO;KAC5B,CAAC,CAAC,CAAC;IACJ,MAAM,QAAQ,GAAG,CAAC,GAAG,cAAc,EAAE,GAAG,aAAa,CAAC,CAAC;IAEvD,yFAAyF;IACzF,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IACtF,MAAM,cAAc,GAAG,IAAI,GAAG,CAC5B,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,WAAW,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAChF,CAAC;IACF,MAAM,UAAU,GAAG,IAAA,2CAA0B,EAC3C,CAAC,GAAG,cAAc,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,wBAAwB;IACnE,aAAa,EACb,QAAQ,CACT,CAAC;IACF,IAAI,UAAU,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC1B,IAAA,4BAAW,EACT,UAAU,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE;YAChC,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CACtB,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAC/D,CAAC,CACF,CAAC;YACF,OAAO,IAAA,2BAAU,EAAC,kCAAW,CAAC,qBAAqB,EAAE;gBACnD,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE;gBACtB,OAAO,EAAE,YAAY,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,WAAW,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE;gBACrE,GAAG,EAAE,SAAS,CAAC,IAAI,CAAC;aACrB,CAAC,CAAC;QACL,CAAC,CAAC,CACH,CAAC;IACJ,CAAC;IAED,qFAAqF;IACrF,MAAM,SAAS,GAAe,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAClD,KAAK,MAAM,OAAO,IAAI,cAAc,EAAE,CAAC;QACrC,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK,EAAE,CAAC;YAC9B,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK;gBAAE,IAAI,CAAC,KAAK,CAAC;oBAAE,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnE,CAAC;IACH,CAAC;IACD,uFAAuF;IACvF,IAAI,aAAa,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC7B,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC;QACrD,KAAK,MAAM,OAAO,IAAI,aAAa,EAAE,CAAC;YACpC,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK,EAAE,CAAC;gBAC9B,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK,EAAE,CAAC;oBAC9B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC;wBAChC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBAChB,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;oBACvB,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,UAAU,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,IAAA,6BAAmB,EAAC,SAAS,CAAC,EAAE,CAAC;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YACzB,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;gBACzB,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;gBAC5B,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACjB,CAAC;QACH,CAAC;IACH,CAAC;IAED,MAAM,cAAc,GAAG,IAAA,mCAAkB,EAAC,KAAK,CAAC,MAAM,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;IAE5E,qFAAqF;IACrF,IAAI,MAAM,GAAyB,IAAI,CAAC;IACxC,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACpC,MAAM,gBAAgB,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,IAAA,6BAAW,EAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7F,KAAK,MAAM,GAAG,IAAI,gBAAgB,EAAE,CAAC;YACnC,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;YAChC,KAAK,MAAM,CAAC,IAAI,MAAM;gBAAE,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBAAE,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAC7D,CAAC;QAED,MAAM,CAAC,GAAG,IAAA,uBAAa,EAAC,KAAK,CAAC,CAAC;QAC/B,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE;YACtC,yFAAyF;YACzF,MAAM,CAAC,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC,OAAO,CAAC;YAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACpC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,SAAS;gBACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBACpC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;wBAAE,IAAA,oBAAU,EAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/D,CAAC;YACH,CAAC;QACH,CAAC,CAAC,CAAC;QACH,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE;YACtB,IAAI,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC;gBAAE,IAAA,oBAAU,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1E,CAAC,CAAC,CAAC;QACH,6FAA6F;QAC7F,KAAK,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,cAAc,EAAE,CAAC;YACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBACtC,IAAI,EAAE,GAAG,CAAC;oBAAE,SAAS;gBACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;oBACtC,IAAI,EAAE,IAAI,EAAE;wBAAE,IAAA,oBAAU,EAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrE,CAAC;YACH,CAAC;QACH,CAAC;QAED,MAAM,MAAM,GAAG,IAAA,uBAAa,EAAC,CAAC,EAAE,sCAAqB,CAAC,CAAC;QACvD,IAAI,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC;YACf,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC;YACrC,IAAA,4BAAW,EAAC;gBACV,IAAA,2BAAU,EAAC,kCAAW,CAAC,kBAAkB,EAAE;oBACzC,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE;oBACrC,GAAG,EAAE,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC;iBACxB,CAAC;aACH,CAAC,CAAC;QACL,CAAC;QACD,MAAM,GAAG,CAAC,CAAC;IACb,CAAC;IAED,OAAO;QACL,KAAK;QACL,SAAS;QACT,WAAW;QACX,YAAY;QACZ,OAAO,EAAE,aAAa;QACtB,UAAU;QACV,YAAY;QACZ,MAAM,EAAE,YAAY;QACpB,QAAQ;QACR,QAAQ;QACR,QAAQ;QACR,IAAI;QACJ,cAAc;QACd,MAAM;QACN,WAAW,EAAE,MAAM,CAAC,CAAC,CAAC,IAAA,qBAAW,EAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;KAC9C,CAAC;AACJ,CAAC,CAAC;AApPW,QAAA,kBAAkB,sBAoP7B"}

package/dist/fea/frame/solver/global-assembly.d.ts

@@ -1,10 +1,21 @@
 import { type Matrix } from "../../../utils/matrix";
+/**
+ * A stiffness contribution to the global system: the nodes it connects and its stiffness over
+ * their DOFs (6 per node, in node order) in global axes. A member is a 2-node element; elements
+ * with more nodes assemble the same way.
+ */
+export interface StiffnessElement {
+    nodes: number[];
+    kGlobal: Matrix;
+}
 export interface AssemblyMember {
     startIndex: number;
     endIndex: number;
     T: Matrix;
     k: Matrix;
 }
+/** Global DOF numbers (6 per node, in node order) of an element connecting `nodes`. */
+export declare const elementDofs: (nodes: readonly number[]) => number[];
 /** Global DOF numbers (6 per node) of a member's 12 end DOFs. */
 export declare const memberDofs: (member: Pick<AssemblyMember, "startIndex" | "endIndex">) => number[];
 /** A member's stiffness in global axes, `Tᵀ·k·T` (N and mm). */

package/dist/fea/frame/solver/global-assembly.js

@@ -1,16 +1,19 @@
 "use strict";
 Object.defineProperty(exports, "__esModule", { value: true });
-exports.scatterMemberForces = exports.memberGlobalStiffness = exports.memberDofs = void 0;
+exports.scatterMemberForces = exports.memberGlobalStiffness = exports.memberDofs = exports.elementDofs = void 0;
 const matrix_1 = require("../../../utils/matrix");
-/** Global DOF numbers (6 per node) of a member's 12 end DOFs. */
-const memberDofs = (member) => {
+/** Global DOF numbers (6 per node, in node order) of an element connecting `nodes`. */
+const elementDofs = (nodes) => {
     const dofs = [];
-    for (const node of [member.startIndex, member.endIndex]) {
+    for (const node of nodes) {
         for (let k = 0; k < 6; k++)
             dofs.push(6 * node + k);
     }
     return dofs;
 };
+exports.elementDofs = elementDofs;
+/** Global DOF numbers (6 per node) of a member's 12 end DOFs. */
+const memberDofs = (member) => (0, exports.elementDofs)([member.startIndex, member.endIndex]);
 exports.memberDofs = memberDofs;
 /** A member's stiffness in global axes, `Tᵀ·k·T` (N and mm). */
 const memberGlobalStiffness = (member) => (0, matrix_1.matMul)((0, matrix_1.transpose)(member.T), (0, matrix_1.matMul)(member.k, member.T));

package/dist/fea/frame/solver/global-assembly.js.map

@@ -1 +1 @@
-{"version":3,"file":"global-assembly.js","sourceRoot":"","sources":["../../../../src/fea/frame/solver/global-assembly.ts"],"names":[],"mappings":";;;AAAA,8CAAmE;AASnE,iEAAiE;AAC1D,MAAM,UAAU,GAAG,CAAC,MAAuD,EAAY,EAAE;IAC9F,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,CAAC,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;IACtD,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC,CAAC;AANW,QAAA,UAAU,cAMrB;AAEF,gEAAgE;AACzD,MAAM,qBAAqB,GAAG,CAAC,MAAsB,EAAU,EAAE,CACtE,IAAA,eAAM,EAAC,IAAA,kBAAS,EAAC,MAAM,CAAC,CAAC,CAAC,EAAE,IAAA,eAAM,EAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;AAD7C,QAAA,qBAAqB,yBACwB;AAE1D;;;;GAIG;AACI,MAAM,mBAAmB,GAAG,CACjC,MAAgB,EAChB,MAA6D,EAC7D,WAAqB,EACrB,IAAY,EACN,EAAE;IACR,MAAM,IAAI,GAAG,IAAA,kBAAU,EAAC,MAAM,CAAC,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;QACtE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,KAAK,CAAC;IAClC,CAAC;AACH,CAAC,CAAC;AAZW,QAAA,mBAAmB,uBAY9B"}
+{"version":3,"file":"global-assembly.js","sourceRoot":"","sources":["../../../../src/fea/frame/solver/global-assembly.ts"],"names":[],"mappings":";;;AAAA,8CAAmE;AAmBnE,uFAAuF;AAChF,MAAM,WAAW,GAAG,CAAC,KAAwB,EAAY,EAAE;IAChE,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;IACtD,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC,CAAC;AANW,QAAA,WAAW,eAMtB;AAEF,iEAAiE;AAC1D,MAAM,UAAU,GAAG,CAAC,MAAuD,EAAY,EAAE,CAC9F,IAAA,mBAAW,EAAC,CAAC,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;AADvC,QAAA,UAAU,cAC6B;AAEpD,gEAAgE;AACzD,MAAM,qBAAqB,GAAG,CAAC,MAAsB,EAAU,EAAE,CACtE,IAAA,eAAM,EAAC,IAAA,kBAAS,EAAC,MAAM,CAAC,CAAC,CAAC,EAAE,IAAA,eAAM,EAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;AAD7C,QAAA,qBAAqB,yBACwB;AAE1D;;;;GAIG;AACI,MAAM,mBAAmB,GAAG,CACjC,MAAgB,EAChB,MAA6D,EAC7D,WAAqB,EACrB,IAAY,EACN,EAAE;IACR,MAAM,IAAI,GAAG,IAAA,kBAAU,EAAC,MAAM,CAAC,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;QACtE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,KAAK,CAAC;IAClC,CAAC;AACH,CAAC,CAAC;AAZW,QAAA,mBAAmB,uBAY9B"}

package/dist/fea/frame/solver/load-case-solver.d.ts

@@ -1,5 +1,5 @@
 import { type LocalMemberLoad } from "../elements/fixed-end-forces";
-import type { LoadCase, MemberLoad, NodalLoad } from "../types/model.type";
+import type { AreaLoad, LoadCase, MemberLoad, NodalLoad } from "../types/model.type";
 import type { FrameSystem } from "./frame-system";
 /**
  * Unrounded result of one load case or combination, in N and mm. Every field is linear in the
@@ -10,6 +10,7 @@ export interface RawState {
     R: number[];
     F: number[];
     memberForces: number[][];
+    shellForces: number[][];
     memberLoads: LocalMemberLoad[][];
 }
 /**
@@ -20,6 +21,14 @@ export interface RawState {
  * @param loadCase - The load case (id and self-weight multiplier)
  * @param loads - The nodal and member loads tagged with this load case
  *
+ * Shell areas take `selfWeight · unitWeight · thickness · modifiers.weight` per unit area along
+ * −Z and their `AreaLoad`s as forces at their corners, each corner the share `∫ Nᵢ dA` of the
+ * element's area on its mean plane (no corner moments, so a warped element's offsets do not load
+ * the rotations auto-restraint holds).
+ *
+ * @param loads - The nodal and member loads tagged with this load case, and its area loads on
+ *   shell areas
+ *
  * @throws {FrameAnalysisError} 304 if a load has a moment about a joint rotation that nothing
  *   resists (an auto-restrained axis), 308 for a load outside its member, 313 if self-weight is
  *   requested for a member without `unitWeight`
@@ -27,6 +36,7 @@ export interface RawState {
 export declare const solveLoadCase: (system: FrameSystem, loadCase: LoadCase, loads: {
     nodalLoads: NodalLoad[];
     memberLoads: MemberLoad[];
+    shellLoads?: AreaLoad[];
 }) => RawState;
 /** Superposes factored states: `Σ factor × state`. Member loads are concatenated, scaled. */
 export declare const combineStates: (terms: {

package/dist/fea/frame/solver/load-case-solver.js

@@ -6,6 +6,7 @@ const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.erro
 const matrix_1 = require("../../../utils/matrix");
 const skyline_1 = require("../../../utils/skyline");
 const fixed_end_forces_1 = require("../elements/fixed-end-forces");
+const shell_forces_1 = require("../post/shell-forces");
 const validate_model_1 = require("../validation/validate-model");
 const global_assembly_1 = require("./global-assembly");
 const GLOBAL_DIRECTIONS = {
@@ -52,6 +53,41 @@ const toLocalLoads = (load, lambda, L) => {
         };
     });
 };
+/**
+ * The part of a member's local loads that falls on one of its pieces, `[offset, offset + L]`,
+ * in the piece's own positions. A distributed load or torque is cut at the piece ends; a point
+ * load goes to the piece it lies in (at an interior node, to the piece that starts there).
+ */
+const clipLoads = (loads, offset, L, last) => {
+    const end = offset + L;
+    const tolerance = 1e-9 * end;
+    const out = [];
+    for (const load of loads) {
+        if (load.type === "point") {
+            const inside = load.a >= offset - tolerance && (load.a < end - tolerance || last);
+            if (inside)
+                out.push({ ...load, a: Math.min(L, Math.max(0, load.a - offset)) });
+            continue;
+        }
+        const from = Math.max(load.start, offset);
+        const to = Math.min(load.end, end);
+        if (to - from <= tolerance)
+            continue;
+        if (load.type === "torque") {
+            out.push({ ...load, start: from - offset, end: to - offset });
+        }
+        else {
+            out.push({
+                ...load,
+                w1: (0, fixed_end_forces_1.loadIntensity)(load, from),
+                w2: (0, fixed_end_forces_1.loadIntensity)(load, to),
+                start: from - offset,
+                end: to - offset,
+            });
+        }
+    }
+    return out;
+};
 /**
  * Solves one load case against the prepared, factorised system.
  *
@@ -60,6 +96,14 @@ const toLocalLoads = (load, lambda, L) => {
  * @param loadCase - The load case (id and self-weight multiplier)
  * @param loads - The nodal and member loads tagged with this load case
  *
+ * Shell areas take `selfWeight · unitWeight · thickness · modifiers.weight` per unit area along
+ * −Z and their `AreaLoad`s as forces at their corners, each corner the share `∫ Nᵢ dA` of the
+ * element's area on its mean plane (no corner moments, so a warped element's offsets do not load
+ * the rotations auto-restraint holds).
+ *
+ * @param loads - The nodal and member loads tagged with this load case, and its area loads on
+ *   shell areas
+ *
  * @throws {FrameAnalysisError} 304 if a load has a moment about a joint rotation that nothing
  *   resists (an auto-restrained axis), 308 for a load outside its member, 313 if self-weight is
  *   requested for a member without `unitWeight`
@@ -67,15 +111,18 @@ const toLocalLoads = (load, lambda, L) => {
 const solveLoadCase = (system, loadCase, loads) => {
     const { members, dofCount, boundary, equation, free, factor, nodeIndex, memberIndex } = system;
     const F = new Array(dofCount).fill(0);
-    const loadsByMember = members.map(() => []);
+    const groups = system.memberGroups;
+    const loadsByMember = groups.map(() => []);
     for (const load of loads.memberLoads) {
         loadsByMember[memberIndex.get(load.memberId)].push(load);
     }
     const selfWeight = loadCase.selfWeight ?? 0;
-    const memberLoads = members.map((member, m) => {
-        const memberLoadsHere = [...loadsByMember[m]];
+    // local loads of each model member, then of each piece it is analysed as
+    const memberLoads = members.map(() => []);
+    groups.forEach((group, g) => {
+        const memberLoadsHere = [...loadsByMember[g]];
         if (selfWeight !== 0) {
-            const { unitWeight, A } = member.section;
+            const { unitWeight, A } = group.section;
             if (unitWeight === undefined)
                 throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.MISSING_UNIT_WEIGHT);
             // kN/m³ · m² = kN/m
@@ -84,19 +131,58 @@ const solveLoadCase = (system, loadCase, loads) => {
                 memberLoadsHere.push({
                     loadCase: loadCase.id,
                     type: "udl",
-                    memberId: member.id,
+                    memberId: group.id,
                     direction: "Z",
                     w,
                 });
             }
         }
-        return memberLoadsHere.flatMap((load) => toLocalLoads(load, member.lambda, member.length));
+        const local = memberLoadsHere.flatMap((load) => toLocalLoads(load, group.lambda, group.length));
+        if (group.pieces.length === 1) {
+            memberLoads[group.pieces[0]] = local;
+            return;
+        }
+        group.pieces.forEach((p, k) => {
+            const piece = members[p];
+            memberLoads[p] = clipLoads(local, piece.offset, piece.length, k === group.pieces.length - 1);
+        });
     });
     const fixed = members.map((member, m) => {
         const condensed = member.condenseFixedEndForces((0, fixed_end_forces_1.fixedEndForces)(memberLoads[m], member.length));
         (0, global_assembly_1.scatterMemberForces)(F, member, condensed, -1);
         return condensed;
     });
+    // Shell surface loads, kN/m² as a global vector per shell, then N at the corners
+    const shellPressure = system.shells.map(() => [0, 0, 0]);
+    for (const load of loads.shellLoads ?? []) {
+        // every element of the area
+        for (const s of system.shellsByArea.get(load.areaId) ?? []) {
+            const normal = system.shells[s].element.axes.lambda[2];
+            const e = load.direction === "3" ? normal : GLOBAL_DIRECTIONS[load.direction];
+            // per projected area → per true area: Aproj / A = |n · e|
+            const q = load.projected
+                ? load.q * Math.abs(normal[0] * e[0] + normal[1] * e[1] + normal[2] * e[2])
+                : load.q;
+            for (let k = 0; k < 3; k++)
+                shellPressure[s][k] += q * e[k];
+        }
+    }
+    system.shells.forEach((shell, s) => {
+        if (selfWeight !== 0) {
+            // kN/m³ · mm · 10⁻³ = kN/m²
+            const { unitWeight = 0, thickness } = shell.section;
+            shellPressure[s][2] -=
+                selfWeight * unitWeight * thickness * (shell.modifiers.weight ?? 1) * 1e-3;
+        }
+        const p = shellPressure[s];
+        if (p[0] === 0 && p[1] === 0 && p[2] === 0)
+            return;
+        shell.nodeIndices.forEach((node, i) => {
+            // kN/m² = 10⁻³ N/mm², times mm²
+            for (let k = 0; k < 3; k++)
+                F[6 * node + k] += p[k] * 1e-3 * shell.cornerAreas[i];
+        });
+    });
     for (const load of loads.nodalLoads) {
         const base = 6 * nodeIndex.get(load.nodeId);
         [load.fx, load.fy, load.fz].forEach((v, k) => (F[base + k] += (v ?? 0) * 1e3));
@@ -134,8 +220,14 @@ const solveLoadCase = (system, loadCase, loads) => {
         (0, global_assembly_1.scatterMemberForces)(KU, member, elastic, 1);
         return elastic.map((v, i) => v + fixed[m][i]);
     });
+    for (const shell of system.shells) {
+        const dofs = (0, global_assembly_1.elementDofs)(shell.nodeIndices);
+        const forces = (0, matrix_1.matVec)(shell.element.kGlobal, dofs.map((d) => U[d]));
+        dofs.forEach((d, a) => (KU[d] += forces[a]));
+    }
     const R = KU.map((v, i) => (boundary.restrained[i] ? v - F[i] : 0));
-    return { U, R, F, memberForces, memberLoads };
+    const shellForces = system.shells.map((shell) => (0, shell_forces_1.shellRawForces)(shell, U));
+    return { U, R, F, memberForces, shellForces, memberLoads };
 };
 exports.solveLoadCase = solveLoadCase;
 const scaleLoad = (load, factor) => {
@@ -160,6 +252,7 @@ const combineStates = (terms) => {
         R: add((s) => s.R, first.R.length),
         F: add((s) => s.F, first.F.length),
         memberForces: first.memberForces.map((_, m) => add((s) => s.memberForces[m], 12)),
+        shellForces: first.shellForces.map((f, sh) => add((s) => s.shellForces[sh], f.length)),
         memberLoads: first.memberLoads.map((_, m) => terms.flatMap(({ state, factor }) => factor === 0 ? [] : state.memberLoads[m].map((l) => scaleLoad(l, factor)))),
     };
 };