Summary
이 연구는 유한 요소 모델을 얻기 위해 상악과 상악치아의 저용량 3 차원 원뿔 빔 기반 환자 이미지를 활용하는 데 필요한 도구를 간략하게 설명합니다. 이 참을성 있는 모형은 정확하게 모든 상악 치아의 CRES를 찾아내기 위하여 이용됩니다.
Abstract
저항의 중심 (CRES)는예측 가능한 치아 운동을위한 기본 기준점으로 간주됩니다. 치아의 CRES를 추정하는 데 사용되는 방법은 전통적인 방사선 및 물리적 측정에서 모델 또는 시체 표본에 대한 시험관 내 분석에 이르기까지 다양합니다. 모델과 단일 치아의 고용량 마이크로 CT 스캔의 유한 요소 분석과 관련된 기술은 많은 약속을 보여 주었지만, 새로운 저용량 및 저해상도 콘 빔 컴퓨터 단층 촬영 (CBCT) 이미지로는 거의 수행되지 않았습니다. 또한, 소수의 선택된 치아(즉, 상악중앙 절개, 개 및 제1 대구치)에 대해서만 CRES가 기술되었다; 나머지는 대부분 무시되었습니다. 또한 CRES를 보다 자세히 결정하는 방법을 설명하여 복제및 구축이 용이해질 수 있습니다.
이 연구는 도구 개발을 위한 일상적인 CBCT 환자 이미지를 사용하였고, 상악치아의 CRES를 찾기 위한 유한 요소 모델을 얻기 위한 워크플로우를 사용했습니다. CBCT 볼륨 이미지는 세분화에 의해 상악치아의 CRES를 결정하는 데 관련된 3차원(3D) 생물학적 구조를 추출하도록 조작하였다. 분할된 오브젝트는 3matic 소프트웨어로 최대 가장자리 길이가 1mm인 사각지대(tet4) 삼각형으로 구성된 가상 메시로 정리및 변환되었습니다. 유한 요소 해석에 사용하기 위해 최대 모서리 길이가 1mm인 테트라헤드론의 솔리드 체적 메시로 모델이 추가로 변환되었습니다. 엔지니어링 소프트웨어인 Abaqus는 모델을 전처리하여 어셈블리를 작성하고 재료 특성, 상호 작용 조건, 경계 조건 및 하중 응용 프로그램을 설정하는 데 사용되었습니다. 하중을 분석할 때 시스템의 응력과 변형을 시뮬레이션하여 CRES를찾는 데 도움을 주어. 이 연구는 치아 움직임을 정확하게 예측하는 첫 번째 단계입니다.
Introduction
치아 또는 치아 세그먼트의 저항 중심 (CRES)은자유 신체의 질량 중심과 유사합니다. 그것은 강체역학의 분야에서 빌린 용어입니다. CRES에서단일 힘이 가해지면, 힘의 작용 선 방향으로 치아의 번역이1,,2가발생합니다. CRES의 위치는 치아의 해부학 및 특성뿐만 아니라 환경 (예 : 치주 인대, 주변 뼈, 인접 치아)에 따라 달라집니다. 치아는 절제된 몸체로, CRES는 자유 신체의 질량 중심과 유사합니다. 가전 제품의 조작에서 대부분의 치열 교정 의사들은 치아 또는 치아 그룹의 CRES에 대한 힘 벡터의 관계를 고려합니다. 실제로, 단일 힘으로 제출될 때 물체가 기울어지거나 신체적 움직임을 표시할지 여부는 주로 물체의 CRES의 위치와 힘 벡터와 CRES사이의 거리에 의해 결정된다. 이것이 정확하게 예측될 수 있는 경우에, 처리 결과는 크게 향상될 것입니다. 따라서, CRES의 정확한 추정은 치열 교정 치아 운동의 효율성을 크게 향상시킬 수 있습니다.
수십 년 동안, 치열 교정 필드는 주어진 치아의 CRES의 위치에 관한 연구를 재검토하고있다, 세그먼트, 또는 아치1,2,2,3,,4,,5,6,7,,88,,10,,911,,12. 그러나, 이러한 연구는 여러 가지 방법으로 그들의 접근에 제한 되었습니다. 대부분의 연구는 단지 몇 치아에 대 한 CRES를 결정, 대부분을 떠나. 예를 들어, 상악 중앙 절개와 상악 절치 세그먼트는 매우 광범위하게 평가되었습니다. 다른 한편으로는, 상악가 개와 첫번째 대구치에 단지 몇몇 연구 결과가 있고 나머지 치아를 위한 아무 것도 없습니다. 또한, 이러한 연구의 대부분은 치아에 대한 일반적인 해부학 적 데이터, 2 차원 (2D) 방사선 사진에서의 측정 및 2D 도면8에 대한 계산을 기반으로 CRES의 위치를결정했습니다. 또한, 현재 문헌의 일부는 인간의 데이터4,,8이아닌 덴티폼 모델의 일반 모델 또는 3차원(3D) 스캔을 사용한다. 치열 교정이 치아 이동을 계획하기 위한 3D 기술로 전환됨에 따라 이 개념을 다시 검토하여 치아 움직임에 대한 3D 과학적 이해를 개발하는 것이 중요합니다.
기술 발전으로 인해 컴퓨팅 성능과 모델링 기능이 향상되면서 보다 복잡한 모델을 만들고 연구할 수 있는 능력이 향상되었습니다. 컴퓨터 단층 촬영 스캐닝 및 콘 빔 컴퓨터 단층 촬영 (CBCT) 스캐닝의 도입은 2D 세계에서 3D로 추력 모델과 계산을 가지고있다. 컴퓨팅 성능과 소프트웨어 의 복잡성의 동시 증가는 연구원이 치아, 뼈, 치주 인대 (PDL) 및 다양한 다른 구조를 분할하는 고급 소프트웨어에 사용하기 위해 정확한 해부학 모델을 추출하기 위해 3D 방사선 그래프를 사용할 수 있게했다7,,8,9,10,,,13,,10,14, 15.14 이러한 분할된 구조는 지정된 힘 또는 변위가 적용될 때 시스템의 응답을 계산하기 위해 엔지니어링 소프트웨어에 사용하기 위해 가상 메시로 변환할 수 있습니다.
이 연구는 살아있는 환자의 CBCT 심상에서 파생된 모형에 적용된 가상 치열 교정 힘 시스템을 검토하기 위하여 이용될 수 있는 특정, 복제 방법론을 제안합니다. 이 방법론을 활용하여, 연구원은 다양한 치아의 CRES를 추정하고 치아 해부학, 뿌리 수 및 3D 공간에서의 뿌리 수 및 배향과 같은 치과 구조의 생물학적 형태를 고려할 수 있습니다. 이 프로세스의 일반적인 개요는 그림 1에나와 있습니다. 이는 CRES를찾기 위한 3D 치아 모델의 생성에 관여하는 논리적 프로세스로 판독기를 지향하는 것이다.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
구강 및 악안면 방사선학 (IRB No. 17-071S-2)의 부서에 보관 된 CBCT 볼륨을 평가하기위한 기관 검토 위원회 면제가 획득되었습니다.
1. 볼륨 선택 및 기준
- 머리와 얼굴(16)의CBCT 이미지를 획득합니다.
- 치아 정렬, 누락된 치아, 복셀 크기, 시야 및 이미지의 전체 품질에 대한 이미지를 검사합니다.
- 복셀 크기가 350μm(0.35mm)를 초과하지 않는지 확인합니다.
2. 치아와 뼈의 분할
- CBCT 이미지의 원시 DICOM 파일을 세분화용 모방 소프트웨어에 로드합니다(그림2). 이미지 > 자르기 프로젝트를클릭합니다. 상악 치아와 상악 치아만 포함하도록 이미지를 자작합니다.
참고: 시야는 상악과 상악 치아를 포착할 수 있을 만큼 커야 합니다. 이미지에는 치아 크라운, 비강 바닥까지의 단단한 구개, 상악 치아의 얼굴 표면 및 경솔한 구개 및 상악점 결절의 후방 범위가 포함되어 있는지 확인하십시오. - 마스크 탭을 마우스 오른쪽 버튼으로 클릭하고 이미지에 대한 새 마스크를 만듭니다. 마스크의 이름을 UL1, UL2, ..., UL7의 왼쪽과 UR1, UR2, ..., 오른쪽의 UR7로 변경합니다.
- 마스크된 CBCT 이미지에서 관심 있는 치아를 식별합니다(보기 참조). 마스크 지우기 도구를 사용하여 마스크를 지웁울 수 있습니다. 두 개의 회색 값이 유사하기 때문에 소프트웨어가 치아와 골격을 구분하지 못할 수 있습니다.
참고: Mimics의 임계값 도구는 치아와 골격을 별도로 분할할 수 없습니다. 따라서 세분화에 대 한 다른 방법이 필요 합니다. - 다중 슬라이스 편집 도구(Ctrl + M)를 클릭합니다. 뷰(축축, 코로나또는 시상)를선택합니다. 필요한 것으로 간주되는 일부 조각을 수동으로 강조 표시(즉, 그리기)합니다.
참고: 슬라이스를 더 강조 표시하려면 구조에 세부 정보가 더 많이 추가됩니다. - 보간 도구를 클릭하여 건너뛴 슬라이스의 볼륨을 채우고 적용합니다.
- 마스크를 마우스 오른쪽 단추로 클릭하고 3D 체적을 계산하는 옵션을 선택하여 치아의 3D 볼륨을 생성합니다.
- 상악아치의 각 치아에 대해 2.2-2.6단계를 반복합니다.
- 모든 3D 상악치아 UL7-UR7을 선택합니다. 마우스 오른쪽 단추를 클릭하여 스무딩을선택합니다. 스무딩 계수를 0.4로 설정하고 반복을 4로 설정합니다.
- 상악골격을 분할하려면 마스크의탭을 마우스 오른쪽 버튼으로 클릭합니다. 이미지에 대한 새 마스크를 만듭니다.
- 미리 정의된 임계값 집합에 대한 드롭다운 메뉴에서 사용자 지정을 선택합니다. 전체 상악 골본을 포함하도록 임계값을 조정합니다. 임계값을 적용하기 전에 구멍 채우기 상자를 선택하십시오.
참고 : 피질 뼈에 ≤1mm의 작은 구멍은 이후 단계에서 쉽게 제거 할 수 있기 때문에 허용됩니다. - 동적 영역 성장 도구를 클릭하여 마스크에 표시되는 큰 구멍을 채웁니다. 다중 레이어 상자를 선택하는 것 외에도 상악 골 마스크를 도구의 대상으로 선택합니다. 최소값은 50, 최대 값은 150을 사용합니다. 마스크에서 강조 표시되지 않은 피질 골격 영역을 클릭하면서 컨트롤 키를 누누적입니다.
- 매끄러운 마스크 기능에 대한 상악 골격 마스크를 마우스 오른쪽 단추로 클릭합니다. 최상의 결과를 얻으려면 이 단계를 3x 반복합니다.
- 마스크를 마우스 오른쪽 단추로 클릭하고 3D 체적을 계산하는 옵션을 선택하여 상악골의 3D 볼륨을 생성합니다.
- 3D 상악골을 선택합니다. 마우스 오른쪽 단추로 클릭하여 스무딩을 선택합니다. 스무딩 계수를 ~0.4로 설정하고 반복을 4로 설정합니다.
- 3D 상악 골격을 선택하고 마우스 오른쪽 단추로 클릭하여 둘러싸기를선택합니다. 가장 작은 디테일은 0.2mm, 간격 닫기 거리는 1mm로 설정합니다. 얇은 벽 보호 옵션을 확인합니다. 확인을누릅니다.
- 3D 상악골의 이름을 "맥실라"로 바꿉니다.
3. 청소 및 메쉬
- 3D 객체를 선택하고 복사(Ctrl + C).
- 3matic 소프트웨어를 열고 선택한 3D 객체를 붙여넣습니다(Ctrl + V). 이러한 개체 트리 및 작업 영역 3D 구조로나타납니다(그림 3).
- 도구 모음에서 수정 탭을 클릭하고 [매끄럽게] 옵션을 사용합니다. 작업 상자에서 원하는 3D 개체 또는 엔터티를 선택하고 기본 매개 변수를 적용합니다.
- 도구 모음에서 완료 탭을 클릭하고 로컬 스무딩 옵션을 사용합니다. 작업 상자에서 원하는 3D 개체 또는 엔터티를 선택합니다. 커서를 사용하여 원하는 영역을 수동으로 부드럽게 합니다.
- 치아를 복제합니다. 개체 트리에서 모든 치아를 선택하고 마우스 오른쪽 단추로 클릭하고 중복을 선택합니다.
- 중복된 모든 치아,그룹을 선택하고 폴더이름을 "그룹 1"으로 지정합니다. 원래 세트는 분석을위한 최종 톱니 역할을합니다.
- 그룹 1의 중복된 톱니를 보려면 곡선 모듈 및 곡선 만들기 옵션을 클릭합니다. 복제된 모든 치아에 대해 시멘로에나멜 접합부(CEJ) 주위에 곡선을 수동으로 그립니다.
- [부드러운 곡선] 옵션 아래에서 곡선, 등고선및 테두리 엔터티를 선택합니다.
- 곡선별 곡면 분할 옵션을 선택하고 선택할 3D 오브젝트를 왼쪽 클릭하여 크라운 및 루트 서피스를 자체 부품으로 분리합니다.
- CEJ에서 치아를 뿌리와 크라운으로 분할하여 치아의 뿌리 구조에서 PDL을 생성합니다.
- 그룹 1(3.6단계에서 생성)의 3D 객체를 그룹 2로 복제합니다. 그룹 2의 경우 개체 트리 상자에서 개체를 클릭합니다. 지표면 목록에서 크라운 서피스를 삭제합니다. 그룹 2의 모든 개체에 대해 이 단계를 수행합니다.
- 그룹 2의 경우 디자인 모듈 > 중공을 클릭합니다. 원하는 매개 변수를 적용합니다(표1).
- 수정 모듈 및 수정 마법사를클릭합니다. 개별 부품을 클릭하고 업데이트한 다음 지정된 지침을 따릅니다.
- 모든 부품에 대해 3.10.3 단계를 반복합니다. 그룹 2의 모든 부분을 "UL1_PDL"에서 "UL7_PDL", "UR1_PDL"을 "UR7_PDL"로 바꿉니다.
- 그룹 1에서 개체 트리 상자에서 개체를 클릭합니다. 지표면 목록에서 루트 서피스를 삭제합니다.
- 구멍 법선 채우기 옵션을 선택하고 등고선 선택을 합니다. 불량 윤곽을 클릭하고 적용합니다. 전체 공간이 채워집니다.
- 설계 모듈 > 로컬 오프셋을 선택하고 전체 크라운 서피스를 선택합니다. 다음 옵션: 방향(외부 선택), 간격띄우기 거리(0.5 선택), 점감 거리(2.0 선택)를 선택합니다. 을 적용합니다.
- 3.13단계를 반복합니다.
- 상악 아치의 각 치아에 대해 3.11-3.14 단계를 반복합니다.
- 리메쉬(그림 3)
- 개체 트리에서 리메시 모듈 > 비매니폴드 어셈블리 > 주 엔티티 > 맥실라 만들기를 클릭합니다. 3.4(원래 톱니)의 모든 객체에 대해 교차 하는 개체를 선택하고 적용을선택합니다.
- 리메시 모듈을 클릭합니다. 비다양형 어셈블리를 분할합니다.
- 교차하는 엔터티를 그룹 1의 모든 개체로 사용하여 3.16.1-3.16.2 단계를 반복하고 적용합니다.
- 선택적 단계로 필요한 경우에만 완료 모듈 > 트림 > 엔티티 > 맥실라를선택합니다. 초과 구조(즉, 노이즈)를 선택하고 적용합니다.
- 수정 모듈을 클릭 > 수정 마법사 > 맥실라 > 업데이트. 주어진 지침을 따르십시오.
- 교차하는 엔터티를 그룹 2의 모든 개체로 사용하고 적용합니다.
- 리메시 모듈 > 적응 형 리메시를클릭합니다. 3.16.6에서 교차하는 모든 엔터티를 선택하고 적용합니다.
- 리메시 모듈 > 분할 비매니폴드 어셈블리를클릭합니다.
- 개체 트리에서 그룹 2에서 Remesh 모듈 > 비다양성 어셈블리 > 주 엔터티 > 개별 개체(PDL) 만들기를 클릭합니다. 교차 하는 개체 를 선택 > 단계 3.4 (해당 치아 유형에 해당)에서 각 개체를 선택 하 고 적용합니다.
- 리메쉬 모듈 > 적응 형 리메시를클릭합니다. 교차 하는 엔터티를 3.16.9에서 선택하고 적용합니다.
- 리메쉬 모듈 > 분할 비매니폴드 어셈블리를클릭합니다.
- 각 치아에 대해 3.16.9-3.16.11 단계를 반복합니다.
- 리메쉬 모듈 > 품질 보존을 클릭하여 삼각형을 줄입니다. 개체 트리에서 모든 엔터티(예: 치아, Pdf 및 맥실라)를 선택하고 적용합니다.
- 리메시 모듈 > 볼륨 메시 만들기 > 엔터티 선택. 메시 매개변수를 선택합니다.
- 모든 엔티티(즉, 치아, PPL 및 맥실라)에 대해 3.18 단계를 반복합니다.
- 입력(.inp) 파일을 3Matic에서 Abaqus로 수동으로 내보냅니다(그림4).
4. 유한 요소 분석
참고 : 모든 사용자 정의 파이썬 스크립트는 보충 첨부 파일에서 찾을 수 있습니다. 그들은 Abaqus에서 매크로 관리자 기능을 사용하여 생성되었습니다.
- 전처리 설정
- 아바쿠스를 열고 표준 모델을선택합니다. 파일 > 작업 디렉토리 설정 > 파일 저장소에 대 한 위치를 선택합니다.
- 파일 > 스크립트를 실행하고 Model_setup_Part1.py를 선택합니다.
- 모델 디렉토리에서 Abaqus에서 .inp 파일을 로드할 파일 경로를 지정합니다.
- 모델 > 시뮬레이션 > 부품 > 맥실라 > 표면을 클릭합니다.
- 대화 상자 "UL1 _socket"에서 표면의 이름을 지정합니다.
- 서피스 영역 선택 에서 각도별로선택합니다. 각도로 "15"를 추가합니다.
- 소켓의 모든 영역이 선택되어 있는지 확인합니다. 완료되면 완료를 누릅니다.
- 개별 소켓에 대해 4.1.4-4.1.7 단계를 반복합니다.
- 모델 > 시뮬레이션 > 부품을클릭합니다. 그런 다음 UL1 > 표면을 선택합니다. 곡면의 이름을 "UL1"로 지정합니다.
- 표면 영역 선택에서 "개별적으로"를 선택합니다. 화면에서 치아를 선택하고 완료를누릅니다.
- 모든 치아에 대해 4.1.9-4.1.10 단계를 반복합니다.
- 모델 > 시뮬레이션 > 부품을클릭합니다. 그런 다음 UL1_PDL > 서피스를 선택합니다. 곡면의 이름을 "UL1_PDL_inner".
- 서피스 영역 선택 에서 각도별로선택합니다. 각도로 "15"를 추가합니다.
참고: 최종 시뮬레이션 중에 오류가 발견되면 각도를 줄이고 표면을 다시 선택합니다. - PDL의 전체 내부 표면적을 선택했는지 확인합니다. 완료되면 완료를 누릅니다.
- UL1_PDL > 서피스를 선택합니다. 서피스의 이름을 "UL1_PDL_outer"으로 지정합니다.
- 서피스 영역 선택 에서 각도별로선택합니다. 각도로 "15"를 추가합니다.
참고: 최종 시뮬레이션 중에 오류가 발견되면 각도를 줄이고 표면을 다시 선택합니다. - PDL의 전체 외부 표면적을 선택했는지 확인합니다. 완료되면 완료를 누릅니다.
- 모든 PPL에 대해 4.13-4.1.19 단계를 반복합니다.
- 파일 > 스크립트를 실행하고 Model_setup_Part2.py를 선택
- 모델 > 시뮬레이션 > BC를클릭합니다. 이름을 BC_all다음 단계를 초기로 선택합니다. 범주에서 "기계"를 선택하고 "선택한 단계 유형"에서 "변위/회전"을 선택합니다. 계속을 누릅니다.
- 경계 조건에 대한 선택 영역 아래에서 각도별로선택합니다. 각도로 "15"를 추가합니다. 세트 만들기를선택합니다. 14개의 톱니에 대한 개별 소켓을 선택합니다. 완료를 누릅니다 .
참고: 이것은 순간적인 치아 운동을 시뮬레이션하는 데 도움이 되었습니다. - 모델 > 시뮬레이션 > 어셈블리 > 세트 > 세트 만들기를클릭합니다. 세트의 이름을 "U1_y_force"으로 지정합니다.
- 집합의 노드 선택에서 개별적으로선택합니다.
참고: 하나의 뉴턴 농축 힘이 양수 Y 방향(분할 힘을 시뮬레이션) 또는 양수 Z 방향(관입력을 시뮬레이션)으로 임의로 선택된 치아 노드에 적용되었습니다. - 상부 중앙 절치(U1)의 볼 표면에서 크라운 의 중심에 있는 노드를 선택하고 완료를 누릅니다.
- 세트 > 세트 만들기를클릭합니다. 세트의 이름을 "U1_z_force".
- 4.1.23-4.1.24 단계를 반복합니다.
- 모든 치아에 대해 4.1.22-4.1.26 단계를 반복합니다.
참고: 4.1.25에서와 같이 특정 치아에 대해 세트가 생성되기 전에 인스턴스 > 해당 치아에 대한 다시 시작으로 이동합니다.
- 모델 설정
- 모델 > 시뮬레이션 > 어셈블리 > 인스턴스를클릭합니다. 모든 인스턴스를 선택하고 다시 시작을 클릭합니다.
- 도구 > 쿼리 > 포인트 / 노드를클릭합니다. 임의로 선택된 중앙 절개기의 중앙에 있는 노드를 선택하고 완료를누릅니다.
- 페이지 하단의 명령 센터 아래에서 4.2.2단계에서 선택한 노드의 X, Y 및 Z 좌표를 복사합니다.
- 수직 도구 모음 아래에서 인스턴스 변환을 선택하고 화면에서 전체 어셈블리(예: 모든 인스턴스)를 선택합니다. 완료를 누릅니다 .
- 번역 벡터의 시작점 선택 상자에서 복사된 좌표를 4.2.3 단계에서 붙여넣거나 X, Y 및 Z 값을 입력합니다. Enter 를 클릭합니다.
- 변환 벡터의 끝점을 선택하거나 X, Y, Z를 입력: 좌표 "0.0", "0.0"및 "0.0"을 입력합니다. Enter 를 클릭합니다.
- 인스턴스 의 위치를보려면 확인을누릅니다.
- 도구 > 쿼리 > 포인트/노드를 클릭하고 중심 절개의 중간선 바로 위에 있는 노드를 선택합니다. 완료를 입력합니다.
- 페이지 하단의 명령 센터 아래에서 4.2.8 단계에서 선택한 노드의 X, Y 및 Z 좌표를 복사합니다.
- 수직 도구 모음 아래에서 인스턴스 변환을 선택하고 화면에서 전체 어셈블리(예: 모든 인스턴스)를 선택합니다. 완료를 입력합니다.
- 복사된 좌표를 변환 벡터의 시작점 선택에 붙여넣거나 X, Y, Z 상자를 입력합니다. Enter 를 클릭합니다.
- 변환 벡터의 끝점 선택에서 X, Y, Z를 입력하거나4.2.9 단계에서 복사된 좌표를 삽입합니다. X 좌표를 0.0으로 변경합니다. Enter 를 클릭합니다.
- 인스턴스 의 위치를보려면 확인을누릅니다.
- 파일 > 스크립트를 실행하고 Model_setup_Part3.py를선택합니다. 재질 특성을 삽입하거나 변경합니다.
- 모델 > 시뮬레이션 > 재질을 클릭하고 뼈 / PDL / 치아를클릭합니다. 조직 특정 특성을 삽입합니다.
- 파일 > 스크립트를 실행하고 Functions.py선택합니다.
- 모델 처리
- 파일 > 스크립트를 실행하고 Job_submission.py를선택합니다.
참고: 작업 모듈은 사용자가 모델에 하나 이상의 작업을 설정하고 작업 관리자는 모델 분석이 시작되고 진행률이 표시되고 완료가 기록되는 위치입니다. - 모두 억제라는제목의 대화 상자에서 구속조건(모델 및 시뮬레이션 > 제약 조건 아래)을 기준으로 치아의 측면(L 또는 R)을 입력합니다. 확인을누릅니다.
- 작업 제출이라는 대화 상자에서 "Y"를 입력하여 지정된 치아/치아에 대한 분석을 실행합니다. 확인을누릅니다.
- 분석 방향이라는 제목의 대화 상자에서 힘 응용 프로그램을 지정하기 위해 "Y"를 입력합니다. 확인을누릅니다.
- 파일 > 스크립트를 실행하고 Job_submission.py를선택합니다.
- C RES 추정을 위한 사후 처리
- 파일 > 스크립트 실행 > Bulk_process.py를선택합니다.
- 여러 작업 분석이라는 대화 상자에서 지정된 치아/치아에 대해 "Y"를 입력합니다. 확인을누릅니다.
- 분석 방향이라는 대화 상자에서 힘 응용 프로그램을 지정하기 위해 "Y"를 입력합니다. 확인을누릅니다.
- 입력 받기라는 제목의 대화 상자에 지정된 인스턴스(예: UL1 또는 UL5 등)로 특정 치아 번호를 입력합니다. 확인을누릅니다.
- 명령 상자에서 점 에 대한 힘 및 예상 위치에 대한 좌표를 선택합니다. 유사하지 않은 경우 4.3.1-4.4.4 단계를 반복합니다.
참고: 각 단계에 대한 작업이 실행된 후 파이썬에서 생성된 사용자 정의 알고리즘이 Abaqus 인터페이스 내에서 실행되어 부하 응용 프로그램의 결과로 생성된 반응력 시스템과 후속 모멘트를 분석했습니다. 알고리즘은 자동으로 하중을 적용하는 새 노드 위치를 제안하여 힘 시스템 내에서 0에 가까운 크기의 모멘트(moment of-zero)가 생성되도록 합니다. 이 반복 프로세스에서 진행 되며, 이를 통해 힘이 적용될 때 0에 가장 가까운 모멘트를 생성하는 노드 위치가 발견되거나 추정될 때까지 진행됩니다. 알고리즘은 토론 섹션에 자세히 설명되어 있습니다.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
절차 섹션(2단계)에 설명된 대로 세분화 및 수동 개요를 확인하기 위해, 마른 두개골에서 상악제 첫 어금니를 추출하고 CBCT 이미지를 촬영했습니다. 이미지 처리 및 편집 소프트웨어 Mimics는 2단계에서 설명한 바와 같이 치아를 수동으로 개략적으로 설명하는 데 사용되었습니다. 그 후 메시를 수행하고 세그먼트 모델을 3matic 소프트웨어로 세척하고 분석을 위해 Abaqus로 가져왔습니다. 우리는 실험실에서 측정된 치아의 FE 모델과 실제 치아에서 측정된 선형 및 체적 측정에서 유의한 차이를 발견하지못했습니다(보충 문서 4).
개체의 CRES를 결정하는 사용자 정의 알고리즘의 유효성을 확인하기 위해, 칼집 내에 둘러싸인 빔의 단순화된 모델이 스크립트 생성의 초기 단계에서 사용되었습니다(그림5A). 강철 케이스는 3자유도의 변위 자유도로 제한되었고 빔/시스 인터페이스의 노드는 함께 묶여 있었습니다. 강제 응용 프로그램에 대한 노드가 임의로 선택되었으며 솔루션이 수렴될 때까지 서브루틴이 반복방식으로 적용되었습니다. 단순화된 모델에서는 길이가 30단위이고 폭이 10단위로 칼집에 둘러싸여 있습니다. 정의된 알고리즘 및 그 계산을 따로 따로 따로 모델 빔의 CRES가 예측되었다(그림5B). 이것은 이론적 계산에 동의했습니다 (보충 문서 3참조). 따라서, 사용자 정의 알고리즘의 유효성은 이 단순화된 모델에서 개발 및 검증되었고, 이어서 상악치아의 CRES의 측정을 위해 구현되었다.
표 2는 구조물에 지정된 재료 특성을 보여줍니다. PDL과 뼈의 재료 특성 모델링의 차이는 치아의 CRES의 최종 위치에 영향을 줄 수 있습니다. 섬유 배향과 관련된 PDL 이방성, 푸아송 비율의 차이, 로딩 패턴 및 크기도 차이를 만들 수 있습니다. PDL은 Ogden 모델(μ1 = 0.07277, α1 = 16.95703, D1= 3 x10-7)22, 23에,23따라 비선형, 초탄성 특성을 할당하였다. 특정 밀도또한 뼈에 대해 = 1.85 g/cm3을 할당했습니다. 치아에 대한 2.02 g /cm3; 및 PDL (즉, 물의 밀도, PDL은 주로 물로 구성되어 있기 때문에)24,,25에대한 1 g /cm3.
힘 벡터를 표준화하고 CRES의위치를 찾기 위해, 카르테시안 좌표계를 구성(X-Y-Z)하고 다음 방향에 의해 정의하였다: Y축(전방 후방 또는 라비오링축)은 후방을 따라 양수 방향으로 후방 봉합사를 따라 배향하고, 수직방향(초조열또는 폐색지각축)의 Z축은 모델의 우수또는 치은부가 양수방향으로, X축은 양수방향으로 횡방향(buccolingual 축)을 가꾸고, X축은 양수방향으로 볼부를갖는다(도 6).
이 좌표계는 두 가지 방법으로 적용되었다: 1) X-Y 평면에서 교차 절개 및 어금니 간 폭을 양분하는 선에 위치한 절개 유두 아래 중앙 절개 의 얼굴 표면 사이에 위치한 기원 (O)와 함께 글로벌 좌표계가 설립되었다; 2) 로컬 좌표계는 모든 치아에 대해 원점으로 'R'으로 구성되었습니다. 모든 치아에 특이적인 'R' 포인트는 크라운의 볼 표면의 기하학적 중심으로 정의되었습니다. 이 사이트는 작업자가 치열 교정 힘을 적용하기 위해 대괄호를 배치할 수 있는 가장 가까운 위치를 근사화하도록 선택되었습니다. 대표적인 결과는 그림 7에나와 있습니다.
전역 및 로컬 좌표계를 기준으로 위치한 CRES는 표 3 및 표 4에나와 있습니다. Y 및 Z 좌표를 따라 힘 계가 적용될 때 X 좌표를 따라 얻어진 CRES의 위치는 서로달랐다(표 5). 그러나, 평균 차이는 작하였다(0.88±0.54 mm).
그림 1: 디자인 흐름도. CRES를찾기 위한 3단계 워크플로. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 2: 세 뷰(X-Y-Z)와 체적 모델에서 상악 치아를 표시하는 모방 소프트웨어의 레이아웃입니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 3: 3matic 소프트웨어의 비 매니폴드 어셈블리를 활용하여 치주 인대(PDL)를 생성하기 위한 단계. 리메쉬모듈(A)비매니폴드 어셈블리 생성,(B)맥실라는 주 엔티티로 설정되고,(C)PDL은 교차하는 엔티티로 설정되고,(D)적응형 리메시,(E)상악골과 PDL을 분할하고,(F)PDL에 대한 단계 B-F를 주 엔티티로,(G)메시를 생성한다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 4: Abaqus 소프트웨어 레이아웃입니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 5: 강철 빔의 단순화 모델. (A)정의된 알고리즘의 정확도를 테스트하는 데 사용되는 강철 외장에 둘러싸인 빔입니다. (B)정의된 알고리즘에 의해 예측된 대로 둘러싸인 빔의 CRES의 위치. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 6: 각 치아에 대한 글로벌 원점(O) 및 로컬 원점(R)을 기준으로 CRES 추정을 위한 좌표계입니다. 이것은 상악 두 번째 초구에 대한 그림입니다. 이 방법은 아치의 모든 치아에 사용되었습니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 7: 상악 치아의 CRES의 3차원 표현. (A)중앙 절개. (B)측면 앞니. (C)개. (D)첫 번째 스모글러. (E)두 번째 소구치. (F)첫 번째 어금니. (G)두 번째 어금니. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
| 중공 유형: | 둘 다 (외부 및 내부) |
| 거리 | 0.2 |
| 가장 작은 세부 사항: | 0.05 |
| 줄일: | 확인 |
| 테두리에서 정리: | 확인 |
| 정리 계수: | 1.1 |
표 1: 중공 공구 매개변수.
| 구조 | 탄성 계수 (MPa) | 푸아송 비율 | 특정 밀도 (g/ cm3) |
| 치아 | 17000 | 0.3 | 2.02 |
| 뼈 | 17000 | 0.3 | 1.85 |
| Pdl | 0.05 | 텍스트 보기 | 1 |
표 2: 유한 요소 모델의 재질 특성입니다.
| 톱니 번호 | 치아 길이 | 루트 길이 | Ⅹ | Y | Z |
| UL1 | 25.2 | 15.1 | 3.4 | 11.0 | 12.9 |
| UL2 | 26.0 | 16.8 | 8.8 | 13.2 | 14.3 |
| UL3 | 29.1 | 19.5 | 15.1 | 18.0 | 15.6 |
| UL4 | 23.8 | 15.7 | 18.4 | 21.5 | 10.6 |
| UL5 | 24.8 | 18.2 | 20.9 | 28.2 | 10.1 |
| UL6 | 22.0 | 16.4 | 25.8 | 38.7 | 11.6 |
| UL7 | 21.4 | 15.0 | 27.4 | 43.2 | 11.4 |
| UR1 | 24.9 | 14.6 | -4.6 | 10.8 | 13.2 |
| UR2 | 26.3 | 16.7 | -9.9 | 13.0 | 13.6 |
| UR3 | 30.9 | 21.1 | -15.6 | 17.7 | 14.2 |
| UR4 | 22.9 | 16.7 | -19.0 | 21.9 | 9.2 |
| UR5 | 23.4 | 16.7 | -21.1 | 29.4 | 8.8 |
| UR6 | 22.2 | 16.3 | -23.9 | 39.6 | 9.8 |
| UR7 | 20.8 | 15.9 | -21.7 | 47.0 | 10.4 |
표 3: 지구점 O와 관련하여 상악치아의 CRES의 3차원(X-Y-Z) 위치.
| 톱니 번호 | 치아 길이 | 루트 길이 | Ⅹ | Y | Z |
| UL1 | 25.2 | 15.1 | -1.1 | 10.9 | 9.4 |
| UL2 | 26.0 | 16.8 | -5.5 | 9.4 | 10.4 |
| UL3 | 29.1 | 19.5 | -5.7 | 9.3 | 13.2 |
| UL4 | 23.8 | 15.7 | -6.4 | 5.7 | 9.0 |
| UL5 | 24.8 | 18.2 | -6.7 | 7.0 | 9.5 |
| UL6 | 22.0 | 16.4 | -6.9 | 8.3 | 10.4 |
| UL7 | 21.4 | 15.0 | -8.6 | 3.3 | 7.3 |
| UR1 | 24.9 | 14.6 | 0.5 | 10.8 | 11.1 |
| UR2 | 26.3 | 16.7 | 5.0 | 10.3 | 9.3 |
| UR3 | 30.9 | 21.1 | 5.7 | 8.5 | 12.0 |
| UR4 | 22.9 | 16.7 | 5.3 | 5.3 | 9.3 |
| UR5 | 23.4 | 16.7 | 5.3 | 6.5 | 9.1 |
| UR6 | 22.2 | 16.3 | 5.6 | 7.8 | 10.1 |
| UR7 | 20.8 | 15.9 | 9.5 | 4.3 | 8.6 |
표 4: CRES가 평가되는 각 치아에 대한 국소 점 R과 관련하여 상악치아의 CRES의 3차원(X-Y-Z) 위치. 여기서 R은 크라운의 볼 표면의 기하학적 중심입니다.
| 톱니 번호 | 년도 | Fz | 차이 |
| UL1 | -1.36 | -0.80 | 0.56 |
| UL2 | -5.73 | -5.23 | 0.5 |
| UL3 | -6.00 | -5.45 | 0.55 |
| UL4 | -6.11 | -6.65 | 0.54 |
| UL5 | -5.95 | -7.40 | 1.46 |
| UL6 | -6.18 | -7.67 | 1.49 |
| UR1 | 0.36 | 0.67 | 0.31 |
| UR2 | 5.23 | 4.77 | 0.46 |
| UR3 | 5.93 | 5.38 | 0.55 |
| UR4 | 4.57 | 6.01 | 1.44 |
| UR5 | 5.88 | 4.69 | 1.91 |
| UR6 | 5.19 | 5.98 | 0.79 |
표 5: Y-(Fy) 및 Z(Fz) 축을 따라 힘이 가해지면 X축을 따라 저항 위치 중심의 변화.
보충 문서 1: FEA에 사용되는 알고리즘의 파이썬 스크립트입니다. 이 파일을 보려면 여기를 클릭하십시오 (다운로드 오른쪽 버튼을 클릭하십시오).
보충 문서 2: 힘 시스템의 해석 개요입니다. 이 파일을 보려면 여기를 클릭하십시오 (다운로드 오른쪽 버튼을 클릭하십시오).
보충 문서 3: 칼집에 싸인 단순 빔의 질량 중심의 이론적 추정. 이 파일을 보려면 여기를 클릭하십시오 (다운로드 오른쪽 버튼을 클릭하십시오).
보충 문서 4: 추출된 상악 첫 번째 어금니의 유한 요소 모델입니다. 이 파일을 보려면 여기를 클릭하십시오 (다운로드 오른쪽 버튼을 클릭하십시오).
Subscription Required. Please recommend JoVE to your librarian.
Discussion
본 연구는 환자의 CBCT 이미지로부터 유래된 상악치아 모델의 유한 요소 분석(FEA)에 대한 일관된 워크플로우를 확립하여 그들의 CRES를결정하는 도구 세트를 보여준다. 임상의의 경우 상악치아의 CRES의 명확하고 간단한 지도는 치아 의 움직임을 계획하고 부작용을 예측하는 귀중한 임상 도구가 될 것입니다. 유한원소법(FEM)은 1973년17년치과생체역학연구에 도입되었으며, 그 이후 폐포지지구조부6,,67,78,89,10,,911,,12의응력 및 변형분야를 분석하기 위해 적용되었다. 워크플로(그림1)에설명된 단계 수에 의해 입증된 바와 같이 유한 요소 모델을 만드는 것은 복잡한 작업입니다. 따라서 방법론의 특정 측면을 단순화해야 했습니다.
첫째, 폐포 소켓에서만 치아 운동이 폐포 뼈의 재흡수 및 apposition가 발생하지 않았다고 가정하여 고려되었습니다. 이러한 유형의 변위는1차 또는 순간치아운동(18)이라고한다. PDL은 즉각적인 치아 변위에서 중요한 개체라는 것이 관찰되었습니다. 뼈와 치아는 치아운동(15)에대한 PDL 응력을 정의하기 위해 경질이라고 합리적으로 가정할 수 있다. 따라서, 본 연구를 위해 응력 분포는 치아 소켓 내에 구속되었다. 경계 조건 작성 도구를 사용하면 모델에 대한 경계 조건을 설정하거나 구속조건을 적용할 수 있습니다. 선택한 점에는 0자유도가 지정되어 모델이 해당 영역에서 견고하게 유지되도록 합니다. 그 결과, 이전 연구에서 수행된 변형된 폐조골 뼈의 고체 원소를 골 변형 을 계산하고 재구성하기 위한 분석 시간이19,,20으로제거되었다.
둘째, 이미지 해상도를 적당한 수준으로 유지하려고 시도했습니다. CBCT 이미지 복셀 크기는 0.27 mm였습니다. 이것은 방사선 투여량을 최소한으로 유지했을 뿐만 아니라 적층 요소에 대한 글로벌 강성 매트릭스를 조립하기 위한 계산 부담을 감소시켰습니다. 그러나 단점은 CBCT 해상도가 스캔에서 PDL을 정확하고 뚜렷하게 캡처하기에 충분하지 않다는 것입니다. 이는 평균 PDL 두께가 약 0.15 mm-0.38 mm(평균: 0.2 mm)21이고 이미지 복셀 크기가 0.27 mm이기 때문입니다. CBCT 스캔의 이러한 단점은 두 가지 문제를 만들었습니다: 1) PDL을 자체적으로 분할할 수 없습니다. 및 2) 임계값을 사용하여 뼈와 치아를 분할하는 것은 둘 사이의 뚜렷한 회색 값 변화의 부족으로 인해 불가능하였다. 그 결과 회색 값이 비슷했기 때문에 소프트웨어가 치아와 뼈를 구분할 수 없었습니다. 즉, 모방은 치아와 뼈를 별도로 분할할 수 없었습니다. 따라서, 다른 세분화 방법이 개발되었다. 모방에서 영역 성장 또는 분할 도구와 같은 수많은 도구를 시도한 후, 치아를 분할하는 가장 좋은 방법은 CBCT의 각 슬라이스에 치아 구조를 수동으로 강조하여 결정되었습니다. 여기서 다중 슬라이스 편집 도구는 효율성 이점을 제공했습니다. 모든 슬라이스를 수동으로 강조 표시하는 대신 일부 조각만 강조 표시하면 됩니다. 이러한 이유로, 그것은 일관된 방식으로 치아의 해부학의 좋은 이미지를 얻기에 가장 큰 정확도를 제공으로, 치아를 분할하기위한 가장 좋은 방법이었다.
모방은 CBCT 이미지의 낮은 해상도로 인해 PDL을 분할 할 수 없었기 때문에 치아의 뿌리 구조로부터 PDL을 성장시킬 필요가 있었습니다. 이를 위해서는 CEJ에서 치아를 뿌리와 크라운으로 분할해야 했습니다. 일단 성장하면, 건설된 PDL은 본질적으로 서로 평행한 두 개의 표면이 0.2 mm 간격으로 떨어져 있었고, 한 표면은 뼈와 밀접하게 접촉하고 다른 표면은 뿌리와 밀접하게 접촉했다. 표면이 유한 요소 분석에서 함께 묶여 치아에 추가된 하중이 PDL을 통해 뼈로 전파되도록 하는 것이 중요했습니다. 엔지니어링 소프트웨어는 표면이 너무 멀리 떨어져 있거나 너무 많이 교차된 모델을 거부했는데, 이는 표면을 연결하는 것이 불가능하고 FEA 모델을 무효화했기 때문에 발생했습니다.
셋째, 모든 모델 표면은 비교적 매끄럽고 작은 표면 지형이 없어 전체 모델 해석에 대해 미미한 부분(예: 볼 피질 표면에서 추가 골격 투영)을 유지했습니다. 해부학의 투영에 미세 요소는 미세 해부학의 복잡한 영역에서 요소의 크기를 감소시켜 최종 모델의 메쉬에 불필요한 합병증을 추가, 따라서 모델의 요소의 수를 증가. 더 작고 많은 요소가 최종 유한 요소 해석에서 컴퓨팅 노력을 증가시다.
Y 및 Z 방향에 힘이 가해졌을 때의 CRES 위치는 X 방향을 따라 위치의 차이로 표현되어 서로 달랐습니다. 그러나, 그 차이는작음(표 5)과임상적으로 통계적으로 미미하였다. 따라서, 한 방향으로 계산된 CRES의 위치는 다른 방향으로 사용될 수 있다. 이전 연구는 또한 CRES에 대한 단일 포인트가10,,26,,27을관찰하지 않는 3D로 평가될 때 나타났다. 따라서 명확한 CRES를 갖는 대신 더 나은 용어가 "저항의 반경"이 될 수 있다고 제안되었습니다. 이러한 차이는 루트 형태, 경계 조건, 재료 특성 및 하중 적용점과 같은 여러 요인에 기인할 수 있습니다.
사용자 지정 알고리즘을 사용한 힘 시스템 분석
치아의 CRES를 찾기 위한 수학적 개념, 이론적 파생, 및 컴퓨터 시뮬레이션은 이전에 상세27,,28,,29,,30에기술되어 있다. 적용된 다양한 하중에 의해 생성된 힘 시스템을 분석하고 치아에 대한 CRES를 예측하기 위해 사용자 지정 알고리즘을 Abaqus 내에서 작성하고 실행했습니다(추가 코딩 파일 참조). 이 알고리즘은 Python을 사용하여 작성되었으며 FEA 소프트웨어 출력 데이터베이스 (.odb 파일)의 데이터를 입력으로 받아들이고 데이터를 처리하며 적용 된 부하에 의해 시스템에서 생성 된 순간에 대한 값을 제공합니다. 또한 노드 위치를 추정하여 시스템 내에서 더 낮은 모멘트 생성을 생성합니다. 이를 통해 사용자는 예측이 단일 위치에 수렴될 때까지 시뮬레이션을 반복적으로 실행할 수 있습니다.
알고리즘은 각 단계에서 적용된 하중의 결과로 절개 좌표, 각 노드의 총 변위 및 각 노드의 반응력에 액세스합니다. 원래 하중 적용과 동일한 방향의 반력과 반대 방향의 반력은 시뮬레이션 중에 치아에 작용하는 골재력 벡터를 결정하기 위해 시스템의 각 노드에서 합산됩니다. 결과 모멘트는 각 노드의 각 반력에 대한 힘 적용 점과 관련하여 계산되며 반력과 동일한 방식으로 합산됩니다. 따라서, 원래 하중 어플리케이션과 동일한 방향으로 골재 력 벡터가 생성되고 그 힘 벡터에 대한 힘 벡터에 의해 생성된 힘 도포는 반대 방향과 그 결과 모멘트에 대한 힘 벡터가 계산된다. 시스템이 정적 평형 상태이기 때문에 모든 힘과 모멘트의 합은 0과 같습니다. 그러나 이러한 방식으로 반력과 모멘트의 고장은 이러한 집계 힘이 시스템의 피벗 포인트역할을 하는 효과적인 위치를 계산할 수 있게 해주며, 이러한 피벗점 사이의 중심점은 CRES에더 가까운 힘 응용 프로그램의 근사치를 제공합니다.
이러한 계산을 수행하기 위해 결과 모멘트의 크기는 해당 힘의 크기로 분할되어 피벗 점에서 힘 응용 지점까지의 거리(R 벡터)의 크기를 지정합니다. R 벡터의 방향은 모멘트 및 힘 벡터의 교차 곱을 통해 결정되며, 여기서 모든 것은 서로 직교되어야 하며, 단위 벡터는 교차 생성물의 크기로 나누어 결정됩니다. 단위 벡터 R은 원래 힘의 점 응용 을 기준으로 각 피벗 점의 좌표의 3D 공간에서 전체 추정을 산출하기 위해 이전에 계산 된 R 벡터의 크기에 곱합니다. 이 두 벡터 사이의 중간점은 다음 반복에서 다음 힘 응용 프로그램의 위치에 대한 추정을 제공합니다. 추가 정보는 보충 문서 2에첨부되어 있습니다.
CRES의 추정은 시스템의 결과 모멘트가 약 0에 추가될 때 결정됩니다. 현재 연구의 경우, 이 결정은 계산된 모멘트의 가장 낮은 양수 및 음수 X-구성요소를 찾아 서 둘을 평균화하여 결정됩니다. 노드의 임의로 생성된 위치 및 임의의 두 노드(0.5 mm) 사이의 고유 거리로 인해, 정확한 제로 모멘트생성이 생성되는 위치를 찾기가어렵다(표 5).
제한
우리의 최선의 노력에도 불구하고,이 연구에 몇 가지 제한이 있습니다. 첫째, PDL은 CBCT상에서 가시화될 수 없기 때문에, 그 자체로 분별로 분할할 수 없었고 0.2 mm의 균일한 두께로 치아의 뿌리 표면에서 생성되었다.30,31 둘째, 정확한 모델을 작성하는 단계의 수는 길어졌습니다. 이것은 모델이 얼마나 빨리 만들어질 수 있는지에 대한 제한이며, 이는 사례별로 환자를위한 개인 치료 계획에 이러한 도구를 사용할 가능성을 제한합니다. 또한 이러한 모델을 생성하는 데 필요한 소프트웨어는 비용이 많이 들며 교육 기관이나 대기업에서 사용할 수 있는 리소스로 제한됩니다. 또한 모델이 만들어지면 FEA를 실행하려면 매우 강력한 컴퓨팅이 필요했습니다. 따라서, 이 방법은 필요한 기술이 널리 이용될 때까지 실행 가능한 치료 계획 도구가 될 수 없다.
미래 연구는 이 모형을 사용하여 상악치아에 유한 요소 분석을 수행하기 위하여 아치및 이의 단을 위한 CRES를 결정하기 위하여 집중해야 합니다, 특히 이의 그 단은 일반적으로 치열 교정에서 조작된, 추출 케이스에 있는 전방 세그먼트 또는 열린 바이트 환자에 침입을 위한 후방 세그먼트와 같은. 이러한 모델에 대해 CRES가 결정되면 기존 데이터에 추가할 추가 CBCT 이미지에서 추가 모델을 개발해야 합니다. CRES 위치의 충분한 데이터 풀을 통해 열 맵을 생성하여 임상의에게 귀중한 참조가 될 수 있는 CRES의 일반적인 위치를 나타낼 수 있습니다.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
저자는 공개 할 것이 없다.
Acknowledgments
저자는 프로젝트를 지원하기위한 찰스 버스톤 재단 상을 인정하고 싶습니다.
Materials
| Name | Company | Catalog Number | Comments |
| 3-matic software | Materialise, Leuven, Belgium. | Cleaning and meshing | |
| Abaqus/CAE software, version 2017 | Dassault Systèmes Simulia Corp., Johnston, RI, USA. | Finite Element Analysis | |
| Mimics software, version 17.0 | Materialise, Leuven, Belgium. | Segmentation of teeth and bone |
References
- Smith, R. J., Burstone, C. J. Mechanics of tooth movement. American Journal of Orthodontics. 85 (4), 294-307 (1984).
- Christiansen, R. L., Burstone, C. J. Centers of rotation within the periodontal space. American Journal of Orthodontics. 55 (4), 353-369 (1969).
- Tanne, K., Nagataki, T., Inoue, Y., Sakuda, M., Burstone, C. J. Patterns of initial tooth displacements associated with various root lengths and alveolar bone heights. American Journal of Orthodontics and Dentofacial Orthopedics. 100 (1), 66-71 (1991).
- Burstone, C. J., Pryputniewicz, R. J. Holographic determination of centers of rotation produced by orthodontic forces. American Journal of Orthodontics. 77 (4), 396-409 (1980).
- Dermaut, L. R., Kleutghen, J. P., De Clerck, H. J. Experimental determination of the Cres of the upper first molar in a macerated, dry human skull submitted to horizontal headgear traction. American Journal of Orthodontics and Dentofacial Orthopedics. 90 (1), 29-36 (1986).
- Tanne, K., Sakuda, M., Burstone, C. J. Three-dimensional finite element analysis for stress in the periodontal tissue by orthodontic forces. American Journal of Orthodontics and Dentofacial Orthopedics. 92 (6), 499-505 (1987).
- Meyer, B. N., Chen, J., Katona, T. R. Does the Cres depend on the direction of tooth movement? American Journal of Orthodontics and Dentofacial Orthopedics. 137 (3), 354-361 (2010).
- Kojima, Y., Fukui, H. A finite element simulation of initial movement, orthodontic movement, and the centre of resistance of the maxillary teeth connected with an archwire. European Journal of Orthodontics. 36 (3), 255-261 (2014).
- Reimann, S., Keilig, L., Jäger, A., Bourauel, C. Biomechanical finite-element investigation of the position of the centre of resistance of the upper incisors. European Journal of Orthodontics. 29 (3), 219-224 (2007).
- Viecilli, R. F., Budiman, A., Burstone, C. J. Axes of resistance for tooth movement: Does the Cres exist in 3-dimensional space? American Journal of Orthodontics and Dentofacial Orthopedics. 143 (2), 163-172 (2013).
- Ammar, H. H., Ngan, P., Crout, R. J., Mucino, V. H., Mukdadi, O. M. Three-dimensional modeling and finite element analysis in treatment planning for orthodontic tooth movement. American Journal of Orthodontics and Dentofacial Orthopedics. 139 (1), 59-71 (2011).
- Sia, S., Koga, Y., Yoshida, N. Determining the center of resistance of maxillary anterior teeth subjected to retraction forces in sliding mechanics. An in vivo study. Angle Orthodontics. 77 (6), 999-1003 (2007).
- Cattaneo, P. M., Dalstra, M., Melsen, B. Moment-to-force ratio, center of rotation, and force level: a finite element study predicting their interdependency for simulated orthodontic loading regimens. American Journal of Orthodontics and Dentofacial Orthopedics. 133 (5), 681-689 (2008).
- Tominaga, J. Y., et al. Effect of play between bracket and archwire on anterior tooth movement in sliding mechanics: A three-dimensional finite element study. Journal of Dental Biomechanics. 3, 1758736012461269 (2012).
- Cai, Y., Yang, X., He, B., Yao, J. Finite element method analysis of the periodontal ligament in mandibular canine movement with transparent tooth correction treatment. BMC Oral Health. 15 (106), (2015).
- Pauwels, R., Araki, K., Siewerdsen, J. H., Thongvigitmanee, S. S. Technical aspects of dental CBCT: state of the art. Dentomaxillofacial Radiology. 44 (1), 20140224 (2015).
- Farah, J. W., Craig, R. G., Sikarskie, D. L. Photoelastic and finite element stress analysis of a restored axisymmetric first molar. Journal of Biomechanics. 6 (5), 511-520 (1973).
- van Driel, W. D., van Leeuwen, E. J., Von den Hoff, J. W., Maltha, J. C., Kuijpers-Jagtman, A. M. Time-dependent mechanical behavior of the periodontal ligament. Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine. 214 (5), 497-504 (2000).
- Bourauel, C., et al. Simulation of orthodontic tooth movements. A comparison of numerical models. Journal of Orofacial Orthopedics. 60 (2), 136-151 (1999).
- Schneider, J., Geiger, M., Sander, F. G. Numerical experiments on longtime orthodontic tooth movement. American Journal of Orthodontics and Dentofacial Orthopedics. 121 (3), 257-265 (2002).
- Ten Cate, A. R. Oral histology, development, structure and function (5th ed). , St. Louis Mosby. (1998).
- McCormack, S. W., Witzel, U., Watson, P. J., Fagan, M. J., Gröning, F. The Biomechanical Function of Periodontal Ligament Fibres in Orthodontic Tooth Movement. PLoS One. 9 (7), e102387 (2014).
- Huang, H., Tang, W., Yan, B., Wu, B., Cao, D. Mechanical responses of the periodontal ligament based on an exponential hyperelastic model: a combined experimental and finite element method. Computer Methods in Biomechanics and Biomedical Engineering. 19 (2), 188-198 (2016).
- Yang, J. A new device for measuring density of jaw bones. Dentomaxillofacial Radiology. 31 (5), 313-316 (2002).
- Gradl, R., et al. Mass density measurement of mineralized tissue with grating-based X-ray phase tomography. PLoS One. 11 (12), e01677979 (2016).
- Jiang, F., Kula, K., Chen, J. Estimating the location of the center of resistance of canines. Angle Orthodontics. 86 (3), 365-371 (2016).
- Nyashin, Y., et al. Center of resistance and center of rotation of a tooth: experimental determination, computer simulation and the effect of tissue nonlinearity. Computer Methods in Biomechanics and Biomedical Engineering. 19, 229-239 (2016).
- Toms, S. R., Eberhardt, A. W. A nonlinear finite element analysis of the periodontal ligament under orthodontic tooth loading. American Journal of Orthodontics and Dentofacial Orthopedics. 123 (6), 657-665 (2003).
- Osipenko, M. A., Nyashin, M. Y., Nyashin, Y. I. Centre of resistance and centre of rotation of a tooth: the definitions, conditions of existence, properties. Russian Journal of Biomechanics. 3 (1), 5-15 (1999).
- Dathe, H., Nägerl, H., Dietmar, K. M. A caveat concerning center of resistance. Journal of Dental Biomechanics. 4, 1758736013499770 (2013).
- Hohmann, A., et al. Influence of different modeling strategies for the periodontal ligament on finite element simulation results. American Journal of Orthodontics and Dentofacial Orthopedics. 139 (6), 775-783 (2011).
