add 250k model

models/PLTtympFranken_250k.py

+#! /usr/bin/env python3
+# -*- coding: utf-8 -*-
+# Dolicorhynchops osborni FHSM VP404
+#   10k faces on the mandible surface
+
+
+def parms(d={}):
+    p = {}
+    path = 'PLTtymFrankencarpus/250k'
+    p['mandible'] = f'{path}/PLTtym_250k_symph_lfs.stl'
+    p['teeth'] = [[-140.616379, 34.634388, 178.193710],
+                [-141.087341, -31.283962, 177.457901]]
+    p['LTMJ'] = [206.616394, -80.626038, 187.317291]
+    p['RTMJ'] = [202.435074, 114.139923, 188.697571]
+    p['muscles'] = [
+        {
+            'file': f'{path}/PLTtym_M_Adductor_externus_L1.stl',
+            'force': 487.38,
+            'focalpt': [146.982380, -77.664399, 278.479435],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Adductor_externus_L2.stl',
+            'force': 487.38,
+            'focalpt': [146.982380, -77.664399, 278.479435],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Adductor_externus_R1.stl',
+            'force': 452.79,
+            'focalpt': [138.878155, 84.493762, 284.866894],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Adductor_externus_R2.stl',
+            'force': 452.79,
+            'focalpt': [138.878155, 84.493762, 284.866894],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Adductor_externus_medialis_L.stl',
+            'force': 362.86,
+            'focalpt': [147.360255, -67.584745, 281.176887],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Adductor_externus_medialis_R.stl',
+            'force': 354.37,
+            'focalpt': [136.040621, 69.415241, 285.431467],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Adductor_mandibulae_profundus_L.stl',
+            'force': 233.63,
+            'focalpt': [190.286665, -73.222918, 239.986014],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Adductor_mandibulae_profundus_R.stl',
+            'force': 225.83,
+            'focalpt': [187.649314, 87.414657, 236.135409],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Pseudotemporalis_profundus_L_lfs.stl',
+            'force': 1537.52,
+            'focalpt': [153.241871, -37.103172, 275.752858],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Pseudotemporalis_profundus_R_lfs.stl',
+            'force': 1458.46,
+            'focalpt': [153.541222, 41.630166, 284.082309],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Pseudotemporalis_superficialis_L.stl',
+            'force': 272.44,
+            'focalpt': [116.915364, -29.987994, 276.369899],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/PLTtym_M_Pseudotemporalis_superficialis_R.stl',
+            'force': 272.42,
+            'focalpt': [116.915364, 29.987994, 276.369899],
+            'method': 'T'                              # 'U', 'T', 'T+N'
+        }
+    ]
+    p['fixations'] = {
+        'teeth': ['z'],
+        'LTMJ': ['x', 'z','y'],
+        'RTMJ': ['x', 'z','y']
+    }
+
+    # material properties
+    p['density'] = 1.850e-9  # [T/mm³] - bone: 1.850 kg/l
+    p['Young'] = 20000.      # [MPa] elastic modulus - bone: 17-20 GPa
+    p['Poisson'] = 0.3       # [-] Poisson's ratio
+
+    p.update(d)
+    return p
+
+
+def getMetafor(p={}):
+    import mandiblemodel as model
+    return model.getMetafor(parms(p))