{"id":289327,"date":"2023-06-19T08:41:36","date_gmt":"2023-06-19T08:41:36","guid":{"rendered":"\/forum\/forums\/topic\/no-contact-penetration-of-conta178-elements-bearing-simulation\/"},"modified":"2023-06-19T08:41:36","modified_gmt":"2023-06-19T08:41:36","slug":"no-contact-penetration-of-conta178-elements-bearing-simulation","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/no-contact-penetration-of-conta178-elements-bearing-simulation\/","title":{"rendered":"No contact penetration of CONTA178 elements – Bearing simulation"},"content":{"rendered":"Hello all,<\/p>\n
I would like to calculate the contact penetration of a tapered roller bearing by using CONTA178 elements between the raceways. So that, I´ve created a model with 30 contact elements including a pre-defined gap along the surface of the raceway and applied a constant pressure around the surface of the outer ring to simulate the pre-load of the bearing. The gap should simulate the profile of the logarithmic raceways.
However, while solving the model I didn´t get any connection between the outer ring and rollers. The CONTA178 elements are still open and there isn´t any contact penetration. I´m not sure whats the reason for that, I tired to vary the boundary conditions a bit, but that doesn´t helped.<\/p>\n
I would be very happy if someone here can help me.<\/p>\n
Thank you.<\/p>\n
import matplotlib.pyplot as plt
import numpy as np
import math as m
import itertools
import pandas<\/p>\n
from ansys.mapdl.core import launch_mapdl<\/p>\n
mapdl = launch_mapdl()<\/p>\n
mapdl.prep7()<\/p>\n
mapdl.units(“SI”)<\/p>\n
mapdl.et(1, “SOLID185”)
mapdl.mp(“EX”, 1, 210e9)
mapdl.mp(“DENS”, 1, 7800)
mapdl.mp(“NUXY”, 1, 0.3)
mapdl.type(1)
mapdl.keyopt(1, 2, 3)<\/p>\n
mapdl.et(2, “CONTA178”)
mapdl.mp(“EX”, 2, 210e9) # Elastic moduli in Pa (kg\/(m*s**2))
mapdl.mp(“DENS”, 2, 7800) # Density in kg\/m3
mapdl.mp(“NUXY”, 2, 0.3) # Poisson’s Ratio
mapdl.type(2)
mapdl.keyopt(2, 3, 1)
mapdl.keyopt(2, 4, 1)
mapdl.keyopt(2, 5, 0)
mapdl.keyopt(2, 10, 0)
mapdl.keyopt(2, 12, 1)<\/p>\n
# Slicing
anzahl_scheiben = 30
anzahl_scheiben_board = 20<\/p>\n
# Geometry
anzahl_rollen = 16<\/p>\n
# Inner ring<\/p>\n
di = 0.065 #Innen-Durchmesser Lagersitz
bi = 0.041 #Breite Lagersitz
α = 7.6 #Druckwinkel
β = 9.6 #Hilfswinkel<\/p>\n
o = 0.009238
p = 0.008282
q = 0.003962
r = 0.000671
s = 0.02539
t = 0.003391
v = 0.005083
w = 0.00079 #0.00086
b = 0.006658<\/p>\n
# Inner ring – Keypoints<\/p>\n
mapdl.k(1, 0, di\/2, 0)
mapdl.k(2, 0, di\/2+o, 0)
mapdl.k(3, 0, di\/2+o, p)
mapdl.k(4, 0, di\/2+o-q, p+r)
mapdl.k(5, 0, di\/2+o-q+t, p+r+s)
mapdl.k(6, 0, di\/2+o-q+t+v, p+r+s-w)
mapdl.k(7, 0, di\/2+o-q+t+v, p+r+s+b)
mapdl.k(8, 0, di\/2, bi)<\/p>\n
mapdl.k(13, 0, 0, 0)
mapdl.k(14, 0, 0, -0.005)
mapdl.k(15, 0, 0, 0.005)<\/p>\n
# Inner ring – Lines<\/p>\n
li_1 = mapdl.l(1, 2)
mapdl.l(2, 3)
mapdl.l(3, 4)
l1_i = mapdl.l(4, 5)
l1_b = mapdl.l(5, 6)
mapdl.l(6, 7)
mapdl.l(7, 8)
li_8 = mapdl.l(8, 1)<\/p>\n
# Inner ring – Surface<\/p>\n
a1 = mapdl.a(1, 2, 3, 4, 5, 6, 7, 8)<\/p>\n
mapdl.lesize(nl1=li_1, ndiv=3) #Boundary Conditions, Definition der Anzahl von Elementen
mapdl.lesize(nl1=li_8, ndiv=10)<\/p>\n
# Outer ring<\/p>\n
ba = 0.032 #Breite Lagersitz
c = 0.00732
e = 0.010942
f1 = 0.032
g = 0.00439
da = di+2*(o+c+e) #Außen-Durchmesser <\/p>\n
# Outer ring – Keypoints<\/p>\n
mapdl.k(9, 0, di\/2+o+c, 0)
mapdl.k(10, 0, da\/2, 0)
mapdl.k(11, 0, da\/2, f1)
mapdl.k(12, 0, da\/2-g, f1)
mapdl.k(16, 0, da\/2-0.00439, f1-0.000102) # Kontaktpunkte Kegelrolle, evtl. über zyl. KOS bestimmen, statt über Winkelfunktion auszurechnen!
mapdl.k(17, 0, da\/2-0.00955, f1-0.025245) # Kontaktpunkte Kegelrolle, evtl. über zyl. KOS bestimmen, statt über Winkelfunktion auszurechnen!<\/p>\n
# Outer ring – Lines<\/p>\n
mapdl.l(9, 10)
l_a_f = mapdl.l(10, 11)
mapdl.l(11, 12)
mapdl.l(12, 16)
l1_a = mapdl.l(16, 17)
mapdl.l(17, 9)<\/p>\n
mapdl.lesize(nl1=l1_a, ndiv=anzahl_scheiben-1) #Scheibenmodell Außenring-Lauffläche
mapdl.lesize(nl1=l1_i, ndiv=anzahl_scheiben-1) #Scheibenmodell Innenring-Lauffläche
mapdl.lesize(nl1=l1_b, ndiv=anzahl_scheiben_board-1) #Scheibenmodell Innenring-Bordkontakt<\/p>\n
mapdl.lesize(nl1=l_a_f, ndiv=10)<\/p>\n
# Outer ring – Surface<\/p>\n
a2 = mapdl.a(9, 10, 11, 12, 16, 17)<\/p>\n
# Extrudieren – Innenring & Außenring<\/p>\n
mapdl.vrotat(a1, pax1=14, pax2=15, arc=360, nseg=anzahl_rollen) #Hier bei nseg noch Anzahl Rollen einfügen
mapdl.vrotat(a2, pax1=14, pax2=15, arc=360, nseg=anzahl_rollen)<\/p>\n
# Roller<\/p>\n
f2 = 0.044114
a1 = 0.004275
a2 = 0.00788
m1 = 0.025241
h = 0.006887
a3 = 0.005862
a4 = 0.006707
a5 = 0.001135
a6 = 0.006322<\/p>\n
# Roller – Keypoints<\/p>\n
mapdl.k(2000000, 0, f2, a2)
mapdl.k(2100000, 0, f2+a1, a2+m1)
mapdl.k(3000000, 0, f2+a3, h)
mapdl.k(3100000, 0, f2+a1+a4, a2+m1-a5)<\/p>\n
# Roller – Lines<\/p>\n
l1_k = mapdl.l(2000000, 3000000)
l2_k = mapdl.l(3000000, 3100000)
mapdl.l(3100000, 2100000)
l3_k = mapdl.l(2000000, 2100000)<\/p>\n
# Roller – Surface<\/p>\n
f_a3 = mapdl.a(2000000, 3000000, 3100000, 2100000)<\/p>\n
mapdl.lesize(nl1=l2_k, ndiv=anzahl_scheiben-1)<\/p>\n
# Extrude – Roller<\/p>\n
mapdl.vsel(“NONE”)
k_rot = mapdl.vrotat(f_a3, pax1=2000000, pax2=2100000, arc=360, nseg=2)<\/p>\n
mapdl.csys(1)
mapdl.vgen(itime=anzahl_rollen, nv1=”ALL”, dy=360\/anzahl_rollen, imove=0)
mapdl.vsel(“ALL”)<\/p>\n
# Meshing<\/p>\n
mapdl.csys(0)
mapdl.vsweep(“ALL”)<\/p>\n
# Creation of contact elements
# Selection of lines<\/p>\n
mapdl.lsel(“S”, “NDIV”, vmin=anzahl_scheiben-1)
#print(mapdl.geometry.lnum)
mapdl.lplot(vtk=True, show_line_numbering=True, background=”white”, color_lines=”black”)<\/p>\n
polar_rollen = 360\/anzahl_rollen
radius_ir = di\/2+o-q+0.001
radius_ri = di\/2+o-q+0.002
radius_ra = f2+a3+0.001
radius_ar = f2+a6+0.0025<\/p>\n
x_ir = []
x_ri = []
x_ra = []
x_ar = []<\/p>\n
y_ir = []
y_ri = []
y_ra = []
y_ar = []<\/p>\n
for i in np.arange(0, 360+polar_rollen, polar_rollen):
x_ir.append(radius_ir*np.sin(np.radians(i)))
x_ri.append(radius_ri*np.sin(np.radians(i)))
x_ra.append(radius_ra*np.sin(np.radians(i)))
x_ar.append(radius_ar*np.sin(np.radians(i)))
y_ir.append(radius_ir*np.cos(np.radians(i)))
y_ri.append(radius_ri*np.cos(np.radians(i)))
y_ra.append(radius_ra*np.cos(np.radians(i)))
y_ar.append(radius_ar*np.cos(np.radians(i)))
<\/p>\n
### Selection of Nodes on selected Lines between future contact (Not shown in this post) ###<\/p>\n
# Contact definition
mapdl.vsel(“ALL”)<\/p>\n
# Inner ring – Roller<\/p>\n
profil_ir = [0.000012544, 0.000010816, 0.000009216, 0.000007744, 0.0000064, 0.000005184, 0.000004096, 0.000003136, 0.000002304, 0.0000016, 0.000001024, 0.000000576, 2.56E-07, 6.4E-08, 0, 6.4E-08, 2.56E-07, 5.76E-07, 0.000001024, 0.0000016, 2.304E-06, 3.136E-06, 0.000004096, 0.000005184, 0.0000064, 0.000007744, 9.216E-06, 0.000010816, 0.000012544, 0.0000144]<\/p>\n
ri_nnum_flip = np.roll(ri_nnum, anzahl_scheiben-1, axis=1)<\/p>\n
mapdl.et(2, “CONTA178”)
mapdl.type(2)
for a in range(0, anzahl_rollen, 1):
for i in range(2, anzahl_scheiben+2, 1):
mapdl.r(i, r1=-210000, r2=profil_ir[i-2])
mapdl.real(i)
mapdl.e(ir_nnum[a][i-2], ri_nnum_flip[a][i-2])<\/p>\n
# Outer ring – Roller<\/p>\n
profil_ar = [0.000012544, 0.000010816, 0.000009216, 0.000007744, 0.0000064, 0.000005184, 0.000004096, 0.000003136, 0.000002304, 0.0000016, 0.000001024, 0.000000576, 2.56E-07, 6.4E-08, 0, 6.4E-08, 2.56E-07, 5.76E-07, 0.000001024, 0.0000016, 2.304E-06, 3.136E-06, 0.000004096, 0.000005184, 0.0000064, 0.000007744, 9.216E-06, 0.000010816, 0.000012544, 0.0000144]<\/p>\n
ar_nnum_flip = np.flip(ar_nnum, axis=1)
ra_nnum_roll = np.roll(ra_nnum, anzahl_scheiben-1, axis=1)<\/p>\n
mapdl.et(2, “CONTA178”)
mapdl.type(2)
for a in range(0, anzahl_rollen, 1):
for i in range(2, anzahl_scheiben+2, 1):
mapdl.r(i, r1=-210000, r2=profil_ar[i-2])
mapdl.real(i)
mapdl.e(ar_nnum_flip[a][i-2], ra_nnum_roll[a][i-2])<\/p>\n
mapdl.eplot(vtk=True, show_edges=True, show_axes=True, line_width=1, show_node_numbering=False)<\/p>\n
# Boundary conditions
# Innenring
# radial
mapdl.allsel()
mapdl.csys(1)
mapdl.nsel(“S”, “LOC”, “X”, vmin=di\/2-0.001, vmax=di\/2+0.001)
mapdl.d(“ALL”, “UX”)
mapdl.d(“ALL”, “UY”)<\/p>\n
# axial
mapdl.allsel()
mapdl.csys(0)
mapdl.asel(“S”, “LOC”, “Z”, vmin=p+r+s+b-0.001, vmax=p+r+s+b+0.001)
mapdl.nsla(“S”, 1)
mapdl.d(“ALL”, “UZ”)<\/p>\n
# Outer ring
# axial
mapdl.allsel()
mapdl.csys(0)
mapdl.asel(“S”, “LOC”, “Z”, vmin=0, vmax=0.001)
mapdl.csys(1)
mapdl.asel(“R”, “LOC”, “X”, vmin=di\/2+o+c, vmax=di\/2+o+c+e)
mapdl.nsla(“S”, 1)
mapdl.csys(0)
mapdl.d(“ALL”, “UZ”)<\/p>\n
# Roller
mapdl.allsel()
mapdl.csys(0)
mapdl.vsel(“S”, “LOC”, “Z”, vmin=p+r+s\/2-0.0015, vmax=p+r+s\/2+0.0015)
mapdl.csys(1)
mapdl.vsel(“S”, “LOC”, “X”, vmin=f2, vmax=f2+0.005)
mapdl.nslv(“S”, 1)
mapdl.csys(1)
mapdl.d(“ALL”, “UX”)
mapdl.d(“ALL”, “UZ”)<\/p>\n
# Forces
mapdl.allsel()
mapdl.csys(1)
mapdl.asel(“S”, “LOC”, “X”, vmin=di\/2+o+c+e-0.001, vmax=di\/2+o+c+e+0.001)
mapdl.sfa(“ALL”, 1, “PRES”, value=2000)<\/p>\n
# Solve
mapdl.run(“\/SOL”)
mapdl.antype(0)
mapdl.outres(“ALL”, “ALL”)
mapdl.run(“\/STATUS,SOLU”)
sol_output = mapdl.solve()
mapdl.finish()
print(sol_output)<\/p>\n
# Plot results
result = mapdl.result
result.plot_principal_nodal_stress(
0,
“SEQV”,
lighting=False,
background=”w”,
show_edges=True,
text_color=”k”,
add_text=False,)<\/p>\n
nnum, stress = result.principal_nodal_stress(0)
von_mises = stress[:, -1] <\/p>\n
max_stress = np.nanmax(von_mises)<\/p>\n
# Contact Penetration
mapdl.allsel()
mapdl.lsel(“S”, “LINE”, vmin=4)
mapdl.nsll(“S”, 1)
mapdl.esln(“S”, 0)
mapdl.esel(“R”, “ENAM”, “”, 178)
mapdl.eplot(vtk=True)
contact_penetration_r = mapdl.post_processing.element_values(item=’nmisc’, comp=3)
print(contact_penetration_r)<\/p>\n
mapdl.post1()
mapdl.set(“last”)
mapdl.allsel()<\/p>\n
mapdl.post_processing.plot_nodal_principal_stress(“1”, smooth_shading=False)<\/p>\n
mapdl.exit()<\/p>\n
![\"v.Mises](\"\/forum\/wp-content\/uploads\/sites\/2\/2023\/06\/19-06-2023-1687164023-PyMAPDL_Bearing.PNG\" "\\"v.Mises")
<\/p>\n","protected":false},"template":"","class_list":["post-289327","topic","type-topic","status-closed","hentry","topic-tag-bearing","topic-tag-CONTA-1","topic-tag-contact","topic-tag-contact-method","topic-tag-mapdl","topic-tag-mechanical-apdl","topic-tag-PyMAPDL-1","topic-tag-python"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_subscription":["255882","88","18229"],"_bbp_author_ip":["96.7.74.173"]," _bbp_last_reply_id":["0"]," _bbp_likes_count":["0"],"_btv_view_count":["913"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["289327"],"_bbp_forum_id":["27791"],"_bbp_engagement":["88","18229","255882"],"_bbp_voice_count":["3"],"_bbp_reply_count":["5"],"_bbp_last_reply_id":["290364"],"_bbp_last_active_id":["290364"],"_bbp_last_active_time":["2023-06-26 13:42:29"]},"test":"tim-hertingersafholland-de"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/289327","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/289327\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=289327"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}