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

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

@@ -0,0 +1,136 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.analyzeFrame = exports.createFrameAnalysis = void 0;
+const distribute_area_loads_1 = require("../area/distribute-area-loads");
+const envelope_1 = require("../post/envelope");
+const present_state_1 = require("../post/present-state");
+const validate_model_1 = require("../validation/validate-model");
+const combinations_1 = require("./combinations");
+const frame_system_1 = require("./frame-system");
+const load_case_solver_1 = require("./load-case-solver");
+/** 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,
+ *   null areas, and combinations
+ * @param options - Diagram station count and detail
+ * @returns `{ result, memberDiagram }`
+ * @throws {FrameAnalysisError} 301–317 for an invalid model (with `details` and every problem in
+ *   `issues`), 304 for an unstable structure
+ */
+const createFrameAnalysis = (model, options = {}) => {
+    const issues = (0, validate_model_1.validateFrameModel)(model);
+    if (issues.length > 0)
+        (0, validate_model_1.throwIssues)(issues);
+    const system = (0, frame_system_1.prepareFrameSystem)(model);
+    const combinations = model.combinations ?? [];
+    const resolved = (0, combinations_1.resolveCombinations)(combinations, new Set(model.loadCases.map((c) => c.id)));
+    const withStations = (options.diagrams ?? "full") === "full";
+    const fromAreas = (0, distribute_area_loads_1.areaMemberLoads)(model);
+    const states = new Map();
+    const cases = {};
+    const nodalByCase = groupByCase(model.nodalLoads);
+    const memberByCase = groupByCase([
+        ...(model.memberLoads ?? []),
+        ...fromAreas.memberLoads,
+    ]);
+    for (const loadCase of model.loadCases) {
+        const state = (0, load_case_solver_1.solveLoadCase)(system, loadCase, {
+            nodalLoads: nodalByCase.get(loadCase.id) ?? [],
+            memberLoads: memberByCase.get(loadCase.id) ?? [],
+        });
+        states.set(loadCase.id, state);
+        cases[loadCase.id] = (0, present_state_1.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 = (0, load_case_solver_1.combineStates)(terms);
+        combinationStates.set(combination.id, state);
+        combinationResults[combination.id] = (0, present_state_1.presentState)(system, state, { id: combination.id, kind: "combination", name: combination.name }, { stations: options.stations, withStations });
+    }
+    const warnings = [
+        ...fromAreas.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 ? (0, envelope_1.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 (0, present_state_1.presentMemberDiagram)(system, state, m, stations ?? options.stations);
+        },
+    };
+};
+exports.createFrameAnalysis = createFrameAnalysis;
+/**
+ * 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.
+ *
+ * Area loads on null areas (`areas`, `areaLoads`) are first turned into member loads on the
+ * members along each area's edges (see `distributeAreaLoads`); the areas add no stiffness.
+ *
+ * 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,
+ *   null areas, 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–317 for an invalid model (with `details` and every problem in
+ *   `issues`), 304 for an unstable structure
+ */
+const analyzeFrame = (model, options = {}) => (0, exports.createFrameAnalysis)(model, options).result;
+exports.analyzeFrame = analyzeFrame;
+//# sourceMappingURL=analyze-frame.js.map

package/dist/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,yEAAgE;AAChE,+CAAiD;AACjD,yDAA2E;AAG3E,iEAA+E;AAC/E,iDAAqD;AACrD,iDAAoD;AACpD,yDAAiF;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;AACI,MAAM,mBAAmB,GAAG,CACjC,KAAiB,EACjB,UAA+B,EAAE,EAClB,EAAE;IACjB,MAAM,MAAM,GAAG,IAAA,mCAAkB,EAAC,KAAK,CAAC,CAAC;IACzC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC;QAAE,IAAA,4BAAW,EAAC,MAAM,CAAC,CAAC;IAE3C,MAAM,MAAM,GAAG,IAAA,iCAAkB,EAAC,KAAK,CAAC,CAAC;IACzC,MAAM,YAAY,GAAG,KAAK,CAAC,YAAY,IAAI,EAAE,CAAC;IAC9C,MAAM,QAAQ,GAAG,IAAA,kCAAmB,EAAC,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;IAC7D,MAAM,SAAS,GAAG,IAAA,uCAAe,EAAC,KAAK,CAAC,CAAC;IAEzC,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;QAC3C,GAAG,CAAC,KAAK,CAAC,WAAW,IAAI,EAAE,CAAC;QAC5B,GAAG,SAAS,CAAC,WAAW;KACzB,CAAC,CAAC;IACH,KAAK,MAAM,QAAQ,IAAI,KAAK,CAAC,SAAS,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,IAAA,gCAAa,EAAC,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,IAAA,4BAAY,EAC/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,IAAA,gCAAa,EAAC,KAAK,CAAC,CAAC;QACnC,iBAAiB,CAAC,GAAG,CAAC,WAAW,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC;QAC7C,kBAAkB,CAAC,WAAW,CAAC,EAAE,CAAC,GAAG,IAAA,4BAAY,EAC/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,SAAS,CAAC,QAAQ;QACrB,GAAG,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,kBAAkB,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAC/E,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;YACrB,GAAG,CAAC;YACJ,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;SAC3F,CAAC,CAAC,CACJ;KACF,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,IAAA,wBAAa,EAAC,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,IAAA,oCAAoB,EAAC,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,QAAQ,IAAI,OAAO,CAAC,QAAQ,CAAC,CAAC;QAC9E,CAAC;KACF,CAAC;AACJ,CAAC,CAAC;AAzFW,QAAA,mBAAmB,uBAyF9B;AAEF;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2BG;AACI,MAAM,YAAY,GAAG,CAAC,KAAiB,EAAE,UAA+B,EAAE,EAAe,EAAE,CAChG,IAAA,2BAAmB,EAAC,KAAK,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC;AADhC,QAAA,YAAY,gBACoB"}

package/dist/fea/frame/solver/boundary-conditions.d.ts

@@ -0,0 +1,14 @@
+import type { Restraint } from "../types/model.type";
+export interface BoundaryConditions {
+    fixed: boolean[];
+    spring: number[];
+    restrained: boolean[];
+}
+/**
+ * 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 declare const buildBoundaryConditions: (restraints: Restraint[], nodeIndex: Map<Restraint["nodeId"], number>, dofCount: number) => BoundaryConditions;

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

@@ -0,0 +1,37 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.buildBoundaryConditions = void 0;
+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
+ */
+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 };
+};
+exports.buildBoundaryConditions = buildBoundaryConditions;
+//# sourceMappingURL=boundary-conditions.js.map

package/dist/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;AACI,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;AAjBW,QAAA,uBAAuB,2BAiBlC"}

package/dist/fea/frame/solver/combinations.d.ts

@@ -0,0 +1,11 @@
+import type { LoadCombination } from "../types/model.type";
+/**
+ * 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 declare const resolveCombinations: (combinations: LoadCombination[], caseIds: Set<string>) => Map<string, Map<string, number>>;

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

@@ -0,0 +1,73 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.resolveCombinations = void 0;
+const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.error");
+/**
+ * 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
+ */
+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 frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.INVALID_LOAD_FACTOR, {
+                combinationId: combination.id,
+            });
+        }
+        for (const term of combination.terms) {
+            if (!Number.isFinite(term.factor))
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.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 frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.UNKNOWN_LOAD_CASE, {
+                    combinationId: combination.id,
+                });
+            if (isCase && !caseIds.has(term.case)) {
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.UNKNOWN_LOAD_CASE, {
+                    combinationId: combination.id,
+                });
+            }
+            if (isCombination && !byId.has(term.combination)) {
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.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 frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.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;
+};
+exports.resolveCombinations = resolveCombinations;
+//# sourceMappingURL=combinations.js.map

package/dist/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,gFAAwF;AAIxF;;;;;;;;GAQG;AACI,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,yCAAkB,CAAC,kCAAW,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,yCAAkB,CAAC,kCAAW,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,yCAAkB,CAAC,kCAAW,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,yCAAkB,CAAC,kCAAW,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,yCAAkB,CAAC,kCAAW,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,yCAAkB,CAAC,kCAAW,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;AAjEW,QAAA,mBAAmB,uBAiE9B"}

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

@@ -0,0 +1,42 @@
+import type { SectionProperties } from "../../../core/types/section-properties.type";
+import type { Matrix } from "../../../utils/matrix";
+import { type SkylineMatrix } from "../../../utils/skyline";
+import type { FrameId, FrameModel, FrameNode } from "../types/model.type";
+import { type BoundaryConditions } from "./boundary-conditions";
+import { type AssemblyMember } from "./global-assembly";
+/** A member with everything that does not depend on the loads. */
+export interface SystemMember extends AssemblyMember {
+    id: FrameId;
+    section: SectionProperties;
+    length: number;
+    lambda: Matrix;
+    released: boolean[];
+    /** Condenses a member's fixed-end forces for its end releases. */
+    condenseFixedEndForces: (fixed: number[]) => number[];
+}
+/** The load-independent part of a frame analysis: geometry, stiffness and its factorisation. */
+export interface FrameSystem {
+    nodes: FrameNode[];
+    nodeIndex: Map<FrameId, number>;
+    memberIndex: Map<FrameId, number>;
+    members: SystemMember[];
+    dofCount: number;
+    boundary: BoundaryConditions;
+    /** Equation number of each global DOF in the factorised system, −1 for a fixed DOF. */
+    equation: Int32Array;
+    /** Global DOF of each equation. */
+    free: number[];
+    /** Skyline Cholesky factor of the free stiffness including springs; null when nothing is free. */
+    factor: SkylineMatrix | null;
+    /** Stored entries of the factor (a measure of bandwidth, memory and work). */
+    profileSize: number;
+}
+/**
+ * Builds the member matrices of a model that passed `validateFrameModel`, assembles the free
+ * stiffness in skyline storage with nodes in reverse Cuthill–McKee order, and factorises it
+ * once for all load cases.
+ *
+ * @throws {FrameAnalysisError} 304 if the structure is a mechanism; `details` names a node and
+ *   a DOF involved in it
+ */
+export declare const prepareFrameSystem: (model: FrameModel) => FrameSystem;

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

@@ -0,0 +1,126 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.prepareFrameSystem = void 0;
+const frame_constant_1 = require("../../../core/constants/frame.constant");
+const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.error");
+const skyline_1 = require("../../../utils/skyline");
+const beam_3d_stiffness_1 = require("../elements/beam-3d-stiffness");
+const static_condensation_1 = require("../elements/static-condensation");
+const coordinate_transform_1 = require("../transformation/coordinate-transform");
+const validate_model_1 = require("../validation/validate-model");
+const boundary_conditions_1 = require("./boundary-conditions");
+const global_assembly_1 = require("./global-assembly");
+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
+ */
+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 = (0, boundary_conditions_1.buildBoundaryConditions)(restraints, nodeIndex, dofCount);
+    const systemMembers = members.map((member) => {
+        const start = nodes[nodeIndex.get(member.startNodeId)];
+        const end = nodes[nodeIndex.get(member.endNodeId)];
+        const { length, lambda } = (0, coordinate_transform_1.memberAxes)([start.x, start.y, start.z], [end.x, end.y, end.z], member.angle ?? 0);
+        const { E, G, A, I22, I33, J } = member.section;
+        const releasedDofs = [
+            ...releasedIndices(member.releases?.start, 0),
+            ...releasedIndices(member.releases?.end, 6),
+        ];
+        const condensed = (0, static_condensation_1.condenseMember)((0, beam_3d_stiffness_1.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: (0, coordinate_transform_1.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 (0, skyline_1.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) => (0, global_assembly_1.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 = (0, skyline_1.createSkyline)(first);
+        systemMembers.forEach((m, index) => {
+            const kGlobal = (0, global_assembly_1.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])
+                        (0, skyline_1.skylineAdd)(K, eqs[a], eqs[b], kGlobal[a][b]);
+                }
+            }
+        });
+        free.forEach((dof, e) => {
+            if (boundary.spring[dof] > 0)
+                (0, skyline_1.skylineAdd)(K, e, e, boundary.spring[dof]);
+        });
+        const result = (0, skyline_1.skylineFactor)(K, frame_constant_1.SINGULARITY_TOLERANCE);
+        if (!result.ok) {
+            const dof = free[result.failedPivot];
+            (0, validate_model_1.throwIssues)([
+                (0, validate_model_1.frameIssue)(frame_analysis_error_1.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 ? (0, skyline_1.skylineSize)(factor) : 0,
+    };
+};
+exports.prepareFrameSystem = prepareFrameSystem;
+//# sourceMappingURL=frame-system.js.map

package/dist/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,uEAA2E;AAC3E,gFAAoE;AAGpE,gDAO4B;AAE5B,qEAAgF;AAChF,yEAAiE;AACjE,iFAA0F;AAE1F,iEAAuE;AACvE,+DAAyF;AACzF,uDAA2F;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;AACI,MAAM,kBAAkB,GAAG,CAAC,KAAiB,EAAe,EAAE;IACnE,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,UAAU,EAAE,GAAG,KAAK,CAAC;IAC7C,MAAM,SAAS,GAAG,IAAI,GAAG,CAAkB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3E,MAAM,WAAW,GAAG,IAAI,GAAG,CAAkB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,QAAQ,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IAClC,MAAM,QAAQ,GAAG,IAAA,6CAAuB,EAAC,UAAU,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;IAE1E,MAAM,aAAa,GAAmB,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QAC3D,MAAM,KAAK,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,WAAW,CAAW,CAAC,CAAC;QACjE,MAAM,GAAG,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAW,CAAC,CAAC;QAC7D,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,IAAA,iCAAU,EACnC,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAC3B,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EACrB,MAAM,CAAC,KAAK,IAAI,CAAC,CAClB,CAAC;QACF,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC;QAChD,MAAM,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,IAAA,oCAAc,EAC9B,IAAA,mDAA+B,EAAC,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,IAAA,2CAAoB,EAAC,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,IAAA,6BAAmB,EAAC,SAAS,CAAC,EAAE,CAAC;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YACzB,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;gBACzB,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;gBAC5B,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACjB,CAAC;QACH,CAAC;IACH,CAAC;IAED,qFAAqF;IACrF,IAAI,MAAM,GAAyB,IAAI,CAAC;IACxC,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACpC,MAAM,eAAe,GAAG,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAA,4BAAU,EAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxF,KAAK,MAAM,GAAG,IAAI,eAAe,EAAE,CAAC;YAClC,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;YAChC,KAAK,MAAM,CAAC,IAAI,MAAM;gBAAE,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBAAE,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAC7D,CAAC;QAED,MAAM,CAAC,GAAG,IAAA,uBAAa,EAAC,KAAK,CAAC,CAAC;QAC/B,aAAa,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE;YACjC,MAAM,OAAO,GAAG,IAAA,uCAAqB,EAAC,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,IAAA,oBAAU,EAAC,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,IAAA,oBAAU,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1E,CAAC,CAAC,CAAC;QAEH,MAAM,MAAM,GAAG,IAAA,uBAAa,EAAC,CAAC,EAAE,sCAAqB,CAAC,CAAC;QACvD,IAAI,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC;YACf,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC;YACrC,IAAA,4BAAW,EAAC;gBACV,IAAA,2BAAU,EAAC,kCAAW,CAAC,kBAAkB,EAAE;oBACzC,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE;oBACrC,GAAG,EAAE,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC;iBACxB,CAAC;aACH,CAAC,CAAC;QACL,CAAC;QACD,MAAM,GAAG,CAAC,CAAC;IACb,CAAC;IAED,OAAO;QACL,KAAK;QACL,SAAS;QACT,WAAW;QACX,OAAO,EAAE,aAAa;QACtB,QAAQ;QACR,QAAQ;QACR,QAAQ;QACR,IAAI;QACJ,MAAM;QACN,WAAW,EAAE,MAAM,CAAC,CAAC,CAAC,IAAA,qBAAW,EAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;KAC9C,CAAC;AACJ,CAAC,CAAC;AA5GW,QAAA,kBAAkB,sBA4G7B"}

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

@@ -0,0 +1,17 @@
+import { type Matrix } from "../../../utils/matrix";
+export interface AssemblyMember {
+    startIndex: number;
+    endIndex: number;
+    T: Matrix;
+    k: Matrix;
+}
+/** Global DOF numbers (6 per node) of a member's 12 end DOFs. */
+export declare const memberDofs: (member: Pick<AssemblyMember, "startIndex" | "endIndex">) => number[];
+/** A member's stiffness in global axes, `Tᵀ·k·T` (N and mm). */
+export declare const memberGlobalStiffness: (member: AssemblyMember) => Matrix;
+/**
+ * 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 declare const scatterMemberForces: (target: number[], member: Pick<AssemblyMember, "startIndex" | "endIndex" | "T">, localForces: number[], sign: 1 | -1) => void;

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

@@ -0,0 +1,33 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.scatterMemberForces = exports.memberGlobalStiffness = exports.memberDofs = void 0;
+const matrix_1 = require("../../../utils/matrix");
+/** Global DOF numbers (6 per node) of a member's 12 end DOFs. */
+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;
+};
+exports.memberDofs = memberDofs;
+/** A member's stiffness in global axes, `Tᵀ·k·T` (N and mm). */
+const memberGlobalStiffness = (member) => (0, matrix_1.matMul)((0, matrix_1.transpose)(member.T), (0, matrix_1.matMul)(member.k, member.T));
+exports.memberGlobalStiffness = memberGlobalStiffness;
+/**
+ * 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.
+ */
+const scatterMemberForces = (target, member, localForces, sign) => {
+    const dofs = (0, exports.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;
+    }
+};
+exports.scatterMemberForces = scatterMemberForces;
+//# sourceMappingURL=global-assembly.js.map

package/dist/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,8CAAmE;AASnE,iEAAiE;AAC1D,MAAM,UAAU,GAAG,CAAC,MAAuD,EAAY,EAAE;IAC9F,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,CAAC,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;IACtD,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC,CAAC;AANW,QAAA,UAAU,cAMrB;AAEF,gEAAgE;AACzD,MAAM,qBAAqB,GAAG,CAAC,MAAsB,EAAU,EAAE,CACtE,IAAA,eAAM,EAAC,IAAA,kBAAS,EAAC,MAAM,CAAC,CAAC,CAAC,EAAE,IAAA,eAAM,EAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;AAD7C,QAAA,qBAAqB,yBACwB;AAE1D;;;;GAIG;AACI,MAAM,mBAAmB,GAAG,CACjC,MAAgB,EAChB,MAA6D,EAC7D,WAAqB,EACrB,IAAY,EACN,EAAE;IACR,MAAM,IAAI,GAAG,IAAA,kBAAU,EAAC,MAAM,CAAC,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;QACtE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,KAAK,CAAC;IAClC,CAAC;AACH,CAAC,CAAC;AAZW,QAAA,mBAAmB,uBAY9B"}

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

@@ -0,0 +1,34 @@
+import { type LocalMemberLoad } from "../elements/fixed-end-forces";
+import type { LoadCase, MemberLoad, NodalLoad } from "../types/model.type";
+import type { FrameSystem } from "./frame-system";
+/**
+ * Unrounded result of one load case or combination, in N and mm. Every field is linear in the
+ * loads, so combinations are formed by factoring and adding states.
+ */
+export interface RawState {
+    U: number[];
+    R: number[];
+    F: number[];
+    memberForces: number[][];
+    memberLoads: LocalMemberLoad[][];
+}
+/**
+ * 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 declare const solveLoadCase: (system: FrameSystem, loadCase: LoadCase, loads: {
+    nodalLoads: NodalLoad[];
+    memberLoads: MemberLoad[];
+}) => RawState;
+/** Superposes factored states: `Σ factor × state`. Member loads are concatenated, scaled. */
+export declare const combineStates: (terms: {
+    state: RawState;
+    factor: number;
+}[]) => RawState;