Uniform format for float in ccx input file

Author: FrozenWinterheart

examples/crankshaft/crankshaft.py

@@ -101,9 +101,13 @@
     )
 
     f_res = 500.
-    for t in np.linspace(0,1,12, endpoint=False):
-        fx, fy = f_res * np.cos(np.pi*t), f_res * np.sin(np.pi*t)
-        # Add step
+    # Apply the force f_res in 12 steps from 0° to 330°
+    n_steps = 12
+    for i in range(n_steps):
+        # calc force components
+        angle = i / n_steps * 2 * np.pi 
+        fx, fy = f_res * np.cos(angle), f_res * np.sin(angle)
+        # make new step
         step = sk.Step(nlgeom=True)  
         load = sk.Cload(ref_load_1 , 1, fx)
         load.add_load(ref_load_1, 2, fy)

pygccx/auxiliary.py

@@ -0,0 +1,5 @@
+from protocols import number
+
+def f2s(x:number) -> str:
+    """Returns a string of the given number in the form 15.7e"""
+    return f'{x:.7e}'

pygccx/mesh/mesh.py

@@ -24,6 +24,7 @@
 from .element import Element
 from .set import Set
 import protocols
+from auxiliary import f2s
 
 @dataclass(repr=False)
 class Mesh:
@@ -322,7 +323,8 @@ def _write_nodes_ccx(self, buffer:list[str]):
         if not self.nodes: return
         buffer += ['*NODE']
         for nid, coords in self.nodes.items():
-            buffer += ['{},{:15.7e},{:15.7e},{:15.7e}'.format(nid, *coords)]
+            buffer += [f'{nid},' + ','.join(map(f2s, coords))]
+            # buffer += ['{},{:15.7e},{:15.7e},{:15.7e}'.format(nid, *coords)]
 
     def _write_elements_ccx(self, buffer:list[str]):
         if not self.elements:return

pygccx/model_keywords/amplitude.py

@@ -23,6 +23,7 @@
 import numpy.typing as npt
 
 from protocols import number
+from auxiliary import f2s
 
 @dataclass
 class Amplitude:
@@ -56,12 +57,12 @@ class Amplitude:
     def __str__(self):
         s = f'*AMPLITUDE,NAME={self.name},'
         if self.use_total_time: s += 'TIME=TOTAL TIME,'
-        if self.shift_x is not None: s += f'SHIFTX={self.shift_x},'
-        if self.shift_y is not None: s += f'SHIFTY={self.shift_y},'
-        s = f'{s[:-1]}\n' # delete last ','
+        if self.shift_x is not None: s += f'SHIFTX={f2s(self.shift_x)},'
+        if self.shift_y is not None: s += f'SHIFTY={f2s(self.shift_y)},'
+        s = s.rstrip(',') + '\n'
 
 
         for t, a in zip(self.times, self.amps):
-            s += f'{t:15.7e},{a:15.7e}\n'
+            s += f'{f2s(t)},{f2s(a)}\n'
 
         return s

pygccx/model_keywords/clearance.py

@@ -19,6 +19,7 @@
 
 from dataclasses import dataclass
 from protocols import ISurface, number
+from auxiliary import f2s
 
 @dataclass
 class Clearance:
@@ -45,4 +46,4 @@ class Clearance:
     """A short description of this Instance. This is written to the ccx input file."""
 
     def __str__(self):
-        return f'*CLEARANCE,MASTER={self.master.name},SLAVE={self.slave.name},VALUE={self.value}\n' 
+        return f'*CLEARANCE,MASTER={self.master.name},SLAVE={self.slave.name},VALUE={f2s(self.value)}\n' 

pygccx/model_keywords/contact_pair.py

@@ -21,6 +21,7 @@
 from typing import Optional, Any
 from enums import EContactTypes, ESetTypes, ESurfTypes
 from protocols import IKeyword, ISet, ISurface, number
+from auxiliary import f2s
 
 @dataclass
 class ContactPair:
@@ -95,7 +96,7 @@ def __str__(self):
         if self.small_sliding: s += f',SMALL SLIDING'
         if self.adjust is not None:
             if isinstance(self.adjust, ISet): s += f',ADJUST={self.adjust.name}'
-            else: s += f',ADJUST={self.adjust}'
+            else: s += f',ADJUST={f2s(self.adjust)}'
         s += '\n'
 
         s += f'{self.dep_surf.name},{self.ind_surf.name}\n'

pygccx/model_keywords/distributing_coupling.py

@@ -21,6 +21,7 @@
 from typing import Any
 from protocols import ISet, number
 from enums import ESetTypes
+from auxiliary import f2s
 
 @dataclass
 class DistribuitingCoupling:
@@ -90,8 +91,8 @@ def __str__(self):
 
         for nid_or_set, weight in self.conditions:
             if isinstance(nid_or_set, ISet):
-                s += f'{nid_or_set.name},{weight:.7e}\n'
+                s += f'{nid_or_set.name},{f2s(weight)}\n'
             else:
-                s += f'{nid_or_set},{weight:.7e}\n'
+                s += f'{nid_or_set},{f2s(weight)}\n'
 
         return s

pygccx/model_keywords/elastic.py

@@ -20,6 +20,7 @@
 from dataclasses import dataclass, field, InitVar
 from enums import EELasticTypes
 from protocols import number
+from auxiliary import f2s
 
 @dataclass
 class Elastic:
@@ -88,7 +89,7 @@ def __str__(self):
         for p in self.elastic_params_for_temps:
             lines = [p[i:i+n] for i in range(0, len(p) or 1, n)]
             for i, line in enumerate(lines):
-                s += ','.join(map(str, line)) 
+                s += ','.join(map(f2s, line)) 
                 s += '\n' if i == len(lines) -1 else ',\n'
 
         return s

pygccx/model_keywords/friction.py

@@ -21,6 +21,7 @@
 from typing import Any
 
 from protocols import number
+from auxiliary import f2s
 @dataclass
 class Friction:
 
@@ -54,5 +55,5 @@ def __setattr__(self, name: str, value: Any) -> None:
 
     def __str__(self):
         s = '*FRICTION\n'
-        s += f'{self.mue},{self.lam}\n'
+        s += f'{f2s(self.mue)},{f2s(self.lam)}\n'
         return s

pygccx/model_keywords/gap.py

@@ -23,6 +23,7 @@
 import numpy.typing as npt
 
 from protocols import ISet, number
+from auxiliary import f2s
 from enums import ESetTypes
 
 @dataclass
@@ -72,9 +73,10 @@ def __setattr__(self, name: str, value: Any) -> None:
     def __str__(self):
         s = f'*GAP,ELSET={self.elset.name}\n'
         nx, ny, nz = self.direction
-        s += f'{self.clearance:.7e},{nx:.7e},{ny:.7e},{nz:.7e},'
+        s += f'{f2s(self.clearance)},{f2s(nx)},{f2s(ny)},{f2s(nz)},'
+        
 
-        s += f',{self.k:.7e}' if self.k is not None else ','
-        s += f',{self.f_inf:.7e}' if self.f_inf is not None else ','
+        s += f',{f2s(self.k)}' if self.k is not None else ','
+        s += f',{f2s(self.f_inf)}' if self.f_inf is not None else ','
 
         return s.rstrip(',') + '\n'