add smilodon

models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py

+#! \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())