@3rav this is interesting because it means the solver type can go with the model, which is handy for us. I understand how to set the *STATIC command, but how do I have a single ccx.exe and 3 links?
Leave SPOOLES, PARDISO and PASTIX installation in the ccx directory, and it is a trivial matter to switch the link in the Tools-> Options->Calculix menu.
FYI, PASTIX has failed on me for thermals involving pipe flow elements, and has diverg…
Check the Components section in Solution after it fails. Sometimes the solution data file has the projected geometry, and you can see what has happened. Other things to try: change the tolerance to 0.07, also try swapping surfaces.
The *TIE command will project nodes from the slave to master surface. This can mess up the shape of your tet10. Make sure your mating surfaces are well aligned, and use as small a number as possible for the position tolerance.
The "don't generate" commands suppress any material commands you have already defined in mecway, since you will be redefining in custom model input. So define the bare minimum just so that mecway doesn't complain, then suppress whatever you have cre…
use ccx, custom model definition.
define minimum properties in Mecway; there are some missing "don't generate" commands, for example you can't "don't generate" *EXPANSION, so only define these in the CCX custom file
see example
The attached example problem demonstrates this use. I just ran it in the latest Mecway 13.1 using CCX 2.17, and it seems to work, but check it yourself to verify.
Use CCX custom commands below. Ramp may not be necessary, just example.
*AMPLITUDE,NAME=Aramp
0,0
.1,1
1,1
*DLOAD,AMPLITUDE=Aramp
rotor,CENTRIF,10000,0,0,0,1,0,0
The pinned support holds the beam at 1 location in X,Y,Z. The beam can still rotate in RX,RY,RZ direction. The force that is applied puts a moment on the beam, rotating about the pin joint which has no torsional stiffness, and it blows up.
The n…
Went through things a little more methodically, and here is my best current recipe on my 8proc laptop.
PASTIX + all DLLS from PARDISO
OMP_NUM_THREADS=6 (8 ran 20% slower)
OPENBLAS_NUM_THREADS=1 (helped by 8%)
PASTIX_MIXED_PRECISION=1 (helped…
We now run with the PaStiX ccx.exe available on the CalculiX discouse blog.
I typically see about 30-35% speedup over PARDISO.
I tried OPENBLAS_NUM=1 and 4, it seemed to run slower. I deleted and it ran the fastest.
I also noticed that OMP_NU…
Go to the Mecway YouTube channel, there are examples from SimCommons, 4P Ingenieria and VMH that demonstrate capability.
SimCommons examples all run in under 1000 nodes.
I will add to what Sergio said - if you do not have standard samples, use some version of the actual part (if available) and load in a similar fashion to the actual application.
For example, an O-ring would be loaded in compression, but an inflatab…
Gluing some "stiff" elements to the load point allows you to distribute the load from a single point. Think of it as a poor-man's rigid region constraint. Victor has added some features recently that may do this for you, but I still use this trick. …
I second Victor's comment about the displacement versus force. If you model a stiff "holder" where the load is currently applied, you might be able to use a single displacement, which is very handy for post-processing.
Convergence is very sensiti…
Those last few posts were a fun read, and good to see a happy ending. I can relate to the trials that we have to go through to make things work sometimes
For doing coefficient fitting with hyperelastic materials, I have found a neat (and free) piec…
I ran on ccx_pardiso_dynamic.exe using 6CPU
machine is 17-4790CPU @3.6GHz
24GB RAM
ran okay in under 2minutes.
https://drive.google.com/file/d/1TgP01C-XBwdcrqXgCE1HR51hhS4WcePh/view?usp=sharing
I use TetGen, it will create meshes from some brutally ugly surfaces, and also has some useful error checking. I believe this is the mesher used in LS-DYNA.
Victor has put a Lab in place at 13.1 that I am finding very useful.
Often we have a model "all set up" and then we want to look at the effect of tolerance or what impact a small modification might make. In these cases, try this:
1- Read in STL she…
Strange. I loaded the tutorial into Mecway, added the missing commands from the original inp file. If I run, Mecway will crash. But if I save the Mecway model as a CCX run file, and then submit directly to CCX, it works. (LEFT is ORIG CCX, RIGHT is …
I have had to reconcile things like the BOLTCALC program, or weld stress calculators, spring stresses, microprocessor junction temperatures. The issue is the same - there is a "handbook" approach and then the FEA result. In most cases I have found…
Here is a "kit" that should get anyone going with the batch trick.
The attached model is a simple contact test that expects a file called input.inp.
The model outputs a file called out.csv
Create a run.bat using the following inputs:
echo off
ech…