2 차 16 진 유한 요소에 8 가우스 포인트가 필요합니까?


10

비 물리적 모드를 도입하지 않고 8 가우스 포인트 미만의 육면체 유한 요소에 대해 2 차 정확도를 얻을 수 있습니까? 단일 중앙 가우스 포인트는 비 물리적 전단 모드를 도입하고 8 가우스 포인트의 표준 대칭 배열은 사면체 이산화에 비해 비쌉니다.

편집 : 누군가가 방정식을 요구했습니다. 내가 관심있는 방정식은 동적 또는 준 정적 비선형 탄성입니다. 준 정적 방정식은

P(ϕ)=0

여기서 , Ω R 3P : R 3 × 3R 3 × 3 은 초 탄성 첫 번째 Piola-Kirchoff 응력 함수입니다. 간단한 예는 압축 형 신형 Hookean입니다. 여기서 P ( F ) = μ ( F - F - T ) + λ F - T log det Fϕ:ΩR3ΩR3P:R3×3R3×3

P(F)=μ(FFT)+λFTlogdetF

정확히 무엇을 시뮬레이션하고 있습니까?
Dan

현재 선형 탄성이지만 문제는 일반적으로 비선형 탄성에 관한 것입니다.
Geoffrey Irving

1
"비 물리적"의 정의는 이들에 의존하기 때문에 관심있는 방정식을 포함시켜야합니다. 또는 최소한 "물리적"인 기능의 공간을 정확하게 정의하십시오.
David Ketcheson

방정식이 추가되었습니다.
Geoffrey Irving

1
dPhi / dx를 사용하면 그라디언트를 의미합니까?
Wolfgang Bangerth

답변:



0

일반적으로 자유도가있는 것보다 셀당 더 적은 직교 점으로 도망 갈 수 없다는 것이 비교적 분명합니다. 3 차원 육면체의 삼선 형 요소의 경우 8 자유도 (정점 당 1 개)가 있으므로 최소 구적 점 수는 8 개가됩니다.

그것은 돌이킬 수 없으며 결과적으로 완전히 쓸모가 없습니다. 그 이유는 1 점 직교 공식이 직교 점에서 동일한 값을 갖는 모든 선형 함수 (시험 공간의 일부)를 구별 할 수 없기 때문입니다. 즉, 중간 점 규칙의 경우 모양 함수 'x'는 함수 '0'이 함수 '-x'와 동일합니다. 다시 말해, 시험 공간에는 정확한 적분이있는 차원 2가 있지만 중간 점 규칙의 경우 자유도가 2 개인데도 공간이 차원 1입니다. 중간 점 규칙의 경우 모양 함수 'x'는 함수 '0'과 동일하고 함수 '-x'와 동일합니다. 다시 말해, 시험 공간에는 정확한 적분이있는 차원 2가 있지만 중간 점 규칙의 경우 자유도가 2 개인데도 공간이 차원 1입니다. 중간 점 규칙의 경우 모양 함수 'x'는 함수 '0'과 동일하고 함수 '-x'와 동일합니다. 다시 말해, 시험 공간에는 정확한 적분이있는 차원 2가 있지만 중간 점 규칙의 경우 자유도가 2 개인데도 공간이 차원 1입니다.


Geoff의 질문이 더 미묘하다고 생각합니다. 잘 형성된 도메인의 사면체에 연속 유한 요소 공간이있는 경우 (예를 들어 격리 된 요소가없는 경우) 명확하게 통합되지 않은 단일 점 구적법으로 벗어날 수 있습니다. 문제는 어떤 방식으로도 육면체 요소와의 통합이 가능한지 여부입니다. 답을 모르지만 직교 점에 추가 메모리 동작이 필요하지 않으므로 얼마나 큰지 잘 모르겠습니다. 유한 요소 잔차 평가를 벡터화하면 메모리에 바인딩되는 것이 일반적이므로 플롭을 사용하는 것이 좋습니다.
Jed Brown

메모리 동작에 대한 좋은 지적.
Geoffrey Irving

1
3×3

1
주석에는 줄 바꿈을 포함 할 수 없다는 불편 함이 있습니다.
Geoffrey Irving

@JedBrown : 좋은 지적입니다. tet에서 선형 함수의 기울기는 상수이므로 질량 행렬에 대해 만든 인수에 따라 단일 직교 점이면 충분합니다 (강성 행렬은 기울기의 질량 행렬 :-). 반면에, 육면체의 삼선 함수의 기울기는 (이방성) 이차 함수이므로 좌표 방향 당 하나 이상의 직교 점 이상이 필요합니다.
Wolfgang Bangerth
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.