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

package/dist/esm/fea/frame/solver/analyze-frame.js

@@ -0,0 +1,123 @@
+import { buildEnvelope } from "../post/envelope.js";
+import { presentMemberDiagram, presentState } from "../post/present-state.js";
+import { throwIssues, validateFrameModel } from "../validation/validate-model.js";
+import { resolveCombinations } from "./combinations.js";
+import { prepareFrameSystem } from "./frame-system.js";
+import { combineStates, solveLoadCase } from "./load-case-solver.js";
+/** Groups loads by their load case id, keeping the input order within each case. */
+const groupByCase = (loads = []) => {
+    const groups = new Map();
+    for (const load of loads) {
+        const group = groups.get(load.loadCase);
+        if (group)
+            group.push(load);
+        else
+            groups.set(load.loadCase, [load]);
+    }
+    return groups;
+};
+/**
+ * Analyses a 3D frame like `analyzeFrame` and keeps the solution, so that full member diagrams
+ * can be requested later one member at a time — e.g. analyse with `diagrams: "extremes"` and
+ * draw a member's diagram when the user selects it. Run it in a Web Worker for large models.
+ *
+ * @param model - Nodes, members, restraints, load cases, loads tagged with their load case,
+ *   and combinations
+ * @param options - Diagram station count and detail
+ * @returns `{ result, memberDiagram }`
+ * @throws {FrameAnalysisError} 301–314 for an invalid model (with `details` and every problem in
+ *   `issues`), 304 for an unstable structure
+ */
+export const createFrameAnalysis = (model, options = {}) => {
+    const issues = validateFrameModel(model);
+    if (issues.length > 0)
+        throwIssues(issues);
+    const system = prepareFrameSystem(model);
+    const combinations = model.combinations ?? [];
+    const resolved = resolveCombinations(combinations, new Set(model.loadCases.map((c) => c.id)));
+    const withStations = (options.diagrams ?? "full") === "full";
+    const states = new Map();
+    const cases = {};
+    const nodalByCase = groupByCase(model.nodalLoads);
+    const memberByCase = groupByCase(model.memberLoads);
+    for (const loadCase of model.loadCases) {
+        const state = solveLoadCase(system, loadCase, {
+            nodalLoads: nodalByCase.get(loadCase.id) ?? [],
+            memberLoads: memberByCase.get(loadCase.id) ?? [],
+        });
+        states.set(loadCase.id, state);
+        cases[loadCase.id] = presentState(system, state, { id: loadCase.id, kind: "case", name: loadCase.name, type: loadCase.type }, { stations: options.stations, withStations });
+    }
+    const combinationStates = new Map();
+    const combinationResults = {};
+    for (const combination of combinations) {
+        // Non-empty: validation rejects combinations without terms
+        const terms = [...resolved.get(combination.id)].map(([caseId, factor]) => ({ state: states.get(caseId), factor }));
+        const state = combineStates(terms);
+        combinationStates.set(combination.id, state);
+        combinationResults[combination.id] = presentState(system, state, { id: combination.id, kind: "combination", name: combination.name }, { stations: options.stations, withStations });
+    }
+    const warnings = [
+        ...Object.values(cases),
+        ...Object.values(combinationResults),
+    ].flatMap((r) => r.warnings.map((w) => ({
+        ...w,
+        message: `${r.kind === "case" ? "Load case" : "Load combination"} "${r.id}": ${w.message}`,
+    })));
+    const combinationList = combinations.map((c) => combinationResults[c.id]);
+    const result = {
+        cases,
+        combinations: combinationResults,
+        envelope: combinationList.length > 0 ? buildEnvelope(combinationList) : undefined,
+        calculationDetails: {
+            nodeCount: system.nodes.length,
+            memberCount: system.members.length,
+            dofCount: system.dofCount,
+            freeDofCount: system.free.length,
+            profileSize: system.profileSize,
+            loadCaseCount: model.loadCases.length,
+            combinationCount: combinations.length,
+        },
+        units: { force: "kN", moment: "kN·m", length: "mm", rotation: "rad" },
+        warnings,
+    };
+    return {
+        result,
+        memberDiagram: (resultId, memberId, stations) => {
+            const state = combinationStates.get(resultId) ?? states.get(resultId);
+            if (!state)
+                throw new Error(`Load case or combination "${resultId}" is not in the model.`);
+            const m = system.memberIndex.get(memberId);
+            if (m === undefined)
+                throw new Error(`Member ${String(memberId)} is not in the model.`);
+            return presentMemberDiagram(system, state, m, stations ?? options.stations);
+        },
+    };
+};
+/**
+ * Analyses a 3D frame by the direct stiffness method (linear elastic, small displacements)
+ * for every load case, then forms every load combination by superposition.
+ *
+ * The model is validated first (`validateFrameModel`). The stiffness matrix is assembled once
+ * in skyline storage with the nodes in reverse Cuthill–McKee order and factorised once; each
+ * load case is a separate solve. A combination is the factored sum of its load cases (nested
+ * combinations are expanded), and its member diagrams are rebuilt from the combined end forces
+ * and factored member loads, so their extremes are exact rather than sums of case extremes.
+ *
+ * Units: node coordinates, section properties and positions in mm and MPa; loads in kN, kN·m
+ * and kN/m; springs in kN/mm and kN·m/rad; unit weight in kN/m³. Results are in kN, kN·m, mm
+ * and rad. Global axes X, Y, Z are right-handed with Z up; members use ETABS-style local axes
+ * 1, 2, 3 (see `memberAxes`). Member end forces and diagrams use right-hand rotations.
+ *
+ * Model and result are plain data (JSON-serialisable), so the analysis can run in a Web Worker.
+ *
+ * @param model - Nodes, members, restraints, load cases, loads tagged with their load case,
+ *   and combinations
+ * @param options - Diagram station count and detail
+ * @returns `cases` and `combinations` keyed by id, `envelope` over the combinations,
+ *   `calculationDetails`, `units` and `warnings`
+ * @throws {FrameAnalysisError} 301–314 for an invalid model (with `details` and every problem in
+ *   `issues`), 304 for an unstable structure
+ */
+export const analyzeFrame = (model, options = {}) => createFrameAnalysis(model, options).result;
+//# sourceMappingURL=analyze-frame.js.map

package/dist/esm/fea/frame/solver/analyze-frame.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"analyze-frame.js","sourceRoot":"","sources":["../../../../../src/fea/frame/solver/analyze-frame.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,EAAE,oBAAoB,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AAG3E,OAAO,EAAE,WAAW,EAAE,kBAAkB,EAAE,MAAM,8BAA8B,CAAC;AAC/E,OAAO,EAAE,mBAAmB,EAAE,MAAM,gBAAgB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,gBAAgB,CAAC;AACpD,OAAO,EAAE,aAAa,EAAE,aAAa,EAAiB,MAAM,oBAAoB,CAAC;AA8BjF,oFAAoF;AACpF,MAAM,WAAW,GAAG,CAAiC,QAAa,EAAE,EAAoB,EAAE;IACxF,MAAM,MAAM,GAAG,IAAI,GAAG,EAAe,CAAC;IACtC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxC,IAAI,KAAK;YAAE,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;;YACvB,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC,CAAC;AAEF;;;;;;;;;;;GAWG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAG,CACjC,KAAiB,EACjB,UAA+B,EAAE,EAClB,EAAE;IACjB,MAAM,MAAM,GAAG,kBAAkB,CAAC,KAAK,CAAC,CAAC;IACzC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC;QAAE,WAAW,CAAC,MAAM,CAAC,CAAC;IAE3C,MAAM,MAAM,GAAG,kBAAkB,CAAC,KAAK,CAAC,CAAC;IACzC,MAAM,YAAY,GAAG,KAAK,CAAC,YAAY,IAAI,EAAE,CAAC;IAC9C,MAAM,QAAQ,GAAG,mBAAmB,CAAC,YAAY,EAAE,IAAI,GAAG,CAAC,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC9F,MAAM,YAAY,GAAG,CAAC,OAAO,CAAC,QAAQ,IAAI,MAAM,CAAC,KAAK,MAAM,CAAC;IAE7D,MAAM,MAAM,GAAG,IAAI,GAAG,EAAoB,CAAC;IAC3C,MAAM,KAAK,GAAoC,EAAE,CAAC;IAClD,MAAM,WAAW,GAAG,WAAW,CAAY,KAAK,CAAC,UAAU,CAAC,CAAC;IAC7D,MAAM,YAAY,GAAG,WAAW,CAAa,KAAK,CAAC,WAAW,CAAC,CAAC;IAChE,KAAK,MAAM,QAAQ,IAAI,KAAK,CAAC,SAAS,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,aAAa,CAAC,MAAM,EAAE,QAAQ,EAAE;YAC5C,UAAU,EAAE,WAAW,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE;YAC9C,WAAW,EAAE,YAAY,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE;SACjD,CAAC,CAAC;QACH,MAAM,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC;QAC/B,KAAK,CAAC,QAAQ,CAAC,EAAE,CAAC,GAAG,YAAY,CAC/B,MAAM,EACN,KAAK,EACL,EAAE,EAAE,EAAE,QAAQ,CAAC,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,CAAC,IAAI,EAAE,IAAI,EAAE,QAAQ,CAAC,IAAI,EAAE,EAC3E,EAAE,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,YAAY,EAAE,CAC7C,CAAC;IACJ,CAAC;IAED,MAAM,iBAAiB,GAAG,IAAI,GAAG,EAAoB,CAAC;IACtD,MAAM,kBAAkB,GAAoC,EAAE,CAAC;IAC/D,KAAK,MAAM,WAAW,IAAI,YAAY,EAAE,CAAC;QACvC,2DAA2D;QAC3D,MAAM,KAAK,GAAG,CAAC,GAAI,QAAQ,CAAC,GAAG,CAAC,WAAW,CAAC,EAAE,CAAyB,CAAC,CAAC,GAAG,CAC1E,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,MAAM,CAAC,GAAG,CAAC,MAAM,CAAa,EAAE,MAAM,EAAE,CAAC,CAC1E,CAAC;QACF,MAAM,KAAK,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC;QACnC,iBAAiB,CAAC,GAAG,CAAC,WAAW,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC;QAC7C,kBAAkB,CAAC,WAAW,CAAC,EAAE,CAAC,GAAG,YAAY,CAC/C,MAAM,EACN,KAAK,EACL,EAAE,EAAE,EAAE,WAAW,CAAC,EAAE,EAAE,IAAI,EAAE,aAAa,EAAE,IAAI,EAAE,WAAW,CAAC,IAAI,EAAE,EACnE,EAAE,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,YAAY,EAAE,CAC7C,CAAC;IACJ,CAAC;IAED,MAAM,QAAQ,GAAa;QACzB,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC;QACvB,GAAG,MAAM,CAAC,MAAM,CAAC,kBAAkB,CAAC;KACrC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CACd,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QACrB,GAAG,CAAC;QACJ,OAAO,EAAE,GAAG,CAAC,CAAC,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,kBAAkB,KAAK,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,OAAO,EAAE;KAC3F,CAAC,CAAC,CACJ,CAAC;IAEF,MAAM,eAAe,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,kBAAkB,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1E,MAAM,MAAM,GAAgB;QAC1B,KAAK;QACL,YAAY,EAAE,kBAAkB;QAChC,QAAQ,EAAE,eAAe,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,SAAS;QACjF,kBAAkB,EAAE;YAClB,SAAS,EAAE,MAAM,CAAC,KAAK,CAAC,MAAM;YAC9B,WAAW,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM;YAClC,QAAQ,EAAE,MAAM,CAAC,QAAQ;YACzB,YAAY,EAAE,MAAM,CAAC,IAAI,CAAC,MAAM;YAChC,WAAW,EAAE,MAAM,CAAC,WAAW;YAC/B,aAAa,EAAE,KAAK,CAAC,SAAS,CAAC,MAAM;YACrC,gBAAgB,EAAE,YAAY,CAAC,MAAM;SACtC;QACD,KAAK,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE;QACrE,QAAQ;KACT,CAAC;IAEF,OAAO;QACL,MAAM;QACN,aAAa,EAAE,CAAC,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,EAAE;YAC9C,MAAM,KAAK,GAAG,iBAAiB,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,MAAM,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YACtE,IAAI,CAAC,KAAK;gBAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,QAAQ,wBAAwB,CAAC,CAAC;YAC3F,MAAM,CAAC,GAAG,MAAM,CAAC,WAAW,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YAC3C,IAAI,CAAC,KAAK,SAAS;gBAAE,MAAM,IAAI,KAAK,CAAC,UAAU,MAAM,CAAC,QAAQ,CAAC,uBAAuB,CAAC,CAAC;YACxF,OAAO,oBAAoB,CAAC,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,QAAQ,IAAI,OAAO,CAAC,QAAQ,CAAC,CAAC;QAC9E,CAAC;KACF,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,KAAiB,EAAE,UAA+B,EAAE,EAAe,EAAE,CAChG,mBAAmB,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC"}

package/dist/esm/fea/frame/solver/boundary-conditions.js

@@ -0,0 +1,33 @@
+const DOF_KEYS = ["ux", "uy", "uz", "rx", "ry", "rz"];
+/** Spring stiffness unit conversion: kN/mm → N/mm, kN·m/rad → N·mm/rad. */
+const SPRING_FACTOR = [1e3, 1e3, 1e3, 1e6, 1e6, 1e6];
+const parse = (value, factor) => {
+    if (value === true)
+        return { fixed: true, spring: 0 };
+    if (typeof value === "number" && value > 0)
+        return { fixed: false, spring: value * factor };
+    return { fixed: false, spring: 0 };
+};
+/**
+ * Collects restraints into per-DOF fixed flags and spring stiffnesses.
+ *
+ * @param restraints - Model restraints
+ * @param nodeIndex - Map from node id to node index (6 DOFs per node)
+ * @param dofCount - Total number of global DOFs
+ */
+export const buildBoundaryConditions = (restraints, nodeIndex, dofCount) => {
+    const fixed = new Array(dofCount).fill(false);
+    const spring = new Array(dofCount).fill(0);
+    for (const restraint of restraints) {
+        const base = 6 * nodeIndex.get(restraint.nodeId);
+        DOF_KEYS.forEach((key, k) => {
+            const parsed = parse(restraint[key], SPRING_FACTOR[k]);
+            if (parsed.fixed)
+                fixed[base + k] = true;
+            spring[base + k] = parsed.spring;
+        });
+    }
+    const restrained = fixed.map((f, i) => f || spring[i] > 0);
+    return { fixed, spring, restrained };
+};
+//# sourceMappingURL=boundary-conditions.js.map

package/dist/esm/fea/frame/solver/boundary-conditions.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"boundary-conditions.js","sourceRoot":"","sources":["../../../../../src/fea/frame/solver/boundary-conditions.ts"],"names":[],"mappings":"AAEA,MAAM,QAAQ,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AAE/D,2EAA2E;AAC3E,MAAM,aAAa,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;AAQrD,MAAM,KAAK,GAAG,CACZ,KAA+B,EAC/B,MAAc,EACsB,EAAE;IACtC,IAAI,KAAK,KAAK,IAAI;QAAE,OAAO,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;IACtD,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,GAAG,CAAC;QAAE,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,GAAG,MAAM,EAAE,CAAC;IAC5F,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;AACrC,CAAC,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAG,CACrC,UAAuB,EACvB,SAA2C,EAC3C,QAAgB,EACI,EAAE;IACtB,MAAM,KAAK,GAAG,IAAI,KAAK,CAAU,QAAQ,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACvD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACnD,KAAK,MAAM,SAAS,IAAI,UAAU,EAAE,CAAC;QACnC,MAAM,IAAI,GAAG,CAAC,GAAI,SAAS,CAAC,GAAG,CAAC,SAAS,CAAC,MAAM,CAAY,CAAC;QAC7D,QAAQ,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,MAAM,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;YACvD,IAAI,MAAM,CAAC,KAAK;gBAAE,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;YACzC,MAAM,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACnC,CAAC,CAAC,CAAC;IACL,CAAC;IACD,MAAM,UAAU,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3D,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC;AACvC,CAAC,CAAC"}

package/dist/esm/fea/frame/solver/combinations.js

@@ -0,0 +1,69 @@
+import { FrameAnalysisError, FrameErrors } from "../../../core/errors/frame-analysis.error.js";
+/**
+ * Expands every combination into net factors on load cases, following nested combinations.
+ *
+ * @param combinations - Model combinations
+ * @param caseIds - Ids of the model's load cases
+ * @returns For each combination id, the net factor on each load case it uses
+ * @throws {FrameAnalysisError} 309 for an unknown or missing reference, 310 for a loop,
+ *   311 for a non-finite factor or an empty combination
+ */
+export const resolveCombinations = (combinations, caseIds) => {
+    const byId = new Map(combinations.map((c) => [c.id, c]));
+    for (const combination of combinations) {
+        if (!Array.isArray(combination.terms) || combination.terms.length === 0) {
+            throw new FrameAnalysisError(FrameErrors.INVALID_LOAD_FACTOR, {
+                combinationId: combination.id,
+            });
+        }
+        for (const term of combination.terms) {
+            if (!Number.isFinite(term.factor))
+                throw new FrameAnalysisError(FrameErrors.INVALID_LOAD_FACTOR, {
+                    combinationId: combination.id,
+                });
+            const isCase = "case" in term && term.case !== undefined;
+            const isCombination = "combination" in term && term.combination !== undefined;
+            if (isCase === isCombination)
+                throw new FrameAnalysisError(FrameErrors.UNKNOWN_LOAD_CASE, {
+                    combinationId: combination.id,
+                });
+            if (isCase && !caseIds.has(term.case)) {
+                throw new FrameAnalysisError(FrameErrors.UNKNOWN_LOAD_CASE, {
+                    combinationId: combination.id,
+                });
+            }
+            if (isCombination && !byId.has(term.combination)) {
+                throw new FrameAnalysisError(FrameErrors.UNKNOWN_LOAD_CASE, {
+                    combinationId: combination.id,
+                });
+            }
+        }
+    }
+    const resolved = new Map();
+    const visiting = new Set();
+    const resolve = (id) => {
+        const done = resolved.get(id);
+        if (done)
+            return done;
+        if (visiting.has(id)) {
+            throw new FrameAnalysisError(FrameErrors.CIRCULAR_COMBINATION, { combinationId: id });
+        }
+        visiting.add(id);
+        const factors = new Map();
+        const addFactor = (caseId, factor) => factors.set(caseId, (factors.get(caseId) ?? 0) + factor);
+        for (const term of byId.get(id).terms) {
+            if ("case" in term) {
+                addFactor(term.case, term.factor);
+            }
+            else {
+                resolve(term.combination).forEach((f, caseId) => addFactor(caseId, term.factor * f));
+            }
+        }
+        visiting.delete(id);
+        resolved.set(id, factors);
+        return factors;
+    };
+    combinations.forEach((c) => resolve(c.id));
+    return resolved;
+};
+//# sourceMappingURL=combinations.js.map

package/dist/esm/fea/frame/solver/combinations.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"combinations.js","sourceRoot":"","sources":["../../../../../src/fea/frame/solver/combinations.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AAIxF;;;;;;;;GAQG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAG,CACjC,YAA+B,EAC/B,OAAoB,EACc,EAAE;IACpC,MAAM,IAAI,GAAG,IAAI,GAAG,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAEzD,KAAK,MAAM,WAAW,IAAI,YAAY,EAAE,CAAC;QACvC,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,WAAW,CAAC,KAAK,CAAC,IAAI,WAAW,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YACxE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,EAAE;gBAC5D,aAAa,EAAE,WAAW,CAAC,EAAE;aAC9B,CAAC,CAAC;QACL,CAAC;QACD,KAAK,MAAM,IAAI,IAAI,WAAW,CAAC,KAAK,EAAE,CAAC;YACrC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC;gBAC/B,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,EAAE;oBAC5D,aAAa,EAAE,WAAW,CAAC,EAAE;iBAC9B,CAAC,CAAC;YACL,MAAM,MAAM,GAAG,MAAM,IAAI,IAAI,IAAI,IAAI,CAAC,IAAI,KAAK,SAAS,CAAC;YACzD,MAAM,aAAa,GAAG,aAAa,IAAI,IAAI,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,CAAC;YAC9E,IAAI,MAAM,KAAK,aAAa;gBAC1B,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,iBAAiB,EAAE;oBAC1D,aAAa,EAAE,WAAW,CAAC,EAAE;iBAC9B,CAAC,CAAC;YACL,IAAI,MAAM,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;gBACtC,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,iBAAiB,EAAE;oBAC1D,aAAa,EAAE,WAAW,CAAC,EAAE;iBAC9B,CAAC,CAAC;YACL,CAAC;YACD,IAAI,aAAa,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,EAAE,CAAC;gBACjD,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,iBAAiB,EAAE;oBAC1D,aAAa,EAAE,WAAW,CAAC,EAAE;iBAC9B,CAAC,CAAC;YACL,CAAC;QACH,CAAC;IACH,CAAC;IAED,MAAM,QAAQ,GAAG,IAAI,GAAG,EAA+B,CAAC;IACxD,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAU,CAAC;IAEnC,MAAM,OAAO,GAAG,CAAC,EAAU,EAAuB,EAAE;QAClD,MAAM,IAAI,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAC9B,IAAI,IAAI;YAAE,OAAO,IAAI,CAAC;QACtB,IAAI,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC;YACrB,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,oBAAoB,EAAE,EAAE,aAAa,EAAE,EAAE,EAAE,CAAC,CAAC;QACxF,CAAC;QACD,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAEjB,MAAM,OAAO,GAAG,IAAI,GAAG,EAAkB,CAAC;QAC1C,MAAM,SAAS,GAAG,CAAC,MAAc,EAAE,MAAc,EAAE,EAAE,CACnD,OAAO,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;QAC3D,KAAK,MAAM,IAAI,IAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAqB,CAAC,KAAK,EAAE,CAAC;YAC3D,IAAI,MAAM,IAAI,IAAI,EAAE,CAAC;gBACnB,SAAS,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;YACpC,CAAC;iBAAM,CAAC;gBACN,OAAO,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,SAAS,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;YACvF,CAAC;QACH,CAAC;QAED,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC;QACpB,QAAQ,CAAC,GAAG,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;QAC1B,OAAO,OAAO,CAAC;IACjB,CAAC,CAAC;IAEF,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3C,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC"}

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

@@ -0,0 +1,122 @@
+import { SINGULARITY_TOLERANCE } from "../../../core/constants/frame.constant.js";
+import { FrameErrors } from "../../../core/errors/frame-analysis.error.js";
+import { createSkyline, reverseCuthillMcKee, skylineAdd, skylineFactor, skylineSize, } from "../../../utils/skyline.js";
+import { beam3dLocalStiffnessRightHanded } from "../elements/beam-3d-stiffness.js";
+import { condenseMember } from "../elements/static-condensation.js";
+import { memberAxes, transformationMatrix } from "../transformation/coordinate-transform.js";
+import { frameIssue, throwIssues } from "../validation/validate-model.js";
+import { buildBoundaryConditions } from "./boundary-conditions.js";
+import { memberDofs, memberGlobalStiffness } from "./global-assembly.js";
+const RELEASE_KEYS = ["u1", "u2", "u3", "r1", "r2", "r3"];
+const DOF_NAMES = ["ux", "uy", "uz", "rx", "ry", "rz"];
+const releasedIndices = (release, offset) => RELEASE_KEYS.flatMap((key, k) => (release?.[key] ? [offset + k] : []));
+/**
+ * Builds the member matrices of a model that passed `validateFrameModel`, assembles the free
+ * stiffness in skyline storage with nodes in reverse Cuthill–McKee order, and factorises it
+ * once for all load cases.
+ *
+ * @throws {FrameAnalysisError} 304 if the structure is a mechanism; `details` names a node and
+ *   a DOF involved in it
+ */
+export const prepareFrameSystem = (model) => {
+    const { nodes, members, restraints } = model;
+    const nodeIndex = new Map(nodes.map((n, i) => [n.id, i]));
+    const memberIndex = new Map(members.map((m, i) => [m.id, i]));
+    const dofCount = 6 * nodes.length;
+    const boundary = buildBoundaryConditions(restraints, nodeIndex, dofCount);
+    const systemMembers = members.map((member) => {
+        const start = nodes[nodeIndex.get(member.startNodeId)];
+        const end = nodes[nodeIndex.get(member.endNodeId)];
+        const { length, lambda } = memberAxes([start.x, start.y, start.z], [end.x, end.y, end.z], member.angle ?? 0);
+        const { E, G, A, I22, I33, J } = member.section;
+        const releasedDofs = [
+            ...releasedIndices(member.releases?.start, 0),
+            ...releasedIndices(member.releases?.end, 6),
+        ];
+        const condensed = condenseMember(beam3dLocalStiffnessRightHanded(E, G, A, I22, I33, J, length), releasedDofs);
+        return {
+            id: member.id,
+            section: member.section,
+            startIndex: nodeIndex.get(member.startNodeId),
+            endIndex: nodeIndex.get(member.endNodeId),
+            T: transformationMatrix(lambda),
+            k: condensed.k,
+            condenseFixedEndForces: condensed.condenseFixedEndForces,
+            released: Array.from({ length: 12 }, (_, i) => releasedDofs.includes(i)),
+            length,
+            lambda,
+        };
+    });
+    // ---- Equation numbering: nodes in RCM order, fixed DOFs skipped ------------------
+    const adjacency = nodes.map(() => []);
+    for (const m of systemMembers) {
+        adjacency[m.startIndex].push(m.endIndex);
+        adjacency[m.endIndex].push(m.startIndex);
+    }
+    const equation = new Int32Array(dofCount).fill(-1);
+    const free = [];
+    for (const node of reverseCuthillMcKee(adjacency)) {
+        for (let k = 0; k < 6; k++) {
+            const dof = 6 * node + k;
+            if (!boundary.fixed[dof]) {
+                equation[dof] = free.length;
+                free.push(dof);
+            }
+        }
+    }
+    // ---- Profile, assembly and factorisation -----------------------------------------
+    let factor = null;
+    if (free.length > 0) {
+        const first = free.map((_, e) => e);
+        const memberEquations = systemMembers.map((m) => memberDofs(m).map((d) => equation[d]));
+        for (const eqs of memberEquations) {
+            const active = eqs.filter((e) => e >= 0);
+            const low = Math.min(...active);
+            for (const e of active)
+                if (low < first[e])
+                    first[e] = low;
+        }
+        const K = createSkyline(first);
+        systemMembers.forEach((m, index) => {
+            const kGlobal = memberGlobalStiffness(m);
+            const eqs = memberEquations[index];
+            // a member's 12 DOFs are distinct, so each unordered pair and each diagonal is added once
+            for (let a = 0; a < 12; a++) {
+                if (eqs[a] < 0)
+                    continue;
+                for (let b = 0; b < 12; b++) {
+                    if (eqs[b] >= eqs[a])
+                        skylineAdd(K, eqs[a], eqs[b], kGlobal[a][b]);
+                }
+            }
+        });
+        free.forEach((dof, e) => {
+            if (boundary.spring[dof] > 0)
+                skylineAdd(K, e, e, boundary.spring[dof]);
+        });
+        const result = skylineFactor(K, SINGULARITY_TOLERANCE);
+        if (!result.ok) {
+            const dof = free[result.failedPivot];
+            throwIssues([
+                frameIssue(FrameErrors.STRUCTURE_UNSTABLE, {
+                    nodeId: nodes[Math.floor(dof / 6)].id,
+                    dof: DOF_NAMES[dof % 6],
+                }),
+            ]);
+        }
+        factor = K;
+    }
+    return {
+        nodes,
+        nodeIndex,
+        memberIndex,
+        members: systemMembers,
+        dofCount,
+        boundary,
+        equation,
+        free,
+        factor,
+        profileSize: factor ? skylineSize(factor) : 0,
+    };
+};
+//# sourceMappingURL=frame-system.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"frame-system.js","sourceRoot":"","sources":["../../../../../src/fea/frame/solver/frame-system.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,qBAAqB,EAAE,MAAM,oCAAoC,CAAC;AAC3E,OAAO,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AAGpE,OAAO,EACL,aAAa,EACb,mBAAmB,EACnB,UAAU,EACV,aAAa,EACb,WAAW,GAEZ,MAAM,oBAAoB,CAAC;AAE5B,OAAO,EAAE,+BAA+B,EAAE,MAAM,+BAA+B,CAAC;AAChF,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,UAAU,EAAE,oBAAoB,EAAE,MAAM,wCAAwC,CAAC;AAE1F,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,uBAAuB,EAA2B,MAAM,uBAAuB,CAAC;AACzF,OAAO,EAAE,UAAU,EAAE,qBAAqB,EAAuB,MAAM,mBAAmB,CAAC;AAE3F,MAAM,YAAY,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AACnE,MAAM,SAAS,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AAEhE,MAAM,eAAe,GAAG,CAAC,OAAqC,EAAE,MAAc,EAAY,EAAE,CAC1F,YAAY,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AA+BzE;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,KAAiB,EAAe,EAAE;IACnE,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,UAAU,EAAE,GAAG,KAAK,CAAC;IAC7C,MAAM,SAAS,GAAG,IAAI,GAAG,CAAkB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3E,MAAM,WAAW,GAAG,IAAI,GAAG,CAAkB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,QAAQ,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IAClC,MAAM,QAAQ,GAAG,uBAAuB,CAAC,UAAU,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;IAE1E,MAAM,aAAa,GAAmB,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QAC3D,MAAM,KAAK,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,WAAW,CAAW,CAAC,CAAC;QACjE,MAAM,GAAG,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAW,CAAC,CAAC;QAC7D,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,UAAU,CACnC,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAC3B,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EACrB,MAAM,CAAC,KAAK,IAAI,CAAC,CAClB,CAAC;QACF,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC;QAChD,MAAM,YAAY,GAAG;YACnB,GAAG,eAAe,CAAC,MAAM,CAAC,QAAQ,EAAE,KAAK,EAAE,CAAC,CAAC;YAC7C,GAAG,eAAe,CAAC,MAAM,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;SAC5C,CAAC;QACF,MAAM,SAAS,GAAG,cAAc,CAC9B,+BAA+B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,CAAC,EAC7D,YAAY,CACb,CAAC;QACF,OAAO;YACL,EAAE,EAAE,MAAM,CAAC,EAAE;YACb,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,UAAU,EAAE,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,WAAW,CAAW;YACvD,QAAQ,EAAE,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAW;YACnD,CAAC,EAAE,oBAAoB,CAAC,MAAM,CAAC;YAC/B,CAAC,EAAE,SAAS,CAAC,CAAC;YACd,sBAAsB,EAAE,SAAS,CAAC,sBAAsB;YACxD,QAAQ,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACxE,MAAM;YACN,MAAM;SACP,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,qFAAqF;IACrF,MAAM,SAAS,GAAe,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAClD,KAAK,MAAM,CAAC,IAAI,aAAa,EAAE,CAAC;QAC9B,SAAS,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;QACzC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC;IAC3C,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,UAAU,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,mBAAmB,CAAC,SAAS,CAAC,EAAE,CAAC;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YACzB,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;gBACzB,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;gBAC5B,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACjB,CAAC;QACH,CAAC;IACH,CAAC;IAED,qFAAqF;IACrF,IAAI,MAAM,GAAyB,IAAI,CAAC;IACxC,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACpC,MAAM,eAAe,GAAG,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxF,KAAK,MAAM,GAAG,IAAI,eAAe,EAAE,CAAC;YAClC,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;YAChC,KAAK,MAAM,CAAC,IAAI,MAAM;gBAAE,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBAAE,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAC7D,CAAC;QAED,MAAM,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC;QAC/B,aAAa,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE;YACjC,MAAM,OAAO,GAAG,qBAAqB,CAAC,CAAC,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,eAAe,CAAC,KAAK,CAAC,CAAC;YACnC,0FAA0F;YAC1F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC5B,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,SAAS;gBACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC5B,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;wBAAE,UAAU,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrE,CAAC;YACH,CAAC;QACH,CAAC,CAAC,CAAC;QACH,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE;YACtB,IAAI,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC;gBAAE,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1E,CAAC,CAAC,CAAC;QAEH,MAAM,MAAM,GAAG,aAAa,CAAC,CAAC,EAAE,qBAAqB,CAAC,CAAC;QACvD,IAAI,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC;YACf,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC;YACrC,WAAW,CAAC;gBACV,UAAU,CAAC,WAAW,CAAC,kBAAkB,EAAE;oBACzC,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE;oBACrC,GAAG,EAAE,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC;iBACxB,CAAC;aACH,CAAC,CAAC;QACL,CAAC;QACD,MAAM,GAAG,CAAC,CAAC;IACb,CAAC;IAED,OAAO;QACL,KAAK;QACL,SAAS;QACT,WAAW;QACX,OAAO,EAAE,aAAa;QACtB,QAAQ;QACR,QAAQ;QACR,QAAQ;QACR,IAAI;QACJ,MAAM;QACN,WAAW,EAAE,MAAM,CAAC,CAAC,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;KAC9C,CAAC;AACJ,CAAC,CAAC"}

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

@@ -0,0 +1,27 @@
+import { matMul, transpose } from "../../../utils/matrix.js";
+/** Global DOF numbers (6 per node) of a member's 12 end DOFs. */
+export const memberDofs = (member) => {
+    const dofs = [];
+    for (const node of [member.startIndex, member.endIndex]) {
+        for (let k = 0; k < 6; k++)
+            dofs.push(6 * node + k);
+    }
+    return dofs;
+};
+/** A member's stiffness in global axes, `Tᵀ·k·T` (N and mm). */
+export const memberGlobalStiffness = (member) => matMul(transpose(member.T), matMul(member.k, member.T));
+/**
+ * Adds `sign · Tᵀ·f` of one member's local end-force vector to a global vector (N and N·mm).
+ * With `sign = −1` and fixed-end forces this adds the equivalent nodal loads; with `sign = +1`
+ * and `k·u` it adds the member's elastic end forces.
+ */
+export const scatterMemberForces = (target, member, localForces, sign) => {
+    const dofs = memberDofs(member);
+    for (let a = 0; a < 12; a++) {
+        let value = 0;
+        for (let j = 0; j < 12; j++)
+            value += member.T[j][a] * localForces[j];
+        target[dofs[a]] += sign * value;
+    }
+};
+//# sourceMappingURL=global-assembly.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"global-assembly.js","sourceRoot":"","sources":["../../../../../src/fea/frame/solver/global-assembly.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,SAAS,EAAe,MAAM,mBAAmB,CAAC;AASnE,iEAAiE;AACjE,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,MAAuD,EAAY,EAAE;IAC9F,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,CAAC,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;IACtD,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC,CAAC;AAEF,gEAAgE;AAChE,MAAM,CAAC,MAAM,qBAAqB,GAAG,CAAC,MAAsB,EAAU,EAAE,CACtE,MAAM,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;AAE1D;;;;GAIG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAG,CACjC,MAAgB,EAChB,MAA6D,EAC7D,WAAqB,EACrB,IAAY,EACN,EAAE;IACR,MAAM,IAAI,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;QACtE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,KAAK,CAAC;IAClC,CAAC;AACH,CAAC,CAAC"}

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

@@ -0,0 +1,119 @@
+import { FrameAnalysisError, FrameErrors } from "../../../core/errors/frame-analysis.error.js";
+import { matVec } from "../../../utils/matrix.js";
+import { skylineSolve } from "../../../utils/skyline.js";
+import { fixedEndForces } from "../elements/fixed-end-forces.js";
+import { memberDofs, scatterMemberForces } from "./global-assembly.js";
+const GLOBAL_DIRECTIONS = {
+    X: [1, 0, 0],
+    Y: [0, 1, 0],
+    Z: [0, 0, 1],
+};
+/** Resolves a member load into local-axis loads in N and mm. */
+const toLocalLoads = (load, lambda, L) => {
+    const components = (direction) => {
+        if (direction === "1" || direction === "2" || direction === "3") {
+            return [{ axis: Number(direction), factor: 1 }];
+        }
+        const v = GLOBAL_DIRECTIONS[direction];
+        return [1, 2, 3]
+            .map((axis) => ({
+            axis,
+            factor: lambda[axis - 1][0] * v[0] + lambda[axis - 1][1] * v[1] + lambda[axis - 1][2] * v[2],
+        }))
+            .filter((c) => Math.abs(c.factor) > 1e-12);
+    };
+    return components(load.direction).map(({ axis, factor }) => load.type === "udl"
+        ? {
+            type: "udl",
+            axis,
+            w: load.w * factor, // kN/m = N/mm
+            start: load.start ?? 0,
+            end: load.end ?? L,
+        }
+        : { type: "point", axis, P: load.P * factor * 1e3, a: load.a });
+};
+/**
+ * Solves one load case against the prepared, factorised system.
+ *
+ * Self-weight is `selfWeight · unitWeight · A` per unit member length along global −Z.
+ *
+ * @param loadCase - The load case (id and self-weight multiplier)
+ * @param loads - The nodal and member loads tagged with this load case
+ *
+ * @throws {FrameAnalysisError} 308 for a load outside its member, 313 if self-weight is
+ *   requested for a member without `unitWeight`
+ */
+export const solveLoadCase = (system, loadCase, loads) => {
+    const { members, dofCount, boundary, equation, free, factor, nodeIndex, memberIndex } = system;
+    const F = new Array(dofCount).fill(0);
+    const loadsByMember = members.map(() => []);
+    for (const load of loads.memberLoads) {
+        loadsByMember[memberIndex.get(load.memberId)].push(load);
+    }
+    const selfWeight = loadCase.selfWeight ?? 0;
+    const memberLoads = members.map((member, m) => {
+        const memberLoadsHere = [...loadsByMember[m]];
+        if (selfWeight !== 0) {
+            const { unitWeight, A } = member.section;
+            if (unitWeight === undefined)
+                throw new FrameAnalysisError(FrameErrors.MISSING_UNIT_WEIGHT);
+            // kN/m³ · m² = kN/m
+            const w = -selfWeight * unitWeight * A * 1e-6;
+            if (w !== 0) {
+                memberLoadsHere.push({
+                    loadCase: loadCase.id,
+                    type: "udl",
+                    memberId: member.id,
+                    direction: "Z",
+                    w,
+                });
+            }
+        }
+        return memberLoadsHere.flatMap((load) => toLocalLoads(load, member.lambda, member.length));
+    });
+    const fixed = members.map((member, m) => {
+        const condensed = member.condenseFixedEndForces(fixedEndForces(memberLoads[m], member.length));
+        scatterMemberForces(F, member, condensed, -1);
+        return condensed;
+    });
+    for (const load of loads.nodalLoads) {
+        const base = 6 * nodeIndex.get(load.nodeId);
+        [load.fx, load.fy, load.fz].forEach((v, k) => (F[base + k] += (v ?? 0) * 1e3));
+        [load.mx, load.my, load.mz].forEach((v, k) => (F[base + 3 + k] += (v ?? 0) * 1e6));
+    }
+    const U = new Array(dofCount).fill(0);
+    if (factor) {
+        const solution = skylineSolve(factor, free.map((dof) => F[dof]));
+        free.forEach((dof) => (U[dof] = solution[equation[dof]]));
+    }
+    // Elastic end forces k·u per member; their global sum is K·U, from which R = K·U − F
+    const KU = new Array(dofCount).fill(0);
+    const memberForces = members.map((member, m) => {
+        const uLocal = matVec(member.T, memberDofs(member).map((d) => U[d]));
+        const elastic = matVec(member.k, uLocal);
+        scatterMemberForces(KU, member, elastic, 1);
+        return elastic.map((v, i) => v + fixed[m][i]);
+    });
+    const R = KU.map((v, i) => (boundary.restrained[i] ? v - F[i] : 0));
+    return { U, R, F, memberForces, memberLoads };
+};
+const scaleLoad = (load, factor) => load.type === "udl" ? { ...load, w: load.w * factor } : { ...load, P: load.P * factor };
+/** Superposes factored states: `Σ factor × state`. Member loads are concatenated, scaled. */
+export const combineStates = (terms) => {
+    const first = terms[0].state;
+    const add = (pick, length) => {
+        const out = new Array(length).fill(0);
+        for (const { state, factor } of terms) {
+            pick(state).forEach((v, i) => (out[i] += factor * v));
+        }
+        return out;
+    };
+    return {
+        U: add((s) => s.U, first.U.length),
+        R: add((s) => s.R, first.R.length),
+        F: add((s) => s.F, first.F.length),
+        memberForces: first.memberForces.map((_, m) => add((s) => s.memberForces[m], 12)),
+        memberLoads: first.memberLoads.map((_, m) => terms.flatMap(({ state, factor }) => factor === 0 ? [] : state.memberLoads[m].map((l) => scaleLoad(l, factor)))),
+    };
+};
+//# sourceMappingURL=load-case-solver.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"load-case-solver.js","sourceRoot":"","sources":["../../../../../src/fea/frame/solver/load-case-solver.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AACxF,OAAO,EAAE,MAAM,EAAe,MAAM,mBAAmB,CAAC;AACxD,OAAO,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AAElD,OAAO,EAAE,cAAc,EAAwB,MAAM,8BAA8B,CAAC;AAGpF,OAAO,EAAE,UAAU,EAAE,mBAAmB,EAAE,MAAM,mBAAmB,CAAC;AAcpE,MAAM,iBAAiB,GAAsD;IAC3E,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CACb,CAAC;AAEF,gEAAgE;AAChE,MAAM,YAAY,GAAG,CAAC,IAAgB,EAAE,MAAc,EAAE,CAAS,EAAqB,EAAE;IACtF,MAAM,UAAU,GAAG,CAAC,SAAwB,EAAyC,EAAE;QACrF,IAAI,SAAS,KAAK,GAAG,IAAI,SAAS,KAAK,GAAG,IAAI,SAAS,KAAK,GAAG,EAAE,CAAC;YAChE,OAAO,CAAC,EAAE,IAAI,EAAE,MAAM,CAAC,SAAS,CAAc,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC;QAC/D,CAAC;QACD,MAAM,CAAC,GAAG,iBAAiB,CAAC,SAAS,CAAC,CAAC;QACvC,OAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAW;aACxB,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,CAAC;YACd,IAAI;YACJ,MAAM,EACJ,MAAM,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;SACvF,CAAC,CAAC;aACF,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC;IAC/C,CAAC,CAAC;IAEF,OAAO,UAAU,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,EAAE,MAAM,EAAE,EAAE,EAAE,CACzD,IAAI,CAAC,IAAI,KAAK,KAAK;QACjB,CAAC,CAAC;YACE,IAAI,EAAE,KAAK;YACX,IAAI;YACJ,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,EAAE,cAAc;YAClC,KAAK,EAAE,IAAI,CAAC,KAAK,IAAI,CAAC;YACtB,GAAG,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC;SACnB;QACH,CAAC,CAAC,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,GAAG,GAAG,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CACjE,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAC3B,MAAmB,EACnB,QAAkB,EAClB,KAA6D,EACnD,EAAE;IACZ,MAAM,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,IAAI,EAAE,MAAM,EAAE,SAAS,EAAE,WAAW,EAAE,GAAG,MAAM,CAAC;IAC/F,MAAM,CAAC,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAE9C,MAAM,aAAa,GAAmB,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC5D,KAAK,MAAM,IAAI,IAAI,KAAK,CAAC,WAAW,EAAE,CAAC;QACrC,aAAa,CAAC,WAAW,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAW,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACrE,CAAC;IAED,MAAM,UAAU,GAAG,QAAQ,CAAC,UAAU,IAAI,CAAC,CAAC;IAC5C,MAAM,WAAW,GAAwB,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACjE,MAAM,eAAe,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,IAAI,UAAU,KAAK,CAAC,EAAE,CAAC;YACrB,MAAM,EAAE,UAAU,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC;YACzC,IAAI,UAAU,KAAK,SAAS;gBAAE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAC5F,oBAAoB;YACpB,MAAM,CAAC,GAAG,CAAC,UAAU,GAAG,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC;YAC9C,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;gBACZ,eAAe,CAAC,IAAI,CAAC;oBACnB,QAAQ,EAAE,QAAQ,CAAC,EAAE;oBACrB,IAAI,EAAE,KAAK;oBACX,QAAQ,EAAE,MAAM,CAAC,EAAE;oBACnB,SAAS,EAAE,GAAG;oBACd,CAAC;iBACF,CAAC,CAAC;YACL,CAAC;QACH,CAAC;QACD,OAAO,eAAe,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,YAAY,CAAC,IAAI,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC;IAC7F,CAAC,CAAC,CAAC;IAEH,MAAM,KAAK,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACtC,MAAM,SAAS,GAAG,MAAM,CAAC,sBAAsB,CAAC,cAAc,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC;QAC/F,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,SAAS,EAAE,CAAC,CAAC,CAAC,CAAC;QAC9C,OAAO,SAAS,CAAC;IACnB,CAAC,CAAC,CAAC;IAEH,KAAK,MAAM,IAAI,IAAI,KAAK,CAAC,UAAU,EAAE,CAAC;QACpC,MAAM,IAAI,GAAG,CAAC,GAAI,SAAS,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAY,CAAC;QACxD,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC/E,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACrF,CAAC;IAED,MAAM,CAAC,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9C,IAAI,MAAM,EAAE,CAAC;QACX,MAAM,QAAQ,GAAG,YAAY,CAC3B,MAAM,EACN,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAC1B,CAAC;QACF,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED,qFAAqF;IACrF,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC/C,MAAM,YAAY,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAC7C,MAAM,MAAM,GAAG,MAAM,CACnB,MAAM,CAAC,CAAC,EACR,UAAU,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CACpC,CAAC;QACF,MAAM,OAAO,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QACzC,mBAAmB,CAAC,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;QAC5C,OAAO,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,CAAC,CAAC,CAAC;IACH,MAAM,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAEpE,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,EAAE,WAAW,EAAE,CAAC;AAChD,CAAC,CAAC;AAEF,MAAM,SAAS,GAAG,CAAC,IAAqB,EAAE,MAAc,EAAmB,EAAE,CAC3E,IAAI,CAAC,IAAI,KAAK,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC;AAE1F,6FAA6F;AAC7F,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,KAA4C,EAAY,EAAE;IACtF,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC;IAC7B,MAAM,GAAG,GAAG,CAAC,IAA+B,EAAE,MAAc,EAAE,EAAE;QAC9D,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC9C,KAAK,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,KAAK,EAAE,CAAC;YACtC,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;QACxD,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC,CAAC;IACF,OAAO;QACL,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;QAClC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;QAClC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;QAClC,YAAY,EAAE,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QACjF,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAC1C,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,EAAE,EAAE,CAClC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,SAAS,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAC1E,CACF;KACF,CAAC;AACJ,CAAC,CAAC"}

package/dist/esm/fea/frame/transformation/coordinate-transform.js

@@ -0,0 +1,71 @@
+import { VERTICAL_MEMBER_TOLERANCE } from "../../../core/constants/frame.constant.js";
+import { FrameAnalysisError, FrameErrors } from "../../../core/errors/frame-analysis.error.js";
+import { zeros } from "../../../utils/matrix.js";
+const cross = (a, b) => [
+    a[1] * b[2] - a[2] * b[1],
+    a[2] * b[0] - a[0] * b[2],
+    a[0] * b[1] - a[1] * b[0],
+];
+const norm = (a) => Math.hypot(a[0], a[1], a[2]);
+const scale = (a, s) => [a[0] * s, a[1] * s, a[2] * s];
+/**
+ * ETABS-style local axes of a member as the rows of a 3×3 direction-cosine matrix.
+ *
+ * Local 1 runs from the start to the end node. For a non-vertical member local 3 is
+ * horizontal (`1 × Z`) and local 2 lies in the vertical plane pointing up; for a vertical
+ * member local 2 is global +X and local 3 is `1 × 2`. `angle` then rotates local 2 and 3
+ * about local 1 by the right-hand rule.
+ *
+ * @param start - Start node coordinates in mm
+ * @param end - End node coordinates in mm
+ * @param angle - Member angle in rad (default 0)
+ * @returns `{ length, lambda }` with `lambda` rows = local 1, 2, 3 in global components
+ * @throws {FrameAnalysisError} 301 if the nodes coincide
+ */
+export const memberAxes = (start, end, angle = 0) => {
+    const d = [end[0] - start[0], end[1] - start[1], end[2] - start[2]];
+    const length = norm(d);
+    if (!(length > 0))
+        throw new FrameAnalysisError(FrameErrors.ZERO_LENGTH_MEMBER);
+    const e1 = scale(d, 1 / length);
+    let e2;
+    let e3;
+    if (Math.hypot(e1[0], e1[1]) < VERTICAL_MEMBER_TOLERANCE) {
+        e2 = [1, 0, 0];
+        e3 = cross(e1, e2);
+    }
+    else {
+        const c = cross(e1, [0, 0, 1]);
+        e3 = scale(c, 1 / norm(c));
+        e2 = cross(e3, e1);
+    }
+    if (angle !== 0) {
+        const cos = Math.cos(angle);
+        const sin = Math.sin(angle);
+        const r2 = [
+            cos * e2[0] + sin * e3[0],
+            cos * e2[1] + sin * e3[1],
+            cos * e2[2] + sin * e3[2],
+        ];
+        const r3 = [
+            -sin * e2[0] + cos * e3[0],
+            -sin * e2[1] + cos * e3[1],
+            -sin * e2[2] + cos * e3[2],
+        ];
+        e2 = r2;
+        e3 = r3;
+    }
+    return { length, lambda: [e1, e2, e3] };
+};
+/** Builds the 12×12 transformation `T` (global → local) from the 3×3 direction cosines. */
+export const transformationMatrix = (lambda) => {
+    const T = zeros(12, 12);
+    for (let block = 0; block < 4; block++) {
+        for (let i = 0; i < 3; i++) {
+            for (let j = 0; j < 3; j++)
+                T[3 * block + i][3 * block + j] = lambda[i][j];
+        }
+    }
+    return T;
+};
+//# sourceMappingURL=coordinate-transform.js.map

package/dist/esm/fea/frame/transformation/coordinate-transform.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"coordinate-transform.js","sourceRoot":"","sources":["../../../../../src/fea/frame/transformation/coordinate-transform.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,yBAAyB,EAAE,MAAM,oCAAoC,CAAC;AAC/E,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AACxF,OAAO,EAAE,KAAK,EAAe,MAAM,mBAAmB,CAAC;AAIvD,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC;IACxC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;CAC1B,CAAC;AACF,MAAM,IAAI,GAAG,CAAC,CAAO,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvD,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAS,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AAE3E;;;;;;;;;;;;;GAaG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,CACxB,KAA+B,EAC/B,GAA6B,EAC7B,KAAK,GAAG,CAAC,EAC2B,EAAE;IACtC,MAAM,CAAC,GAAS,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1E,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QAAE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,kBAAkB,CAAC,CAAC;IAEhF,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;IAChC,IAAI,EAAQ,CAAC;IACb,IAAI,EAAQ,CAAC;IACb,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,yBAAyB,EAAE,CAAC;QACzD,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACf,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACrB,CAAC;SAAM,CAAC;QACN,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC/B,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3B,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACrB,CAAC;IAED,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC5B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAS;YACf,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;YACzB,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;YACzB,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;SAC1B,CAAC;QACF,MAAM,EAAE,GAAS;YACf,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;YAC1B,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;YAC1B,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;SAC3B,CAAC;QACF,EAAE,GAAG,EAAE,CAAC;QACR,EAAE,GAAG,EAAE,CAAC;IACV,CAAC;IAED,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC;AAC1C,CAAC,CAAC;AAEF,2FAA2F;AAC3F,MAAM,CAAC,MAAM,oBAAoB,GAAG,CAAC,MAAc,EAAU,EAAE;IAC7D,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;QACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7E,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC"}

package/dist/esm/fea/frame/types/model.type.js

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=model.type.js.map

package/dist/esm/fea/frame/types/model.type.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"model.type.js","sourceRoot":"","sources":["../../../../../src/fea/frame/types/model.type.ts"],"names":[],"mappings":""}

package/dist/esm/fea/frame/types/results.type.js

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=results.type.js.map