From 00be3897ffd5aad440a0ffcdceda00150393ce50 Mon Sep 17 00:00:00 2001
From: Romain Boman <r.boman@uliege.be>
Date: Tue, 3 Jan 2023 15:30:01 +0100
Subject: [PATCH] add smilodon

---
 ...venile_LACMHC2002LR3_mandible_500k_lfs.ply |  3 +
 .../500k/m. masseter left.stl                 |  3 +
 .../500k/m. masseter right.stl                |  3 +
 .../500k/m. pterygoideus left.stl             |  3 +
 .../500k/m. pterygoideus right.stl            |  3 +
 .../500k/m. temporalis left.stl               |  3 +
 .../500k/m. temporalis right.stl              |  3 +
 .../LACMHC_2002_L_R_3_30.py                   | 80 +++++++++++++++++++
 8 files changed, 101 insertions(+)
 create mode 100644 models/Smilodon_fatalis_juvenile/500k/Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply
 create mode 100644 models/Smilodon_fatalis_juvenile/500k/m. masseter left.stl
 create mode 100644 models/Smilodon_fatalis_juvenile/500k/m. masseter right.stl
 create mode 100644 models/Smilodon_fatalis_juvenile/500k/m. pterygoideus left.stl
 create mode 100644 models/Smilodon_fatalis_juvenile/500k/m. pterygoideus right.stl
 create mode 100644 models/Smilodon_fatalis_juvenile/500k/m. temporalis left.stl
 create mode 100644 models/Smilodon_fatalis_juvenile/500k/m. temporalis right.stl
 create mode 100644 models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py

diff --git a/models/Smilodon_fatalis_juvenile/500k/Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply b/models/Smilodon_fatalis_juvenile/500k/Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply
new file mode 100644
index 0000000..704b629
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:81f9e34fce122fb420ccc282b69be73238a179eb52706bc5e1aa5ea5d849b1f9
+size 9499786
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. masseter left.stl b/models/Smilodon_fatalis_juvenile/500k/m. masseter left.stl
new file mode 100644
index 0000000..3d1b61b
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. masseter left.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:766f2501ed240819add0958920c104c4f7627140348aeed3a0ace4b7bccad8ae
+size 2099797
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. masseter right.stl b/models/Smilodon_fatalis_juvenile/500k/m. masseter right.stl
new file mode 100644
index 0000000..e030bd4
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. masseter right.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:96adcd27e29d766622f54f4a9cda724be8db454120b99461f9c73b4c50d4fe03
+size 2442162
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus left.stl b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus left.stl
new file mode 100644
index 0000000..96341ff
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus left.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0fd6e556989938bc7fb061246119894b5f09451444135b8ebdf472f17618e906
+size 425077
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus right.stl b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus right.stl
new file mode 100644
index 0000000..652b4f6
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus right.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4dab0b4312d219619e69abd75085e7f55ce73d907397fbfb0e0c9760a1ff343d
+size 351754
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. temporalis left.stl b/models/Smilodon_fatalis_juvenile/500k/m. temporalis left.stl
new file mode 100644
index 0000000..617c9e3
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. temporalis left.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:199b75986ecae6a42fe432fd213323ea8b7956913fb20a678ddad1037776d915
+size 2761401
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. temporalis right.stl b/models/Smilodon_fatalis_juvenile/500k/m. temporalis right.stl
new file mode 100644
index 0000000..7e1e68d
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. temporalis right.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5c189fc79791e9bd2cd7aac6e8580e948e587be313c882ddd31a72b27e33f20b
+size 3375261
diff --git a/models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py b/models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py
new file mode 100644
index 0000000..2105643
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py
@@ -0,0 +1,80 @@
+#! \usr\bin\env python3
+# -*- coding: utf-8 -*-
+# Smilodon fatalis juvenile L-R-3
+
+
+def parms(d={}):
+    p = {}
+    import os
+    path = os.path.join(os.path.dirname(__file__),'500k')
+    p['bone'] = f'{path}\Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply'
+    p['contact_pts'] = [ 
+                   [-3.128700, 44.665501, -9.424970] ,
+                   [3.667650, 45.900398, -8.284160],
+                   [9.765930, 44.752102, -6.164700]]
+    p['axis_pt1'] = [-27.858801, -35.784000, -1.227750]
+    p['axis_pt2'] = [22.619600, -43.196201, 4.224330] 
+    p['muscles'] = [
+        {
+            'file': f'{path}\m. masseter left.stl', 
+            'force': 60.7224,
+            'focalpt': [-41.7800, -33.8660, 14.6611],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. masseter right.stl',
+            'force': 75.40299,
+            'focalpt': [41.2037, -39.5222, 14.6911],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. pterygoideus left.stl',
+            'force': 12.51531,
+            'focalpt': [-12.7555, -31.5549, 12.3858],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. pterygoideus right.stl',
+            'force': 9.85383,
+            'focalpt': [10.5836, -32.2269, 12.1060],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. temporalis left.stl',
+            'force': 86.78298,
+            'focalpt': [-22.6034, -79.5285, 22.4788],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. temporalis right.stl',
+            'force': 104.6819,
+            'focalpt': [15.2406, -81.9674, 20.9924],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        }
+    ]
+    p['fixations'] = {
+        'contact_pts': ['z'],
+        'axis_pt1': ['z','y'],
+        'axis_pt2': ['z','y']
+    }
+
+    # material properties
+    p['density'] = 1.662e-9  # [T\mm³]
+    p['Young'] = 18000.      # [MPa]
+    p['Poisson'] = 0.3       # [-]
+
+    p['dt0'] = 0.3       # [s]
+    p['tolNR'] = 1e-3    # [-]
+
+    p.update(d)
+    return p
+
+
+def getMetafor(p={}):
+    import bonemodel as model
+    return model.getMetafor(parms(p))
+
+
+if __name__ == "__main__":
+    import models.bonemodel2 as model
+    model.solve(parms())
-- 
GitLab