Home

CalculiX USER`S MANUAL - CalculiX GraphiX, Version 2.9 -

Contents

1. off bg Pw k 95 Commands to control the graphic output This command is intended for batch mode See also Viewing for the menu controlled functions e fill Element faces are filled e line Elements are displayed as fireframes e point Element edges are displayed as points e edge triggers the display of the model edges were e elem triggers the display of the element edges They are deselected with the additional parameter off e surf and volu are used to display the structure either only by it s outer skin surf or by drawing all elements volu e front and back define which side of the structure should be drawn Either the side which faces the user or the back side If the back side is displayed then internal structures are visible e vec triggers the vector mode All vector entities like displacements are displayed with arrows pointing in the direction of the vector and with a length proportional to the value of the vector See Toggle Vector Plot for a detailed description of the equivalent menu function e disp will show the deformed structure based on a formerly selected dis placement dataset no entity must be selected See Toggle Add Displacement for a detailed description of the equivalent menu function e bg Without second parameter toggles the background colour The sec ond parameter w forces white while the
2. bia lt line gt lt bias gt bia lt set gt lt bias gt l L mult div lt factor gt This keyword is used to define the bias of a single line or of a set of lines see qadd The bias defines the ratio of the length of the first element to the length of the last element For example bia all 4 5 will force a ratio in which the last element is 4 5 times bigger than the first one Real numbers are permitted since version 1 5 see also qbia To convert from pre 1 5 versions start the program with the oldbias option 8 5 body body lt name char lt 9 gt lt surf1 gt lt surf2 gt lt surf1 gt lt surf2 gt lt surf3 gt lt surf4 gt gt lt surf5 gt lt surf6 gt lt surf6 gt lt set gt This keyword is used to define or redefine a volume body Each body must have five six or seven surfaces to be mesh able with hexaeder elements oth erwise it can only meshed with tets if NETGEN 4 is installed However it is sufficient to specify just the top and the bottom surfaces But if sur faces with 3 or 5 edges are involved then this surfaces have to be the top and bottom surfaces This is also true if surfaces have different line divisions at opposite edges The missing surfaces between the top and bottom surfaces will be created automatically if they do not already exist they will always have 4 edges with the same division on opposide edges But all need
3. dsl provides the same functionality 7 4 Frame Adjusts the drawing box 7 5 Zoom Use this command to zoom into a rectangular section of the window After this option is chosen use the mouse to select the opposite corners of a rectangle The display will zoom in on the rectangular area Note the rectangle is never shown on the screen see also Zoom 7 6 Center Used to choose a new center point for the structure After this option is chosen pick either a node a point or the corner of an entity To easily find the element corners the function Toggle Element Edges is triggered automatically see also aent 7 7 Enquire Used to investigate parameters like the value and the position of a certain node of the model Pick a node after this option is chosen To easily find the element corners the function Toggle Element Edges is triggered automatically see also RENj 21 7 8 Cut Used to cut elements and to create a section of new elements and nodes Either pick three nodes or in case a dataset entity of a vector was already selected use the menu entry vector and select just one node The cutting plane is then determined by the direction of the vector displacements worstPS To eas ily find the element corners the function Toggle Element Edgesl is triggered automatically see also gcut and ut 7 9 Graph Used to generate a 2D plot The option Length will provide a plot value ove
4. LINE L009 P006 P009 pO 4 LINE LOOA P009 POOA pO 4 LINE LOOC POOA pO 8 LINE LOOG POOA POOG pO 4 LINE LOOI POOG POOI pO 4 GSUR A001 BLEND L003 L002 L001 L004 GSUR A002 BLEND L007 L006 L005 L003 AS 128 Py Figure 13 disc made of four 90 degree segments GSUR A003 BLEND LOOC LOOA LOO9 L007 GSUR A004 BLEND L004 LOOI LOOG LOOC ELTY all QU4 C 2 Cylinder PNT po 0 00000 0 00000 0 00000 PNT py 0 00000 1 00000 0 00000 PNT p1 1 00000 0 00000 0 00000 PNT P001 1 00000 1 00000 0 00000 PNT P002 0 00000 0 00000 1 00000 PNT P003 0 00000 1 00000 1 00000 PNT P006 1 00000 0 00000 0 00000 PNT P007 1 00000 1 00000 0 00000 PNT POOA 0 00000 0 00000 1 00000 PNT POOC 0 00000 1 00000 1 00000 LINE L001 pi P001 2 129 LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE GSUR GSUR GSUR GSUR ELTY L002 L003 L004 L005 L006 L007 L008 L009 LOOA LOOI LOOJ A001 A002 A003 A005 Figure 14 cylinder made of four 90 degree segments P002 P003 2 pi P002 po 8 P001 P006 P002 P003 POOA P006 P007 POOA POOC BLEND BLEND BLEND BLEND all QU4 P003 py 8 P007 2 P006 po 8 P007 py 8 POOC 2 POOA pO 8 POOC py 8 pi po 8 P001 py 8 L001 L002 L005 L008 L003 L006 L009 LOOI 130 L002 L005 L008 L001 L004 L007 LOOA LOOJ C 3 Sphere Figure 15 Segment of a Sphere PNT py
5. Shaded surfaces plot si all can only be displayed if the interiour was previously calculated which is done with the command rep The attribute t applies only to nodes and will display only the ones which have attached texts plot nt all will display only the nodes which have attached texts out of the set all They are created with qadd aeng or Tatxt The attribute q works only for elements With this attribute only elements which do not pass the element quality check are displayed plot eq all The threshold for the element quality is defined with Tegall The attribute v works only for nodes plot nv setl will display the nodes in set setl with their values This is only meaningful in post processing mode see how to assign values to nodes To plot additional entities see plus 8 61 plus plus Dn e f p 2 Ig pb S L sh si 1 gt 52 Pa rg p gq 4 ry gt lt set gt w kK Yr g pb y m n 4 This keyword is used to display the entities of an additional set after a command was used see also minus Further details are explained in section 8 62 pnt pnt lt name char lt 9 gt lt x gt lt y gt lt z gt lt line gt lt ratio gt lt times gt lt P1 gt lt P2 gt lt ratio gt lt times gt lt setname containing nodes gt This keyword is used to define or redefine a point There are four possibil
6. The data structures for the mesh and datasets are available The default function calculates the hydrodynamic stresses with the command call hydro See JUser Functions for details 8 8 cntr entr lt pnt nod gt Defines a new center point or node This command is intended for batch mode See also qent for the cursor controlled command 8 9 comp comp lt set chars gt u d e 27 This keyword is used to add all entities to the specified set see seta which de pend on the already included entities u up or to include all entities necessary to describe the already included entities d down For example the set lines stores lines and should also include all dependent points comp lines do Or the set lines should also include all surfs and bodies which depend on the lines comp lines up In some cases you will need only the end points of lines With the option e edges comp lines e only end points are included in the set One exception to this logic was in troduced for convenience comp nodes do will add all faces described by the nodes in set nodes despite the fact that faces are made from nodes Wildcards can be used to search for setnames of a certain expression comp E do will complete all sets starting with E 8 10 copy copy lt set gt lt new_set gt scal lt fx gt lt fy gt lt fz gt lt pnt gt a tra lt dx gt lt d
7. WN S NS NS TA Figure 17 Airfoil for cfd codes 24375 0 00427 0 23750 0 00565 0 23125 0 00653 0 22500 0 00730 0 21250 0 00875 0 20000 0 00993 0 18750 0 01070 0 17500 0 01115 0 15000 0 01120 0 12500 0 01043 0 10000 0 00917 0 07500 0 00750 0 05000 0 00540 0 02500 0 00308 0 01250 0 00175 0 00000 0 00000 0 50000 0 50000 0 134 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 W g or UV UV UT FO UV uo UT TU UV TV UT UV UT UT TU UT TU UV UT TU UV UT UT TU UT UT TU UT UT UT UT UT UT uv H H 2422224222422 2422 242 242 2 242 242 242 2 2424222422242 2422 2422 242224232 2423422222422 2422422424244 Cir DA DA Da DA a DD DA a Dr DN Dr Dr Dr Dr Dr Dr Dr Dr Dr Dr Dr DD qr Dr Dr Dr Dr Dr Dr Dr Dr Dr Dr DD Dr Dr Dr Dr Dr Dr Dr Dt H P2 p3 p4 PO1A PO1B P046 P059 P049 POSA P04C PO4D PO5C P04G PO4H PO4I POSD PO4L PO4M PO5F PO4P PO4R POAT PO5G PO4W PO4X P04Z PO5H P052 P054 P056 P058 POSI P05J PO5L PO5M POSN POSS POSV POON PO2M PO3B P03C PO3E PO3P pli pl2 50000 50000 50000 18162 18180 27025 26599 25144 35589 22636 20128 38027 12604 07588 02572 31932 27652 27025 20962 25373 22676 20124 01132 12604 07588
8. 02572 00065 18074 18133 00118 00010 01244 00610 50000 50000 50000 50000 50000 23448 23471 24164 23405 24536 24464 00000 00000 I I o OGOOGO oo0DO0OO0O0O0Oo0oOooOoOo 50000 50000 50000 01898 01094 01256 28688 02439 17566 03241 03643 00958 03833 03616 03231 21136 00154 00803 27840 01567 02143 02394 29145 02508 02519 02355 30887 03754 02465 02891 02250 50000 50000 34112 29206 05780 05314 00217 01345 00608 02804 02029 00794 00400 00000 00000 135 oo0DOOOO0OO0O0O0O0O0O0O0O0O0O0O0O0O0O0O0O000000000000000000000000o0o00o 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 SEQA SEQA SEQA SEQA SEQA SEQA SEQA SEQA SEQA SEQA SEQA SEQA SEQA LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE S006 S007 SOOW SOOR SooS SOOT S001 SOOE SOOL SOOA SOOX S002 SOOP L003 LOOC LOOE L004 LO5F LOSS
9. LO5C L006 L007 L008 LOOG LOON LO3R LOOI LO3S LO4V LO4W LO4X L04Z L050 L052 L054 L056 L058 L059 LOSA LO5D LOSV LOSI LO5L LO5Z LO6C LO6F pnt pnt pnt pnt pnt pnt pnt pnt pnt pnt pnt pnt pnt PO1A PO1A P019 PO1B P019 P058 PO5S PO3B P052 PO5V P058 PO1A P019 PO3B P052 P058 PO5F PO1A PO3E P052 PO2M P03C POON PO3P P052 POOA P017 P049 P04G P04Z PO5D P009 P002 PO4D POOT PO4R P004 POOR 910 POOC P016 P046 PO4H PO4X PO5C P008 P000 PO4C POOV PO4T P003 POOS P019 S006 120 PO1B S007 120 p4 P1 150 P2 PO5L 204 PO5L PO5S 912 P1 PO5J 120 P054 P056 P019 PO5V PO3C P056 PO5N P054 910 910 910 210 SOOW SOOR 210 5005 p4 PO59 204 P059 P054 PO5F PO5F P059 PO5F PO5H P056 P052 PO5F P1 20 P059 PO5G PO5H P058 POSG 912 912 4 SOOT 120 120 912 912 p3 POSI 130 POSI P05J P056 PO5M PO5N P058 P019 PO5M p4 12 P2 13 PO5N p3 20 PO5M PO5S PO5V PO5G 0 0 130 4 912 130 130 130 130 120 120 150 POOE P015 PO4L PO4I PO4W PO5A P007 POOP PO3B POOW P054 PO3E PO2M 136 POOG P014 PO4M P056 P054 P059 P006 PO3P POOX POOI POOJ POOK POOL P019 P013 P012 P011 P010 PO1B PO4P PO3C P005 POON P00Z PO1B LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LINE LCMB LCMB GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR GSUR SETA SETA SET
10. Or the point can be moved to the position of a second point for this go to the second point which coordinates should be used and press p again If the coordinates of a node should be used press n instead Also a normal projection to a nurbs related surface from cad systems is feasible by choosing the target surface with the s key or a NURBS surface with the S key either before or after the point was marked Press q to quit the command or u to undo the last action If you picked the wrong point the one which should be moved just pick the same again and pick then the correct one It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 65 8 82 qnod gt gnod RETURN gt w Pp PR Pu This keyword is used to move nodes see also node To move a node move the mouse pointer to the desired node and press p pick then go to the desired position and press m move See galiJor the sectionlOrientation how to rotate the model into a certain position Press q to quit the command or u to undo the last action It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 83 qrem gt grem l
11. The eq and ne compare strings and should not be used for nummerical values since no conver sion to a common format is done Two strings are equal if they have the same length and all characters are equal See also while valul Btacil and How to run cgx in batch model 8 39 Icmb 1cmb lt name char lt 9 gt lt line gt lt line gt gt lt line gt up to 14 lines ADD lt line gt lt line gt gt 42 lt line gt up to 14 lines This keyword is used to define extend or redefine a combined line Iemb Com bined lines are necessary if the edge of a surface should be made of more than one line Usually the user does not create lemb s directly They are created automatically during the process of defining a surface with the command asur There is no limitation to the number of lines in a combined line However with one command not more than 14 lines can be specified at a time To specify more than that or to extend an existing lemb a modify command has to follow For example lemb U260 U249 U248 U247 U243 U237 U236 U231 U219 defines the lemb U260 with 8 lines and their orientation in the lemb The following command lemb U260 ADD U218 U217 extends the lemb U260 by two additional lines 8 40 length length lt set gt This keyword is used to calculate the length of all lines stored in a set The command writes to th
12. This is just for visualization Meshing takes place in NETGEN 4 see the read command to handle edges generated by NETGEN It is hoped by the author that common CAD formats will be supported by stand alone interfaces which translate into fbd commands So far a vda to fbd interface is available Tet meshes can be generated based on the resulting fbd files The following file formats are available to write a mesh and certain boundary conditions Abaqus which is also used by the CalculiX solver cecx Ansys most boundary conditions available Code_Aster mesh and sets of nodes and elements are available Samcef mesh and sets of nodes and elements are available dolfyn a free cfd code 5 duns a free cfd code 6 isaac a free cfd code 7 OpenFOAM a free cfd code 8 only 8 noded brick elements are sup ported Nastran most boundary conditions available tochnog a free fem code 9 only 8 noded brick elements are supported The following solver input file formats can be read to check the mesh sets and certain boundary conditions Abaqus this is also used by the CalculiX solver ccx Netgen read Netgen native format vol The following file formats are available to read solver results frd format files of this format are used to read results of previous calcula tions like displacements and stresses This format is described in section Result Format It is also used by cecx duns a free
13. ine or in combination with the graph command if values along a path should be displayed Use prnt sq to list all existing sequences and use qseq or seqaj to create them Up to now they are only used to store nodes and points As sets they will be written to a file if fbd format is specified In this case also node sequences can be stored B 3 How to enquire node numbers and values at certain locations A very common problem is how to get the exact value on a node position during post processing To actually get the value add the element edges to the view see Toggle Element Edges and type geng and press the RETURN key during typing the mouse pointer MUST stay in the main window do NOT move the pointer into the xterm Then move the mouse pointer over the location of a node and press the n key The node number the value at that node and the location will be written in the xterm from which the cgx was started See also the feng command for batch controlled value extraction B 4 How to write values to a file When you need to write certain values see Tvalul to a file you may use the echo system command e Choose the relevant dataset menu or ds command e Add the relevant nodes to a set qadd setname or seta setname n Open the stack stack on e Place node values on the stack prnt se setname Store the node valu nod1 pop e Store the value valu vall pop e Write to file sys echo Node nod1 Value v
14. inp the parameter inp might be omitted An openFoam file 8 can be read in a similar way read CaseDirName foam as an ccx input file Sets will be created if defined in the openFoam file Results can be used for mapping purposes For further details see How to map loads in the appendix The user might read a result file in frd format during runtime If a mesh exists and should not be overwritten just add the parameter add read result file add this forces the program to change the numbers which are already in use to unused ones Existing datasets will be extended by the new node numbers and their data The option ext read result file ext will also extend the existing datasets but in this case the nodes and elements are updated modified If no parameter follows then existing nodes and elements are updated and the new datasets will be appended to existing ones It is also possible to read files written in the result format frd during run time only to define sets of nodes or elements without changing the definitions of them The command read result file setname will create a set of the name setname and all nodes and elements listed in the file mesh frd will be added to that set But this will NOT create or modify the nodes and elements All nodes and elements must exist and will not be changed An stl file can be read with read file stl If a file with the numbers of nodes or elements is not available
15. kz amp lt o o F Ju E Lo ER O Zz E 2 5 L JO 2 E gt L JO o L JO i i i fi i i i o o o o o o o o o o o o o o O o o o o o O 00 N co LO Y mM N S l Figure 10 example of a time history plot B 7 How the mesh is related to the geometry Very often you need the embedded nodes element faces or elements of geometric entities to apply boundary conditions If you understand the underlying concept you will be able to get them more easily Geometric entities are the mothers of nodes faces and elements so to say and will remember them In turn if an entity is not the mother of a given mesh entity it will not remember it Therefore a body will only know about nodes which are not placed on surfaces lines or points A surface will only know about nodes which are not placed on lines or points And so on Therefore if you need the nodes on a surface and not only the ones just inside 115 the surface then create a set with this surface and do a comp set do This will add the lines and points together with their nodes to the set On the other hand if you have only a node and you need the geometric entity in which the node is embedded you might also type comp set do B 8 How to change the order of elements Use the command mids to change from second order to first order or vice versa In this case the amount of elements will not change Or use the command end with the parameter quadlin to change from se
16. merg s all Creation of the LOAD set and ANCHORAG one gadd LOAD use rr keys to select the surface R qadd ANCHORAG use rr keys to select the surface 126 NOTA It s easy to verify the different sets for example K plot b all you can see all the volumes K plus s LOAD you can see the set LOAD K plus s ANCHORAG K plus s SYMETRY 9 Mesh K plot ld all K div all mult 2 K elty all HE20 to specify HEXAHEDRA with 20 nodes K elty LOAD qu8 to mesh the set LOAD otherwise no quads will be created K elty ANCHORAG qu8 K elty SYMETRY K mesh all to mesh the part with all dat name K send all sam to export the mesh into Samcef format K send LOAD sam nam to export groups into Samcef format K send ANCHORAG sam nam see previous remark K send SYMETRY sam nam 10 Modifications It s possible now to make some modifications a open all dat file with your favorite text editor Vi for me b open ANCHORAG nam amp the SYMETRIC nam files and do the same as previously c concatenate under Linux the files using the following schema cat all dat LOAD nam gt s1 m cat s1 m ANCHORAG nam gt s2 m cat s2 m SYMETRY nam gt part dat all the sx m files will be erased afterward d open PART dat file and go to the end gt then add RETURN e the mesh file now works with Samcef Another interesting way add for each nam file an input in your bank file input part dat inpu
17. ry send sur abq slide ry or use ty if the nodes should move only tangential to the y axis In both cases the axial movement is still permitted The parameter sur will write the surfaces of the mesh either in abaqus format for the SURFACE command or in frd format as shell elements For example send top abq sur will write the elements in set top together with the face nr Important Only 79 free surfaces of the mesh are regarded internal faces are unknown and can not be identified The front or rear side of the face pos or neg is selected with the or parameter send surset abq sur will write all surfaces in set surset in abaqus format but all surfaces which have a potentially negative side will be written in that manner The parameter pres is used to assign pressure values to element faces For example send surf abq pres 0 05 will assign the pressure value 0 05 to all element faces in set surf If a dataset with nodal values is available it is also possible to use this values instead of using a uniform value for all faces For example send surf abq pres lc1 el will specify the dataset nr 1 and the entity nr 1 to be used instead of a uni form value The front or rear side of the face pos or neg is selected with the or parameter it has to be the last parameter The parameter film is used to assign free stream temperatures and ther mal heat coefficients to element faces I
18. will not be deleted It is just not longer a member of that set To remove entities from the set dummy type setr dummy p pl p2 This will remove the points pl and p2 from the set The following entities are known Nodes n Elements e Points p Lines 1 Surfaces s Bodies b Nurb Lines L Nurb Surfaces S Values v other sets se and shapes sh The program will automatically determine the type of the entities if not specified but then the names must be unique Wildcards can be used to search for setnames of a certain expression 8 105 shpe gt shpe lt name gt pln lt P1 gt lt P2 gt lt P3 gt cy1 lt P1 gt lt P2 gt lt R1 gt con lt P1 gt lt P2 gt lt R1 gt lt R2 gt sph lt P1 gt lt R1 gt This keyword is used to create a shape which can be used to define the interiour of surfaces or to be used as a target for projections see proj or to split entities see split A plane shape is defined with the parameter pln followed by the names of three points shpe H001 pln P1 P2 P3 A cylinder can be defined by shpe H001 cyl P1 P2 10 5 were 10 5 is the radius A cone can be defined by shpe H001 con P1 P2 10 5 2 were 10 5 is the radius at point P1 and 2 at point P2 A sphere can be defined by shpe H001 sph P1 10 5 were 10 5 is the radius at point P1 If automatic name generation is desired then 89 use instead of a name See also qshp for the mouse controlled defini
19. worstPS Be aware of the key u undo to return to the un cutted structure See also cut for the command line function It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 71 qdel qdel RETURN wW Pe PP PEP PPPPPELPPPE PEL PHR g 58 LC39 STRESS Tim 13 000000 Entity Mises 33e 03 54e 01 e max min 21e 02 98e 02 74e 02 51e 02 27e 02 04e 02 30e 02 57e 02 33e 02 10e 02 856e 02 63e 02 39e 02 16e 02 92e 02 69e 02 46e 02 22e 02 86e 01 51e 01 17e 01 82e 01 NON ORPRRPPRPRNNNNW WDD na A DR KR W latfmtmf2 frd Figure 5 qcut A section through a model defined by three nodes This keyword is used to delete entities see also del Higher entities depending ones will be deleted to To delete more than one entity with one stroke type a before and create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the rectangle To delete only one entity type i before Press q to quit the command It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 72 qdis qdis RETURN gt w PE Pp lg This k
20. 0 00000 1 00000 0 00000 PNT pi 1 00000 0 00000 0 00000 PNT P001 0 70711 0 00000 0 70711 PNT P003 0 00000 0 00000 1 00000 PNT P006 0 70711 0 50000 0 50000 PNT P008 0 00000 0 70711 0 70711 PNT POOC 0 70711 0 00000 0 00000 PNT POOK 0 70711 0 70711 0 00000 PNT POOL 0 00000 0 00000 0 00000 PNT POON 0 00000 1 00000 0 00000 LINE L001 pi P001 POOL 8 LINE L002 P001 P003 POOL 8 LINE L003 pi P006 POOL 8 LINE L004 P006 P008 POOL 8 LINE L006 P001 P006 POOC 8 131 LINE L008 P003 P008 POOL 8 LINE LOOA pi POOK POOL 8 LINE LOOC POOK POON POOL 8 LINE LOOG P006 POOK POOC 8 LINE LOOJ P008 POON POOL 8 GSUR A005 BLEND L003 LOO1 L006 GSUR A002 B beginfappendix LEND L002 L006 L004 L008 GSUR A006 BLEND L003 LOOG LOOA GSUR A004 BLEND L004 LOOG LOOC LOOJ ELTY all QU4 C 4 Sphere Volume Figure 16 Segment of a Sphere Volume PNT py 0 00000 1 00000 0 00000 PNT p1 1 00000 0 00000 0 00000 132 PNT P006 0 70711 0 50000 0 50000 PNT P008 0 00000 0 70711 0 70711 PNT POOC 0 70711 0 00000 0 00000 PNT POOK 0 70711 0 70711 0 00000 PNT POOL 0 00000 0 00000 0 00000 PNT POON 0 00000 1 00000 0 00000 LINE L001 pi POOL 8 LINE L002 POOL P008 8 LINE L003 pi P006 POOL 8 LINE L004 P006 P008 POOL 8 LINE L005 POOL POON 8 LINE LOOA pi POOK POOL 8 LINE LOOC POOK POON POOL 8 LINE LOOG P006 POOK POOC 8 8 LINE LOOJ P008 POON POOL GSUR A001 BLEND L003 LOO1 L002 L004 GSUR
21. 1 gt 2 gt 1 gt 1 abs max min pow gt gt sqr sin cos tan asin acos atan gt int float l exp gt name name An entity which stores a string of characters Most characters are valid but no white spaces are accepted from the command line The command parser will scan each parameter of each command for expressions which match the name of a value and will replace the parameter by the content of the value For example pnt PO xvalue 0 0 uses the value xvalue If the user has previously defined the value with valu xvalue 1 24 then the command parser will replace xvalue by 1 24 in the pnt command ATTENTION The commands prnt all commands starting with set except seta and setr and valu itself are not parsed for substitution To suppress the substitution the value has to be masked by a leading For example the command del v value will delete the value itself The content of a value can be defined by the user valu argl 1 24 or derived from the stack previously filled by a command valu argl pop if the stack was activated A certain number on the stack can be addressed with valu argl pop 2 were the 2 means that 2 successive pops are executed In this way the second value on the stack can be loaded at once Values can be written to the stack when the push parameter is used valu argl pu
22. 2 3 2 2 o Ogg Oyy 0220 0550 yy OyyOzz O22T gg Ory Oya 2 2 2 2 pel 0730 yyO 22 20 74042020 un Trax O yz OyyO zz T O 220 zy 0 11 3 maxShear Stresses Entity name maxShear This is the maximum shear stress maxShear 0 5 01 03 11 4 Cylindrical Stresses Entity names SXX STT SRR SXT STR SRX The Cylindrical Stresses are truly cartesin stresses in a cylindrical system re garding the node position relative to the axis of the cylindrical system The stress tensor is rotated individually for each node They are calculated on de mand see trfm 12 User Functions The user might define his own functions to manipulate the mesh or the re sults with the user function stored in file userFunction c See the command Icalll how to call a user function The user can generate new nodes elements or datasets or extend or manipulate existing datasets or interfaces to other software The file userFunction c includes an example which calculates the hydrostatic stress and stores the result in a new dataset There the user can see how to deal with datasets A Known Problems A 1 Program is not responding If the program seems to hang then leave the window with the mouse pointer and go in again If that does not help then probably a command is waiting for input Stay with the mouse pointer inside the window and press the q key several times Another very popular error is to move the
23. A002 BLEND L005 LOO1 LOOA LOOC GSUR A006 BLEND L003 LOOG LOOA GSUR A004 BLEND L004 LOOG LOOC LOOJ GSUR A003 BLEND L002 LOOJ L005 GBOD BOO1 NORM A006 A003 A004 A002 A001 ELTY all HE20 C 5 Airfoil for cfd codes All surfaces must be oriented in the same way The sets are used to define areas for the boundary conditions PNT P002 0 24688 0 00667 0 00000 PNT P003 0 24375 0 00903 0 00000 PNT P004 0 23750 0 01228 0 00000 PNT P005 0 23125 0 01450 0 00000 PNT P006 0 22500 0 01608 0 00000 PNT P007 0 21250 0 01798 0 00000 PNT P008 0 20000 0 01875 0 00000 PNT P009 0 18750 0 01900 0 00000 PNT POOA 0 17500 0 01888 0 00000 PNT POOC 0 15000 0 01785 0 00000 PNT POOE 0 12500 0 01602 0 00000 PNT POOG 0 10000 0 01368 0 00000 PNT POOL 0 07500 0 01090 0 00000 PNT POOJ 0 05000 0 00770 0 00000 PNT POOK 0 02500 0 00420 0 00000 PNT POOL 0 01250 0 00230 0 00000 PNT POOO 0 25000 0 00000 0 00000 PNT POOP 0 24688 0 00308 0 00000 133 DDD ED DD DI SS DD IS SIS POOR POOS POOT POOV POOW POOX P00Z P010 P011 P012 P013 P014 P015 P016 P017 P019 P1 0 0 0 0 0 0 0 Q sQ Q 0 0 0 0 0 0 wS TT TT ARATE ASST UE A S 7 LG A LE es ll Lo LL LLE a Upp Ms FEY S i Perry ni QUE ne Jet LLIT TIC haa i ZI Cy GH 7 CI UTG hh RS RS V L WA A NO SS 7 Wu
24. Dataset 1 are written For frd format the ds pa rameter w o any further values forces cgx to write all datasets but w o the internally calculated ones send all frd ds If a sign appears at the end of the ds parameter send all frd ds or send all frd ds1 then also the internally calculated values like vMises stress are written Ad ditionally the entities might be specified with the e parameter followed by the numbers to write The numbers have to be separated by a sign or in case of a range by a sign send all frd dsl e9 25 30 The bin parameter has the same meaning as the ds parameter but in this case the result is written in the binary form of the frd format always all datasets In case data should be written for the crack analysis tool Darwin send all dar ds will write all datasets Since Darwin changed the format from version 7 1 on it is possible to change the format with the parameter v7 1 as the last parameter If the set is of the ordered type and includes nodes see fasegj then the data are written in tabular form for use in a 1D crack prop Darwin analysis In certain circumstances the user needs an easy solver independent format to write node and element numbers This is provided by the Ist parameter send setname Ist will write the node and element numbers in lines of seven space separated rows Missing lower entities nodes points lines etc will be a
25. LL LL LL LL LL B 15 Remarks Concerning Ansyg LL 00000 2G B 17 Remarks Concerning dolfyn 2 0004 97 97 98 104 105 105 105 106 107 107 110 110 110 111 111 111 111 112 112 B 20 Remarks Concerning NETGEN 125 ns 4G Pie oe Sie O 125 cani hue eae eas 126 128 ul des p mere ean E 128 ee eee ee ee eee eee 129 Tents sua eins NA fl puke ly Gee aut gs A hare gece 6 131 CA Sphere Volume o 132 C5 Arrtoil tor Cid coded s sra dE eo se ts OY Bo ee ea 133 C 6 Tf and while demad o ss di s es a Ga A ego OS ee Gr a 138 1 Introduction This document is the description of CalculiX GraphiX cgx This program is designed to generate and display finite elements FE and results coming from CalculiX CrunchiX ccx Tf you have any problems using cgx this document should solve them If not you might send an email to the author 3 The Con cept and File Format sections give some background on functionality and mesher capabilities The Getting Started section describes how to run the verification examples you should have obtained along with the code of the program You might use this section to check whether you installed CalculiX correctly Then a detailed overview is given of the menu and all the available keywords in al phabetical order in the Menu and Commands sections respectively Finally the User s Manual ends with th
26. This keyword is used to close an open set Without parameter setc will close the at last opened set 8 101 sete sete lt set gt in e p 1 1 c 1 8 1 b 8 gt L 1 30 sh uN max min strict This keyword is used to enquire other sets which have entities in common or which have an identical content Example sete blade p min searches for sets who store at least the same points as blade always the set all The option max searches for sets whose entities are completely included in the specified set The option strict searches for identical sets regarding the specified type of entity 8 102 seti gseti lt set gt ar Pe Pp PE e PSs PEPE PE Pse eh gt lt set gt lt set gt This keyword is used to generate an intersection of certain sets what do the specified sets have in common of the specified entity type Example seti intersectSet p setl set2 set3 generates a set intersectSet with points which are stored in each of the sets setl set2 set3 8 103 seto seto lt set gt This keyword is used to enquire open sets seto Or to mark a set as open seto set All newly defined or redefined entities will be members of all open sets See sete how to close an open set 88 8 104 setr gsetr lt set gt mn e p PI sP PE 1S v se sh lt name chars gt lt name chars gt This keyword is used to remove entities from a set see also The entity
27. be moved to their nearest neighbour in set trgt set as long as the neighbour is not more than 0 01 units away move slave set equ trgt set 0 01 mir will mirror the set The mirror plane is placed normal to the direction running from P1 to P2 and placed at P2 or defined by a point and a main axis x y z as shown in the following example move part mir P y 8 54 movi movi delay lt sec gt start stop make lt pic nr gt lt pic nr gt lt prolog gif gt clean This keyword is used to start or stop the recording of a movie After start all frames will be stored in single gif files until the stop command is issued See also the menu options Start Recording Gif Movie Use the option make to assemble the movie from the individual files An existing movie will be copied in front of a range of frames if it is specified as the last argument The range consists of the nr of the first and last picture to be used With the option delay a time delay in seconds between frames can be specified With the option clean all single gif files will be erased Here is an example command sequences movi delay 0 01 movi start movi stop movi make 0 500 movi make 0 500 prolog gif movi clean Further remarks in How to change the format of the movie file 49 8 55 msg msg onloff This keyword is used to enable or disable full printout during runtime This is useful for debugg
28. between int float and exponent See the command ffont on how to change the font size To catch more than one entity with one stroke type a all before and create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the rectangle To catch only the entity which is closest to the user type i before DAT3 STRESS Time 100 000000 Entity Mises 71e 02 32e 00 71e 02 25e 02 79e 02 33e 02 87e 02 40e 02 94e 02 48e 02 02e 02 56e 02 09e 02 63e 02 17e 02 71e 02 25e 02 78e 02 32e 02 86e 02 40e 02 37e 01 75e 01 32e 00 max min e A DV PRR WW BP KUOOTATATWANN OW WW FO latimsta frd Figure 8 qtxt Node attached texts at two locations one with node nr and value in exp form and one just with his value in floating point form The font with gives biggest numbers was used 8 89 quit quit This keyword is used to terminate the program without save 69 8 90 read read lt geo file gt add lt ccx file gt inp add ext nom lt openFoam file gt foam add ext nom lt result file frd gt add ext nom lt setname gt lt stl file gt add ext nom lt list file gt n e lt column gt lt edge file gt edg lt netgen file gt ng lt pixel file gt l
29. but instead of releasing the left mouse button over a Dataset continue to the bottom of the sub menu to Entity Continue from that item to the right and release the mouse button when the pointer is over an entity Now the data is displayed in the drawing area 17 7 2 Viewing In the following sections changing properties and styles of the displayed struc ture are explained See the command view to control the functions with the command line 7 2 1 Show Elements With Light This is the default view of the mesh if the program was started in viewing mode If used any animation will be interrupted and no values are displayed 7 2 2 Show Bad Elements This option presents elements which have a negative Jacobian value at least at one integration point The solver cecx can not deal with those elements So far only TET and HEX elements are checked These elements are stored in the set called NJBY See also the command Boel 7 2 3 Fill This is the default mode and forces the element faces to be rendered 7 2 4 Lines The edges of the element faces are displayed This is especially useful to see into the structure to find hot spots in the displayed field With Toggle Move Z Zoom and gcutl a more detailed analysis can follow For very dense meshes switch to Dots 7 2 5 Dots The corners of the element faces will be displayed This is especially useful if values inside the structure need checked 7 2 6 Toggle
30. c d code 6 isaac a free cfd code 7 e OpenFOAM a free cfd code 8 e Nastran the f06 file can be read sf only CHEXA displacements and stresses For a more detailed description on how to use cgx to read this formats see Program Parameters and the program specific ips and Hints sections See he t command for how to write them from cgx 4 Getting Started For installation help see Calculix cgx X X INSTALL After the program is installed on your machine you should check the functionality by running the examples included in the distribution The examples are located in Cal culix cgx X X examples Begin with a result file called result frd Just type cgx result frd and some information is echoed in the xterm and a new window called main window appears on the screen The name conventions used for the different areas in the main window are explained in figure Now you should move the mouse pointer into the menu area and press the left mouse button Keep it pressed and continue over the menu item Dataset to Disp There you re lease the button Then press the left button again and continue over Dataset and Entity to D1 For background informations look into the subsection Datasets and Entity which explains how to display results After seeing the values you might play around a bit with the Menu Before going further you should read the section Input Devices
31. deleted The nodes are stored in a new set called delete and the user might ani this set 8 50 min min lt value gt f i e Pr Pu This keyword is used to define the lower value in the scale in the menu area see figure 2 The number representation can be changed between int float exp For example max 1100 i will set the upper value to 1100 and the representation to integer A third ga parameter 1 lock or u unlock can be provided which locks the scale to certain max or min values The selection of a different dataset will not change the scale 8 51 minus minus mn e Pp PIs Ib 17S 1 L sh lt set gt This keyword is used to remove entities of a set from the screen see also plus The following entities are known Nodes n Elements e Points p Lines l Surfaces s Bodies b Nurb Surfaces S Nurb Lines L and Shapes sh Only the set which was used to display the entities can be used to remove them AT 8 52 mm mm lt value gt PP Pi Pel 1 u This keyword combines the functionallity of the commands max and fminjin one command The minimum value is set to max 8 53 move move lt set gt scal lt fx gt lt fy gt lt fz gt lt pnt gt tra lt dx gt lt dy gt lt dz gt Prot lt pi gt lt p2 gt 1 x 1 y 1 z Px ylz lt alfa gt lt alfa1 gt lt ax1 gt lt alfa2 gt lt ax2 gt Prad lt pi gt lt p2
32. in the result format frd then a so called list file could be read instead In such a file either numbers of nodes or elements can be defined The numbers found in a spec ified column is interpreted as a node or element number If n is specified the numbers are interpreted as nodes and if e is specified the numbers are interpreted as elements For example 71 read list file e3 will define a set storing names of elements from the third column of file list file txt The created set name is always the name of the file In case NETGEN is used for meshing then the edges which are generated by NETGEN can be included for modifications read edge file edg will create beam elements based on the defined edges When the modifications are done the beam elements can be exported in the NETGEN edge format with the send setname stl command and used for meshing in netgen see also How to deal with cad geometryl The full netgen mesh can be imported with read netgen file ng Be aware that all netgen generated element types are imported 1D 2D and 3D elements This elements are distributed in various sets for further opera tions If a pixel file in xwd format is available it can be included as a background picture The user can then create geometry based on this picture read pixelfile 2 4 here the picture will be scaled by a factor of 2 in x and 4 in y direction Delete it with del pic The picture can
33. lines lemb c A for surfaces s B for Bodies b Q for nurb lines nurl L N for nurb surfaces nurs S A for sets se and H for Shapes sh For example asgn p U will assign the character U as the leading character to all newly created names of points The automatically created names of geometric entities use 4 charac ters If all possible names with the chosen leading letter are in use then the next alphabetical letter is chosen as a leading letter so after PZZZ follows Q000 If no more letter follow then the amount of letters per name is increased The maximum number is 8 Each entity has its own name space Different entities might use the same name Remark Currently nurbs lines are automatically used to create splines sharing the same name Nurbs lines can not be used for other purposes than to be displayed and so far they can not be written to a file The command is also used to controll the behaviour of the unstructured triangulator This unstructured mesher 14 uses the tree parameters alpha beta nadapt for mesh control Current default is 0 4 for alpha and beta and 4 for nadapt In case Nastran input should be generated it is possible to switch from MPCs to RBEs when using the send command in combination with the areampc op tion The value after rbe represents the thermal expansion coefficient of this elements 25 asgn rbe 0 5e 6 It should be noted that coincident nodes are connected by MPCs either way 8 4 bia
34. lt set gt lt target set gt lt shpe gt gt tra lt dx gt lt dy gt lt dz gt lt offset gt rot lt p1 gt lt p2 gt lt offset gt Prot x y z lt offset gt rad lt p1 gt lt p2 gt lt offset gt Prad x y z lt offset gt nor lt offset gt lt tol gt measures distances between entities of two sets For example between points or nodes stored in one set to surfaces or shapes stored in a second set The average maximum and minimum distance is determined The distance is measured normal rotational radial or translatoric The command works analogous to the proj command Please look there for details The command writes to the tack 31 8 16 div div lt line gt lt division gt lt set gt lt division gt mult div lt factor div gt lt factor bias gt auto lt node dist gt lt angle gt lt elem ratio gt This keyword is used to define the division of line s The div keyword also works on a set of lines seelgadd The division controls the number of nodes created when the geometry is meshed see and mesh For example div all 4 attaches the division of 4 to all lines With the keyword mult or div in combi nation with a value it is possible to multiply or divide already assigned divisions div all mult 2 Or in case you need a starting point for the individual divisions you can use the option aut
35. nodes from the set in dep In this case the equations will be written in the nastran format nas in nastran called MPC The parameter rx12 defines the displacement system as rectangular r the rotational axis is x and the 12 defines the num ber of segments in 360 deg therefore the angle of the segment is 360 deg 12 Attention The sign of the number of segments must be negative if the angle between the independent side and the dependent side is negative When the nr of segments is omitted the value is calculated individual for each node The c triggers the correction of the position of the dependent nodes to a position defined by the angle of the segment highly recommended u would prevent the correction The 1 will be the identifier for the equations if the format is nas nastran In case the format would be ans ansys then the 1 would be the number of the first equation No number is required for abq abaqus and calculix Tf the coordinate system of the displacements is cylindrical rtz then the example would be 82 send dep indep nas cyempce cx12 cl Only the r from rx12 is changed to c A thermal connection is cre ated with send dep indep nas cycmpc tx c The thermal connection is triggered by the t a pressure connection with p See comments above for the single parameters In case of cfd elements the syntax is similar but only the function na
36. not be modified in cgx Only scal ing during reading with either a global factor or separate scaling in x and y direction is supported Other modifications have to be made with external software 8 91 rep rep lt set gt This keyword is used to refresh entities of a certain set This is done automati cally but with one exception The rendering of the interior definition of surfaces must be triggered manually with this command The rendered surfaces can be displayed with plot or plus using the type si 8 92 rnam rnam lt set gt lt set gt This keyword is used to rename a set 72 8 93 rot rot u d r 1 c lt angle gt x P 3 xo y y 17 2 7 This keyword is used to rotate the view direction to the model For example rot u 10 will rotate the view direction 10 degrees upwards The meaning of the other letters is down d right r left 1 clockwise c The remaining letters will rotate the view into a pre defined direction For example the user will look in x direction after the command rot x and against the x direction with the command rot x and so on 8 94 save save This keyword is used to save the geometry if any to a file named as the input file but with the extension fbd If a file with that name exists already then this file will be saved with the new file extension fbb as a backup 8 95 scal scal s v d lt value gt This k
37. png format The program convert must be installed 7 11 5 Start Recording Gif Movie All frames during an animation are stored The recording ends after the right mouse button is pressed while in the menu area Finally all frames are combined in the file movie gif which can be displayed with various tools Firefox I2 or realplay Tf the animation is stopped with the middle mouse button while in the menu area then the movie stops until it is released by pressing the middle mouse button again See movil for the keyboard options Further remarks in How to change the format of the movie file 7 12 Help Starts the html help and displays this document It only works if the specified html viewer is available The default is Firefox P2 but this can be changed in the cgx h file The search path for the documentation is also defined in the cgx h file Please make sure that the documentation is in the speci fied location or change the path in the cgx h file and recompile the pro gram after the object files are deleted The default location for the html help is CalculiX cgx X X doc cgx and CalculiX ccx_X X doc ccx for cgx and cecx respectively The html files must be downloaded directly or compiled from the latex source for this function to work properly The INSTALL file tells how to compile the latex code to html The INSTALL file is located Cal culiX cgx_X X and CalculiX ccx_X X for cgx and cex respectivel
38. strained in the related DOF number e dep BOU Created if EQUATION is found All dependent nodes are stored in dep e ind BOU Analog to dep e clo nr Created if CLOAD is found The nr counts the number of appearance e dlo nr Created if DLOAD is found e mpc nr Created if MPC is found e rad nr Created if MPC is found e fim nr Created if FILM is found e cflx nr Created if CFLUX is found e dflx nr Created if DFLUX is found e tie l 2 Created if TTE is found The set tiel stores the slave entities and tie2 the master The single sets which define the contact areas are linked together so if identified with the qenq command the referenced linked opposing set is listed in the xterm e trans Created if TRANSFORM is found The name consists of the definition of the transformation e elty Created for all known element types like C3D20 Check the element quality with Show Bad Elements or use egall to set thresholds and plot to actually plot the affected elements Forces regard their referenced coordinate system TRANSFORM The val ues are automatically transformed into the global cartesian system so that the vectors point in the correct direction Multiple load definitions inside one STEP on a single entity will sum up ccx compatible This applies to cflux dflux cload dload 122 B 15 Remarks Concerning Ansys The cgx is capable to write the fo
39. such a way you can create and modify geometry or meshes in batch mode or evaluate results in batch mode But you need a graphic capable computer because cgx needs it also if no graphic output is requested Nevertheless the popup of the graphic window can be suppressed if cgx was started with the bg parameter Some commands make only sense in a batch file like iff and while Please see also the demo in Land while demol on how to work with this commands See also THow to write values to a file It should be noted that a successive save or exit command will overwrite the batch file if it has the file extension fbd To prevent this the user should use a different extension like fbl B 13 How to deal with cad geometry In general hexahedra elements perform better than tets but if the mesh should be derived from a cad geometry it is often more convenient to create a tetrahedra mesh as to modify or rebuild the geometry to make it meshable with hexahedra elements The following section gives some hints which alternatives are available e A simple step reader is integrated in cgx cgx step filename It can deal with points and lines This is sufficient to create simple hexaeder elements for axis symmetric structures like a rotor but the experienced user might use the following approach also for more complex geometries You might start with a file containing a 2D section subdivide it in mesh able surfaces and sweep it i
40. the user wants to replace it by a spline line then the user defines or selects the new spline and then moves the mouse in an area were one of the old lines the new one might be present as well is located and presses the e key The definition of the lemb will be changed without destroying the surface definition Press q to quit the command or u to delete the last created line It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 79 qmsh qmsh RETURN f v gt n gt 1 c s m 1d 1 q A very usefull command to optimize surface meshes as a basis for tet meshes Based on identified lines I or surfaces s meshes can be deleted d and generated m Surfaces can be combined c to get rid of very small ones usually first identified by a common line with I Example Identify a line with key I The related surfaces are selected e Then press key d if the mesh retated to this surfaces should be deleted e The line divisions might be changed by pressing the key with the desired number This function of qmsh works like the Tgdiv command e Press key m to re mesh the surfaces e To combine two selected surfaces press key c after selecting surfaces individually or by selecting a line 64 A sur
41. to be selected For example you type send sur abq slide s and this should attach the node with the node nr 1 to a plane parallel to the surface of the corresponding element face This element uses the node 1 2 3 4 5 6 7 and 8 where 1 2 3 and 4 are on the surface of the mesh The program will detect this element face and constrain the node only then if the nodes 2 3 and 4 were also included in the set sur The detected elements or faces of volume elements will be assigned to the set So it is possible to check the detected elements Instead of providing a set with the necessary nodes you might instead provide a set with shell elements or and faces of volume elements The nodes will be derived from them In most cases this is the most convenient way A special case which works probably only in ccx abaqus might give wrong results is triggered by the c parameter Tt behaves as for the s parameter but works for cycsym calcs send cycsur abq slide c A new set with new nodes will be created with the name lt set gt COPY which has to be used in the cycsym definition in ccx instead of lt set gt In addi tion equations are specified which connect the new nodes with the original ones in the normal direction of the element faces In this way two faces used in a cyclic symmetry connection will slide at each other instead of beeing glued in all directions If nodes should slide radially to the y axis then use
42. types Values NODE XX Where NODE XX node number or blank if KEY 2 data 109 In case of material dependant data REMARK Implemented only for NMATS 1 first line Short Format 1X 12 415 Long Format 1X 12 110 315 Values KEY NODENR NMATS Where KEY 1 NODENR Node number NMATS Number of different materials at this node unused second and following lines Short Format 1X 12 15 6E12 5 Long Format 1X 12 110 6E12 5 Values KEY MAT XX YY ZZ XY YZ ZX Where KEY 2 MAT material property number if KEY 2 unused XX data Last Record only FORMAT 0 1 ascii omitted for FORMAT 2 Format 1X 3 Values KEY 11 Pre defined Calculations Listing of the automatically calculated additional results 11 1 Von Mises Equivalent Stress Entity name Mises 1 ToM Fy los Oy oy Oz 0z 09 672 672 672 11 2 Principal Stresses Entity names P1 P2 P3 worstPS The principal stresses o are named P1 P2 P3 From the three principal stresses o the absolute maximum value will be calculated and named worstPS For ex ample if a node has the three values 100 0 and 110 MPa then 110 MPa would be shown The three principal stresses o4 09 03 are derived from the following 110 equation Ors A Oxy Gee nx 0 yx Oyy A Oyz ny 0 Ose Oxy Or X nz 0 They are given by the three roots of the equation stress tensor is symmetric Ory Oya etc
43. use this command mencoder lt inputfile gt ovc lavc lavcopts vcodec flv of lavf In general To convert a video file from one format into the other use the program mencoder It comes with the mplayer packages mplayer movie wmv B 2 How to define a set of entities Some knowledge is necessary to efficiently select entities nodes points Asa golden rule in complicated situations never try to create a set just by adding the entities to the set see qadd It is much better to catch a bigger group which includes certainly the wanted ones Then display the set with the command and remove all unwanted entities with the qrem command Entities are selected if at least one pixel of it is inside the selection rectangle But there 112 is one exception Only the lower left corner of any text names of points etc can be selected Surfaces can also be selected by picking its shaded interiour see Tep You might add all lower entities as points for example by typing comp set do This is necessary for example if you had selected surfaces and you want to move them in space Only points have locations and therefore nothing will happen unless you completed the set by the related points with comp command There is another type of set called sequence The data structure is the same but with one exception the data keep their order in which they were selected This type of set is used for splines see glin and
44. with 3 amplitudes and 3 phase angles in with 6 amplitudes and 6 phase angles in ICIND1 sub component index or row number 108 ICIND2 column number for ICTYPE 4 IEXIST 0 data are provided 1 data are to be calculated by predefined functions not used 2 as 0 but flaged by cgx ICNAME Name of the predefined calculation not used ALL calculate the total displacement if ICTYPE 2 This record must be repeated for each entity 4 Type of Record not used This record will be necessary in combination with the request for predefined calculations This type of record is not allowed in combination with binary coding of data Format 1X 12 215 2013 Values KEY IRECTY NUMCPS LSTCPS I I 1 NUMCPS Where KEY 6 IRECTY Record variant identification number NUMCPS Number of components LSTCPS For each variant component the position of the corresponding component in attribute definition 5 Type of Record The following records are data records and the format is repeated for each node In case of material independent data ascii coding Following records ascci FORMAT 0 1 Short Format 1X 12 15 6E12 5 Long Format 1X 12 110 6E12 5 Values KEY NODE XX Where KEY 1 if its the first line of data for a given node 2 if its a continuation line node number or blank if KEY 2 data NODE XX binary coding Following records ascci FORMAT 2 int NCOMPS xfloat int and float are ansi c data
45. 0 For dolfyn duns and isaac the same syntax has to be used send all duns viscous wall profil subsonic inflow in subsonic outflow out 84 send all isaac WALL profil FARFIELD far EXTRAPOLATE out send all dolfyn INLE ingang OUTL uitgang etc To define so called gap elements and related control commands These elements will connect parts if they are closer as a certain distance For example if the distance is zero contact The command send dep indep abq gap 1 0 0 will connect the nodes in the set dep and indep with gap elements but only if they match each other in the direction 1 8 97 sega sega lt seq gt nod pnt lt name gt lt gt afte befo lt name gt nod pnt lt name gt lt gt l end nod pnt lt name gt lt gt This keyword is used to create or redefine a set marked as a sequential set This set is used for spline definitions see ine With the command qlin such a sequential set is automatically created To begin such a set type for example SEQA Q003 PNT P004 P005 P006 POOM POON The program will create or overwrite the set Q003 The command will con tinue in the next line if the sign is found SEQA Q003 PNT P004 P005 P006 P007 P008 P009 The parameter AFTE will insert additional points after the first specified point in the existing sequence The parameter BEFO will insert additional points be fore the first specifie
46. 34 gsur gsur lt name char lt 9 gt BLEND lt nurbs gt lt shape gt gt lt line lcmb gt lt line lcmb gt 3 5 times This keyword is used to define or redefine a surface in the most basic way see also qsur Each surface which is intended for a regular mesh must have three to five edges which might consist of a single line or a combined line see Icmbj to be mesh able However the recommend amount of edges is four For example gsur 5004 BLEND L002 LOOE L006 LOOC will create the surface S004 with a mathematically positive orientation indicated by the sign after the surface name The keyword BLEND indicates that the interior of the surface will be defined according to Coons or a NURBS surface or shape is referenced It should be mentioned that only nurbs or shape related surfaces can be meshed with unstructured triangles Use a or in front of the lines or lembs to indicate the orientation These signs will be corrected automatically if necessary If the surface is intended for an unstructured mesh it is not necessary to care about the correct orientation of the lines and the number of edges is not limmited If automatic name generation is desired then use instead of a name 8 35 gtol gtol lt geometric tol gt This keyword is used to enquire the default geometric tolerance gtol Or it is used to set the default geometric tolerance gtol
47. A SETA SETA SETA SETA SETA SETA SETA SETA SETA SETA SETA SETA SETA SETA LO6H LO6I LO6J L001 LOOA LOOK L009 L000 L002 LOOL PO5G PO5L PO5H PO1A PO3E P052 PO2M P03C POON PO3P PO5 PO5 POS POO P03 P03 PO1 POS P03 P02 I H J N P B B 4 E M cl pli p12 120 L005 P000 P019 C001 C004 A001 A002 A003 A004 A005 AOOI A00J AO0OK AOOL AOON A000 AOOP AOOR A0O0S AOOT wall wall wall wall prof prof prof prof prof prof prof in 1 out out out out out E L BL BL BL BL BL BL BL BL BL BL BL BL BL BL BL LO LO LO LO PRPPe EE EE A do do do te t tt il ik il il il il il I L004 1 LO5 1 LO5 1 LOO 1 LOO 1 LO5 1 1 1 I 1 1 001 OOK END END END END END END END END END END END END END END END 5C 59 SA SD LOOC LOOE L001 LOOA L009 L002 LOOL F S G I I 204 130 204 S001 SOOE SOOL SOOA S00X S002 SOOP 120 L002 LOON L003 L006 LOOC L008 LO6C LO4W LO4V L04Z L06J LO4V LO6H LO4W LO4X L058 L056 210 120 110 210 110 205 205 LOOA L000 C004 LO3S LOO7 LO5Z LOOL L050 L004 LO5C LO5D L052 LO6F LO3R L054 LO5V LO6I 137 LOOL L006 L008 LO3R LOOG LOSV LO4X L04Z L06J LO5F LO6H LOSI L058 L056 LO5L LOSS L009 cool LOOE L003 LO6C L007 C004 L050 L054 LOGI LO5A L059
48. CalculiX USER S MANUAL CalculiX GraphiX Version 2 9 Klaus Wittig November 14 2015 Figure 1 A complex model made from scratch using second order brick elements Contents 2 Concept 6 Input Device co E 7 2 8 Toggle Element Edges 7 2 9 Toggle Surfaces Volumes 7 2 10 Toggle Move Z Zoo 7 3 3 Steps per Period 7 3 4 Time per Period 8 27 fipa 37 8 28 fon 37 92 8 116ucut 93 94 95 96 8 120while 96 9 Element Type 10 Result Forma 10 1 Model Header Record 12 User Function A Known Problemas A 1 Program is not responding B_Tips and Hint B 1 How to change the format of the movie file B 3 How to enquire node numbers and values at certain locations B 5 How to select only nodes on the surface B 12 How to run cex in batch mode B 14 How to check an input file for ccx B 16 Remarks Concerning Code Aste B 18 Remarks Concerning Duns and Isaac B 19 Remarks Concerning Nastra 10 2 User Header Record o o o o o A 2 During Meshing se s e sa ana s da a B 2 How to define a set of entitieg B 4 How to write values to a filg B 6 How to generate a time history plot B 11 How to map load B 13 How to deal with cad seometryl
49. Culling Back Front This removes the faces of volume elements for all elements or for the surface of the structure depending on the state of Toggle Surfaces Volumes With this option the user can visualize internal structures like cracks or a core of a hollow structure For shell elements the behavior is different Initially only the front faces are illuminated and the back faces are dark This is helpful to determine the orienta tion of the elements If you want to see all faces of the shell elements illuminated regardless of the orientation then use this option If you want to change the bJ orientation use the command Jafipp 18 7 2 7 Toggle Model Edges Per default all free element edges are shown The user can remove show them with this option 7 2 8 Toggle Element Edges Per default just the free element edges are shown The user might add all edges to the structure with that option 7 2 9 Toggle Surfaces Volumes This switches the way each volume elements are displayed Either all faces of the elements or just the element faces on the surface of the structure are displayed Depending on the state of Toggle Culling Back Front either the faces pointing to the user or the faces pointing away are displayed The default is just to show the surface pointing to the user In the lower left corner of the drawing area see figure P a character is printed indicating the program is in the surface mode s or in th
50. L Like Rendering Toolkit from Brian Paul http www mesa3d org OpenGL Utility Toolkit GLUT from Mark J Kilgard CalculiX GraphiX cgx from Klaus Wittig klaus wittigQ calculix de NETGEN unstructured mesher http www hpfem jku at netgen dolfyn Open Source CFD code http www dolfyn net Duns a two and three dimensional cfd code http sourceforge net projects duns ISAAC a two and three dimensional cfd code http isaac cfd sourceforge net OpenFOAM a three dimensional cfd code http www opencfd co uk Tochnog a free fem code http tochnog sourceforge net Tutorial for CalculiX from Dr Guido Dhondt http www dhondt de tutorial html ImageMagick 5 1 0 00 01 01 Q 8 cristyg mystic es dupont com Copy right Copyright C 2000 ImageMagick Studio Mozilla Foundation http www firefox com S A Coons Surfaces for computer aided design of space forms Project MAC MIT 1964 Revised to MAC TR 41 1967 mesh2d unstructured 2D mesher from B Kaan Karamete Ph D No URL Available Paul Dierckx Curve and Surface Fitting with Splines Oxford University Press 1993 140
51. L052 LO3S LO6F LO5Z SETA out 1 LO5L C 6 If and while demo Paos P201 P3014 P401 P411 144 Pza P3t P paz P2214 P321 P421 paza P2314 P331 paor 300 P400 pioo P200 P P110 P210 P310 P410 P120 P220 P320 P420 pazo P230 P330 P430 batchfile fbl if amp value amp while demo Figure 18 Result of If and while demo The if and while commands can be nested A demo which produces some points on the window follows text if amp value amp while demo def the leading letter of point names valu vp P define the initial x value valu vx 0 define parameters valu v2 4 valu v3 1 start loop while vx lt v2 valu vy 0 138 valu vz 0 seto S1 while vy lt v2 define the pnt coordinates valu vy vy v3 valu vy int vy valu vx int vx valu vz int vz define the pnt name valu pi vp vy valu pi amp pl vx valu pi amp pl vz generate the pnt pnt pl vx vy vz endwhile setc valu vy 0 valu vz 1 seto S2 while vy lt v2 valu vy vy v3 valu vy int vy valu vx int vx valu vz int vz valu p2 amp vp vy valu p2 amp p2 vx valu p2 amp p2 vz pnt p2 vx vy vz endwhile setc valu vx vx v3 if vx plot pa Sir else plus pa 52 b endif endwhile send all fbd demo on how to write the content of variables to writedemo txt valu vx int vx valu vy int vy sys echo VX vx VY vx VZ vz gt writedemo txt 139 References 1 10 11 12 13 14 15 OpenG
52. See also the commands steps max min or the combination of max and min Tm and Iscall which might be used to modify the colour representation of the displayed values For example type min 0 to set the lower value of the colour bar to zero Now you should study the following interactive commands Use qera to enquire values at nodes Use gtxt to generate node attached texts showing their number and value Use Jqcut to generate a section through the model And use leraph to generate a 2D time history plot for results with several time steps or a 2D plot of values along a sequence of nodes see Jasedj Watch out if you type a command the cgx window MUST stay active and not the xterm from which the program was started It is better to stay with the mouse pointer in the cgx window Next Quit the program and type cgx b geometry fbd in the xterm The program starts again but now you see only a wire frame of the geometry Move the mouse pointer into the new window and type Imeshi all The mouse pointer MUST stay in this window during typing and NOT in the xterm from which the program was started After you see ready in the 10 parent xterm the mesh is created To actually see it type plus ea all Now you see the mesh in green color To see the mesh as a wire frame choose in the main menu Viewing and continue to the entry Toggle Element Edges and then again in Viewing choose Dots To se
53. ad Nastran f06 file ng read Netgen native format with surface domains step read an ascii step file points and lines only stepsplit read step and write its parts to the filesystem in separate directories stl read an ascii stl file triangles v default read a result file in frd format and optional a solver input file ccx in addition which provides the sets and loads used in the calculation special purpose options mksets make node sets from DLOAD values setname _ lt value gt read forces the program to read the complete result file at startup If no option is provided then a result file frd is assumed see Result Format A file containing commands or geometric informations is assumed if the option b is specified Such a file will be created if you use exit or Bavel after you have interactively created geometry Option a awaits the same format as option b but merging defining of line divisions and the calculation of the interior of the surfaces is done automatically and the illuminated structure is presented after startup This should be used if the commandfile was generated by an interface program which convertes cad data to cgx format for example vda2fbd With option a and b the program will start also if no file is specified An input file for the solver can be read with option c Certain key words are known and the affected nodes or elements are stored in sets For example the def
54. all gt file txt See also if and while and see the demo in T and while demo on how to work with this commands They are usefull for automatic post processing 113 B 5 How to select only nodes on the surface Some times you need to select nodes only on the surface of the mesh This can be done when the mesh is displayed in the surface mode see Toggle Surfaces Volumes using the menu entry Show Elements With Light To find node positions more easily add the element edges to the view see Toggle Element Edges A different way uses the element faces see plot with parameter f Store the desired faces in a new set qadd and complete the set downwards comp set do Check the selected nodes with plot n setname B 6 How to generate a time history plot So called time history plots can be created based on a previous displayed se quence of data sets see Toggle Dataset Sequence with the command graph It is also possible to use only the command line For example edsle3 e ds 1210 e graph set t will produce a time history plot for the nodes stored in set over all loadcases from 1 to 10 for the entity nr 3 Or e graph set t DISP D1 will display the displacement in direction 1 for all loadcases For more details and other options look into the graph command description 114 40 FP AAA AAA AE A TEA EE I ON Lo Le Jig N Il o he e Z o F qa EL 48 g E E
55. ance value which restricts the deviation of the normal vectors of faces from the selected face As long as the deviation is below the specified value in de grees all adjacent faces will be selected in a loop qadd areaset t25 IF pit Press q to quit the command It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 66 qali qali RETURN gt w Pp PR Pq This keyword is used to align a plane defined by three points or nodes with the screen working plane This is useful if a point or a node should be moved manually along a defined plane see qpnt To define the plane move the mouse pointer over the first entity and press n if its a node or a p if its a point Then define the next two entities in the same way Press q to quit the command It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 67 qbia qbia RETURN aw P Pir Pepee 10 gt 99 q This keyword is used to change the bias of a line The bias defines a coarsening or refinement of the mesh along a line The number defines the ratio of the length of the first element to the length of the last element at a given line It works by pressing a nu
56. are NOT marked in the dataset as 19 vectors can be displayed with vector needles This command line approach with ds is the only way to display duns cfd results with vector needles See also the command cal how to manipulate the length of the vectors 7 2 13 Toggle Add Displacement It is possible to display results on the deformed structure For example you can display a stress field on the deformed structure If you know a suitable amplification factor for your displacements then use the call command to issue this value but this can also be done later Of course displacements for the Loadcase must be available 7 3 Animate This option allows the animation of displacements See also Tanim ds and call to use this functionality with the command line It is possible to create this sequence from just one Dataset see IStartl This is useful for displaying mode shapes See also Toggle Dataset Sequence to create a sequence from multiple Datasets to visualize dynamic responses 7 3 1 Start Creates a sequence of display lists to visualize displacements for example mode shapes The program recognizes displacements just by the name of the dataset This name must start with the letters DISP otherwise the animation will not start see Nodal Results Block 7 3 2 Tune Value Controls the amplitude of the animation If Toggle Real Displacements was chosen before the tune value is equivalent to the amplifi
57. ault set s bou dof store nodes which are restricted in the corresponding degree of freedom and the set s dep dof and ind dof store dependent and independent nodes used in equations A special case is OpenFOAM The results are organized in a directory struc 14 ture consisting of a case containing time directories in which the result files are stored The user must call cgx using the case directory cgx foam case The program will then search the time directories The time directories must con tain a time file to be recognized Or in other words each directory in this level containing a time file is regarded as a result directory 6 Input Devices 6 1 Mouse The mouse is used to manipulate the view point and scaling of the object inside the drawing area figure 2 Rotation of the object is controlled by the left mouse button zoom in and out by the middle mouse button and translation of the object is controlled by the right mouse button Inside the menu area the mouse triggers the main menu with the left button In addition the mouse controls the of nodal values The animation will stop if the mouse pointer is not in the drawing area but will start again if the pointer enters the drawing area This can be prevented by pressing the middle mouse button while the mouse pointer is in the menu area Pressing the right button will release the next frame A frozen animation can be released by pressing the middle button The previous fra
58. c e Equations or RBEs areampc to glue components A previous command fasgn defines if mpcs or rbes will be created e RBE2 mpc for rbe spiders e Pressure pres so far only CHEXA8 e Temperatures ds all datasets with just one entity will be written as temperatures The resulting files have to be combined with the help of an editor and extended by material data and the necessary controll commands The f06 file with results can be read so far only CHEXA with displacements and stresses B 20 Remarks Concerning NETGEN It is not necessary to write the mesh in abaqus format if you use netgen as a mesher The native netgen format vol can be read by cgx cgx ng file vol as well The netgen mesh format vol includes also the surface patches which were defined by the edges of the model and used for the generation of the volume mesh This patches can be used to define boundary conditions or loads The nodes and faces of this patches are stored in sets named set nr To get an overview over the patches type prnt se To see were the patches are lo cated type fplotjf all and use the PAGE DOWN key to scan through all sets A netgen surface mesh can be written based on the faces of elements The faces of hex tet quad and tria elements are triangulated and written in the stand alone netgen mesher format file ng This mesher can be found in the netgen sub directory nglib and is named ng vol It will create t
59. cation of the animation 7 3 3 Steps per Period Determines how many display lists for one period of animation will be used If Toggle Dataset Sequence was chosen then these number of display lists will be interpreted as one period see Time per Period 7 3 4 Time per Period Determines how many seconds per period 7 3 5 Toggle Real Displacements To see the correct displacement of each node The animation can be controlled with the help of the Mouse 20 7 3 6 The user can switch on additional undeformed model edges This is usefull for hardcopies were this edges give a reference to the undeformed shape 7 3 7 The user can switch on additional undeformed element edges This is usefull for hardcopies were this edges give a reference to the undeformed shape 7 3 8 Toggle Dataset Sequence Creates a sequence of display lists to visualize values of a sequence of Datasets The Datasets must use the same type for example only displacements or only stresses To activate the animation after you have selected Toggle Dataset Sequence choose the first Dataset to be displayed then the second and then the last one Finally choose the entity The first two datasets define the spacing between the requested datasets and the third one defines the last dataset to be displayed The last two selections of datasets can be omitted Then all datasets which use the same name starting from the selected one will be used The command
60. ces or a single surface In case of a single surface with the parameter auto all surfaces will be flipped in a way that they share a common direction in or out of a volume See also qfip 8 27 fipc gt flpc This command changes the colour of the scale Initialy the default is red for high values and blue for low values The command will invert the current state 8 28 font font d 1 lt value 1 6 gt 37 This command changes the drawing font Six fonts of different heigh are avail able font d 6 selects the highest font for the drawing area were font 11 selects the smallest one for the legend The default is selected if no number is specified font d or font 1 8 29 frame frame This command fits the model into the drawing area 8 30 gbod gbod lt name char lt 9 gt NORM lt surf gt lt surf gt gt 5 7 surfaces This keyword is used to define or redefine a body in the most basic way see also qbod Each body must have five to seven surfaces to be mesh able However the number of recommended surfaces is six The first two surfaces should be the top and the bottom surfaces For example gbod B001 NORM S001 5002 S005 S004 5003 S006 will create a body B001 The keyword NORM is a necessary placeholder for future functionallity but has no actual meaning Next follow the surfaces with a sign or in front that
61. cond order to first order But in this case each second order element will be splitted in 8 first order elements B 9 How to connect independent meshes Sometimes it is advisable to glue independent meshes together instead of trying to create one big seamless mesh Or you just want to create a contact formulation were you know that no separation will happen In this situations you might use equations also called mpc s which connect one dependent node with one or more independent nodes The independent side should be coarser than the dependent side to avoid gaps in the connection See the command fsend in combination with the option areampc how to create such connections Sliding and fixed contact as well as press fits can be modelled The necessary sets should be defined based on geometry not on the mesh In this case the definition will be stored with the save or exit command and can be used after the next program call Background The movement of each location in an element or on its surface is fully described by the movement of its associated nodes and its shape function In an iterative algorithm element coordinates are varied until the real world coordinates of a dependent node are matched Based on the element coordinates the shape function gives the participation weight factors of the independent nodes the coefficients in the mpc s This approach delivers valid results as long as the dependent node is not located outsid
62. connecting lines and surfaces are created and at last the bodies The divisions of the new lines between setl and set2 is specified with the parameter div or the default is used Existing sets are extended by the copied entities if the last parameter a append is provided Several transformations are available For example scal creates a scaled copy the scaling factors fx fy fz could be chosen independently swep partl part2 scal 2 swep partl part2 scal 1 1 2 PO tra will create a copy and will move it away by the vector dx dy dz and the optional parameter a will assign the new entities to sets were the mother of each entity is included swep setl set2 tra 10 20 30 rot will create a copy and will move it around the axis defined by the points pl and p2 by alfa degrees the connecting lines will be arcs swep set 1 set2 rot p0 px 20 rad will create a copy and will move it radially to the x y or z axis or will create a spherical section if just a single point is defined swep cylinderl cylinder2 rad x 20 swep spherel sphere2 rad pPO 10 mir will create a mirrored copy The mirror plane is placed normal to the 91 direction running from P1 to P2 and placed at P2 swep sectionl section2 mir pl p2 8 111 sys sys lt shell command parameters gt This keyword is used to issue any shell command unix or dos shell For example to move files created by the send command to certain file names or and locations an
63. d in fixed format Node element definitions and results might be in ascii or binary coding The ascii format is able to store element and node numbers up to 99999 in the short form or up to 9999999999 in the long form An example for the short form is shown below 1Ctest 1221C defines a new calc of name test 1UDATE 26 01 2000 1U stores user job informations 2C 1222C starts a block of node coordinates 1 1 0 00000E 00 0 00000E 00 0 00000E 00 1 node A 2 0 10000E 01 0 00000E 00 0 00000E 00 2 node 3 lend of the current block 3C 17230 starts a block of elem definitions 1 1 4 0 Olfirst elem type of that elem is 4 he20 2 1 2 3 4 13 14 15 16 5 6 7 Bo iss 2 12 17 18 19 20 twenty nodes defining that element D y ED ques 3 lend of the current block 1PHID 10 Idefines a parameter HID value 10 100CL101 100C starts a user defined result block 4 DISP 3 1 lAttribute Header Record Dataset 5 D1 1 2 1 O Component Def Record Entity 5 D2 1 2 2 0 5 D3 1 2 3 0 VV VV VV MV MM V V MV V V VV V V V I 1 0 00000E 00 1 00000E 00 1 00000E 00 Nodal Values 104 gt 1 2 1 00000E 00 0 00000E 00 0 00000E 00 gt 3 lend of the current block gt 9999 lend of data The binary format applys only for data lines and the end of the current block line 3 which is omitted All other lines are the same 10 1 Model Header Record Purpose Defines the name of the mod
64. d point and the parameter END will add additional points to a sequence 8 98 seql seql lt set gt lt nr gt Makes splines from all sorts of lines The nr of new created inner points is defined by the parameter nr Also existing splines will be redefined 85 8 99 seta seta lt set gt 207 en le 9p 217 1267 128 b REPS Pvt se S sh 1d lt div gt lt name chars gt n e lt name gt lt name gt lt steps gt This keyword is used to create or redefine a set see alsojgadd All entities like points or bodies and so on must be stored at least in one set to be reachable The set all is created automatically at startup and will be open see Seto all the time unless explicitly closed see Setc To add points to the set dummy type seta dummy p pl p2 This will add the points pl and p2 to the set dummy The following entities are known Nodes n Elements e Faces f Points p Lines 1 Surfaces s Bodies b Nurb Lines L Nurb Surfaces S Values v names of other sets se or shapes sh If the entity of the specified type does not exits a set of that name is assumed and if existing then all it s entities of the specified type are appended seta setl n set2 will append only the nodes in set2 to setl If the type se is used then the full content of set2 is appended Wildcards can be used to search for set names of a certain expression The program will automatically
65. d to start the analysis On certain platforms cgx will not wait for the completition of the command if the amp key was provided as the last argument of the command Otherwhise cgx waits until the command was completed 8 112 test test vw 1 1 n gt e p 1 1 c 8 1 b S L se sh lt entity gt Tests the existence of a given entity Returns TRUE or FALSE and writes it on the stack 8 113 text text lt string gt This keyword is used to define a comment This commend will show up in the menu area of the main window below the file name 8 114 tra tra P P Pwe PdP I Pr lt relative distance gt This keyword is used to move the model in the window For example trau l will move the model 0 1 times the model dimensions upwards The meaning of the other letters is forward f down d right r left 1 8 115 trfm trim rec cyl x y 172 S lt first Dataset Nr gt lt last Dataset Nr gt 92 Changes dataset entities from one coordinate system to another The option cyl transforms the global results to cylindrical and rec from cylindrical to global cartesian In both cases the axis of the cylindrical system must be pro vided Optionally the first and the last dataset of a range of a unique type can be specified The current dataset is selected if no dataset is specified The transformation into the cylindrical system takes place in a way that tensors and vectors are transfo
66. dded to the set be fore the set is written if the parameter comp follows For example geometry like bodies send set fbd comp 76 will be extended by surfaces lines and points or send set frd comp will include all nodes used by the elements to the set elem and will then write the file The parameter quadlin forces the conversion from second order elements into single order elements were each element is subdivided into 8 single order ele ments This takes place during writing and will not change the current state of the mesh in cgx Up to now it is only available for the abq format send all abq quadlin If the parameter tm follows the definition of a solver format send set abq tmf then the mesh plus all temperatures with the necessary solver context for tmf calculations is written If the parameter sta follows send set abq sta 2900 then the mesh plus all temperatures with the necessary solver context for a static calculations is written Here the value for speed is a reference value for eventual numerical values in the TEXT Block of the Datasets in the result file frd file see Nodal Results Block A scaling factor will be determined based on the reference value here a ref speed of 2900 and the individual TEXT Block values factor speed 2 refspeed 2 This factor is used in AMPLITUDE commands which will also be written and will be used to scale static loads which are referenced in the STEP data I
67. determine the type of the entities if not specified but then the names must be unique More than one name can be specified A minus sign between two numbers of nodes or elements specifies a range of entities with steps of steps seta setl n 1001 1100 12 If the sign is specified instead of a setname then the program generates automatically sets with system defined setnames and stores entities in it This can be used is to separate independent meshes and line loops The single inde pendent meshes are then referenced by new setnames for example seta all will determine all separate disjunct meshes in set all and store them in sets called CF lt nr gt Figurel8 99 illustrates this function Other entities have to be specified as an additional argument after the sign 86 seta 1 all will determine all separate disjunct line loops Other entieties are not im plemented so far Values valu can be used as arguments It is necessary to mask a value with a leading V In cases were the value should not be replaces by its content seta set v vall val2 will not substitute the vall and val2 but add the values itself to the set Lines with a division higher than specified with seta setname 1d50 can be stored in a set Here all lines with a division above 50 are stored in setname Figure 9 All disjunct meshes in the jet engine are successive coloured 87 8 100 setc setc lt set gt
68. ds for the mapping process 2D to 2D is listed below Define a set with the slave nodes or faces For example with qadd slave comp slave do to extend the set by the referenced faces or nodes Then open a set and read the master model seto new read result frd add setc If necessary move neu in space to match the position of the slave see move Add the faces to the new set only nodes and elements are already stored in new comp new do plot f new Define the master set qadd master catch the right faces comp master do Map the values of dataset 1 with map slave master surf dsi Check the mapping with ds lel plot fv slave Write the file with the mapped values send slave abq pres dsi el Check the file slave_dsie1 dlo with an editor and use it in an inp file If the master values are not available in either inp or frd format but in any other format which can be read by cgx isaac openFoam it is then possible to write them in frd format with the send command 118 B 12 How to run cgx in batch mode Most commands can be executed in batch mode Actually if you read a file with geometry fbd file then you run cgx in batch mode already You just have to add exit or quit at the end of the file and it will be obvious You might use the read command in such a command file to reference other command files or to read files with mesh entities for whatever purposes In
69. e Etack 8 41 line line lt name char lt 9 gt lt p1 gt lt p2 gt lt cplseq gt lt div gt lt bias gt This keyword is used to define or redefine a line A line depends on points A line can only be defined if the necessary points are already defined There are three different types of lines available The straight line line 11 pl p2 4 is defined by its name 11 the name could have up to 8 characters by the points pl and p2 and optionally by the division The arc line pl p2 cp 4 needs a center point cp The radius changes linear from pl to p2 if the center point cp is excentric The name is chosen automatically triggered by the char acter The spline line 11 pl p2 seq 4 43 4 El needs a so called sequential set seg use the command Fegal or aseq to define such a set This set seq stores the spline points between the end points in the right order The spline function is described in 15 Usually a line is defined interactively with alin 8 42 Inor lnor lt name gt lt p1 gt lt p2 gt lt p3 gt lt length gt A new line normal to a plane defined by three points and a length is created It starts at the last point 8 43 mata mata lt Material Nr gt lt set gt This keyword is used to assign a material number to a set of elements Currently this feature is only useful if the NASTRAN format is used The material number is a numeric attribute assigned to each elemen
70. e written The sliding condition parameter slide were nodes are bound to the element faces to which they belong can be used for structural calculations or for fluid calculations In the later case it forces the flow to follow the surface of the elements This is necessary in inviscid calculations to prevent the fluid to penetrate the walls Results from a previous calculation can be written in the solver input format for further calculations parameter ds tmf etc For example temperatures for thermal mechanic fatigue analysis or velocities pressures and thermals for a restart of a cfd calculation B 11 How to map loads Values at nodes can be mapped interpolated from a second mesh with the map command This second mesh data source or master is usually avail able as a result file from a previous calculation The read command with the add parameter can be used to include this file in the current model Solver input files inp and result files frd can be used as a data source The command 117 will add an offset to the nodes and elements so that existing nodes and elements will not be overwritten The original mesh is the target or slave of the map ping process The mapping process will add the values to the dataset s The mapped values can be visualized by selecting the mapped dataset and entity and then by plotting the slave elements or faces with plot fv slave An example of the necessary comman
71. e appendix and some references used while writing the code 2 Concept This program uses the openGL library for visualization and the glut library 2 for window management and event handling This results in very high speed if a hardware accelerated openGL library is available and still high speed for software rendering MesaGL 1 The cgx has pre and post processor capabilities It is able to generate and display beam shell and brick elements in its linear and quadratic form fig HJ It can display but not create pentahedra and tetrahedra elements But tets can be generated from within cgx if NETGEN 4 is installed and the program ng vol part of NETGEN is accessible see also How to deal with cad geometry The built in mesher creates a structured mesh based on a description of the geometry For example it uses lines for beam elements surfaces for shell elements and volumes bodies for brick elements The program distinguishes between the mesh and the underlying geometry Elements are made from faces and faces are made from nodes If you move a node the corresponding face s and element s will follow The geometry behaves according to the mesh Lines are made from points surfaces are made from lines and bodies are made of surfaces Surfaces might have 3 to 5 edges and bodies might have 5 to 7 surfaces As a result if you modify the position of a point all related geometry will follow In other words if the location of geomet
72. e of the independent element Therefore the location of the dependent node has to be modified in such cases The dependent will be moved to the surface of the independent element again decribed by its shape function and the procedure is repeated B 10 How to define loads and constraints Loads and constraints are not saved in any database They are just created and written to a file with the send command But the sets which are needed for the definition are stored together with the geometry if you type the Bavej command Of course the sets must have been defined based on geometry and not based on mesh entities like nodes because mesh entities are not saved with the save command You must know that geometry sets know also their mesh entities after a Imesh command If you store your commands to write the boundary conditions in a file you can easily repeat them by just reading this 116 Element Face Dependent Node Four Independent Nodes At The Element Corners Figure 11 Dependent node on element face to create a connection command file see read Several types of loads can be applied That is forces pressure thermal heat coefficients and radiation etc For unsupported loads write a set of element node or face labels parameter names and apply the desired load to that set in the solver input file Other boundary conditions like single point constraints for structural thermal and fluid calculations can b
73. e the mesh illuminated chose in the main menu Viewing and continue to the entry Show Elements With Light To see it filled choose in the main menu Viewing and continue to the entity TEM Most of the time it is sufficient to see the surface elements only For this purpose choose in the main menu Viewing and continue to the entry Toggle Surfaces Volumes If you start cgx in the post processor mode as you did in the first example cgx result frd the surface mode is automatically set To see the interior of the structure choose in the main menu Viewing and continue to the entity Toggle Culling Back Front To save the mesh in the format used by the solver type send all abq To store the mesh in the result format type send all frd To create a new model start the cgx by typing cgx b file where file will be the name of the new model if you later exit the program with the command exit The way to create a model from scratch is roughly as follows create e points with apnt or pat e lines with ala e surfaces with qsur e Bodies with Jabod If possible create higher geometry by sweeping or copying geometry with EWepj or copy You might move or scale your model with the command Movel The commands require sets to work with Sets reference entities like bodies or nodes They are usefull because you can deal with a bunch of entities at once See the section How to define a set of entitie
74. e the mouse pointer over the line and press the s key Press q to quit the command It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 87 qsur qsur lt name gt optional RETURN we a bh 21 3 Yr 12 29 rg ig y Ig 9 aq This keyword is used to create or change a surface see also gsur The user might specify a name in the command line or by picking an existing surface with the s key Tf the surface is supposed to be meshed with unstructured triangles element type tr3u or tr6u it is sufficient to select all lines with the I key lowercase L This can be done either in the all mode press a or individual mode press i which is default Then generate the surface by pressing the g key Such surfaces allow holes in it If the surface is not plane it will need to reference 67 an embedded NURBS surface see Inurs or a shape see shpe Despite it is quite an exception it should be noted here that this kind of mesh might be extruded in the 3rd dimension by using the Bwep command to generate penta elements for the cfd solver Only external meshers might use the triangles as a geometry description see How to deal with cad geometry For the definition of a regular meshable surface
75. e volume mode v 7 2 10 Toggle Move Z Zoom Instead of zooming in with the help of the middle mouse button it is also possible to move a clipping plane through the structure to get a view of the inside The clipping plane is parallel to the screen and will be moved in the direction to and from the user by pressing the middle mouse button and moving the pointer up and down while inside the drawing area Usually it needs some mouse movements until the clipping plane has reached the structure Depending on hardware this functionality could be slow After zooming in consider using the Jplot and plus commands to customize your view 7 2 11 Toggle Background Color With this option it is possible to switch between a black and a white back ground 7 2 12 Toggle Vector Plot It is possible to add small needles to the plot which point with their heads in the direction of the vectors Only entities which are marked in the database as vectors will be affected See Nodal Results Block for information on how entities are marked as vectors Internally calculated vector results like the worst principal stress are marked automatically If one component or the value of a vector is selected then the option takes immediate effect This option can be used in combination with Animate Toggle Dataset Sequence See also the keyboard command lds how to select datasets and entities with the keyboard In this case entities which
76. ed It makes sense for volumes were a vector pointing in or out can be determined for all surfaces In this case all related surfaces and embedded elements are oriented in the same way as the selected one It works only in situations were only two surfaces share a common edge This is the case for volumes without inner surfaces or a 2D model In case of a volume all surface normals will point either inwards or outwards depending on the orientation of the selected surface Press q to quit the command It is also possible to measure distances between two pixels on the screen Just press the key vv on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 77 qint qint RETURN w 1212 tg This keyword is used to extend two lines to their intersecting location The end points of the two lines are moved to the intersecting location but they are not merged Use the merg command to merge duplicate points Only straight lines are supported Select lines with the 1 key Press q to quit the command 8 78 qlin qlin lt name gt optional RETURN w b c e g 1 m p gt grs ex This keyword is used to create a sequence of lines or just one see also lline The user might specify a name in the command line or by picking an existing line with the key 1 In this case an existing line might be redefined without destroying related geo
77. ed lines must exist More precisely only single lines or existing combined lines can be detected The user must define the missing surface if just a chain of lines and no lemb is defined between two corner points of the top and bottom sur faces before he can successfully use the body command It is a more convenient way to define a body than the command gbod but exactly 2 or all surfaces must be specified otherwise the body will not be created The most convenient way to define bodies is to use the command qbod For example body b1 s1 s2 26 will look for the missing surfaces and if necessary create them if all lines between the corner points of sl and s2 are defined the result is the creation of body bl Or for example body s1 s2 s3 s4 s5 will create a body and a new name for it The new name is triggered by the sign Here the body is based on 5 surfaces If the surfaces are not connected the body is not mesh able In case a body should only be meshable with tets it can be composed of more than 7 surfaces The definition can be provided by a set of surfaces body surfset will create a body based on the surfaces referenced by surfset 8 6 break break This keyword is used to end the interpretation of a command file The program returns to the interactive mode 8 7 call call lt parameters gt This keyword is used to allow the user to control his own functionality in the file userFunction c
78. ed on this line which in turn are used to redefine the selected line as a sequence line which uses the new spline points Sequences can be shown with the prnt sq command Node sequences can be written to a file and read again with the commands send seqname fbd and read seqname fbd 8 85 qshp qshp RETURN 2w1 n gt p gt 1 g8 1 h 1 s S c 1g This keyword is used to define a shape A shape needs 3 points or nodes points at their location will be generated and used for its definition The points are selected with the p nodes n key After the points are selected it will be generated by a g keystroke The shape can be assigned to a surface by selecting first either a shape h or a nurbs S and then the surface with the s key The shape will then define the interiour of this surface The selected shape or nurbs stay selected until cleared with c It can also be used for projection purposes see proj or splitting see Em 8 86 qspl qspl RETURN w s q This keyword is used to split one or more lines at a certain position can also be done with alin A point is created at the splitting position the original line is deleted and two new lines will appear instead All lines running through the selected location will be split at once and the newly created splitting points will be merged to one if they are closer to each other than defined by fgtoll To split a line mov
79. el Format 1X 1 C A6 Values KEY CODE NAME 10 2 User Header Record Purpose Stores additional user informations regarding the job ie user id creation date model informations Format 1X 1 U A66 Values KEY CODE STRING 10 3 Nodal Point Coordinate Block Purpose Defines the nodal coordinates 1 Record Format 1X 2 C 18X 112 37X 11 Values KEY CODE NUMNOD FORMAT Where KEY 2 CODE C NUMNOD Number of nodes in this block FORMAT Format indicator 0 short format 1 long format 2 binary format coordinates float 3 binary format coordinates double Following records ascci FORMAT 0 1 Short Format 1X 1 15 3E12 5 Long Format 1X 1 110 3E12 5 Values KEY NODE X Y Z Where KEY 1 NODE node number Kes coordinates 105 Following records binary FORMAT 2 Format int NCOMPS float int and float are ansi c data types Values NODE X Y Z Where NODE X node number coordinates Following records binary FORMAT 3 Format int NCOMPS double int and double are ansi c data types Values NODE X Y Z Where NODE Koz node number coordinates Last Record only FORMAT 0 amp 1 ascii omitted for FORMAT 2 amp 3 Format 1X 3 Values KEY 10 4 Element Definition Block Purpose Defines the topology of the elements 1 Record Format 1X 3 C 18X 112 37X 11 Values KEY CODE NUMELEM FORMAT Whe
80. es gt This keyword is used to define a nurbs shape Surfaces might use nurbs to define the interior geometry There are two possible ways of definition The first is using predefined point names and the second is by specifying the coordinates explicitly When the coordinates are defined then the parameter COMPACT 50 must be used as shown above but when point names are used then COM PACT must be omitted CAD interfaces might use this functionality REMARK The knot vector has to have a multiplicity of degree 1 There is also a small nurbs building capability in cgx It is possible to use existing surfaces with 4 edges which do not already reference a given nurbs The new nurbs will follow the Coons algorithm but can be modified by moving the control points The surfaces must be stored in a set For example nurs surfaceSet will define nurbs for all surfaces stored in the set surfaceSet This nurbs can be used to define the interiour of other surfaces This is necessary if tr3u ele ments unstructured triangles should be used and if the surface is not related to a given shape Note SU offers another option to create nurbs related sur faces by associating existing surfaces to an overlapping existing NURBS The interiour of the surface is then defined by the NURBS 8 59 ori ori lt set gt This keyword is used to trigger the orientation of the entities This is done automatically and it should
81. es on which the values are known were included from an external result file with the read command with the parameter add For further details see How to map loads in the appendix 8 45 mats mats This keyword is used to get an overview of all material numbers defined New sets will be created for all material numbers of elements This feature is usually used if a mesh with mateial numbers was read For example if a native netgen file 4 was read then this command will generate sets which can be used to asign boundary conditions 8 46 max max lt value gt DI 4 2e 1 au This keyword is used to define the upper value of the scale in the menu area see figure The number representation can be changed between int float exp For example max 1100 i will set the upper value to 1100 and the representation to integer A third parameter 1 lock or u unlock can be provided which locks the scale to certain max or min values The selection of a different dataset will not change the scale 8 47 merg merg n e p 1 c s lt set gt lt gtol gt nolock This keyword is used to merge close points and nodes or equally defined enti ties The following entities are known Nodes n Elements e Points p Lines 1 Combined Lines c Surfaces s For example to merge points included in the set point set type merg p point set 45 Only entities included in the se
82. es the 2D meshing regions as separate sets which can be used to apply bound ary conditions cgx ng file vol See the Program Parameters section on how to do that e The CAD model is available in VDA format In this case the user can mesh the surfaces with structured elty setname tr3 etc or unstruc tured triangles 14 elty setname tr3u etc and use them to be replaced by tets created by an external auto mesher called from within cgx The cgx uses the tet mesher from NETGEN 4 for this task Therefore cgx can only generate a tet mesh if NETGEN is installed and the program ng_vol part of NETGEN is accessible For the full functionallity replace the original ng vol in the NETGEN package with the modified program ng vol from the cgx distribution and build it again The following section describes the process to generate a tet mesh with cgx based on a cad model in VDA file format Figure 12 CAD geometry meshed with tet elements e Use the interface program to convert the cad format to fbd format You find example data vda in CalculiX cgx 2 7 examples cad 120 ie vda2fbd halter vda gt halter fbd In some cases the header of the vda file causes trouble and must be mod ified to make the program run Start cgx with that file cgx a halter fbd Usually some warnings appear on the screen but the program will fix that automatically triggered by the parameter a After all geometry has been read the progra
83. et2 nor 0 7 If a point does not hit any surface from the target set then it will not be moved A tolerance can be specified for the projections No projection takes place if the target surface is farther away as this value defines proj setl set2 nor 0 7 10 8 65 qadd gadd lt set seq gt t lt value gt RETURN gt LU Pa PE Prae eE PpP gt a PE PE PE PE P leag gt e at eu gt This keyword is used to add entities to a set See also seta But a set will not keep the sequence in which the entities were selected UseBeqajor the command asegj if the order of the selected entities has to be kept After an entitiy was selected you get certain informations about the entity If the node which belongs to the maximum or minimum value in a certain area has to be stored in a set you might use the same key strokes as described for the command aeng To catch more than one entity with one stroke type a all at first Then create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the selection rectangle To catch only the entity which is closest to the user type i before Then move the mouse pointer over the entity s and press one of the following keys for Nodes n Elements e Faces f Points p Lines l Surfaces s Bodies b Nurb Surfaces S node attached texts t and for Nurb Lines L 56 If faces f of a certain area have to be selected the user might specify a toler
84. et4 elements which use and keep the shape of the provided tri3 elements B 21 Remarks Concerning OpenFOAM The mesh can be written in OpenFOAM polyMesh format with the command Bendl If you work in the polyMesh directory of the OpenFOAM case then all mesh related files will be already in place So far the physical type is not written in the boundary file as it is not mandatory The results of an OpenFOAM calculation can be viewed by specifying the parameter foam and the case the relative or absolute path including the directory name of the case See also Program Parameters Imesh 125 B 22 Remarks Concerning Samcef From Paul CARRICO 2006 04 17 BASIC TUTORIAL FOR HOWTO USE THE SAMCEF EXPORT FORMAT ANNAN In K Definition of the points pnt p0 0 0 0 pnt pi 100 pnt p2010 pnt p3 210 plot pa all Definition of the lines plus 1 all qlin link the points pO pl p3 p2 p2 p0 Creation of the first surface qsur Creation of 10 between pO amp p1 and 11 between pi amp p3 sets gadd lo gadd 11 Creation of the 2 other surfaces swep 10 11 tra 0 3 0 swep 11 lib tra 3 0 0 Creation of the SYMETRY set plot s all qadd SYMETRY use both a and rr keys to select all the surfaces Creation of the volumes swep SYMETRY si swep tra 0 0 1 all the volume will be automatically created Looking for common points lines and surfaces the order merg p all merg 1 all
85. eyword is used to measure distances between two nodes or points Move the mouse pointer over one entity and press the following key for a node n for a point p or for a center point c If the key c was pressed then the distance between the two nodes or points are given in cylindrical coordinates Here lcir is the length of the arc da the angle dr is r2 r1 If no center point was chosen then the distance and its xyz components are given in the Cartesian system But in 59 addition the cylindrical distances around the origin and around the xyz axis are 99 NU also given Press q to quit the command e dist global distance e dx dy dz distance in the three Cartesian directions da global angle dax day daz angle around x y z dr global radius difference e drx dry drz radius difference around x y z It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 73 qdiv gdiv RETURN Iw a l i 1 79 10 99 1 q This keyword is used to change the division of a line by pressing a number between 1 and 9 when the mouse pointer is at the position of a line see also div To define numbers between 10 and 99 press the space bar when the mouse pointer is at the position of a line and then the number To select more than one line with one stroke type a before a
86. eyword is used to scale the graphic presentation of values but no values of entities The scaling of the value itself can be done with the command ds For example the command scal s 0 5 will scale the range of values presented in the color plot by a factor of 0 5 The command scal v 5 will scale the length of vectors by a factor of 5 The command 73 scal d 2 will scale the deformed shape by a factor of 2 Without a value it restores the default value See also Toggle Vector Plot Toggle Add Displacementy Datasets 8 96 send send lt set gt send lt set gt fbd c e gt stl gt abq ans ast dar duns frd gmp gt 1st nas ng pat tcg lt RETURN gt bin ds lt nr gt e lt nr gt lt nr gt P tmf sta lt refspeed gt crp lt timefact gt lt refspeed gt lt writefreq gt comp mpc lt rotation gt v lt node gt lt v1 gt lt v2 gt lt v3 gt n lt node gt names spc lt dofs 1 6 tlp gt lt value gt gt lt dofs 1 6 tlp gt ds lt nr gt e lt nr gt gt nor lt faci gt lt fac2 gt lt fac3 gt spct lt dofs 1 6 tlp gt lt value gt gt lt dofs 1 6 1t p gt ds lt nr gt e lt nr gt P slide c s rx gt ry l rz tz ty tz pre
87. f the parameter crp follows send set abq crp 1 2900 1 then the mesh plus all temperatures with the necessary solver context for creep calculations is written Here the value for timefact scales the time values of the Datasets speed is a reference value for scaling static loads see option sta and writefreq limits the output to the result file If the parameter names follows the solver type then just the element or node numbers are written This is useful for defining sets for the solvers The com mand send steel abq names TT will write the node and element numbers included in steel This file could be used to define a set for material assignments boundary conditions or loads If all sets should be written at once use send all abq nam If the solver format duns is used then related numbers of surfaces 3D or lines 2D will be written This information is necessary to apply boundary conditions to duns The numbers are used in the connectivity file duns conn which will be created by the mesh command The user has to refer this numbers in the duns script2 file when assigning boundary conditions It is also possible to define and send some constraints to the file system To constrain degrees of freedom dof of selected nodes use the parameter spc in combination with the numbers of the constrained dofs and optionally a forced deflection For example send hinge abq spc 123 will co
88. face without holes can be splitted by a new line e First create the desired splitting line for example with falinj e Convert the line to a spline with qseq optional e Project the line points to the surface with qpnt or proj optional e Type qmsh and identify the line with key P y y a Identify the surface with key s e Press key b to break split the surface in two 66 99 This keyword is also used to manually define elements The key n selects nodes and the key s or v generate a surface or volume element based on the node selection Mostly used to repair stl files manually 8 80 qnor qnor RETURN w p A new line normal to a plane defined by three points and a length is created It starts at the last point 8 81 qpnt gt qpnt lt name gt optional RETURN gt w Pp Pg PePe PEPE Pre This keyword is used to create or move points see also pnt The user might specify a name in the command line if a certain name should be used To create a point move the mouse pointer to the desired location and press the key g generate or over an existing node and press n uses then the node coordinates After a point was selected with p it can be moved in different ways Either in the screen plane for this go to the desired position and press m move see gal or the section Orientation how to rotate the model into a certain position
89. gt lt line lcmb gt lt line 1cmb gt This keyword is used to define or redefine a surface It is a more convenient way to define a surface than the command EST Either individual lines or a set can be specified To be meshable with structured elements it requires 3 4 or 5 edges lines or lembs To be meshable with tr3u or tr6u it requires either to be plane or to reference a shape or nurbs which can not be provided with this 90 command see fgsur The most convenient way to define surfaces is to use the command RSU If the name of the surface should be automatically generated then just type instead of a name 8 110 swep swep lt set gt lt set gt scal lt fx gt lt fy gt lt fz gt lt PO gt lt div gt a tra lt dx gt lt dy gt lt dz gt lt div gt a rot lt p1 gt lt p2 gt lt alfa gt lt div gt a JI Prot x y z p lt pnt gt lt dr gt lt div gt a JI rad lt p1 gt lt p2 gt lt alfa gt lt div gt a JI Prad x y z p lt pnt gt lt dr gt lt div gt a JI mir lt P1 gt lt P2 gt lt div gt a This keyword is used to sweep entities into the next higher dimension Sweeping a point will create a line sweeping a line will create a surface and sweeping a surface will create a body At first a copy of the first set see Seta and Copy will be created The copy of the master set is included in the second set Then the
90. gt 1 x 1 y 1 2 1 1 IP x ly ez gt p lt pnt gt lt dr gt lt dri gt lt ax1 gt lt dr2 gt lt ax2 gt equ lt trgt set gt lt tol gt mir lt P1 gt lt P2 gt This keyword is used to move nodes or points which are stored in a set Re lated entities will be moved as well For example to move a line it is necessary to include their points in a set seegomp Several transformations are available For example scal will scale the entities of the set the scaling factors fx fy fz can be chosen independently and a reference point can be used move part scal 2 move part scal 1 1 2 move part scal 2 PO move part scal 1 1 2 PO tra will move it away by the vector dx dy dz move all tra 10 20 30 rot will move it around the axis defined by the points pl and p2 or the axis x y z by alfa degrees move all rot p0 px 20 rad will move it radially to the x y or z axis or two points as above or to a single point move cylinder rad x 20 move sphere rad pPO 10 48 The axis for the rad or rot commands can also be specified by one point and on main axis x y z as shown in the following example move all rot P y 20 The change in length or angle might be interpolated for the rad or rot cases move set rad x 1 120 2 140 The number 1 specifies the radial change around x at 120 length units along the x axis and 2 is the change at 140 length units With parameter equ points will
91. he normal direction of the independent faces works also with option f Additional nodes are generated and can be used to request the reaction forces on the dependent nodes They are stored in a set named N dep set name ind set name Two files are produced The one with the equations has to be included in the model definition section and the boundary file in the step section Special cases are the cfd solvers Duns Isaac and OpenFoam The boundary patches are an integral part of the mesh So it is necessary to specify all bound ary patches when writing the mesh All free surfaces of the mesh must be specified This is an example for OpenFoam send all foam cyclic cycl cyc2 cx patch in patch out wall wall will write the so called polyMesh description to the file system After the send command the set all with the mesh is specified then the format foam then cyclic boundary conditions cyclic cyc1 cyc2 cx between set cycl and cyc2 of the axi symmetric case around x cx then boundary conditions of type patch for set in patch in then boundary conditions of type patch for set out patch out and boundary conditions of type wall for set wall wall wall The symmetric boundary conditions base type cyclic can be axi symmetric c around x y z or rectangular r in direction of x y z Only for OpenFoam and in the rectangular Cartesian case also a vector pointing in the direction of the symmetry can be specified ie 1 1
92. he only exception are the element types tr3u and tr6u This types are generated in an unstructured way using the mesher from 14 In this case holes are permitted and a surface coming from a cad system should be meshable These surfaces MUST reference a given NURBS surface or a shape to be meshable if they are plane a shape will be genrated atomatically This is usually the case if the data is derived from a cad system with the interface program vda2fbd If NETGEN 4 is installed and if the elements form a closed volume they can be used to generate a tet mesh mesh set with trias tet or mesh set with trias tet lt element target size gt 46 This is a second method to generate tets The other one is to use Teltyl to asign tet elements to bodies 8 49 mids mids lt set gt lin gen rem This keyword is used correct the midside node position of higher order elements stored in a set It is performed automatically if a new mesh is created or if nodes are projected to target surfaces The correction will use a circular track defined by the corner nodes and the midside node With the parameter lin the corrected position is halfway between the corner nodes for all inner nodes Nodes on the surface are not affected by the lin option The gen option will generate midside nodes for linear elements like he8 or te4 The rem option will remove midside nodes from the element formulation but the nodes are not
93. hen the 1 would be the number of the first equation No number is required for abq abaqus and calculix see also How to connect independent meshes There is also the slide option in combination with the areampc option For example if the mesh of a turbine blade and a disk should be connected with each other in a simpler but realistic way then a sliding condition between this parts can be established The command send dep indep abq areampc slide will connect the nodes in the set dep to element faces described by nodes in cluded in the set indep but only in the direction perpendicular to a plane defined by nodes of the dep set Therefore all dep nodes and all indep nodes must lie 83 in the same plane and will slide in the same plane Another case is considered with the presfit option in combination with the areampc option For example if a cylindrical press fit should be simulated then a forced displacement between the two intersecting surfaces is necessary This forces the dependent nodes to move to the independent face Two modes are available send dep indep abq areampc presfit f simulates sticking friction and with the option s send dep indep abq areampc presfit s works for sliding conditions The user might request a certain value for the press fit if the overlapping of the mesh do not represent the necessary distance send dep indep abq areampc presfit s0 06 will move the dep nodes 0 06 in t
94. hermal and thermomechanical calculations The Thermal solver per forms linear and non linear calculations for pure thermal but also for thermomechanical simulations e Input interface EFICAS is the input interface for coding the input file but it s not a pre processing as you can find in many commercial code 123 e Tools Many tools are provided with Code Aster HOMARD for mesh refinement GIBI and GMSH for post processing CGX to ASTER export format i e HOWTO use this export format The export format allows to export meshes and sets from CGX to Code Aster for linear and quadratic 1D 2D and 3D elements For this just type send all aster for exporting a complete mesh The file will have the mail extension Type send set aster nam for exporting the GROUP NO node group and the GROUP MA mesh group which compose the set The file will have nam extension The later item is particularly useful to apply boundary conditions onto DOF pressure force displacement temperature etc Since of these bound ary conditions are coded with EFICAS in the input file I thought it was not necessary to develop another features than the 2 previous ones So far results can not be read B 17 Remarks Concerning dolfyn Some support for dolfyn a free cfd code was provided by Runar Tenfjord In the CalculiX cgx nr examples dolfyn directory is a patch for the dolfyn source code included This patch enables dolfyn to write frd
95. hout destroying the definition of related bodies or other geometries After selecting a surface with the s key the referenced lines might be replaced by selecting new ones The previous selected lines are not longer referenced Replace assignments by the blended function with the b key To quit the command use the q key or use u to undo the last action It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 88 qtxt qtxt RETURN g a 3 p ip 2m an ry f 2q q This keyword is used to move node attached texts showing node number and value see figure 8 to certain positions in the drawing area or to manipulate them otherwhise They are created with the key g while the mouse pointer is over the node to which it should be attached Attention They are only visible after the command plus nt all was used but they are immediatelly visible when the the commands REM or gadd were used instead of qtxt To move a text pick it at the lower left corner with the key p and place it with the key m in the new location To move it back to its node use b Delete 68 them with d The node nr and the value in the text can be switched on and off with the n and v key The f key will toggle the format of the value
96. im tune lt value gt steps lt value gt time lt value gt gt real on off model on off gt elem on off This keyword is used to manipulate the animation of displacements See also Ids and beal The amplification is controlled with tune steps defines the number of frames over one periode time controlls the duration of one periode real switches of the automatic amplification and the real displacements are 24 used instead In addition the displacements of the negative part of the periode is set to zero model switches the static model undeformed edges on or of elem does this for the element edges 8 2 area area lt set gt This keyword is used to calculate the area and the center of gravity of a set of shell elements or surfaces of volume elements If a dataset is active then an averaged value is calculated The command writes to the tack 8 3 asgn asgn Pn le lp ll Pes eb 17S L se gt sh alpha beta nadapt lt value gt gt rbe lt value gt mpc This keyword is used to manipulate the behaviour of successive commands For example to define the first node or element number which will be used for the next mesh generation And it is used to redefine the leading character of new entities The default is D for points p L for lines 1 C for combined
97. indicates the orientation of each surface These signs will be corrected automatically if necessary If the an automatically generated name is desired then type instead of a name 8 31 gonly gonly on off This keyword is used to erase the contents of the menu area Sometimes this is useful for hardcopies 38 8 32 graph graph lt amplitude chars gt amp 1 lt material chars gt mat 1 lt seq gt length lt set gt p nr time descr lt parameter name gt gt lt Dataset name gt lt entity name gt lt parameter name gt gt lt first Dataset Nr gt lt last Dataset Nr gt 1 This keyword is used to create time history plots of nodal values The values of the pre selected nodes stored in the referenced set will be written to a file called graph out A gnuplot command file will be written called graph gnu and executed The resulting file graph ps will be automatically displayed with the postscript viewer The default is ghostview but this can be changed in the cgx h file before the program is build See also How to generate a time history plot for further details A set with the nodes of interest must be created usually with gadd be fore a graph can be drawn If the parameter 1 length is used to plot values along the length of a set of nodes then this set has to be of the sequence type usually created wi
98. ines surfaces But surfaces can only be identified if the command rep lt setname gt was issued before In case of nodes with related values results it will also determine the highes or lowest value in the specified range or all nodes above or below a certain value The coordinates might be taken from one node or point in a given 3rd set or given in cartesian coordinates option rec or cylindrical coordinates option cyl Some coordinates might be omitted to specify an infinite range The J key has to be used in this case The mode is defined by the keys i individual a all h high V low were h and 1 will search the highest or lowest value in range This value will also be written to a file for automatic processing purposes For example enq all newset rec 10 100 0 1 i will search for entities in set all at the given location x 10 y infinite range triggered by the key z 100 with a tolerance of 0 1 and only the closest entity of all kinds of entities are stored in set newset triggered by the i key With the a key all in range would be stored The key h or 1 would trigger a search for the highest or lowest value in the specified range In combination with a certain value enq all newset rec 10 100 0 1 h 1013 all nodes with a value above 1013 would be stored in set newset The command enq domaini newset cyl x 100 10 0 1 a will search for entities in set doma
99. ing purposes The default is off 8 56 node node lt nr gt lt x gt lt y gt lt z gt This keyword is used to define or redefine a node For example node 23580 10 0 1 defines a node with the number 23580 at the position x 10 y 0 z 1 8 57 nurl nurl lt name char9 gt DEFINE COMPACT gt lt pstart gt lt pend gt lt deg gt lt npnt gt lt nknt gt lt div gt l CONTROL lt index gt lt pnt gt lt x y z gt lt weight gt KNOT lt index gt lt value gt P END This keyword is used to define a nurbs line So far this command is only used to read a nurbs line definition Nurbs lines are converted automatically into a spline Nurbs lines can be displayed but not saved There are two possible ways of definition Either by using predefined point names or by specifying the coordinates explicitly When the coordinates are defined the parameter COM PACT must be used as shown above When the point names are used then COMPACT must be omitted CAD interfaces might use this functionality 8 58 nurs nurs lt name char9 gt DEFINE COMPACT gt lt u deg gt lt v deg gt lt u npnt gt lt v npnt gt lt u nknt gt lt v nknt gt CONTROL lt u index gt lt v index gt lt pnt gt lt x y z gt gt lt weight gt P KNOT lt U gt lt V gt lt index gt lt value gt END lt gt lt setname containing surfac
100. inl at radius 100 around x at x 10 with a tolerance of 0 1 All entities in range are stored in set newset The following example illustrates the use of enq to find the highest value close to a location given by a node number activate dataset 3 with entity 7 ds 3 e 7 H read a file vith one node to define a set read posi frd posi search the highest value around a radius of 1 enq all t set posi 1 h H move the file with the search result to a meaningfull name sys mv enq_1c3_e7_1 0ut pos1 out The command writes to the stack 36 8 24 eqal eqal jbir aspr mca lt value gt The keyword without parameters lists the current element quality thresholds 0O off e JBIR The maximum ratio of jacobian determinants xsj between all gaus points Were xsj Vworldcoordinates Velementcoordinates and jbir xsj max xsj min e ASPR Element aspect ratio between all element sides e MCA Maximum corner angle between all element sides With parameter value sets the value of the threshold 8 25 exit exit This command terminates the program and saves the geometry if any to a file named as the input file but with the extension fbd If a file with that name exists already then this file will be saved with the new file extension fbb as a backup see also pave 8 26 flip flip lt setname gt lt surf gt auto This command changes the orientation of a set of surfa
101. instead of a uniform value The front or rear side of the face pos or neg is selected with the or parameter it has to be the last parameter The parameter dflux is used to assign an energy stream to element faces It works similar to the pres option For example send surf abq dflux 0 5 will assign the flux of 0 5 to all element faces in set surf If a dataset with nodal values is available it is also possible to use this values instead of using a uniform value for all faces For example send surf abq dflux lei el will specify the dataset nr 1 and the entity nr 1 to be used instead of a uni form value The front or rear side of the face pos or neg is selected with the or parameter it has to be the last parameter The parameter cflux is used to assign an energy stream to nodes For example send surf abq cflux 0 5 will assign the flux of 0 5 to all nodes in set surf The parameter force is used to assign force values to nodes For example send nodes abq force 1 20 0 will assign the specified forces to all nodes of the set nodes The parameter mpc is either used to create input for the user subroutine umpc which forces all nodes from set to rotate by an average value speci fied with rotation around the vector v or to create a rigid body In case of a rotation the value has to be in degree were 90 degree is orthogonal For example send nodes abq mpc 4 1 0 0 81
102. isplay purposes The Dataset might be a positive number which has to match the nr in the Dataset menu or an V lower case L which is interpreted as the last available Dataset or a negative number Then it is interpreted as the last ds minus the specified number For example ds lel will display the first entity of the first Dataset dslel will display the last Dataset To start the animation of the second to last Dataset assuming displacements right away ds 1 a Also sequences can be defined ds 24 10 e 1 Here the 1st entity of each second Dataset is selected The selection starts at the second and ends at the 10th dataset If more than one entity is defined then a vector plot will be displayed If a 4th entity is defined then this entity will be used for the basic color plot ds 2 4 10 e 12 13 14 15 In case the deformed shape should be shown together with the fringe plot then the e parameter has to be replaced by an a character ds 2 4 10 a 12 REMARK Vector plots can not use the deformed shape Therefore only one entity is supported In addition it is possible to scale the entities of the specified datasets ds 1 s 1 2 will scale all entities of dataset 1 by a factor of 1 2 33 8 18 elem elem lt nr gt lt firstNode gt lt lastNode gt be2 be3 tr3 tr6 gt qu que he8 he20 This keyword is used to define elements based on nodes and its type see
103. ities to define a point To define a point just with coordinates pnt pl 11 1 2 34 or pnt 11 1 2 34 where the name is chosen automatically It is also possible to create points on a line or in the direction from P1 to P2 by defining a spacing ratio and number of points pnt L1 0 25 3 or pnt P1 P2 0 25 3 will create 3 new points at the positions 0 25 0 5 and 0 75 times the length of the line or the distance from P1 to P2 and it is also possible to create points on the positions of existing nodes The command pnt set will create new points on the positions of the nodes included in the specified set Usually when points are defined interactive the command is used 8 63 prnt prnt se sq lt RETURN set chars gt range n e lt set chars gt range 53 PaPe PEP PP PP PE rP IPL HS I Py lt entity gt amp lt RETURN amplitude chars gt mat lt RETURN material chars gt l par lt RETURN parameter gt st usr This keyword is used to print entity definitions The following entities are known Nodes n Elements e Faces f Points p Lines l Surfaces s Bodies b Nurb Lines L Nurb Surfaces S Values v Sets se and Sequences sq To see all known sets type prnt se Or type prnt sq to see all known sequential sets sequences Wildcards can be used to search for setnames of a certain expression In this case all se
104. l complete the set downwards before writing or send set fbd e will write all model edges as small lines running from node to node and a surface description with triangles stl send all stl This triangles are based on elements which were created by meshing surfaces or by automatically triangulated element faces of all types of supported elements If be2 elements are included in the mesh meshed lines a so called edge file for NETGEN will also be created if stl is written NETGEN provides also a stand alone mesher called ng vol The mesher can be found in the netgen sub directory nglib As for the stl format triangles can be written for this mesher with send set ng It has the advantage that the triangles are directly used to define tetras and not as with stl are only used to define the outer shape of the body The following mesh related formats are known Femview and CADfix frd Nastran nas Abaqus abq Ansys ans CodeAster ast Samcef sam Dar win dar patran pat only sets gagemap gmp only sets and Tochnog tcg but only ccx Abaqus is fully supported Good support is also provided for ansys darwin and nastran If no further parameter follows then just the defini tion of the nodes and elements will be written send all abq will write the mesh in the format used by Abaqus and the CalculiX solver If the parameter ds is provided followed by the dataset nr 75 send set abq dsl then the values stored in
105. llowing mesh entities to files see send the necessary key parameters are listed in brackets below e Nodes and Elements e Sets of nodes and elements nam Single point constraints spc Equations areampc e Pressure pres Temperatures ds all datasets with just one entity will be written as temperatures The resulting files have to be combined with the help of an editor and extended by material data and the necessary controll commands So far results can not be read B 16 Remarks Concerning Code Aster From Paul CARRICO 2005 02 12 Brief presentation of CODE ASTER Code Aster is an implicit solver under GPL licence from the French company EDF Electricit de France Code Aster and its documentations is downloadable at the following address http www code aster org NOTA the documentation is in French at the moment but many users are translating it in English and in German The following list is not exhaustive but it briefly presents the capabilities of the solver e Mechanical modeling Code aster allows linear and non linear calculations in static dynamic acoustic etc Many mechanical laws are implemented in the solver such as damage fatigue creep viscosities elastic plastic etc For isotropic and non isotropic materials orthotropic ones for example Because of EDF fields the materials used in Code Aster are essentially Metallic ones and Geo material one but there are probably some others e T
106. m will merge points and lines to close the volume Then all nurbs related surfaces might be trimmed optional This has to be done manually by typing rep all Finally with the command plus si all all surfaces are rendered and the result is displayed in the main window The whole process might need some time If problems occur the user might start the program with option b instead of a cgx b halter fbd Then no automatic post processing is done and the user has to prepare the geometry manually See the commands If a surface points inwards it is not illuminated and appears only dark grey Flip the surfaces in a way that its normal direction points outwards use afipl Press the a key to selet the all mode and then select a correctly oriented surface by pressing key s all other surfaces will immediatelly use the same orientation Define the element type see elty In most cases it is advisable to use triangles first and not directly tetraeders because this has the advantage that the surface mesh can interactively improved before the tets are gen erated Remark If you choose tets from the beginning you must first create a body with body all before you can use the elty command Type the command mesh all to generate the surface mesh see mesh Visualize and check the mesh plot f all plus Ip all plus ld all If some surfaces can not be meshed then modify the divisio
107. mber between 1 and 9 when the mouse pointer is at the position of a line To define numbers between 10 and 99 press the space bar when the mouse pointer is at the position of a line and then the number To select more than one line with one stroke type a before and create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the rectangle To select only one line type i before Press c to change the direction of the bias Press q to quit the command see alse bia 8 68 qbod gt qbod lt name gt optional RETURN w b a i r s g 2 2 2 2 qg Pu 57 This keyword is used to create a body see also and body The user might specify a name in the command line or by picking an existing body with the key b Otherwise the program chooses an unused name It is possible to create the body out of five to seven surfaces which are needed to define a body or just of two opposite surfaces but then these two surfaces must be connected on their corner points by lines To be more precise only single lines or existing combined lines will be detected If a combined line would be necessary but does not exist then the user should define a surface using this lines which will create the necessary combined line Other missed surfaces will be created automatically To catch more than one surface with one stroke type a before and create a rectang
108. me can be reloaded by pressing the middle mouse button twice and the right button once while the mouse is in the menu area 6 2 Keyboard The Keyboard is used for command line input and specifying the type of entities when selecting them with the mouse pointer The command line is preferable in situations where pure mouse operation is not convenient i e to define a certain value or for batch controlled operations Therefore most commands are only available over the command line The stream coming from the keyboard is echoed in the parent xterm but during typing the mouse pointer must stay inside the main window Otherwise the commands will not be recognized by the program The following special keys are used Special Keys ARROW_UP previous command ARROW_DOWN next command PAGE_UP entities of previous set if the last command was plot or plus or the previous Loadcase PAGE_DOWN entities of next set if the last command was plot or plus or the next Loadcase 15 MENU AREA DRAWING AREA COORDINATE SYSTEM DRAWING MODE SURFACE VOLUME FILENAME 4340109 Figure 2 structure of the main window 7 Menu The main menu pops up when pressing the left mouse button inside the menu area figure 3 It should be noted that there are equivalent command line functions for most of the menu functions This can be used for batch controlled post processing Next the entities inside the main menu will be explai
109. me itself differs send dep indep abq cycmpcf rx12 cl The function name for cfd elements is cycmpcf instead of cyempc In addition it is possible to glue independent meshes together For this purpose the dependent nodes are tied to independent elements by equations Choose the finer mesh for the dependent side The equations are based on the shape functions of the element types For example send dep indep nas areampc 123 cl will connect the nodes in the set dep to element faces described by nodes in cluded in the set indep The set dep must contain all nodes which should be glued and the set indep should contain all nodes of the elements surfaces to which the dep nodes should be glued The numbers 123 are the degrees of freedom which will be connected t will create a thermal connection p a pressure connection The c triggers the correction of the position of the de pendent nodes to a position on the surface of the independent elements highly recommended and default u would prevent the correction f forces the dependend node away from the independent face Of course the mesh has to be written after the use of such a command otherwhise the corrected node positions would not be regarded and the equations would lead to increased stiff ness and decreased accuracy The 1 will be the identifier for the equations if the format is nas nastran In case the format would be ans ansys t
110. metries To start a sequence of lines move the mouse pointer over a point and press the key b begin To define just a straight line go to the end point and press g generate This point is also the starting point for the next line So no b is necessary for the next line 63 If the line is a sequential line spline then the points of the sequence must be defined with the key t in the correct order but the last point must be selected with the key g If the line is an arc then the center point must be selected with the c key or one point on the line between the end points must be selected with the m key but the last point must be selected with the key g In case you need to split a line during the process you can do that by pressing the s key When you want to modify the length then pick an existing line with the key 1 and then the endpoint you want to move with the key p The displacement is requested and the point is moved in the direction of the line by the specified displacement The x key will transform the line into a straight line One quite usefull function for cad based surfaces is triggered by the e key exchange This allows to modify an existing surface gsur With that command a line or a line sequence in a lemb can be exchanged by the previous selected 1 key or generated line For example one edge of an already existing surface uses a lemb but
111. mouse pointer into the xterm were the typed commands and the protocol is visible Please the mouse pointer MUST stay in the main window during typing 111 A 2 During Meshing If a body can not be meshed although no obvious error is in its definition then the reason could be that the shape of the body is very far from being brick like Subdivide the body to improve the shapes of the single ones For the moment there are some restrictions for the definition of five or seven sided bodies For example the first two surfaces in the body topology see gbod have to be defined in the same order That means the first line of the first surface has to be connected with the first line in the second surface by one of the remaining surfaces This is always the case if the body is a product of a Bwepj command A 3 Program generates a segmentation fault Write a mail to the author and if possible add the input file 3 B Tips and Hints The following collection will give you background information and procedures to deal with common situations B 1 How to change the format of the movie file Use a shell command like convert to split up movie gif with multiple layers into multiple frames as seperate jpeg files convert movie gif 4d jpg Convert any series of multiple jpeg files into WMV format with mencoder mf jpg mf type jpg fps 25 ovc lavc lavcopts vcodec wmv2 o movie wmv To convert a DivX video file into Flash Movie format
112. n the 3rd dimension to create your geometry see Bwep If features exist in the 3rd dimension then this features must be included in the 2D section Then sweep the 2D section to the location were the feature starts then right to the end and at last to the end of the geometry Then delete the unwanted bodies before and after the feature You might also project the sweped sections to target surfaces if the feature is shaped in the 3rd dimension see proj Sometimes the geometry consists of several parts which are arranged on different positions in space This is called assembly When cgx is started with the parameter stepsplit instead of step it will write the single parts to separate directories using their original coordinate systems The user can prepare meshable geometry or meshes based on them But he must use the filenames which he finds there The final meshed assembly can be build afterwards by calling the fbl file which was written by cgx This fbl file contains the original step commands to position and eventu ally duplicate the single geometries meshes from the subdirectories 119 e For CAD models in STEP or IGES format you should consider to use a tet mesher like NETGEN 4 which generates quite nice tet meshes You can read this meshes with cgx and combine them with cgx geometry and meshes Then create your boundary conditions etc You might read the native netgen format vol instead of abaqus format because this includ
113. nd create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the rectangle To select only one line type i before Press q to quit the command General rules are described in divi It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 74 qenq r ho 2 2 2 3 2 gqgenq RETURN wali PePPe Pere 11718 I gt Ps PL lent PE eu v n e mu ev lt ae This keyword is used to gain information about entities It is especially useful to get the values on particular nodes If the maximum or minimum value in a certain area has to be searched type first m to go in the max min mode and create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the selection rectangle Then the key h high to search the node with the maximum value The minimum is searched with I low The search result is then shown in the xterm and also attached to the node With the u key the 60 last search result can be deleted undo The t key changes into the qtxt mode The atxt functionallity is now available which allows to manipulate the node attached string and its position on the screen Use gadd instead of qenq if you need t
114. ned 7 1 Datasets Datasets are selected with the menu item Dataset A dataset is a block of nodal values These could be displacements due to a linear analysis or for a specific time step during a nonlinear analysis It could also contain other values like stresses strains temperatures or something else To select a dataset make sure that the mouse pointer is inside the menu area Then press the left mouse button and move the mouse pointer over the menu entry Dataset then continue to the right A sub menu pops up showing all available datasets with a leading number and sometimes followed by a dataset value usually time or frequency and a dataset description Move the mouse pointer over a dataset you are interested in and release the left mouse button The dataset is now selected A results Entity must be chosen to see the values in the drawing 16 MENU AREA DRAWING AREA Hardcopy DIPs Hardcopy HELP wel Hardcopy beamp frd Figure 3 structure of the main menu area This Dataset might also contain automatically calculated values like the v Mises stress and the maximum principal stress see Pre defined Calculations and Result Format See also the command ds to control the functionality with the command line 7 1 1 Entity To view data from the dataset its also necessary to specify the entity i e dx for a displacement Dataset It works in the same way as for selecting the dataset
115. never be necessary to use it manually 8 60 plot plot n gt e gt 1 gt p 1 1 8s b S L sh si gt P ra Ppr bt lv lt set gt PE Pr Pe PE Py PE a t This keyword is used to display the entities of a set Entities already visible will be erased The following types of entities are known Nodes n Elements e Faces f Points p Lines l Surfaces s Bodies b Nurb Surfaces S Nurb Lines L Shapes sh and the shaded illuminated surfaces si The entities can be displayed in the following colors White w Black k Red r Green g Blue b Yellow y Magenta m Nickel n metallic grey and turquois t To display the entities with attributes use the type in combination with an attribute second letter For example 51 plot la all will display all lines with their names The attribute d works only for lines plot ld all will display all lines with their division and bias see Dia The division is given by the last two numbers 1 99 and the bias by the leading number If there is more than one number in front of the division the number has to be divided by a factor of ten to get the bias 10130 means a bias of 10 1 and a div of 30 Also the attribute p works only for lines In this case the lines with its end points are drawn plot lp all This is useful to detect the begin and end of all lines If end points are deleted the line is also deleted Therefore special care with end points is necessary
116. ns of the lines of this surfaces see adiv or Jamsbl or chose the structured element type tr6 for this surfaces Or define a certain element size transition with feltyl See especially the command gmshj to fix parts of the surface mesh in a convenient way Generate the tet mesh mesh all tet A target size for the tet elements can be specified if ng vol from the cgx distribution was installed for example mesh all tet 2 5 will define a maximum element size of 2 5 If the model consists of several unconnected parts separate them in single sets and mesh one after the other see seta farem comp Create sets of nodes faces or elements to create boundary conditions if you had not created them already based on the geometry and export them with Bend 121 e Create the input deck for ccx with an editor and start the calculation e Look at the results with cgx cgx results frd B 14 How to check an input file for ccx A quick check of a ccx input file can be done with cgx by reading the file with the option c cgx c file inp After startup all defined node element and surface sets are availabe together with internal sets which group togeter certain entities according to their purposes The following internal sets will be created if appropriate data were found e bou DOF Created if BOUNDARY is found All affected nodes are stored in bou were bou DOF store just the nodes which are con
117. nstrain the translational degrees of freedom of the nodes in set hinge Or send hinge abq spc 12356 will leave just one rotational degree of freedom unconstrained And send move nodes abq spc 1 0 1 will move the nodes in set move nodes by 0 1 in direction 1 send move nodes abq spc 1 ds1 el will use the values from dataset nr 1 and the entity nr 1 for the forced dis placement in direction 1 and send move nodes abq spc nor 0 1 will force the nodes included in set move nodes by 0 1 in the normal direc tion normal to the element faces The components of the vectors at the node positions can be scaled individually facl 3 A subsequent calculation with ccx will move the mesh accordingly After solving the deformed mesh can be used as a new mesh for further calculations With this procedure geometric variations are possible without manipulating the geometry morphing For cfd there is a spc function available which works on faces instead on nodes send Inlet abq spcf 123 78 The degees of freedom 123 define velocities and not deflections Its called spcf and works otherwhise accordingly to the spc command The parameter slide in combination with one additional parameter will create equations that force the selected nodes to move only in a plane If the additional parameter s is specified then the specified nodes will be attached to their element faces All nodes of the affected element faces have
118. nt definitions If the geometry was already meshed the mesh will NOT be deleted If the mesh command is executed again after new assign ments has taken place additional elements will be created elty all he20 assigns 20 node brick elements to all bodies in the set all elty part1 he8 redefines that definition for all bodies in the set part1 elty part2 tr6u assigns 6 node unstructured triangle elements to all surfaces in set part2 elty part2 tr6u 0 5 will do the same but specifies a mesh refinement factor of 0 5 gt 1 coarser than the average boundary spacing lt 1 denser elty part3 te10 assigns 10 node elements to all bodies in set part3 But this works only if NETGEN 4 is installed and the program ng vol is accessible elty part3 tel0 3 5 will do the same but specifies a target size for the elements In this case the modified program ng vol from the cgx distribution must be available Replace the original ng_vol in the NETGEN package and build it again 8 23 enq enq lt set gt lt set gt set lt setname gt gt rec lt x value gt _ lt y value gt _ lt z value gt _ JJ gt 35 cyl x y z lt r value gt _ lt axis value gt _ gt lt tol gt i a h 1 lt value gt This command is used to locate entities from a certain set first provided set and stores them in the second set The following entities are handled nodes points l
119. o save the node in a set for further use To catch more than one entity with one stroke type a all before and create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the rectangle To catch only the entity which is closest to the user type i before Then move the mouse pointer over the entity s and press one of the following keys for Nodes n Elements e Faces f Points p Lines 1 Surfaces s Bodies b Nurb Surfaces S and for Nurb Lines L The position of nodes or points are given in Cartesian and cylindrical co ordinates see figure 6 axyz are the 3 angles around x y and z rxyz are the 3 radii around x y and z In a second row the sets to which the picked entity belongs are listed Press q to quit the command It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry rz X Figure 6 qenq Definition of the cylindrical system 61 8 75 qfil qfil lt radius gt RETURN w 12 tg The command creates a fillet line a radius between two lines who intersect After the command name the value for the fillet radius has to be specified Then for each pair of selected lines a fillet of this value will be created Select lines with the 1 key Press qu to quit the command Currently it works only for
120. o with the optional parameters node dist and angle Node dist is the maximum allowed distance between nodes and angle is the maximum al lowed angle defined by three sequential nodes If one parameter is not fulfilled then the division is halved until the requirements are fulfilled Default values are defined in the file cgx h and can be listed with div without parameters div all auto uses the defaults But they can be overruled div all auto 2 10 0 5 will use a maximum element lenght of 2 the angle between successive nodes is less than 10 degree and the minimum element is only half of the maximum length Some rules must be fulfilled before a geometry is meshable For linear elements ie qu or he8 the sum of all divisions of each surface must be even In case of quadratic elements ie que or he20 this sum must be divisible by 4 without residue Opposite edges of a given surface might have different divi sions For example on the left side of a given surface the division is 8 and on the right side it is only 4 But only two opposite surfaces of a body can use this feature These surfaces are called top and bottom surfaces All other surfaces of this body must have unique divisions on opposite edges 32 8 17 ds ds lt 1 Dataset Nr gt lt 2 Dataset Nr gt lt n Dataset Nr gt a e gt lt entity nr gt up to 4 times ll Ps lt value gt This keyword is used to define one or more and one or more Entity for d
121. offset gt lt to1 gt Prot x y z lt offset gt lt tol gt rad lt p1 gt lt p2 gt lt offset gt lt tol gt Prad Py lz lt offset gt lt tol gt nor lt offset gt lt tol gt This keyword is used to project points with all related geometry or nodes onto a set of surfaces or element faces Alternative a shape can be specified as the target In consequence the set with the entities which should be moved must contain nodes and or points and the target set must contain surfaces and or element faces Several transformations are available For example tra will move points in the direction of the vector dx dy dz onto elements or surfaces included in set2 Alternatively an offset could be specified as well proj setl set2 tra 0 0 5 0 7 55 rot will move points around the axis defined by the points p1 and p2 or around the x y z axis onto elements or surfaces included in set2 proj setl set2 rot p0 px rad will move points radial to the axis defined by the points pl and p2 or radial to the x y or z axis onto elements or surfaces included in set2 Alterna tively a set of lines could be used instead of surfaces as the target set Then the geometry will be moved onto an imaginary rotational surface defined by these lines proj setl set2 rad x nor will move points in a direction normal to the target surface onto surfaces included in set2 An offset might be specified proj setl s
122. p0 px 20 rad will create a copy and will move it radially to the x y or z axis or will create a spherical section if just a single point is defined copy cylinderl cylinder2 rad x 20 copy spherel sphere2 rad pPO 10 mir will create a mirrored copy The mirror plane is placed normal to the direction running from P1 to P2 and placed at P2 copy sectionl section2 mir P1 P2 8 11 corrad gt corrad lt set gt This is a very special command to adjust improperly defined arc lines like in fillets The center points of arc lines included in the set are moved in a way that each arc line will run tangentially into a connected straight line But because the end points of the arc lines are not moved only one side of each arc line will run into a connected line The other side is not controlled and might end in a sharp corner Therefore for each arc line exactly one connected straight line must be included into the set figure 29 Line L002 and L003 are members of set LIN Point D004 is the center point of arc line L003 and will be moved Line L003 runs now tangential Command into line L002 corrad LIN 02 02 L003 L003 D004 ne Figure 4 Effect of the corrad command 8 12 csysa gt csysa lt sysNr gt lt set gt Specifies the displacement coordinate system for each node Nastran only 8 13 cut cut lt pnt nod gt lt pnt nod gt lt pnt nod gt This keyword is used to define a cutting
123. parameter k forces black as background colour 8 119 volu volu lt set gt This keyword is used to calculate the volume and the center of gravity of a set of volume elements If an ccx input file with density data was read then the mass will be also calculated If a dataset is active then an averaged value is calculated The command writes to the tack 8 120 while while lt value gt eq ne lt gt gt lt value gt 96 A command to compare two values If the compare is True the following com mands are executed until the endwhile command is found This procedure is repeated until the compare is False while argl arg2 will repeat the commands between while and endwhile until the nummeri cal value stored in argl is not equal to the nummerical value stored in value arg2 The values are locally converted to float format for the nummerical comparison The eq and ne compare strings and should not be used for num merical values since no conversion to a common format is done Two strings are equal if they have the same length and all characters are equal See also El valu Stack and How to run cgx in batch model 8 121 zap zap lt set gt This keyword is used to delete all entities of a set and the set itself All depending entities will be deleted as well 8 122 zoom zoom lt scale gt lt pix gt lt ply gt lt p2x gt lt p2y gt Thi
124. plane through elements to visualize internal results The plane is either defined by three nodes or points or in case a dataset entity of a vector was already selected by just one node or point The cutting plane is then determined by the direction of the vector displace ments worstPS The menu option Show Elements With Light or the com mands Tucutl Tvievlsur or view volu will display the whole model again and will delete the plane This command is intended for batch mode See qcut for the cursor controlled command 8 14 del del p gt 1 gt 1 gt 10 c 1 s b t S L se sh lt entity gt se0 mesh pic 30 This keyword is used to delete entities the whole mesh see also qdel or a background picture For example del se part will delete the set part but all included entities are untouched The following entities are known Points p Lines 1 Combined Lines c Surfaces s Bodies b Node Texts t Nurb Surfaces S Nurb Lines L Sets se and Shapes sh When an entity is deleted all dependent higher entities are deleted as well Special cases are del 10 set I zero were all lines with zero length in set set are deleted and del sel will delete all empty sets If a background picture was loaded with the read command it can be deleted with del pic See also api on how to delete a set with all its referenced entities 8 15 dist dist
125. r It is advisable to save your work from time to time without exiting the pro gram This is done with the command or with Exit will write all geometry to an fbd file and if a file of this name exists already then the extension of this file will be renamed from fbd to fbb quit closes the program without saving A solver input file can be written with the help of an editor emacs nedit etc If you write a ccx command file then include the mesh the boundary conditions etc with the ccx command INCLUDE After you finished your input file for the solver ccx you might read it by calling the program again with cgx c solverfile inp 12 for a final check All predefined sets are available together with automati cally generated sets which store boundaries equations and more These sets start with the sign For example the set bou stores all constrained nodes where the set boul bou2 bou3 store the constraints for the individual di rections Further the set dep and ind store the dependent and independent nodes involved in equations etc See which sets are defined with the command e prnt se Each line starts with the set index then the set name followed by the number of all referenced entities The sets can be specified by index or name For example if the index of set blade is 5 the following commands are equivalent e plot p 5 e plot p blade The use of wildcards is pos
126. r distance between nodes The option Datasets will provide a plot value over Dataset nr and the option Time will provide value over Time For the later two options it is necessary to first create an animation with either the command ds or the menu option Toggle Dataset Sequence see also graphl and How to generate a time history plot To easily find the element corners the function Toggle Element Edges is triggered automatically 7 10 Orientation 7 10 1 x View To look along the x axis 7 10 2 x View To look against the x axis 7 10 3 y View To look along the y axis 7 10 4 y View To look against the y axis 7 10 5 2 View To look along the z axis 7 10 6 z View To look against the z axis 22 7 11 Hardcopy To create a hard copy during animation it is useful to stop the animation first with the middle mouse button while inside the menu area of the main window and then release one picture after the other with the right button until the desired amplitude or step is reached 7 11 1 Tga Hardcopy To create a window dump in tga format You might use the program convert to convert this format to others 7 11 2 Ps Hardcopy To create a window dump in postscript format The program convert must be installed 7 11 3 Gif Hardcopy To create a window dump in gif format The program convert must be installed 7 11 4 Png Hardcopy To create a wndow dump in
127. re KEY 3 CODE C NUMELEM Number of elements in this block FORMAT Format indicator O short format 1 long format 2 binary format Following records ascci FORMAT 0 1 The following block of records must be repeated for each element The first record initializes an element definition Short Format 1X 1 15 315 Long Format 1X 1 110 315 Values KEY ELEMENT TYPE GROUP MATERIAL Where KEY 1 NODE element number TYPE element type see section Element Types GROUP element group number see command grps 106 MATERIAL element material number see command mats Then the nodes in the correct order have to follow Short Format 1X 2 1515 Long Format 1X 2 10110 Values KEY NODE NODE NODE NODE Where KEY 2 NODE node number Additional lines must follow if more nodes are used Following records binary FORMAT 2 Format 4 int nodes int int and float are ansi c data types Values ELEMENT TYPE GROUP MATERIAL NODE NODE NODE NODE Where ELEMENT element number TYPE element type see section Element Types GROUP element group number see command grps MATERIAL element material number see command mats NODE node number Last Record only FORMAT 0 1 ascii omitted for FORMAT 2 Format 1X 3 Values KEY 10 5 Parameter Header Record Purpose Stores informations related to data
128. result files which can be visualized with cgx There is also an example which allows cgx to be used as an pre processor for dolfyn The mesh can be written in dolfyn format with the command send B 18 Remarks Concerning Duns and Isaac If you intent to create a 2D mesh for the cfd code duns or isaac you have to watch out that all surfaces are created in the same order That means that all surfaces must be defined clockwise or counter clockwise For isaac they have to be counter clockwise if you look in z direction The block structure can be detected only in this case You can check the mesh by simply mesh it with linear shell elements and display them All elements must be either dark or illuminated A later qflp or flip command will not cure wrong oriented surfaces since it only changes the sign in their definition and not the basic edge sequence which is necessary here The results of a calculation can be opened by specifying the parameter duns2d duns3d isaac2d isaac3d and the filenmae without any extention cgx isaac2d RAE2822 See also Program Parameters mesh Bend Watch also the airfoil example in the distrib B 19 Remarks Concerning Nastran The cgx is capable to write the following mesh entities to files see send the necessary key parameters are listed in brackets below e Nodes and Elements 124 e node displacement coordinate system see command fesysa e Single point constraints sp
129. ric entities is changed it is necessary to move the points on which the entities rely It should be noted that faces exist only on free surfaces of the model In addition entities can be grouped together to make sets Sets are useful to handle parts of a model For example sets can be used to manipulate or display a few entities at a time see also How to define a set of entities Also a simple but powerful entity which can store values character strings is available This values can be derived from previous commands or calculated results by using an internal stack Simple calculations can be performed The values can be used to substitute parameters of subsequent commands The user might measure a distance or calculate a distance and use this value to move a part of the mesh Together with a while loop an if case distinguishing com mand and the possibility to use system calls via the sys command elaborated batch files can be written After a mesh is created in cgx it needs written to a file for use with the solver Likewise several boundary conditions and loads can be written to files see also How to connect independent meshes How to define loads and constraints and lsend These files need to be added into the control file for later use in cecx Additional commands material description and so on must be added with the help of an external editor After the analysis is completed the results can be vi
130. rmed into a new local cartesian system which is alligned with the directions of a true cylindrical system In this way the dimensions are maintained for example the displacement in angular direction is not transformed into an angle but into a displacement in tangential direction The transformation from a cylindrical into a cartesian system works accord ingly Therefore successive cyl and rec commands are permitted This command sequence can be used to rotate the model with its datasets in the correct way which means that all results are also rotated Example Choose the desired dataset and an entity with the menu or the ds command Then type trfm cyl z to transform the dataset from a rectangular system into a cylindrical around the global z axis Type tr m rec z to transform from cylindrical which exists after the first call to the rectan gular system which re produces the original values To transform several datasets of the same type at once trfm rec z 1 10000 This command transforms all datasets starting with the first to the last if the last dataset has a number below 10001 but only the ones of the same type as the 1st 8 116 ucut ucut If a section through the mesh was created with the icut or qcut command then this command will delete the cut and display the un cutted structure 93 8 117 valu valu lt name gt lt value gt push pop nr gt P 19
131. s lt value gt ds lt nr gt e lt nr gt film lt temp gt ds lt nr gt e lt nr gt gt lt coeff gt ds lt nr gt e lt nr gt rad lt temp gt ds lt nr gt e lt nr gt gt lt emis gt ds lt nr gt e lt nr gt dflux lt load gt ds lt nr gt e lt nr gt cflux lt load gt force lt f1 gt lt f2 gt lt f3 gt quadlin sur l J dolfyn duns isaac foam lt type gt lt set gt gt cyclic lt set gt lt set gt tx ty tz rx ry rz ex cy cz lt vx vy vz gt send lt dep set gt lt indep set gt nas abq ans cycmpc cycmpcf gt px py pz tx ty tz gt re Pry P rz cx cy cz lt segments gt c u lt NR gt gt l areampc lt dofs t p 1 6 gt c u lt Nr gt lt value gt gt areampc gt slide gt 74 areampc presfit s lt value gt gt P P gap lt vx gt lt vy gt lt vz gt lt tol gt This keyword is used to send data to the file system The following formats are known but not all options for all formats are fully supported so far The geometry format fbd send set fbd or send set fbd c wil
132. s about how to create them You can write a file with basic commands like pnt to create the basis for your construction and read it with the read command Most commands can be used in batch mode This allows the user to write a command file for repeated actions The interactive commands start with the letter q Please make yourself familiar with all of them before you start to model complex geometry After the geometry is created the divisions of the lines can be changed to control the density of the elements Display the lines and their divisions with e plot id all To change the element division use qdivy 11 The default division is 4 With a division of 4 a line will have 6 nodes and will therefore be the edge of two element of the quadratic type Next the type of the elements must be defined This can be done for each of the different sets A new assignment will replace a previous one Delete all previous assignments with e eltylall and assign new types with e eltylall he20 If a mesh is already defined type e del mesh and mesh again with e mesh all Then choose the menu entity Viewing Show Elements With Light to see the mesh lighted Lastly export the mesh in the calculix solver format with e Isend all abd With the send command it is also possible to write boundary conditions loads and equations to files The equations are useful to glue parts togethe
133. s keyword is used to scale the model in the window For example zoom 2 will increase the size of the representation of the model by a factor of 2 A certain region of the model can be specified with two corner points of an imagi nary rectangle The coordinates are relative to the graphic window which has its origin at the left lower corner and as a fraction of the edge lengths For example zoom 0 0 0 5 0 5 will display the third quadrant of the window scaled by a factor of 2 97 9 Element Types Node numbering of the elements and the type numbers used in the Result For mat frd file The solvers might use different node numbering rules 1 2 node beam element be2 type 11 98 3 node beam element be3 type 12 3 node shell element tr3 tr3u type 7 6 node shell element tr6 type 8 4 node shell element que type 9 100 4 O O 8 o O O 1 5 8 node shell element qu8 type 10 4 1 2 4 node tet element type 3 101 8 10 3 5 6 2 10 node tet element type 6 8 7 8 node brick element he8 type 1 102 20 A 7 5 V A 18 6 16 15 l 14 4 11 12 dp o Dra ana de 3 lo 10 9 2 20 node brick element he20 type 4 6 node penta element pe6 type 2 103 15 node penta element pel5 type 5 10 Result Format Listing of the mesh and the nodal results format The data are store
134. section Element Typeslin the appendix for the correct node order For example elem 1 1 2 3 4 qu4 creates a shell element with four nodes 8 19 else else See the if command 8 20 endif endif See the if command 8 21 endwhile endwhile See the while command 8 22 elty elty lt set gt be2 be2d be3 gt tr3 tr3u tr3e tr3s tr3c gt tro tr6u tr6e tr6s tr6c gt qua quee qu4s quec gt que qu8e qu8s quec gt he8 he8f he8i he8r gt he20 he20r gt gt te4 te10 lt parameter gt This keyword is used to assign a specific element type to a set of entities see section in the appendix The element name is composed of the following parts The leading two letters define the shape be beam tr triangle qu quadrangle he hexahedra then the number of nodes and at last an attribute describing the mathematical formulation or other features c axisymmetric e plain strain s plain stress u unstructured mesh r reduced integration i incompatible modes f fluid element for ccx If the element type 34 is omitted the assignment is deleted If all parameters are omitted the actual assignments are posted elty will print only the sets with assigned elements Multiple definitions are pos sible For example elty all deletes all eleme
135. sets ie bondary conditions and loads They should consist of a keyword and a value Format 1X 1 P A66 Values KEY CODE STRING Where KEY 1 CODE P STRING Keyword Value ie FORCE 1000 10 6 Nodal Results Block Purpose Stores values on node positions 1 Record Format 1X 100 C 6A1 E12 5 112 20A1 12 15 10A1 12 107 Values KEY CODE SETNAME VALUE NUMNOD TEXT ICTYPE NUMSTP ANALYS FORMAT Where KEY 100 CODE C SETNAME Name not used VALUE Could be frequency time or any numerical value NUMNOD Number of nodes in this nodal results block TEXT Any text ICTYPE Analysis type O static 1 time step 2 frequency 3 load step 4 0 1 2 Record user named NUMSTP Step number ANALYS Type of analysis description FORMAT Format indicator short format long format 2 binary format Format 1X 12 2X 8A1 215 Values KEY NAME NCOMPS Where KEY 4 Dataset name to be used in the menu NAME NCOMPS IRTYPE 1 3 Type of Record IRTYPE Number of entities Nodal data material independent 2 Nodal data material dependant 3 Element data at nodes not used Format 1X 12 2X 8A1 515 8A1 Values KEY NAME MENU ICTYPE ICIND1 ICIND2 IEXIST ICNAME Where KEY 5 NAME MENU 1 ICTYPE 1 2 4 12 14 scalar vector matrix vector degree tensor degree Entity name to be used in the menu for this comp Type of entity with 3 components
136. sh 94 A white space separated string stored in a valu will be written to the stack in separate sub strings With this method single pieces of a white space separated string can be splitted and stored in separate values stack on opens the stack valu complicatedString push splits the string stored in complicatedString and writes the single pieces to the stack The command valu subString pop 3 loads the 3rd substring into the varible subString The command is able to perform simple calculations like valu result argl arg2 result will store the product from argl and arg2 It is important to un derstand that the command deals only with other values That means that a direct multiplication of two numbers is not possible The number has to be stored in a value before it can be used valu argl 1 24 assigns the string 1 24 to the value argl During the calulation the strings are converted to double precision numbers and the result is stored as a string representing an exponential number The int float and exp convert between integer floating point and exponential format The operator is used to concatenate two strings valu string3 amp string string2 See also Iprntl and How to write values to a file 8 118 view view fi11 line point edge off elem off surf volu gt front back vec off disp
137. sible to search for a certain expression e pinl se Now all sets starting with a in their names will be listed Predefined loads are stored as IDatasetsl to be visualized Sets with the name of the load type CLOAD DLOAD store the related nodes faces or elements Use the command or to visualize entities of sets Then run the input file with ccx The result file frd can be visualized with egx filename frd filename inp were the solver input file filename inp is optional With this file the sets boundary conditions and loads used in the calculation are available together with the results If you have problems doing the above or if you want to learn more and in more detail about the cgx continue with the tutorial and look in the appendix section lips and Hints and Known Problems 5 Program Parameters usage cgx al b bg c duns2d duns3d isaac2d isaac3d foam ng 13 stepl stl filename ccxfile a automatic build mode geometry file derived from a cad file is expected b build mode geometry file in fbd format is expected bg bacEground suppress creation of graphic output otherwhise as b geometry command file must be provided c read an solver input file ccx Abaqus duns2d read duns result files 2D duns3d read duns result files 3D isaac2d read isaac result files 2D isaac3d read isaac result files 3D foam read the OpenFOAM result directory structure 06 re
138. straight lines A curved line can be split and the part which should be used for the fillet can be transformed to a straight form with the ali command s and x key strokes Figure 7 shows on the right side the initial situation and on the left side the created fillet The command has also shifted point P002 to the left Always the end point of the first selected line closest to the intersection is moved The second line gets a new end point 002 02 Figure 7 qfil Based on two intersecting lines a fillet is created 8 76 qflp qflp RETURN w Pe Pe Pra ri Pa 62 This keyword is used to invert the outer and inner side the orientation of shell elements and surfaces The orientation of shell elements or surfaces can be seen by the interior color The outer face reflects light if the face is dark grey it is the back side If Toggle Culling Back Front was selected before then both sides are illuminated To flip the orientation of a surface and all related shell elements select either a shell element with the e key or if the illuminated surfaces are displayed see Fep how to do that select the surface with the s key To see the effect on the elements immediatelly they must have been displayed with the plot command If only elements are in the database flip them with the e key The a key has a different meaning than usually If pressed before a selection then a so called auto mode is activat
139. sualized by calling the cgx program again in an independent session The program is primary controlled by the keyboard with individual commands for each function Only a subset of commands which are most important for post processing is also available through a pop up menu Shaded animations of static and dynamic results the common color plots and time history plots can be created Also a cut through the model can be done which creates a section and it is possible to zoom through the model Skilled users might include their own functions For example someone may need his own functions to manipulate the result data or he may need an interface to read or write his own results format see also cal Both the pre and post processing can be automated in batch mode see also How to run cgx in batch mode 3 File Formats The following file formats are available to write w and or read r geometric entities e fbd format r w this format consists of a collection of commands ex plained in the section Commands and it is mainly used to store geo metrical information like points lines surfaces and bodies All geometry generated by the user is stored in this format But it can also be used to define a batch job which uses the available commands step format r reverse engineered based on some cad files Only points and certain types of lines are supported currently stl format r w this format describes a shape using only triangles
140. t LOAD nam input ANCHORAG nam etc 11 IMPORTANT REMARK After it s possible to modify the mesh into BACON extrusions etc that s why the element hypothesis is not added at the end of the 127 file gt you must define the element definition AFTER the last mesh modification HYP MINDLIN 12 Comments if you ve any remark or any comment or any suggestion to improve this export format please send a mail to paul carrico_at_free fr So far results can not be read C Simple Examples The following listings show simple geometry input files The pictures show this geometry together with their labels and the generated mesh The models were made based on three points Two points defined one axis of rotation and one was the basis of several swepj and merg operations In case of the sphere the surfaces on the pole had to be redefined using only three lines per surface C 1 Disc PNT py 0 00000 1 00000 0 00000 PNT pO 0 00000 0 00000 0 00000 PNT P001 0 70711 0 00000 0 70711 PNT P003 0 00000 0 00000 1 00000 PNT P005 0 70711 0 00000 0 70711 PNT P006 1 00000 0 00000 0 00000 PNT P009 0 70711 0 00000 0 70711 PNT POOA 0 00000 0 00000 1 00000 PNT POOG 0 70711 0 00000 0 70711 PNT POOL 1 00000 0 00000 0 00000 LINE L001 POOI P001 pO LINE L002 P001 P003 po 4 LINE L003 P003 pO 8 LINE L004 pO POOI 8 LINE L005 P003 P005 pO 4 LINE L006 P005 P006 pO 4 LINE L007 P006 pO 8
141. t are considered The value gtol determines the maximum distance between merged nodes and points The parameter nolock will force merging even if the dependent entities like lines are degenerated afterwards For example a degenerate line will have two equal points 8 48 mesh mesh lt set gt fast tet lt size gt l block gt lonly l nolength noangle length angle This keyword is used to start the meshing of the model before using the mesh command the element types must be defined with the elty command Existing elements will not be deleted Therefore its possible to start cgx in the viewing mode v with a mesh alone and then create bodies and fill them with addi tional elements To delete a mesh use the command dellmesh See the command lsend to describe areas for boundary conditions In case a blocked grid for cfd calculations should be generated use the ad ditional parameter block mesh setname block see also the section Remarks Concerning Duns and Isaac The mesh optimizer for structured elements is controlled with the additional parameters nolength and noangle These parameters switch off length and angle optimizations of elements The parameter fast suppresses the Nurbs based meshing and uses the fast coons algorithm for surface meshing Usually all elements are generated in a structured way That means that no holes in a surface or volume are permitted T
142. t of the mesh and will be stored with the mesh in the frd format or nastran format see send For example mata 7 part assigns the material number 7 to all elements in the set part Elements can have just one material number The default number is 1 These numbers will be saved with the mesh if the database is written to the file system with the command send all frd or send all nas If the frd file is used later the material number s are available immediately 8 44 map map lt slave set gt lt master set gt gt gt surf volu xy Paz trx gt ery rz gt ds lt nr gt This keyword is used to map or interpolate values from one mesh to another For example map slave master surf dsl will map the values from Dataset 1 to the nodes of set slave All available 44 datasets will be mapped if no number follows the ds parameter The param eter surf is used for mapping of values from surface to surface 2D to 2D A typical application would be the mapping of pressure The parameter volu triggers the mapping from a 3D mesh to another 3D mesh of the same shape i e temperatures The parameters x y z are used for mapping from 2D to 2D 3D in the indicated direction The parameters rx ry rz are used for ro tational mapping of 2D to 2D 3D around the x y z axis The master nodes must belong to elements Unconnected nodes are not allowed Usually the mas ter nod
143. t set gt RETURN tw 2a 121 Pere Pe PE Pp PI 12s gt 1b gt q This keyword is used to erase entities from a set see also setr To remove entities move the mouse pointer over the entity and press the following keys for Nodes n Elements e Faces f Points p Lines 1 Surfaces s Bodies b Nurb Surfaces S and for Nurb Lines L To catch more than one entity with one stroke type a before and create a rectangular picking area by pressing two times the r key Both strokes define opposite corners of the rectangle To catch only one entity type i before Press q to quit the command It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 84 qseq qseg lt set gt RETURN n p RETURN nr This keyword is used to define a sequence The sequence might consist of nodes or points A sequence is nothing else than a set which keeps the selection order of the entities qseq nodeset will store all selected nodes selected with the n key in the order in which they were selected This set might be used for a subsequent graph command Point sequences are used to define splines If qseq is typed without a setname then an existing line has to be selected with a nummerical keystroke or the space bar followed by two numbers New points 66 will be creat
144. t works similar to the pres option For example send surf abq film 1200 0 5 will assign the film temperature of 1200 and the coefficient 0 5 to all element faces in set surf If datasets with nodal values are available it is also possible to use this values instead of using a uniform value for all faces For example send surf abq film dsl el 0 5 will specify the dataset nr 1 and the entity nr 1 to be used for the tempera ture instead of a uniform value or send surf abq rad ds1 el Ic1 e2 will specify also the dataset nr 1 with the entity nr 2 to be used for the film coefficients instead of a uniform value The front or rear side of the face pos or neg is selected with the or parameter it has to be the last parameter The parameter rad is used to assign a sink temperatures and the emissiv ity to element faces It works similar to the pres option For example send surf abq rad 1200 0 5 80 will assign the sink temperature of 1200 and the emissivity 0 5 to all element faces in set surf If datasets with nodal values are available it is also possible to use this values instead of using a uniform value for all faces For example send surf abq rad lc1 el 1200 will specify the dataset nr 1 and the entity nr 1 to be used for the tempera ture instead of a uniform value or send surf abq rad lc1 el lc1 e2 will specify also the dataset nr 1 with the entity nr 2 to be used for the emissiv ity
145. t zoom gt lt x_zoom gt lt y_zoom gt This keyword is used to read files or commands from a file Most commands can be read from a file but not all of them In general all basic commands to create geometry are understood and up to now this is the only way to read pre defined geometry during run time To read commands from a file like pnt or Seta and so on type read geo file this will eventually replace entities if their names were already in use To pre vent this type read geo file add this forces the program to change the names which are already in use to unused ones Therefore no existing geometry will be overwritten If an ccx input file with loads boundary conditions and sets was read read ccx file inp then the loads and boundary conditions are stored in automatically created sets which start with a leading For example the ccx command BOUND ARY will trigger the creation of the set bou The ccx sets are stored in cgx sets of the same name The load values are stored in Datasets The option nom no mesh read ccx file inp nom is used to suppress the reading of nodes and elements This speeds up the reading and makes sense if the mesh exists and only the set definitions and loads should be read If the parameter add is used read cex file inp add 70 then existing node or element numbers are not overwritten and the program choses new numbers If the ccx input file ends with
146. th qseqj graph seq length The length is calculated as the sum of the distances between successive nodes Since version 1 8 there are two ways of operation One way is to specify the Datasets and the entity as parameters of the command graph set time DISP D1 which will display the displacement D1 of nodes in set vs the dataset time of all DISP Datasets The first and the last dataset can be specified optionally For the traditional way a sequence of datasets must have been selected see Toggle Dataset Sequence After the selection of the datasets an Entity must be specified Then the user could use the graph command to generate the history plot of this pre defined sequence The command graph set nr will display the values on the node positions vs the dataset numbers The command graph set time will display the values on the node positions vs the dataset values usually time or frequency seelNodal Results Block parameter VALUE and the command 39 graph set descr will display the values on the node positions vs the dataset description only makes sense if the description is a numerical value See Nodal Results Block parameter TEXT for the location of this data in the frd file In addition a second data file is written graph out2 which stores the node number and the node coordinates and all values at this node position for all datasets This file can be used to display
147. three to five edges must be defined To create a surface move the mouse pointer over the first line of the first edge and select it with the 1 key number one Tf more than one line is necessary to define the first edge select the following ones one after the other with the 1 key Each marked line is listed in the mother xterm If all lines of the first edge are selected select the lines of the second edge by pressing the 2 key then the third and eventually the fourth and fifth The last selected line must match the first All lines defining an edge will create a combined line seellemb and this combined line will show up in the definition of the surface Then generate the surface by pressing the g key and you might continue with the next surface By default the interior of the surface is defined according to Coons which is a blended function But it can also follow an embedded NURBS surface or a shape To relate a surface to an existing NURBS surface select another surface which already uses the desired NURBS with the S key uppercase or to reference the shape with the h key Then select the target surface with s The target surface is converted to a NURBS surface which is necessary to be meshable with unstructured triangles tr3u tr6u and to be used to generate tet elements The selected shape or nurbs will stay selected until cleared with c An existing surface might be completly redefined wit
148. tion of shapes 8 106 split split lt set1 gt lt set2 gt This keyword is used to split lines surfaces and elements The split command will not generate new surfaces but it will create points along the intersection between the surfaces in setl and the surfaces or shapes or faces in set2 So far only tetraeder elements in set1 can be splitted by surfaces faces or shapes in set2 In this case the splitted tetraeder elements will be replaced by new tetraeder and pentaeder elements It can be used to cut out parts of the mesh to be either used directly in other models or which can be remeshed with tets based on an stl file derived from them Or it can be used to apply damages to a structure See How to map loads in the appendix on how to map loads to the new created nodes from the original mesh 8 107 stack stack onloff free Several functions which return values to the user might now place them also on a software stack used by farealdistldsferra etol length prntfvalulvolu The valul command is able to read this values from the stack and place them in variables The free parameter will empty the stack but does not open or close it See also and while 8 108 steps steps lt value gt This keyword is used to define the number of different colors in the post processor mode The default is 21 8 109 surf surf lt name char lt 9 gt lt set gt gt lt line lcmb gt lt line lcmb
149. ts matching the expression will be listed prnt se N lists all sets starting with N To see the contents of a specific set type prnt se setName In case the stack was activated the value and the coordinates of stored nodes are written to the stack The content of the stack can be listed with prnt st The index and the value are written The value with the highest index is addressed next To print the definition of a line type prnt lineName For elements and nodes is an additional parameter range availabe In this case the range of node or element numbers will be displayed max and min nr and in addition holes in the numbering are detected 54 prnt n setname range If an ccx or abaqus input file was read then it is also possible to print the amplitudes AMPLITUDE in ccx or the material propperties MATERIAL in ccx wildcards can be used prnt amp amplitude name prnt mat material name If an ccx result file was read then the user headers meta data prnt usr and parameters of each dataset can be listed and written to the stack with either prnt par to list all parameters of the active dataset or a single one like for example prnt par STEP can be listed and written to the stack see stack See also fif and while 8 64 proj proj lt set gt lt target set gt lt shpe gt gt tra lt dx gt lt dy gt lt dz gt lt offset gt lt to1 gt rot lt p1 gt lt p2 gt lt
150. ular picking area by pressing two times the r key Both strokes define opposite corners of the rectangle Type s to select surfaces To catch surfaces individually type i before its also the default After selecting exactly six or two opposite surfaces press g to generate the body Press q to quit the command or u to undo the last action It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 69 qcnt gcnt RETURN w n p This keyword is used to define a new center point or node by pressing n or p when the mouse pointer is at the position of a node or a point It is also possible to measure distances between two pixels on the screen Just press the key w on the positions of the two pixels The distance is calculated in the scale of the displayed geometry 8 70 qcut qcut RETURN w q l n p u v This keyword is used to define a cutting plane trough elements to visualize internal results see figure The plane is defined either by picking three nodes select with key n or points select with key p or in case a dataset entity of a vector was already selected by just one node select with key v The cutting plane is then determined by the direction of the vector displacements
151. value The command will recognize points or nodes as equally located when the distance between them is less than gtol The command writes to the tack 8 36 hcpy hcpy ps xwd gif png make 1s 1 clean def xwd 41 This keyword creates a hardcopy in one of the above formats Default is xwd Use the program convert to convert to other formats if needed hcpy ps will create a ps file hcpy make ls will bundle all ps files created in one file using the landscape ls format The ls parameter may be omitted A maximum of six pictures is placed on each page hcpy clean will remove all hcpy files See also the commands Text rot ds max and min which might be used in hcpy batch files These commands should be used in the shown order before using the hcpy command 8 37 help help This keyword prints a short overview of all commands 8 38 if if lt value gt eq ne lt gt lt value gt This keyword is used to compare two values If the compare is True the fol lowing commands are executed If the compare is False the code after else is executed Normal operation continues after endif if argl arg2 will skip successive commands when the nummerical value stored in argl is not equal to the nummerical value stored in value arg2 The values are locally converted to float format for the nummerical comparison
152. values over node positions etc It is up to the user to generate a suitable plot out of this data The graph files can be edited and combined For example to plot one type of value vs another type of value If an cecx or Abaqus input file was read then it is also possible to create time history plots of the amplitudes AMPLITUDE in ccx or the material properties can be displayed MATERIAL in ccx graph amplitude name amp graph material name mat Wildcards can be used to search for amplitudes or materials of a certain expression They are listed with all values if the parameter I follows after amp or mat graph N amp 1 lists all amplitudes which contain the N key in its name with all values The parameters related to datasets can also be displayed with the graph command graph p time DISP HID This command plots the nodal diameters used in a cyclic symmetry calcula tion over the dataset value Only few parameters are written from ccx in frd format Other applications might define its own parameters in frd format see Parameter Header Record for the format of a parameter record 8 33 grps 2 grps 2 The result format allows to assign elements to a group One element can only belong to one group This command creates sets and stores all elements of a 40 certain group in a certain set called grp nr See Element Definition Block how this applies to the result format frd 8
153. will assign the nodes of the set nodes to the user subroutine umpc and will force them to rotate by 4 degree around the x axis Two files are produced The one with the equations has to be included in the model definition section and the boundary file in the step section In case of a rigid body request there are two other parameters in combina tion with mpc available The rotation value is replaced by a nodenr of an independent node which will be created based on either provided coordinates send nodes nas mpc v4711 1 0 0 or by averaged coordinates based on the specified set here nodes which are interpreted as dependent nodes send nodes nas mpc n4711 In both cases an RBE2 nastran or RIGID BODY REF NODE nodenr ccx abq element will be created which connect all this nodes See also asgn on how to predefine a thermal expansion coefficient alfa for this element only nastran Another useful method are so called cyclic symmetry equations These equations are used when just a section of a rotation symmetric part like a disk is modeled These equations force the two cutting planes of such a section to move exactly equal in the cylindrical system If the coordinate system of the displacements for the solver is rectangular xyz then the syntax is send dep indep nas cyempc rx12 cl Here dep is the set containing the nodes of the dependent side These nodes will be replaced by the solver with the independent
154. y 23 7 13 Quit This terminates the program without a save 8 Commands This section is a reference to all commands and their parameters in alphabetic order If a command is typed the mouse pointer must be in the main window figure 2 Only the echo of the input stream is visible in the parent xterm The keywords are not case sensitive but all command parameters are case sensitive Each reference starts with a short description of the command The following syntax is used for these descriptions Known commands and syntax AS Keyword either uppercase or lowercase lt gt Parameter case sensitive Esc ls combination of parameters or optional parameter Remark OR amp AND from to gt command continues in the next line RETURN press the RETURN key Entities with the exception of nodes and elements are referenced by names which can contain letters and numbers Usually one to four characters is rec ommended If a new entity uses an existing name the old definition will be overwritten To overcome this problem alias names can be used An alias name is defined with the sign in front An already defined alias name can be referenced by placing the sign in front For example LINE L1 4P1 4P2 SET will create a line with the alias name L1 and will use the alias names P1 and P2 to define the end points and uses the set SET to define the point sequence between the end points 8 1 anim an
155. y gt lt dz gt a rot lt p1 gt lt p2 gt lt alfa gt a JI Prot x y z lt alfa gt a 1I Prot lt p1 gt x y z lt alfa gt a rad lt p1 gt lt p2 gt lt dr gt a Prad x y z p lt PNT gt lt dr gt a Prad lt p1 gt x y z lt dr gt a mir lt P1 gt lt P2 gt a This keyword is used to create a copy of a set see seta about sets Geometry nodes and elements with their results can be copied The copy of results is usefull to evaluate additional sectors in case of a cyclic symmetric calculation 28 The copy is included in the new set Existing sets are extended by the copied entities if the last parameter a append is provided Several transformations are available For example scal creates a scaled copy the scaling factors fx fy fz can be chosen independently Several transformations are available For example scal creates a scaled copy the scaling factors fx fy fz can be chosen independently copy partl part2 scal 2 PO copy partl part2 scal 1 1 2 PO tra will create a copy and will move it away by the vector dx dy dz and the optional parameter a will assign the new entities to sets were the mother of each entity is included copy setl set2 tra 10 20 30 a rot will create a copy and will move it around the axis defined by the points pl and p2 by alfa degrees copy setl set2 rot

Download Pdf Manuals

Related Search

Related Contents

Zeitschaltuhr „EMT707RCC“  Home Decorators Collection BSZ9113AS Instructions / Assembly  Zenpro SmartControl  Kodak EKTAPRO User's Manual  CORTADORA DE CÉSPED ESTIMADO CLIENTE  BigDog® Mowers A Series General Service Manual  POWER FEED® 10M WIRE FEEDER  Avaya 4600 Series IP Telephones Application Note    Home Dynamix 32RN-8302-450 Use and Care Manual  

Copyright © All rights reserved. Failed to retrieve file