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

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;

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

@@ -0,0 +1,124 @@
+"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);
+    };
+    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`
+ */
+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) => load.type === "udl" ? { ...load, w: load.w * factor } : { ...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,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;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,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;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,112 @@
+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;
+}
+export type MemberLoad = UniformMemberLoad | PointMemberLoad;
+/** 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[];
+    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":""}