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

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

@@ -0,0 +1,267 @@
+import { AREA_GEOMETRY_TOLERANCE } from "../../../core/constants/frame.constant.js";
+import { throwIssues, validateFrameModel } from "../validation/validate-model.js";
+import { areaPlane, clipPolygon, edgeVector, resolveAreaSupports, } from "./area-geometry.js";
+import { edgeMemberChain, memberAdjacency } from "./edge-members.js";
+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 = 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 = clipPolygon(region, [ni[0] - nj[0], ni[1] - nj[1]], dot2(ni, corners[i]) - dot2(nj, cj));
+    });
+    const e = 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 = edgeVector(corners, i);
+    const L = Math.hypot(...e);
+    const sin = Math.abs(e[0] * span[1] - e[1] * span[0]) / L;
+    if (sin <= 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)) <= 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.
+ */
+export 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 = memberAdjacency(model.members);
+    const prepared = new Map();
+    for (const area of areas) {
+        if (!loadedAreas.has(area.id))
+            continue;
+        const ids = area.nodeIds;
+        const plane = areaPlane(ids.map(position));
+        const supports = resolveAreaSupports(plane, area.distribution, area.supportEdges);
+        const tolerance = AREA_GEOMETRY_TOLERANCE * plane.size;
+        const n = ids.length;
+        const edges = [];
+        for (const i of supports.loadedEdges) {
+            const chain = edgeMemberChain(ids[i], ids[(i + 1) % n], position, adjacency);
+            const L = Math.hypot(...edgeVector(plane.corners, i));
+            if (area.distribution === "cantilever") {
+                const depth = Math.hypot(...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]) > 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 = 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) > 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 };
+};
+/**
+ * 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`
+ */
+export const distributeAreaLoads = (model) => {
+    const issues = validateFrameModel(model);
+    if (issues.length > 0)
+        throwIssues(issues);
+    return areaMemberLoads(model);
+};
+//# sourceMappingURL=distribute-area-loads.js.map

package/dist/esm/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,OAAO,EAAE,uBAAuB,EAAE,MAAM,oCAAoC,CAAC;AAI7E,OAAO,EAAE,WAAW,EAAE,kBAAkB,EAAE,MAAM,8BAA8B,CAAC;AAC/E,OAAO,EACL,SAAS,EACT,WAAW,EACX,UAAU,EACV,mBAAmB,GAIpB,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,eAAe,EAAE,eAAe,EAAoB,MAAM,gBAAgB,CAAC;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,UAAU,CAAC,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,WAAW,CAClB,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,UAAU,CAAC,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,UAAU,CAAC,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,uBAAuB;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,uBAAuB,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;AACH,MAAM,CAAC,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,eAAe,CAAC,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,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAc,CAAC;QACxD,MAAM,QAAQ,GAAG,mBAAmB,CAAC,KAAK,EAAE,IAAI,CAAC,YAAY,EAAE,IAAI,CAAC,YAAY,CAAE,CAAC;QACnF,MAAM,SAAS,GAAG,uBAAuB,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,eAAe,CAAC,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,UAAU,CAAC,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,UAAU,CAAC,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,uBAAuB,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,uBAAuB,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,uBAAuB,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;AAEF;;;;;;;;;;;;;GAaG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAG,CACjC,KAAiB,EACkC,EAAE;IACrD,MAAM,MAAM,GAAG,kBAAkB,CAAC,KAAK,CAAC,CAAC;IACzC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC;QAAE,WAAW,CAAC,MAAM,CAAC,CAAC;IAC3C,OAAO,eAAe,CAAC,KAAK,CAAC,CAAC;AAChC,CAAC,CAAC"}

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

@@ -0,0 +1,62 @@
+import { AREA_GEOMETRY_TOLERANCE } from "../../../core/constants/frame.constant.js";
+/** Indexes members by the nodes they connect. */
+export 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;
+};
+/**
+ * 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 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 = 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;
+};
+//# sourceMappingURL=edge-members.js.map

package/dist/esm/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,OAAO,EAAE,uBAAuB,EAAE,MAAM,oCAAoC,CAAC;AAgB7E,iDAAiD;AACjD,MAAM,CAAC,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;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,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,uBAAuB,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"}

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

@@ -1,5 +1,7 @@
 import { FrameAnalysisError, FrameErrors } from "../../../core/errors/frame-analysis.error.js";
-// 4-point Gauss–Legendre rule, exact for the cubic integrands used here
+/** Intensity of a distributed load at `x` (inside its extent), in N/mm. */
+export const loadIntensity = (load, x) => load.w1 + ((load.w2 - load.w1) * (x - load.start)) / (load.end - load.start);
+// 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,
 ];
@@ -54,7 +56,17 @@ export const fixedEndForces = (loads, L) => {
             equivalent[i] += unit[i] * factor;
     };
     for (const load of loads) {
-        if (load.type === "point") {
+        if (load.type === "torque") {
+            if (!(load.start >= 0 && load.end <= L && load.start < load.end)) {
+                throw new FrameAnalysisError(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 FrameAnalysisError(FrameErrors.LOAD_OUTSIDE_MEMBER);
             add(unitEquivalentLoad(load.axis, load.a, L), load.P);
@@ -66,7 +78,8 @@ export const fixedEndForces = (loads, L) => {
             const half = (load.end - load.start) / 2;
             const mid = (load.end + load.start) / 2;
             GAUSS_POINTS.forEach((p, index) => {
-                add(unitEquivalentLoad(load.axis, mid + half * p, L), load.w * half * GAUSS_WEIGHTS[index]);
+                const x = mid + half * p;
+                add(unitEquivalentLoad(load.axis, x, L), loadIntensity(load, x) * half * GAUSS_WEIGHTS[index]);
             });
         }
     }

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

@@ -1 +1 @@
-{"version":3,"file":"fixed-end-forces.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/fixed-end-forces.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AAOxF,wEAAwE;AACxE,MAAM,YAAY,GAAG;IACnB,CAAC,kBAAkB,EAAE,CAAC,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CACjF,CAAC;AACF,MAAM,aAAa,GAAG;IACpB,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CAC/E,CAAC;AAEF;;;;GAIG;AACH,MAAM,kBAAkB,GAAG,CAAC,IAAe,EAAE,CAAS,EAAE,CAAS,EAAY,EAAE;IAC7E,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QAChB,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,OAAO,GAAG,CAAC;IACb,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IAC5C,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACrC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IACtC,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;IACf,CAAC;SAAM,CAAC;QACN,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QACb,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC;IAChB,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,KAAwB,EAAE,CAAS,EAAY,EAAE;IAC9E,MAAM,UAAU,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,GAAG,GAAG,CAAC,IAAc,EAAE,MAAc,EAAE,EAAE;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,UAAU,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;IACjE,CAAC,CAAC;IAEF,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;YAC1B,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;gBAC/B,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,GAAG,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACjE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,CAAC;YACD,MAAM,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACxC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE;gBAChC,GAAG,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,GAAG,GAAG,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;YAC9F,CAAC,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,OAAO,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AACnC,CAAC,CAAC"}
+{"version":3,"file":"fixed-end-forces.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/fixed-end-forces.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AAYxF,2EAA2E;AAC3E,MAAM,CAAC,MAAM,aAAa,GAAG,CAC3B,IAAuD,EACvD,CAAS,EACD,EAAE,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;AAE1F,0FAA0F;AAC1F,MAAM,YAAY,GAAG;IACnB,CAAC,kBAAkB,EAAE,CAAC,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CACjF,CAAC;AACF,MAAM,aAAa,GAAG;IACpB,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CAC/E,CAAC;AAEF;;;;GAIG;AACH,MAAM,kBAAkB,GAAG,CAAC,IAAe,EAAE,CAAS,EAAE,CAAS,EAAY,EAAE;IAC7E,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QAChB,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,OAAO,GAAG,CAAC;IACb,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IAC5C,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACrC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IACtC,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;IACf,CAAC;SAAM,CAAC;QACN,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QACb,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC;IAChB,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,KAAwB,EAAE,CAAS,EAAY,EAAE;IAC9E,MAAM,UAAU,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,GAAG,GAAG,CAAC,IAAc,EAAE,MAAc,EAAE,EAAE;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,UAAU,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;IACjE,CAAC,CAAC;IAEF,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,CAAC,IAAI,KAAK,QAAQ,EAAE,CAAC;YAC3B,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACjE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,CAAC;YACD,2EAA2E;YAC3E,MAAM,SAAS,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;YACnD,MAAM,QAAQ,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC7C,UAAU,CAAC,CAAC,CAAC,IAAI,SAAS,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC;YAChD,UAAU,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;QAC9C,CAAC;aAAM,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;YACjC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;gBAC/B,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,GAAG,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACjE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,CAAC;YACD,MAAM,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACxC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE;gBAChC,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,GAAG,CAAC,CAAC;gBACzB,GAAG,CACD,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,EACnC,aAAa,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,aAAa,CAAC,KAAK,CAAC,CACrD,CAAC;YACJ,CAAC,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,OAAO,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AACnC,CAAC,CAAC"}

package/dist/esm/fea/frame/index.js

@@ -1,5 +1,6 @@
 export { analyzeFrame, createFrameAnalysis } from "./solver/analyze-frame.js";
 export { validateFrameModel } from "./validation/validate-model.js";
+export { distributeAreaLoads } from "./area/distribute-area-loads.js";
 export { beam3dLocalStiffness, beam3dLocalStiffnessRightHanded, } from "./elements/beam-3d-stiffness.js";
 export { memberAxes } from "./transformation/coordinate-transform.js";
 export { memberDemandEnvelope } from "./post/demand-envelope.js";

package/dist/esm/fea/frame/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../src/fea/frame/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,YAAY,EAAE,mBAAmB,EAAE,MAAM,wBAAwB,CAAC;AAE3E,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,OAAO,EACL,oBAAoB,EACpB,+BAA+B,GAChC,MAAM,8BAA8B,CAAC;AACtC,OAAO,EAAE,UAAU,EAAE,MAAM,uCAAuC,CAAC;AACnE,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../src/fea/frame/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,YAAY,EAAE,mBAAmB,EAAE,MAAM,wBAAwB,CAAC;AAE3E,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,OAAO,EAAE,mBAAmB,EAAE,MAAM,8BAA8B,CAAC;AACnE,OAAO,EACL,oBAAoB,EACpB,+BAA+B,GAChC,MAAM,8BAA8B,CAAC;AACtC,OAAO,EAAE,UAAU,EAAE,MAAM,uCAAuC,CAAC;AACnE,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC"}

package/dist/esm/fea/frame/post/member-diagram.js

@@ -7,6 +7,9 @@ const DEFLECTION_GRID_INTERVALS = 40;
 /** Golden-section steps refining each deflection extreme between its neighbouring grid points. */
 const REFINE_STEPS = 30;
 const GOLDEN = (Math.sqrt(5) - 1) / 2;
+// 3-point Gauss–Legendre rule on [−1, 1], exact up to degree 5
+const GAUSS3_POINTS = [-Math.sqrt(3 / 5), 0, Math.sqrt(3 / 5)];
+const GAUSS3_WEIGHTS = [5 / 9, 8 / 9, 5 / 9];
 const FORCE_DECIMALS = 4;
 const DISPLACEMENT_DECIMALS = 5;
 /**
@@ -18,11 +21,17 @@ const DISPLACEMENT_DECIMALS = 5;
  */
 const internalForces = (startForce, loads, x, side) => {
     const [f1, f2, f3, m1, m2, m3] = startForce;
-    // Q[axis]: resultant of loads in [0, x]; I[axis]: ∫(s − x)·q ds over [0, x]
+    // Q[axis]: resultant of loads in [0, x]; I[axis]: ∫(s − x)·q ds over [0, x]; Qt: torque in [0, x]
     const Q = [0, 0, 0, 0];
     const I = [0, 0, 0, 0];
+    let Qt = 0;
     for (const load of loads) {
-        if (load.type === "point") {
+        if (load.type === "torque") {
+            const m = Math.min(x, load.end);
+            if (m > load.start)
+                Qt += load.t * (m - load.start);
+        }
+        else if (load.type === "point") {
             const included = side === "left" ? load.a < x : load.a <= x;
             if (included) {
                 Q[load.axis] += load.P;
@@ -32,8 +41,12 @@ const internalForces = (startForce, loads, x, side) => {
         else {
             const m = Math.min(x, load.end);
             if (m > load.start) {
-                Q[load.axis] += load.w * (m - load.start);
-                I[load.axis] += (load.w * ((m - x) ** 2 - (load.start - x) ** 2)) / 2;
+                // q(s) = w1 + k·(s − start); in u = s − x over [a, b]: q = w1 + k·(u − a)
+                const k = (load.w2 - load.w1) / (load.end - load.start);
+                const a = load.start - x;
+                const b = m - x;
+                Q[load.axis] += load.w1 * (b - a) + (k * (b - a) ** 2) / 2;
+                I[load.axis] += ((load.w1 - k * a) * (b ** 2 - a ** 2)) / 2 + (k * (b ** 3 - a ** 3)) / 3;
             }
         }
     }
@@ -41,7 +54,7 @@ const internalForces = (startForce, loads, x, side) => {
         N: -f1 - Q[1],
         V2: -f2 - Q[2],
         V3: -f3 - Q[3],
-        T: -m1,
+        T: -m1 - Qt,
         M2: -m2 - x * f3 + I[3],
         M3: -m3 + x * f2 - I[2],
     };
@@ -52,8 +65,8 @@ const internalForces = (startForce, loads, x, side) => {
  *
  * Each field is `a + b·x + x·∫₀ˣκ − ∫₀ˣ s·κ`, with the running integrals accumulated once along
  * the sorted grid. The grid contains every load breakpoint, so each interval lies inside one
- * smooth piece where N is at most linear and M at most quadratic, and one Simpson panel per
- * interval integrates the (at most cubic) integrands exactly.
+ * smooth piece where N is at most quadratic and M at most cubic, and one 3-point Gauss panel per
+ * interval integrates the (at most quartic) integrands exactly.
  *
  * The constants come from the member's own end displacements: end translations first, an end
  * rotation only where a translation is released. A released rotation (hinge) never enters, so
@@ -64,22 +77,29 @@ const displacementsOnGrid = (member, startForce, endDisplacements, loads, grid)
     const L = member.length;
     const n = grid.length;
     // integrands: [N/EA, M3/EI33, s·M3/EI33, −M2/EI22, −s·M2/EI22]
-    const integrands = (s, side) => {
-        const f = internalForces(startForce, loads, s, side);
+    const integrands = (s) => {
+        const f = internalForces(startForce, loads, s, "right");
         const k3 = f.M3 / (E * I33);
         const k2 = -f.M2 / (E * I22);
         return [f.N / (E * A), k3, s * k3, k2, s * k2];
     };
+    /** ∫ₐᵇ of every integrand: one 3-point Gauss panel, interior points only (no load side). */
+    const panel = (a, b) => {
+        const half = (b - a) / 2;
+        const mid = (a + b) / 2;
+        const out = [0, 0, 0, 0, 0];
+        GAUSS3_POINTS.forEach((p, k) => {
+            const g = integrands(mid + half * p);
+            for (let c = 0; c < 5; c++)
+                out[c] += half * GAUSS3_WEIGHTS[k] * g[c];
+        });
+        return out;
+    };
     const cumulative = Array.from({ length: 5 }, () => new Float64Array(n));
     for (let i = 1; i < n; i++) {
-        const a = grid[i - 1];
-        const b = grid[i];
-        const ga = integrands(a, "right");
-        const gm = integrands((a + b) / 2, "right");
-        const gb = integrands(b, "left");
-        for (let c = 0; c < 5; c++) {
-            cumulative[c][i] = cumulative[c][i - 1] + ((b - a) / 6) * (ga[c] + 4 * gm[c] + gb[c]);
-        }
+        const g = panel(grid[i - 1], grid[i]);
+        for (let c = 0; c < 5; c++)
+            cumulative[c][i] = cumulative[c][i - 1] + g[c];
     }
     const [strain, k3, sk3, k2, sk2] = cumulative;
     const last = n - 1;
@@ -121,7 +141,7 @@ const displacementsOnGrid = (member, startForce, endDisplacements, loads, grid)
         u2: plane2.a + plane2.b * x + P3(i),
         u3: plane3.a + plane3.b * x + P2(i),
     }));
-    /** Displacement anywhere: the grid value plus one exact Simpson panel to x. */
+    /** Displacement anywhere: the grid value plus one exact Gauss panel to x. */
     const at = (x) => {
         let lo = 0;
         let hi = last;
@@ -135,10 +155,8 @@ const displacementsOnGrid = (member, startForce, endDisplacements, loads, grid)
         const a = grid[lo];
         if (x - a <= 0)
             return onGrid[lo];
-        const ga = integrands(a, "right");
-        const gm = integrands((a + x) / 2, "right");
-        const gx = integrands(x, "left");
-        const c = cumulative.map((col, k) => col[lo] + ((x - a) / 6) * (ga[k] + 4 * gm[k] + gx[k]));
+        const g = panel(a, x);
+        const c = cumulative.map((col, k) => col[lo] + g[k]);
         return {
             u1: axialStart + c[0],
             u2: plane2.a + plane2.b * x + (x * c[1] - c[2]),
@@ -194,16 +212,40 @@ export const memberDiagram = (member, startForce, endDisplacements, loads, stati
                 out.push(v);
         return out;
     };
-    // Zero-shear points: V is linear inside each piece between breakpoints, so the root is exact
+    // Inside each piece between breakpoints N and V are at most quadratic (linear loads), so the
+    // quadratic through both ends and the midpoint is exact. Its roots (zero shear: the moment
+    // extremes) and its vertex (zero load: the N and V extremes) are added to the grid.
     const roots = [];
     for (let i = 0; i < sortedBreaks.length - 1; i++) {
         const xa = sortedBreaks[i];
         const xb = sortedBreaks[i + 1];
         const fa = internalForces(startForce, loads, xa, "right");
+        const fm = internalForces(startForce, loads, (xa + xb) / 2, "right");
         const fb = internalForces(startForce, loads, xb, "left");
-        for (const key of ["V2", "V3"]) {
-            if (fa[key] * fb[key] < 0)
-                roots.push(xa + ((xb - xa) * fa[key]) / (fa[key] - fb[key]));
+        const inside = (t) => t > 0 && t < 1 && roots.push(xa + (xb - xa) * t);
+        for (const key of ["N", "V2", "V3"]) {
+            // f(t) = A·t² + B·t + C on t ∈ [0, 1]
+            const C = fa[key];
+            const A = 2 * fa[key] - 4 * fm[key] + 2 * fb[key];
+            const B = 4 * fm[key] - 3 * fa[key] - fb[key];
+            const scale = Math.max(Math.abs(fa[key]), Math.abs(fm[key]), Math.abs(fb[key]));
+            const curved = Math.abs(A) > 1e-12 * scale;
+            if (curved)
+                inside(-B / (2 * A));
+            if (key === "N")
+                continue;
+            if (!curved) {
+                if (Math.abs(B) > 1e-12 * scale)
+                    inside(-C / B);
+            }
+            else {
+                const disc = B * B - 4 * A * C;
+                if (disc >= 0) {
+                    const sq = Math.sqrt(disc);
+                    inside((-B - sq) / (2 * A));
+                    inside((-B + sq) / (2 * A));
+                }
+            }
         }
     }
     const coarse = Array.from({ length: DEFLECTION_GRID_INTERVALS + 1 }, (_, i) => (L * i) / DEFLECTION_GRID_INTERVALS);
@@ -258,8 +300,9 @@ export const memberDiagram = (member, startForce, endDisplacements, loads, stati
         }
     }
     // ---- Extremes -------------------------------------------------------------------
-    // N, V and T are linear and M quadratic inside each piece, so their extremes lie at a
-    // breakpoint (either side) or a zero-shear point — all on the grid, so these are exact.
+    // N and V are at most quadratic, T linear and M at most cubic inside each piece, so their
+    // extremes lie at a breakpoint (either side), a zero-load or a zero-shear point — all on the
+    // grid, so these are exact.
     // Deflections are taken on the grid, then refined by golden-section search.
     const extremes = {};
     const best = {};