@theerapat-s28/ts-structural-eng-tools 3.3.0 → 3.4.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/CHANGELOG.md +18 -1
- package/README.md +65 -11
- package/dist/core/constants/fea.constant.d.ts +21 -0
- package/dist/core/constants/fea.constant.js +26 -0
- package/dist/core/constants/fea.constant.js.map +1 -0
- package/dist/core/constants/frame.constant.d.ts +0 -34
- package/dist/core/constants/frame.constant.js +2 -35
- package/dist/core/constants/frame.constant.js.map +1 -1
- package/dist/core/constants/rc.constant.d.ts +8 -0
- package/dist/core/constants/rc.constant.js +10 -1
- package/dist/core/constants/rc.constant.js.map +1 -1
- package/dist/core/constants/shell.constant.d.ts +23 -0
- package/dist/core/constants/shell.constant.js +28 -0
- package/dist/core/constants/shell.constant.js.map +1 -0
- package/dist/core/errors/fea.error.d.ts +5 -0
- package/dist/core/errors/fea.error.js +9 -1
- package/dist/core/errors/fea.error.js.map +1 -1
- package/dist/core/index.d.ts +4 -2
- package/dist/core/index.js +16 -10
- package/dist/core/index.js.map +1 -1
- package/dist/esm/core/constants/fea.constant.js +23 -0
- package/dist/esm/core/constants/fea.constant.js.map +1 -0
- package/dist/esm/core/constants/frame.constant.js +1 -34
- package/dist/esm/core/constants/frame.constant.js.map +1 -1
- package/dist/esm/core/constants/rc.constant.js +9 -0
- package/dist/esm/core/constants/rc.constant.js.map +1 -1
- package/dist/esm/core/constants/shell.constant.js +25 -0
- package/dist/esm/core/constants/shell.constant.js.map +1 -0
- package/dist/esm/core/errors/fea.error.js +9 -1
- package/dist/esm/core/errors/fea.error.js.map +1 -1
- package/dist/esm/core/index.js +4 -2
- package/dist/esm/core/index.js.map +1 -1
- package/dist/esm/fea/analysis/mesh/expand-model.js +1 -1
- package/dist/esm/fea/analysis/mesh/expand-model.js.map +1 -1
- package/dist/esm/fea/analysis/mesh/member-connectivity.js +2 -1
- package/dist/esm/fea/analysis/mesh/member-connectivity.js.map +1 -1
- package/dist/esm/fea/analysis/post/envelope.js +70 -0
- package/dist/esm/fea/analysis/post/envelope.js.map +1 -1
- package/dist/esm/fea/analysis/post/present-state.js +9 -2
- package/dist/esm/fea/analysis/post/present-state.js.map +1 -1
- package/dist/esm/fea/analysis/post/wall-forces.js +299 -0
- package/dist/esm/fea/analysis/post/wall-forces.js.map +1 -0
- package/dist/esm/fea/analysis/post/wall-geometry.js +173 -0
- package/dist/esm/fea/analysis/post/wall-geometry.js.map +1 -0
- package/dist/esm/fea/analysis/solver/analysis-system.js +1 -1
- package/dist/esm/fea/analysis/solver/analysis-system.js.map +1 -1
- package/dist/esm/fea/analysis/solver/analyze.js +13 -5
- package/dist/esm/fea/analysis/solver/analyze.js.map +1 -1
- package/dist/esm/fea/analysis/solver/load-case-solver.js +4 -3
- package/dist/esm/fea/analysis/solver/load-case-solver.js.map +1 -1
- package/dist/esm/fea/analysis/validation/validate-model.js +20 -2
- package/dist/esm/fea/analysis/validation/validate-model.js.map +1 -1
- package/dist/esm/fea/analysis/validation/validate-walls.js +226 -0
- package/dist/esm/fea/analysis/validation/validate-walls.js.map +1 -0
- package/dist/esm/fea/frame/area/area-geometry.js +1 -1
- package/dist/esm/fea/frame/area/area-geometry.js.map +1 -1
- package/dist/esm/fea/frame/area/area-member-loads.js +1 -1
- package/dist/esm/fea/frame/area/area-member-loads.js.map +1 -1
- package/dist/esm/fea/frame/area/edge-members.js +1 -1
- package/dist/esm/fea/frame/area/edge-members.js.map +1 -1
- package/dist/esm/fea/frame/post/member-diagram.js +12 -0
- package/dist/esm/fea/frame/post/member-diagram.js.map +1 -1
- package/dist/esm/fea/frame/solver/nodes-along-members.js +1 -1
- package/dist/esm/fea/frame/solver/nodes-along-members.js.map +1 -1
- package/dist/esm/fea/frame/transformation/coordinate-transform.js +1 -1
- package/dist/esm/fea/frame/transformation/coordinate-transform.js.map +1 -1
- package/dist/esm/fea/frame/validation/validate-null-areas.js +1 -1
- package/dist/esm/fea/frame/validation/validate-null-areas.js.map +1 -1
- package/dist/esm/fea/shared/solver/auto-restraint.js +1 -1
- package/dist/esm/fea/shared/solver/auto-restraint.js.map +1 -1
- package/dist/esm/fea/shared/validation/issues.js +8 -0
- package/dist/esm/fea/shared/validation/issues.js.map +1 -1
- package/dist/esm/fea/shell/elements/membrane-quad.js +1 -1
- package/dist/esm/fea/shell/elements/shell-axes.js +2 -1
- package/dist/esm/fea/shell/elements/shell-axes.js.map +1 -1
- package/dist/esm/fea/shell/mesh/mesh-shells.js +1 -1
- package/dist/esm/fea/shell/mesh/mesh-shells.js.map +1 -1
- package/dist/esm/fea/shell/post/strip-forces.js +341 -0
- package/dist/esm/fea/shell/post/strip-forces.js.map +1 -0
- package/dist/esm/fea/shell/solver/shell-loads.js +3 -0
- package/dist/esm/fea/shell/solver/shell-loads.js.map +1 -1
- package/dist/esm/fea/shell/validation/validate-shells.js +9 -2
- package/dist/esm/fea/shell/validation/validate-shells.js.map +1 -1
- package/dist/esm/fea/shell/validation/validate-strips.js +99 -0
- package/dist/esm/fea/shell/validation/validate-strips.js.map +1 -0
- package/dist/esm/index.js +1 -1
- package/dist/esm/index.js.map +1 -1
- package/dist/esm/rc/rc-section-properties.js +1 -1
- package/dist/esm/rc/rc-section-properties.js.map +1 -1
- package/dist/fea/analysis/index.d.ts +2 -2
- package/dist/fea/analysis/mesh/expand-model.js +3 -3
- package/dist/fea/analysis/mesh/expand-model.js.map +1 -1
- package/dist/fea/analysis/mesh/member-connectivity.js +2 -1
- package/dist/fea/analysis/mesh/member-connectivity.js.map +1 -1
- package/dist/fea/analysis/post/envelope.js +70 -0
- package/dist/fea/analysis/post/envelope.js.map +1 -1
- package/dist/fea/analysis/post/present-state.d.ts +7 -0
- package/dist/fea/analysis/post/present-state.js +10 -3
- package/dist/fea/analysis/post/present-state.js.map +1 -1
- package/dist/fea/analysis/post/wall-forces.d.ts +65 -0
- package/dist/fea/analysis/post/wall-forces.js +305 -0
- package/dist/fea/analysis/post/wall-forces.js.map +1 -0
- package/dist/fea/analysis/post/wall-geometry.d.ts +73 -0
- package/dist/fea/analysis/post/wall-geometry.js +188 -0
- package/dist/fea/analysis/post/wall-geometry.js.map +1 -0
- package/dist/fea/analysis/solver/analysis-system.js +2 -2
- package/dist/fea/analysis/solver/analysis-system.js.map +1 -1
- package/dist/fea/analysis/solver/analyze.d.ts +6 -3
- package/dist/fea/analysis/solver/analyze.js +13 -5
- package/dist/fea/analysis/solver/analyze.js.map +1 -1
- package/dist/fea/analysis/solver/load-case-solver.d.ts +1 -0
- package/dist/fea/analysis/solver/load-case-solver.js +5 -4
- package/dist/fea/analysis/solver/load-case-solver.js.map +1 -1
- package/dist/fea/analysis/types/model.type.d.ts +58 -1
- package/dist/fea/analysis/types/results.type.d.ts +125 -1
- package/dist/fea/analysis/validation/validate-model.d.ts +1 -1
- package/dist/fea/analysis/validation/validate-model.js +20 -2
- package/dist/fea/analysis/validation/validate-model.js.map +1 -1
- package/dist/fea/analysis/validation/validate-walls.d.ts +35 -0
- package/dist/fea/analysis/validation/validate-walls.js +230 -0
- package/dist/fea/analysis/validation/validate-walls.js.map +1 -0
- package/dist/fea/frame/area/area-geometry.js +5 -5
- package/dist/fea/frame/area/area-geometry.js.map +1 -1
- package/dist/fea/frame/area/area-member-loads.js +7 -7
- package/dist/fea/frame/area/area-member-loads.js.map +1 -1
- package/dist/fea/frame/area/edge-members.js +2 -2
- package/dist/fea/frame/area/edge-members.js.map +1 -1
- package/dist/fea/frame/post/member-diagram.d.ts +9 -0
- package/dist/fea/frame/post/member-diagram.js +14 -1
- package/dist/fea/frame/post/member-diagram.js.map +1 -1
- package/dist/fea/frame/solver/nodes-along-members.js +2 -2
- package/dist/fea/frame/solver/nodes-along-members.js.map +1 -1
- package/dist/fea/frame/transformation/coordinate-transform.js +2 -2
- package/dist/fea/frame/transformation/coordinate-transform.js.map +1 -1
- package/dist/fea/frame/validation/validate-null-areas.js +2 -2
- package/dist/fea/frame/validation/validate-null-areas.js.map +1 -1
- package/dist/fea/shared/solver/auto-restraint.js +3 -3
- package/dist/fea/shared/solver/auto-restraint.js.map +1 -1
- package/dist/fea/shared/validation/issues.js +8 -0
- package/dist/fea/shared/validation/issues.js.map +1 -1
- package/dist/fea/shell/elements/membrane-quad.js +2 -2
- package/dist/fea/shell/elements/shell-axes.js +6 -5
- package/dist/fea/shell/elements/shell-axes.js.map +1 -1
- package/dist/fea/shell/index.d.ts +2 -2
- package/dist/fea/shell/mesh/mesh-shells.js +2 -2
- package/dist/fea/shell/mesh/mesh-shells.js.map +1 -1
- package/dist/fea/shell/post/strip-forces.d.ts +69 -0
- package/dist/fea/shell/post/strip-forces.js +349 -0
- package/dist/fea/shell/post/strip-forces.js.map +1 -0
- package/dist/fea/shell/solver/shell-loads.d.ts +3 -1
- package/dist/fea/shell/solver/shell-loads.js +3 -0
- package/dist/fea/shell/solver/shell-loads.js.map +1 -1
- package/dist/fea/shell/types/model.type.d.ts +31 -0
- package/dist/fea/shell/types/results.type.d.ts +54 -0
- package/dist/fea/shell/validation/validate-shells.d.ts +7 -2
- package/dist/fea/shell/validation/validate-shells.js +9 -2
- package/dist/fea/shell/validation/validate-shells.js.map +1 -1
- package/dist/fea/shell/validation/validate-strips.d.ts +11 -0
- package/dist/fea/shell/validation/validate-strips.js +103 -0
- package/dist/fea/shell/validation/validate-strips.js.map +1 -0
- package/dist/index.d.ts +2 -2
- package/dist/index.js +6 -2
- package/dist/index.js.map +1 -1
- package/dist/rc/rc-section-properties.js +4 -4
- package/dist/rc/rc-section-properties.js.map +1 -1
- package/package.json +1 -1
|
@@ -0,0 +1,299 @@
|
|
|
1
|
+
import { VERTICAL_MEMBER_TOLERANCE } from "../../../core/constants/fea.constant.js";
|
|
2
|
+
import { memberInternalForces } from "../../frame/post/member-diagram.js";
|
|
3
|
+
import { elementDofs } from "../../shared/solver/global-assembly.js";
|
|
4
|
+
import { roundToDecimalPlaces } from "../../../utils/math.js";
|
|
5
|
+
import { cross, dot, fromNegativeSide, sub, wallAxes, wallCutPositions, wallExtent, wallTolerance, } from "./wall-geometry.js";
|
|
6
|
+
/** The resultants, in the order of the raw values. */
|
|
7
|
+
export const WALL_COMPONENTS = ["P", "V2", "V3", "T", "M2", "M3"];
|
|
8
|
+
const FORCE_DECIMALS = 4;
|
|
9
|
+
const LENGTH_DECIMALS = 3;
|
|
10
|
+
const round = (value, decimals = FORCE_DECIMALS) => roundToDecimalPlaces(value, decimals) + 0; // + 0: no −0
|
|
11
|
+
const roundPoint = (p) => [round(p[0], 3), round(p[1], 3), round(p[2], 3)];
|
|
12
|
+
const range = (values) => ({ min: Math.min(...values), max: Math.max(...values) });
|
|
13
|
+
/**
|
|
14
|
+
* One cut: the shell elements and member pieces on its side that reach it, the members that
|
|
15
|
+
* cross it, its section (a leg per area, a point per member) and the section's centroid.
|
|
16
|
+
*/
|
|
17
|
+
const prepareCut = (system, shape, e1, s, face, tol) => {
|
|
18
|
+
const negative = fromNegativeSide(face);
|
|
19
|
+
const sign = negative ? 1 : -1;
|
|
20
|
+
const at = (n) => {
|
|
21
|
+
const node = system.nodes[n];
|
|
22
|
+
return [node.x, node.y, node.z];
|
|
23
|
+
};
|
|
24
|
+
const shells = [];
|
|
25
|
+
const legs = [];
|
|
26
|
+
for (const { areaId, shells: list } of shape.areas) {
|
|
27
|
+
const onCut = [];
|
|
28
|
+
let normal;
|
|
29
|
+
let thickness = 0;
|
|
30
|
+
for (const index of list) {
|
|
31
|
+
const shell = system.shells[index];
|
|
32
|
+
const points = shell.nodeIndices.map(at);
|
|
33
|
+
const r = range(points.map((p) => dot(p, e1)));
|
|
34
|
+
if (r.max - r.min <= tol)
|
|
35
|
+
continue;
|
|
36
|
+
if (Math.abs((negative ? r.max : r.min) - s) > tol)
|
|
37
|
+
continue;
|
|
38
|
+
const corners = points.flatMap((p, i) => (Math.abs(dot(p, e1) - s) <= tol ? [i] : []));
|
|
39
|
+
shells.push({ shell: index, corners, sign });
|
|
40
|
+
onCut.push(...corners.map((i) => points[i]));
|
|
41
|
+
normal ?? (normal = shell.element.axes.lambda[2]);
|
|
42
|
+
thickness = shell.section.thickness;
|
|
43
|
+
}
|
|
44
|
+
if (!normal)
|
|
45
|
+
continue;
|
|
46
|
+
const d = cross(e1, normal);
|
|
47
|
+
const dn = Math.hypot(d[0], d[1], d[2]);
|
|
48
|
+
if (dn === 0)
|
|
49
|
+
continue;
|
|
50
|
+
// legs run up (spandrel), else toward +X (+Y for a leg parallel to Y)
|
|
51
|
+
let along = [d[0] / dn, d[1] / dn, d[2] / dn];
|
|
52
|
+
const flip = Math.abs(along[2]) > VERTICAL_MEMBER_TOLERANCE
|
|
53
|
+
? along[2] < 0
|
|
54
|
+
: Math.abs(along[0]) > VERTICAL_MEMBER_TOLERANCE
|
|
55
|
+
? along[0] < 0
|
|
56
|
+
: along[1] < 0;
|
|
57
|
+
if (flip)
|
|
58
|
+
along = [-along[0], -along[1], -along[2]];
|
|
59
|
+
let start = onCut[0];
|
|
60
|
+
let end = onCut[0];
|
|
61
|
+
for (const p of onCut) {
|
|
62
|
+
if (dot(p, along) < dot(start, along))
|
|
63
|
+
start = p;
|
|
64
|
+
if (dot(p, along) > dot(end, along))
|
|
65
|
+
end = p;
|
|
66
|
+
}
|
|
67
|
+
const length = dot(sub(end, start), along);
|
|
68
|
+
if (length > tol)
|
|
69
|
+
legs.push({ areaId, start, end, length, thickness });
|
|
70
|
+
}
|
|
71
|
+
const ends = [];
|
|
72
|
+
const crossings = [];
|
|
73
|
+
const members = [];
|
|
74
|
+
for (const { memberId, pieces } of shape.members) {
|
|
75
|
+
for (const piece of pieces) {
|
|
76
|
+
const member = system.members[piece];
|
|
77
|
+
const p0 = at(member.startIndex);
|
|
78
|
+
const p1 = at(member.endIndex);
|
|
79
|
+
const s0 = dot(p0, e1);
|
|
80
|
+
const s1 = dot(p1, e1);
|
|
81
|
+
if (Math.abs(s1 - s0) <= tol)
|
|
82
|
+
continue; // lies in the cut plane
|
|
83
|
+
const r = range([s0, s1]);
|
|
84
|
+
if (r.min < s - tol && r.max > s + tol) {
|
|
85
|
+
const t = (s - s0) / (s1 - s0);
|
|
86
|
+
const point = [
|
|
87
|
+
p0[0] + t * (p1[0] - p0[0]),
|
|
88
|
+
p0[1] + t * (p1[1] - p0[1]),
|
|
89
|
+
p0[2] + t * (p1[2] - p0[2]),
|
|
90
|
+
];
|
|
91
|
+
crossings.push({ piece, x: t * member.length, point, sign: s1 > s0 ? 1 : -1 });
|
|
92
|
+
members.push({ memberId, point, A: member.section.A });
|
|
93
|
+
break;
|
|
94
|
+
}
|
|
95
|
+
if (Math.abs((negative ? r.max : r.min) - s) <= tol) {
|
|
96
|
+
// the piece's end on the cut: the one with the larger s from below, the smaller from above
|
|
97
|
+
const end = negative === s1 > s0 ? 1 : 0;
|
|
98
|
+
const point = end === 1 ? p1 : p0;
|
|
99
|
+
ends.push({ piece, end, point, sign });
|
|
100
|
+
members.push({ memberId, point, A: member.section.A });
|
|
101
|
+
break;
|
|
102
|
+
}
|
|
103
|
+
}
|
|
104
|
+
}
|
|
105
|
+
// centroid: legs by length × thickness at their midpoints, members by A
|
|
106
|
+
let area = 0;
|
|
107
|
+
const sum = [0, 0, 0];
|
|
108
|
+
const addWeighted = (p, w) => {
|
|
109
|
+
area += w;
|
|
110
|
+
for (let k = 0; k < 3; k++)
|
|
111
|
+
sum[k] += w * p[k];
|
|
112
|
+
};
|
|
113
|
+
for (const leg of legs) {
|
|
114
|
+
const mid = [0, 1, 2].map((k) => (leg.start[k] + leg.end[k]) / 2);
|
|
115
|
+
addWeighted(mid, leg.length * leg.thickness);
|
|
116
|
+
}
|
|
117
|
+
for (const m of members)
|
|
118
|
+
addWeighted(m.point, m.A);
|
|
119
|
+
let centroid;
|
|
120
|
+
if (area > 0) {
|
|
121
|
+
centroid = sum.map((v) => v / area);
|
|
122
|
+
}
|
|
123
|
+
else {
|
|
124
|
+
// only corners touch the cut: the mean of the points on it
|
|
125
|
+
const points = [
|
|
126
|
+
...shells.flatMap((t) => t.corners.map((i) => at(system.shells[t.shell].nodeIndices[i]))),
|
|
127
|
+
...ends.map((t) => t.point),
|
|
128
|
+
...crossings.map((t) => t.point),
|
|
129
|
+
];
|
|
130
|
+
centroid = [0, 1, 2].map((k) => points.reduce((acc, p) => acc + p[k], 0) / Math.max(points.length, 1));
|
|
131
|
+
}
|
|
132
|
+
return {
|
|
133
|
+
face,
|
|
134
|
+
centroid,
|
|
135
|
+
section: {
|
|
136
|
+
legs: legs.map((leg) => ({
|
|
137
|
+
areaId: leg.areaId,
|
|
138
|
+
start: roundPoint(leg.start),
|
|
139
|
+
end: roundPoint(leg.end),
|
|
140
|
+
length: round(leg.length, LENGTH_DECIMALS),
|
|
141
|
+
thickness: leg.thickness,
|
|
142
|
+
})),
|
|
143
|
+
members: members.map((m) => ({ memberId: m.memberId, point: roundPoint(m.point), A: m.A })),
|
|
144
|
+
length: round(legs.reduce((acc, leg) => acc + leg.length, 0), LENGTH_DECIMALS),
|
|
145
|
+
area: round(area, LENGTH_DECIMALS),
|
|
146
|
+
centroid: roundPoint(centroid),
|
|
147
|
+
},
|
|
148
|
+
shells,
|
|
149
|
+
ends,
|
|
150
|
+
crossings,
|
|
151
|
+
};
|
|
152
|
+
};
|
|
153
|
+
/**
|
|
154
|
+
* Prepares the piers and spandrels of a model that passed `validateModel`: their axes, cuts
|
|
155
|
+
* (see `wallCutPositions`) and, for each cut, the shell elements and member pieces whose nodal
|
|
156
|
+
* forces make up its resultant.
|
|
157
|
+
*
|
|
158
|
+
* @param walls - The model's piers and spandrels
|
|
159
|
+
*/
|
|
160
|
+
export const prepareWalls = (walls, system) => {
|
|
161
|
+
const at = (n) => {
|
|
162
|
+
const node = system.nodes[n];
|
|
163
|
+
return [node.x, node.y, node.z];
|
|
164
|
+
};
|
|
165
|
+
const prepare = (kind, wall) => {
|
|
166
|
+
const areas = (wall.areaIds ?? []).map((areaId) => ({
|
|
167
|
+
areaId,
|
|
168
|
+
shells: system.shellsByArea.get(areaId) ?? [],
|
|
169
|
+
}));
|
|
170
|
+
const members = (wall.memberIds ?? []).map((memberId) => ({
|
|
171
|
+
memberId,
|
|
172
|
+
pieces: system.memberGroups[system.memberIndex.get(memberId)].pieces,
|
|
173
|
+
}));
|
|
174
|
+
const shape = {
|
|
175
|
+
areas: areas.map(({ areaId, shells }) => ({
|
|
176
|
+
areaId,
|
|
177
|
+
elements: shells.map((s) => system.shells[s].nodeIndices.map(at)),
|
|
178
|
+
})),
|
|
179
|
+
members: members.map(({ memberId, pieces }) => ({
|
|
180
|
+
memberId,
|
|
181
|
+
pieces: pieces.map((m) => [
|
|
182
|
+
at(system.members[m].startIndex),
|
|
183
|
+
at(system.members[m].endIndex),
|
|
184
|
+
]),
|
|
185
|
+
})),
|
|
186
|
+
};
|
|
187
|
+
const axes = wallAxes(kind, shape, {
|
|
188
|
+
angle: wall.angle,
|
|
189
|
+
reverse: wall.reverse,
|
|
190
|
+
});
|
|
191
|
+
const e1 = axes[0];
|
|
192
|
+
const extent = wallExtent(shape, e1);
|
|
193
|
+
const tol = wallTolerance(extent);
|
|
194
|
+
const offset = kind === "pier" ? 0 : extent.min;
|
|
195
|
+
const extra = (wall.cuts ?? []).map((c) => c + offset);
|
|
196
|
+
const cuts = wallCutPositions(kind, shape, e1, extra).map(({ s, face }) => ({
|
|
197
|
+
position: s - offset,
|
|
198
|
+
...prepareCut(system, { areas, members }, e1, s, face, tol),
|
|
199
|
+
}));
|
|
200
|
+
return { id: wall.id, axes, length: extent.max - extent.min, cuts };
|
|
201
|
+
};
|
|
202
|
+
return {
|
|
203
|
+
piers: (walls.piers ?? []).map((p) => prepare("pier", p)),
|
|
204
|
+
spandrels: (walls.spandrels ?? []).map((s) => prepare("spandrel", s)),
|
|
205
|
+
};
|
|
206
|
+
};
|
|
207
|
+
/**
|
|
208
|
+
* The unrounded global force (N) and moment about the cut's centroid (N·mm) on the +1 face of
|
|
209
|
+
* the negative part of one cut: the nodal forces `k·u − f` of the shell elements on the cut's
|
|
210
|
+
* side at their corners on it, the end forces of the member pieces ending on it from that side
|
|
211
|
+
* (both negated when the side is the positive one), and the internal forces of the members
|
|
212
|
+
* crossing it. Linear in the state.
|
|
213
|
+
*/
|
|
214
|
+
const cutRawForces = (system, state, cut) => {
|
|
215
|
+
const F = [0, 0, 0];
|
|
216
|
+
const M = [0, 0, 0];
|
|
217
|
+
const c = cut.centroid;
|
|
218
|
+
const addAt = (p, f, m, sign) => {
|
|
219
|
+
const r = sub(p, c);
|
|
220
|
+
const rf = cross(r, f);
|
|
221
|
+
for (let k = 0; k < 3; k++) {
|
|
222
|
+
F[k] += sign * f[k];
|
|
223
|
+
M[k] += sign * (m[k] + rf[k]);
|
|
224
|
+
}
|
|
225
|
+
};
|
|
226
|
+
for (const { shell: index, corners, sign } of cut.shells) {
|
|
227
|
+
const shell = system.shells[index];
|
|
228
|
+
const dofs = elementDofs(shell.nodeIndices);
|
|
229
|
+
const u = dofs.map((d) => state.U[d]);
|
|
230
|
+
const p = state.shellPressures[index];
|
|
231
|
+
for (const i of corners) {
|
|
232
|
+
const f = [0, 0, 0, 0, 0, 0];
|
|
233
|
+
for (let r = 0; r < 6; r++) {
|
|
234
|
+
const row = shell.element.kGlobal[6 * i + r];
|
|
235
|
+
let v = 0;
|
|
236
|
+
for (let j = 0; j < 24; j++)
|
|
237
|
+
v += row[j] * u[j];
|
|
238
|
+
f[r] = v;
|
|
239
|
+
}
|
|
240
|
+
for (let k = 0; k < 3; k++)
|
|
241
|
+
f[k] -= p[k] * 1e-3 * shell.cornerAreas[i];
|
|
242
|
+
const node = system.nodes[shell.nodeIndices[i]];
|
|
243
|
+
addAt([node.x, node.y, node.z], [f[0], f[1], f[2]], [f[3], f[4], f[5]], sign);
|
|
244
|
+
}
|
|
245
|
+
}
|
|
246
|
+
const toGlobal = (lambda, local) => [
|
|
247
|
+
lambda[0][0] * local[0] + lambda[1][0] * local[1] + lambda[2][0] * local[2],
|
|
248
|
+
lambda[0][1] * local[0] + lambda[1][1] * local[1] + lambda[2][1] * local[2],
|
|
249
|
+
lambda[0][2] * local[0] + lambda[1][2] * local[1] + lambda[2][2] * local[2],
|
|
250
|
+
];
|
|
251
|
+
for (const { piece, end, point, sign } of cut.ends) {
|
|
252
|
+
const { lambda } = system.members[piece];
|
|
253
|
+
const f = state.memberForces[piece].slice(6 * end, 6 * end + 6);
|
|
254
|
+
addAt(point, toGlobal(lambda, f.slice(0, 3)), toGlobal(lambda, f.slice(3, 6)), sign);
|
|
255
|
+
}
|
|
256
|
+
for (const { piece, x, point, sign } of cut.crossings) {
|
|
257
|
+
const { lambda } = system.members[piece];
|
|
258
|
+
const f = memberInternalForces(state.memberForces[piece].slice(0, 6), state.memberLoads[piece], x);
|
|
259
|
+
addAt(point, toGlobal(lambda, f.slice(0, 3)), toGlobal(lambda, f.slice(3, 6)), sign);
|
|
260
|
+
}
|
|
261
|
+
return { F, M };
|
|
262
|
+
};
|
|
263
|
+
/** The presented resultant of one cut in the wall's axes (kN, kN·m). */
|
|
264
|
+
const cutResultants = (system, state, axes, cut) => {
|
|
265
|
+
const { F, M } = cutRawForces(system, state, cut);
|
|
266
|
+
return {
|
|
267
|
+
P: round(dot(F, axes[0]) * 1e-3),
|
|
268
|
+
V2: round(dot(F, axes[1]) * 1e-3),
|
|
269
|
+
V3: round(dot(F, axes[2]) * 1e-3),
|
|
270
|
+
T: round(dot(M, axes[0]) * 1e-6),
|
|
271
|
+
M2: round(dot(M, axes[1]) * 1e-6),
|
|
272
|
+
M3: round(dot(M, axes[2]) * 1e-6),
|
|
273
|
+
};
|
|
274
|
+
};
|
|
275
|
+
const presentAxes = (axes) => axes.map((e) => e.map((v) => round(v, 9)));
|
|
276
|
+
/** The presented forces of one pier in a raw state. */
|
|
277
|
+
export const presentPierForces = (system, state, pier) => ({
|
|
278
|
+
pierId: pier.id,
|
|
279
|
+
axes: presentAxes(pier.axes),
|
|
280
|
+
cuts: pier.cuts.map((cut) => ({
|
|
281
|
+
z: round(cut.position, LENGTH_DECIMALS),
|
|
282
|
+
face: cut.face,
|
|
283
|
+
...cutResultants(system, state, pier.axes, cut),
|
|
284
|
+
section: cut.section,
|
|
285
|
+
})),
|
|
286
|
+
});
|
|
287
|
+
/** The presented forces of one spandrel in a raw state. */
|
|
288
|
+
export const presentSpandrelForces = (system, state, spandrel) => ({
|
|
289
|
+
spandrelId: spandrel.id,
|
|
290
|
+
axes: presentAxes(spandrel.axes),
|
|
291
|
+
length: round(spandrel.length, LENGTH_DECIMALS),
|
|
292
|
+
cuts: spandrel.cuts.map((cut) => ({
|
|
293
|
+
x: round(cut.position, LENGTH_DECIMALS),
|
|
294
|
+
face: cut.face,
|
|
295
|
+
...cutResultants(system, state, spandrel.axes, cut),
|
|
296
|
+
section: cut.section,
|
|
297
|
+
})),
|
|
298
|
+
});
|
|
299
|
+
//# sourceMappingURL=wall-forces.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"wall-forces.js","sourceRoot":"","sources":["../../../../../src/fea/analysis/post/wall-forces.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,yBAAyB,EAAE,MAAM,kCAAkC,CAAC;AAC7E,OAAO,EAAE,oBAAoB,EAAE,MAAM,oCAAoC,CAAC;AAE1E,OAAO,EAAE,WAAW,EAAE,MAAM,wCAAwC,CAAC;AACrE,OAAO,EAAE,oBAAoB,EAAE,MAAM,iBAAiB,CAAC;AAYvD,OAAO,EACL,KAAK,EACL,GAAG,EACH,gBAAgB,EAChB,GAAG,EACH,QAAQ,EACR,gBAAgB,EAChB,UAAU,EACV,aAAa,GAMd,MAAM,iBAAiB,CAAC;AAEzB,sDAAsD;AACtD,MAAM,CAAC,MAAM,eAAe,GAAkC,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;AAEjG,MAAM,cAAc,GAAG,CAAC,CAAC;AACzB,MAAM,eAAe,GAAG,CAAC,CAAC;AA6C1B,MAAM,KAAK,GAAG,CAAC,KAAa,EAAE,QAAQ,GAAG,cAAc,EAAE,EAAE,CACzD,oBAAoB,CAAC,KAAK,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,aAAa;AAC1D,MAAM,UAAU,GAAG,CAAC,CAAO,EAAQ,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAEvF,MAAM,KAAK,GAAG,CAAC,MAAgB,EAAE,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC;AAE7F;;;GAGG;AACH,MAAM,UAAU,GAAG,CACjB,MAAsB,EACtB,KAGC,EACD,EAAQ,EACR,CAAS,EACT,IAAiB,EACjB,GAAW,EACoB,EAAE;IACjC,MAAM,QAAQ,GAAG,gBAAgB,CAAC,IAAI,CAAC,CAAC;IACxC,MAAM,IAAI,GAAW,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,MAAM,EAAE,GAAG,CAAC,CAAS,EAAQ,EAAE;QAC7B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC7B,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC;IAEF,MAAM,MAAM,GAAgB,EAAE,CAAC;IAC/B,MAAM,IAAI,GAA2B,EAAE,CAAC;IACxC,KAAK,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;QACnD,MAAM,KAAK,GAAW,EAAE,CAAC;QACzB,IAAI,MAAwB,CAAC;QAC7B,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,KAAK,MAAM,KAAK,IAAI,IAAI,EAAE,CAAC;YACzB,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;YACnC,MAAM,MAAM,GAAG,KAAK,CAAC,WAAW,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;YACzC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG;gBAAE,SAAS;YACnC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG;gBAAE,SAAS;YAC7D,MAAM,OAAO,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACvF,MAAM,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC;YAC7C,KAAK,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7C,MAAM,KAAN,MAAM,GAAK,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAS,EAAC;YAChD,SAAS,GAAG,KAAK,CAAC,OAAO,CAAC,SAAS,CAAC;QACtC,CAAC;QACD,IAAI,CAAC,MAAM;YAAE,SAAS;QACtB,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,MAAM,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,IAAI,EAAE,KAAK,CAAC;YAAE,SAAS;QACvB,sEAAsE;QACtE,IAAI,KAAK,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACpD,MAAM,IAAI,GACR,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,yBAAyB;YAC5C,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC;YACd,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,yBAAyB;gBAC9C,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC;gBACd,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACrB,IAAI,IAAI;YAAE,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACpD,IAAI,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACrB,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACnB,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC;YACtB,IAAI,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC;gBAAE,KAAK,GAAG,CAAC,CAAC;YACjD,IAAI,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC;gBAAE,GAAG,GAAG,CAAC,CAAC;QAC/C,CAAC;QACD,MAAM,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,EAAE,KAAK,CAAC,CAAC;QAC3C,IAAI,MAAM,GAAG,GAAG;YAAE,IAAI,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE,SAAS,EAAE,CAAC,CAAC;IACzE,CAAC;IAED,MAAM,IAAI,GAAc,EAAE,CAAC;IAC3B,MAAM,SAAS,GAAgB,EAAE,CAAC;IAClC,MAAM,OAAO,GAA8B,EAAE,CAAC;IAC9C,KAAK,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,IAAI,KAAK,CAAC,OAAO,EAAE,CAAC;QACjD,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;YAC3B,MAAM,MAAM,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;YACjC,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;YAC/B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YACvB,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,IAAI,GAAG;gBAAE,SAAS,CAAC,wBAAwB;YAChE,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YAC1B,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC;gBACvC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;gBAC/B,MAAM,KAAK,GAAS;oBAClB,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBAC3B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBAC3B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;iBAC5B,CAAC;gBACF,SAAS,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBAC/E,OAAO,CAAC,IAAI,CAAC,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC;gBACvD,MAAM;YACR,CAAC;YACD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;gBACpD,2FAA2F;gBAC3F,MAAM,GAAG,GAAU,QAAQ,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChD,MAAM,KAAK,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;gBAClC,IAAI,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,CAAC,CAAC;gBACvC,OAAO,CAAC,IAAI,CAAC,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC;gBACvD,MAAM;YACR,CAAC;QACH,CAAC;IACH,CAAC;IAED,wEAAwE;IACxE,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,MAAM,GAAG,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5B,MAAM,WAAW,GAAG,CAAC,CAAO,EAAE,CAAS,EAAE,EAAE;QACzC,IAAI,IAAI,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,CAAC,CAAC;IACF,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE,CAAC;QACvB,MAAM,GAAG,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAS,CAAC;QAChF,WAAW,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC;IAC/C,CAAC;IACD,KAAK,MAAM,CAAC,IAAI,OAAO;QAAE,WAAW,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,IAAI,QAAc,CAAC;IACnB,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;QACb,QAAQ,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAS,CAAC;IAC9C,CAAC;SAAM,CAAC;QACN,2DAA2D;QAC3D,MAAM,MAAM,GAAG;YACb,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzF,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;YAC3B,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;SACjC,CAAC;QACF,QAAQ,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CACtB,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CACrE,CAAC;IACZ,CAAC;IAED,OAAO;QACL,IAAI;QACJ,QAAQ;QACR,OAAO,EAAE;YACP,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;gBACvB,MAAM,EAAE,GAAG,CAAC,MAAM;gBAClB,KAAK,EAAE,UAAU,CAAC,GAAG,CAAC,KAAK,CAAC;gBAC5B,GAAG,EAAE,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC;gBACxB,MAAM,EAAE,KAAK,CAAC,GAAG,CAAC,MAAM,EAAE,eAAe,CAAC;gBAC1C,SAAS,EAAE,GAAG,CAAC,SAAS;aACzB,CAAC,CAAC;YACH,OAAO,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,KAAK,EAAE,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC3F,MAAM,EAAE,KAAK,CACX,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,EAC9C,eAAe,CAChB;YACD,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,eAAe,CAAC;YAClC,QAAQ,EAAE,UAAU,CAAC,QAAQ,CAAC;SAC/B;QACD,MAAM;QACN,IAAI;QACJ,SAAS;KACV,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAC1B,KAAyD,EACzD,MAAsB,EACgC,EAAE;IACxD,MAAM,EAAE,GAAG,CAAC,CAAS,EAAQ,EAAE;QAC7B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC7B,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC;IACF,MAAM,OAAO,GAAG,CAAC,IAAc,EAAE,IAA6B,EAAgB,EAAE;QAC9E,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,OAAO,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;YAClD,MAAM;YACN,MAAM,EAAE,MAAM,CAAC,YAAY,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,EAAE;SAC9C,CAAC,CAAC,CAAC;QACJ,MAAM,OAAO,GAAG,CAAC,IAAI,CAAC,SAAS,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC;YACxD,QAAQ;YACR,MAAM,EAAE,MAAM,CAAC,YAAY,CAAC,MAAM,CAAC,WAAW,CAAC,GAAG,CAAC,QAAQ,CAAW,CAAC,CAAC,MAAM;SAC/E,CAAC,CAAC,CAAC;QACJ,MAAM,KAAK,GAAc;YACvB,KAAK,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,CAAC;gBACxC,MAAM;gBACN,QAAQ,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;aAClE,CAAC,CAAC;YACH,OAAO,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,CAAC;gBAC9C,QAAQ;gBACR,MAAM,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAgB,EAAE,CAAC;oBACtC,EAAE,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC;oBAChC,EAAE,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;iBAC/B,CAAC;aACH,CAAC,CAAC;SACJ,CAAC;QACF,MAAM,IAAI,GAAG,QAAQ,CAAC,IAAI,EAAE,KAAK,EAAE;YACjC,KAAK,EAAG,IAAiB,CAAC,KAAK;YAC/B,OAAO,EAAG,IAAqB,CAAC,OAAO;SACxC,CAAC,CAAC;QACH,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,MAAM,GAAG,UAAU,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACrC,MAAM,GAAG,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;QAClC,MAAM,MAAM,GAAG,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;QAChD,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;QACvD,MAAM,IAAI,GAAG,gBAAgB,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,CAAC;YAC1E,QAAQ,EAAE,CAAC,GAAG,MAAM;YACpB,GAAG,UAAU,CAAC,MAAM,EAAE,EAAE,KAAK,EAAE,OAAO,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,CAAC;SAC5D,CAAC,CAAC,CAAC;QACJ,OAAO,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,IAAI,EAAE,CAAC;IACtE,CAAC,CAAC;IACF,OAAO;QACL,KAAK,EAAE,CAAC,KAAK,CAAC,KAAK,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;QACzD,SAAS,EAAE,CAAC,KAAK,CAAC,SAAS,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;KACtE,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,YAAY,GAAG,CAAC,MAAsB,EAAE,KAAe,EAAE,GAAgB,EAAE,EAAE;IACjF,MAAM,CAAC,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC;IACvB,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,CAAO,EAAE,IAAY,EAAE,EAAE;QACxD,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACpB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACpB,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAChC,CAAC;IACH,CAAC,CAAC;IAEF,KAAK,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,IAAI,GAAG,CAAC,MAAM,EAAE,CAAC;QACzD,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,WAAW,CAAC,KAAK,CAAC,WAAW,CAAC,CAAC;QAC5C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,KAAK,CAAC,cAAc,CAAC,KAAK,CAAC,CAAC;QACtC,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,GAAG,GAAG,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC7C,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;oBAAE,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChD,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACX,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;YACvE,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC;YAChD,KAAK,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;QAChF,CAAC;IACH,CAAC;IAED,MAAM,QAAQ,GAAG,CAAC,MAAkB,EAAE,KAAe,EAAQ,EAAE,CAAC;QAC9D,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QAC3E,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QAC3E,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;KAC5E,CAAC;IACF,KAAK,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC;QACnD,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,KAAK,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QAChE,KAAK,CAAC,KAAK,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACvF,CAAC;IACD,KAAK,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,GAAG,CAAC,SAAS,EAAE,CAAC;QACtD,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,oBAAoB,CAC5B,KAAK,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACrC,KAAK,CAAC,WAAW,CAAC,KAAK,CAAC,EACxB,CAAC,CACF,CAAC;QACF,KAAK,CAAC,KAAK,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACvF,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAClB,CAAC,CAAC;AAEF,wEAAwE;AACxE,MAAM,aAAa,GAAG,CACpB,MAAsB,EACtB,KAAe,EACf,IAAc,EACd,GAAgB,EACA,EAAE;IAClB,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,YAAY,CAAC,MAAM,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;IAClD,OAAO;QACL,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAChC,EAAE,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACjC,EAAE,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACjC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAChC,EAAE,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACjC,EAAE,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;KAClC,CAAC;AACJ,CAAC,CAAC;AAEF,MAAM,WAAW,GAAG,CAAC,IAAc,EAAsB,EAAE,CACzD,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAS,CAAuB,CAAC;AAE3E,uDAAuD;AACvD,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAC/B,MAAsB,EACtB,KAAe,EACf,IAAkB,EACN,EAAE,CAAC,CAAC;IAChB,MAAM,EAAE,IAAI,CAAC,EAAE;IACf,IAAI,EAAE,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC;IAC5B,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAC5B,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,QAAQ,EAAE,eAAe,CAAC;QACvC,IAAI,EAAE,GAAG,CAAC,IAA0C;QACpD,GAAG,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,IAAI,EAAE,GAAG,CAAC;QAC/C,OAAO,EAAE,GAAG,CAAC,OAAO;KACrB,CAAC,CAAC;CACJ,CAAC,CAAC;AAEH,2DAA2D;AAC3D,MAAM,CAAC,MAAM,qBAAqB,GAAG,CACnC,MAAsB,EACtB,KAAe,EACf,QAAsB,EACN,EAAE,CAAC,CAAC;IACpB,UAAU,EAAE,QAAQ,CAAC,EAAE;IACvB,IAAI,EAAE,WAAW,CAAC,QAAQ,CAAC,IAAI,CAAC;IAChC,MAAM,EAAE,KAAK,CAAC,QAAQ,CAAC,MAAM,EAAE,eAAe,CAAC;IAC/C,IAAI,EAAE,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAChC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,QAAQ,EAAE,eAAe,CAAC;QACvC,IAAI,EAAE,GAAG,CAAC,IAA8C;QACxD,GAAG,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,QAAQ,CAAC,IAAI,EAAE,GAAG,CAAC;QACnD,OAAO,EAAE,GAAG,CAAC,OAAO;KACrB,CAAC,CAAC;CACJ,CAAC,CAAC"}
|
|
@@ -0,0 +1,173 @@
|
|
|
1
|
+
import { VERTICAL_MEMBER_TOLERANCE, WALL_CUT_TOLERANCE } from "../../../core/constants/fea.constant.js";
|
|
2
|
+
export const sub = (p, q) => [p[0] - q[0], p[1] - q[1], p[2] - q[2]];
|
|
3
|
+
export const dot = (p, q) => p[0] * q[0] + p[1] * q[1] + p[2] * q[2];
|
|
4
|
+
export const cross = (p, q) => [
|
|
5
|
+
p[1] * q[2] - p[2] * q[1],
|
|
6
|
+
p[2] * q[0] - p[0] * q[2],
|
|
7
|
+
p[0] * q[1] - p[1] * q[0],
|
|
8
|
+
];
|
|
9
|
+
const norm = (p) => Math.hypot(p[0], p[1], p[2]);
|
|
10
|
+
const unit = (p) => {
|
|
11
|
+
const n = norm(p);
|
|
12
|
+
return n > 0 ? [p[0] / n, p[1] / n, p[2] / n] : [0, 0, 0];
|
|
13
|
+
};
|
|
14
|
+
/** Unit normal of a polygon by Newell's method (zero for a degenerate polygon). */
|
|
15
|
+
export const polygonNormal = (points) => {
|
|
16
|
+
const n = [0, 0, 0];
|
|
17
|
+
points.forEach((p, i) => {
|
|
18
|
+
const q = points[(i + 1) % points.length];
|
|
19
|
+
n[0] += (p[1] - q[1]) * (p[2] + q[2]);
|
|
20
|
+
n[1] += (p[2] - q[2]) * (p[0] + q[0]);
|
|
21
|
+
n[2] += (p[0] - q[0]) * (p[1] + q[1]);
|
|
22
|
+
});
|
|
23
|
+
return unit(n);
|
|
24
|
+
};
|
|
25
|
+
/** The horizontal direction in a vertical plane of normal `n`, toward +X (+Y if along Y). */
|
|
26
|
+
export const horizontalInPlane = (n) => {
|
|
27
|
+
const d = unit([-n[1], n[0], 0]);
|
|
28
|
+
const flip = Math.abs(d[0]) > VERTICAL_MEMBER_TOLERANCE ? d[0] < 0 : d[1] < 0;
|
|
29
|
+
return flip ? [-d[0], -d[1], -d[2]] : d;
|
|
30
|
+
};
|
|
31
|
+
/**
|
|
32
|
+
* Axes of a pier or spandrel (see `WallPier`, `WallSpandrel`).
|
|
33
|
+
*
|
|
34
|
+
* - Pier: 1 = +Z; 2 = global X rotated by `angle` about +Z, else the horizontal direction of
|
|
35
|
+
* the area with the largest horizontal extent (the first wins a tie), else global X; 3 = 1 × 2.
|
|
36
|
+
* - Spandrel: 1 = horizontal in the plane of its first area, toward +X (+Y for a plane parallel
|
|
37
|
+
* to Y), reversed by `reverse`; 2 = +Z; 3 = 1 × 2.
|
|
38
|
+
*/
|
|
39
|
+
export const wallAxes = (kind, shape, options) => {
|
|
40
|
+
const up = [0, 0, 1];
|
|
41
|
+
if (kind === "pier") {
|
|
42
|
+
let e2 = [1, 0, 0];
|
|
43
|
+
if (options.angle !== undefined) {
|
|
44
|
+
e2 = [Math.cos(options.angle), Math.sin(options.angle), 0];
|
|
45
|
+
}
|
|
46
|
+
else {
|
|
47
|
+
let best = 0;
|
|
48
|
+
for (const area of shape.areas) {
|
|
49
|
+
if (area.elements.length === 0)
|
|
50
|
+
continue;
|
|
51
|
+
const d = horizontalInPlane(polygonNormal(area.elements[0]));
|
|
52
|
+
const s = area.elements.flat().map((p) => dot(p, d));
|
|
53
|
+
const extent = Math.max(...s) - Math.min(...s);
|
|
54
|
+
if (extent > best) {
|
|
55
|
+
best = extent;
|
|
56
|
+
e2 = d;
|
|
57
|
+
}
|
|
58
|
+
}
|
|
59
|
+
}
|
|
60
|
+
return [up, e2, cross(up, e2)];
|
|
61
|
+
}
|
|
62
|
+
const first = shape.areas.find((a) => a.elements.length > 0);
|
|
63
|
+
let e1 = first ? horizontalInPlane(polygonNormal(first.elements[0])) : [1, 0, 0];
|
|
64
|
+
if (options.reverse)
|
|
65
|
+
e1 = [-e1[0], -e1[1], -e1[2]];
|
|
66
|
+
return [e1, up, cross(e1, up)];
|
|
67
|
+
};
|
|
68
|
+
/** Smallest and largest coordinate along `e1` of every element corner and member end. */
|
|
69
|
+
export const wallExtent = (shape, e1) => {
|
|
70
|
+
const s = [
|
|
71
|
+
...shape.areas.flatMap((a) => a.elements.flat()),
|
|
72
|
+
...shape.members.flatMap((m) => m.pieces.flat()),
|
|
73
|
+
].map((p) => dot(p, e1));
|
|
74
|
+
return { min: Math.min(...s), max: Math.max(...s) };
|
|
75
|
+
};
|
|
76
|
+
/** The distance within which a point counts as on a cut (see `WALL_CUT_TOLERANCE`). */
|
|
77
|
+
export const wallTolerance = (extent) => WALL_CUT_TOLERANCE * Math.max(extent.max - extent.min, 1);
|
|
78
|
+
const range = (points, e1) => {
|
|
79
|
+
const s = points.map((p) => dot(p, e1));
|
|
80
|
+
return { min: Math.min(...s), max: Math.max(...s) };
|
|
81
|
+
};
|
|
82
|
+
/** Whether anything of the shape lies on the negative / positive side touching the cut at s. */
|
|
83
|
+
const touching = (shape, e1, s, tol) => {
|
|
84
|
+
let negative = false;
|
|
85
|
+
let positive = false;
|
|
86
|
+
const check = (points) => {
|
|
87
|
+
const r = range(points, e1);
|
|
88
|
+
if (r.max - r.min <= tol)
|
|
89
|
+
return; // lies in the cut plane
|
|
90
|
+
if (Math.abs(r.max - s) <= tol)
|
|
91
|
+
negative = true;
|
|
92
|
+
if (Math.abs(r.min - s) <= tol)
|
|
93
|
+
positive = true;
|
|
94
|
+
};
|
|
95
|
+
shape.areas.forEach((a) => a.elements.forEach(check));
|
|
96
|
+
shape.members.forEach((m) => m.pieces.forEach(check));
|
|
97
|
+
return { negative, positive };
|
|
98
|
+
};
|
|
99
|
+
/**
|
|
100
|
+
* The cuts of a pier or spandrel, in order along axis 1.
|
|
101
|
+
*
|
|
102
|
+
* - Pier: at each elevation where one of its areas or members starts or ends, a `"top"` cut
|
|
103
|
+
* when something of the pier ends there from below and a `"bottom"` cut when something starts
|
|
104
|
+
* there going up (the top first); then each extra cut not at such an elevation, as `"cut"`.
|
|
105
|
+
* - Spandrel: `"left"` at its smallest coordinate, `"right"` at its largest, and each extra cut
|
|
106
|
+
* strictly between them as `"cut"`.
|
|
107
|
+
*
|
|
108
|
+
* @param extra - The extra cuts as coordinates along axis 1 (mm)
|
|
109
|
+
*/
|
|
110
|
+
export const wallCutPositions = (kind, shape, e1, extra = []) => {
|
|
111
|
+
const extent = wallExtent(shape, e1);
|
|
112
|
+
const tol = wallTolerance(extent);
|
|
113
|
+
const auto = [];
|
|
114
|
+
const addAuto = (s) => {
|
|
115
|
+
if (!auto.some((a) => Math.abs(a - s) <= tol))
|
|
116
|
+
auto.push(s);
|
|
117
|
+
};
|
|
118
|
+
if (kind === "pier") {
|
|
119
|
+
for (const area of shape.areas) {
|
|
120
|
+
if (area.elements.length === 0)
|
|
121
|
+
continue;
|
|
122
|
+
const r = range(area.elements.flat(), e1);
|
|
123
|
+
addAuto(r.min);
|
|
124
|
+
addAuto(r.max);
|
|
125
|
+
}
|
|
126
|
+
for (const member of shape.members) {
|
|
127
|
+
const r = range(member.pieces.flat(), e1);
|
|
128
|
+
addAuto(r.min);
|
|
129
|
+
addAuto(r.max);
|
|
130
|
+
}
|
|
131
|
+
}
|
|
132
|
+
else {
|
|
133
|
+
addAuto(extent.min);
|
|
134
|
+
addAuto(extent.max);
|
|
135
|
+
}
|
|
136
|
+
auto.sort((a, b) => a - b);
|
|
137
|
+
const cuts = [];
|
|
138
|
+
for (const s of auto) {
|
|
139
|
+
if (kind === "spandrel") {
|
|
140
|
+
cuts.push({ s, face: s === extent.min ? "left" : "right" });
|
|
141
|
+
continue;
|
|
142
|
+
}
|
|
143
|
+
const side = touching(shape, e1, s, tol);
|
|
144
|
+
if (side.negative)
|
|
145
|
+
cuts.push({ s, face: "top" });
|
|
146
|
+
if (side.positive)
|
|
147
|
+
cuts.push({ s, face: "bottom" });
|
|
148
|
+
}
|
|
149
|
+
extra.forEach((s, cut) => {
|
|
150
|
+
if (auto.some((a) => Math.abs(a - s) <= tol))
|
|
151
|
+
return;
|
|
152
|
+
if (cuts.some((c) => c.face === "cut" && Math.abs(c.s - s) <= tol))
|
|
153
|
+
return;
|
|
154
|
+
cuts.push({ s, face: "cut", cut });
|
|
155
|
+
});
|
|
156
|
+
// stable: a "top" stays before the "bottom" at the same level
|
|
157
|
+
return cuts.sort((a, b) => a.s - b.s);
|
|
158
|
+
};
|
|
159
|
+
/** Whether a face takes the elements on the negative side of its cut (below, or to the left). */
|
|
160
|
+
export const fromNegativeSide = (face) => face === "top" || face === "right" || face === "cut";
|
|
161
|
+
/** Areas with a shell element on both sides of the cut at s (no mesh line there). */
|
|
162
|
+
export const areasAcrossCut = (shape, e1, s, tol) => shape.areas
|
|
163
|
+
.filter((a) => a.elements.some((el) => {
|
|
164
|
+
const r = range(el, e1);
|
|
165
|
+
return r.min < s - tol && r.max > s + tol;
|
|
166
|
+
}))
|
|
167
|
+
.map((a) => a.areaId);
|
|
168
|
+
/** Whether some element or member of the shape reaches the cut at s. */
|
|
169
|
+
export const reachesCut = (shape, e1, s, tol) => [...shape.areas.flatMap((a) => a.elements), ...shape.members.flatMap((m) => m.pieces)].some((points) => {
|
|
170
|
+
const r = range(points, e1);
|
|
171
|
+
return r.min <= s + tol && r.max >= s - tol;
|
|
172
|
+
});
|
|
173
|
+
//# sourceMappingURL=wall-geometry.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"wall-geometry.js","sourceRoot":"","sources":["../../../../../src/fea/analysis/post/wall-geometry.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,yBAAyB,EAAE,kBAAkB,EAAE,MAAM,kCAAkC,CAAC;AA6BjG,MAAM,CAAC,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;AACvF,MAAM,CAAC,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;AACjF,MAAM,CAAC,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC;IAC/C,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;CAC1B,CAAC;AACF,MAAM,IAAI,GAAG,CAAC,CAAO,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvD,MAAM,IAAI,GAAG,CAAC,CAAO,EAAQ,EAAE;IAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IAClB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAC5D,CAAC,CAAC;AAEF,mFAAmF;AACnF,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,MAAuB,EAAQ,EAAE;IAC7D,MAAM,CAAC,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACtB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;QAC1C,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxC,CAAC,CAAC,CAAC;IACH,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC;AACjB,CAAC,CAAC;AAEF,6FAA6F;AAC7F,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,CAAO,EAAQ,EAAE;IACjD,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,yBAAyB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC9E,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1C,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,QAAQ,GAAG,CACtB,IAAc,EACd,KAAgB,EAChB,OAA8C,EACpC,EAAE;IACZ,MAAM,EAAE,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;QACpB,IAAI,EAAE,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACzB,IAAI,OAAO,CAAC,KAAK,KAAK,SAAS,EAAE,CAAC;YAChC,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7D,CAAC;aAAM,CAAC;YACN,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,KAAK,MAAM,IAAI,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;gBAC/B,IAAI,IAAI,CAAC,QAAQ,CAAC,MAAM,KAAK,CAAC;oBAAE,SAAS;gBACzC,MAAM,CAAC,GAAG,iBAAiB,CAAC,aAAa,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC7D,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBACrD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC/C,IAAI,MAAM,GAAG,IAAI,EAAE,CAAC;oBAClB,IAAI,GAAG,MAAM,CAAC;oBACd,EAAE,GAAG,CAAC,CAAC;gBACT,CAAC;YACH,CAAC;QACH,CAAC;QACD,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IACjC,CAAC;IACD,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC7D,IAAI,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,iBAAiB,CAAC,aAAa,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAU,CAAC;IAC3F,IAAI,OAAO,CAAC,OAAO;QAAE,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;AACjC,CAAC,CAAC;AAEF,yFAAyF;AACzF,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAgC,EAAE;IACrF,MAAM,CAAC,GAAG;QACR,GAAG,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,EAAE,CAAC;QAChD,GAAG,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC;KACjD,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACzB,OAAO,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;AACtD,CAAC,CAAC;AAEF,uFAAuF;AACvF,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,MAAoC,EAAE,EAAE,CACpE,kBAAkB,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;AAE5D,MAAM,KAAK,GAAG,CAAC,MAAuB,EAAE,EAAQ,EAAE,EAAE;IAClD,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACxC,OAAO,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;AACtD,CAAC,CAAC;AAEF,gGAAgG;AAChG,MAAM,QAAQ,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAE,CAAS,EAAE,GAAW,EAAE,EAAE;IACtE,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,MAAM,KAAK,GAAG,CAAC,MAAuB,EAAE,EAAE;QACxC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;QAC5B,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG;YAAE,OAAO,CAAC,wBAAwB;QAC1D,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,GAAG;YAAE,QAAQ,GAAG,IAAI,CAAC;QAChD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,GAAG;YAAE,QAAQ,GAAG,IAAI,CAAC;IAClD,CAAC,CAAC;IACF,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;IACtD,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;IACtD,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,CAAC;AAChC,CAAC,CAAC;AAEF;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAC9B,IAAc,EACd,KAAgB,EAChB,EAAQ,EACR,QAA2B,EAAE,EACV,EAAE;IACrB,MAAM,MAAM,GAAG,UAAU,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IAClC,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,MAAM,OAAO,GAAG,CAAC,CAAS,EAAE,EAAE;QAC5B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9D,CAAC,CAAC;IACF,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;QACpB,KAAK,MAAM,IAAI,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;YAC/B,IAAI,IAAI,CAAC,QAAQ,CAAC,MAAM,KAAK,CAAC;gBAAE,SAAS;YACzC,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,CAAC;YAC1C,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACf,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,CAAC;QACD,KAAK,MAAM,MAAM,IAAI,KAAK,CAAC,OAAO,EAAE,CAAC;YACnC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,CAAC;YAC1C,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACf,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,CAAC;IACH,CAAC;SAAM,CAAC;QACN,OAAO,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;QACpB,OAAO,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACtB,CAAC;IACD,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAE3B,MAAM,IAAI,GAAsB,EAAE,CAAC;IACnC,KAAK,MAAM,CAAC,IAAI,IAAI,EAAE,CAAC;QACrB,IAAI,IAAI,KAAK,UAAU,EAAE,CAAC;YACxB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,KAAK,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;YAC5D,SAAS;QACX,CAAC;QACD,MAAM,IAAI,GAAG,QAAQ,CAAC,KAAK,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QACzC,IAAI,IAAI,CAAC,QAAQ;YAAE,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,CAAC;QACjD,IAAI,IAAI,CAAC,QAAQ;YAAE,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,CAAC,CAAC;IACtD,CAAC;IACD,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE;QACvB,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAAE,OAAO;QACrD,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAAE,OAAO;QAC3E,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC;IACrC,CAAC,CAAC,CAAC;IACH,8DAA8D;IAC9D,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AACxC,CAAC,CAAC;AAEF,iGAAiG;AACjG,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAAC,IAAiB,EAAE,EAAE,CACpD,IAAI,KAAK,KAAK,IAAI,IAAI,KAAK,OAAO,IAAI,IAAI,KAAK,KAAK,CAAC;AAEvD,qFAAqF;AACrF,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAE,CAAS,EAAE,GAAW,EAAW,EAAE,CAC5F,KAAK,CAAC,KAAK;KACR,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CACZ,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE;IACrB,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;AAC5C,CAAC,CAAC,CACH;KACA,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;AAE1B,wEAAwE;AACxE,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAE,CAAS,EAAE,GAAW,EAAW,EAAE,CACxF,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,EAAE,GAAG,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CACzF,CAAC,MAAM,EAAE,EAAE;IACT,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IAC5B,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC;AAC9C,CAAC,CACF,CAAC"}
|
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
import { SINGULARITY_TOLERANCE } from "../../../core/constants/
|
|
1
|
+
import { SINGULARITY_TOLERANCE } from "../../../core/constants/fea.constant.js";
|
|
2
2
|
import { FeaErrors } from "../../../core/errors/fea.error.js";
|
|
3
3
|
import { memberGlobalStiffness } from "../../frame/solver/member-assembly.js";
|
|
4
4
|
import { prepareSystemMembers, } from "../../frame/solver/system-member.js";
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"analysis-system.js","sourceRoot":"","sources":["../../../../../src/fea/analysis/solver/analysis-system.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,qBAAqB,EAAE,MAAM,
|
|
1
|
+
{"version":3,"file":"analysis-system.js","sourceRoot":"","sources":["../../../../../src/fea/analysis/solver/analysis-system.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,qBAAqB,EAAE,MAAM,kCAAkC,CAAC;AACzE,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,uCAAuC,CAAC;AAC9E,OAAO,EACL,oBAAoB,GAGrB,MAAM,qCAAqC,CAAC;AAC7C,OAAO,EACL,kBAAkB,EAClB,0BAA0B,GAE3B,MAAM,uCAAuC,CAAC;AAC/C,OAAO,EACL,uBAAuB,GAExB,MAAM,4CAA4C,CAAC;AACpD,OAAO,EAAE,WAAW,EAAyB,MAAM,wCAAwC,CAAC;AAE5F,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,mCAAmC,CAAC;AAC5E,OAAO,EAAE,mBAAmB,EAAoB,MAAM,oCAAoC,CAAC;AAE3F,OAAO,EACL,aAAa,EACb,mBAAmB,EACnB,UAAU,EACV,aAAa,EACb,WAAW,GAEZ,MAAM,oBAAoB,CAAC;AAI5B,MAAM,SAAS,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AA8BhE;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,qBAAqB,GAAG,CAAC,KAAe,EAAkB,EAAE;IACvE,MAAM,EAAE,KAAK,EAAE,OAAO,GAAG,EAAE,EAAE,UAAU,EAAE,GAAG,KAAK,CAAC;IAClD,iEAAiE;IACjE,MAAM,MAAM,GAAG,CAAC,KAAK,CAAC,MAAM,IAAI,EAAE,CAAuC,CAAC;IAC1E,MAAM,SAAS,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACzE,MAAM,WAAW,GAAG,IAAI,GAAG,CAAgB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7E,MAAM,UAAU,GAAG,IAAI,GAAG,CAAgB,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7E,MAAM,QAAQ,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IAClC,MAAM,QAAQ,GAAG,uBAAuB,CAAC,UAAU,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;IAE1E,MAAM,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,YAAY,EAAE,GAAG,oBAAoB,CAC1E,OAAO,EACP,KAAK,EACL,SAAS,CACV,CAAC;IACF,MAAM,YAAY,GAAG,mBAAmB,CAAC,MAAM,EAAE,KAAK,EAAE,SAAS,CAAC,CAAC;IAEnE,MAAM,YAAY,GAAG,IAAI,GAAG,EAAmB,CAAC;IAChD,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE;QAC7B,MAAM,IAAI,GAAG,YAAY,CAAC,GAAG,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;QACzC,IAAI,IAAI;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;;YAClB,YAAY,CAAC,GAAG,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACxC,CAAC,CAAC,CAAC;IAEH,yFAAyF;IACzF,MAAM,cAAc,GAAuB,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QACnE,KAAK,EAAE,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,CAAC,QAAQ,CAAC;QACjC,OAAO,EAAE,qBAAqB,CAAC,CAAC,CAAC;KAClC,CAAC,CAAC,CAAC;IACJ,MAAM,aAAa,GAAuB,YAAY,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAClE,KAAK,EAAE,EAAE,CAAC,WAAW;QACrB,OAAO,EAAE,EAAE,CAAC,OAAO,CAAC,OAAO;KAC5B,CAAC,CAAC,CAAC;IACJ,MAAM,QAAQ,GAAG,CAAC,GAAG,cAAc,EAAE,GAAG,aAAa,CAAC,CAAC;IAEvD,yFAAyF;IACzF,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IACtF,MAAM,cAAc,GAAG,IAAI,GAAG,CAC5B,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,WAAW,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAChF,CAAC;IACF,MAAM,UAAU,GAAG,0BAA0B,CAC3C,CAAC,GAAG,cAAc,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,wBAAwB;IACnE,aAAa,EACb,QAAQ,CACT,CAAC;IACF,IAAI,UAAU,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC1B,WAAW,CACT,UAAU,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE;YAChC,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CACtB,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAC/D,CAAC,CACF,CAAC;YACF,OAAO,UAAU,CAAC,SAAS,CAAC,qBAAqB,EAAE;gBACjD,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE;gBACtB,OAAO,EAAE,YAAY,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,WAAW,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE;gBACrE,GAAG,EAAE,SAAS,CAAC,IAAI,CAAC;aACrB,CAAC,CAAC;QACL,CAAC,CAAC,CACH,CAAC;IACJ,CAAC;IAED,qFAAqF;IACrF,MAAM,SAAS,GAAe,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAClD,KAAK,MAAM,OAAO,IAAI,cAAc,EAAE,CAAC;QACrC,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK,EAAE,CAAC;YAC9B,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK;gBAAE,IAAI,CAAC,KAAK,CAAC;oBAAE,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnE,CAAC;IACH,CAAC;IACD,uFAAuF;IACvF,IAAI,aAAa,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC7B,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC;QACrD,KAAK,MAAM,OAAO,IAAI,aAAa,EAAE,CAAC;YACpC,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK,EAAE,CAAC;gBAC9B,KAAK,MAAM,CAAC,IAAI,OAAO,CAAC,KAAK,EAAE,CAAC;oBAC9B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC;wBAChC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBAChB,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;oBACvB,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,UAAU,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,MAAM,IAAI,IAAI,mBAAmB,CAAC,SAAS,CAAC,EAAE,CAAC;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YACzB,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;gBACzB,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;gBAC5B,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACjB,CAAC;QACH,CAAC;IACH,CAAC;IAED,MAAM,cAAc,GAAG,kBAAkB,CAAC,KAAK,CAAC,MAAM,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;IAE5E,qFAAqF;IACrF,IAAI,MAAM,GAAyB,IAAI,CAAC;IACxC,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACpC,MAAM,gBAAgB,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,WAAW,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7F,KAAK,MAAM,GAAG,IAAI,gBAAgB,EAAE,CAAC;YACnC,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;YAChC,KAAK,MAAM,CAAC,IAAI,MAAM;gBAAE,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBAAE,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAC7D,CAAC;QAED,MAAM,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC;QAC/B,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE;YACtC,yFAAyF;YACzF,MAAM,CAAC,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC,OAAO,CAAC;YAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACpC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,SAAS;gBACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBACpC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;wBAAE,UAAU,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/D,CAAC;YACH,CAAC;QACH,CAAC,CAAC,CAAC;QACH,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE;YACtB,IAAI,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC;gBAAE,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1E,CAAC,CAAC,CAAC;QACH,6FAA6F;QAC7F,KAAK,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,cAAc,EAAE,CAAC;YACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBACtC,IAAI,EAAE,GAAG,CAAC;oBAAE,SAAS;gBACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;oBACtC,IAAI,EAAE,IAAI,EAAE;wBAAE,UAAU,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrE,CAAC;YACH,CAAC;QACH,CAAC;QAED,MAAM,MAAM,GAAG,aAAa,CAAC,CAAC,EAAE,qBAAqB,CAAC,CAAC;QACvD,IAAI,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC;YACf,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC;YACrC,WAAW,CAAC;gBACV,UAAU,CAAC,SAAS,CAAC,kBAAkB,EAAE;oBACvC,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE;oBACrC,GAAG,EAAE,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC;iBACxB,CAAC;aACH,CAAC,CAAC;QACL,CAAC;QACD,MAAM,GAAG,CAAC,CAAC;IACb,CAAC;IAED,OAAO;QACL,KAAK;QACL,SAAS;QACT,WAAW;QACX,YAAY;QACZ,OAAO,EAAE,aAAa;QACtB,UAAU;QACV,YAAY;QACZ,MAAM,EAAE,YAAY;QACpB,QAAQ;QACR,QAAQ;QACR,QAAQ;QACR,IAAI;QACJ,cAAc;QACd,MAAM;QACN,WAAW,EAAE,MAAM,CAAC,CAAC,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;KAC9C,CAAC;AACJ,CAAC,CAAC"}
|
|
@@ -1,7 +1,9 @@
|
|
|
1
1
|
import { areaMemberLoads } from "../../frame/area/area-member-loads.js";
|
|
2
|
+
import { prepareStrips } from "../../shell/post/strip-forces.js";
|
|
2
3
|
import { throwIssues } from "../../shared/validation/issues.js";
|
|
3
4
|
import { roundToDecimalPlaces } from "../../../utils/math.js";
|
|
4
5
|
import { buildEnvelope } from "../post/envelope.js";
|
|
6
|
+
import { prepareWalls } from "../post/wall-forces.js";
|
|
5
7
|
import { presentMemberDiagram, presentState } from "../post/present-state.js";
|
|
6
8
|
import { expandModel } from "../mesh/expand-model.js";
|
|
7
9
|
import { validateModel } from "../validation/validate-model.js";
|
|
@@ -40,7 +42,7 @@ const groupByCase = (loads = []) => {
|
|
|
40
42
|
* case, null areas, and combinations
|
|
41
43
|
* @param options - Diagram station count and detail
|
|
42
44
|
* @returns `{ result, memberDiagram }`
|
|
43
|
-
* @throws {FeaError} 301–317, 326 and
|
|
45
|
+
* @throws {FeaError} 301–317, 326, 327 and 329–331 for an invalid model (with `details` and
|
|
44
46
|
* every problem in `issues`), 328 for a translation nothing resists at a node connected only
|
|
45
47
|
* to shells, 304 for an unstable structure
|
|
46
48
|
*/
|
|
@@ -54,6 +56,8 @@ export const createAnalysis = (model, options = {}) => {
|
|
|
54
56
|
const resolved = resolveCombinations(combinations, new Set(model.loadCases.map((c) => c.id)));
|
|
55
57
|
const withStations = (options.diagrams ?? "full") === "full";
|
|
56
58
|
const fromAreas = areaMemberLoads(model);
|
|
59
|
+
const strips = prepareStrips(model.strips ?? [], system.shells);
|
|
60
|
+
const walls = prepareWalls(model, system);
|
|
57
61
|
const states = new Map();
|
|
58
62
|
const cases = {};
|
|
59
63
|
const nodalByCase = groupByCase(model.nodalLoads);
|
|
@@ -70,7 +74,7 @@ export const createAnalysis = (model, options = {}) => {
|
|
|
70
74
|
shellLoads: shellByCase.get(loadCase.id) ?? [],
|
|
71
75
|
});
|
|
72
76
|
states.set(loadCase.id, state);
|
|
73
|
-
cases[loadCase.id] = presentState(system, state, { id: loadCase.id, kind: "case", name: loadCase.name, type: loadCase.type }, { stations: options.stations, withStations });
|
|
77
|
+
cases[loadCase.id] = presentState(system, state, { id: loadCase.id, kind: "case", name: loadCase.name, type: loadCase.type }, { stations: options.stations, withStations }, strips.strips, walls);
|
|
74
78
|
}
|
|
75
79
|
const combinationStates = new Map();
|
|
76
80
|
const combinationResults = {};
|
|
@@ -79,11 +83,12 @@ export const createAnalysis = (model, options = {}) => {
|
|
|
79
83
|
const terms = [...resolved.get(combination.id)].map(([caseId, factor]) => ({ state: states.get(caseId), factor }));
|
|
80
84
|
const state = combineStates(terms);
|
|
81
85
|
combinationStates.set(combination.id, state);
|
|
82
|
-
combinationResults[combination.id] = presentState(system, state, { id: combination.id, kind: "combination", name: combination.name }, { stations: options.stations, withStations });
|
|
86
|
+
combinationResults[combination.id] = presentState(system, state, { id: combination.id, kind: "combination", name: combination.name }, { stations: options.stations, withStations }, strips.strips, walls);
|
|
83
87
|
}
|
|
84
88
|
const warnings = [
|
|
85
89
|
...fromAreas.warnings,
|
|
86
90
|
...expanded.warnings,
|
|
91
|
+
...strips.warnings,
|
|
87
92
|
...[...Object.values(cases), ...Object.values(combinationResults)].flatMap((r) => r.warnings.map((w) => ({
|
|
88
93
|
...w,
|
|
89
94
|
message: `${r.kind === "case" ? "Load case" : "Load combination"} "${r.id}": ${w.message}`,
|
|
@@ -157,7 +162,10 @@ export const createAnalysis = (model, options = {}) => {
|
|
|
157
162
|
* Area loads on null areas (`areas`, `areaLoads`) are first turned into member loads on the
|
|
158
163
|
* members along each area's edges (see `distributeAreaLoads`); the areas add no stiffness.
|
|
159
164
|
* Shell areas (`shells`) are four-node flat shell elements with their own stiffness; area loads
|
|
160
|
-
* and self-weight on them act as forces at their corners.
|
|
165
|
+
* and self-weight on them act as forces at their corners. Strips (`strips`) integrate the shell
|
|
166
|
+
* forces across cuts at stations along them (`stripForces`), for strip design of slabs.
|
|
167
|
+
* Piers and spandrels (`piers`, `spandrels`) sum the nodal forces of their shell elements and
|
|
168
|
+
* members on cuts across them (`pierForces`, `spandrelForces`), for wall design.
|
|
161
169
|
*
|
|
162
170
|
* Model and result are plain data (JSON-serialisable), so the analysis can run in a Web Worker.
|
|
163
171
|
*
|
|
@@ -166,7 +174,7 @@ export const createAnalysis = (model, options = {}) => {
|
|
|
166
174
|
* @param options - Diagram station count and detail
|
|
167
175
|
* @returns `cases` and `combinations` keyed by id, `envelope` over the combinations,
|
|
168
176
|
* `calculationDetails`, `units` and `warnings`
|
|
169
|
-
* @throws {FeaError} 301–317, 326 and
|
|
177
|
+
* @throws {FeaError} 301–317, 326, 327 and 329–331 for an invalid model (with `details` and
|
|
170
178
|
* every problem in `issues`), 328 for a translation nothing resists at a node connected only
|
|
171
179
|
* to shells, 304 for an unstable structure
|
|
172
180
|
*/
|