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

package/dist/fea/frame/elements/static-condensation.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"static-condensation.js","sourceRoot":"","sources":["../../../../src/fea/frame/elements/static-condensation.ts"],"names":[],"mappings":";;;AAAA,gFAAwF;AACxF,8CAAmE;AAQnE;;;;;;;;;;;GAWG;AACI,MAAM,cAAc,GAAG,CAAC,CAAS,EAAE,QAAkB,EAAmB,EAAE;IAC/E,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC1B,OAAO,EAAE,CAAC,EAAE,sBAAsB,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAC9D,CAAC;IACD,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9F,MAAM,GAAG,GAAG,CAAC,IAAc,EAAE,IAAc,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE1F,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,yBAAyB;IACzB,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/F,MAAM,MAAM,GAAG,IAAA,0BAAiB,EAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC3C,IAAI,CAAC,MAAM;QAAE,MAAM,IAAI,yCAAkB,CAAC,kCAAW,CAAC,eAAe,CAAC,CAAC;IACvE,MAAM,EAAE,GAAG,QAAQ,CAAC,MAAM,CAAC;IAE3B,MAAM,KAAK,GAAW,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvD,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;QACzB,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;YACzB,IAAI,KAAK,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YACxB,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;gBACxB,KAAK,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACtC,CAAC,CAAC,CAAC;YACH,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtB,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,kDAAkD;IAClD,MAAM,SAAS,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CACvC,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAC/F,CAAC;IAEF,OAAO;QACL,CAAC,EAAE,KAAK;QACR,sBAAsB,EAAE,CAAC,KAAK,EAAE,EAAE;YAChC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC5C,MAAM,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACzC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;gBACzB,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5E,CAAC,CAAC,CAAC;YACH,OAAO,KAAK,CAAC;QACf,CAAC;KACF,CAAC;AACJ,CAAC,CAAC;AA5CW,QAAA,cAAc,kBA4CzB"}

package/dist/fea/frame/index.d.ts

@@ -0,0 +1,9 @@
+export { analyzeFrame, createFrameAnalysis } from "./solver/analyze-frame";
+export type { AnalyzeFrameOptions, FrameAnalysis } from "./solver/analyze-frame";
+export { validateFrameModel } from "./validation/validate-model";
+export { beam3dLocalStiffness, beam3dLocalStiffnessRightHanded, } from "./elements/beam-3d-stiffness";
+export { memberAxes } from "./transformation/coordinate-transform";
+export { memberDemandEnvelope } from "./post/demand-envelope";
+export type { MemberDemandEnvelope, DemandEnvelopeOptions } from "./post/demand-envelope";
+export type { FrameId, FrameNode, DofRestraint, Restraint, MemberEndRelease, Member, NodalLoad, LoadDirection, UniformMemberLoad, PointMemberLoad, MemberLoad, LoadCaseType, LoadCase, CombinationTerm, LoadCombination, FrameModel, } from "./types/model.type";
+export type { NodalDisplacement, NodalReaction, MemberEndForce, MemberEndForces, MemberStation, StationComponent, ExtremeValue, MemberDiagram, GoverningValue, GoverningExtremeValue, Bound, FrameCaseResult, FrameEnvelope, FrameResult, } from "./types/results.type";

package/dist/fea/frame/index.js

@@ -0,0 +1,16 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.memberDemandEnvelope = exports.memberAxes = exports.beam3dLocalStiffnessRightHanded = exports.beam3dLocalStiffness = exports.validateFrameModel = exports.createFrameAnalysis = exports.analyzeFrame = void 0;
+var analyze_frame_1 = require("./solver/analyze-frame");
+Object.defineProperty(exports, "analyzeFrame", { enumerable: true, get: function () { return analyze_frame_1.analyzeFrame; } });
+Object.defineProperty(exports, "createFrameAnalysis", { enumerable: true, get: function () { return analyze_frame_1.createFrameAnalysis; } });
+var validate_model_1 = require("./validation/validate-model");
+Object.defineProperty(exports, "validateFrameModel", { enumerable: true, get: function () { return validate_model_1.validateFrameModel; } });
+var beam_3d_stiffness_1 = require("./elements/beam-3d-stiffness");
+Object.defineProperty(exports, "beam3dLocalStiffness", { enumerable: true, get: function () { return beam_3d_stiffness_1.beam3dLocalStiffness; } });
+Object.defineProperty(exports, "beam3dLocalStiffnessRightHanded", { enumerable: true, get: function () { return beam_3d_stiffness_1.beam3dLocalStiffnessRightHanded; } });
+var coordinate_transform_1 = require("./transformation/coordinate-transform");
+Object.defineProperty(exports, "memberAxes", { enumerable: true, get: function () { return coordinate_transform_1.memberAxes; } });
+var demand_envelope_1 = require("./post/demand-envelope");
+Object.defineProperty(exports, "memberDemandEnvelope", { enumerable: true, get: function () { return demand_envelope_1.memberDemandEnvelope; } });
+//# sourceMappingURL=index.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/fea/frame/index.ts"],"names":[],"mappings":";;;AAAA,wDAA2E;AAAlE,6GAAA,YAAY,OAAA;AAAE,oHAAA,mBAAmB,OAAA;AAE1C,8DAAiE;AAAxD,oHAAA,kBAAkB,OAAA;AAC3B,kEAGsC;AAFpC,yHAAA,oBAAoB,OAAA;AACpB,oIAAA,+BAA+B,OAAA;AAEjC,8EAAmE;AAA1D,kHAAA,UAAU,OAAA;AACnB,0DAA8D;AAArD,uHAAA,oBAAoB,OAAA"}

package/dist/fea/frame/post/demand-envelope.d.ts

@@ -0,0 +1,36 @@
+import type { FrameId } from "../types/model.type";
+import type { FrameResult } from "../types/results.type";
+export interface MemberDemandEnvelope {
+    memberId: FrameId;
+    MuPositive: number;
+    MuNegative: number;
+    Vu: number;
+    Nu: number;
+    VuRange: {
+        from: number;
+        to: number;
+    };
+    /** Load case or combination id producing each demand; null when that demand is zero. */
+    governing: Record<"MuPositive" | "MuNegative" | "Vu" | "Nu", string | null>;
+}
+export interface DemandEnvelopeOptions {
+    /** Combination ids to envelope. Default: every combination, unless `cases` is given. */
+    combinations?: string[];
+    /** Load case ids to envelope as well (or instead). Default: none. */
+    cases?: string[];
+    /** mm. Take Vu only between this distance from each end, e.g. the effective depth d. Default 0. */
+    shearCriticalDistance?: number;
+}
+/**
+ * Reduces a member's force diagrams to the demands used in beam design — Mu⁺, Mu⁻, Vu, Nu —
+ * enveloped over load combinations (and optionally load cases), with the governing id of each.
+ * Works on plain numbers and does not depend on the `rc` module.
+ *
+ * @param result - Result of `analyzeFrame`
+ * @param memberId - Member to reduce
+ * @param options - Combinations and/or cases to envelope, `shearCriticalDistance` in mm
+ * @returns Demands in kN and kN·m with their governing load case or combination
+ * @throws {Error} if the member or a requested id is not in the result, if nothing is selected,
+ *   or if `shearCriticalDistance` is used on a result analysed with `diagrams: "extremes"`
+ */
+export declare const memberDemandEnvelope: (result: FrameResult, memberId: FrameId, options?: DemandEnvelopeOptions) => MemberDemandEnvelope;

package/dist/fea/frame/post/demand-envelope.js

@@ -0,0 +1,101 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.memberDemandEnvelope = void 0;
+const math_1 = require("../../../utils/math");
+/** Demands of one diagram; Vu between `from` and `to`. */
+const diagramDemands = (diagram, from, to) => {
+    const { N, M3, V2 } = diagram.extremes;
+    const moments = {
+        MuPositive: Math.max(M3.max.value, 0),
+        MuNegative: Math.max(-M3.min.value, 0),
+        Nu: Math.abs(N.max.value) >= Math.abs(N.min.value) ? N.max.value : N.min.value,
+    };
+    if (diagram.stations.length === 0) {
+        // analysed with diagrams: "extremes" — Vu over the whole length comes from the extremes
+        if (from > 0) {
+            throw new Error('shearCriticalDistance needs member stations: analyse with diagrams: "full", or pass a diagram from createFrameAnalysis(...).memberDiagram.');
+        }
+        return { ...moments, Vu: Math.max(Math.abs(V2.max.value), Math.abs(V2.min.value)) };
+    }
+    // V2 is piecewise linear with stations at every breakpoint, so interpolation is exact.
+    // Stations at a point load come as a (left, right) pair sharing the same x.
+    const v2At = (x, side) => {
+        const stations = diagram.stations;
+        const hits = stations.filter((s) => s.x === x);
+        if (hits.length > 0)
+            return side === "left" ? hits[0].V2 : hits[hits.length - 1].V2;
+        for (let i = 0; i < stations.length - 1; i++) {
+            const a = stations[i];
+            const b = stations[i + 1];
+            if (a.x < x && x < b.x)
+                return a.V2 + ((b.V2 - a.V2) * (x - a.x)) / (b.x - a.x);
+        }
+        return stations[stations.length - 1].V2;
+    };
+    const candidates = [v2At(from, "right"), v2At(to, "left")];
+    for (const s of diagram.stations) {
+        if (s.x >= from && s.x <= to)
+            candidates.push(s.V2);
+    }
+    return { ...moments, Vu: (0, math_1.roundToDecimalPlaces)(Math.max(...candidates.map(Math.abs)), 4) + 0 };
+};
+/**
+ * Reduces a member's force diagrams to the demands used in beam design — Mu⁺, Mu⁻, Vu, Nu —
+ * enveloped over load combinations (and optionally load cases), with the governing id of each.
+ * Works on plain numbers and does not depend on the `rc` module.
+ *
+ * @param result - Result of `analyzeFrame`
+ * @param memberId - Member to reduce
+ * @param options - Combinations and/or cases to envelope, `shearCriticalDistance` in mm
+ * @returns Demands in kN and kN·m with their governing load case or combination
+ * @throws {Error} if the member or a requested id is not in the result, if nothing is selected,
+ *   or if `shearCriticalDistance` is used on a result analysed with `diagrams: "extremes"`
+ */
+const memberDemandEnvelope = (result, memberId, options = {}) => {
+    const pick = (source, ids, label) => ids.map((id) => {
+        const r = source[id];
+        if (!r)
+            throw new Error(`${label} "${id}" is not in the frame result.`);
+        return r;
+    });
+    const combinationIds = options.combinations ?? (options.cases ? [] : Object.keys(result.combinations));
+    const selected = [
+        ...pick(result.combinations, combinationIds, "Load combination"),
+        ...pick(result.cases, options.cases ?? [], "Load case"),
+    ];
+    if (selected.length === 0) {
+        throw new Error("No load combinations or cases selected; the result has no combinations, so pass `cases`.");
+    }
+    const diagrams = selected.map((r) => {
+        const diagram = r.memberDiagrams.find((d) => d.memberId === memberId);
+        if (!diagram)
+            throw new Error(`Member ${String(memberId)} is not in the frame result.`);
+        return { id: r.id, diagram };
+    });
+    const L = diagrams[0].diagram.length;
+    const from = Math.min(Math.max(options.shearCriticalDistance ?? 0, 0), L / 2);
+    const to = L - from;
+    const envelope = {
+        MuPositive: 0,
+        MuNegative: 0,
+        Vu: 0,
+        Nu: 0,
+        governing: { MuPositive: null, MuNegative: null, Vu: null, Nu: null },
+    };
+    for (const { id, diagram } of diagrams) {
+        const d = diagramDemands(diagram, from, to);
+        for (const key of ["MuPositive", "MuNegative", "Vu"]) {
+            if (d[key] > envelope[key]) {
+                envelope[key] = d[key];
+                envelope.governing[key] = id;
+            }
+        }
+        if (Math.abs(d.Nu) > Math.abs(envelope.Nu)) {
+            envelope.Nu = d.Nu;
+            envelope.governing.Nu = id;
+        }
+    }
+    return { memberId, ...envelope, VuRange: { from, to } };
+};
+exports.memberDemandEnvelope = memberDemandEnvelope;
+//# sourceMappingURL=demand-envelope.js.map

package/dist/fea/frame/post/demand-envelope.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"demand-envelope.js","sourceRoot":"","sources":["../../../../src/fea/frame/post/demand-envelope.ts"],"names":[],"mappings":";;;AAAA,0CAAuD;AAyBvD,0DAA0D;AAC1D,MAAM,cAAc,GAAG,CAAC,OAAsB,EAAE,IAAY,EAAE,EAAU,EAAE,EAAE;IAC1E,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,OAAO,CAAC,QAAQ,CAAC;IACvC,MAAM,OAAO,GAAG;QACd,UAAU,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC;QACrC,UAAU,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC;QACtC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK;KAC/E,CAAC;IACF,IAAI,OAAO,CAAC,QAAQ,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAClC,wFAAwF;QACxF,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;YACb,MAAM,IAAI,KAAK,CACb,4IAA4I,CAC7I,CAAC;QACJ,CAAC;QACD,OAAO,EAAE,GAAG,OAAO,EAAE,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;IACtF,CAAC;IAED,uFAAuF;IACvF,4EAA4E;IAC5E,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,IAAsB,EAAU,EAAE;QACzD,MAAM,QAAQ,GAAG,OAAO,CAAC,QAAQ,CAAC;QAClC,MAAM,IAAI,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;QAC/C,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC;YAAE,OAAO,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC;QACpF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7C,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;YACtB,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;gBAAE,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAClF,CAAC;QACD,OAAO,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC;IAC1C,CAAC,CAAC;IAEF,MAAM,UAAU,GAAa,CAAC,IAAI,CAAC,IAAI,EAAE,OAAO,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,MAAM,CAAC,CAAC,CAAC;IACrE,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,QAAQ,EAAE,CAAC;QACjC,IAAI,CAAC,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC,CAAC,CAAC,IAAI,EAAE;YAAE,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACtD,CAAC;IACD,OAAO,EAAE,GAAG,OAAO,EAAE,EAAE,EAAE,IAAA,2BAAoB,EAAC,IAAI,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;AAChG,CAAC,CAAC;AAEF;;;;;;;;;;;GAWG;AACI,MAAM,oBAAoB,GAAG,CAClC,MAAmB,EACnB,QAAiB,EACjB,UAAiC,EAAE,EACb,EAAE;IACxB,MAAM,IAAI,GAAG,CAAC,MAAuC,EAAE,GAAa,EAAE,KAAa,EAAE,EAAE,CACrF,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE;QACb,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC;QACrB,IAAI,CAAC,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,GAAG,KAAK,KAAK,EAAE,+BAA+B,CAAC,CAAC;QACxE,OAAO,CAAC,CAAC;IACX,CAAC,CAAC,CAAC;IACL,MAAM,cAAc,GAClB,OAAO,CAAC,YAAY,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAC;IAClF,MAAM,QAAQ,GAAG;QACf,GAAG,IAAI,CAAC,MAAM,CAAC,YAAY,EAAE,cAAc,EAAE,kBAAkB,CAAC;QAChE,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,OAAO,CAAC,KAAK,IAAI,EAAE,EAAE,WAAW,CAAC;KACxD,CAAC;IACF,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC1B,MAAM,IAAI,KAAK,CACb,0FAA0F,CAC3F,CAAC;IACJ,CAAC;IAED,MAAM,QAAQ,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;QAClC,MAAM,OAAO,GAAG,CAAC,CAAC,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,QAAQ,CAAC,CAAC;QACtE,IAAI,CAAC,OAAO;YAAE,MAAM,IAAI,KAAK,CAAC,UAAU,MAAM,CAAC,QAAQ,CAAC,8BAA8B,CAAC,CAAC;QACxF,OAAO,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,OAAO,EAAE,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC;IACrC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,qBAAqB,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IAC9E,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC;IAEpB,MAAM,QAAQ,GAAuD;QACnE,UAAU,EAAE,CAAC;QACb,UAAU,EAAE,CAAC;QACb,EAAE,EAAE,CAAC;QACL,EAAE,EAAE,CAAC;QACL,SAAS,EAAE,EAAE,UAAU,EAAE,IAAI,EAAE,UAAU,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,EAAE;KACtE,CAAC;IACF,KAAK,MAAM,EAAE,EAAE,EAAE,OAAO,EAAE,IAAI,QAAQ,EAAE,CAAC;QACvC,MAAM,CAAC,GAAG,cAAc,CAAC,OAAO,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC;QAC5C,KAAK,MAAM,GAAG,IAAI,CAAC,YAAY,EAAE,YAAY,EAAE,IAAI,CAAU,EAAE,CAAC;YAC9D,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC;gBAC3B,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBACvB,QAAQ,CAAC,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;YAC/B,CAAC;QACH,CAAC;QACD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC;YAC3C,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;YACnB,QAAQ,CAAC,SAAS,CAAC,EAAE,GAAG,EAAE,CAAC;QAC7B,CAAC;IACH,CAAC;IAED,OAAO,EAAE,QAAQ,EAAE,GAAG,QAAQ,EAAE,OAAO,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,CAAC;AAC1D,CAAC,CAAC;AAvDW,QAAA,oBAAoB,wBAuD/B"}

package/dist/fea/frame/post/envelope.d.ts

@@ -0,0 +1,8 @@
+import type { FrameCaseResult, FrameEnvelope } from "../types/results.type";
+/**
+ * Max/min of every presented result across combinations, each with the combination that
+ * governs it. Combinations must come from the same model (same nodes and members, same order).
+ *
+ * @param results - Combination results, in model order
+ */
+export declare const buildEnvelope: (results: FrameCaseResult[]) => FrameEnvelope;

package/dist/fea/frame/post/envelope.js

@@ -0,0 +1,63 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.buildEnvelope = void 0;
+const DISPLACEMENT_KEYS = ["ux", "uy", "uz", "rx", "ry", "rz"];
+const REACTION_KEYS = ["fx", "fy", "fz", "mx", "my", "mz"];
+const COMPONENTS = ["N", "V2", "V3", "T", "M2", "M3", "u1", "u2", "u3"];
+/** Keeps the larger max and smaller min; the first combination wins a tie. */
+const update = (bound, max, min) => {
+    if (!bound)
+        return { max, min };
+    return {
+        max: max.value > bound.max.value ? max : bound.max,
+        min: min.value < bound.min.value ? min : bound.min,
+    };
+};
+/**
+ * Max/min of every presented result across combinations, each with the combination that
+ * governs it. Combinations must come from the same model (same nodes and members, same order).
+ *
+ * @param results - Combination results, in model order
+ */
+const buildEnvelope = (results) => {
+    const first = results[0];
+    const nodalDisplacements = first.nodalDisplacements.map((d, n) => {
+        const entry = { nodeId: d.nodeId };
+        for (const key of DISPLACEMENT_KEYS) {
+            let bound;
+            for (const r of results) {
+                const g = { value: r.nodalDisplacements[n][key], combinationId: r.id };
+                bound = update(bound, g, g);
+            }
+            entry[key] = bound;
+        }
+        return entry;
+    });
+    const reactions = first.reactions.map((d, n) => {
+        const entry = { nodeId: d.nodeId };
+        for (const key of REACTION_KEYS) {
+            let bound;
+            for (const r of results) {
+                const g = { value: r.reactions[n][key], combinationId: r.id };
+                bound = update(bound, g, g);
+            }
+            entry[key] = bound;
+        }
+        return entry;
+    });
+    const members = first.memberDiagrams.map((d, m) => {
+        const extremes = {};
+        for (const c of COMPONENTS) {
+            let bound;
+            for (const r of results) {
+                const e = r.memberDiagrams[m].extremes[c];
+                bound = update(bound, { ...e.max, combinationId: r.id }, { ...e.min, combinationId: r.id });
+            }
+            extremes[c] = bound;
+        }
+        return { memberId: d.memberId, extremes };
+    });
+    return { combinationIds: results.map((r) => r.id), nodalDisplacements, reactions, members };
+};
+exports.buildEnvelope = buildEnvelope;
+//# sourceMappingURL=envelope.js.map

package/dist/fea/frame/post/envelope.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"envelope.js","sourceRoot":"","sources":["../../../../src/fea/frame/post/envelope.ts"],"names":[],"mappings":";;;AASA,MAAM,iBAAiB,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AACxE,MAAM,aAAa,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AACpE,MAAM,UAAU,GAAuB,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;AAE5F,8EAA8E;AAC9E,MAAM,MAAM,GAAG,CACb,KAA2B,EAC3B,GAAM,EACN,GAAM,EACI,EAAE;IACZ,IAAI,CAAC,KAAK;QAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAChC,OAAO;QACL,GAAG,EAAE,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG;QAClD,GAAG,EAAE,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG;KACnD,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;GAKG;AACI,MAAM,aAAa,GAAG,CAAC,OAA0B,EAAiB,EAAE;IACzE,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IAEzB,MAAM,kBAAkB,GAAG,KAAK,CAAC,kBAAkB,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QAC/D,MAAM,KAAK,GAAG,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,EAAiD,CAAC;QAClF,KAAK,MAAM,GAAG,IAAI,iBAAiB,EAAE,CAAC;YACpC,IAAI,KAAwC,CAAC;YAC7C,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;gBACxB,MAAM,CAAC,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,aAAa,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC;gBACvE,KAAK,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC9B,CAAC;YACD,KAAK,CAAC,GAAG,CAAC,GAAG,KAA8B,CAAC;QAC9C,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC,CAAC,CAAC;IAEH,MAAM,SAAS,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QAC7C,MAAM,KAAK,GAAG,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,EAAwC,CAAC;QACzE,KAAK,MAAM,GAAG,IAAI,aAAa,EAAE,CAAC;YAChC,IAAI,KAAwC,CAAC;YAC7C,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;gBACxB,MAAM,CAAC,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,aAAa,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC;gBAC9D,KAAK,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC9B,CAAC;YACD,KAAK,CAAC,GAAG,CAAC,GAAG,KAA8B,CAAC;QAC9C,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC,CAAC,CAAC;IAEH,MAAM,OAAO,GAAG,KAAK,CAAC,cAAc,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QAChD,MAAM,QAAQ,GAAG,EAA4D,CAAC;QAC9E,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE,CAAC;YAC3B,IAAI,KAA+C,CAAC;YACpD,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;gBACxB,MAAM,CAAC,GAAG,CAAC,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;gBAC1C,KAAK,GAAG,MAAM,CAAC,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,aAAa,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,aAAa,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;YAC9F,CAAC;YACD,QAAQ,CAAC,CAAC,CAAC,GAAG,KAAqC,CAAC;QACtD,CAAC;QACD,OAAO,EAAE,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,QAAQ,EAAE,CAAC;IAC5C,CAAC,CAAC,CAAC;IAEH,OAAO,EAAE,cAAc,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,kBAAkB,EAAE,SAAS,EAAE,OAAO,EAAE,CAAC;AAC9F,CAAC,CAAC;AA3CW,QAAA,aAAa,iBA2CxB"}

package/dist/fea/frame/post/member-diagram.d.ts

@@ -0,0 +1,26 @@
+import type { SectionProperties } from "../../../core/types/section-properties.type";
+import type { Matrix } from "../../../utils/matrix";
+import type { LocalMemberLoad } from "../elements/fixed-end-forces";
+import type { FrameId } from "../types/model.type";
+import type { MemberDiagram } from "../types/results.type";
+/** The member data a diagram needs. */
+export interface DiagramMember {
+    id: FrameId;
+    length: number;
+    section: Pick<SectionProperties, "E" | "A" | "I22" | "I33">;
+    lambda: Matrix;
+    released: boolean[];
+}
+/**
+ * Builds the force and displacement diagrams of one member from its end forces, its end
+ * displacements and its loads.
+ *
+ * @param member - Member geometry, section, axes and end releases
+ * @param startForce - Force the start node exerts on the member, right-hand local, N and N·mm
+ * @param endDisplacements - The 12 local right-hand displacements of the member's end nodes, mm and rad
+ * @param loads - Member loads in local axes, N and mm
+ * @param stationCount - Number of evenly spaced stations (at least 2); load breakpoints and
+ *   zero-shear points are added, with both sides reported at a point load
+ * @param withStations - false returns the extremes only, with an empty `stations` array
+ */
+export declare const memberDiagram: (member: DiagramMember, startForce: number[], endDisplacements: number[], loads: LocalMemberLoad[], stationCount?: number, withStations?: boolean) => MemberDiagram;

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

@@ -0,0 +1,336 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.memberDiagram = void 0;
+const math_1 = require("../../../utils/math");
+const FORCE_COMPONENTS = ["N", "V2", "V3", "T", "M2", "M3"];
+const DISPLACEMENT_COMPONENTS = ["u1", "u2", "u3"];
+const COMPONENTS = [...FORCE_COMPONENTS, ...DISPLACEMENT_COMPONENTS];
+/** Evenly spaced grid for the deflection extremes (force extremes need no grid: see below). */
+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;
+const FORCE_DECIMALS = 4;
+const DISPLACEMENT_DECIMALS = 5;
+/**
+ * Internal forces at distance `x` from the member start, in N and N·mm. Taken from the
+ * free body of the segment [0, x]; N is positive in tension, `M3` is positive for sagging
+ * when local 2 points up. At a point load, `side` picks the value just before or after it.
+ *
+ * @param startForce - Right-hand local force the start node exerts on the member: [F1, F2, F3, M1, M2, M3]
+ */
+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]
+    const Q = [0, 0, 0, 0];
+    const I = [0, 0, 0, 0];
+    for (const load of loads) {
+        if (load.type === "point") {
+            const included = side === "left" ? load.a < x : load.a <= x;
+            if (included) {
+                Q[load.axis] += load.P;
+                I[load.axis] += load.P * (load.a - x);
+            }
+        }
+        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;
+            }
+        }
+    }
+    return {
+        N: -f1 - Q[1],
+        V2: -f2 - Q[2],
+        V3: -f3 - Q[3],
+        T: -m1,
+        M2: -m2 - x * f3 + I[3],
+        M3: -m3 + x * f2 - I[2],
+    };
+};
+/**
+ * Local displacements u1, u2, u3 (mm) at every grid point, from Euler–Bernoulli kinematics:
+ * `u1' = N/EA`, `u2'' = M3/EI33` (R3 = u2'), `u3'' = −M2/EI22` (R2 = −u3').
+ *
+ * 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.
+ *
+ * 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
+ * the curve is correct at hinged ends whose rotation differs from the node's.
+ */
+const displacementsOnGrid = (member, startForce, endDisplacements, loads, grid) => {
+    const { E, A, I22, I33 } = member.section;
+    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 k3 = f.M3 / (E * I33);
+        const k2 = -f.M2 / (E * I22);
+        return [f.N / (E * A), k3, s * k3, k2, s * k2];
+    };
+    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 [strain, k3, sk3, k2, sk2] = cumulative;
+    const last = n - 1;
+    /** a, b of `u = a + b·x + P(x)` from the first independent pair of end conditions. */
+    const planeConstants = (translation, rotation, rotationSign, PL, // P(L)
+    dPL) => {
+        const { released } = member;
+        // rows [ca, cb | rhs] of ca·a + cb·b = rhs
+        const conditions = [];
+        if (!released[translation])
+            conditions.push([1, 0, endDisplacements[translation]]);
+        if (!released[6 + translation]) {
+            conditions.push([1, L, endDisplacements[6 + translation] - PL]);
+        }
+        if (!released[rotation])
+            conditions.push([0, 1, rotationSign * endDisplacements[rotation]]);
+        if (!released[6 + rotation]) {
+            conditions.push([0, 1, rotationSign * endDisplacements[6 + rotation] - dPL]);
+        }
+        for (let i = 0; i < conditions.length; i++) {
+            for (let j = i + 1; j < conditions.length; j++) {
+                const [a1, b1, r1] = conditions[i];
+                const [a2, b2, r2] = conditions[j];
+                const det = a1 * b2 - a2 * b1;
+                if (Math.abs(det) > 1e-12 * Math.max(1, L)) {
+                    return { a: (r1 * b2 - r2 * b1) / det, b: (a1 * r2 - a2 * r1) / det };
+                }
+            }
+        }
+        return { a: 0, b: 0 }; // unreachable: such releases are rejected as a mechanism (307)
+    };
+    const P3 = (i) => grid[i] * k3[i] - sk3[i];
+    const P2 = (i) => grid[i] * k2[i] - sk2[i];
+    const plane2 = planeConstants(1, 5, 1, P3(last), k3[last]); // u2, with u2' = R3
+    const plane3 = planeConstants(2, 4, -1, P2(last), k2[last]); // u3, with u3' = −R2
+    const axialStart = member.released[0] ? endDisplacements[6] - strain[last] : endDisplacements[0];
+    const onGrid = grid.map((x, i) => ({
+        u1: axialStart + strain[i],
+        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. */
+    const at = (x) => {
+        let lo = 0;
+        let hi = last;
+        while (lo < hi) {
+            const mid = (lo + hi + 1) >> 1;
+            if (grid[mid] <= x)
+                lo = mid;
+            else
+                hi = mid - 1;
+        }
+        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]));
+        return {
+            u1: axialStart + c[0],
+            u2: plane2.a + plane2.b * x + (x * c[1] - c[2]),
+            u3: plane3.a + plane3.b * x + (x * c[3] - c[4]),
+        };
+    };
+    return { onGrid, at };
+};
+const SCALE = {
+    N: 1e-3,
+    V2: 1e-3,
+    V3: 1e-3,
+    T: 1e-6,
+    M2: 1e-6,
+    M3: 1e-6,
+    u1: 1,
+    u2: 1,
+    u3: 1,
+};
+const present = (c, value) => (0, math_1.roundToDecimalPlaces)(value * SCALE[c], c === "u1" || c === "u2" || c === "u3" ? DISPLACEMENT_DECIMALS : FORCE_DECIMALS) + 0;
+/**
+ * Builds the force and displacement diagrams of one member from its end forces, its end
+ * displacements and its loads.
+ *
+ * @param member - Member geometry, section, axes and end releases
+ * @param startForce - Force the start node exerts on the member, right-hand local, N and N·mm
+ * @param endDisplacements - The 12 local right-hand displacements of the member's end nodes, mm and rad
+ * @param loads - Member loads in local axes, N and mm
+ * @param stationCount - Number of evenly spaced stations (at least 2); load breakpoints and
+ *   zero-shear points are added, with both sides reported at a point load
+ * @param withStations - false returns the extremes only, with an empty `stations` array
+ */
+const memberDiagram = (member, startForce, endDisplacements, loads, stationCount = 21, withStations = true) => {
+    const L = member.length;
+    const tolerance = 1e-9 * L;
+    const breakpoints = new Set([0, L]);
+    for (const load of loads) {
+        if (load.type === "point")
+            breakpoints.add(load.a);
+        else {
+            breakpoints.add(load.start);
+            breakpoints.add(load.end);
+        }
+    }
+    const sortedBreaks = [...breakpoints].sort((a, b) => a - b);
+    const isBreak = (x) => sortedBreaks.some((b) => Math.abs(b - x) < tolerance);
+    /** Sorted, with points closer than the tolerance merged (the first one kept). */
+    const mergeSorted = (values) => {
+        const sorted = [...values].sort((a, b) => a - b);
+        const out = [];
+        for (const v of sorted)
+            if (out.length === 0 || v - out[out.length - 1] >= tolerance)
+                out.push(v);
+        return out;
+    };
+    // Zero-shear points: V is linear inside each piece between breakpoints, so the root is exact
+    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 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 coarse = Array.from({ length: DEFLECTION_GRID_INTERVALS + 1 }, (_, i) => (L * i) / DEFLECTION_GRID_INTERVALS);
+    const count = Math.max(2, Math.floor(stationCount));
+    const uniform = withStations
+        ? Array.from({ length: count }, (_, i) => (L * i) / (count - 1))
+        : [];
+    // One grid for everything: breakpoints, zero-shear points, stations and the deflection grid
+    const grid = mergeSorted([...sortedBreaks, ...roots, ...uniform, ...coarse]);
+    const field = displacementsOnGrid(member, startForce, endDisplacements, loads, grid);
+    const displacements = field.onGrid;
+    const gridIndex = (x) => {
+        let lo = 0;
+        let hi = grid.length - 1;
+        while (lo < hi) {
+            const mid = (lo + hi) >> 1;
+            if (grid[mid] < x - tolerance)
+                lo = mid + 1;
+            else
+                hi = mid;
+        }
+        return lo;
+    };
+    // ---- Stations -------------------------------------------------------------------
+    const stations = [];
+    if (withStations) {
+        const [l1, l2, l3] = member.lambda;
+        const toStation = (x, f, u) => {
+            const station = { x: (0, math_1.roundToDecimalPlaces)(x, 3) };
+            for (const c of FORCE_COMPONENTS)
+                station[c] = present(c, f[c]);
+            for (const c of DISPLACEMENT_COMPONENTS)
+                station[c] = present(c, u[c]);
+            // global = Σ local component × its axis direction
+            const g = [0, 1, 2].map((k) => u.u1 * l1[k] + u.u2 * l2[k] + u.u3 * l3[k]);
+            station.ux = (0, math_1.roundToDecimalPlaces)(g[0], DISPLACEMENT_DECIMALS) + 0;
+            station.uy = (0, math_1.roundToDecimalPlaces)(g[1], DISPLACEMENT_DECIMALS) + 0;
+            station.uz = (0, math_1.roundToDecimalPlaces)(g[2], DISPLACEMENT_DECIMALS) + 0;
+            return station;
+        };
+        const stationXs = mergeSorted([...uniform, ...sortedBreaks, ...roots]);
+        for (const x of stationXs) {
+            const u = displacements[gridIndex(x)];
+            const right = internalForces(startForce, loads, x, "right");
+            if (isBreak(x) && x > 0 && x < L) {
+                const left = internalForces(startForce, loads, x, "left");
+                if (FORCE_COMPONENTS.some((c) => Math.abs(left[c] - right[c]) > 1e-9)) {
+                    stations.push(toStation(x, left, u));
+                }
+            }
+            stations.push(toStation(x, right, u));
+        }
+    }
+    // ---- 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.
+    // Deflections are taken on the grid, then refined by golden-section search.
+    const extremes = {};
+    const best = {};
+    const consider = (c, value, x, i) => {
+        const b = best[c];
+        if (!b)
+            best[c] = { max: [value, x, i], min: [value, x, i] };
+        else {
+            if (value > b.max[0])
+                b.max = [value, x, i];
+            if (value < b.min[0])
+                b.min = [value, x, i];
+        }
+    };
+    grid.forEach((x, i) => {
+        const sides = isBreak(x) && x > 0 && x < L ? ["right", "left"] : ["right"];
+        for (const side of sides) {
+            const f = internalForces(startForce, loads, x, side);
+            for (const c of FORCE_COMPONENTS)
+                consider(c, f[c], x, i);
+        }
+        for (const c of DISPLACEMENT_COMPONENTS)
+            consider(c, displacements[i][c], x, i);
+    });
+    const refine = (c, i, sign) => {
+        // maximise sign·u between the neighbouring grid points
+        let a = grid[Math.max(0, i - 1)];
+        let b = grid[Math.min(grid.length - 1, i + 1)];
+        const value = (x) => sign * field.at(x)[c];
+        let x1 = b - GOLDEN * (b - a);
+        let x2 = a + GOLDEN * (b - a);
+        let f1 = value(x1);
+        let f2 = value(x2);
+        for (let step = 0; step < REFINE_STEPS; step++) {
+            if (f1 < f2) {
+                a = x1;
+                x1 = x2;
+                f1 = f2;
+                x2 = a + GOLDEN * (b - a);
+                f2 = value(x2);
+            }
+            else {
+                b = x2;
+                x2 = x1;
+                f2 = f1;
+                x1 = b - GOLDEN * (b - a);
+                f1 = value(x1);
+            }
+        }
+        return f1 >= f2 ? [sign * f1, x1] : [sign * f2, x2];
+    };
+    for (const c of DISPLACEMENT_COMPONENTS) {
+        const b = best[c];
+        const [maxValue, maxX] = refine(c, b.max[2], 1);
+        if (maxValue > b.max[0])
+            b.max = [maxValue, maxX, b.max[2]];
+        const [minValue, minX] = refine(c, b.min[2], -1);
+        if (minValue < b.min[0])
+            b.min = [minValue, minX, b.min[2]];
+    }
+    for (const c of COMPONENTS) {
+        const make = ([value, x]) => ({
+            value: present(c, value),
+            x: (0, math_1.roundToDecimalPlaces)(x, 3),
+        });
+        extremes[c] = { max: make(best[c].max), min: make(best[c].min) };
+    }
+    return { memberId: member.id, length: L, stations, extremes };
+};
+exports.memberDiagram = memberDiagram;
+//# sourceMappingURL=member-diagram.js.map