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

package/dist/esm/core/errors/frame-analysis.error.js

@@ -0,0 +1,70 @@
+export class FrameAnalysisError extends Error {
+    constructor(std = { code: 0, message: "Unknown error" }, details = {}, issues) {
+        super(std.message);
+        this.code = std.code;
+        this.message = std.message;
+        this.details = details;
+        this.issues = issues ?? [{ code: std.code, message: std.message, details }];
+        Object.setPrototypeOf(this, FrameAnalysisError.prototype);
+    }
+}
+export const FrameErrors = {
+    // 3XX : Frame analysis errors
+    ZERO_LENGTH_MEMBER: {
+        code: 301,
+        message: "A member's start and end nodes coincide; the member has zero length.",
+    },
+    UNKNOWN_NODE_REFERENCE: {
+        code: 302,
+        message: "A member, restraint or load refers to a node or member id that is not in the model.",
+    },
+    DUPLICATE_ID: {
+        code: 303,
+        message: "Node ids, member ids and restrained node ids must each be unique, and load case and combination ids must be unique across both lists.",
+    },
+    STRUCTURE_UNSTABLE: {
+        code: 304,
+        message: "The restrained structure is unstable (singular stiffness matrix); add supports or members to prevent a mechanism.",
+    },
+    INVALID_SECTION_PROPERTIES: {
+        code: 305,
+        message: "Member section properties E, G, A, I22, I33 and J must all be positive, and unitWeight, when given, must not be negative.",
+    },
+    NO_RESTRAINTS: {
+        code: 306,
+        message: "The model has no restrained degrees of freedom.",
+    },
+    INVALID_RELEASE: {
+        code: 307,
+        message: "The end releases on a member make it a mechanism (rigid-body motion within the member).",
+    },
+    LOAD_OUTSIDE_MEMBER: {
+        code: 308,
+        message: "A member load lies outside the member, or its start is not before its end.",
+    },
+    UNKNOWN_LOAD_CASE: {
+        code: 309,
+        message: "A load or a load combination term refers to a load case or combination that is not in the model, or a term names neither.",
+    },
+    CIRCULAR_COMBINATION: {
+        code: 310,
+        message: "Load combinations reference each other in a loop.",
+    },
+    INVALID_LOAD_FACTOR: {
+        code: 311,
+        message: "A combination factor or self-weight multiplier is not a finite number, or a load combination has no terms.",
+    },
+    NO_LOAD_CASES: {
+        code: 312,
+        message: "The model has no load cases; define at least one in loadCases.",
+    },
+    MISSING_UNIT_WEIGHT: {
+        code: 313,
+        message: "A load case requests self-weight but a member section has no unitWeight; set unitWeight (0 for a weightless member).",
+    },
+    INVALID_INPUT_VALUE: {
+        code: 314,
+        message: "An input value is invalid: a number that is not finite (NaN, Infinity or not a number), a negative spring stiffness, or a value outside its allowed options.",
+    },
+};
+//# sourceMappingURL=frame-analysis.error.js.map

package/dist/esm/core/errors/frame-analysis.error.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"frame-analysis.error.js","sourceRoot":"","sources":["../../../../src/core/errors/frame-analysis.error.ts"],"names":[],"mappings":"AAcA,MAAM,OAAO,kBAAmB,SAAQ,KAAK;IAO3C,YACE,MAA8B,EAAE,IAAI,EAAE,CAAC,EAAE,OAAO,EAAE,eAAe,EAAE,EACnE,UAA6B,EAAE,EAC/B,MAA0B;QAE1B,KAAK,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QACnB,IAAI,CAAC,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;QACrB,IAAI,CAAC,OAAO,GAAG,GAAG,CAAC,OAAO,CAAC;QAC3B,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,MAAM,GAAG,MAAM,IAAI,CAAC,EAAE,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC,CAAC;QAC5E,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,kBAAkB,CAAC,SAAS,CAAC,CAAC;IAC5D,CAAC;CACF;AAcD,MAAM,CAAC,MAAM,WAAW,GAA2C;IACjE,8BAA8B;IAC9B,kBAAkB,EAAE;QAClB,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,sEAAsE;KAChF;IACD,sBAAsB,EAAE;QACtB,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,qFAAqF;KAC/F;IACD,YAAY,EAAE;QACZ,IAAI,EAAE,GAAG;QACT,OAAO,EACL,uIAAuI;KAC1I;IACD,kBAAkB,EAAE;QAClB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,mHAAmH;KACtH;IACD,0BAA0B,EAAE;QAC1B,IAAI,EAAE,GAAG;QACT,OAAO,EACL,2HAA2H;KAC9H;IACD,aAAa,EAAE;QACb,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,iDAAiD;KAC3D;IACD,eAAe,EAAE;QACf,IAAI,EAAE,GAAG;QACT,OAAO,EACL,yFAAyF;KAC5F;IACD,mBAAmB,EAAE;QACnB,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,4EAA4E;KACtF;IACD,iBAAiB,EAAE;QACjB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,2HAA2H;KAC9H;IACD,oBAAoB,EAAE;QACpB,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,mDAAmD;KAC7D;IACD,mBAAmB,EAAE;QACnB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,4GAA4G;KAC/G;IACD,aAAa,EAAE;QACb,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,gEAAgE;KAC1E;IACD,mBAAmB,EAAE;QACnB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,sHAAsH;KACzH;IACD,mBAAmB,EAAE;QACnB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,8JAA8J;KACjK;CACF,CAAC"}

package/dist/esm/core/errors/rc-design.error.js

@@ -0,0 +1,64 @@
+export class RCDesignError extends Error {
+    constructor(std = { code: 0, message: "Unknown error" }) {
+        super(std.message);
+        this.code = std.code;
+        this.message = std.message;
+        Object.setPrototypeOf(this, RCDesignError.prototype);
+    }
+}
+export const Errors = {
+    // 1XX : RC Beam Design Errors
+    /**
+     * @deprecated Never thrown since 1.4.0. A section whose compression steel
+     * lies below the neutral axis is valid — that steel acts in tension and
+     * `rectBeamMomentCapacity` now returns a reduced Mn with a warning instead.
+     * Retained so existing `Errors.COMP_STEEL_DISTANCE_OVER_AC` references and
+     * `code === 101` checks keep compiling; safe to delete in the next major.
+     */
+    COMP_STEEL_DISTANCE_OVER_AC: {
+        code: 101,
+        message: "d' of compressive steel is greater than c, consider using singly reinforced section.",
+    },
+    TENSILE_STEEL_NOT_YIELDING: {
+        code: 102,
+        message: "Out of tension controlled region, tensile steel is not yielding.",
+    },
+    SECTION_TOO_NARROW_FOR_BAR: {
+        code: 103,
+        message: "Section width is too narrow to fit even one bar of the given diameter with the required cover and spacing.",
+    },
+    BAR_LAYOUT_EXCEEDS_SECTION_HEIGHT: {
+        code: 104,
+        message: "Combined top and bottom rebar layers exceed the section height h.",
+    },
+    SHEAR_STRENGTH_EXCEEDS_MAX: {
+        code: 105,
+        message: "Vs exceeds 0.66*sqrt(fc')*bw*d (ACI318-19, 22.5.1.2); section is too small for shear, increase bw or d.",
+    },
+    // 2XX : Strengthening / plate-to-concrete connection errors
+    NON_POSITIVE_SHEAR_FLOW: {
+        code: 201,
+        message: "Interface shear flow requires a positive shear force V and at least one jacketing plate with non-zero area.",
+    },
+    INVALID_BOLT_CAPACITY: {
+        code: 202,
+        message: "Resolved bolt shear capacity is not positive; check Ase, futa, or the supplied phiVboltOverride.",
+    },
+    PLATE_OUTSIDE_SECTION: {
+        code: 203,
+        message: "Inconsistent jacketed section geometry; plate and section dimensions must be positive and d must lie within h.",
+    },
+    CRACKED_NA_NOT_FOUND: {
+        code: 204,
+        message: "No valid cracked neutral axis was found within the concrete depth; check the section and plate properties.",
+    },
+    SIDE_PLATE_GEOMETRY_INVALID: {
+        code: 205,
+        message: "Inconsistent side plate geometry; thickness, engaged depth, fy and sides must be positive, the depth must fit within the section, and a stripped layout requires a positive width and spacing.",
+    },
+    STRENGTHENED_SHEAR_NO_HEADROOM: {
+        code: 206,
+        message: "Existing Vs already reaches 0.66*sqrt(fc')*bw*d (ACI318-19, 22.5.1.2); the web crushes before any added side plate force can be developed, so shear strengthening cannot be realised on this section.",
+    },
+};
+//# sourceMappingURL=rc-design.error.js.map

package/dist/esm/core/errors/rc-design.error.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"rc-design.error.js","sourceRoot":"","sources":["../../../../src/core/errors/rc-design.error.ts"],"names":[],"mappings":"AAAA,MAAM,OAAO,aAAc,SAAQ,KAAK;IAItC,YAAY,MAAyB,EAAE,IAAI,EAAE,CAAC,EAAE,OAAO,EAAE,eAAe,EAAE;QACxE,KAAK,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QACnB,IAAI,CAAC,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;QACrB,IAAI,CAAC,OAAO,GAAG,GAAG,CAAC,OAAO,CAAC;QAC3B,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,aAAa,CAAC,SAAS,CAAC,CAAC;IACvD,CAAC;CACF;AAOD,MAAM,CAAC,MAAM,MAAM,GAAsC;IACvD,8BAA8B;IAE9B;;;;;;OAMG;IACH,2BAA2B,EAAE;QAC3B,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,sFAAsF;KAChG;IACD,0BAA0B,EAAE;QAC1B,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,kEAAkE;KAC5E;IACD,0BAA0B,EAAE;QAC1B,IAAI,EAAE,GAAG;QACT,OAAO,EACL,4GAA4G;KAC/G;IACD,iCAAiC,EAAE;QACjC,IAAI,EAAE,GAAG;QACT,OAAO,EAAE,mEAAmE;KAC7E;IACD,0BAA0B,EAAE;QAC1B,IAAI,EAAE,GAAG;QACT,OAAO,EACL,yGAAyG;KAC5G;IAED,4DAA4D;IAC5D,uBAAuB,EAAE;QACvB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,6GAA6G;KAChH;IACD,qBAAqB,EAAE;QACrB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,kGAAkG;KACrG;IACD,qBAAqB,EAAE;QACrB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,gHAAgH;KACnH;IACD,oBAAoB,EAAE;QACpB,IAAI,EAAE,GAAG;QACT,OAAO,EACL,4GAA4G;KAC/G;IACD,2BAA2B,EAAE;QAC3B,IAAI,EAAE,GAAG;QACT,OAAO,EACL,gMAAgM;KACnM;IACD,8BAA8B,EAAE;QAC9B,IAAI,EAAE,GAAG;QACT,OAAO,EACL,uMAAuM;KAC1M;CACF,CAAC"}

package/dist/esm/core/index.js

@@ -0,0 +1,9 @@
+export { isSinglyReinforced } from "./types/rc-beam.type.js";
+export { hasTopPlate, hasBottomPlate, hasSidePlate } from "./types/plate-jacketing.type.js";
+// Constants
+export { FLEXURAL_STRENGTH_REDUCTION_FACTOR, CONCRETE_ULTIMATE_STRAIN, SHEAR_STRENGTH_REDUCTION_FACTOR, ANCHOR_SHEAR_STRENGTH_REDUCTION_FACTOR, MAX_PLATE_BOLT_SPACING, PLATE_SHEAR_YIELD_COEFFICIENT, SIDE_PLATE_EFFECTIVE_STRAIN_LIMIT, SIDE_PLATE_REDUCTION_FACTOR_PSI, } from "./constants/rc.constant.js";
+export { CONCRETE_POISSON_RATIO, CONCRETE_UNIT_WEIGHT, BEAM_INERTIA_MODIFIER, COLUMN_INERTIA_MODIFIER, VERTICAL_MEMBER_TOLERANCE, SINGULARITY_TOLERANCE, EQUILIBRIUM_TOLERANCE, } from "./constants/frame.constant.js";
+// Errors
+export { RCDesignError, Errors } from "./errors/rc-design.error.js";
+export { FrameAnalysisError, FrameErrors } from "./errors/frame-analysis.error.js";
+//# sourceMappingURL=index.js.map

package/dist/esm/core/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/core/index.ts"],"names":[],"mappings":"AAWA,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAW1D,OAAO,EAAE,WAAW,EAAE,cAAc,EAAE,YAAY,EAAE,MAAM,8BAA8B,CAAC;AAKzF,YAAY;AACZ,OAAO,EACL,kCAAkC,EAClC,wBAAwB,EACxB,+BAA+B,EAC/B,sCAAsC,EACtC,sBAAsB,EACtB,6BAA6B,EAC7B,iCAAiC,EACjC,+BAA+B,GAChC,MAAM,yBAAyB,CAAC;AAEjC,OAAO,EACL,sBAAsB,EACtB,oBAAoB,EACpB,qBAAqB,EACrB,uBAAuB,EACvB,yBAAyB,EACzB,qBAAqB,EACrB,qBAAqB,GACtB,MAAM,4BAA4B,CAAC;AAEpC,SAAS;AACT,OAAO,EAAE,aAAa,EAAE,MAAM,EAAE,MAAM,0BAA0B,CAAC;AACjE,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,+BAA+B,CAAC"}

package/dist/esm/core/types/output-message.type.js

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

package/dist/esm/core/types/output-message.type.js.map

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

package/dist/esm/core/types/plate-jacketing.type.js

@@ -0,0 +1,18 @@
+/**
+ * A plate counts as present only when both its width and thickness are positive,
+ * so a bottom-only jacket may omit the top plate fields or pass them as zero.
+ */
+export function hasTopPlate(props) {
+    return (props.topSteelWidth ?? 0) > 0 && (props.topSteelThickness ?? 0) > 0;
+}
+export function hasBottomPlate(props) {
+    return (props.bottomSteelWidth ?? 0) > 0 && (props.bottomSteelThickness ?? 0) > 0;
+}
+/**
+ * A side plate counts as present only when both its thickness and engaged depth
+ * are positive.
+ */
+export function hasSidePlate(props) {
+    return props.thickness > 0 && props.depth > 0;
+}
+//# sourceMappingURL=plate-jacketing.type.js.map

package/dist/esm/core/types/plate-jacketing.type.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"plate-jacketing.type.js","sourceRoot":"","sources":["../../../../src/core/types/plate-jacketing.type.ts"],"names":[],"mappings":"AAWA;;;GAGG;AACH,MAAM,UAAU,WAAW,CAAC,KAAyB;IACnD,OAAO,CAAC,KAAK,CAAC,aAAa,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,iBAAiB,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;AAC9E,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,KAAyB;IACtD,OAAO,CAAC,KAAK,CAAC,gBAAgB,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,oBAAoB,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;AACpF,CAAC;AA6DD;;;GAGG;AACH,MAAM,UAAU,YAAY,CAAC,KAA0B;IACrD,OAAO,KAAK,CAAC,SAAS,GAAG,CAAC,IAAI,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC;AAChD,CAAC"}

package/dist/esm/core/types/rc-beam.type.js

@@ -0,0 +1,11 @@
+/**
+ * A section counts as singly reinforced when it carries no compression steel.
+ *
+ * `d_` is deliberately not sign-checked: it is measured from the extreme
+ * concrete compression fibre, so compression steel above that face — a steel
+ * plate bonded to the top of the beam — legitimately has `d_ <= 0`.
+ */
+export function isSinglyReinforced(section) {
+    return !("As_" in section) || !("d_" in section) || section.As_ <= 0;
+}
+//# sourceMappingURL=rc-beam.type.js.map

package/dist/esm/core/types/rc-beam.type.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"rc-beam.type.js","sourceRoot":"","sources":["../../../../src/core/types/rc-beam.type.ts"],"names":[],"mappings":"AAoDA;;;;;;GAMG;AACH,MAAM,UAAU,kBAAkB,CAAC,OAAwB;IACzD,OAAO,CAAC,CAAC,KAAK,IAAI,OAAO,CAAC,IAAI,CAAC,CAAC,IAAI,IAAI,OAAO,CAAC,IAAI,OAAO,CAAC,GAAG,IAAI,CAAC,CAAC;AACvE,CAAC"}

package/dist/esm/core/types/section-properties.type.js

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

package/dist/esm/core/types/section-properties.type.js.map

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

package/dist/esm/core/types/unit.type.js

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

package/dist/esm/core/types/unit.type.js.map

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

package/dist/esm/fea/frame/elements/beam-3d-stiffness.js

@@ -0,0 +1,84 @@
+import { zeros } from "../../../utils/matrix.js";
+/**
+ * Local DOF indices (0-based) whose reference-matrix rotation sign is opposite to the
+ * right-hand rule: R2 and R3 at the start and end nodes.
+ */
+const FLIPPED_DOFS = [4, 5, 10, 11];
+/**
+ * Local stiffness matrix of a 3D prismatic beam element, exactly as in the reference
+ * "Complete Beam Stiffness Matrix".
+ *
+ * DOF order for node i then node j: [U1, U2, U3, R1, R2, R3]. Local 1 runs along the member.
+ * The reference matrix measures the bending rotations R2 and R3 opposite to the right-hand
+ * rule; use `beam3dLocalStiffnessRightHanded` for the right-handed form.
+ *
+ * The reference uses one `EI`; here the U2–R3 terms use `E·I33` and the U3–R2 terms use
+ * `E·I22`, which reduces to the reference when `I22 = I33`.
+ *
+ * @param E - Young's modulus in MPa
+ * @param G - Shear modulus in MPa
+ * @param A - Area in mm²
+ * @param I22 - Second moment about local 2 in mm⁴
+ * @param I33 - Second moment about local 3 in mm⁴
+ * @param J - Torsion constant in mm⁴
+ * @param L - Member length in mm
+ * @returns 12×12 symmetric matrix in N and mm
+ */
+export const beam3dLocalStiffness = (E, G, A, I22, I33, J, L) => {
+    const k = zeros(12, 12);
+    // 1-based row/column, upper triangle as written in the reference; mirrored for symmetry.
+    const set = (row, col, value) => {
+        k[row - 1][col - 1] = value;
+        k[col - 1][row - 1] = value;
+    };
+    const ax = (E * A) / L;
+    const tor = (G * J) / L;
+    // U2–R3 bending (I33)
+    const a3 = (12 * E * I33) / L ** 3;
+    const b3 = (6 * E * I33) / L ** 2;
+    const c3 = (4 * E * I33) / L;
+    const d3 = (2 * E * I33) / L;
+    // U3–R2 bending (I22)
+    const a2 = (12 * E * I22) / L ** 3;
+    const b2 = (6 * E * I22) / L ** 2;
+    const c2 = (4 * E * I22) / L;
+    const d2 = (2 * E * I22) / L;
+    set(1, 1, ax);
+    set(1, 7, -ax);
+    set(7, 7, ax);
+    set(2, 2, a3);
+    set(2, 6, -b3);
+    set(2, 8, -a3);
+    set(2, 12, -b3);
+    set(6, 6, c3);
+    set(6, 8, b3);
+    set(6, 12, d3);
+    set(8, 8, a3);
+    set(8, 12, b3);
+    set(12, 12, c3);
+    set(3, 3, a2);
+    set(3, 5, b2);
+    set(3, 9, -a2);
+    set(3, 11, b2);
+    set(5, 5, c2);
+    set(5, 9, -b2);
+    set(5, 11, d2);
+    set(9, 9, a2);
+    set(9, 11, -b2);
+    set(11, 11, c2);
+    set(4, 4, tor);
+    set(4, 10, -tor);
+    set(10, 10, tor);
+    return k;
+};
+/**
+ * Converts a local matrix between the reference rotation convention and the right-hand rule
+ * (`D·k·D`, with `D` flipping the sign of R2 and R3 at both ends). The conversion is its own inverse.
+ */
+export const flipBendingRotationSigns = (k) => {
+    const flipped = new Set(FLIPPED_DOFS);
+    return k.map((row, i) => row.map((value, j) => (flipped.has(i) !== flipped.has(j) ? -value : value)));
+};
+/** Same element as `beam3dLocalStiffness` with all rotations following the right-hand rule. */
+export const beam3dLocalStiffnessRightHanded = (E, G, A, I22, I33, J, L) => flipBendingRotationSigns(beam3dLocalStiffness(E, G, A, I22, I33, J, L));
+//# sourceMappingURL=beam-3d-stiffness.js.map

package/dist/esm/fea/frame/elements/beam-3d-stiffness.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"beam-3d-stiffness.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/beam-3d-stiffness.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAe,MAAM,mBAAmB,CAAC;AAEvD;;;GAGG;AACH,MAAM,YAAY,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;AAEpC;;;;;;;;;;;;;;;;;;;GAmBG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAG,CAClC,CAAS,EACT,CAAS,EACT,CAAS,EACT,GAAW,EACX,GAAW,EACX,CAAS,EACT,CAAS,EACD,EAAE;IACV,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,yFAAyF;IACzF,MAAM,GAAG,GAAG,CAAC,GAAW,EAAE,GAAW,EAAE,KAAa,EAAE,EAAE;QACtD,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QAC5B,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;IAC9B,CAAC,CAAC;IAEF,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IAExB,sBAAsB;IACtB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,sBAAsB;IACtB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAE7B,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IAEd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAEhB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAEhB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACjB,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IAEjB,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,wBAAwB,GAAG,CAAC,CAAS,EAAU,EAAE;IAC5D,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,YAAY,CAAC,CAAC;IACtC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CACtB,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAC5E,CAAC;AACJ,CAAC,CAAC;AAEF,+FAA+F;AAC/F,MAAM,CAAC,MAAM,+BAA+B,GAAG,CAC7C,CAAS,EACT,CAAS,EACT,CAAS,EACT,GAAW,EACX,GAAW,EACX,CAAS,EACT,CAAS,EACD,EAAE,CAAC,wBAAwB,CAAC,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC"}

package/dist/esm/fea/frame/elements/fixed-end-forces.js

@@ -0,0 +1,75 @@
+import { FrameAnalysisError, FrameErrors } from "../../../core/errors/frame-analysis.error.js";
+// 4-point Gauss–Legendre rule, exact for the cubic integrands used here
+const GAUSS_POINTS = [
+    -0.8611363115940526, -0.3399810435848563, 0.3399810435848563, 0.8611363115940526,
+];
+const GAUSS_WEIGHTS = [
+    0.3478548451374538, 0.6521451548625461, 0.6521451548625461, 0.3478548451374538,
+];
+/**
+ * Equivalent nodal loads of a unit transverse/axial load at distance `x`, from the
+ * Euler–Bernoulli shape functions: the exact fixed-end solution for a prismatic member.
+ * Returned per unit of load, right-hand local DOF order [i: U1..R3, j: U1..R3].
+ */
+const unitEquivalentLoad = (axis, x, L) => {
+    const xi = x / L;
+    const out = new Array(12).fill(0);
+    if (axis === 1) {
+        out[0] = 1 - xi;
+        out[6] = xi;
+        return out;
+    }
+    const n1 = 1 - 3 * xi ** 2 + 2 * xi ** 3;
+    const n2 = L * (xi - 2 * xi ** 2 + xi ** 3);
+    const n3 = 3 * xi ** 2 - 2 * xi ** 3;
+    const n4 = L * (-(xi ** 2) + xi ** 3);
+    if (axis === 2) {
+        // load along local 2: force U2, moment about 3 (R3 = +dv/dx)
+        out[1] = n1;
+        out[5] = n2;
+        out[7] = n3;
+        out[11] = n4;
+    }
+    else {
+        // load along local 3: force U3, moment about 2 (R2 = −dw/dx)
+        out[2] = n1;
+        out[4] = -n2;
+        out[8] = n3;
+        out[10] = -n4;
+    }
+    return out;
+};
+/**
+ * Fixed-end forces of a member under local loads: the end forces (right-hand local axes,
+ * N and N·mm) that hold both ends fully fixed. Equivalent nodal loads are the negative.
+ *
+ * @param loads - Loads in local axes, N and mm
+ * @param L - Member length in mm
+ * @throws {FrameAnalysisError} 308 if a load lies outside the member
+ */
+export const fixedEndForces = (loads, L) => {
+    const equivalent = new Array(12).fill(0);
+    const add = (unit, factor) => {
+        for (let i = 0; i < 12; i++)
+            equivalent[i] += unit[i] * factor;
+    };
+    for (const load of loads) {
+        if (load.type === "point") {
+            if (!(load.a >= 0 && load.a <= L))
+                throw new FrameAnalysisError(FrameErrors.LOAD_OUTSIDE_MEMBER);
+            add(unitEquivalentLoad(load.axis, load.a, L), load.P);
+        }
+        else {
+            if (!(load.start >= 0 && load.end <= L && load.start < load.end)) {
+                throw new FrameAnalysisError(FrameErrors.LOAD_OUTSIDE_MEMBER);
+            }
+            const half = (load.end - load.start) / 2;
+            const mid = (load.end + load.start) / 2;
+            GAUSS_POINTS.forEach((p, index) => {
+                add(unitEquivalentLoad(load.axis, mid + half * p, L), load.w * half * GAUSS_WEIGHTS[index]);
+            });
+        }
+    }
+    return equivalent.map((v) => -v);
+};
+//# sourceMappingURL=fixed-end-forces.js.map

package/dist/esm/fea/frame/elements/fixed-end-forces.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"fixed-end-forces.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/fixed-end-forces.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AAOxF,wEAAwE;AACxE,MAAM,YAAY,GAAG;IACnB,CAAC,kBAAkB,EAAE,CAAC,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CACjF,CAAC;AACF,MAAM,aAAa,GAAG;IACpB,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CAC/E,CAAC;AAEF;;;;GAIG;AACH,MAAM,kBAAkB,GAAG,CAAC,IAAe,EAAE,CAAS,EAAE,CAAS,EAAY,EAAE;IAC7E,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QAChB,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,OAAO,GAAG,CAAC;IACb,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IAC5C,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACrC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IACtC,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;IACf,CAAC;SAAM,CAAC;QACN,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QACb,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC;IAChB,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,KAAwB,EAAE,CAAS,EAAY,EAAE;IAC9E,MAAM,UAAU,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,GAAG,GAAG,CAAC,IAAc,EAAE,MAAc,EAAE,EAAE;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,UAAU,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;IACjE,CAAC,CAAC;IAEF,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;YAC1B,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;gBAC/B,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,GAAG,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACjE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,CAAC;YACD,MAAM,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACxC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE;gBAChC,GAAG,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,GAAG,GAAG,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;YAC9F,CAAC,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,OAAO,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AACnC,CAAC,CAAC"}

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

@@ -0,0 +1,55 @@
+import { FrameAnalysisError, FrameErrors } from "../../../core/errors/frame-analysis.error.js";
+import { solveLinearSystem } from "../../../utils/matrix.js";
+/**
+ * Statically condenses released local DOFs out of a member, so those DOFs carry no force.
+ *
+ * With released set r and retained set c:
+ * `k* = k_cc − k_cr·k_rr⁻¹·k_rc`, `f* = f_c − k_cr·k_rr⁻¹·f_r`.
+ * The stiffness is condensed once; the returned function condenses the fixed-end forces of
+ * each load case.
+ *
+ * @param k - 12×12 local stiffness (right-hand)
+ * @param released - Local DOF indices (0–11) to release
+ * @throws {FrameAnalysisError} 307 if the released set makes the member a mechanism
+ */
+export const condenseMember = (k, released) => {
+    if (released.length === 0) {
+        return { k, condenseFixedEndForces: (fixed) => [...fixed] };
+    }
+    const retained = Array.from({ length: 12 }, (_, i) => i).filter((i) => !released.includes(i));
+    const sub = (rows, cols) => rows.map((r) => cols.map((c) => k[r][c]));
+    const krr = sub(released, released);
+    const krc = sub(released, retained);
+    const kcr = sub(retained, released);
+    const kcc = sub(retained, retained);
+    // [k_rr⁻¹·k_rc | k_rr⁻¹]
+    const rhs = released.map((_, i) => [...krc[i], ...released.map((__, j) => (i === j ? 1 : 0))]);
+    const solved = solveLinearSystem(krr, rhs);
+    if (!solved)
+        throw new FrameAnalysisError(FrameErrors.INVALID_RELEASE);
+    const nc = retained.length;
+    const kStar = k.map((row) => row.map(() => 0));
+    retained.forEach((r, ri) => {
+        retained.forEach((c, ci) => {
+            let value = kcc[ri][ci];
+            released.forEach((_, j) => {
+                value -= kcr[ri][j] * solved[j][ci];
+            });
+            kStar[r][c] = value;
+        });
+    });
+    // k_cr·k_rr⁻¹, applied to f_r for every load case
+    const kcrKrrInv = retained.map((_, ri) => released.map((__, m) => released.reduce((s, ___, j) => s + kcr[ri][j] * solved[j][nc + m], 0)));
+    return {
+        k: kStar,
+        condenseFixedEndForces: (fixed) => {
+            const fStar = new Array(12).fill(0);
+            const fr = released.map((r) => fixed[r]);
+            retained.forEach((r, ri) => {
+                fStar[r] = fixed[r] - kcrKrrInv[ri].reduce((s, v, m) => s + v * fr[m], 0);
+            });
+            return fStar;
+        },
+    };
+};
+//# sourceMappingURL=static-condensation.js.map

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

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

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

@@ -0,0 +1,6 @@
+export { analyzeFrame, createFrameAnalysis } from "./solver/analyze-frame.js";
+export { validateFrameModel } from "./validation/validate-model.js";
+export { beam3dLocalStiffness, beam3dLocalStiffnessRightHanded, } from "./elements/beam-3d-stiffness.js";
+export { memberAxes } from "./transformation/coordinate-transform.js";
+export { memberDemandEnvelope } from "./post/demand-envelope.js";
+//# sourceMappingURL=index.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../src/fea/frame/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,YAAY,EAAE,mBAAmB,EAAE,MAAM,wBAAwB,CAAC;AAE3E,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,OAAO,EACL,oBAAoB,EACpB,+BAA+B,GAChC,MAAM,8BAA8B,CAAC;AACtC,OAAO,EAAE,UAAU,EAAE,MAAM,uCAAuC,CAAC;AACnE,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC"}

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

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

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

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