add pathera pardus model

models/Panthera_pardus_RC_30.py

+#! /usr/bin/env python3
+# -*- coding: utf-8 -*-
+# Dolicorhynchops osborni FHSM VP404
+#   10k faces on the mandible surface
+
+
+def parms(d={}):
+    p = {}
+    path = 'Panthera_pardus/150k'
+    p['mandible'] = f'{path}/Panthera_pardus_AMNH_113745_mandible.stl'
+    p['teeth'] = [16.4418, 17.3567, -53.2842]
+    p['LTMJ'] = [-37.433147, -7.689619, 48.622185]
+    p['RTMJ'] = [37.795021,	-7.337307, 49.882557]
+    
+    p['muscles'] = [
+        {
+            'file': f'{path}/m. masseter left.stl', 
+            'force': 150.9284,
+            'focalpt': [-49.8903, 18.6, 30.6564],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/m. masseter right.stl',
+            'force': 230.87,
+            'focalpt': [52.1813, 12.4872, 34.5441],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/m. pterygoideus left.stl',
+            'force': 25.7054,
+            'focalpt': [-10.1567, -2.1217, 34.2502],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/m. pterygoideus right.stl',
+            'force': 35.3892,
+            'focalpt': [11.0602, -3.4042, 33.662],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/m. temporalis left.stl',
+            'force': 221.773,
+            'focalpt': [-20.3207, 11.5864, 82.6092],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}/m. temporalis right.stl',
+            'force': 366.678,
+            'focalpt': [22.9732, 8.3983, 82.0678],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        }
+    ]
+    p['fixations'] = {
+        'teeth': ['x', 'z','y'],
+        'LTMJ': ['z','y'],
+        'RTMJ': ['z','y']
+    }
+
+    # material properties
+    p['density'] = 1.662e-9  # [T/mm³] - bone: 1.850 kg/l
+    p['Young'] = 18000.     # [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))