From 0428759b1369bbaccc4433f8623c2e2ddc6183e0 Mon Sep 17 00:00:00 2001 From: "KOKO\\Mimi" Date: Wed, 12 Aug 2026 02:36:46 +0900 Subject: [PATCH] docs: rerun MITC4 numerical review --- .../linear-static-mitc4-shell-review.md | 176 +++++++++++++----- 1 file changed, 134 insertions(+), 42 deletions(-) diff --git a/docs/numerical-reviews/linear-static-mitc4-shell-review.md b/docs/numerical-reviews/linear-static-mitc4-shell-review.md index fb97a65..3c632fd 100644 --- a/docs/numerical-reviews/linear-static-mitc4-shell-review.md +++ b/docs/numerical-reviews/linear-static-mitc4-shell-review.md @@ -7,24 +7,34 @@ - source_research: `docs/research/linear-static-mitc4-shell-research.md` - source_formulation: `docs/formulations/mitc4-shell-formulation.md` - source_commits: requirements `5c3025a`, research `02680dc`, formulation revision - `6cb0302`, prior numerical review `0a5aad4` + `6cb0302`, prior numerical review `ebb2657`, I/O definition `73dc761`, reference + metadata policy `2bd297d`, S4 legacy artifacts `57122b0`, observed S4R artifacts + `a28dba4` - source_formulation_sha256: `6DD807E7D9D02C2454CAAF0289E391F5812CC309CDF1C00E92ACEAA97841957D` +- source_io_definition: `docs/io-definitions/linear-static-mitc4-shell-io.md` +- source_io_definition_sha256: `081AB60E504382F2A52C3170FF1F599A1B3E312143065E80F7AD467558958BA3` - status: `needs-reference-model` -- review_iteration: `2` +- review_iteration: `3` - owner_agent: `numerical-review-agent` -- date: `2026-08-11` +- date: `2026-08-12` - implementation_planning_authorized: `false` - review_scope: `현재 선형 정적 formulation 전체 및 미래 기하비선형 tangent 절의 수학 구조` -- reference_artifact_state: rerun workspace에서 `reference/shell/` 미관찰; - reference artifact는 실행하거나 수정하지 않음 +- reference_artifact_state: `reference/shell/` S4 4 files and `reference/shellR/` + S4R 4 files observed read-only; both lack `metadata.json`, which is + `absent-allowed`; both also lack the generally required bundle `README.md`, and + the feature-specific Reference Model Contract is absent ## Review Verdict - verdict: `needs-reference-model` - primary_reason: 1차 리뷰의 confirmed formulation defect `NR-C01`-`NR-C05`와 결정 `NR-D01`-`NR-D02`는 revision `6cb0302`에 수학적으로 일관되게 반영되었다. - 그러나 Implementation Planning 전에 승인해야 하는 drilling/geometry/director - calibration 및 Abaqus reference-model evidence가 아직 없다. + S4/S4R artifact 존재성도 현재 확인되었다. 그러나 두 bundle은 동일한 + 평면 정규 `6 x 6` mesh의 단일 두께 사례이며, 승인된 MITC4 Reference + Model Contract, provenance/unit/step-frame 계약, U/UR tolerance 제안, + drilling sweep, geometry/director separation sweep이 없다. 따라서 존재하는 + artifact는 source-type presence만 닫고 Implementation Planning 수치 gate를 + 닫지 못한다. - secondary_dependencies: `needs-numerical-calibration` - current_linear_kernel_assessment: MITC4의 20-DOF 물리 kinematics, covariant shear tying, plane-stress constitutive law, `2 x 2 x 2` integration, 24-DOF 변환 및 recovery의 @@ -35,9 +45,11 @@ 미래 기하비선형 구현은 계속 금지한다. 현재 선형 정적 정식화에서 새 confirmed mathematical defect는 발견되지 않았다. -그럼에도 요구조건이 Implementation Planning 전 승인을 요구하는 수치 상수와 reference -evidence가 없으므로 `implementation_planning_authorized=false`를 유지한다. 다음 rerun은 -Reference Model과 수치 calibration evidence가 준비된 뒤 수행해야 한다. +다만 `pass-for-implementation-planning`은 미정 상수를 임의 default로 넘겨도 +된다는 뜻이 아니다. 요구조건 `036`, `056`, `058`-`062`와 formulation +Section 21.1이 명시적으로 요구하는 evidence가 없으므로 +`implementation_planning_authorized=false`를 유지한다. 다음 rerun은 승인된 +Reference Model Contract와 수치 calibration evidence가 준비된 뒤 수행해야 한다. ## 1. Review Method and Evidence Boundary @@ -50,7 +62,12 @@ Reference Model과 수치 calibration evidence가 준비된 뒤 수행해야 한 독립적으로 재유도했다. 4. 1차 리뷰의 geometry collapse, drilling 박판 극한, 혼합 DOF scaling 및 nonlinear coordinate-map 반례를 수정본에 다시 적용했다. -5. 확인된 결함의 해소 여부, 알려진 위험, 근거가 부족한 open decision을 분리했다. +5. I/O 계약과 project-wide optional `metadata.json` policy를 현재 artifact + inventory에 대조했다. +6. 두 input의 키워드, mesh/property/load 동일성, CSV row count, + finite/unique node identity, hash와 displacement component scale을 read-only로 + 확인했다. +7. 확인된 결함의 해소 여부, 알려진 위험, 근거가 부족한 open decision을 분리했다. 코드, Harness, Abaqus 및 reference solver는 실행하지 않았고 formulation 문서는 수정하지 않았다. 이 리포트는 reference comparison이나 release 승인이 아니다. @@ -61,13 +78,21 @@ Reference Model과 수치 calibration evidence가 준비된 뒤 수행해야 한 `Required revision`은 formulation commit `67251e0`에 대한 역사적 기록이며, 현재 판정은 각 항목 끝의 `Rerun disposition`이 지배한다. -| finding | rerun status | revision evidence | remaining boundary | +| finding/blocker | rerun status | current evidence | remaining boundary | | --- | --- | --- | --- | | `NR-C01` | resolved | Sections 9.2-9.3, 16.2, 17.1/17.3 | fixed measure threshold calibration | | `NR-C02` | resolved | Sections 12.2, 17.4 | drilling family/plateau calibration | | `NR-C03` | resolved | Sections 12.3/12.5, 17.1 | numerical rank/conditioning evidence | | `NR-C04` | resolved | Sections 7.1-7.2 | none for current linear formulation | -| `NR-C05` | resolved at scope boundary | Sections 15.1-15.5, 17.6 | actual nonlinear chart/update remains future-only | +| `NR-C05` | resolved | Sections 15.1-15.5, 17.6 | current linear scope only; actual nonlinear chart/update remains future-only | +| `NR-O01` | open | dimensionless `rho_d` contract only | family, sweep, stable plateau, nominal coefficient, scaled rank/conditioning bound | +| `NR-O02` | open | separate energy and zero-denominator classification only | warning threshold and physical-response contamination bound | +| `NR-O03` | open | current artifacts have a flat common normal field | curved-mesh incident-normal distribution and hard `theta_smooth` acceptance | +| `NR-O04` | open | fixed geometry-measure inventory and one regular planar mesh | valid distortion/warp versus collapsed/inverted separation sweep and thresholds | +| `NR-O05` | partially resolved | S4/S4R U/UR rows and component scales are observable | approved unit/provenance contract, mesh/convergence evidence, U mixed tolerance, UR warning threshold | +| reference source-type presence | resolved | S4 and S4R inputs plus displacement CSVs are both present | presence is not usability or numerical calibration | +| Reference Model Contract/artifact usability | open | no feature-specific contract; both bundle READMEs are absent; S4R path is noncanonical under current I/O contract | approved alias/canonical disposition, required README handling, provenance, units, step/frame, schema, tolerance | +| `metadata.json` policy | resolved | absent in both bundles and optional by project policy | absence is `absent-allowed`; only a present/conflicting file would require action | ### NR-C01 — Resolved — Jacobian geometry measure completeness @@ -482,35 +507,71 @@ $$ ### 4.2 Open decisions — blocking evidence gaps -#### NR-O01 — Drilling reference family, nominal coefficient, and plateau +#### NR-O01 — Open — Drilling reference family, nominal coefficient, and plateau NR-C02/NR-C03을 반영한 normalized sweep 결과가 없으므로 후보 A/B/physical-block 중 어느 것도 승인하지 않는다. Exact-arithmetic positivity만으로 practical factorization과 physical contamination을 동시에 보장할 수 없다. -#### NR-O02 — Drilling energy warning ratio +#### NR-O02 — Open — Drilling energy warning ratio Source는 static drill-specific acceptable ratio를 제공하지 않는다. Nominal plateau에서 `U/N/M/Q` 민감도와 함께 calibrate해야 하며, 현 단계에서 `1%`, `5%` 등의 임의 값을 승인하지 않는다. `Ephys`가 zero/near-zero이면 ratio를 만들지 않고 두 energy와 case classification을 보고하는 formulation 원칙은 승인한다. -#### NR-O03 — Smooth-director angle +#### NR-O03 — Open — Smooth-director angle `20 degrees`는 usable initial candidate이나 Abaqus grouping heuristic을 FESA hard -rejection으로 바꾸는 결정이다. Pinched cylinder와 LE3의 coarse-to-fine mesh에서 incident -normal deviation 분포와 false-fold/false-smoothing case를 확인하기 전에는 승인하지 -않는다. +rejection으로 바꾸는 결정이다. 현재 S4/S4R artifact는 모두 평면이어서 +incident-normal deviation이 이 경계를 시험하지 않는다. Pinched cylinder와 LE3의 +coarse-to-fine mesh에서 incident normal deviation 분포와 +false-fold/false-smoothing case를 확인하기 전에는 승인하지 않는다. -#### NR-O04 — Geometry thresholds +#### NR-O04 — Open — Geometry thresholds -NR-C01의 measure inventory는 고정되었다. 이제 valid distortion/warp sweep과 collapsed -negative sequence가 pass/fail 사이에 분리 영역을 제공해야 threshold를 승인할 수 있다. +NR-C01의 measure inventory는 고정되었다. 현재 artifact는 평면 정규 mesh의 +한 점만 제공하므로 threshold 근거가 아니다. Valid distortion/warp sweep과 +collapsed/inverted negative sequence가 pass/fail 사이에 분리 영역을 제공해야 +threshold를 승인할 수 있다. -#### NR-O05 — U tolerance and UR warning threshold +#### NR-O05 — Partially resolved — U tolerance and UR warning threshold -이 값은 Reference Model과 공동 소유이며 artifact/mesh evidence가 없다. U-only blocking, -UR warning-only라는 정책은 승인하지만 수치값은 승인하지 않는다. +이 값은 Reference Model과 공동 소유이다. 현재 artifact로 S4/S4R 각 49개 +U/UR node row와 component scale은 관찰할 수 있어 presence 부분은 해소되었다. +그러나 unit label, approved provenance, mesh/convergence uncertainty, tolerance proposal가 +없으므로 수치값은 승인하지 않는다. U-only blocking, UR warning-only 정책만 +유지한다. + +### 4.3 Read-only reference artifact re-inspection + +두 bundle은 수정 없이 다음과 같이 관찰되었다. + +| model candidate | input identity | observed result inventory | numerical relevance | contract status | +| --- | --- | --- | --- | --- | +| `reference/shell/` | `shell.inp`, SHA-256 `4005851E1AB22FD3A16AC17A8D5DA3E051233F69F37419079F3553AD134ECFCF`, `TYPE=S4` | displacement 49 rows, reaction 49 rows, stress 288 rows | S4 source presence and one regular flat-plate U/UR scale | I/O-approved S4 legacy alias, but feature Reference Model Contract and bundle README are absent | +| `reference/shellR/` | `shellR.inp`, SHA-256 `1325940FB42B78961CF25E84379BF2693846FAD22473E7688AC5456B37B18CB4`, `TYPE=S4R` | displacement 49 rows, reaction 49 rows, stress 288 rows | S4R source presence and paired same-model U/UR scale | observed noncanonical bundle; no approved canonical/legacy-alias disposition and no bundle README | + +`shell.inp` and `shellR.inp`는 job name과 source `TYPE`만 다르다. 두 모델은 +49 nodes, 36 elements의 평면 정규 `6 x 6` mesh, `t=0.5`, `E=2.1e11`, +`nu=0.3`, 중앙 절점 `-100000` global-Z CLOAD와 전체 둘레 6-DOF 구속을 +공유한다. 중앙 절점 `U3`는 S4 `-2.37408203e-5`, S4R +`-2.45598239e-5`이다. 각 displacement/reaction CSV의 node key는 49개이며 +중복이 없고 수치값은 유한하다. + +이 inventory는 새로 중요한 정보이지만 calibration portfolio는 아니다. 두 +deck은 단일 두께의 regular planar case이므로 shear/membrane locking trend, +distortion/warp rejection, curved director angle, drilling plateau 및 mesh convergence를 +판정할 수 없다. S4/S4R Abaqus 차이는 두 Abaqus formulation이 다름을 보여줄 +뿐 FESA drilling/geometry 상수를 선택하지 않는다. + +두 bundle 모두 `metadata.json`이 없다. 이는 project-wide 정책상 +`absent-allowed`이며 불완전 판정 근거가 아니다. 다만 선택 JSON의 부재가 +Reference Model Contract이 소유해야 할 provenance, units, coordinate system, +step/frame, schema, tolerance를 없애지는 않는다. 현재는 그 contract 자체가 +없으며, 일반 bundle 정책이 요구하는 `README.md`도 두 경로 모두에 없다. +`metadata.json` 부재와 이 두 필수 계약 공백은 서로 구분한다. 따라서 현재 두 +bundle은 numerical acceptance evidence로 사용할 수 없다. ## 5. Numerical Risk Assessment @@ -558,8 +619,9 @@ UR warning-only라는 정책은 승인하지만 수치값은 승인하지 않는 - patch_tests: Section 7.2의 membrane, bending, shear, twist 및 rigid/drill states. - mms_or_mes: CLOAD-only 범위에서는 prescribed polynomial element state를 이용한 manufactured element solution을 사용하며 distributed-load full-domain MMS는 `N/A`다. -- benchmark_reference_comparison: Section 7.3의 pinched cylinder, LE3 및 S4/S4R - displacement evidence. +- benchmark_reference_comparison: Section 7.3의 pinched cylinder, LE3 및 승인된 + S4/S4R displacement evidence. 현재 평면 정규 쌍은 source-type presence 자료이며 + 단독 benchmark portfolio로는 불충분하다. - missing_evidence: Section 7.5의 drilling, geometry, director-angle, reference 및 future nonlinear evidence. @@ -593,7 +655,8 @@ Distributed load가 없어도 prescribed nodal field와 element action으로 pat - regular/distorted pinched cylinder sequence; - NAFEMS LE3 coarse-to-fine sequence와 incident-normal angle inventory; - 필요 시 equivalent nodal loads를 명시한 Scordelis-Lo 보조 sequence; -- S4 및 S4R source artifact 각 하나 이상. +- S4 및 S4R source artifact 각 하나 이상. 현재 두 regular-flat artifact는 + presence를 닫지만 thin/thick/distorted/curved coverage를 대체하지 않는다. [[Shell Locking Phenomenon]]에 따라 한 개 mesh의 displacement 일치만으로 locking absence를 주장하지 않는다. [[MITC4 Shell Element]]의 patch 및 Scordelis-Lo 근거는 @@ -617,7 +680,14 @@ Requirements/Formulation으로 되돌려 재검토해야 한다. - dimensionally normalized drilling sweep 결과; - fixed geometry metric의 valid/invalid separation data; - curved benchmark별 incident-normal angle distribution; -- reference `S4`/`S4R` artifact inventory와 U/UR scale. +- approved MITC4 Reference Model Contract; +- S4R noncanonical bundle의 승인된 canonical/legacy-alias disposition; +- 두 bundle의 required `README.md` 처리; +- reference provenance, units, coordinate system, step/frame, schema, tolerance; +- U mixed tolerance와 UR warning threshold의 mesh/convergence 근거. + +S4/S4R artifact presence와 raw U/UR component scale 관찰은 더 이상 missing이 아니다. +`metadata.json`도 missing evidence가 아니며 부재를 `absent-allowed`로 기록한다. Future nonlinear chart/update 및 objective drilling evidence도 아직 없지만 현재 linear gate의 dependency는 아니다. 이를 근거로 기하비선형 구현을 시작해서는 안 된다. @@ -629,8 +699,11 @@ gate의 dependency는 아니다. 이를 근거로 기하비선형 구현을 시 | current-linear formulation consistency | pass | 추가 formulation revision 불필요 | | `NR-C01`-`NR-C05` disposition | resolved | 1차 confirmed defect 종료 | | `NR-D01`-`NR-D02` integration | pass | 수식/normalized metric 유지 | -| drilling/geometry/director calibration | pending | Implementation Planning 차단 | -| reference model and U/UR thresholds | pending | primary verdict `needs-reference-model` | +| `NR-O01`-`NR-O04` drilling/geometry/director calibration | open | Implementation Planning 차단 | +| S4/S4R artifact presence | resolved | 두 source type의 read-only bundle 존재 확인 | +| `NR-O05` U/UR scale observation | partially resolved | raw rows/scales는 존재하나 threshold 승인 불가 | +| reference model contract and artifact usability | open | primary verdict `needs-reference-model` | +| optional `metadata.json` absence | resolved | `absent-allowed`; gate 차단 아님 | | future geometric-nonlinear executability | out of current scope | 별도 future formulation/review 전 구현 금지 | ## 8. Required Revisions @@ -653,20 +726,31 @@ gate의 dependency는 아니다. 이를 근거로 기하비선형 구현을 시 ### Reference Model Agent -- 현재 workspace에 관찰되지 않는 `reference/shell/`을 artifact가 다시 보일 때 - read-only로 inventory하고 S4R 선언 artifact의 실제 파일, provenance, input subset 및 - displacement row를 확인한다. - `docs/reference-models/linear-static-mitc4-shell-reference-models.md`에 canonical portfolio와 artifact/schema 계약을 작성한다. -- 별도 S4 artifact를 준비한다. +- `reference/shell/`을 exact S4 legacy alias로 등록하고, 관찰된 + `reference/shellR/` noncanonical bundle을 수정·rename하지 않은 채 승인된 + legacy alias로 취급할지 아니면 canonical S4R artifact가 별도 필요한지 + 명시적으로 판정한다. +- 두 bundle에 일반 artifact 정책상 필요한 `README.md`가 없음을 기록하고, + reference artifact 변경 권한이 있는 별도 승인 단계로 처리한다. Numerical + Review가 이 파일들을 생성하거나 legacy 예외를 새로 승인하지 않는다. +- 두 bundle의 provenance, units, coordinate system, model/step/frame, schema, + tolerance를 계약에 고정한다. `metadata.json`은 생성 요구하지 않으며 + 부재를 `absent-allowed`로 기록한다. - LE3/pinched-cylinder mesh마다 incident-normal angle, mesh size 및 U reference scale을 - metadata/evidence로 남긴다. + Reference Model Contract/evidence로 남긴다. +- flat/thin/thick/distorted/curved 포트폴리오와 valid/invalid geometry sequence를 + 정의해 현재 단일 regular-flat pair의 coverage 한계를 닫는다. - U mixed tolerance와 nonblocking UR warning threshold를 제안한다. ### Numerical Review Rerun - Formulation revision `6cb0302`의 defect correction은 이번 rerun에서 확인했다. -- 다음 rerun에서 geometry negative sequence와 drilling normalized sweep evidence를 +- S4/S4R presence와 optional-metadata policy는 이번 rerun에서 현재 inventory로 + 다시 판정했다. +- 다음 rerun에서 승인된 Reference Model Contract, geometry negative + sequence, drilling normalized sweep, director-angle distribution 및 tolerance 제안을 검토한다. - Section 4.2의 open decision을 닫거나 명시적으로 다시 반환한다. - 그 전에는 `pass-for-implementation-planning`을 부여하지 않는다. @@ -682,16 +766,20 @@ gate의 dependency는 아니다. 이를 근거로 기하비선형 구현을 시 ### I/O Definition Agent -- formal I/O gate handoff는 이번 verdict로 승인하지 않는다. 비권위 schema 메모를 - 작성하더라도 geometry threshold, - drilling coefficient/energy warning, U/UR tolerance를 implementation default로 - 선점하지 않는다. +- 현재 I/O definition은 이 rerun의 reviewed input이지만 Numerical Review 통과를 + 의미하지 않는다. Geometry threshold, drilling coefficient/energy warning, + U/UR tolerance를 implementation default로 선점하지 않는다. +- Section 7.3의 “S4R bundle missing” 서술은 현재 물리 inventory와 다르다. + 다만 `reference/shellR/`의 noncanonical naming과 contract 부재 때문에 “usable + canonical S4R evidence missing” dependency는 남는다. 이 구분을 다음 I/O 메타데이터 + refresh에 반영한다. - NR-D01의 normalized drilling-load projection과 exact-zero handling을 diagnostic contract에 반영한다. ### Reference Model Agent -- Section 8의 read-only inventory 및 benchmark evidence를 준비한다. +- Section 4.3/8의 read-only inventory를 계약으로 승인하고 누락 benchmark + evidence를 준비한다. - Abaqus S4/S4R formulation equivalence를 주장하지 않고 global U만 blocking evidence로 사용한다. @@ -709,6 +797,10 @@ gate의 dependency는 아니다. 이를 근거로 기하비선형 구현을 시 - `docs/formulations/mitc4-shell-formulation.md` - `docs/research/linear-static-mitc4-shell-research.md` - `docs/requirements/linear-static-mitc4-shell.md` +- `docs/io-definitions/linear-static-mitc4-shell-io.md` +- `docs/reference-models/README.md` +- `reference/shell/` S4 input/displacement/reaction/stress bundle (read-only) +- `reference/shellR/` S4R input/displacement/reaction/stress bundle (read-only) - `docs/reference-papers/MITC4/AContinuumMechanicsBasedFourNodeShell/` - `docs/reference-papers/MITC4/FourNodeQuadrilateralShellElementMITC4/` - `docs/reference-papers/MITC4/유한요소해석법을이용한쉘구조물의동적좌굴해석/`