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

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

@@ -0,0 +1,138 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.combineStates = exports.solveLoadCase = void 0;
+const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.error");
+const matrix_1 = require("../../../utils/matrix");
+const skyline_1 = require("../../../utils/skyline");
+const fixed_end_forces_1 = require("../elements/fixed-end-forces");
+const global_assembly_1 = require("./global-assembly");
+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);
+    };
+    if (load.type === "torque") {
+        // kN·m/m = 10³ N·mm/mm
+        return [{ type: "torque", t: load.t * 1e3, start: load.start ?? 0, end: load.end ?? L }];
+    }
+    return components(load.direction).map(({ axis, factor }) => {
+        if (load.type === "point")
+            return { type: "point", axis, P: load.P * factor * 1e3, a: load.a };
+        const [w1, w2] = load.type === "udl" ? [load.w, load.w] : [load.w1, load.w2];
+        return {
+            type: "distributed",
+            axis,
+            w1: w1 * factor, // kN/m = N/mm
+            w2: w2 * factor,
+            start: load.start ?? 0,
+            end: load.end ?? L,
+        };
+    });
+};
+/**
+ * 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`
+ */
+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 frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.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((0, fixed_end_forces_1.fixedEndForces)(memberLoads[m], member.length));
+        (0, global_assembly_1.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 = (0, skyline_1.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 = (0, matrix_1.matVec)(member.T, (0, global_assembly_1.memberDofs)(member).map((d) => U[d]));
+        const elastic = (0, matrix_1.matVec)(member.k, uLocal);
+        (0, global_assembly_1.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 };
+};
+exports.solveLoadCase = solveLoadCase;
+const scaleLoad = (load, factor) => {
+    if (load.type === "distributed")
+        return { ...load, w1: load.w1 * factor, w2: load.w2 * factor };
+    if (load.type === "torque")
+        return { ...load, t: load.t * factor };
+    return { ...load, P: load.P * factor };
+};
+/** Superposes factored states: `Σ factor × state`. Member loads are concatenated, scaled. */
+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)))),
+    };
+};
+exports.combineStates = combineStates;
+//# sourceMappingURL=load-case-solver.js.map

package/dist/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,gFAAwF;AACxF,8CAAwD;AACxD,gDAAkD;AAElD,mEAAoF;AAGpF,uDAAoE;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,IAAI,IAAI,CAAC,IAAI,KAAK,QAAQ,EAAE,CAAC;QAC3B,uBAAuB;QACvB,OAAO,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,GAAG,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,IAAI,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IAC3F,CAAC;IACD,OAAO,UAAU,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,EAAE,MAAM,EAAE,EAAmB,EAAE;QAC1E,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO;YAAE,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,GAAG,GAAG,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;QAC/F,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAC7E,OAAO;YACL,IAAI,EAAE,aAAa;YACnB,IAAI;YACJ,EAAE,EAAE,EAAE,GAAG,MAAM,EAAE,cAAc;YAC/B,EAAE,EAAE,EAAE,GAAG,MAAM;YACf,KAAK,EAAE,IAAI,CAAC,KAAK,IAAI,CAAC;YACtB,GAAG,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC;SACnB,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC;AAEF;;;;;;;;;;GAUG;AACI,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,yCAAkB,CAAC,kCAAW,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,IAAA,iCAAc,EAAC,WAAW,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC;QAC/F,IAAA,qCAAmB,EAAC,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,IAAA,sBAAY,EAC3B,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,IAAA,eAAM,EACnB,MAAM,CAAC,CAAC,EACR,IAAA,4BAAU,EAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CACpC,CAAC;QACF,MAAM,OAAO,GAAG,IAAA,eAAM,EAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QACzC,IAAA,qCAAmB,EAAC,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;AArEW,QAAA,aAAa,iBAqExB;AAEF,MAAM,SAAS,GAAG,CAAC,IAAqB,EAAE,MAAc,EAAmB,EAAE;IAC3E,IAAI,IAAI,CAAC,IAAI,KAAK,aAAa;QAAE,OAAO,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,MAAM,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,MAAM,EAAE,CAAC;IAChG,IAAI,IAAI,CAAC,IAAI,KAAK,QAAQ;QAAE,OAAO,EAAE,GAAG,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC;IACnE,OAAO,EAAE,GAAG,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC;AACzC,CAAC,CAAC;AAEF,6FAA6F;AACtF,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;AApBW,QAAA,aAAa,iBAoBxB"}

package/dist/fea/frame/transformation/coordinate-transform.d.ts

@@ -0,0 +1,21 @@
+import { type Matrix } from "../../../utils/matrix";
+/**
+ * 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 declare const memberAxes: (start: [number, number, number], end: [number, number, number], angle?: number) => {
+    length: number;
+    lambda: Matrix;
+};
+/** Builds the 12×12 transformation `T` (global → local) from the 3×3 direction cosines. */
+export declare const transformationMatrix: (lambda: Matrix) => Matrix;

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

@@ -0,0 +1,76 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.transformationMatrix = exports.memberAxes = void 0;
+const frame_constant_1 = require("../../../core/constants/frame.constant");
+const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.error");
+const matrix_1 = require("../../../utils/matrix");
+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
+ */
+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 frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.ZERO_LENGTH_MEMBER);
+    const e1 = scale(d, 1 / length);
+    let e2;
+    let e3;
+    if (Math.hypot(e1[0], e1[1]) < frame_constant_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] };
+};
+exports.memberAxes = memberAxes;
+/** Builds the 12×12 transformation `T` (global → local) from the 3×3 direction cosines. */
+const transformationMatrix = (lambda) => {
+    const T = (0, matrix_1.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;
+};
+exports.transformationMatrix = transformationMatrix;
+//# sourceMappingURL=coordinate-transform.js.map

package/dist/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,uEAA+E;AAC/E,gFAAwF;AACxF,8CAAuD;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;AACI,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,yCAAkB,CAAC,kCAAW,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,0CAAyB,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;AAvCW,QAAA,UAAU,cAuCrB;AAEF,2FAA2F;AACpF,MAAM,oBAAoB,GAAG,CAAC,MAAc,EAAU,EAAE;IAC7D,MAAM,CAAC,GAAG,IAAA,cAAK,EAAC,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;AARW,QAAA,oBAAoB,wBAQ/B"}

package/dist/fea/frame/types/model.type.d.ts

@@ -0,0 +1,170 @@
+import type { SectionProperties } from "../../../core/types/section-properties.type";
+export type FrameId = string | number;
+export interface FrameNode {
+    id: FrameId;
+    x: number;
+    y: number;
+    z: number;
+}
+/** A restrained DOF: `true` is fixed, a number is a spring stiffness, omitted/`false` is free. */
+export type DofRestraint = boolean | number;
+export interface Restraint {
+    nodeId: FrameId;
+    ux?: DofRestraint;
+    uy?: DofRestraint;
+    uz?: DofRestraint;
+    rx?: DofRestraint;
+    ry?: DofRestraint;
+    rz?: DofRestraint;
+}
+/** Local end DOFs of a member that transmit no force (hinges, rollers…). */
+export interface MemberEndRelease {
+    u1?: boolean;
+    u2?: boolean;
+    u3?: boolean;
+    r1?: boolean;
+    r2?: boolean;
+    r3?: boolean;
+}
+export interface Member {
+    id: FrameId;
+    startNodeId: FrameId;
+    endNodeId: FrameId;
+    section: SectionProperties;
+    angle?: number;
+    releases?: {
+        start?: MemberEndRelease;
+        end?: MemberEndRelease;
+    };
+}
+export interface NodalLoad {
+    loadCase: string;
+    nodeId: FrameId;
+    fx?: number;
+    fy?: number;
+    fz?: number;
+    mx?: number;
+    my?: number;
+    mz?: number;
+}
+/** Load direction: local member axes "1", "2", "3" or global axes "X", "Y", "Z". */
+export type LoadDirection = "1" | "2" | "3" | "X" | "Y" | "Z";
+/** Uniform load per unit length of member, along the whole member or from `start` to `end`. */
+export interface UniformMemberLoad {
+    loadCase: string;
+    type: "udl";
+    memberId: FrameId;
+    direction: LoadDirection;
+    w: number;
+    start?: number;
+    end?: number;
+}
+export interface PointMemberLoad {
+    loadCase: string;
+    type: "point";
+    memberId: FrameId;
+    direction: LoadDirection;
+    P: number;
+    a: number;
+}
+/**
+ * Linearly varying load per unit length of member, from `w1` at `start` to `w2` at `end`
+ * (triangular when one of them is 0).
+ */
+export interface TrapezoidalMemberLoad {
+    loadCase: string;
+    type: "trapezoidal";
+    memberId: FrameId;
+    direction: LoadDirection;
+    w1: number;
+    w2: number;
+    start?: number;
+    end?: number;
+}
+/**
+ * Uniform torque per unit length of member, about local axis 1 (right-hand rule), along the
+ * whole member or from `start` to `end`.
+ */
+export interface TorqueMemberLoad {
+    loadCase: string;
+    type: "torque";
+    memberId: FrameId;
+    t: number;
+    start?: number;
+    end?: number;
+}
+export type MemberLoad = UniformMemberLoad | PointMemberLoad | TrapezoidalMemberLoad | TorqueMemberLoad;
+/**
+ * How a null area sends its load to the members along its edges:
+ * - `"twoWay"`: each point to the nearest edge (corner bisectors); every edge needs members.
+ * - `"oneWay"`: strips spanning between `supportEdges`, parallel to the unsupported edges.
+ * - `"cantilever"`: everything to the one edge in `supportEdges`, as a line load plus a torque.
+ */
+export type AreaDistribution = "twoWay" | "oneWay" | "cantilever";
+/**
+ * A load-distribution-only area: it adds no stiffness and only turns area loads into member
+ * loads on the members along its edges. Edge i runs from `nodeIds[i]` to `nodeIds[i + 1]`
+ * (the last back to the first); local axis 3 is the normal by the right-hand rule around the
+ * nodes, so nodes listed counter-clockwise seen from above give local 3 = +Z.
+ */
+export interface NullArea {
+    id: FrameId;
+    nodeIds: FrameId[];
+    distribution: AreaDistribution;
+    /** Edge indices that carry the load: one-way supports, or the one cantilever support. */
+    supportEdges?: number[];
+}
+/** Area load direction: the area's local normal "3" or a global axis. */
+export type AreaLoadDirection = "3" | "X" | "Y" | "Z";
+export interface AreaLoad {
+    loadCase: string;
+    areaId: FrameId;
+    direction: AreaLoadDirection;
+    q: number;
+}
+/** Descriptive category of a load case; it labels the case and does not change the analysis. */
+export type LoadCaseType = "dead" | "superDead" | "live" | "roofLive" | "snow" | "wind" | "seismic" | "other";
+/**
+ * One independently analysed set of loads, e.g. dead, superimposed dead or live. The loads
+ * themselves are listed in `FrameModel.nodalLoads` and `FrameModel.memberLoads`, each tagged
+ * with its `loadCase`.
+ */
+export interface LoadCase {
+    id: string;
+    name?: string;
+    type?: LoadCaseType;
+    /**
+     * Self-weight multiplier (1 = full self-weight). Applies `selfWeight · unitWeight · A` per unit
+     * length of every member along global −Z. Default 0.
+     */
+    selfWeight?: number;
+}
+/** One term of a load combination: a factor times a load case or another combination. */
+export type CombinationTerm = {
+    case: string;
+    factor: number;
+} | {
+    combination: string;
+    factor: number;
+};
+/** A linear combination (sum of factored terms) of load cases and other combinations. */
+export interface LoadCombination {
+    id: string;
+    name?: string;
+    terms: CombinationTerm[];
+}
+export interface FrameModel {
+    nodes: FrameNode[];
+    members: Member[];
+    restraints: Restraint[];
+    loadCases: LoadCase[];
+    /** Nodal loads of every load case, each tagged with its `loadCase`. */
+    nodalLoads?: NodalLoad[];
+    /** Member loads of every load case, each tagged with its `loadCase`. */
+    memberLoads?: MemberLoad[];
+    /** Null areas that distribute `areaLoads` to the members along their edges. */
+    areas?: NullArea[];
+    /** Area loads of every load case, each tagged with its `loadCase`. */
+    areaLoads?: AreaLoad[];
+    combinations?: LoadCombination[];
+}

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

@@ -0,0 +1,3 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+//# sourceMappingURL=model.type.js.map

package/dist/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/fea/frame/types/results.type.d.ts

@@ -0,0 +1,134 @@
+import type { Warnings } from "../../../core/types/output-message.type";
+import type { FrameId, LoadCaseType } from "./model.type";
+export interface NodalDisplacement {
+    nodeId: FrameId;
+    ux: number;
+    uy: number;
+    uz: number;
+    rx: number;
+    ry: number;
+    rz: number;
+}
+export interface NodalReaction {
+    nodeId: FrameId;
+    fx: number;
+    fy: number;
+    fz: number;
+    mx: number;
+    my: number;
+    mz: number;
+}
+/** Force the node exerts on one member end, in local axes 1, 2, 3 (right-hand rotations). */
+export interface MemberEndForce {
+    f1: number;
+    f2: number;
+    f3: number;
+    m1: number;
+    m2: number;
+    m3: number;
+}
+export interface MemberEndForces {
+    memberId: FrameId;
+    start: MemberEndForce;
+    end: MemberEndForce;
+}
+/**
+ * Internal forces and displacements at a distance `x` from the member start. Forces act on
+ * the positive (+1) face of the cut; N is positive in tension. Displacements include the
+ * rigid-body motion of the member, so at the ends they equal the node displacements (except
+ * across a released translation).
+ */
+export interface MemberStation {
+    x: number;
+    N: number;
+    V2: number;
+    V3: number;
+    T: number;
+    M2: number;
+    M3: number;
+    u1: number;
+    u2: number;
+    u3: number;
+    ux: number;
+    uy: number;
+    uz: number;
+}
+/** Station components with extremes: internal forces and local displacements. */
+export type StationComponent = "N" | "V2" | "V3" | "T" | "M2" | "M3" | "u1" | "u2" | "u3";
+export interface ExtremeValue {
+    value: number;
+    x: number;
+}
+export interface MemberDiagram {
+    memberId: FrameId;
+    length: number;
+    stations: MemberStation[];
+    extremes: Record<StationComponent, {
+        max: ExtremeValue;
+        min: ExtremeValue;
+    }>;
+}
+/** Max or min of a result across combinations, with the combination that produces it. */
+export interface GoverningValue {
+    value: number;
+    combinationId: string;
+}
+export interface GoverningExtremeValue extends GoverningValue {
+    x: number;
+}
+export type Bound<T> = {
+    max: T;
+    min: T;
+};
+/** Results of one load case or one load combination. */
+export interface FrameCaseResult {
+    id: string;
+    kind: "case" | "combination";
+    name?: string;
+    type?: LoadCaseType;
+    nodalDisplacements: NodalDisplacement[];
+    reactions: NodalReaction[];
+    memberEndForces: MemberEndForces[];
+    memberDiagrams: MemberDiagram[];
+    equilibriumResidual: {
+        force: number;
+        moment: number;
+    };
+    warnings: Warnings;
+}
+/** Max/min of every result across a set of combinations, with the governing combination. */
+export interface FrameEnvelope {
+    combinationIds: string[];
+    nodalDisplacements: ({
+        nodeId: FrameId;
+    } & Record<"ux" | "uy" | "uz" | "rx" | "ry" | "rz", Bound<GoverningValue>>)[];
+    reactions: ({
+        nodeId: FrameId;
+    } & Record<"fx" | "fy" | "fz" | "mx" | "my" | "mz", Bound<GoverningValue>>)[];
+    members: {
+        memberId: FrameId;
+        extremes: Record<StationComponent, Bound<GoverningExtremeValue>>;
+    }[];
+}
+export interface FrameResult {
+    cases: Record<string, FrameCaseResult>;
+    combinations: Record<string, FrameCaseResult>;
+    /** Envelope over all combinations; absent when the model has none. */
+    envelope?: FrameEnvelope;
+    calculationDetails: {
+        nodeCount: number;
+        memberCount: number;
+        dofCount: number;
+        freeDofCount: number;
+        profileSize: number;
+        loadCaseCount: number;
+        combinationCount: number;
+    };
+    units: {
+        force: "kN";
+        moment: "kN·m";
+        length: "mm";
+        rotation: "rad";
+    };
+    warnings: Warnings;
+}

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

@@ -0,0 +1,3 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+//# sourceMappingURL=results.type.js.map

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

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

package/dist/fea/frame/validation/validate-model.d.ts

@@ -0,0 +1,18 @@
+import { type FrameErrorDetails, type FrameModelIssue } from "../../../core/errors/frame-analysis.error";
+import type { FrameModel } from "../types/model.type";
+/** Builds an issue whose message names where the problem is. */
+export declare const frameIssue: (std: {
+    code: number;
+    message: string;
+}, details?: FrameErrorDetails) => FrameModelIssue;
+/** Throws the first issue as a `FrameAnalysisError` carrying all of them. */
+export declare const throwIssues: (issues: FrameModelIssue[]) => never;
+/**
+ * Checks a frame model without analysing it and returns every problem found, each with the
+ * node, member, load case, combination and field it concerns. An empty array means the
+ * model can be analysed, apart from an unstable structure (304), which only the solve detects.
+ *
+ * @param model - The model to check
+ * @returns Every problem found (codes 301–303 and 305–317), in model order
+ */
+export declare const validateFrameModel: (model: FrameModel) => FrameModelIssue[];