sorotokicode - Sorotoki - Stable (v2.09.26)

• Improved stability for the Fem.addTendon() for 3D simulations. The FEM solver now ensure rotations are being followed more accurately.

• Fixed issues with LocalsFastYH and LocalsFastMN mex files. the cpp compiller complained about matrix(vector) instead of vector(vector) indexing for the deformationGradient computation. Weirdly the same code does not result in an issue for LocalsFastNH (Neohookean). Issue seems to be unique to Matlab 2022 release, and not Matlab 2021 or earlier...

- MATLAB
Published by BJCaasenbrood over 3 years ago

sorotokicode - Sorotoki - Stable (v2.09.15.2)

Sorotoki - Stable (v2.09.15) Latest

• Minor update to scripts and Readme.md
• mex compile error for LocalsNHFast.m has been fixed (tested for Matlab 2022b).

- MATLAB
Published by BJCaasenbrood over 3 years ago

sorotokicode - Sorotoki - Stable (v2.09.15)

• Minor update to scripts and Readme.md ? Apparently there is a mex compile error for LocalsNHFast.m... Issue remains unclear, perhaps a Matlab versioning issue.

- MATLAB
Published by BJCaasenbrood over 3 years ago

sorotokicode - Sorotoki - Stable (v2.08.17)

DOI release with Zenodo.

- MATLAB
Published by BJCaasenbrood over 3 years ago

sorotokicode - Sorotoki - Stable (v2.07.25)

SOROTOKI - Alpha - 2.07.25 - July 25 - 2022

• Documentation update to the examples Modeling.
• Added minor description in auxiliary functions in /src/base

Fem.m

• Generated new mex files for the Local element computation. Now each matrix computation in during the fem simulation calls a .mex64 file. This significantly improves performance ~40% compare to .m files. During install, Sorotoki will ask to install the new mex files.

• Added (proper) friction effects to the contact simulations. A good example is fembouncingball.m which initially started spinning due to improper implementation of the friction forces during contact. The interaction friction can be modified under Material.Cfr (default = 1e-6 for most sample materials). Added example to show the friction effects is: ./scripts/fem/2D/dynamics/femcrawler.m

? High interaction friction forces lead to instabilities, implement with caution. A fix will be issued later.

• Added initial conditions to the dynamic Fem simulations. Now Fem.Utmp can be overwritten before calling fem.simulate. This can be extremely handy if we like to simulate from a quasi-static equilibrium, e.g., gravity deflection. One can first call fem.solve(), get the displacement through U0 = fem.Log.U(end,:), then a new simulation:

fem.reset(); fem.set('Utmp',U0); fem.simulate();

We can also ensure there is no ramping of the gravity force by setting the following fem.set('GravityRamp',false);

Shapes.m + Major update to the class. Now, Shapes.reconstruct() will produce better POD bases from the finite element simulations. As example:

shp = Shapes(fem,Modal,'NNode',100,'L0',120,'FilterRadius',[15,15]); shp = shp.reference([0,0],[119,0]); shp = shp.reconstruct();

This code will generate a basis from the data in Fem, and project its dynamics onto a curve spanned by the points: (0,0) -> (119,0). During the reconstruction, all system tensors are build accordingly.

• Shapes.show() is implemented and will show the strain basis.

Model.m - Fixed a minus error in the gravitational force computation that lead to incorrect orientation of the gravity vector.

- MATLAB
Published by BJCaasenbrood over 3 years ago