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

package/dist/esm/fea/frame/elements/shell/flat-shell.js

@@ -0,0 +1,111 @@
+import { zeros } from "../../../../utils/matrix.js";
+import { dkqPlateStiffness } from "./dkq-plate.js";
+import { membraneQuadStiffness } from "./membrane-quad.js";
+import { shellAxes } from "./shell-axes.js";
+/** Positions of the membrane DOFs [u, v, θz] and the plate DOFs [w, θx, θy] among a node's six. */
+const MEMBRANE_DOFS = [0, 1, 5];
+const PLATE_DOFS = [2, 3, 4];
+/**
+ * Local stiffness of a flat four-node shell: the drilling membrane (`membraneQuadStiffness`) on
+ * u₁, u₂, r₃ and the DKQ plate (`dkqPlateStiffness`) on u₃, r₁, r₂ of every corner. For a flat
+ * element of uniform section the two do not couple. `"membrane"` and `"plate"` keep only their
+ * own part, leaving the other DOFs without stiffness.
+ *
+ * @param corners - The four corners in the element plane (mm), counter-clockwise
+ * @param properties - Material, thickness, type and stiffness modifiers
+ * @returns The 24×24 stiffness over `[u₁ u₂ u₃ r₁ r₂ r₃]` per corner (N and mm)
+ */
+export const flatShellLocalStiffness = (corners, properties) => {
+    const { E, nu, t, type } = properties;
+    const k = zeros(24, 24);
+    const scatter = (part, dofs) => {
+        for (let a = 0; a < 12; a++) {
+            const i = 6 * Math.floor(a / 3) + dofs[a % 3];
+            for (let b = 0; b < 12; b++)
+                k[i][6 * Math.floor(b / 3) + dofs[b % 3]] += part[a][b];
+        }
+    };
+    if (type !== "plate") {
+        const membrane = { E, nu, t, modifiers: properties.membraneModifiers };
+        scatter(membraneQuadStiffness(corners, membrane), MEMBRANE_DOFS);
+    }
+    if (type !== "membrane") {
+        const plate = { E, nu, t, modifiers: properties.bendingModifiers };
+        scatter(dkqPlateStiffness(corners, plate), PLATE_DOFS);
+    }
+    return k;
+};
+/**
+ * Per corner, the 6×6 transformation from the global DOFs at the real node to the local DOFs at
+ * its projection on the mean plane: the direction cosines λ on translations and rotations, plus
+ * a rigid link across the offset h along local 3, `u_proj = u + θ × (−h·e₃)`, which adds
+ * `−h·r₂` to u₁ and `+h·r₁` to u₂. For a flat element it is just λ.
+ */
+export const shellNodeTransforms = ({ lambda, offsets }) => offsets.map((h) => {
+    const T = zeros(6, 6);
+    for (let i = 0; i < 3; i++) {
+        for (let j = 0; j < 3; j++) {
+            T[i][j] = lambda[i][j];
+            T[3 + i][3 + j] = lambda[i][j];
+        }
+    }
+    for (let j = 0; j < 3; j++) {
+        T[0][3 + j] = -h * lambda[1][j];
+        T[1][3 + j] = h * lambda[0][j];
+    }
+    return T;
+});
+/**
+ * Transforms a 24×24 local shell stiffness to global DOFs block by block,
+ * `k_ab = T_aᵀ·k_local,ab·T_b`, without forming the 24×24 transformation.
+ *
+ * @param kLocal - Stiffness over the local DOFs of the projected corners
+ * @param transforms - The four corner transformations from `shellNodeTransforms`
+ */
+export const shellGlobalStiffness = (kLocal, transforms) => {
+    const k = zeros(24, 24);
+    for (let a = 0; a < 4; a++) {
+        for (let b = 0; b < 4; b++) {
+            const Ta = transforms[a];
+            const Tb = transforms[b];
+            // block·T_b
+            const right = zeros(6, 6);
+            for (let i = 0; i < 6; i++) {
+                for (let m = 0; m < 6; m++) {
+                    const v = kLocal[6 * a + i][6 * b + m];
+                    if (v === 0)
+                        continue;
+                    for (let j = 0; j < 6; j++)
+                        right[i][j] += v * Tb[m][j];
+                }
+            }
+            // T_aᵀ·(block·T_b)
+            for (let i = 0; i < 6; i++) {
+                for (let m = 0; m < 6; m++) {
+                    const v = Ta[m][i];
+                    if (v === 0)
+                        continue;
+                    for (let j = 0; j < 6; j++)
+                        k[6 * a + i][6 * b + j] += v * right[m][j];
+                }
+            }
+        }
+    }
+    return k;
+};
+/**
+ * A flat four-node shell element (drilling membrane plus DKQ plate) built on the mean plane of
+ * its corners, linked rigidly to the real corners when they are not coplanar.
+ *
+ * @param points - The four corners in global axes (mm), in order around the element; local 3
+ *   follows the right-hand rule around them
+ * @param properties - Material, thickness, type and stiffness modifiers
+ * @param angle - Rotation of local 1 and 2 about local 3 in rad (default 0)
+ */
+export const flatShellElement = (points, properties, angle = 0) => {
+    const axes = shellAxes(points, angle);
+    const transforms = shellNodeTransforms(axes);
+    const kLocal = flatShellLocalStiffness(axes.corners, properties);
+    return { axes, transforms, kLocal, kGlobal: shellGlobalStiffness(kLocal, transforms) };
+};
+//# sourceMappingURL=flat-shell.js.map

package/dist/esm/fea/frame/elements/shell/flat-shell.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"flat-shell.js","sourceRoot":"","sources":["../../../../../../src/fea/frame/elements/shell/flat-shell.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAe,MAAM,mBAAmB,CAAC;AAEvD,OAAO,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AAChD,OAAO,EAAE,qBAAqB,EAAE,MAAM,iBAAiB,CAAC;AAGxD,OAAO,EAAE,SAAS,EAA6B,MAAM,cAAc,CAAC;AA4BpE,mGAAmG;AACnG,MAAM,aAAa,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAChC,MAAM,UAAU,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAE7B;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAG,CACrC,OAAwB,EACxB,UAA+B,EACvB,EAAE;IACV,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,UAAU,CAAC;IACtC,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,MAAM,OAAO,GAAG,CAAC,IAAY,EAAE,IAAc,EAAE,EAAE;QAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC9C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvF,CAAC;IACH,CAAC,CAAC;IACF,IAAI,IAAI,KAAK,OAAO,EAAE,CAAC;QACrB,MAAM,QAAQ,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,SAAS,EAAE,UAAU,CAAC,iBAAiB,EAAE,CAAC;QACvE,OAAO,CAAC,qBAAqB,CAAC,OAAO,EAAE,QAAQ,CAAC,EAAE,aAAa,CAAC,CAAC;IACnE,CAAC;IACD,IAAI,IAAI,KAAK,UAAU,EAAE,CAAC;QACxB,MAAM,KAAK,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,SAAS,EAAE,UAAU,CAAC,gBAAgB,EAAE,CAAC;QACnE,OAAO,CAAC,iBAAiB,CAAC,OAAO,EAAE,KAAK,CAAC,EAAE,UAAU,CAAC,CAAC;IACzD,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;;;GAKG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAG,CAAC,EAAE,MAAM,EAAE,OAAO,EAAa,EAAY,EAAE,CAC9E,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;IAChB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACvB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,CAAC;IACH,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAChC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACjC,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC,CAAC;AAEL;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAG,CAAC,MAAc,EAAE,UAAoB,EAAU,EAAE;IACnF,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;YACzB,MAAM,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;YACzB,YAAY;YACZ,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;oBACvC,IAAI,CAAC,KAAK,CAAC;wBAAE,SAAS;oBACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;wBAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1D,CAAC;YACH,CAAC;YACD,mBAAmB;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACnB,IAAI,CAAC,KAAK,CAAC;wBAAE,SAAS;oBACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;wBAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzE,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;;;;;;GAQG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAC9B,MAAuB,EACvB,UAA+B,EAC/B,KAAK,GAAG,CAAC,EACS,EAAE;IACpB,MAAM,IAAI,GAAG,SAAS,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;IACtC,MAAM,UAAU,GAAG,mBAAmB,CAAC,IAAI,CAAC,CAAC;IAC7C,MAAM,MAAM,GAAG,uBAAuB,CAAC,IAAI,CAAC,OAAO,EAAE,UAAU,CAAC,CAAC;IACjE,OAAO,EAAE,IAAI,EAAE,UAAU,EAAE,MAAM,EAAE,OAAO,EAAE,oBAAoB,CAAC,MAAM,EAAE,UAAU,CAAC,EAAE,CAAC;AACzF,CAAC,CAAC"}

package/dist/esm/fea/frame/elements/shell/membrane-quad.js

@@ -0,0 +1,107 @@
+import { DRILLING_PENALTY_RATIO } from "../../../../core/constants/frame.constant.js";
+import { zeros } from "../../../../utils/matrix.js";
+import { bilinearShape, cartesianDerivatives, GAUSS_3X3, midsideShape, quadJacobian, } from "./quad-geometry.js";
+import { isotropicRigidity } from "./rigidity.js";
+/**
+ * Strain–displacement rows of the four-node membrane with drilling rotations at (ξ, η).
+ *
+ * Displacements are bilinear plus Allman-type edge terms: for the edge from corner i to corner
+ * j (counter-clockwise) with mid-side function N_k,
+ * `u += N_k·(y_j − y_i)/8·(θ_j − θ_i)` and `v += N_k·(x_i − x_j)/8·(θ_j − θ_i)`, a quadratic
+ * bulge along the edge's outward normal that matches the difference of the corner rotations.
+ *
+ * @param corners - The four corners in the element plane (mm), counter-clockwise
+ * @param xi - Natural coordinate ξ
+ * @param eta - Natural coordinate η
+ */
+export const membraneStrainRows = (corners, xi, eta) => {
+    const jacobian = quadJacobian(corners, xi, eta);
+    const corner = bilinearShape(xi, eta);
+    const cornerD = cartesianDerivatives(jacobian, corner);
+    const edgeD = cartesianDerivatives(jacobian, midsideShape(xi, eta));
+    // u,x  u,y  v,x  v,y over the 12 DOFs
+    const ux = new Array(12).fill(0);
+    const uy = new Array(12).fill(0);
+    const vx = new Array(12).fill(0);
+    const vy = new Array(12).fill(0);
+    const rotation = new Array(12).fill(0);
+    for (let i = 0; i < 4; i++) {
+        ux[3 * i] = cornerD.dX[i];
+        uy[3 * i] = cornerD.dY[i];
+        vx[3 * i + 1] = cornerD.dX[i];
+        vy[3 * i + 1] = cornerD.dY[i];
+        rotation[3 * i + 2] = corner.N[i];
+    }
+    for (let k = 0; k < 4; k++) {
+        const i = k;
+        const j = (k + 1) % 4;
+        // (y_j − y_i, x_i − x_j) is the edge length times its outward normal
+        const cu = (corners[j][1] - corners[i][1]) / 8;
+        const cv = (corners[i][0] - corners[j][0]) / 8;
+        for (const [node, sign] of [
+            [j, 1],
+            [i, -1],
+        ]) {
+            const theta = 3 * node + 2;
+            ux[theta] += sign * cu * edgeD.dX[k];
+            uy[theta] += sign * cu * edgeD.dY[k];
+            vx[theta] += sign * cv * edgeD.dX[k];
+            vy[theta] += sign * cv * edgeD.dY[k];
+        }
+    }
+    for (let a = 0; a < 12; a++)
+        rotation[a] -= 0.5 * (vx[a] - uy[a]);
+    return { B: [ux, vy, uy.map((v, a) => v + vx[a])], rotation, detJ: jacobian.detJ };
+};
+/**
+ * Plane-stress membrane rigidity `E·t/(1 − ν²)·[[1, ν, 0], [ν, 1, 0], [0, 0, (1 − ν)/2]]`
+ * (N/mm), with the stiffness modifiers applied (see `isotropicRigidity`).
+ */
+export const membraneRigidity = ({ E, nu, t, modifiers }) => isotropicRigidity((E * t) / (1 - nu * nu), nu, modifiers);
+/**
+ * Local stiffness of the four-node membrane with drilling rotations (after Ibrahimbegovic,
+ * Taylor & Wilson 1990): bilinear displacements plus Allman-type edge terms, and the
+ * Hughes–Brezzi penalty `½·γ·t·∫(θ − ω)² dA` with `γ = penaltyRatio · G`.
+ *
+ * - The strain energy is integrated with 3×3 Gauss points. With 2×2 the edge terms leave a
+ *   spurious zero-energy mode in rectangular elements (alternating corner rotations whose
+ *   normal strains vanish exactly at the 2×2 points).
+ * - The penalty is integrated at the centre (one point). With 2×2 it over-constrains distorted
+ *   elements and they lock as γ grows; at the centre it only ties the element's mean rotation
+ *   to the material, so results hardly change with γ.
+ * - Stiffness modifiers scale the rigidity only; the penalty keeps the unmodified G.
+ *
+ * The DOFs are `[u₁ v₁ θ₁ … u₄ v₄ θ₄]`: in-plane displacements along the element's x and y and
+ * the rotation about its normal, right-handed, at corners listed counter-clockwise about that
+ * normal. The quad must be convex (det J > 0 throughout), which validation guarantees.
+ *
+ * @param corners - The four corners in the element plane (mm)
+ * @param properties - E (MPa), ν and thickness t (mm)
+ * @param penaltyRatio - γ / G (default `DRILLING_PENALTY_RATIO`)
+ * @returns The symmetric 12×12 stiffness in N and mm (N/mm, N, N·mm/rad); exactly three
+ *   zero-energy modes, the in-plane rigid-body motions
+ */
+export const membraneQuadStiffness = (corners, properties, penaltyRatio = DRILLING_PENALTY_RATIO) => {
+    const D = membraneRigidity(properties);
+    const k = zeros(12, 12);
+    for (const { xi, eta, weight } of GAUSS_3X3) {
+        const { B, detJ } = membraneStrainRows(corners, xi, eta);
+        const scale = detJ * weight;
+        // Bᵀ·D·B via D·B
+        const DB = D.map((row) => B[0].map((_, a) => row[0] * B[0][a] + row[1] * B[1][a] + row[2] * B[2][a]));
+        for (let a = 0; a < 12; a++) {
+            for (let b = 0; b < 12; b++) {
+                k[a][b] += (B[0][a] * DB[0][b] + B[1][a] * DB[1][b] + B[2][a] * DB[2][b]) * scale;
+            }
+        }
+    }
+    // one-point rule at the centre: weight 4
+    const { rotation, detJ } = membraneStrainRows(corners, 0, 0);
+    const penalty = ((penaltyRatio * properties.E * properties.t) / (2 * (1 + properties.nu))) * detJ * 4;
+    for (let a = 0; a < 12; a++) {
+        for (let b = 0; b < 12; b++)
+            k[a][b] += penalty * rotation[a] * rotation[b];
+    }
+    return k;
+};
+//# sourceMappingURL=membrane-quad.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"membrane-quad.js","sourceRoot":"","sources":["../../../../../../src/fea/frame/elements/shell/membrane-quad.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,sBAAsB,EAAE,MAAM,oCAAoC,CAAC;AAC5E,OAAO,EAAE,KAAK,EAAe,MAAM,mBAAmB,CAAC;AAEvD,OAAO,EACL,aAAa,EACb,oBAAoB,EACpB,SAAS,EACT,YAAY,EACZ,YAAY,GAEb,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,iBAAiB,EAA0B,MAAM,YAAY,CAAC;AAoBvE;;;;;;;;;;;GAWG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAChC,OAAwB,EACxB,EAAU,EACV,GAAW,EACS,EAAE;IACtB,MAAM,QAAQ,GAAG,YAAY,CAAC,OAAO,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IAChD,MAAM,MAAM,GAAG,aAAa,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IACtC,MAAM,OAAO,GAAG,oBAAoB,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC;IACvD,MAAM,KAAK,GAAG,oBAAoB,CAAC,QAAQ,EAAE,YAAY,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;IAEpE,sCAAsC;IACtC,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACzC,MAAM,QAAQ,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC1B,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC1B,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC9B,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC9B,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACpC,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACtB,qEAAqE;QACrE,MAAM,EAAE,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC/C,MAAM,EAAE,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC/C,KAAK,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI;YACzB,CAAC,CAAC,EAAE,CAAC,CAAC;YACN,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;SACR,EAAE,CAAC;YACF,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YAC3B,EAAE,CAAC,KAAK,CAAC,IAAI,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YACrC,EAAE,CAAC,KAAK,CAAC,IAAI,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YACrC,EAAE,CAAC,KAAK,CAAC,IAAI,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YACrC,EAAE,CAAC,KAAK,CAAC,IAAI,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACvC,CAAC;IACH,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;QAAE,QAAQ,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAElE,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,EAAE,IAAI,EAAE,QAAQ,CAAC,IAAI,EAAE,CAAC;AACrF,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,SAAS,EAAsB,EAAU,EAAE,CACtF,iBAAiB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,SAAS,CAAC,CAAC;AAE5D;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,MAAM,CAAC,MAAM,qBAAqB,GAAG,CACnC,OAAwB,EACxB,UAA8B,EAC9B,YAAY,GAAG,sBAAsB,EAC7B,EAAE;IACV,MAAM,CAAC,GAAG,gBAAgB,CAAC,UAAU,CAAC,CAAC;IACvC,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,KAAK,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,SAAS,EAAE,CAAC;QAC5C,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,kBAAkB,CAAC,OAAO,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QACzD,MAAM,KAAK,GAAG,IAAI,GAAG,MAAM,CAAC;QAC5B,iBAAiB;QACjB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CACvB,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAC3E,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC5B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;YACpF,CAAC;QACH,CAAC;IACH,CAAC;IAED,yCAAyC;IACzC,MAAM,EAAE,QAAQ,EAAE,IAAI,EAAE,GAAG,kBAAkB,CAAC,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7D,MAAM,OAAO,GACX,CAAC,CAAC,YAAY,GAAG,UAAU,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;IACxF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,OAAO,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IAC9E,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC"}

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

@@ -0,0 +1,127 @@
+/** 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 η.
+ */
+export const 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 η.
+ */
+export const 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. */
+export 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 };
+};
+/**
+ * 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.
+ */
+export 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],
+});
+/**
+ * 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.
+ */
+export 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 = midsideShape(xi, eta);
+    return { N: [...N, ...mid.N], dXi: [...dXi, ...mid.dXi], dEta: [...dEta, ...mid.dEta] };
+};
+/**
+ * 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
+ */
+export const quadJacobian = (corners, xi, eta) => {
+    const { dXi, dEta } = 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],
+        ],
+    };
+};
+/**
+ * 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`
+ */
+export 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]),
+});
+/**
+ * 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²
+ */
+export const cornerAreas = (corners) => {
+    const areas = [0, 0, 0, 0];
+    for (const { xi, eta, weight } of GAUSS_2X2) {
+        const { detJ } = quadJacobian(corners, xi, eta);
+        bilinearShape(xi, eta).N.forEach((n, i) => (areas[i] += n * detJ * weight));
+    }
+    return areas;
+};
+//# sourceMappingURL=quad-geometry.js.map

package/dist/esm/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;AACH,MAAM,CAAC,MAAM,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;AACH,MAAM,CAAC,MAAM,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;AACrE,MAAM,CAAC,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;AAEF;;;GAGG;AACH,MAAM,CAAC,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;AAEH;;;;GAIG;AACH,MAAM,CAAC,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,YAAY,CAAC,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;AAQF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,OAAwB,EAAE,EAAU,EAAE,GAAW,EAAgB,EAAE;IAC9F,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,aAAa,CAAC,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;AAEF;;;;GAIG;AACH,MAAM,CAAC,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;AAEH;;;;;;;GAOG;AACH,MAAM,CAAC,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,SAAS,EAAE,CAAC;QAC5C,MAAM,EAAE,IAAI,EAAE,GAAG,YAAY,CAAC,OAAO,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QAChD,aAAa,CAAC,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"}

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

@@ -0,0 +1,19 @@
+/**
+ * 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 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],
+    ];
+};
+//# sourceMappingURL=rigidity.js.map

package/dist/esm/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;AACH,MAAM,CAAC,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"}

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

@@ -0,0 +1,97 @@
+import { AREA_GEOMETRY_TOLERANCE, SHELL_WARPING_LIMIT, VERTICAL_MEMBER_TOLERANCE, } from "../../../../core/constants/frame.constant.js";
+import { quadJacobian } from "./quad-geometry.js";
+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)
+ */
+export 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]) < 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)),
+    };
+};
+/**
+ * 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 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 = 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 = shellAxes(points);
+    for (const [xi, eta] of [
+        [-1, -1],
+        [1, -1],
+        [1, 1],
+        [-1, 1],
+    ]) {
+        if (!(quadJacobian(axes.corners, xi, eta).detJ > (AREA_GEOMETRY_TOLERANCE * area) / 4)) {
+            return false;
+        }
+    }
+    return Math.max(...axes.offsets.map(Math.abs)) <= SHELL_WARPING_LIMIT * Math.sqrt(area);
+};
+//# sourceMappingURL=shell-axes.js.map

package/dist/esm/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,OAAO,EACL,uBAAuB,EACvB,mBAAmB,EACnB,yBAAyB,GAC1B,MAAM,oCAAoC,CAAC;AAG5C,OAAO,EAAE,YAAY,EAAa,MAAM,iBAAiB,CAAC;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;AACH,MAAM,CAAC,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,yBAAyB,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;AAEF;;;;GAIG;AACH,MAAM,CAAC,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,uBAAuB,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,SAAS,CAAC,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,YAAY,CAAC,IAAI,CAAC,OAAO,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,uBAAuB,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,mBAAmB,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAC1F,CAAC,CAAC"}

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

@@ -0,0 +1,129 @@
+import { AREA_GEOMETRY_TOLERANCE } from "../../../core/constants/frame.constant.js";
+import { shellAxes } from "../elements/shell/shell-axes.js";
+import { meshShells, pointOnSegment } from "./mesh-shells.js";
+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) ||
+                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 = shellAxes(points);
+        const size = Math.max(...axes.corners.map(([x, y]) => Math.hypot(x, y)));
+        const tolerance = 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 = 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 = 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.
+ */
+export 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 = 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),
+    };
+};
+//# sourceMappingURL=expand-model.js.map