@theerapat-s28/ts-structural-eng-tools 3.1.0 → 3.2.0

package/dist/fea/frame/solver/{global-assembly.d.ts → member-assembly.d.ts}

@@ -1,21 +1,10 @@
 import { type Matrix } from "../../../utils/matrix";
-/**
- * A stiffness contribution to the global system: the nodes it connects and its stiffness over
- * their DOFs (6 per node, in node order) in global axes. A member is a 2-node element; elements
- * with more nodes assemble the same way.
- */
-export interface StiffnessElement {
-    nodes: number[];
-    kGlobal: Matrix;
-}
 export interface AssemblyMember {
     startIndex: number;
     endIndex: number;
     T: Matrix;
     k: Matrix;
 }
-/** Global DOF numbers (6 per node, in node order) of an element connecting `nodes`. */
-export declare const elementDofs: (nodes: readonly number[]) => number[];
 /** 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). */

package/dist/fea/frame/solver/{global-assembly.js → member-assembly.js}

@@ -1,19 +1,10 @@
 "use strict";
 Object.defineProperty(exports, "__esModule", { value: true });
-exports.scatterMemberForces = exports.memberGlobalStiffness = exports.memberDofs = exports.elementDofs = void 0;
+exports.scatterMemberForces = exports.memberGlobalStiffness = exports.memberDofs = void 0;
+const global_assembly_1 = require("../../shared/solver/global-assembly");
 const matrix_1 = require("../../../utils/matrix");
-/** Global DOF numbers (6 per node, in node order) of an element connecting `nodes`. */
-const elementDofs = (nodes) => {
-    const dofs = [];
-    for (const node of nodes) {
-        for (let k = 0; k < 6; k++)
-            dofs.push(6 * node + k);
-    }
-    return dofs;
-};
-exports.elementDofs = elementDofs;
 /** Global DOF numbers (6 per node) of a member's 12 end DOFs. */
-const memberDofs = (member) => (0, exports.elementDofs)([member.startIndex, member.endIndex]);
+const memberDofs = (member) => (0, global_assembly_1.elementDofs)([member.startIndex, member.endIndex]);
 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));
@@ -33,4 +24,4 @@ const scatterMemberForces = (target, member, localForces, sign) => {
     }
 };
 exports.scatterMemberForces = scatterMemberForces;
-//# sourceMappingURL=global-assembly.js.map
+//# sourceMappingURL=member-assembly.js.map

package/dist/fea/frame/solver/member-assembly.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"member-assembly.js","sourceRoot":"","sources":["../../../../src/fea/frame/solver/member-assembly.ts"],"names":[],"mappings":";;;AAAA,4EAAqE;AACrE,8CAAmE;AASnE,iEAAiE;AAC1D,MAAM,UAAU,GAAG,CAAC,MAAuD,EAAY,EAAE,CAC9F,IAAA,6BAAW,EAAC,CAAC,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;AADvC,QAAA,UAAU,cAC6B;AAEpD,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/member-loads.d.ts

@@ -0,0 +1,21 @@
+import type { FeaId } from "../../shared/types/model.type";
+import { type LocalMemberLoad } from "../elements/fixed-end-forces";
+import type { MemberLoad } from "../types/model.type";
+import type { MemberGroup, SystemMember } from "./system-member";
+/**
+ * The local loads (N and mm) on every member piece in one load case: each member load, and the
+ * self-weight `selfWeight · unitWeight · A` per unit length along global −Z, resolved in the
+ * whole member's local axes, then cut at its pieces.
+ *
+ * @param groups - Model members and the pieces each is analysed as
+ * @param pieces - Every member piece
+ * @param memberIndex - Index into `groups` of each member id
+ * @param loads - The load case's member loads
+ * @param loadCaseId - The load case, which tags the self-weight load
+ * @param selfWeight - The load case's self-weight multiplier
+ * @returns The local loads of each piece, in the order of `pieces`
+ * @throws {FeaError} 313 if self-weight is requested for a member without `unitWeight`
+ */
+export declare const memberPieceLoads: (groups: MemberGroup[], pieces: SystemMember[], memberIndex: Map<FeaId, number>, loads: MemberLoad[], loadCaseId: string, selfWeight: number) => LocalMemberLoad[][];
+/** A local member load scaled by `factor`, e.g. for a load combination. */
+export declare const scaleMemberLoad: (load: LocalMemberLoad, factor: number) => LocalMemberLoad;

package/dist/fea/frame/solver/member-loads.js

@@ -0,0 +1,146 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.scaleMemberLoad = exports.memberPieceLoads = void 0;
+const fea_error_1 = require("../../../core/errors/fea.error");
+const fixed_end_forces_1 = require("../elements/fixed-end-forces");
+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 }];
+    }
+    // per projected length → per member length: Lproj / L = |local 1 × d| = √(1 − (local 1 · d)²)
+    let projection = 1;
+    if (load.type !== "point" && load.projected) {
+        const v = GLOBAL_DIRECTIONS[load.direction];
+        const c = lambda[0][0] * v[0] + lambda[0][1] * v[1] + lambda[0][2] * v[2];
+        projection = Math.sqrt(Math.max(0, 1 - c * c));
+    }
+    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 * projection, // kN/m = N/mm
+            w2: w2 * factor * projection,
+            start: load.start ?? 0,
+            end: load.end ?? L,
+        };
+    });
+};
+/**
+ * The part of a member's local loads that falls on one of its pieces, `[offset, offset + L]`,
+ * in the piece's own positions. A distributed load or torque is cut at the piece ends; a point
+ * load goes to the piece it lies in (at an interior node, to the piece that starts there).
+ */
+const clipLoads = (loads, offset, L, last) => {
+    const end = offset + L;
+    const tolerance = 1e-9 * end;
+    const out = [];
+    for (const load of loads) {
+        if (load.type === "point") {
+            const inside = load.a >= offset - tolerance && (load.a < end - tolerance || last);
+            if (inside)
+                out.push({ ...load, a: Math.min(L, Math.max(0, load.a - offset)) });
+            continue;
+        }
+        const from = Math.max(load.start, offset);
+        const to = Math.min(load.end, end);
+        if (to - from <= tolerance)
+            continue;
+        if (load.type === "torque") {
+            out.push({ ...load, start: from - offset, end: to - offset });
+        }
+        else {
+            out.push({
+                ...load,
+                w1: (0, fixed_end_forces_1.loadIntensity)(load, from),
+                w2: (0, fixed_end_forces_1.loadIntensity)(load, to),
+                start: from - offset,
+                end: to - offset,
+            });
+        }
+    }
+    return out;
+};
+/**
+ * The local loads (N and mm) on every member piece in one load case: each member load, and the
+ * self-weight `selfWeight · unitWeight · A` per unit length along global −Z, resolved in the
+ * whole member's local axes, then cut at its pieces.
+ *
+ * @param groups - Model members and the pieces each is analysed as
+ * @param pieces - Every member piece
+ * @param memberIndex - Index into `groups` of each member id
+ * @param loads - The load case's member loads
+ * @param loadCaseId - The load case, which tags the self-weight load
+ * @param selfWeight - The load case's self-weight multiplier
+ * @returns The local loads of each piece, in the order of `pieces`
+ * @throws {FeaError} 313 if self-weight is requested for a member without `unitWeight`
+ */
+const memberPieceLoads = (groups, pieces, memberIndex, loads, loadCaseId, selfWeight) => {
+    const loadsByMember = groups.map(() => []);
+    for (const load of loads) {
+        loadsByMember[memberIndex.get(load.memberId)].push(load);
+    }
+    // local loads of each model member, then of each piece it is analysed as
+    const memberLoads = pieces.map(() => []);
+    groups.forEach((group, g) => {
+        const memberLoadsHere = [...loadsByMember[g]];
+        if (selfWeight !== 0) {
+            const { unitWeight, A } = group.section;
+            if (unitWeight === undefined)
+                throw new fea_error_1.FeaError(fea_error_1.FeaErrors.MISSING_UNIT_WEIGHT);
+            // kN/m³ · m² = kN/m
+            const w = -selfWeight * unitWeight * A * 1e-6;
+            if (w !== 0) {
+                memberLoadsHere.push({
+                    loadCase: loadCaseId,
+                    type: "udl",
+                    memberId: group.id,
+                    direction: "Z",
+                    w,
+                });
+            }
+        }
+        const local = memberLoadsHere.flatMap((load) => toLocalLoads(load, group.lambda, group.length));
+        if (group.pieces.length === 1) {
+            memberLoads[group.pieces[0]] = local;
+            return;
+        }
+        group.pieces.forEach((p, k) => {
+            const piece = pieces[p];
+            memberLoads[p] = clipLoads(local, piece.offset, piece.length, k === group.pieces.length - 1);
+        });
+    });
+    return memberLoads;
+};
+exports.memberPieceLoads = memberPieceLoads;
+/** A local member load scaled by `factor`, e.g. for a load combination. */
+const scaleMemberLoad = (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 };
+};
+exports.scaleMemberLoad = scaleMemberLoad;
+//# sourceMappingURL=member-loads.js.map

package/dist/fea/frame/solver/member-loads.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"member-loads.js","sourceRoot":"","sources":["../../../../src/fea/frame/solver/member-loads.ts"],"names":[],"mappings":";;;AAAA,0DAAiE;AAIjE,mEAAmF;AAInF,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,8FAA8F;IAC9F,IAAI,UAAU,GAAG,CAAC,CAAC;IACnB,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO,IAAI,IAAI,CAAC,SAAS,EAAE,CAAC;QAC5C,MAAM,CAAC,GAAG,iBAAiB,CAAC,IAAI,CAAC,SAA4B,CAAC,CAAC;QAC/D,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1E,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACjD,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,GAAG,UAAU,EAAE,cAAc;YAC5C,EAAE,EAAE,EAAE,GAAG,MAAM,GAAG,UAAU;YAC5B,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;;;;GAIG;AACH,MAAM,SAAS,GAAG,CAChB,KAAwB,EACxB,MAAc,EACd,CAAS,EACT,IAAa,EACM,EAAE;IACrB,MAAM,GAAG,GAAG,MAAM,GAAG,CAAC,CAAC;IACvB,MAAM,SAAS,GAAG,IAAI,GAAG,GAAG,CAAC;IAC7B,MAAM,GAAG,GAAsB,EAAE,CAAC;IAClC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;YAC1B,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,IAAI,MAAM,GAAG,SAAS,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,SAAS,IAAI,IAAI,CAAC,CAAC;YAClF,IAAI,MAAM;gBAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC;YAChF,SAAS;QACX,CAAC;QACD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;QAC1C,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACnC,IAAI,EAAE,GAAG,IAAI,IAAI,SAAS;YAAE,SAAS;QACrC,IAAI,IAAI,CAAC,IAAI,KAAK,QAAQ,EAAE,CAAC;YAC3B,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,EAAE,KAAK,EAAE,IAAI,GAAG,MAAM,EAAE,GAAG,EAAE,EAAE,GAAG,MAAM,EAAE,CAAC,CAAC;QAChE,CAAC;aAAM,CAAC;YACN,GAAG,CAAC,IAAI,CAAC;gBACP,GAAG,IAAI;gBACP,EAAE,EAAE,IAAA,gCAAa,EAAC,IAAI,EAAE,IAAI,CAAC;gBAC7B,EAAE,EAAE,IAAA,gCAAa,EAAC,IAAI,EAAE,EAAE,CAAC;gBAC3B,KAAK,EAAE,IAAI,GAAG,MAAM;gBACpB,GAAG,EAAE,EAAE,GAAG,MAAM;aACjB,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAEF;;;;;;;;;;;;;GAaG;AACI,MAAM,gBAAgB,GAAG,CAC9B,MAAqB,EACrB,MAAsB,EACtB,WAA+B,EAC/B,KAAmB,EACnB,UAAkB,EAClB,UAAkB,EACG,EAAE;IACvB,MAAM,aAAa,GAAmB,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC3D,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,aAAa,CAAC,WAAW,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAW,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACrE,CAAC;IAED,yEAAyE;IACzE,MAAM,WAAW,GAAwB,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC9D,MAAM,CAAC,OAAO,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;QAC1B,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,KAAK,CAAC,OAAO,CAAC;YACxC,IAAI,UAAU,KAAK,SAAS;gBAAE,MAAM,IAAI,oBAAQ,CAAC,qBAAS,CAAC,mBAAmB,CAAC,CAAC;YAChF,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,UAAU;oBACpB,IAAI,EAAE,KAAK;oBACX,QAAQ,EAAE,KAAK,CAAC,EAAE;oBAClB,SAAS,EAAE,GAAG;oBACd,CAAC;iBACF,CAAC,CAAC;YACL,CAAC;QACH,CAAC;QACD,MAAM,KAAK,GAAG,eAAe,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,YAAY,CAAC,IAAI,EAAE,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;QAChG,IAAI,KAAK,CAAC,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC9B,WAAW,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;YACrC,OAAO;QACT,CAAC;QACD,KAAK,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YAC5B,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACxB,WAAW,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,KAAK,EAAE,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,EAAE,CAAC,KAAK,KAAK,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAC/F,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,OAAO,WAAW,CAAC;AACrB,CAAC,CAAC;AA5CW,QAAA,gBAAgB,oBA4C3B;AAEF,2EAA2E;AACpE,MAAM,eAAe,GAAG,CAAC,IAAqB,EAAE,MAAc,EAAmB,EAAE;IACxF,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;AAJW,QAAA,eAAe,mBAI1B"}

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

@@ -0,0 +1,50 @@
+import type { SectionProperties } from "../../../core/types/section-properties.type";
+import type { FeaId, FeaNode } from "../../shared/types/model.type";
+import type { Matrix } from "../../../utils/matrix";
+import { type CondensedMember } from "../elements/static-condensation";
+import type { Member } from "../types/model.type";
+import type { AssemblyMember } from "./member-assembly";
+/**
+ * A member, or a piece of one, with everything that does not depend on the loads. A member is
+ * cut into pieces where generated shell-mesh nodes lie inside it; otherwise it is one piece.
+ */
+export interface SystemMember extends AssemblyMember {
+    id: FeaId;
+    /** Distance of the piece's start from the member's start (mm); 0 for a whole member. */
+    offset: number;
+    section: SectionProperties;
+    length: number;
+    lambda: Matrix;
+    /** True when any end DOF is released or tied by a spring. */
+    hasReleases: boolean;
+    /** Condenses a member's fixed-end forces for its end releases. */
+    condenseFixedEndForces: CondensedMember["condenseFixedEndForces"];
+    /** The member's own end displacements from its end nodes' (differ across releases). */
+    memberEndDisplacements: CondensedMember["memberEndDisplacements"];
+}
+/** A model member and the pieces it is analysed as (one unless the shell mesh splits it). */
+export interface MemberGroup {
+    id: FeaId;
+    /** Indices of its pieces among all pieces, from the member's start to its end. */
+    pieces: number[];
+    section: SectionProperties;
+    length: number;
+    lambda: Matrix;
+}
+/**
+ * Builds the pieces of every member of a model that passed `validateModel`: local axes,
+ * transformation and local stiffness condensed for the end releases, which apply only at the
+ * member's own ends. A member is one piece unless mesh nodes lie inside it.
+ *
+ * @param members - Members, in model order; after `expandModel` each lists the generated nodes
+ *   it passes through, in order from its start
+ * @param nodes - The analysed nodes, generated mesh nodes included
+ * @param nodeIndex - Map from node id to index into `nodes`
+ * @returns `pieces` to assemble, and one `group` per member indexing its pieces
+ */
+export declare const prepareSystemMembers: (members: (Member & {
+    through?: FeaId[];
+})[], nodes: FeaNode[], nodeIndex: Map<FeaId, number>) => {
+    pieces: SystemMember[];
+    groups: MemberGroup[];
+};

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

@@ -0,0 +1,70 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.prepareSystemMembers = void 0;
+const beam_3d_stiffness_1 = require("../elements/beam-3d-stiffness");
+const end_releases_1 = require("../elements/end-releases");
+const static_condensation_1 = require("../elements/static-condensation");
+const coordinate_transform_1 = require("../transformation/coordinate-transform");
+/**
+ * Builds the pieces of every member of a model that passed `validateModel`: local axes,
+ * transformation and local stiffness condensed for the end releases, which apply only at the
+ * member's own ends. A member is one piece unless mesh nodes lie inside it.
+ *
+ * @param members - Members, in model order; after `expandModel` each lists the generated nodes
+ *   it passes through, in order from its start
+ * @param nodes - The analysed nodes, generated mesh nodes included
+ * @param nodeIndex - Map from node id to index into `nodes`
+ * @returns `pieces` to assemble, and one `group` per member indexing its pieces
+ */
+const prepareSystemMembers = (members, nodes, nodeIndex) => {
+    const point = (id) => {
+        const n = nodes[nodeIndex.get(id)];
+        return [n.x, n.y, n.z];
+    };
+    const systemMembers = [];
+    const memberGroups = members.map((member) => {
+        // pieces between the member's ends and the mesh nodes inside it (see `expandModel`)
+        const chain = [member.startNodeId, ...(member.through ?? []), member.endNodeId];
+        const whole = (0, coordinate_transform_1.memberAxes)(point(member.startNodeId), point(member.endNodeId), member.angle ?? 0);
+        const { E, G, A, I22, I33, J } = member.section;
+        const pieces = [];
+        let offset = 0;
+        for (let k = 0; k + 1 < chain.length; k++) {
+            const { length, lambda } = chain.length === 2
+                ? whole
+                : (0, coordinate_transform_1.memberAxes)(point(chain[k]), point(chain[k + 1]), member.angle ?? 0);
+            // releases only at the member's own ends
+            const releases = (0, end_releases_1.memberEndReleases)({
+                start: k === 0 ? member.releases?.start : undefined,
+                end: k === chain.length - 2 ? member.releases?.end : undefined,
+            });
+            const condensed = (0, static_condensation_1.condenseMember)((0, beam_3d_stiffness_1.beam3dLocalStiffnessRightHanded)(E, G, A, I22, I33, J, length), releases);
+            pieces.push(systemMembers.length);
+            systemMembers.push({
+                id: member.id,
+                offset,
+                section: member.section,
+                startIndex: nodeIndex.get(chain[k]),
+                endIndex: nodeIndex.get(chain[k + 1]),
+                T: (0, coordinate_transform_1.transformationMatrix)(lambda),
+                k: condensed.k,
+                hasReleases: releases.length > 0,
+                condenseFixedEndForces: condensed.condenseFixedEndForces,
+                memberEndDisplacements: condensed.memberEndDisplacements,
+                length,
+                lambda,
+            });
+            offset += length;
+        }
+        return {
+            id: member.id,
+            pieces,
+            section: member.section,
+            length: whole.length,
+            lambda: whole.lambda,
+        };
+    });
+    return { pieces: systemMembers, groups: memberGroups };
+};
+exports.prepareSystemMembers = prepareSystemMembers;
+//# sourceMappingURL=system-member.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"system-member.js","sourceRoot":"","sources":["../../../../src/fea/frame/solver/system-member.ts"],"names":[],"mappings":";;;AAIA,qEAAgF;AAChF,2DAA6D;AAC7D,yEAAuF;AACvF,iFAA0F;AAiC1F;;;;;;;;;;GAUG;AACI,MAAM,oBAAoB,GAAG,CAClC,OAA2C,EAC3C,KAAgB,EAChB,SAA6B,EACsB,EAAE;IACrD,MAAM,KAAK,GAAG,CAAC,EAAS,EAA4B,EAAE;QACpD,MAAM,CAAC,GAAG,KAAK,CAAC,SAAS,CAAC,GAAG,CAAC,EAAE,CAAW,CAAC,CAAC;QAC7C,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC;IACF,MAAM,aAAa,GAAmB,EAAE,CAAC;IACzC,MAAM,YAAY,GAAkB,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QACzD,oFAAoF;QACpF,MAAM,KAAK,GAAG,CAAC,MAAM,CAAC,WAAW,EAAE,GAAG,CAAC,MAAM,CAAC,OAAO,IAAI,EAAE,CAAC,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC;QAChF,MAAM,KAAK,GAAG,IAAA,iCAAU,EAAC,KAAK,CAAC,MAAM,CAAC,WAAW,CAAC,EAAE,KAAK,CAAC,MAAM,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC;QAChG,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC;QAChD,MAAM,MAAM,GAAa,EAAE,CAAC;QAC5B,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GACtB,KAAK,CAAC,MAAM,KAAK,CAAC;gBAChB,CAAC,CAAC,KAAK;gBACP,CAAC,CAAC,IAAA,iCAAU,EAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC;YAC1E,yCAAyC;YACzC,MAAM,QAAQ,GAAG,IAAA,gCAAiB,EAAC;gBACjC,KAAK,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC,SAAS;gBACnD,GAAG,EAAE,CAAC,KAAK,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC,CAAC,SAAS;aAC/D,CAAC,CAAC;YACH,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,QAAQ,CACT,CAAC;YACF,MAAM,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,CAAC;YAClC,aAAa,CAAC,IAAI,CAAC;gBACjB,EAAE,EAAE,MAAM,CAAC,EAAE;gBACb,MAAM;gBACN,OAAO,EAAE,MAAM,CAAC,OAAO;gBACvB,UAAU,EAAE,SAAS,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAW;gBAC7C,QAAQ,EAAE,SAAS,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAW;gBAC/C,CAAC,EAAE,IAAA,2CAAoB,EAAC,MAAM,CAAC;gBAC/B,CAAC,EAAE,SAAS,CAAC,CAAC;gBACd,WAAW,EAAE,QAAQ,CAAC,MAAM,GAAG,CAAC;gBAChC,sBAAsB,EAAE,SAAS,CAAC,sBAAsB;gBACxD,sBAAsB,EAAE,SAAS,CAAC,sBAAsB;gBACxD,MAAM;gBACN,MAAM;aACP,CAAC,CAAC;YACH,MAAM,IAAI,MAAM,CAAC;QACnB,CAAC;QACD,OAAO;YACL,EAAE,EAAE,MAAM,CAAC,EAAE;YACb,MAAM;YACN,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,MAAM,EAAE,KAAK,CAAC,MAAM;YACpB,MAAM,EAAE,KAAK,CAAC,MAAM;SACrB,CAAC;IACJ,CAAC,CAAC,CAAC;IACH,OAAO,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,YAAY,EAAE,CAAC;AACzD,CAAC,CAAC;AAzDW,QAAA,oBAAoB,wBAyD/B"}

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

@@ -11,7 +11,7 @@ import { type Matrix } from "../../../utils/matrix";
  * @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
+ * @throws {FeaError} 301 if the nodes coincide
  */
 export declare const memberAxes: (start: [number, number, number], end: [number, number, number], angle?: number) => {
     length: number;

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

@@ -2,7 +2,7 @@
 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 fea_error_1 = require("../../../core/errors/fea.error");
 const matrix_1 = require("../../../utils/matrix");
 const cross = (a, b) => [
     a[1] * b[2] - a[2] * b[1],
@@ -23,13 +23,13 @@ const scale = (a, s) => [a[0] * s, a[1] * s, a[2] * s];
  * @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
+ * @throws {FeaError} 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);
+        throw new fea_error_1.FeaError(fea_error_1.FeaErrors.ZERO_LENGTH_MEMBER);
     const e1 = scale(d, 1 / length);
     let e2;
     let e3;

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

@@ -1 +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"}
+{"version":3,"file":"coordinate-transform.js","sourceRoot":"","sources":["../../../../src/fea/frame/transformation/coordinate-transform.ts"],"names":[],"mappings":";;;AAAA,uEAA+E;AAC/E,0DAAiE;AACjE,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,oBAAQ,CAAC,qBAAS,CAAC,kBAAkB,CAAC,CAAC;IAEpE,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

@@ -1,37 +1,5 @@
 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;
-/**
- * Fixed DOFs along a straight line: every node on the segment between two nodes, including the
- * nodes the engine generates when it meshes shell areas, e.g. a wall base fixed along its
- * length. `true` fixes the DOF; there are no line springs. Combined with point restraints at
- * the same node, a fixed DOF stays fixed.
- */
-export interface LineRestraint {
-    nodeIds: [FrameId, FrameId];
-    ux?: boolean;
-    uy?: boolean;
-    uz?: boolean;
-    rx?: boolean;
-    ry?: boolean;
-    rz?: boolean;
-}
-export interface Restraint {
-    nodeId: FrameId;
-    ux?: DofRestraint;
-    uy?: DofRestraint;
-    uz?: DofRestraint;
-    rx?: DofRestraint;
-    ry?: DofRestraint;
-    rz?: DofRestraint;
-}
+import type { FeaId } from "../../shared/types/model.type";
 /**
  * How a member end DOF is tied to its node: `true` (or 0) is a full release that carries no
  * force, a positive number is a spring stiffness in series with the member end (partial
@@ -48,9 +16,9 @@ export interface MemberEndRelease {
     r3?: DofRelease;
 }
 export interface Member {
-    id: FrameId;
-    startNodeId: FrameId;
-    endNodeId: FrameId;
+    id: FeaId;
+    startNodeId: FeaId;
+    endNodeId: FeaId;
     section: SectionProperties;
     angle?: number;
     releases?: {
@@ -58,16 +26,6 @@ export interface Member {
         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";
 /**
@@ -79,7 +37,7 @@ export type LoadDirection = "1" | "2" | "3" | "X" | "Y" | "Z";
 export interface UniformMemberLoad {
     loadCase: string;
     type: "udl";
-    memberId: FrameId;
+    memberId: FeaId;
     direction: LoadDirection;
     w: number;
     start?: number;
@@ -89,7 +47,7 @@ export interface UniformMemberLoad {
 export interface PointMemberLoad {
     loadCase: string;
     type: "point";
-    memberId: FrameId;
+    memberId: FeaId;
     direction: LoadDirection;
     P: number;
     a: number;
@@ -102,7 +60,7 @@ export interface PointMemberLoad {
 export interface TrapezoidalMemberLoad {
     loadCase: string;
     type: "trapezoidal";
-    memberId: FrameId;
+    memberId: FeaId;
     direction: LoadDirection;
     w1: number;
     w2: number;
@@ -117,7 +75,7 @@ export interface TrapezoidalMemberLoad {
 export interface TorqueMemberLoad {
     loadCase: string;
     type: "torque";
-    memberId: FrameId;
+    memberId: FeaId;
     t: number;
     start?: number;
     end?: number;
@@ -137,132 +95,9 @@ export type AreaDistribution = "twoWay" | "oneWay" | "cantilever";
  * nodes, so nodes listed counter-clockwise seen from above give local 3 = +Z.
  */
 export interface NullArea {
-    id: FrameId;
-    nodeIds: FrameId[];
+    id: FeaId;
+    nodeIds: FeaId[];
     distribution: AreaDistribution;
     /** Edge indices that carry the load: one-way supports, or the one cantilever support. */
     supportEdges?: number[];
 }
-/**
- * How a shell area carries load: in its plane and in bending (`"shell"`), in its plane only
- * (`"membrane"`: u1, u2 and the rotation about local 3), or in bending only (`"plate"`: u3 and
- * the rotations about local 1 and 2).
- */
-export type ShellType = "shell" | "membrane" | "plate";
-/** Material and thickness of a shell area. */
-export interface ShellSection {
-    thickness: number;
-    E: number;
-    nu: number;
-    unitWeight?: number;
-}
-/**
- * Stiffness modifiers of a shell area along its local axes 1 and 2 (default 1 each), e.g. for
- * cracked walls and slabs: f11, f22, f12 scale the membrane rigidity and m11, m22, m12 the
- * bending rigidity (each matrix becomes S·D·S, S = diag(√f11, √f22, √f12)); weight scales the
- * self-weight only.
- */
-export interface ShellModifiers {
-    f11?: number;
-    f22?: number;
-    f12?: number;
-    m11?: number;
-    m22?: number;
-    m12?: number;
-    weight?: number;
-}
-/**
- * How the engine meshes a shell area: `maxSize` (mm) is the largest element side wanted;
- * `divisions` the least number of elements along edges 0–1 and 1–2 (from `nodeIds[0]` to
- * `nodeIds[1]` and from `nodeIds[1]` to `nodeIds[2]`). Counts may grow so that areas sharing an
- * edge, or meeting it at a node, share their mesh nodes.
- */
-export interface ShellMesh {
-    maxSize?: number;
-    divisions?: [number, number];
-}
-/**
- * A four-node flat shell area: a drilling-rotation membrane and a thin (DKQ) plate, built on
- * the mean plane of its corners. Without `mesh` it is one element, unless a neighbour's mesh or
- * a node lying on its edge forces divisions. Local axis 3 is the normal by the
- * right-hand rule around `nodeIds`; local 2 is +Z projected onto the plane (upward on a wall)
- * and local 1 = 2 × 3 is horizontal, except on a horizontal shell, where local 1 = +X.
- * `angle` rotates local 1 and 2 about local 3. Its id shares one namespace with null areas, so
- * `AreaLoad` can load it.
- */
-export interface ShellArea {
-    id: FrameId;
-    nodeIds: FrameId[];
-    type?: ShellType;
-    section: ShellSection;
-    angle?: number;
-    modifiers?: ShellModifiers;
-    mesh?: ShellMesh;
-}
-/** Area load direction: the area's local normal "3" or a global axis. */
-export type AreaLoadDirection = "3" | "X" | "Y" | "Z";
-/**
- * Uniform load on a null area or a shell area. `q` is per unit true area, or with
- * `projected: true` (global directions only) per unit area projected onto the plane normal to
- * `direction`, e.g. per plan area on a ramp; it is then applied as `q · |n · direction|` on
- * the true area. A null area sends it to the members along its edges; a shell area takes it
- * as forces at its corners, each corner the share `∫ Nᵢ dA` of its area.
- */
-export interface AreaLoad {
-    loadCase: string;
-    areaId: FrameId;
-    direction: AreaLoadDirection;
-    q: number;
-    projected?: boolean;
-}
-/** 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 and `selfWeight · unitWeight · thickness` (times `modifiers.weight`)
-     * per unit area of every shell, 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[];
-    /** Fixed DOFs of every node along straight lines, generated mesh nodes included. */
-    lineRestraints?: LineRestraint[];
-    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[];
-    /** Shell areas: four-node flat shell elements with membrane and/or plate stiffness. */
-    shells?: ShellArea[];
-    /** Area loads of every load case, each tagged with its `loadCase`. */
-    areaLoads?: AreaLoad[];
-    combinations?: LoadCombination[];
-}