Hindawi

BioMed Research International

Volume 2017, Article ID 2827953, 7 pages
https://doi.org/10.1155/2017/2827953

Research Article
Influence of Deformation and Stress between Bone and Implant
from Various Bite Forces by Numerical Simulation Analysis

Hsin-Chung Cheng,1,2 Boe-Yu Peng,1,2 May-Show Chen,2,3 Chiung-Fang Huang,4,5

Yi Lin,6 and Yung-Kang Shen4
1
School of Dentistry, College of Oral Medicine, Taipei Medical University, Taipei, Taiwan
2
Department of Dentistry, Taipei Medical University Hospital, Taipei, Taiwan
3
School of Oral Hygiene, College of Oral Medicine, Taipei Medical University, Taipei, Taiwan
4
School of Dental Technology, College of Oral Medicine, Taipei Medical University, Taipei, Taiwan
5
Division of Family Dentistry, Department of Dentistry, Taipei Medical University Hospital, Taipei, Taiwan
6
Department of Business Administration, Takming University of Science and Technology, Taipei, Taiwan

Correspondence should be addressed to Yung-Kang Shen; ykshen@tmu.edu.tw

Received 18 January 2017; Accepted 27 March 2017; Published 28 May 2017

Academic Editor: Ren-Yeong Huang

Copyright © 2017 Hsin-Chung Cheng et al. This is an open access article distributed under the Creative Commons Attribution
License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly
cited.

Endosseous oral implant is applied for orthodontic anchorage in subjects with multiple tooth agenesis. Its effectiveness under
orthodontic loading has been demonstrated clinically and experimentally. This study investigates the deformation and stress on
the bone and implant for different bite forces by three-dimensional (3D) finite element (FE) methods. A numerical simulation of
deformation and stress distributions around implants was used to estimate the survival life for implants. The model was applied to
determine the pattern and distribution of deformations and stresses within the endosseous implant and on supporting tissues when
the endosseous implant is used for orthodontic anchorage. A threaded implant was placed in an edentulous segment of a human
mandible with cortical and cancellous bone. Analytical results demonstrate that maximum stresses were always located around the
implant neck in marginal bone. The results also reveal that the stress for oblique force has the maximum value followed by the
horizontal force; the vertical force causes the stress to have the minimum value between implant and bone. Thus, this area should
be preserved clinically to maintain the structure and function of a bone implant.

1. Introduction that maximum stresses happened at the location of the first

screw hole (closest to the condoyle) of an implant. The highest
The lifetime of a buried implant is evaluated by two phases microstrain was observed in bone adjacent to the first screw
as masticatory forces are applied: the unloaded healing phase hole.
and functional phase. Implants may fail during either phase, An applied mechanical force generates stress and strain in
typically for different reasons. Failure during the first phase bone, deforming its structural arrangement. Van Oosterwyck
happens within a short time after implant placement and is et al. [6] showed that a dehiscence can happen when a
associated primarily with inflammation [1]. Failure during narrow alveolar bone ridge is served with an oral implant.
the second phase takes place after implant loading and The presence of buccal and/or lingual dehiscence leads to
is associated mainly with direction of a load incorrectly a marked increase in marginal bone strain at the implant’s
oriented along the long axis of an implant [2, 3]. Implant mesial and distal sides. Faegh and Muftu [7] revealed that
size affects the area of possible retention in bone; factors the highest continuous interfacial stresses exist in the region
such as occlusion, masticatory force, number of implants, where an implant collar engages the cortical region and near
and implant position within prosthesis affect forces acting the implant apex in the subcortical region. Detolla et al. [8]
on bone adjacent to implants [4]. Kashi et al. [5] showed took advantage of computer-aided design (CAD) and FEA
2 BioMed Research International

to evaluate stresses on implant surfaces and in surrounding

bone. Menicucci et al. [9], who used two-dimensional (2D)
and three-dimensional (3D) finite element analysis (FEA) to
determine peri-implant stress occurring during tooth load-
ing, revealed that the static load was more harmful for peri-
implant bone than transitional load. Tada et al. [10] applied
three-dimensional finite element analysis to evaluate the
effect of implant type and length. The results demonstrated
that cancellous bone of higher rather than lower density
might ensure a better biomechanical environment of implant. CT image Integrated STL file
Geng et al. [11] reviewed the application of FEA in dental Figure 1: 3D model creation for numerical simulation.
implant studies. Based on clinical observations, some dentists
have declared that marginal bone loss around the implant
neck is approximately 0.5–1 mm or as high as 1.5 mm during
the first year after implant loading. Subsequently, the loss rate
of bone is considered either stable or significantly reduced
(bone loss of approximately 0.1 mm), or resort of the bone First molar
crest continues and the implant is ultimately lost. These find-
ings are in accordance with those gained 3D mathematical
models of dental implants under nonaxisymmetrical loading,
demonstrating that maximum stress happened around the
implant neck.
This study compared the effect of various bite forces (ver-
tical force, horizontal force, and oblique force) on deforma-
tion and stress distributions around the implant. The vertical
force is the normal chewing ability (masticatory force). The Figure 2: Mandible model for numerical simulation.
horizontal force means the molar force when human being
is sleeping. The oblique force is as the human being eats the
irregular food (special masticatory force). Thus, single cylin-
drical dental implant placed vertically into the molar region selected based on those of the implant most frequently used
of the mandible was modeled using 3D graphics. Deforma- in Taiwan (data were obtained from implant distributors).
tion and stress distribution in the bone socket after loading The mesial and distal borders of the end of the modeled
by bite force in different directions was computed by FEA. mandible section were restricted, such that displacement
of nodes in all directions was zero. The whole bone was
2. Numerical Method supposed as a homogenous and isotropic material with the
characteristics of cancellous bone. This was realized as spongy
A computer tomography (CT) scanner of dentistry (i-CAT bone that changed its structure after successful implant
3D; Imaging Sciences International, Inc., USA) was applied to placement and the implant interface changes into increasing
scan a patient’s mandible. The patient was a 35-year-old male similar to cortical bone. The tooth is zirconia material. The
recruited from Taipei Medical University Hospital, Taiwan. most frequently used implants in Taiwan are 4.5 mm in
Mill and Amira V3.1.1 software program was taken advantage diameter, 0.7 mm in pitch, and 8 mm in length. Relying on
of to integrate the dental CT scan and 3D model. Finally, this implant size, the models consisted of 27855 elements and
study employs the Solidworks software to get the STL file of 47646 nodes (Figure 4).
the 3D model. This model can be applied as the initial model Geometric three-dimensional models of the implant
for numerical simulation using ANSYS software. Figure 1 and molar region of the mandible and material properties
shows the process of this schedule. Finite element analysis (previously mentioned) of bone were simplified to decrease
is commonly applied for deformation and stress analysis of computing time and memory consumption. The authors
nonrigid bodies. The area of interest, the model (an implant suppose these geometric modifications do not influence
and part of a mandible list in Figure 2), is divided into element computational accuracy in terms of deformation and stress
meshes. The material properties of implant system bone are distributions. The simplifications applied in this study did not
shown in Table 1. Figure 3 shows the model of implant system affect analytical results, because all models had undergone
and cancellous bone. The physical properties of each element the same simplification. Model simplification has a basis;
are fixed. Boundary conditions, restrictions to physical prop- important factors are taken into account; so the differences
erties stemming from the physical behavior of the patient, with the real results are not.
are used for outer model elements. Deformation or stress, for To yield the 3D model, pre- and postprocessor ANSYS
example, was determined within each element in the model. computer-aided engineering (CAE) in FEA software (ANSYS
In this study, FEA was used to analyze deformations and V12) was applied. This pre- and postprocessor offers a para-
stresses created by different bite forces around cylindrical metric definition of geometry and the FE mesh (ANSYS,
dental implants. Dimensions of the reference implant were V12).
BioMed Research International 3

Solidworks before the combination

Artificial
ceramic tooth

Abutment

Implant

Cancellous bone

(a) (b)

(c)
Figure 3: Implant system: (a) part, (b) assembly, and (c) cross section.

Table 1: Material properties of implant system.

Tooth (zirconia) Abutment (titanium) Implant (titanium) Cancellous bone

Young’s modules (Pa) 2 × 1011 9.6 × 1010 9.6 × 1010 1.37 × 109
Poisson’s ratio 0.3 0.3 0.3 0.3

Loading of implants in three dimensions with various bite The authors had done the test of element quality; element
forces (masticatory force (vertical force), molar force (hori- quality is mostly close to 1. So the results are credible. The
zontal force), and special masticatory force (oblique force)) deformation and stress (Pa) at the implant-bone interface
(Figure 5), the values of the bite forces used are 500 N, 1000 N, were computed using FEA software.
1500 N, 2000 N, and 2500 N. The special masticatory force
in different angles (30∘ , 45∘ , and 60∘ ) is used relative to 3. Results and Discussion
the occlusal plane. This 3D loading acted on the center of
the upper surface of the tooth. Force magnitudes, as well Figure 6 shows the deformation distributions of various
as the acting point, were chosen based on previous work masticatory forces (500 N–1000 N). Figures 6(a) and 6(b)
at Taipei Medical University Hospital. Assuming that the show that the maximum deformation value occurs on the top
implant binds to the bone, the interfacial condition is bonded of the tooth. Deformation value decreases as the deformation
[12]. Pre- and postprocessing were conducted on a personal position goes to the bottom of implant. Figures 6(c) and
computer (PC). The CPU of PC had an Intel type (Core 6(d) indicate that the maximum deformation value occurs
2 Q8200) and 3.5 GB RAM. Computation time for each on the top of the tooth. Deformation value decreases as
simplified single-size implant was approximately one and half the deformation position goes to the implant abutment. The
hours. All computations were fulfilled for the 3D models. deformation distribution of implant presents the symmetry
4 BioMed Research International

(a)

0 0.0025 0.005

0.0013 0.0037
(m)

(b)

Figure 4: Mesh and screw type of implant system: (a) mesh and (b) screw type (diameter = 4.5 mm, pitch = 0.7 mm).

휃 = 30∘ , 45∘ , 60∘

(a) (b) (c)

Figure 5: Various bite forces: (a) vertical force, (b) horizontal force, and (c) oblique force.
BioMed Research International 5

(a) (b) (c)

(d) (e) (f)

(g) (h)
Figure 6: Deformation distributions for various forces: (a) 500 N, vertical force; (b) 1000 N, vertical force; (c) 500 N, horizontal force; (d)
1000 N, horizontal force; (e) 500 N, oblique force (45∘ ); (f) 1000 N, oblique force (45∘ ); (g) 500 N, oblique force (60∘ ); (h) 1000 N, oblique force
(60∘ ).

by means of the central axis of the tooth. The deformation tooth. The deformation distribution of implant system by
distribution of implant system by horizontal force has the oblique force (60∘ ) has the difference among it by vertical
difference among it by vertical force. Figures 6(e) and 6(f) force. Figure 7 shows the stress distributions of various mas-
reveal that the maximum deformation value takes place on ticatory forces (500 N–1000 N). A study goal was to identify
the top of the tooth. Deformation value decreases as the the exact location of the implant abutment. Analytical results
deformation position goes to the implant abutment. The for Figures 7(a) and 7(b) demonstrate that maximum stress
deformation distribution by oblique force (45∘ ) shows the was on the neck of the implant abutment. Additionally, the
symmetry by means of the central axis of the tooth. The neck of the implant abutment was the weak point. Analytical
deformation distribution of implant system by vertical force results for Figures 7(c) and 7(d) demonstrate that maximum
demonstrates the oblique shape on the tooth. Figures 6(g) stress was on the neck of the implant abutment. Additionally,
and 6(h) show that the maximum deformation value occurs the neck of the implant abutment was the weak point. The
on the top of the tooth. Deformation value decreases as stress distribution of the implant system decreases with
the deformation position goes to the implant abutment. The the central axis of the implant. The stress distribution of
deformation distribution of implant by oblique force (60∘ ) implant system also presents decrease from neck of implant
presents the symmetry by means of the central axis of the abutment to the top of the tooth. Analytical results for
6 BioMed Research International

(a) (b) (c)

(d) (e) (f)

(g) (h)
Figure 7: Stress distributions for various forces: (a) 500 N, vertical force; (b) 1000 N, vertical force; (c) 500 N, horizontal force; (d) 1000 N,
horizontal force; (e) 500 N, oblique force (45∘ ); (f) 1000 N, oblique force (45∘ ); (g) 500 N, oblique force (60∘ ); (h) 1000 N, oblique force (60∘ ).

Figures 7(e) and 7(f) demonstrate that maximum stress was Figure 8(a) reveals the maximum deformation of implant
on the neck of the implant abutment. Additionally, the neck system on different bite forces (various values). The deforma-
of the implant abutment was the weak point. The stress tion of the implant system has the maximum value by hori-
distribution of the implant system decreases with the central zontal force; then it is by oblique force 30∘ , 45∘ , and 60∘ , and
axis of the implant. The stress distribution of implant system it by vertical force earns the minimum value. The maximum
also presents decrease from neck of implant abutment to value of deformation increases as the bite force increases.
the top of the tooth. Analytical results for Figures 7(g) and Figure 8(b) reveals the maximum stress of implant system
7(h) indicate that maximum stress was on the abutment. A on different bite forces (various values). The maximum stress
study goal was to identify the exact location of the implant value of the implant system occurs by horizontal force; then
abutment. Analytical results demonstrate that maximum it is by oblique force 30∘ , 45∘ , and 60∘ ; and stress by vertical
stress was on the neck of the implant abutment. Additionally, force earns the minimum value. The reason is that horizontal
the neck of the implant abutment was the weak point. The force can generate a larger moment than other loading on
stress distribution of the implant system decreases with the the neck of the implant abutment. Vertical force acts directly
central axis of the implant. The stress distribution of implant on the abutment of implant. This phenomenon gets smaller
system also presents decrease from neck of implant abutment moment on the implant abutment. The maximum value of
to the top of the tooth. stress increases as the bite force increases.
BioMed Research International 7

250.0 system is the weak point as the bite forces come from the
human being.
200.0
Deformation (∗10 −3 mm)

Conflicts of Interest
150.0
The authors declare that they have no conflicts of interest.

100.0
References
[1] A. Sethi, T. Kaus, and P. Scohor, “The use of angulated
50.0
abutments in implant dentistry: five year clinical results of an
ongoing prospective study,” International Journal of Oral and
0.0 Maxillofacial Implants, vol. 15, no. 6, pp. 801–810, 2000.
500 1000 1500 2000 2500 [2] I. Akpinar, F. Demirel, L. Parnas, and S. Sahin, “A comparison
Force (N) of stress and strain distribution characteristics of two different
rigid implant designs for distal-extension fixed prostheses,”
Horizontal forces Oblique forces (60∘ )
Quintessence International, vol. 27, no. 1, pp. 11–17, 1996.
Oblique forces (30∘ ) Vertical forces
Oblique forces (45∘ ) [3] H. W. Anselm Wiskott and U. C. Belser, “Lack of integration
of smooth titanium surfaces: a working hypothesis based
(a) on strains generated in the surrounding bone,” Clinical Oral
Implants Research, vol. 10, no. 6, pp. 429–444, 1999.
6.E + 09
[4] D. R. Carter, M. C. Van der Meulen, and G. S. Beaupre,
“Mechanical factors in bone growth and development,” Bone,
5.E + 09
vol. 18, no. 1, supplement 1, pp. S5–S10, 1996.
[5] A. Kashi, A. R. Chowdhury, and S. Saha, “Finite element analysis
4.E + 09
of a TMJ implant,” Journal of Dental Research, vol. 89, no. 3, pp.
Stress (Pa)

241–245, 2010.
3.E + 09 [6] H. Van Oosterwyck, J. Duyck, J. Vander Sloten, G. Van Der
Perre, and I. Naert, “Peri-implant bone tissue strains in cases
2.E + 09 of dehiscence: a finite element study,” Clinical Oral Implants
Research, vol. 13, no. 3, pp. 327–333, 2002.
1.E + 09 [7] S. Faegh and S. Muftu, “Load transfer along the bone-dental
implant interface,” Journal of Biomechanics, vol. 43, no. 9, pp.
0.E + 00 1761–1770, 2010.
500 1000 1500 2000 2500 [8] D. H. Detolla, S. Andreana, A. Patra, R. Buhite, and B. Comella,
Force (N) “Role of the finite element model in dental implants,” The
Journal of Oral Implantology, vol. 26, no. 2, pp. 77–81, 2000.
Horizontal forces Oblique forces (60∘ )
Oblique forces (30∘ ) Vertical forces [9] G. Menicucci, A. Mossolov, M. Mozzati, M. Lorenzetti, and
Oblique forces (45∘ ) G. Preti, “Tooth-implant connection: some biomechanical
aspects based on finite element analyses,” Clinical Oral Implants
(b) Research, vol. 13, no. 3, pp. 334–341, 2002.
Figure 8: Maximum deformation and stress value on various bite [10] S. Tada, R. Stegaroiu, E. Kitamura, O. Miyakawa, and H.
forces: (a) deformation and (b) stress. Kusakari, “Influence of implant design and bone quality
on stress/strain distribution in bone around implants: a 3-
dimensional finite element analysis,” International Journal of
Oral and Maxillofacial Implants, vol. 18, pp. 357–368, 2003.
4. Conclusions [11] J. P. Geng, K. B. Tan, and G. R. Liu, “Application of finite element
analysis in implant dentistry: a review of the literature,” The
The deformation of the implant system has the maximum Journal of Prosthetic Dentistry, vol. 85, pp. 585–598, 2001.
value by horizontal force, then it is by oblique force 30∘ , [12] W. Xu, J. H. Wang, J. Geng, and H. M. Huang, “Application of
45∘ , and 60∘ ; and it by vertical force gets the minimum commercial FEA software,” in Application of the Finite Element
value. The maximum stress value of the implant system Method in Implant Density, pp. 92–126, Springer, 2008.
occurs by horizontal force; then it is by oblique force 30∘ ,
45∘ , and 60∘ and stress by vertical force gets the minimum
value. The deformation and stress distribution by molar
force (horizontal force) show the maximum value for dental
implant system. The molar force for the patient is key factor
after dental implantation clinic. The results also show that
the distribution of deformation and stress increases as the
bite force increases. This study shows the neck of the implant
 MEDIATORS of

INFLAMMATION

The Scientific Gastroenterology Journal of

World Journal
Hindawi Publishing Corporation
Research and Practice
Hindawi Publishing Corporation
Hindawi Publishing Corporation
Diabetes Research
Hindawi Publishing Corporation
Disease Markers
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

Journal of International Journal of

Immunology Research
Hindawi Publishing Corporation
Endocrinology
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

Submit your manuscripts at

https://www.hindawi.com

BioMed
PPAR Research
Hindawi Publishing Corporation
Research International
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

Journal of
Obesity

Evidence-Based
Journal of Stem Cells Complementary and Journal of
Ophthalmology
Hindawi Publishing Corporation
International
Hindawi Publishing Corporation
Alternative Medicine
Hindawi Publishing Corporation Hindawi Publishing Corporation
Oncology
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

Parkinson’s
Disease

Computational and
Mathematical Methods
in Medicine
Behavioural
Neurology
AIDS
Research and Treatment
Oxidative Medicine and
Cellular Longevity
Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014