@theerapat-s28/ts-structural-eng-tools 1.5.0 → 2.1.0

package/dist/fea/frame/area/distribute-area-loads.js

@@ -0,0 +1,272 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.distributeAreaLoads = exports.areaMemberLoads = void 0;
+const frame_constant_1 = require("../../../core/constants/frame.constant");
+const validate_model_1 = require("../validation/validate-model");
+const area_geometry_1 = require("./area-geometry");
+const edge_members_1 = require("./edge-members");
+const GLOBAL_AXES = [
+    ["X", 0],
+    ["Y", 1],
+    ["Z", 2],
+];
+const dot2 = (a, b) => a[0] * b[0] + a[1] * b[1];
+/** Linear interpolation of a profile at `s` (0 outside it). */
+const valueAt = (profile, s) => {
+    for (let k = 1; k < profile.length; k++) {
+        const a = profile[k - 1];
+        const b = profile[k];
+        if (s >= a.s && s <= b.s)
+            return b.s > a.s ? a.w + ((b.w - a.w) * (s - a.s)) / (b.s - a.s) : b.w;
+    }
+    return 0;
+};
+/**
+ * Width of a convex polygon, given in (s, d) coordinates, across the s axis: the extent in d
+ * at each vertex position s. Between vertices the width is linear, so this is exact.
+ */
+const widthProfile = (polygon, tolerance) => {
+    const positions = [];
+    for (const s of polygon.map((p) => p[0]).sort((a, b) => a - b)) {
+        if (positions.length === 0 || s - positions[positions.length - 1] > tolerance)
+            positions.push(s);
+    }
+    return positions.map((s) => {
+        let min = Infinity;
+        let max = -Infinity;
+        polygon.forEach((p, i) => {
+            const q = polygon[(i + 1) % polygon.length];
+            const lo = Math.min(p[0], q[0]) - tolerance;
+            const hi = Math.max(p[0], q[0]) + tolerance;
+            if (s < lo || s > hi)
+                return;
+            const ds = q[0] - p[0];
+            const values = Math.abs(ds) <= tolerance ? [p[1], q[1]] : [p[1] + ((s - p[0]) / ds) * (q[1] - p[1])];
+            for (const d of values) {
+                min = Math.min(min, d);
+                max = Math.max(max, d);
+            }
+        });
+        return { s, w: max - min };
+    });
+};
+/** Unit inward normal of edge i of a counter-clockwise polygon (the area is on its left). */
+const inwardNormal = (corners, i) => {
+    const e = (0, area_geometry_1.edgeVector)(corners, i);
+    const L = Math.hypot(...e);
+    return [-e[1] / L, e[0] / L];
+};
+/** A profile pinned to the edge ends: clamped to [0, L], with breakpoints at 0 and L. */
+const pinToEdge = (profile, L) => {
+    const inside = profile.filter((p) => p.s > 0 && p.s < L);
+    return [{ s: 0, w: valueAt(profile, 0) }, ...inside, { s: L, w: valueAt(profile, L) }];
+};
+/**
+ * Two-way: every point of the area goes to the nearest edge line, so edge i takes the area
+ * clipped by `dᵢ(p) ≤ dⱼ(p)` for every other edge j (the corner-bisector, or straight-skeleton,
+ * split of a convex polygon), carried straight across onto the edge.
+ */
+const twoWayProfile = (plane, i, tolerance) => {
+    const { corners } = plane;
+    const ni = inwardNormal(corners, i);
+    let region = corners;
+    corners.forEach((cj, j) => {
+        if (j === i)
+            return;
+        const nj = inwardNormal(corners, j);
+        // dᵢ(p) ≤ dⱼ(p)  ⇔  (nᵢ − nⱼ)·p ≤ nᵢ·cᵢ − nⱼ·cⱼ
+        region = (0, area_geometry_1.clipPolygon)(region, [ni[0] - nj[0], ni[1] - nj[1]], dot2(ni, corners[i]) - dot2(nj, cj));
+    });
+    const e = (0, area_geometry_1.edgeVector)(corners, i);
+    const L = Math.hypot(...e);
+    const t = [e[0] / L, e[1] / L];
+    const local = region.map((p) => {
+        const r = [p[0] - corners[i][0], p[1] - corners[i][1]];
+        return [dot2(r, t), dot2(r, ni)];
+    });
+    return pinToEdge(widthProfile(local, tolerance), L);
+};
+/**
+ * One-way: strips parallel to `span` run across the area and half of each strip goes to each
+ * edge it ends on. A strip of width db meets edge i (at angle θ to the span) over ds = db/sin θ,
+ * so the edge takes `chord(b) / 2 · sin θ` per unit length; edges parallel to the span take none.
+ */
+const oneWayProfile = (plane, i, span, tolerance) => {
+    const { corners } = plane;
+    const across = [-span[1], span[0]];
+    const chord = widthProfile(corners.map((p) => [dot2(p, across), dot2(p, span)]), tolerance);
+    const e = (0, area_geometry_1.edgeVector)(corners, i);
+    const L = Math.hypot(...e);
+    const sin = Math.abs(e[0] * span[1] - e[1] * span[0]) / L;
+    if (sin <= frame_constant_1.AREA_GEOMETRY_TOLERANCE)
+        return null;
+    const b0 = dot2(corners[i], across);
+    const b1 = b0 + dot2(e, across);
+    const at = (b) => (valueAt(chord, b) / 2) * sin;
+    const inside = chord
+        .filter((p) => (p.s - b0) * (p.s - b1) < 0)
+        .map((p) => ({ s: ((p.s - b0) / (b1 - b0)) * L, w: at(p.s) }));
+    const profile = [{ s: 0, w: at(b0) }, ...inside, { s: L, w: at(b1) }];
+    return profile.sort((a, b) => a.s - b.s);
+};
+/**
+ * Cuts an edge profile at the members along the edge and maps each piece onto its member, as
+ * positions from the member start. Pieces that carry nothing are dropped.
+ */
+const segmentPieces = (profile, segment, memberLength, tolerance) => {
+    const { from, to, reversed } = segment;
+    const cuts = [
+        from,
+        ...profile.map((p) => p.s).filter((s) => s > from + tolerance && s < to - tolerance),
+        to,
+    ];
+    const x = (s) => {
+        const along = ((s - from) / (to - from)) * memberLength;
+        return Math.min(memberLength, Math.max(0, reversed ? memberLength - along : along));
+    };
+    const peak = Math.max(...profile.map((p) => Math.abs(p.w)));
+    const pieces = [];
+    for (let k = 1; k < cuts.length; k++) {
+        const wa = valueAt(profile, cuts[k - 1]);
+        const wb = valueAt(profile, cuts[k]);
+        if (Math.max(Math.abs(wa), Math.abs(wb)) <= frame_constant_1.AREA_GEOMETRY_TOLERANCE * peak)
+            continue;
+        const [xa, xb] = [x(cuts[k - 1]), x(cuts[k])];
+        pieces.push(xa < xb ? { start: xa, end: xb, w1: wa, w2: wb } : { start: xb, end: xa, w1: wb, w2: wa });
+    }
+    return pieces;
+};
+/**
+ * Turns the area loads of a model into member loads, assuming the model is valid (as
+ * `validateFrameModel` checks): each area's edge profiles are built once, then scaled for each
+ * of its loads and split into global X, Y and Z trapezoidal loads (plus torques for a cantilever).
+ * Everything is linear in `q`, so load-case superposition still holds.
+ */
+const areaMemberLoads = (model) => {
+    const areas = model.areas ?? [];
+    const areaLoads = model.areaLoads ?? [];
+    const memberLoads = [];
+    const warnings = [];
+    if (areas.length === 0 || areaLoads.length === 0)
+        return { memberLoads, warnings };
+    const nodeById = new Map(model.nodes.map((n) => [n.id, n]));
+    const position = (id) => {
+        const node = nodeById.get(id);
+        return [node.x, node.y, node.z];
+    };
+    const lengthOf = (member) => {
+        const [a, b] = [position(member.startNodeId), position(member.endNodeId)];
+        return Math.hypot(b[0] - a[0], b[1] - a[1], b[2] - a[2]);
+    };
+    const loadedAreas = new Set(areaLoads.map((l) => l.areaId));
+    const adjacency = (0, edge_members_1.memberAdjacency)(model.members);
+    const prepared = new Map();
+    for (const area of areas) {
+        if (!loadedAreas.has(area.id))
+            continue;
+        const ids = area.nodeIds;
+        const plane = (0, area_geometry_1.areaPlane)(ids.map(position));
+        const supports = (0, area_geometry_1.resolveAreaSupports)(plane, area.distribution, area.supportEdges);
+        const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * plane.size;
+        const n = ids.length;
+        const edges = [];
+        for (const i of supports.loadedEdges) {
+            const chain = (0, edge_members_1.edgeMemberChain)(ids[i], ids[(i + 1) % n], position, adjacency);
+            const L = Math.hypot(...(0, area_geometry_1.edgeVector)(plane.corners, i));
+            if (area.distribution === "cantilever") {
+                const depth = Math.hypot(...(0, area_geometry_1.edgeVector)(plane.corners, (i + 1) % n));
+                edges.push({
+                    chain,
+                    profile: [
+                        { s: 0, w: depth },
+                        { s: L, w: depth },
+                    ],
+                    torque: depth ** 2 / 2,
+                });
+            }
+            else if (area.distribution === "twoWay") {
+                edges.push({ chain, profile: twoWayProfile(plane, i, tolerance) });
+            }
+            else {
+                const profile = oneWayProfile(plane, i, supports.span, tolerance);
+                if (profile)
+                    edges.push({ chain, profile });
+                else {
+                    warnings.push({
+                        reference: "Null area, one-way distribution",
+                        message: `Area "${String(area.id)}": support edge ${i} runs parallel to the span and receives no load.`,
+                    });
+                }
+            }
+        }
+        prepared.set(area.id, { plane, edges });
+    }
+    for (const load of areaLoads) {
+        const area = prepared.get(load.areaId);
+        if (!area || load.q === 0)
+            continue;
+        const { plane, edges } = area;
+        const e = load.direction === "3"
+            ? plane.normal
+            : GLOBAL_AXES.map(([axis]) => (axis === load.direction ? 1 : 0));
+        const components = GLOBAL_AXES.filter(([, k]) => Math.abs(e[k]) > frame_constant_1.AREA_GEOMETRY_TOLERANCE);
+        // Lever arm (inward, in plane) × load: only the normal part of the load twists the edge,
+        // about +edge direction; validation rejects a load along a cantilever's support edge.
+        const normalPart = e[0] * plane.normal[0] + e[1] * plane.normal[1] + e[2] * plane.normal[2];
+        const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * plane.size;
+        for (const { chain, profile, torque } of edges) {
+            for (const segment of chain) {
+                const { member } = segment;
+                const L = lengthOf(member);
+                for (const piece of segmentPieces(profile, segment, L, tolerance)) {
+                    for (const [axis, k] of components) {
+                        memberLoads.push({
+                            loadCase: load.loadCase,
+                            type: "trapezoidal",
+                            memberId: member.id,
+                            direction: axis,
+                            w1: load.q * e[k] * piece.w1 * 1e-3, // kN/m² · mm → kN/m
+                            w2: load.q * e[k] * piece.w2 * 1e-3,
+                            start: piece.start,
+                            end: piece.end,
+                        });
+                    }
+                }
+                if (torque !== undefined && Math.abs(normalPart) > frame_constant_1.AREA_GEOMETRY_TOLERANCE) {
+                    // member local 1 is ±edge direction
+                    const sign = segment.reversed ? -1 : 1;
+                    memberLoads.push({
+                        loadCase: load.loadCase,
+                        type: "torque",
+                        memberId: member.id,
+                        t: load.q * normalPart * torque * sign * 1e-6, // kN/m² · mm² → kN·m/m
+                    });
+                }
+            }
+        }
+    }
+    return { memberLoads, warnings };
+};
+exports.areaMemberLoads = areaMemberLoads;
+/**
+ * The member loads that a model's null areas send to the members along their edges — what
+ * `analyzeFrame` adds to `memberLoads` — so they can be checked or previewed.
+ *
+ * Two-way areas send each point to the nearest edge (corner-bisector split); one-way areas
+ * split strips parallel to their unsupported edges half to each end; cantilever areas send
+ * `q · c` and a torque `q · c² / 2` per unit length to their support edge (c = depth). Edges
+ * split by intermediate nodes are cut at those nodes. Loads are per unit true area and come
+ * out as trapezoidal loads in global X, Y and Z (and torques), each tagged with its load case.
+ *
+ * @param model - The frame model, validated first
+ * @returns `memberLoads` and `warnings` (a one-way support edge parallel to the span)
+ * @throws {FrameAnalysisError} 301–317 for an invalid model, with every problem in `issues`
+ */
+const distributeAreaLoads = (model) => {
+    const issues = (0, validate_model_1.validateFrameModel)(model);
+    if (issues.length > 0)
+        (0, validate_model_1.throwIssues)(issues);
+    return (0, exports.areaMemberLoads)(model);
+};
+exports.distributeAreaLoads = distributeAreaLoads;
+//# sourceMappingURL=distribute-area-loads.js.map

package/dist/fea/frame/area/distribute-area-loads.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"distribute-area-loads.js","sourceRoot":"","sources":["../../../../src/fea/frame/area/distribute-area-loads.ts"],"names":[],"mappings":";;;AAAA,uEAA6E;AAI7E,iEAA+E;AAC/E,mDAQyB;AACzB,iDAAoF;AAepF,MAAM,WAAW,GAAgC;IAC/C,CAAC,GAAG,EAAE,CAAC,CAAC;IACR,CAAC,GAAG,EAAE,CAAC,CAAC;IACR,CAAC,GAAG,EAAE,CAAC,CAAC;CACT,CAAC;AAEF,MAAM,IAAI,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAE7D,+DAA+D;AAC/D,MAAM,OAAO,GAAG,CAAC,OAAgB,EAAE,CAAS,EAAU,EAAE;IACtD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACrB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;YACtB,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3E,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,YAAY,GAAG,CAAC,OAAe,EAAE,SAAiB,EAAW,EAAE;IACnE,MAAM,SAAS,GAAa,EAAE,CAAC;IAC/B,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QAC/D,IAAI,SAAS,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,GAAG,SAAS,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,SAAS;YAC3E,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IACD,OAAO,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;QACzB,IAAI,GAAG,GAAG,QAAQ,CAAC;QACnB,IAAI,GAAG,GAAG,CAAC,QAAQ,CAAC;QACpB,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YACvB,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC;YAC5C,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;YAC5C,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;YAC5C,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE;gBAAE,OAAO;YAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACvB,MAAM,MAAM,GACV,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACxF,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;gBACvB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;gBACvB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YACzB,CAAC;QACH,CAAC,CAAC,CAAC;QACH,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;IAC7B,CAAC,CAAC,CAAC;AACL,CAAC,CAAC;AAEF,6FAA6F;AAC7F,MAAM,YAAY,GAAG,CAAC,OAAe,EAAE,CAAS,EAAQ,EAAE;IACxD,MAAM,CAAC,GAAG,IAAA,0BAAU,EAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IACjC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AAC/B,CAAC,CAAC;AAEF,yFAAyF;AACzF,MAAM,SAAS,GAAG,CAAC,OAAgB,EAAE,CAAS,EAAW,EAAE;IACzD,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACzD,OAAO,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,OAAO,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;AACzF,CAAC,CAAC;AAEF;;;;GAIG;AACH,MAAM,aAAa,GAAG,CAAC,KAAgB,EAAE,CAAS,EAAE,SAAiB,EAAW,EAAE;IAChF,MAAM,EAAE,OAAO,EAAE,GAAG,KAAK,CAAC;IAC1B,MAAM,EAAE,GAAG,YAAY,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IACpC,IAAI,MAAM,GAAG,OAAO,CAAC;IACrB,OAAO,CAAC,OAAO,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE;QACxB,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO;QACpB,MAAM,EAAE,GAAG,YAAY,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;QACpC,gDAAgD;QAChD,MAAM,GAAG,IAAA,2BAAW,EAClB,MAAM,EACN,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAC9B,IAAI,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,CACpC,CAAC;IACJ,CAAC,CAAC,CAAC;IACH,MAAM,CAAC,GAAG,IAAA,0BAAU,EAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IACjC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAQ,EAAE;QACnC,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7D,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACnC,CAAC,CAAC,CAAC;IACH,OAAO,SAAS,CAAC,YAAY,CAAC,KAAK,EAAE,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC;AACtD,CAAC,CAAC;AAEF;;;;GAIG;AACH,MAAM,aAAa,GAAG,CACpB,KAAgB,EAChB,CAAS,EACT,IAAU,EACV,SAAiB,EACD,EAAE;IAClB,MAAM,EAAE,OAAO,EAAE,GAAG,KAAK,CAAC;IAC1B,MAAM,MAAM,GAAS,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,YAAY,CACxB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAQ,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,MAAM,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAC1D,SAAS,CACV,CAAC;IACF,MAAM,CAAC,GAAG,IAAA,0BAAU,EAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IACjC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC1D,IAAI,GAAG,IAAI,wCAAuB;QAAE,OAAO,IAAI,CAAC;IAChD,MAAM,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IAChC,MAAM,EAAE,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;IACxD,MAAM,MAAM,GAAG,KAAK;SACjB,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;SAC1C,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACjE,MAAM,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACtE,OAAO,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC3C,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,aAAa,GAAG,CACpB,OAAgB,EAChB,OAAoB,EACpB,YAAoB,EACpB,SAAiB,EACyC,EAAE;IAC5D,MAAM,EAAE,IAAI,EAAE,EAAE,EAAE,QAAQ,EAAE,GAAG,OAAO,CAAC;IACvC,MAAM,IAAI,GAAG;QACX,IAAI;QACJ,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,SAAS,IAAI,CAAC,GAAG,EAAE,GAAG,SAAS,CAAC;QACpF,EAAE;KACH,CAAC;IACF,MAAM,CAAC,GAAG,CAAC,CAAS,EAAE,EAAE;QACtB,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,YAAY,CAAC;QACxD,OAAO,IAAI,CAAC,GAAG,CAAC,YAAY,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,YAAY,GAAG,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;IACtF,CAAC,CAAC;IACF,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,MAAM,MAAM,GAAG,EAAE,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACrC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,IAAI,wCAAuB,GAAG,IAAI;YAAE,SAAS;QACrF,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,MAAM,CAAC,IAAI,CACT,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAC1F,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC,CAAC;AAEF;;;;;GAKG;AACI,MAAM,eAAe,GAAG,CAC7B,KAAiB,EACkC,EAAE;IACrD,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,IAAI,EAAE,CAAC;IAChC,MAAM,SAAS,GAAG,KAAK,CAAC,SAAS,IAAI,EAAE,CAAC;IACxC,MAAM,WAAW,GAAiB,EAAE,CAAC;IACrC,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,SAAS,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,EAAE,WAAW,EAAE,QAAQ,EAAE,CAAC;IAEnF,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,MAAM,QAAQ,GAAG,CAAC,EAAW,EAAQ,EAAE;QACrC,MAAM,IAAI,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC;QAC/B,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC;IACF,MAAM,QAAQ,GAAG,CAAC,MAAc,EAAE,EAAE;QAClC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,WAAW,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;QAC1E,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3D,CAAC,CAAC;IAEF,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;IAC5D,MAAM,SAAS,GAAG,IAAA,8BAAe,EAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IACjD,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAsD,CAAC;IAC/E,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;YAAE,SAAS;QACxC,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC;QACzB,MAAM,KAAK,GAAG,IAAA,yBAAS,EAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAc,CAAC;QACxD,MAAM,QAAQ,GAAG,IAAA,mCAAmB,EAAC,KAAK,EAAE,IAAI,CAAC,YAAY,EAAE,IAAI,CAAC,YAAY,CAAE,CAAC;QACnF,MAAM,SAAS,GAAG,wCAAuB,GAAG,KAAK,CAAC,IAAI,CAAC;QACvD,MAAM,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC;QAErB,MAAM,KAAK,GAAiB,EAAE,CAAC;QAC/B,KAAK,MAAM,CAAC,IAAI,QAAQ,CAAC,WAAW,EAAE,CAAC;YACrC,MAAM,KAAK,GAAG,IAAA,8BAAe,EAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,EAAE,SAAS,CAAE,CAAC;YAC9E,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,IAAA,0BAAU,EAAC,KAAK,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC;YAEtD,IAAI,IAAI,CAAC,YAAY,KAAK,YAAY,EAAE,CAAC;gBACvC,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,IAAA,0BAAU,EAAC,KAAK,CAAC,OAAO,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBACpE,KAAK,CAAC,IAAI,CAAC;oBACT,KAAK;oBACL,OAAO,EAAE;wBACP,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE;wBAClB,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE;qBACnB;oBACD,MAAM,EAAE,KAAK,IAAI,CAAC,GAAG,CAAC;iBACvB,CAAC,CAAC;YACL,CAAC;iBAAM,IAAI,IAAI,CAAC,YAAY,KAAK,QAAQ,EAAE,CAAC;gBAC1C,KAAK,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,OAAO,EAAE,aAAa,CAAC,KAAK,EAAE,CAAC,EAAE,SAAS,CAAC,EAAE,CAAC,CAAC;YACrE,CAAC;iBAAM,CAAC;gBACN,MAAM,OAAO,GAAG,aAAa,CAAC,KAAK,EAAE,CAAC,EAAE,QAAQ,CAAC,IAAK,EAAE,SAAS,CAAC,CAAC;gBACnE,IAAI,OAAO;oBAAE,KAAK,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,OAAO,EAAE,CAAC,CAAC;qBACvC,CAAC;oBACJ,QAAQ,CAAC,IAAI,CAAC;wBACZ,SAAS,EAAE,iCAAiC;wBAC5C,OAAO,EAAE,SAAS,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,mBAAmB,CAAC,kDAAkD;qBACxG,CAAC,CAAC;gBACL,CAAC;YACH,CAAC;QACH,CAAC;QACD,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC,CAAC;IAC1C,CAAC;IAED,KAAK,MAAM,IAAI,IAAI,SAAS,EAAE,CAAC;QAC7B,MAAM,IAAI,GAAG,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACvC,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,CAAC,KAAK,CAAC;YAAE,SAAS;QACpC,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,GAAG,IAAI,CAAC;QAC9B,MAAM,CAAC,GACL,IAAI,CAAC,SAAS,KAAK,GAAG;YACpB,CAAC,CAAC,KAAK,CAAC,MAAM;YACd,CAAC,CAAE,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,KAAK,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAU,CAAC;QAC/E,MAAM,UAAU,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,wCAAuB,CAAC,CAAC;QAC3F,yFAAyF;QACzF,sFAAsF;QACtF,MAAM,UAAU,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC5F,MAAM,SAAS,GAAG,wCAAuB,GAAG,KAAK,CAAC,IAAI,CAAC;QAEvD,KAAK,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,MAAM,EAAE,IAAI,KAAK,EAAE,CAAC;YAC/C,KAAK,MAAM,OAAO,IAAI,KAAK,EAAE,CAAC;gBAC5B,MAAM,EAAE,MAAM,EAAE,GAAG,OAAO,CAAC;gBAC3B,MAAM,CAAC,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC;gBAC3B,KAAK,MAAM,KAAK,IAAI,aAAa,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC,EAAE,SAAS,CAAC,EAAE,CAAC;oBAClE,KAAK,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,UAAU,EAAE,CAAC;wBACnC,WAAW,CAAC,IAAI,CAAC;4BACf,QAAQ,EAAE,IAAI,CAAC,QAAQ;4BACvB,IAAI,EAAE,aAAa;4BACnB,QAAQ,EAAE,MAAM,CAAC,EAAE;4BACnB,SAAS,EAAE,IAAqB;4BAChC,EAAE,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,IAAI,EAAE,oBAAoB;4BACzD,EAAE,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,IAAI;4BACnC,KAAK,EAAE,KAAK,CAAC,KAAK;4BAClB,GAAG,EAAE,KAAK,CAAC,GAAG;yBACf,CAAC,CAAC;oBACL,CAAC;gBACH,CAAC;gBACD,IAAI,MAAM,KAAK,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,GAAG,wCAAuB,EAAE,CAAC;oBAC3E,oCAAoC;oBACpC,MAAM,IAAI,GAAG,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACvC,WAAW,CAAC,IAAI,CAAC;wBACf,QAAQ,EAAE,IAAI,CAAC,QAAQ;wBACvB,IAAI,EAAE,QAAQ;wBACd,QAAQ,EAAE,MAAM,CAAC,EAAE;wBACnB,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,UAAU,GAAG,MAAM,GAAG,IAAI,GAAG,IAAI,EAAE,uBAAuB;qBACvE,CAAC,CAAC;gBACL,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,EAAE,WAAW,EAAE,QAAQ,EAAE,CAAC;AACnC,CAAC,CAAC;AA5GW,QAAA,eAAe,mBA4G1B;AAEF;;;;;;;;;;;;;GAaG;AACI,MAAM,mBAAmB,GAAG,CACjC,KAAiB,EACkC,EAAE;IACrD,MAAM,MAAM,GAAG,IAAA,mCAAkB,EAAC,KAAK,CAAC,CAAC;IACzC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC;QAAE,IAAA,4BAAW,EAAC,MAAM,CAAC,CAAC;IAC3C,OAAO,IAAA,uBAAe,EAAC,KAAK,CAAC,CAAC;AAChC,CAAC,CAAC;AANW,QAAA,mBAAmB,uBAM9B"}

package/dist/fea/frame/area/edge-members.d.ts

@@ -0,0 +1,25 @@
+import type { FrameId, Member } from "../types/model.type";
+import type { Vec3 } from "./area-geometry";
+/** The members meeting at each node, with the node at their other end. */
+export type MemberAdjacency = Map<FrameId, {
+    member: Member;
+    other: FrameId;
+}[]>;
+/** One member covering part of an area edge, with positions along the edge from its start. */
+export interface EdgeSegment {
+    member: Member;
+    from: number;
+    to: number;
+    reversed: boolean;
+}
+/** Indexes members by the nodes they connect. */
+export declare const memberAdjacency: (members: Member[]) => MemberAdjacency;
+/**
+ * The chain of collinear members covering an area edge from `start` to `end`, found by walking
+ * the members at each node towards `end`: each step takes the member whose far node lies on the
+ * edge nearest ahead, so intermediate nodes (beams framing in) split the edge into segments.
+ *
+ * @param position - Node coordinates in mm
+ * @returns The segments in order along the edge, or `null` if a gap leaves part of it uncovered
+ */
+export declare const edgeMemberChain: (start: FrameId, end: FrameId, position: (id: FrameId) => Vec3, adjacency: MemberAdjacency) => EdgeSegment[] | null;

package/dist/fea/frame/area/edge-members.js

@@ -0,0 +1,67 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.edgeMemberChain = exports.memberAdjacency = void 0;
+const frame_constant_1 = require("../../../core/constants/frame.constant");
+/** Indexes members by the nodes they connect. */
+const memberAdjacency = (members) => {
+    const adjacency = new Map();
+    const link = (node, entry) => {
+        const list = adjacency.get(node);
+        if (list)
+            list.push(entry);
+        else
+            adjacency.set(node, [entry]);
+    };
+    for (const member of members) {
+        link(member.startNodeId, { member, other: member.endNodeId });
+        link(member.endNodeId, { member, other: member.startNodeId });
+    }
+    return adjacency;
+};
+exports.memberAdjacency = memberAdjacency;
+/**
+ * The chain of collinear members covering an area edge from `start` to `end`, found by walking
+ * the members at each node towards `end`: each step takes the member whose far node lies on the
+ * edge nearest ahead, so intermediate nodes (beams framing in) split the edge into segments.
+ *
+ * @param position - Node coordinates in mm
+ * @returns The segments in order along the edge, or `null` if a gap leaves part of it uncovered
+ */
+const edgeMemberChain = (start, end, position, adjacency) => {
+    const A = position(start);
+    const B = position(end);
+    const d = [B[0] - A[0], B[1] - A[1], B[2] - A[2]];
+    const L = Math.hypot(...d);
+    const t = [d[0] / L, d[1] / L, d[2] / L];
+    const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * L;
+    const segments = [];
+    let current = start;
+    let s = 0;
+    while (current !== end) {
+        let next;
+        for (const { member, other } of adjacency.get(current) ?? []) {
+            const P = position(other);
+            const r = [P[0] - A[0], P[1] - A[1], P[2] - A[2]];
+            const along = r[0] * t[0] + r[1] * t[1] + r[2] * t[2];
+            const off = Math.hypot(r[0] - along * t[0], r[1] - along * t[1], r[2] - along * t[2]);
+            if (off > tolerance || along <= s + tolerance || along > L + tolerance)
+                continue;
+            if (!next || along < next.s)
+                next = { member, other, s: along };
+        }
+        if (!next)
+            return null;
+        const to = next.other === end ? L : next.s;
+        segments.push({
+            member: next.member,
+            from: s,
+            to,
+            reversed: next.member.startNodeId !== current,
+        });
+        current = next.other;
+        s = to;
+    }
+    return segments;
+};
+exports.edgeMemberChain = edgeMemberChain;
+//# sourceMappingURL=edge-members.js.map

package/dist/fea/frame/area/edge-members.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"edge-members.js","sourceRoot":"","sources":["../../../../src/fea/frame/area/edge-members.ts"],"names":[],"mappings":";;;AAAA,uEAA6E;AAgB7E,iDAAiD;AAC1C,MAAM,eAAe,GAAG,CAAC,OAAiB,EAAmB,EAAE;IACpE,MAAM,SAAS,GAAoB,IAAI,GAAG,EAAE,CAAC;IAC7C,MAAM,IAAI,GAAG,CAAC,IAAa,EAAE,KAAyC,EAAE,EAAE;QACxE,MAAM,IAAI,GAAG,SAAS,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACjC,IAAI,IAAI;YAAE,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;;YACtB,SAAS,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC;IACpC,CAAC,CAAC;IACF,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE,CAAC;QAC7B,IAAI,CAAC,MAAM,CAAC,WAAW,EAAE,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC;QAC9D,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,CAAC,WAAW,EAAE,CAAC,CAAC;IAChE,CAAC;IACD,OAAO,SAAS,CAAC;AACnB,CAAC,CAAC;AAZW,QAAA,eAAe,mBAY1B;AAEF;;;;;;;GAOG;AACI,MAAM,eAAe,GAAG,CAC7B,KAAc,EACd,GAAY,EACZ,QAA+B,EAC/B,SAA0B,EACJ,EAAE;IACxB,MAAM,CAAC,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;IACxB,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxD,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC/C,MAAM,SAAS,GAAG,wCAAuB,GAAG,CAAC,CAAC;IAE9C,MAAM,QAAQ,GAAkB,EAAE,CAAC;IACnC,IAAI,OAAO,GAAG,KAAK,CAAC;IACpB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,OAAO,OAAO,KAAK,GAAG,EAAE,CAAC;QACvB,IAAI,IAA+D,CAAC;QACpE,KAAK,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,IAAI,SAAS,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,EAAE,EAAE,CAAC;YAC7D,MAAM,CAAC,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;YAC1B,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACxD,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACtD,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtF,IAAI,GAAG,GAAG,SAAS,IAAI,KAAK,IAAI,CAAC,GAAG,SAAS,IAAI,KAAK,GAAG,CAAC,GAAG,SAAS;gBAAE,SAAS;YACjF,IAAI,CAAC,IAAI,IAAI,KAAK,GAAG,IAAI,CAAC,CAAC;gBAAE,IAAI,GAAG,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC;QAClE,CAAC;QACD,IAAI,CAAC,IAAI;YAAE,OAAO,IAAI,CAAC;QACvB,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;QAC3C,QAAQ,CAAC,IAAI,CAAC;YACZ,MAAM,EAAE,IAAI,CAAC,MAAM;YACnB,IAAI,EAAE,CAAC;YACP,EAAE;YACF,QAAQ,EAAE,IAAI,CAAC,MAAM,CAAC,WAAW,KAAK,OAAO;SAC9C,CAAC,CAAC;QACH,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC;QACrB,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC;AAtCW,QAAA,eAAe,mBAsC1B"}

package/dist/fea/frame/elements/beam-3d-stiffness.d.ts

@@ -0,0 +1,29 @@
+import { type Matrix } from "../../../utils/matrix";
+/**
+ * Local stiffness matrix of a 3D prismatic beam element, exactly as in the reference
+ * "Complete Beam Stiffness Matrix".
+ *
+ * DOF order for node i then node j: [U1, U2, U3, R1, R2, R3]. Local 1 runs along the member.
+ * The reference matrix measures the bending rotations R2 and R3 opposite to the right-hand
+ * rule; use `beam3dLocalStiffnessRightHanded` for the right-handed form.
+ *
+ * The reference uses one `EI`; here the U2–R3 terms use `E·I33` and the U3–R2 terms use
+ * `E·I22`, which reduces to the reference when `I22 = I33`.
+ *
+ * @param E - Young's modulus in MPa
+ * @param G - Shear modulus in MPa
+ * @param A - Area in mm²
+ * @param I22 - Second moment about local 2 in mm⁴
+ * @param I33 - Second moment about local 3 in mm⁴
+ * @param J - Torsion constant in mm⁴
+ * @param L - Member length in mm
+ * @returns 12×12 symmetric matrix in N and mm
+ */
+export declare const beam3dLocalStiffness: (E: number, G: number, A: number, I22: number, I33: number, J: number, L: number) => Matrix;
+/**
+ * Converts a local matrix between the reference rotation convention and the right-hand rule
+ * (`D·k·D`, with `D` flipping the sign of R2 and R3 at both ends). The conversion is its own inverse.
+ */
+export declare const flipBendingRotationSigns: (k: Matrix) => Matrix;
+/** Same element as `beam3dLocalStiffness` with all rotations following the right-hand rule. */
+export declare const beam3dLocalStiffnessRightHanded: (E: number, G: number, A: number, I22: number, I33: number, J: number, L: number) => Matrix;

package/dist/fea/frame/elements/beam-3d-stiffness.js

@@ -0,0 +1,90 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.beam3dLocalStiffnessRightHanded = exports.flipBendingRotationSigns = exports.beam3dLocalStiffness = void 0;
+const matrix_1 = require("../../../utils/matrix");
+/**
+ * Local DOF indices (0-based) whose reference-matrix rotation sign is opposite to the
+ * right-hand rule: R2 and R3 at the start and end nodes.
+ */
+const FLIPPED_DOFS = [4, 5, 10, 11];
+/**
+ * Local stiffness matrix of a 3D prismatic beam element, exactly as in the reference
+ * "Complete Beam Stiffness Matrix".
+ *
+ * DOF order for node i then node j: [U1, U2, U3, R1, R2, R3]. Local 1 runs along the member.
+ * The reference matrix measures the bending rotations R2 and R3 opposite to the right-hand
+ * rule; use `beam3dLocalStiffnessRightHanded` for the right-handed form.
+ *
+ * The reference uses one `EI`; here the U2–R3 terms use `E·I33` and the U3–R2 terms use
+ * `E·I22`, which reduces to the reference when `I22 = I33`.
+ *
+ * @param E - Young's modulus in MPa
+ * @param G - Shear modulus in MPa
+ * @param A - Area in mm²
+ * @param I22 - Second moment about local 2 in mm⁴
+ * @param I33 - Second moment about local 3 in mm⁴
+ * @param J - Torsion constant in mm⁴
+ * @param L - Member length in mm
+ * @returns 12×12 symmetric matrix in N and mm
+ */
+const beam3dLocalStiffness = (E, G, A, I22, I33, J, L) => {
+    const k = (0, matrix_1.zeros)(12, 12);
+    // 1-based row/column, upper triangle as written in the reference; mirrored for symmetry.
+    const set = (row, col, value) => {
+        k[row - 1][col - 1] = value;
+        k[col - 1][row - 1] = value;
+    };
+    const ax = (E * A) / L;
+    const tor = (G * J) / L;
+    // U2–R3 bending (I33)
+    const a3 = (12 * E * I33) / L ** 3;
+    const b3 = (6 * E * I33) / L ** 2;
+    const c3 = (4 * E * I33) / L;
+    const d3 = (2 * E * I33) / L;
+    // U3–R2 bending (I22)
+    const a2 = (12 * E * I22) / L ** 3;
+    const b2 = (6 * E * I22) / L ** 2;
+    const c2 = (4 * E * I22) / L;
+    const d2 = (2 * E * I22) / L;
+    set(1, 1, ax);
+    set(1, 7, -ax);
+    set(7, 7, ax);
+    set(2, 2, a3);
+    set(2, 6, -b3);
+    set(2, 8, -a3);
+    set(2, 12, -b3);
+    set(6, 6, c3);
+    set(6, 8, b3);
+    set(6, 12, d3);
+    set(8, 8, a3);
+    set(8, 12, b3);
+    set(12, 12, c3);
+    set(3, 3, a2);
+    set(3, 5, b2);
+    set(3, 9, -a2);
+    set(3, 11, b2);
+    set(5, 5, c2);
+    set(5, 9, -b2);
+    set(5, 11, d2);
+    set(9, 9, a2);
+    set(9, 11, -b2);
+    set(11, 11, c2);
+    set(4, 4, tor);
+    set(4, 10, -tor);
+    set(10, 10, tor);
+    return k;
+};
+exports.beam3dLocalStiffness = beam3dLocalStiffness;
+/**
+ * Converts a local matrix between the reference rotation convention and the right-hand rule
+ * (`D·k·D`, with `D` flipping the sign of R2 and R3 at both ends). The conversion is its own inverse.
+ */
+const flipBendingRotationSigns = (k) => {
+    const flipped = new Set(FLIPPED_DOFS);
+    return k.map((row, i) => row.map((value, j) => (flipped.has(i) !== flipped.has(j) ? -value : value)));
+};
+exports.flipBendingRotationSigns = flipBendingRotationSigns;
+/** Same element as `beam3dLocalStiffness` with all rotations following the right-hand rule. */
+const beam3dLocalStiffnessRightHanded = (E, G, A, I22, I33, J, L) => (0, exports.flipBendingRotationSigns)((0, exports.beam3dLocalStiffness)(E, G, A, I22, I33, J, L));
+exports.beam3dLocalStiffnessRightHanded = beam3dLocalStiffnessRightHanded;
+//# sourceMappingURL=beam-3d-stiffness.js.map

package/dist/fea/frame/elements/beam-3d-stiffness.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"beam-3d-stiffness.js","sourceRoot":"","sources":["../../../../src/fea/frame/elements/beam-3d-stiffness.ts"],"names":[],"mappings":";;;AAAA,8CAAuD;AAEvD;;;GAGG;AACH,MAAM,YAAY,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;AAEpC;;;;;;;;;;;;;;;;;;;GAmBG;AACI,MAAM,oBAAoB,GAAG,CAClC,CAAS,EACT,CAAS,EACT,CAAS,EACT,GAAW,EACX,GAAW,EACX,CAAS,EACT,CAAS,EACD,EAAE;IACV,MAAM,CAAC,GAAG,IAAA,cAAK,EAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,yFAAyF;IACzF,MAAM,GAAG,GAAG,CAAC,GAAW,EAAE,GAAW,EAAE,KAAa,EAAE,EAAE;QACtD,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QAC5B,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;IAC9B,CAAC,CAAC;IAEF,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IAExB,sBAAsB;IACtB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,sBAAsB;IACtB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAE7B,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IAEd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAEhB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAEhB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACjB,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IAEjB,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AA7DW,QAAA,oBAAoB,wBA6D/B;AAEF;;;GAGG;AACI,MAAM,wBAAwB,GAAG,CAAC,CAAS,EAAU,EAAE;IAC5D,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,YAAY,CAAC,CAAC;IACtC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CACtB,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAC5E,CAAC;AACJ,CAAC,CAAC;AALW,QAAA,wBAAwB,4BAKnC;AAEF,+FAA+F;AACxF,MAAM,+BAA+B,GAAG,CAC7C,CAAS,EACT,CAAS,EACT,CAAS,EACT,GAAW,EACX,GAAW,EACX,CAAS,EACT,CAAS,EACD,EAAE,CAAC,IAAA,gCAAwB,EAAC,IAAA,4BAAoB,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AARxE,QAAA,+BAA+B,mCAQyC"}

package/dist/fea/frame/elements/fixed-end-forces.d.ts

@@ -0,0 +1,36 @@
+/**
+ * A member load resolved into local axes, in N and mm. A distributed load varies linearly
+ * from `w1` at `start` to `w2` at `end` (uniform when they are equal); a torque is uniform
+ * about local 1.
+ */
+export type LocalMemberLoad = {
+    type: "distributed";
+    axis: 1 | 2 | 3;
+    w1: number;
+    w2: number;
+    start: number;
+    end: number;
+} | {
+    type: "point";
+    axis: 1 | 2 | 3;
+    P: number;
+    a: number;
+} | {
+    type: "torque";
+    t: number;
+    start: number;
+    end: number;
+};
+/** Intensity of a distributed load at `x` (inside its extent), in N/mm. */
+export declare const loadIntensity: (load: Extract<LocalMemberLoad, {
+    type: "distributed";
+}>, x: number) => number;
+/**
+ * Fixed-end forces of a member under local loads: the end forces (right-hand local axes,
+ * N and N·mm) that hold both ends fully fixed. Equivalent nodal loads are the negative.
+ *
+ * @param loads - Loads in local axes, N and mm
+ * @param L - Member length in mm
+ * @throws {FrameAnalysisError} 308 if a load lies outside the member
+ */
+export declare const fixedEndForces: (loads: LocalMemberLoad[], L: number) => number[];

package/dist/fea/frame/elements/fixed-end-forces.js

@@ -0,0 +1,93 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.fixedEndForces = exports.loadIntensity = void 0;
+const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.error");
+/** Intensity of a distributed load at `x` (inside its extent), in N/mm. */
+const loadIntensity = (load, x) => load.w1 + ((load.w2 - load.w1) * (x - load.start)) / (load.end - load.start);
+exports.loadIntensity = loadIntensity;
+// 4-point Gauss–Legendre rule, exact up to degree 7: cubic shape functions × linear loads
+const GAUSS_POINTS = [
+    -0.8611363115940526, -0.3399810435848563, 0.3399810435848563, 0.8611363115940526,
+];
+const GAUSS_WEIGHTS = [
+    0.3478548451374538, 0.6521451548625461, 0.6521451548625461, 0.3478548451374538,
+];
+/**
+ * Equivalent nodal loads of a unit transverse/axial load at distance `x`, from the
+ * Euler–Bernoulli shape functions: the exact fixed-end solution for a prismatic member.
+ * Returned per unit of load, right-hand local DOF order [i: U1..R3, j: U1..R3].
+ */
+const unitEquivalentLoad = (axis, x, L) => {
+    const xi = x / L;
+    const out = new Array(12).fill(0);
+    if (axis === 1) {
+        out[0] = 1 - xi;
+        out[6] = xi;
+        return out;
+    }
+    const n1 = 1 - 3 * xi ** 2 + 2 * xi ** 3;
+    const n2 = L * (xi - 2 * xi ** 2 + xi ** 3);
+    const n3 = 3 * xi ** 2 - 2 * xi ** 3;
+    const n4 = L * (-(xi ** 2) + xi ** 3);
+    if (axis === 2) {
+        // load along local 2: force U2, moment about 3 (R3 = +dv/dx)
+        out[1] = n1;
+        out[5] = n2;
+        out[7] = n3;
+        out[11] = n4;
+    }
+    else {
+        // load along local 3: force U3, moment about 2 (R2 = −dw/dx)
+        out[2] = n1;
+        out[4] = -n2;
+        out[8] = n3;
+        out[10] = -n4;
+    }
+    return out;
+};
+/**
+ * Fixed-end forces of a member under local loads: the end forces (right-hand local axes,
+ * N and N·mm) that hold both ends fully fixed. Equivalent nodal loads are the negative.
+ *
+ * @param loads - Loads in local axes, N and mm
+ * @param L - Member length in mm
+ * @throws {FrameAnalysisError} 308 if a load lies outside the member
+ */
+const fixedEndForces = (loads, L) => {
+    const equivalent = new Array(12).fill(0);
+    const add = (unit, factor) => {
+        for (let i = 0; i < 12; i++)
+            equivalent[i] += unit[i] * factor;
+    };
+    for (const load of loads) {
+        if (load.type === "torque") {
+            if (!(load.start >= 0 && load.end <= L && load.start < load.end)) {
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.LOAD_OUTSIDE_MEMBER);
+            }
+            // linear torsion shape functions 1 − ξ and ξ, integrated over [start, end]
+            const resultant = load.t * (load.end - load.start);
+            const centroid = (load.start + load.end) / 2;
+            equivalent[3] += resultant * (1 - centroid / L);
+            equivalent[9] += (resultant * centroid) / L;
+        }
+        else if (load.type === "point") {
+            if (!(load.a >= 0 && load.a <= L))
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.LOAD_OUTSIDE_MEMBER);
+            add(unitEquivalentLoad(load.axis, load.a, L), load.P);
+        }
+        else {
+            if (!(load.start >= 0 && load.end <= L && load.start < load.end)) {
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.LOAD_OUTSIDE_MEMBER);
+            }
+            const half = (load.end - load.start) / 2;
+            const mid = (load.end + load.start) / 2;
+            GAUSS_POINTS.forEach((p, index) => {
+                const x = mid + half * p;
+                add(unitEquivalentLoad(load.axis, x, L), (0, exports.loadIntensity)(load, x) * half * GAUSS_WEIGHTS[index]);
+            });
+        }
+    }
+    return equivalent.map((v) => -v);
+};
+exports.fixedEndForces = fixedEndForces;
+//# sourceMappingURL=fixed-end-forces.js.map