@theerapat-s28/ts-structural-eng-tools 2.2.0 → 3.0.0

package/dist/fea/frame/elements/shell/quad-geometry.js

@@ -0,0 +1,136 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.cornerAreas = exports.cartesianDerivatives = exports.quadJacobian = exports.serendipityShape = exports.midsideShape = exports.bilinearShape = exports.GAUSS_3X3 = exports.GAUSS_2X2 = void 0;
+/** Natural coordinates (ξ, η) of the four corners, counter-clockwise about the normal. */
+const CORNER_XI = [-1, 1, 1, -1];
+const CORNER_ETA = [-1, -1, 1, 1];
+const GAUSS = 1 / Math.sqrt(3);
+/**
+ * 2×2 Gauss rule over the parent square −1 ≤ ξ, η ≤ 1: points (±1/√3, ±1/√3), weight 1 each.
+ * Exact for polynomials up to cubic in each of ξ and η.
+ */
+exports.GAUSS_2X2 = [
+    { xi: -GAUSS, eta: -GAUSS, weight: 1 },
+    { xi: GAUSS, eta: -GAUSS, weight: 1 },
+    { xi: GAUSS, eta: GAUSS, weight: 1 },
+    { xi: -GAUSS, eta: GAUSS, weight: 1 },
+];
+const GAUSS_3 = Math.sqrt(0.6);
+const GAUSS_3_POINTS = [
+    { at: -GAUSS_3, weight: 5 / 9 },
+    { at: 0, weight: 8 / 9 },
+    { at: GAUSS_3, weight: 5 / 9 },
+];
+/**
+ * 3×3 Gauss rule over the parent square: points 0 and ±√0.6 in each direction, weights the
+ * products of 5/9, 8/9, 5/9. Exact for polynomials up to quintic in each of ξ and η.
+ */
+exports.GAUSS_3X3 = GAUSS_3_POINTS.flatMap((p) => GAUSS_3_POINTS.map((q) => ({ xi: p.at, eta: q.at, weight: p.weight * q.weight })));
+/** Bilinear corner functions `Nᵢ = ¼(1 + ξᵢξ)(1 + ηᵢη)`, i = 1…4. */
+const bilinearShape = (xi, eta) => {
+    const N = [];
+    const dXi = [];
+    const dEta = [];
+    for (let i = 0; i < 4; i++) {
+        const a = CORNER_XI[i];
+        const b = CORNER_ETA[i];
+        N.push(0.25 * (1 + a * xi) * (1 + b * eta));
+        dXi.push(0.25 * a * (1 + b * eta));
+        dEta.push(0.25 * b * (1 + a * xi));
+    }
+    return { N, dXi, dEta };
+};
+exports.bilinearShape = bilinearShape;
+/**
+ * Mid-side functions N₅…N₈, one per edge 1→2, 2→3, 3→4, 4→1 (index k = edge from corner k to
+ * corner k + 1): 1 at the edge's midpoint and 0 at every corner and every other midpoint.
+ */
+const midsideShape = (xi, eta) => ({
+    N: [
+        0.5 * (1 - xi * xi) * (1 - eta),
+        0.5 * (1 + xi) * (1 - eta * eta),
+        0.5 * (1 - xi * xi) * (1 + eta),
+        0.5 * (1 - xi) * (1 - eta * eta),
+    ],
+    dXi: [-xi * (1 - eta), 0.5 * (1 - eta * eta), -xi * (1 + eta), -0.5 * (1 - eta * eta)],
+    dEta: [-0.5 * (1 - xi * xi), -(1 + xi) * eta, 0.5 * (1 - xi * xi), -(1 - xi) * eta],
+});
+exports.midsideShape = midsideShape;
+/**
+ * 8-node serendipity functions: the corner functions
+ * `Nᵢ = ¼(1 + ξᵢξ)(1 + ηᵢη)(ξᵢξ + ηᵢη − 1)`, i = 1…4, followed by the mid-side functions
+ * N₅…N₈ of `midsideShape`. Each is 1 at its own node and 0 at the other seven.
+ */
+const serendipityShape = (xi, eta) => {
+    const N = [];
+    const dXi = [];
+    const dEta = [];
+    for (let i = 0; i < 4; i++) {
+        const a = CORNER_XI[i];
+        const b = CORNER_ETA[i];
+        N.push(0.25 * (1 + a * xi) * (1 + b * eta) * (a * xi + b * eta - 1));
+        dXi.push(0.25 * a * (1 + b * eta) * (2 * a * xi + b * eta));
+        dEta.push(0.25 * b * (1 + a * xi) * (a * xi + 2 * b * eta));
+    }
+    const mid = (0, exports.midsideShape)(xi, eta);
+    return { N: [...N, ...mid.N], dXi: [...dXi, ...mid.dXi], dEta: [...dEta, ...mid.dEta] };
+};
+exports.serendipityShape = serendipityShape;
+/**
+ * Jacobian of the bilinear map from the parent square to a quadrilateral.
+ *
+ * @param corners - The four corners in the element plane (mm), counter-clockwise
+ * @param xi - Natural coordinate ξ
+ * @param eta - Natural coordinate η
+ * @returns `det J` and `J⁻¹`; `det J > 0` everywhere for a convex, counter-clockwise quad
+ */
+const quadJacobian = (corners, xi, eta) => {
+    const { dXi, dEta } = (0, exports.bilinearShape)(xi, eta);
+    let xXi = 0;
+    let yXi = 0;
+    let xEta = 0;
+    let yEta = 0;
+    for (let i = 0; i < 4; i++) {
+        xXi += dXi[i] * corners[i][0];
+        yXi += dXi[i] * corners[i][1];
+        xEta += dEta[i] * corners[i][0];
+        yEta += dEta[i] * corners[i][1];
+    }
+    const detJ = xXi * yEta - yXi * xEta;
+    return {
+        detJ,
+        inverse: [
+            [yEta / detJ, -yXi / detJ],
+            [-xEta / detJ, xXi / detJ],
+        ],
+    };
+};
+exports.quadJacobian = quadJacobian;
+/**
+ * Converts ξ, η derivatives of shape functions to x, y derivatives: `[N,x; N,y] = J⁻¹·[N,ξ; N,η]`.
+ *
+ * @returns `dX[i]` and `dY[i]` for every function i of `shape`
+ */
+const cartesianDerivatives = ({ inverse }, shape) => ({
+    dX: shape.dXi.map((d, i) => inverse[0][0] * d + inverse[0][1] * shape.dEta[i]),
+    dY: shape.dXi.map((d, i) => inverse[1][0] * d + inverse[1][1] * shape.dEta[i]),
+});
+exports.cartesianDerivatives = cartesianDerivatives;
+/**
+ * Integrals `∫ Nᵢ dA` of the bilinear corner functions over a quad (2×2 Gauss, exact): the
+ * share of a uniform surface load each corner takes. They sum to the quad's area; each is
+ * A/4 for a parallelogram.
+ *
+ * @param corners - The four corners in the element plane (mm), counter-clockwise
+ * @returns Four areas in mm²
+ */
+const cornerAreas = (corners) => {
+    const areas = [0, 0, 0, 0];
+    for (const { xi, eta, weight } of exports.GAUSS_2X2) {
+        const { detJ } = (0, exports.quadJacobian)(corners, xi, eta);
+        (0, exports.bilinearShape)(xi, eta).N.forEach((n, i) => (areas[i] += n * detJ * weight));
+    }
+    return areas;
+};
+exports.cornerAreas = cornerAreas;
+//# sourceMappingURL=quad-geometry.js.map

package/dist/fea/frame/elements/shell/quad-geometry.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"quad-geometry.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/shell/quad-geometry.ts"],"names":[],"mappings":";;;AAGA,0FAA0F;AAC1F,MAAM,SAAS,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAU,CAAC;AAC1C,MAAM,UAAU,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAU,CAAC;AAE3C,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;AAE/B;;;GAGG;AACU,QAAA,SAAS,GAA2D;IAC/E,EAAE,EAAE,EAAE,CAAC,KAAK,EAAE,GAAG,EAAE,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE;IACtC,EAAE,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE;IACrC,EAAE,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE;IACpC,EAAE,EAAE,EAAE,CAAC,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE;CACtC,CAAC;AAEF,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;AAC/B,MAAM,cAAc,GAAG;IACrB,EAAE,EAAE,EAAE,CAAC,OAAO,EAAE,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE;IAC/B,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE;IACxB,EAAE,EAAE,EAAE,OAAO,EAAE,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE;CAC/B,CAAC;AAEF;;;GAGG;AACU,QAAA,SAAS,GACpB,cAAc,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAC3B,cAAc,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,CAClF,CAAC;AASJ,qEAAqE;AAC9D,MAAM,aAAa,GAAG,CAAC,EAAU,EAAE,GAAW,EAAe,EAAE;IACpE,MAAM,CAAC,GAAa,EAAE,CAAC;IACvB,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACxB,CAAC,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC5C,GAAG,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACrC,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,EAAE,CAAC;AAC1B,CAAC,CAAC;AAZW,QAAA,aAAa,iBAYxB;AAEF;;;GAGG;AACI,MAAM,YAAY,GAAG,CAAC,EAAU,EAAE,GAAW,EAAe,EAAE,CAAC,CAAC;IACrE,CAAC,EAAE;QACD,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QAC/B,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;QAChC,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QAC/B,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;KACjC;IACD,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;IACtF,IAAI,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC;CACpF,CAAC,CAAC;AATU,QAAA,YAAY,gBAStB;AAEH;;;;GAIG;AACI,MAAM,gBAAgB,GAAG,CAAC,EAAU,EAAE,GAAW,EAAe,EAAE;IACvE,MAAM,CAAC,GAAa,EAAE,CAAC;IACvB,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACxB,CAAC,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACrE,GAAG,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC5D,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD,MAAM,GAAG,GAAG,IAAA,oBAAY,EAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IAClC,OAAO,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,IAAI,EAAE,CAAC,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;AAC1F,CAAC,CAAC;AAbW,QAAA,gBAAgB,oBAa3B;AAQF;;;;;;;GAOG;AACI,MAAM,YAAY,GAAG,CAAC,OAAwB,EAAE,EAAU,EAAE,GAAW,EAAgB,EAAE;IAC9F,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,IAAA,qBAAa,EAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IAC7C,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9B,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9B,IAAI,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC;IACD,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC;IACrC,OAAO;QACL,IAAI;QACJ,OAAO,EAAE;YACP,CAAC,IAAI,GAAG,IAAI,EAAE,CAAC,GAAG,GAAG,IAAI,CAAC;YAC1B,CAAC,CAAC,IAAI,GAAG,IAAI,EAAE,GAAG,GAAG,IAAI,CAAC;SAC3B;KACF,CAAC;AACJ,CAAC,CAAC;AApBW,QAAA,YAAY,gBAoBvB;AAEF;;;;GAIG;AACI,MAAM,oBAAoB,GAAG,CAClC,EAAE,OAAO,EAAgB,EACzB,KAAwC,EACR,EAAE,CAAC,CAAC;IACpC,EAAE,EAAE,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9E,EAAE,EAAE,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;CAC/E,CAAC,CAAC;AANU,QAAA,oBAAoB,wBAM9B;AAEH;;;;;;;GAOG;AACI,MAAM,WAAW,GAAG,CAAC,OAAwB,EAAY,EAAE;IAChE,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,KAAK,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,iBAAS,EAAE,CAAC;QAC5C,MAAM,EAAE,IAAI,EAAE,GAAG,IAAA,oBAAY,EAAC,OAAO,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QAChD,IAAA,qBAAa,EAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC;IAC9E,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC,CAAC;AAPW,QAAA,WAAW,eAOtB"}

package/dist/fea/frame/elements/shell/rigidity.d.ts

@@ -0,0 +1,14 @@
+import type { Matrix } from "../../../../utils/matrix";
+/** Stiffness modifiers on the components 11, 22 and 12 of a rigidity matrix (1 = unmodified). */
+export type RigidityModifiers = readonly [number, number, number];
+/**
+ * Isotropic rigidity `c·[[1, ν, 0], [ν, 1, 0], [0, 0, (1 − ν)/2]]` with stiffness modifiers applied
+ * as `S·D·S`, `S = diag(√f₁₁, √f₂₂, √f₁₂)`: D₁₁, D₂₂ and D₃₃ scale by their own modifier and the
+ * coupling D₁₂ by √(f₁₁·f₂₂), so D stays symmetric and positive definite and equal modifiers
+ * scale it uniformly.
+ *
+ * @param c - E·t/(1 − ν²) for a membrane, E·t³/(12(1 − ν²)) for a plate
+ * @param nu - Poisson's ratio
+ * @param modifiers - Positive modifiers [f₁₁, f₂₂, f₁₂] (default all 1)
+ */
+export declare const isotropicRigidity: (c: number, nu: number, modifiers?: RigidityModifiers) => Matrix;

package/dist/fea/frame/elements/shell/rigidity.js

@@ -0,0 +1,23 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.isotropicRigidity = void 0;
+/**
+ * Isotropic rigidity `c·[[1, ν, 0], [ν, 1, 0], [0, 0, (1 − ν)/2]]` with stiffness modifiers applied
+ * as `S·D·S`, `S = diag(√f₁₁, √f₂₂, √f₁₂)`: D₁₁, D₂₂ and D₃₃ scale by their own modifier and the
+ * coupling D₁₂ by √(f₁₁·f₂₂), so D stays symmetric and positive definite and equal modifiers
+ * scale it uniformly.
+ *
+ * @param c - E·t/(1 − ν²) for a membrane, E·t³/(12(1 − ν²)) for a plate
+ * @param nu - Poisson's ratio
+ * @param modifiers - Positive modifiers [f₁₁, f₂₂, f₁₂] (default all 1)
+ */
+const isotropicRigidity = (c, nu, modifiers = [1, 1, 1]) => {
+    const [f11, f22, f12] = modifiers;
+    return [
+        [c * f11, c * nu * Math.sqrt(f11 * f22), 0],
+        [c * nu * Math.sqrt(f11 * f22), c * f22, 0],
+        [0, 0, ((c * (1 - nu)) / 2) * f12],
+    ];
+};
+exports.isotropicRigidity = isotropicRigidity;
+//# sourceMappingURL=rigidity.js.map

package/dist/fea/frame/elements/shell/rigidity.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"rigidity.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/shell/rigidity.ts"],"names":[],"mappings":";;;AAKA;;;;;;;;;GASG;AACI,MAAM,iBAAiB,GAAG,CAC/B,CAAS,EACT,EAAU,EACV,YAA+B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAChC,EAAE;IACV,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,SAAS,CAAC;IAClC,OAAO;QACL,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;QAC3C,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QAC3C,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;KACnC,CAAC;AACJ,CAAC,CAAC;AAXW,QAAA,iBAAiB,qBAW5B"}

package/dist/fea/frame/elements/shell/shell-axes.d.ts

@@ -0,0 +1,43 @@
+import type { Matrix } from "../../../../utils/matrix";
+import { type Vec2 } from "./quad-geometry";
+/** A point or vector in global axes: X, Y, Z in mm. */
+export type Vec3 = [number, number, number];
+/** A four-node shell element's local axes and its corners on the mean plane. */
+export interface ShellAxes {
+    /** Direction cosines, rows = local 1, 2, 3 in global components. */
+    lambda: Matrix;
+    /** Centroid of the four corners (mm), the origin of the local coordinates. */
+    origin: Vec3;
+    /** Corners projected onto the mean plane, as local 1–2 coordinates (mm). */
+    corners: Vec2[];
+    /**
+     * Distance of each corner from the mean plane along local 3 (mm): +h, −h, +h, −h, with
+     * h = ¼(x₁ − x₂ + x₃ − x₄)·n. All zero for a flat element.
+     */
+    offsets: number[];
+}
+/**
+ * Local axes of a four-node shell element, ETABS-style, and its corners on the mean plane.
+ *
+ * - Local 3 is the unit normal `n ∥ d₁ × d₂` of the diagonals `d₁ = x₃ − x₁`, `d₂ = x₄ − x₂`: the
+ *   right-hand rule around the nodes, as for null areas (for four nodes Newell's vector area is
+ *   exactly `½ d₁ × d₂`).
+ * - Non-horizontal element: local 2 is +Z projected onto the plane (upward on a wall) and
+ *   local 1 = 2 × 3 is horizontal.
+ * - Horizontal element (the normal's horizontal part below the vertical-member tolerance):
+ *   local 1 = +X and local 2 = 3 × 1, so local 2 = +Y when local 3 = +Z.
+ * - `angle` then rotates local 1 and 2 about local 3 by the right-hand rule.
+ *
+ * The mean plane passes through the centroid with normal n; a warped element's corners lie
+ * alternately h above and below it.
+ *
+ * @param points - The four corners in global axes (mm), in order around the element
+ * @param angle - Rotation of local 1 and 2 about local 3 in rad (default 0)
+ */
+export declare const shellAxes: (points: readonly Vec3[], angle?: number) => ShellAxes;
+/**
+ * Whether four corners make a valid shell element: no repeated corner, a non-degenerate
+ * convex quadrilateral in order around it (det J > 0 at every corner of the bilinear map), and
+ * no corner farther than `SHELL_WARPING_LIMIT · √A` from the mean plane.
+ */
+export declare const isValidShellQuad: (points: readonly Vec3[]) => boolean;

package/dist/fea/frame/elements/shell/shell-axes.js

@@ -0,0 +1,102 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.isValidShellQuad = exports.shellAxes = void 0;
+const frame_constant_1 = require("../../../../core/constants/frame.constant");
+const quad_geometry_1 = require("./quad-geometry");
+const sub = (a, b) => [a[0] - b[0], a[1] - b[1], a[2] - b[2]];
+const dot = (a, b) => a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
+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 unit = (a) => {
+    const n = Math.hypot(a[0], a[1], a[2]);
+    return [a[0] / n, a[1] / n, a[2] / n];
+};
+/**
+ * Local axes of a four-node shell element, ETABS-style, and its corners on the mean plane.
+ *
+ * - Local 3 is the unit normal `n ∥ d₁ × d₂` of the diagonals `d₁ = x₃ − x₁`, `d₂ = x₄ − x₂`: the
+ *   right-hand rule around the nodes, as for null areas (for four nodes Newell's vector area is
+ *   exactly `½ d₁ × d₂`).
+ * - Non-horizontal element: local 2 is +Z projected onto the plane (upward on a wall) and
+ *   local 1 = 2 × 3 is horizontal.
+ * - Horizontal element (the normal's horizontal part below the vertical-member tolerance):
+ *   local 1 = +X and local 2 = 3 × 1, so local 2 = +Y when local 3 = +Z.
+ * - `angle` then rotates local 1 and 2 about local 3 by the right-hand rule.
+ *
+ * The mean plane passes through the centroid with normal n; a warped element's corners lie
+ * alternately h above and below it.
+ *
+ * @param points - The four corners in global axes (mm), in order around the element
+ * @param angle - Rotation of local 1 and 2 about local 3 in rad (default 0)
+ */
+const shellAxes = (points, angle = 0) => {
+    const e3 = unit(cross(sub(points[2], points[0]), sub(points[3], points[1])));
+    let e1;
+    let e2;
+    if (Math.hypot(e3[0], e3[1]) < frame_constant_1.VERTICAL_MEMBER_TOLERANCE) {
+        e1 = [1, 0, 0];
+        e2 = cross(e3, e1);
+    }
+    else {
+        e2 = unit([-e3[2] * e3[0], -e3[2] * e3[1], 1 - e3[2] * e3[2]]); // Z − (Z·n)·n
+        e1 = cross(e2, e3);
+    }
+    if (angle !== 0) {
+        const c = Math.cos(angle);
+        const s = Math.sin(angle);
+        const r1 = [c * e1[0] + s * e2[0], c * e1[1] + s * e2[1], c * e1[2] + s * e2[2]];
+        const r2 = [-s * e1[0] + c * e2[0], -s * e1[1] + c * e2[1], -s * e1[2] + c * e2[2]];
+        e1 = r1;
+        e2 = r2;
+    }
+    const origin = [0, 1, 2].map((k) => (points[0][k] + points[1][k] + points[2][k] + points[3][k]) / 4);
+    const relative = points.map((p) => sub(p, origin));
+    return {
+        lambda: [e1, e2, e3],
+        origin,
+        corners: relative.map((r) => [dot(r, e1), dot(r, e2)]),
+        offsets: relative.map((r) => dot(r, e3)),
+    };
+};
+exports.shellAxes = shellAxes;
+/**
+ * Whether four corners make a valid shell element: no repeated corner, a non-degenerate
+ * convex quadrilateral in order around it (det J > 0 at every corner of the bilinear map), and
+ * no corner farther than `SHELL_WARPING_LIMIT · √A` from the mean plane.
+ */
+const isValidShellQuad = (points) => {
+    const size = Math.max(...points.map((p, i) => {
+        const q = points[(i + 1) % 4];
+        return Math.hypot(q[0] - p[0], q[1] - p[1], q[2] - p[2]);
+    }));
+    const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * size;
+    for (let i = 0; i < 4; i++) {
+        for (let j = i + 1; j < 4; j++) {
+            const [p, q] = [points[i], points[j]];
+            if (!(Math.hypot(q[0] - p[0], q[1] - p[1], q[2] - p[2]) > tolerance))
+                return false;
+        }
+    }
+    const d1 = [0, 1, 2].map((k) => points[2][k] - points[0][k]);
+    const d2 = [0, 1, 2].map((k) => points[3][k] - points[1][k]);
+    const area = Math.hypot(d1[1] * d2[2] - d1[2] * d2[1], d1[2] * d2[0] - d1[0] * d2[2], d1[0] * d2[1] - d1[1] * d2[0]) / 2;
+    if (!(area > tolerance * size))
+        return false;
+    const axes = (0, exports.shellAxes)(points);
+    for (const [xi, eta] of [
+        [-1, -1],
+        [1, -1],
+        [1, 1],
+        [-1, 1],
+    ]) {
+        if (!((0, quad_geometry_1.quadJacobian)(axes.corners, xi, eta).detJ > (frame_constant_1.AREA_GEOMETRY_TOLERANCE * area) / 4)) {
+            return false;
+        }
+    }
+    return Math.max(...axes.offsets.map(Math.abs)) <= frame_constant_1.SHELL_WARPING_LIMIT * Math.sqrt(area);
+};
+exports.isValidShellQuad = isValidShellQuad;
+//# sourceMappingURL=shell-axes.js.map

package/dist/fea/frame/elements/shell/shell-axes.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"shell-axes.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/shell/shell-axes.ts"],"names":[],"mappings":";;;AAAA,uEAI4C;AAG5C,mDAA0D;AAK1D,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAChF,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1E,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,EAAQ,EAAE;IAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACxC,CAAC,CAAC;AAiBF;;;;;;;;;;;;;;;;;GAiBG;AACI,MAAM,SAAS,GAAG,CAAC,MAAuB,EAAE,KAAK,GAAG,CAAC,EAAa,EAAE;IACzE,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7E,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,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,cAAc;QAC9E,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACrB,CAAC;IACD,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC1B,MAAM,EAAE,GAAS,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACvF,MAAM,EAAE,GAAS,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1F,EAAE,GAAG,EAAE,CAAC;QACR,EAAE,GAAG,EAAE,CAAC;IACV,CAAC;IAED,MAAM,MAAM,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAChC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAC/D,CAAC;IACV,MAAM,QAAQ,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;IACnD,OAAO;QACL,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;QACpB,MAAM;QACN,OAAO,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAQ,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC5D,OAAO,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;KACzC,CAAC;AACJ,CAAC,CAAC;AA9BW,QAAA,SAAS,aA8BpB;AAEF;;;;GAIG;AACI,MAAM,gBAAgB,GAAG,CAAC,MAAuB,EAAW,EAAE;IACnE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CACnB,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACrB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3D,CAAC,CAAC,CACH,CAAC;IACF,MAAM,SAAS,GAAG,wCAAuB,GAAG,IAAI,CAAC;IACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;gBAAE,OAAO,KAAK,CAAC;QACrF,CAAC;IACH,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,MAAM,IAAI,GACR,IAAI,CAAC,KAAK,CACR,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAC9B,GAAG,CAAC,CAAC;IACR,IAAI,CAAC,CAAC,IAAI,GAAG,SAAS,GAAG,IAAI,CAAC;QAAE,OAAO,KAAK,CAAC;IAC7C,MAAM,IAAI,GAAG,IAAA,iBAAS,EAAC,MAAM,CAAC,CAAC;IAC/B,KAAK,MAAM,CAAC,EAAE,EAAE,GAAG,CAAC,IAAI;QACtB,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACR,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACP,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;KACR,EAAE,CAAC;QACF,IAAI,CAAC,CAAC,IAAA,4BAAY,EAAC,IAAI,CAAC,OAAO,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,wCAAuB,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACvF,OAAO,KAAK,CAAC;QACf,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,oCAAmB,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAC1F,CAAC,CAAC;AAnCW,QAAA,gBAAgB,oBAmC3B"}

package/dist/fea/frame/index.d.ts

@@ -6,5 +6,5 @@ export { beam3dLocalStiffness, beam3dLocalStiffnessRightHanded, } from "./elemen
 export { memberAxes } from "./transformation/coordinate-transform";
 export { memberDemandEnvelope } from "./post/demand-envelope";
 export type { MemberDemandEnvelope, DemandEnvelopeOptions } from "./post/demand-envelope";
-export type { FrameId, FrameNode, DofRestraint, Restraint, DofRelease, MemberEndRelease, Member, NodalLoad, LoadDirection, UniformMemberLoad, PointMemberLoad, TrapezoidalMemberLoad, TorqueMemberLoad, MemberLoad, AreaDistribution, NullArea, AreaLoadDirection, AreaLoad, LoadCaseType, LoadCase, CombinationTerm, LoadCombination, FrameModel, } from "./types/model.type";
-export type { NodalDisplacement, NodalReaction, MemberEndForce, MemberEndForces, ReleaseDeformation, MemberReleaseDeformations, AutoRestrainedDof, MemberStation, StationComponent, ExtremeValue, MemberDiagram, GoverningValue, GoverningExtremeValue, Bound, FrameCaseResult, FrameEnvelope, FrameResult, } from "./types/results.type";
+export type { FrameId, FrameNode, DofRestraint, Restraint, LineRestraint, DofRelease, MemberEndRelease, Member, NodalLoad, LoadDirection, UniformMemberLoad, PointMemberLoad, TrapezoidalMemberLoad, TorqueMemberLoad, MemberLoad, AreaDistribution, NullArea, ShellType, ShellSection, ShellModifiers, ShellArea, ShellMesh, AreaLoadDirection, AreaLoad, LoadCaseType, LoadCase, CombinationTerm, LoadCombination, FrameModel, } from "./types/model.type";
+export type { NodalDisplacement, NodalReaction, MemberEndForce, MemberEndForces, ReleaseDeformation, MemberReleaseDeformations, AutoRestrainedDof, MemberStation, StationComponent, ExtremeValue, MemberDiagram, ShellCornerForces, ShellForceComponent, ShellElementForces, FrameMesh, GoverningValue, GoverningExtremeValue, Bound, FrameCaseResult, FrameEnvelope, FrameResult, } from "./types/results.type";

package/dist/fea/frame/mesh/expand-model.d.ts

@@ -0,0 +1,29 @@
+import type { Warnings } from "../../../core/types/output-message.type";
+import type { FrameId, FrameModel, Member } from "../types/model.type";
+import type { FrameMesh } from "../types/results.type";
+import { type ShellElement } from "./mesh-shells";
+/** A member with the generated mesh nodes lying inside it, where the solver splits it. */
+export interface SplitMember extends Member {
+    /** Generated node ids strictly between its ends, in order from the start. */
+    through?: FrameId[];
+}
+/** The model the solver analyses: shells meshed, line restraints spread to their nodes. */
+export interface ExpandedModel {
+    /**
+     * Nodes include the generated ones; `shells` are the elements; members list the generated
+     * nodes they pass through; restraints include the lines.
+     */
+    model: FrameModel & {
+        shells: ShellElement[];
+        members: SplitMember[];
+    };
+    /** The shell mesh for the result; undefined without shells. */
+    mesh: FrameMesh | undefined;
+    warnings: Warnings;
+}
+/**
+ * Meshes a valid model's shell areas, spreads its line restraints and marks the members that
+ * generated nodes split (see `meshShells`).
+ * Without shells and line restraints the model is returned as it is.
+ */
+export declare const expandModel: (model: FrameModel) => ExpandedModel;

package/dist/fea/frame/mesh/expand-model.js

@@ -0,0 +1,133 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.expandModel = void 0;
+const frame_constant_1 = require("../../../core/constants/frame.constant");
+const shell_axes_1 = require("../elements/shell/shell-axes");
+const mesh_shells_1 = require("./mesh-shells");
+const DOF_KEYS = ["ux", "uy", "uz", "rx", "ry", "rz"];
+/**
+ * Restraints of every node on each line restraint's segment (ends included), merged with the
+ * point restraints: a DOF fixed by either is fixed.
+ */
+const spreadLineRestraints = (model, nodes) => {
+    const lines = model.lineRestraints ?? [];
+    if (lines.length === 0)
+        return model.restraints;
+    const position = new Map(nodes.map((n) => [n.id, [n.x, n.y, n.z]]));
+    const merged = new Map(model.restraints.map((r) => [r.nodeId, { ...r }]));
+    for (const line of lines) {
+        const [P, Q] = line.nodeIds.map((id) => position.get(id));
+        for (const node of nodes) {
+            const on = line.nodeIds.includes(node.id) ||
+                (0, mesh_shells_1.pointOnSegment)([node.x, node.y, node.z], P, Q, true) !== null;
+            if (!on)
+                continue;
+            const restraint = merged.get(node.id) ?? { nodeId: node.id };
+            for (const key of DOF_KEYS)
+                if (line[key] === true)
+                    restraint[key] = true;
+            merged.set(node.id, restraint);
+        }
+    }
+    return [...merged.values()];
+};
+/**
+ * Warnings for member nodes lying inside a shell area but not on its mesh: such a node is not
+ * connected to the shell.
+ */
+const disconnectedNodes = (model, meshNodes) => {
+    const warnings = [];
+    const shells = model.shells ?? [];
+    const position = new Map(model.nodes.map((n) => [n.id, [n.x, n.y, n.z]]));
+    const memberNodes = new Set(model.members.flatMap((m) => [m.startNodeId, m.endNodeId]));
+    for (const shell of shells) {
+        const points = shell.nodeIds.map((id) => position.get(id));
+        const axes = (0, shell_axes_1.shellAxes)(points);
+        const size = Math.max(...axes.corners.map(([x, y]) => Math.hypot(x, y)));
+        const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * size;
+        for (const id of memberNodes) {
+            if (meshNodes.has(id))
+                continue;
+            const p = position.get(id);
+            const r = [p[0] - axes.origin[0], p[1] - axes.origin[1], p[2] - axes.origin[2]];
+            const [e1, e2, e3] = axes.lambda;
+            const z = r[0] * e3[0] + r[1] * e3[1] + r[2] * e3[2];
+            if (Math.abs(z) > Math.max(tolerance, ...axes.offsets.map(Math.abs)))
+                continue;
+            const x = r[0] * e1[0] + r[1] * e1[1] + r[2] * e1[2];
+            const y = r[0] * e2[0] + r[1] * e2[1] + r[2] * e2[2];
+            const inside = axes.corners.every(([ax, ay], k) => {
+                const [bx, by] = axes.corners[(k + 1) % 4];
+                return (bx - ax) * (y - ay) - (by - ay) * (x - ax) > tolerance * size;
+            });
+            if (inside) {
+                warnings.push({
+                    reference: "Shell area meshing",
+                    message: `Node "${String(id)}" lies inside shell area "${String(shell.id)}" but not on its mesh, so its members are not connected to the shell; divide the area through the node.`,
+                });
+            }
+        }
+    }
+    return warnings;
+};
+/**
+ * The generated nodes lying strictly inside each member (a beam along a meshed edge), so the
+ * member can be split there and share them with the shells.
+ */
+const splitMembers = (model, generated) => {
+    if (generated.length === 0)
+        return model.members;
+    const position = new Map(model.nodes.map((n) => [n.id, [n.x, n.y, n.z]]));
+    return model.members.map((member) => {
+        const [P, Q] = [position.get(member.startNodeId), position.get(member.endNodeId)];
+        const lo = [0, 1, 2].map((k) => Math.min(P[k], Q[k]));
+        const hi = [0, 1, 2].map((k) => Math.max(P[k], Q[k]));
+        const slack = frame_constant_1.AREA_GEOMETRY_TOLERANCE * Math.hypot(Q[0] - P[0], Q[1] - P[1], Q[2] - P[2]);
+        const through = [];
+        for (const node of generated) {
+            const x = [node.x, node.y, node.z];
+            if (x.some((v, k) => v < lo[k] - slack || v > hi[k] + slack))
+                continue;
+            const t = (0, mesh_shells_1.pointOnSegment)(x, P, Q);
+            if (t !== null)
+                through.push({ id: node.id, t });
+        }
+        return through.length === 0
+            ? member
+            : { ...member, through: through.sort((a, b) => a.t - b.t).map((h) => h.id) };
+    });
+};
+/**
+ * Meshes a valid model's shell areas, spreads its line restraints and marks the members that
+ * generated nodes split (see `meshShells`).
+ * Without shells and line restraints the model is returned as it is.
+ */
+const expandModel = (model) => {
+    const shells = model.shells ?? [];
+    if (shells.length === 0) {
+        return {
+            model: { ...model, shells: [], restraints: spreadLineRestraints(model, model.nodes) },
+            mesh: undefined,
+            warnings: [],
+        };
+    }
+    const meshed = (0, mesh_shells_1.meshShells)(model.nodes, shells);
+    const nodes = [...model.nodes, ...meshed.nodes];
+    const meshNodes = new Set(meshed.elements.flatMap((e) => e.nodeIds));
+    return {
+        model: {
+            ...model,
+            nodes,
+            members: splitMembers(model, meshed.nodes),
+            shells: meshed.elements,
+            restraints: spreadLineRestraints(model, nodes),
+        },
+        mesh: {
+            nodes: meshed.nodes,
+            shells: meshed.elements.map(({ id, areaId, nodeIds }) => ({ id, areaId, nodeIds })),
+        },
+        warnings: disconnectedNodes(model, meshNodes),
+    };
+};
+exports.expandModel = expandModel;
+//# sourceMappingURL=expand-model.js.map

package/dist/fea/frame/mesh/expand-model.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"expand-model.js","sourceRoot":"","sources":["../../../../src/fea/frame/mesh/expand-model.ts"],"names":[],"mappings":";;;AAAA,uEAA6E;AAG7E,6DAAoE;AAGpE,+CAA8E;AAE9E,MAAM,QAAQ,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AAoB/D;;;GAGG;AACH,MAAM,oBAAoB,GAAG,CAAC,KAAiB,EAAE,KAA0B,EAAe,EAAE;IAC1F,MAAM,KAAK,GAAG,KAAK,CAAC,cAAc,IAAI,EAAE,CAAC;IACzC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC,UAAU,CAAC;IAChD,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,MAAM,MAAM,GAAG,IAAI,GAAG,CAAqB,KAAK,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC9F,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,CAAC;QAC3D,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;YACzB,MAAM,EAAE,GACN,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;gBAC9B,IAAA,4BAAc,EAAC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,KAAK,IAAI,CAAC;YAChE,IAAI,CAAC,EAAE;gBAAE,SAAS;YAClB,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC;YAC7D,KAAK,MAAM,GAAG,IAAI,QAAQ;gBAAE,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,IAAI;oBAAE,SAAS,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;YAC1E,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,SAAS,CAAC,CAAC;QACjC,CAAC;IACH,CAAC;IACD,OAAO,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC;AAC9B,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,iBAAiB,GAAG,CAAC,KAAiB,EAAE,SAAuB,EAAY,EAAE;IACjF,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,EAAE,CAAC;IAClC,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,WAAW,EAAE,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IACxF,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,IAAA,sBAAS,EAAC,MAAM,CAAC,CAAC;QAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACzE,MAAM,SAAS,GAAG,wCAAuB,GAAG,IAAI,CAAC;QACjD,KAAK,MAAM,EAAE,IAAI,WAAW,EAAE,CAAC;YAC7B,IAAI,SAAS,CAAC,GAAG,CAAC,EAAE,CAAC;gBAAE,SAAS;YAChC,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC;YAC5B,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;YACtF,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;YACjC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACrD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;gBAAE,SAAS;YAC/E,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACrD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACrD,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE;gBAChD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC3C,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,SAAS,GAAG,IAAI,CAAC;YACxE,CAAC,CAAC,CAAC;YACH,IAAI,MAAM,EAAE,CAAC;gBACX,QAAQ,CAAC,IAAI,CAAC;oBACZ,SAAS,EAAE,oBAAoB;oBAC/B,OAAO,EAAE,SAAS,MAAM,CAAC,EAAE,CAAC,6BAA6B,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC,yGAAyG;iBACnL,CAAC,CAAC;YACL,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,YAAY,GAAG,CAAC,KAAiB,EAAE,SAAsB,EAAiB,EAAE;IAChF,IAAI,SAAS,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC,OAAO,CAAC;IACjD,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,OAAO,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QAClC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,WAAW,CAAE,EAAE,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAE,CAAC,CAAC;QACpF,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtD,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtD,MAAM,KAAK,GAAG,wCAAuB,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1F,MAAM,OAAO,GAAiC,EAAE,CAAC;QACjD,KAAK,MAAM,IAAI,IAAI,SAAS,EAAE,CAAC;YAC7B,MAAM,CAAC,GAAS,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;YACzC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;gBAAE,SAAS;YACvE,MAAM,CAAC,GAAG,IAAA,4BAAc,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAClC,IAAI,CAAC,KAAK,IAAI;gBAAE,OAAO,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;QACnD,CAAC;QACD,OAAO,OAAO,CAAC,MAAM,KAAK,CAAC;YACzB,CAAC,CAAC,MAAM;YACR,CAAC,CAAC,EAAE,GAAG,MAAM,EAAE,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;IACjF,CAAC,CAAC,CAAC;AACL,CAAC,CAAC;AAEF;;;;GAIG;AACI,MAAM,WAAW,GAAG,CAAC,KAAiB,EAAiB,EAAE;IAC9D,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,EAAE,CAAC;IAClC,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACxB,OAAO;YACL,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,MAAM,EAAE,EAAE,EAAE,UAAU,EAAE,oBAAoB,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,CAAC,EAAE;YACrF,IAAI,EAAE,SAAS;YACf,QAAQ,EAAE,EAAE;SACb,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,IAAA,wBAAU,EAAC,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IAChD,MAAM,SAAS,GAAG,IAAI,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;IACrE,OAAO;QACL,KAAK,EAAE;YACL,GAAG,KAAK;YACR,KAAK;YACL,OAAO,EAAE,YAAY,CAAC,KAAK,EAAE,MAAM,CAAC,KAAK,CAAC;YAC1C,MAAM,EAAE,MAAM,CAAC,QAAQ;YACvB,UAAU,EAAE,oBAAoB,CAAC,KAAK,EAAE,KAAK,CAAC;SAC/C;QACD,IAAI,EAAE;YACJ,KAAK,EAAE,MAAM,CAAC,KAAK;YACnB,MAAM,EAAE,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,CAAC;SACpF;QACD,QAAQ,EAAE,iBAAiB,CAAC,KAAK,EAAE,SAAS,CAAC;KAC9C,CAAC;AACJ,CAAC,CAAC;AA1BW,QAAA,WAAW,eA0BtB"}

package/dist/fea/frame/mesh/mesh-shells.d.ts

@@ -0,0 +1,38 @@
+import { type Vec3 } from "../elements/shell/shell-axes";
+import type { FrameId, FrameNode, ShellArea } from "../types/model.type";
+/** A shell element of the analysed model: a shell area, or one element of its mesh. */
+export interface ShellElement extends ShellArea {
+    /** The shell area it belongs to (its own id when the area is one element). */
+    areaId: FrameId;
+}
+/** The result of meshing a model's shell areas. */
+export interface ShellMeshResult {
+    /** Generated nodes, in global coordinates (mm). */
+    nodes: FrameNode[];
+    /** Every shell element, area by area in model order. */
+    elements: ShellElement[];
+    /** Areas that could not be meshed (no conforming division counts, a bad generated element, or a generated id already in use). */
+    problems: FrameId[];
+}
+/** Parameter t ∈ (0, 1) of `x` on the segment p → q when it lies strictly inside it, else null. */
+export declare const pointOnSegment: (x: Vec3, p: Vec3, q: Vec3, inclusive?: boolean) => number | null;
+/**
+ * Meshes the shell areas of a model whose shells and nodes are valid.
+ *
+ * - Each area edge is cut into segments at the model nodes lying on it (hard points), so an
+ *   edge shared whole by two areas, or partly where areas meet at a node (a T-junction), shares
+ *   its segments.
+ * - Every segment gets a division count: at least 1, at least `length / maxSize` for an area
+ *   with `mesh.maxSize`, and raised until in every area opposite edges add up to the same
+ *   count (and at least `mesh.divisions`), adding each division where the elements are longest.
+ * - Segment nodes are evenly spaced; interior nodes come from transfinite interpolation of the
+ *   four edges, so a planar area stays planar.
+ *
+ * An area that ends with 1 × 1 divisions is its own element, with its own id. Generated ids:
+ * `"<a>~<b>:<k>"` for the k-th node from model node a to b on a segment, `"<areaId>:<i>,<j>"`
+ * for interior nodes and for elements.
+ *
+ * @param nodes - Model nodes
+ * @param shells - Shell areas, each with 4 valid corner nodes
+ */
+export declare const meshShells: (nodes: FrameNode[], shells: ShellArea[]) => ShellMeshResult;