Introduction

Module Master is the preprocessor of POLIGONSOFT designated for setup of three-dimensional models (GM) for simulations in the processor modules (solvers) Euler, Fourier, Hooke etc.

Preparation to the simulation consists of:

  • preparation of the computational area (or GM);
  • preparation of material properties for all domains of the model;
  • preparation of boundary conditions (heat transfer parameters, etc.);
  • preparation of other data depending on the task being solved.

All these actions are performed in the modules of POLIGONSOFT unified in the preprocessors group.


POLIGONSOFT Preprocessors

This chapter describes the process of simulation preparation depending on the case. The Preparing the Modeling Area section describes the general principles of GM preparation, the possibilities of Master module are discussed in details. Calculation of temperature fields in a casting and a mold is the main task of modeling, regardless of the casting technology and the formulation of the problem. Therefore, the steps described in this section are required to prepare any calculation. The next sections describes how to set-up the cases of filling, radiation heat transfer, simulation of the shrinkage porosity and stress.

Next there is the description of the irregular modelling tasks, such as: exothermic sleeves, virtual mold and simulation with the non-coincident meshes and the description of the task settings of the simulation of a special casting technologies (centrifugal casting, high pressure die casting (HPDC), continuous casting and so on).

At the end of the chapter, there is the detailed description of the work with POLIGONSOFT databases.

To start the calculation, the user should perform the following actions in the Master module:

  • import the finite element model of the computational domain, created in the mesh generator;
  • set the right orientation of the model;
  • indicate the belonging of each body of the model to a casting or mold;
  • set casting alloy and mold materials properties;
  • describe the boundary conditions, for example, set the parameters of heat transfer between the bodies of the model and with the environment;
  • determine the initial conditions for the calculation, for example, the initial temperatures of all bodies;
  • set other parameters in accordance with the simulated technology: the location of the melt supply and its speed, the movement or rotation of bodies, etc.;
  • run calculation.

The module starts from the main window with the Master button, or the same command from menu Preprocessor.


Start of the Master Module

The main part of the module window (see the figure below) is occupied by the working area, where the GM is displayed (1). On the left is the Model Tree (2), which displays a list of GM volumes. The Model Tree allows you to control the visualization of the GM and set the most important parameters of the computational domain. Toolbars (3) are installed on all sides of the module window. The menu bar (4) is also installed at the top of the window.

When working with the program, you can use both the commands of the main menu in the upper part of the window, and the buttons on the toolbars. Almost all items of the main menu are duplicated by the corresponding buttons on the toolbars. Further, when describing the functions of the module, their call will be indicated only using the buttons on the toolbars, implying that there are corresponding menu items, as well as hotkeys (for some commands) indicated in the names of menu items and button tips.


Master module window

Basic concepts: Volumes, Types, Indices

The main task of module Master is to provide to the user the possibility to designate the groups of the finite elements with the different properties and their faces with different characteristics. Such identification of the faces and elements is the base of simulation process.

VOLUMES. When importing a finite element mesh of the computational domain, the Master module groups the finite elements of the mesh into volumes. The division of the computational domain into volumes (often called bodies) is usually already indicated in the imported mesh and, as a rule, is a consequence of dividing the original CAD model into volumes (i.e., even before dividing it into finite elements). The finite elements (and their faces) combined into volumes are automatically assigned some indices and types. When importing a mesh, the assignment of indices and types to finite elements will be done automatically, but often not in the way required for a particular task. Therefore, one of the user's tasks is to change the automatically specified parameters to those that correspond to the modeling task. It is also possible to create new volumes while editing the mesh. For more information on creating new volumes, see the Editing Element Indices in the Arbitrary Areas section.

VOLUME TYPES. The finite elements included in the casting meshes have special properties and are elements of the Casting type, while the finite elements included in the mold meshes have other properties and are of the Mold type. This division is due to the fact that in addition to calculating the temperature fields, which is performed in all volumes of the calculation area, additional problems are solved in the Casting type volumes: melt flow, phase transition (solidification), formation of shrinkage cavities, macro- and microporosity, structure, etc.

When visualizing volume types, red is used for casting bodies, and blue for the mold (see the figure below). In addition, another type can be used in the Master module - Empty. This type is used as a buffer for elements that you want to temporarily exclude from the calculation area. Volumes of this type are not visualized, but gray is used for color indexing of the Empty type in the dialog boxes of the Master module. When saving a geometric model file, all elements of the Empty type are deleted.


Volumes Types

When importing a finite element model into the Master module, the first volume recorded in the file receives the Casting type, all other volumes - the Mold type. How to change volume types is described in the Editing Volume Types and Indices section.

VOLUME INDICES. In addition, both the casting and the mold can be non-uniform, i.e. not all finite elements entering into them necessarily have the same properties. Thus, several materials can be used simultaneously in the molds; they can have the cores, cooling plates, etc in their composition. In it’s turn the casting being monolithic and homogeneous by the physical properties, may have the areas within itself that are distinguished, for example, by the value of the initial temperature. In order to take these differences into account in the calculation, elements with the same conditions or properties are combined into groups and assigned numbers or indices.

There are 9 element indices in total, and they are numbered from 1 to 9. To visualize volumes consisting of elements with different indices, POLIGONSOFT uses a specific color scheme (see the figure below).


Volumes indices

The element index (and its color) means that some condition or property (physical parameters, temperature, etc.) will be assigned to this group of elements. These conditions and properties must also be defined in the Master module by connecting the appropriate databases, templates and process conditions. For example, for volumes of the "mold" type, the index of the element always determines the its material. The properties of this material will be assigned to this volume when connected mold material properties to the model. At this stage, it is advisable to connect one of the technological process templates to the Master module in order to view ready-made sets of mold materials and determine volume indices. This is done in the Start Simulation and Parameters window on the General tab in the Mold Properties line or on the Project Overview tab in Master module window.

Accordingly, if the same properties or conditions are assigned to elements with different indices, then these elements will behave in the same way during simulation. Such an indexing system makes it possible to economically and efficiently set a wide variety of simulation conditions, for example, two domains of identical materials can easily be assigned the same thermophysical properties, but different initial temperatures.

When the model is first loaded into the Master module, the volumes are consecutively numbered, i.e. the first volume gets index 1, the second one gets index 2, and so on. If there are more than 9 volumes in the model, the tenth volume will again receive index 1, and so on. How to assign the element indices is described in the section on editing the indices of volumes and arbitrary areas.

BOUNDARY INDICES. The faces of finite elements can be in different conditions: one thing is the face of an element that is on the boundary between the casting and the mold, and quite another thing is the face of an element that is on the outer boundary of the casting, i.e. between the casting and the environment. For this reason, the faces of elements that form boundaries (surfaces) with certain conditions are assigned the same index. Each specified boundary index reflects the processes that occur during the transition from one volume of the computational domain to another (most often, this is heat transfer, for example, from the casting to the mold). Indices are assigned only to the outer faces of the elements of the volumes of the Casting and Mold types, i.e. the Casting-Mold, Casting-Environment and Mold-Environment boundaries. The faces of elements that are in contact with faces of the same type (i.e. Casting-Casting and Mold-Mold) may not be specially indexed, since at the boundary of their contact, no jump of the modeled function (for example, temperature distribution) is assumed and, therefore, there is no need to specify the characteristics of the process of transition through such edges. Such edges are automatically assigned the index 0, they are not available for assigning a boundary index. How to assign boundary indices is described in the section on editing the boundary indices.

There are 19 boundary indices in total, and they are numbered from 0 to 18. POLIGONSOFT has some rules for assigning boundary indices that cannot be violated. In particular, zero indices are assigned to symmetry boundaries, internal boundaries, and surfaces with zero flows (from the point of view of thermal, diffusion, filtration, electrical, etc. flows, all these are boundaries of the same type).

Indices from 1 to 9 can be assigned to both external boundaries and boundaries between volumes (e.g., between a casting and a mold). In this case, even if there is actually nothing beyond such a boundary, solvers will assume that some conventional “semi-infinite” mold is implied behind this boundary. The material properties of such a virtual mold will be used in accordance with the index that this boundary has. This rule is valid for both casting boundaries and mold boundaries – in any case, the virtual volume will be considered a Mold. The problem of modeling with a virtual mold is described in more detail in the Virtual Mold section. This indices also can be assigned to the boundaries between the casting and the ambient with assigned temperature.

Indexes from 10 to 18 can be assigned only to boundaries with the environment. However, if there is actually a body behind such a boundary, the solvers will still consider it to be a boundary with the environment (this can be used if for some reason it is necessary to "disable" the connection between the casting and the mold in the calculation).

Just like for element indices, boundary indices correspond to some physical parameters, for example, heat transfer parameters. They are assigned when connecting one of the process templates in the Start Simulation and Parameters window on the General tab in the Heat Transfer Parameters line. Just like for element indices, a certain color scheme is adopted for edges, used when visualizing boundaries (see the figure below).


Boundaries Indices

It should be noted that at the boundary between two elements of different types there are always two faces – the face of one element and the face of the other element (see the figure below).


Boundary between the bodies (elements)

In general, each of these faces can have its own index, for which unique conditions can be set. But in normal conditions, when modeling casting technologies, the same parameters are set for both faces. Therefore, either such faces have the same boundary index, or they have different indexes, but the parameters assigned to these indices are similar.

When a mesh model is first loaded into the Master module, all faces are assigned default indices: the Casting-Mold boundaries are assigned index 1 on the casting side and on the mold side. The Casting-Environment boundaries are assigned index 17, and the Mold-Environment boundaries are assigned index 18. If one of the process templates is connected to the Master module, the user can see which parameters are assigned to each boundary index. The template is connected in the Start Simulation and Parameters window on the General tab in the Heat Transfer Parameters line. The template is configured in such a way that when using the default boundary indices, the user in most cases receives a fully configured heat transfer model for the selected process.

If necessary, you can change the boundary indices to create additional conditions on the boundaries. For example, some of the outer boundaries that are considered "by default" as boundaries with the environment may actually be symmetry boundaries and should be assigned the index 0.

When first loaded, the first volume of the mesh model is assigned the Casting type, and all other volumes are assigned the Mold type. This is often not the case, and the user assigns volume types independently. After editing volume types, you will need to reassign the default boundary indices (for more information, see the Editing Boundary Indices section).

Toolbars

This section provides a complete listing of toolbars commands . As mentioned above, almost all the menu items are duplicated by the corresponding buttons on the toolbars. Each command is given a short description and a link to the corresponding section.


Standard Toolbar


The most common functions (loading GM, recording GM, etc.) are concentrated on the Standard toolbar in the upper part of the preprocessor window (see the figure below)


Standard Toolbar



The Open Geometry (Ctrl+O) button is used to load the mesh model in different formats. See Loading and saving the model section for more details.



The Save Geometry (Ctrl+S) button saves the geometric model with all the assigned conditions and calculation parameters in the POLIGONSOFT format with the g3d extension. See Loading and saving the model section for more details.



The Model Information (Ctrl+I) button opens a window that contains information about the loaded geometric model. See Model Information section for more details.



The Start Simulation and Parameters (Ctrl+R) button opens a window where all the required simulation data is set and the calculation starts. See Simulation chapter for more details.



The View Results button launches the postprocessor to view the calculation results. See Viewing Results chapter for more details.



The Model Tree button enables displaying the model tree on the left side of the Master module window. For information on using the model tree, see Editing Volume Types and Indices section.



The Language (F3) button cyclically changes the interface language. Press it several times to switch to the desired language.



The Special Colors button opens the Master module settings window, where you can set the background color of the workspace, casting mesh, mold mesh, section plane color, border and element mark colors. See Special Colors. section for more details.

The Settings (Ctrl+Alt+S) button opens a window that contains all the parameters that control the preprocessor work and its individual functions. See Module Settings section for more details.


Process Toolbar


Functions and tools associated with the case setup (materials, process parameters, boundary and initial conditions) are called by buttons:

  • on the Process toolbar (see the figure below), which by default is located along the top border of the Master module window;
  • on the Process group on the Main tab of the ribbon menu.

Process Toolbar



The Material Properties button opens the editor where the materials used in the calculation are created and edited. See Material Properties section for more details.



The Inlets button opens a dialog where you can specify the locations of the melt inlets into the mold, their sizes and parameters. See Melt Inlet section for more details.



The Heat Transfer button opens the editor where heat transfer parameters are created and edited at the boundaries between volumes (body-body) and at external boundaries (body-environment). See Heat Transfer Parameters section for more details.



The Temperature Modes button opens the editor where the temperature modes of the casting and mold volumes are created and edited. See Temperature Modes section for more details.



The Rotation button opens a dialog where you can set the speed and direction of rotation of the model during centrifugal casting. See Centrifugal Casting section for more details.

The Translation button opens an editor where modes for moving casting and mold volumes are created and edited. See the Translation Modes section for more details.

The Dynamic Environment оbutton opens an editor where you can create and edit an environment consisting of several layers with different parameters. See the Dynamic Environment section for more details.

The Nucleation button opens an editor where the parameters for the nucleation of the grain structure at the casting boundaries are created and edited. See the Surface Nucleation section for more details.

The Inheritance button opens an editor where the macrostructure inheritance parameters at the boundaries of the calculation area are created and edited. See the Structure Inheritance section for more details.


Edit Toolbar


The largest number of special functions of the Master module are called by the buttons on the Edit toolbar, which by default is located along the right border of the Master module window.


Edit Toolbar



The Create Shell button opens a dialog for creating a shell of a given thickness around the selected bodies. See Shells Creation section for more details.



The Extrude Mesh button opens a dialog for creating new volumes by extruding a surface mesh at a specified distance. See New Volumes Creation (Extrude Mesh) section for more details.



The Select Contiguous button turns on the surface selection mode. See Editing Boundary Indices section for more details.



The Select Elements by Frame button turns on the frame selection mode. See Editing the indices of the elements in the arbitrary areas section for more details.



The Set Default Boundary Indices button assigns default boundaries for the entire GM. See Editing Boundary Indices section for more details.



The Select Boundaries button opens a dialog for selecting boundaries by index. See Editing Boundary Indices section for more details.



The Deselect All Boundaries (Ctrl+G) button deselects all the selected boundaries. See Editing Boundary Indices section for more details.



The Assign Boundary Index button opens a dialog for assigning an index to the selected boundaries. See Editing Boundary Indices section for more details.



The Select Elements button opens a dialog for selecting elements by index. See Editing the indices of the elements in the arbitrary areas section for more details.



The Deselect All Elements (Ctrl+E) button deselects all selected elements. See Editing the indices of the elements in the arbitrary areas section for more details.



The Assign Volume Index button opens a dialog for assigning an index to the selected elements. See Editing the indices of the elements in the arbitrary areas section for more details.



The Check Mesh Quality button checks the mesh quality according to various criteria and opens a dialog for viewing the report and correcting bad elements. See Mesh Quality Control section for more details.



The Refine Mesh button performs a shift of mesh nodes to improve its quality. See Mesh Quality Control section for more details.



The Split Elements button opens a dialog for choosing an algorithm for remeshing the selected mesh fragment. See Mesh Quality Control section for more details.



The Mark Parameters button opens the preprocessor settings window for selecting parameters and display styles for selected elements and boundaries. See Mark Parameters section for more details.

The Mesh Quality button shows a histogram of the deviation of mesh elements from an ideal tetrahedron. See Mesh Quality Control section for more details.


Transformation Toolbar


Commands changing the geometric model: its orientation, shape, size, position, etc., are located on the Transformation toolbar (see the figure below).


Transformation Toolbar



The Fix Rotation button rotates the coordinate system to its base position (front view), and changes the orientation of the model in space.



The Reflect by X button flips the geometry around the X axis.



The Reflect by Y button flips the geometry around the Y axis.



The Reflect by Z button flips the geometry around the Z axis.



The Scale button resizes the model.



The Copy/Move button copies or moves the visible bodies of the model in the given direction.

The Measure (F4) button opens a tool for measuring the distance between two given points.

For more details about each command, see the Change Model Orientation, Position and Size section.

Module Settings

Settings button on the Standard toolbar opens a window that contains all the parameters that control the work of the preprocessor as a whole and its individual functions. To open the settings window, you can use the key combination Ctrl + Alt + S (see the figure below).


Master module Settings window

The settings are grouped in tabs, according to the list on the left side of the window. To take effect, you should click the Apply button. However, as a rule, these settings do not require constant adjustment. The default values set the maximum comfort when working with the geometric model. The models settings and settings for some special functions specified in this window will be described further in the corresponding sections. This section describes settings that affect the preprocessor operation.

To return to the default module settings, click the Restore Defaults button in the lower left corner of the window and then confirm the selection with the OK button. The settings for the active tab are reset. Selecting Reset All resets all module settings.

General


The General tab contains user interface settings (see the figure below).

User Interface Group

Sets the general appearance of the user interface - tabbed ribbon menu or classic menu. In both cases, you can additionally use toolbars. Default value: Ribbon and Toolbars.

  • Toolbar Icons Size Group. Specifies the size of buttons on toolbars.
    Default value: 16px.


  • Ribbon Tabs Position Group. Specifies where the ribbon tabs will appear above or below the ribbon.
    Default value: Top Position.


  • Language. Determines the interface language. It can be changed "on the fly" by pressing the F3 key on the keyboard, and the language will change to the next one in the list.

General Tab

Simulation Parameters Group

  • Save Calculation State. The parameter determines whether calculation state files with the dump extension will be created and updated in the DUMP folder. These files are updated every time results are written and can take up significant disk space. With their help, you can continue an interrupted calculation from the moment of the last entry. When this option is disabled, the DUMP folder and the files in it are not created. In this case, the calculation can be continued only from the moment the current solver is launched (from the beginning of filling, solidification, etc.). For more details, see the Continuing an Interrupted Simulation section.
    Parameter is enabled by default.

Special Colors


The Special Colors button on the Standard toolbar opens the Settings window on the corresponding tab where the user can set the background color, casting mesh color, mold mesh color, section plane color, edges and element label colors (see the figure below).


Special Colors Tab

To set a color, click on the corresponding colored square. The standard MS Windows Color Picker window will open. Having set the desired color in it, you should click OK button to return to the window for setting Special Colors. After making all the changes, click the OK button.

The section plane color has a checkbox Auto (it is checked by default). This means that when setting up the section plane, its color will be automatically set to black or white depending on the set background color of the window. If the checkbox is disabled, the color specified by the user will be used.


Selection


The visualization of the selected (marked) elements or their external faces is performed using marks - highlighting the edges of the faces and / or a point in the center with color (see the figure below).


The style of marks can be changed by the user. To do this, you can use the Selection Settings button on the Edit toolbar or the Settings tab of the ribbon menu. The Settings window will open, where the Selection tab contains the settings that determine the appearance of the element and face marks (see the figure below).


Boundary selection result
Element selection result
An example of visualization of marked elements and faces


The color of element and face marks is set in the Settings window, which can be opened by pressing the Special Colors button on the Standard panel (for more details, see the Special Colors section).


Mark Parameters Tab

At the top of the tab there is an option Select boundaries together with opposite. If the option is enabled, then when selecting internal mated boundaries of the "cast-mold" or "mold-to-mold" type, a new boundary index is assigned simultaneously to both of its sides. For example, if you assign a new index to a casting-chill boundary on the chill side, the casting-side boundary will automatically get the same index. When this option is off, you can assign different indexes to mating boundaries. Option is enabled by default.

Below are two groups of parameters: Face Marks and Element Marks. With their help, you can achieve the most convenient marks visualization.





Central Point options draw a point at the center of the marked triangle or element. The parameter is disabled by default.




The Offset options control the visual gap between selected objects. When the gap is zero, the edges of the object are painted directly, and when the gap is increased, the edges of the marks are drawn inside the selected object. Specified as a percentage.



The Surface Angle parameter specifies the criterion for determining whether adjacent feature faces belong to the same surface. If the dihedral angle between adjacent faces is equal to or greater than the specified value, the faces are considered to be on the same surface. This allows you to select faces not element-by-element, but immediately with surfaces. For example, if you set the parameter equal to 180°, the faces that lie in the same plane will be selected, and if the parameter value is 0°, all the edges of the object are selected. The parameter is set in degrees from 0 to 180. The default value is 150°.

The Element Critical Angles parameter specifies the minimum and maximum dihedral angles in degrees. These values are used when analyzing the finite element model for quality. The default limit angles are set to 6.5° and 175°, which usually do not need to be changed. For more information on checking the mesh quality, see Mesh Quality Control section.

Preparing the Modeling Area

Loading and saving the model


The creation of three-dimensional geometric model and its dividing to the finite elements are the complex processes requiring the specialized software. The programs performing the dividing operation save the results of its work in the original file formats. That’s why all models prepared in the external mesh generators should be converted from their formats to the internal format of POLIGONSOFT.

Initially, the working field of the program is free, and most of the buttons on the toolbars are unavailable. And this is natural, since the source file containing the geometric model (GM) has not yet been loaded into the program. Recall that from the point of view of numerical simulations by the finite element method, the geometric model is a finite element mesh.


In order to load the GM, you should click the Open Geometry, button, usually located in the upper left corner of the window on the Working commands toolbar. You can also use the drop-down menu of the "File" menu item or the Ctrl+O key combination.

Note that the menu items duplicate all buttons on the toolbars. In the description we will go from the buttons of the toolbars, without mentioning each time about the presence of a similar item in the submenu and the hot keys noted there.

After clicking the Open Geometry File button, the Loading Options window will open.


Loading Options dialog

In this window, use the radio buttons to select the file type (1) that you want to download. The possible file types cover formats of almost all modern programs intended for dividing the geometric models into the finite elements: FEMAP, HyperMesh, NASTRAN, PATRAN, ANSYS, CADDS, NX, CATIA, ProCAST and SALOME (supplied with POLIGONSOFT). Selecting the POLIGONSOFT type means that it will load a file that has already been saved in POLIGONSOFT format.

Then you can specify the units in which the coordinates of the nodes in the downloaded file (2) are recorded or specify the unique scale factor. In various CAD-systems and mesh-generator programs, the coordinates of points of the geometric objects are specified in different units of length: millimeters, inches, meters, etc. POLIGONSOFT relies on the SI system, and measures all coordinates in meters. Therefore, in order to observe the true size of the geometric model, the corresponding coordinate resimulation should be performed. For example, in the transition from millimeters to meters, all coordinates should be multiplied by 0.001. To change the current scale factor, you should normally set the text cursor in the input field and enter the corresponding number on the keyboard. To check how the program responded to the introduction of the scale factor, you can check the size of the model on the Information Panel. (It is called by a special button, which will be described below.)

After clicking the OK button, Open Geometry File window will open, in which you can find the file to be loaded on the disks.

If you select the Define by Extension radio button in the Loading Options window, it will be possible to load any file with a finite element mesh, since the filter All files (*. *) will be used. When loading, the units of length previously set for the selected file type will be used. If you try to download a file that does not match any of the used formats, an error message will be displayed.

Files can also be loaded by dragging them onto the working area of the preprocessor window.

In POLIGONSOFT there is a possibility of inheriting the parameters specified during the preparation of the previous simulation. This is very convenient when you need to simulate the same casting, but with a different gating-feeding system under the same (or close to the same) conditions. In such cases, most of the parameters such as: casting alloy, mold materials, heat transfer parameters, initial temperatures, etc. remain unchanged and can be translated to a new geometric model. Therefore, after loading the mesh, a message is displayed (see the figure below), in which it is proposed to apply the parameters of the previous calculation to the loaded geometry (Yes button), or load the default calculation parameters (No button).


Selecting the option for loading simulation parameters

After loading the model, an image of the modeled object will appear in the workspace, and most toolbar buttons will become active. By default, the geometric model is displayed in casting and mold surface mode, or as configured by the user during the previous preprocessor run. Managing geometric model visualization is discussed in the Common Elements of 3D Modules chapter.

In the preprocessor Master it is possible to add other GM from previously created geometry files to the loaded model. This can be convenient if some of the GM bodies involved in the calculation do not change in the process of technology development. For example, when simulating vacuum casting, a furnace with heaters is involved in the calculation. This GM can be added to any model of the casting block that is planned to be filled in this furnace (see figure below). Of course, both GMs must be performed in a single coordinate system. To add a GM, select the Append Geometry ... command from the File menu and load the required file.


Example of adding a furnace model to a casting model

A list of previously opened geometry files is stored in the File-Recently Open Files menu.


The Save Geometry command will save the geometric model with all the assigned conditions and calculation parameters in the POLIGONSOFT file with the g3d extension, in which it is ready for calculation. To call the command, you can use the key combination Ctrl+S.

After clicking the Save Geometry button, a standard dialog will open, in which you can specify a file name and specify its location on disk. By default, in the Save Geometry dialog, in the File Type field, the "Geometry Files (* .g3d)" filter is set. When saving with this filter, a geometry file will be saved in the specified location and a DATA folder will be created, into which files with additional data necessary to start the simulation will be written. If you need to run the calculation on another computer or in a different folder, copy the g3d file and the DATA folder to a new location, open the geometry file in the Master module and all the simulation parameters will be updated.

It is also possible to save the geometry file without additional data, i.e. only a finite element model without reference to the given simulation parameters. To do this, in the save dialog, change the file type to "Mesh only (* .g3d)". In this case, only the g3d file will be saved, no DATA folder will be created.


Change Model Orientation, Position and Size


In POLIGONSOFT, the gravity vector is always directed against the Y-axis and this is important for the correct calculation of shrinkage porosity. This is not always taken into account when creating a CAD model of a casting and a mold, therefore, after loading the mesh model into the Master module, it may be necessary to change its orientation. In this case, it is necessary to fix the visual angle of rotation, i.e. enter it into the real coordinates of the GM. The corresponding commands are located on the Transformation Toolbar.

Remember that when you rotate the object (s) in the working area (see the Common Elements of 3D Modules section), the view or the viewing angle of the model changes, but the orientation in the coordinate system remains the same, because the coordinate system rotates with the object.

Click the Reset Rotations button on the Rotations toolbar to set the model to a base position: the X-axis is horizontal, the Y-axis is vertical, and the Z-axis is perpendicular to the plane of the screen. This position corresponds to the front view.

In the base position, it is clear how the model is oriented in the POLIGONSOFT coordinate system. For example, on Fig. below (left) you can see that the loaded model is oriented incorrectly. For convenience, you can turn off the display of form volumes in the model tree. Use the buttons on the Rotations toolbar to set the correct position of the model. For example, for the case shown in Fig. below (left), you need to double-click the Rotate Across X 90° button to rotate the model 180° around the X axis.

Do not use the mouse to rotate the model! Adjusting the exact position with the mouse "by eye" is not possible. Use the buttons on the Rotations toolbar.

After the model is in the correct position on the screen, click the Fix Rotation button on the Transformation toolbar.

When you click on the Fix Rotation button, the position of the object in the working area will remain unchanged, but the coordinates of the mesh nodes will be recalculated, and the coordinate system will unfold to the base "front view" position (see the figure below on the right). Now the direction of the gravity vector (against the Y-axis) corresponds to the technological process.


Front view before rotation
Front view after rotation
The result of the coordinate system rotation

In some cases, additional editing of the model may be required, for example, its scaling or mirroring about one of the axes.


Mirroring the Model


In some cases, it is necessary to mirror the geometric model along an axis. This means that the coordinates of all points along this axis will be reversed. To do this, use a group of three buttons on the Transformation toolbar.




The Reflect by X command flips the geometry around the X axis.



The Reflect by Y command flips the geometry around the Y axis.

The Reflect by Z command flips the geometry around the Z axis.

Press the desired button to get the result. If the wrong button was selected by mistake, you must press it again to return the model to its original position. In Fig. below is the result of flipping the model along the Y axis.


Before reflection
After reflection
The result of flipping the model along the Y axis


Scaling the Model


Sometimes you need to change the dimensions of the model. Most often, such a need arises if the units of measurement were specified incorrectly when importing the mesh (for more information, see the Loading and saving the model section).For the dimensions of a model, see the Model Information section.


The Scale button on the Transformation toolbar opens a dialog where you can set the scaling factor of the model (see the figure below).


Scale dialog

Enter the desired value in the Scaling factor field and click the OK button. You can control the result in the Model Information window.


Copying and moving model volumes


You can create copies of one or more model volumes or move these volumes in a given direction. This can be useful, for example, when modeling the filling of several identical molds standing side by side. In this case, there is no need to create mesh models for all molds.


The Copy/Move button on the Transformations toolbar opens a dialog where you can set the direction of copying or moving the visible volumes of the model (see figure below).


Copying the casting and shell volumes with a shift along the x and z axes

In the left part of the dialog, select the action to be performed. In the right part of the dialog, you need to enter the components of the displacement vector in mm. Clicking the OK button will complete the operation.

Note that the action is performed only for visible objects (visibility is set in the model tree). It should be noted that the intersection of grids is not allowed.

Measurements


The Measure (F4) button opens a tool for measuring the distance between two given points (see figure below).


Measuring the distance between two points

To put a point, you need to right-click while holding down the Ctrl key. After two points are set, the distance between them is shown in the drawing area. The dialog displays the length of the segment between two points and the length of its projections on the coordinate axes.

To remove a point, you need to click the left mouse button while holding down the Ctrl key.


Use the Reset button (Ctrl+D) to remove all points and reset the dialog state.


Materials Assignment


At the beginning of the loaded model set-up for calculation, it is advisable to load all the necessary materials into the project: casting material and mold materials. This can be done in the Start Simulation and Parameters window on the General tab by filling in the lines Alloy Properties and Mold Properties (see the figure below).


Start Simulation and Parameters window

At the end of each line there are icons for necessary operations are performed.




The button is located in the Mold Properties line and loads files of the pmms type (or the obsolete bdf type) into the project, which can contain up to 9 mold materials used for a particular casting technology. Such files can be loaded from the POLIGONSOFT database and edited by the user (see later in this section).



The button is located in the "Alloy properties" line and opens the material editor, where the user can load the material from the POLIGONSOFT database or enter his own data. See the Material Properties section for details.

The button removes data from the project and clears the line.

Another way to assign materials is the Materials tree, which is located on the Project Overview tab (1) in the Master module main window (see the figure below). To set the casting material, press  button in the corresponding line (2). The material editor will open and you can import the required material from the database (see the Material Properties section for more details).


Material tree on the Project Overview tab

МThe mold materials can be specified in the same way as the casting material, i.e. importing each material for each mold volume index directly from the material database. To do this, successively press  buttons opposite each desired volume index (3) and load the material through the Material Properties editor. Volume indexes assigned to mold bodies are marked with a "(+)" at the beginning of the line. If desired, you can save the loaded set of mold materials to a pmms (POLIGONSOFT Mold Materials Set) file and then use it in other calculations, to do this use the button (4). Such files, created for different cast technologies, are in the POLIGONSOFT database, and they can be loaded with a button (loading of obsolete bdf files is supported).


Editing Volume Types and Indices


As described in Basic Concepts: Volumes, Types, Indexes section, in the imported mesh, finite elements are grouped in volumes according to the casting and mold design. Preparing the calculation, all volumes of the model should be assigned one of two types: "casting" or "mold". In addition, each volume is assigned a number or "index", which means that this volume is assigned some condition or property (physical parameters, temperature, etc.). These conditions and properties are set through the connection of the corresponding databases, templates and process conditions. For example, for "mold" volumes , the volume index always determines its material. The properties of this material will be assigned to the volume by the index number from the loaded mold materials property set. Below, in the table, there is a list of properties and conditions that can be assigned to the volumes of the model through the volume index. Note that this list is not the same for casting and mold volumes.


Property / Condition Casting Mold

Material Properties

 

+

Temperature Mode

+ +

Translation

+ +

Separation of the volume mesh

+ +

Filled volume

+

 

Radiation Area

+ +

Empty for stress

 

+

Removable element of the
gating system

+

 


Positions that are linked only to the volume index and do not depend on its type are highlighted in green. This means that by assigning the same volume indices to the casting and mold volumes, you can set them the same conditions, for example, the temperature mode or translation. This can have undesirable consequences, so you should carefully consider what conditions are associated with the selected index.

When importing a mesh created in a mesh generator, the first volume recorded in the file gets the "casting" type by default, all other volumes - the "mold" type. At the same time, the volumes receive continuous numbering, i.e. the first volume gets index 1, the second - 2, and so on. If there are more than 9 volumes in the model, the tenth volume will get index 1 again, and so on. In most cases, the "default" settings do not correspond to the real situation and you need to manually specify the Type and Index parameters for each model volume.

Before you start editing the types and indices of volumes, you should connect one of the templates of the technological process (for example, mold materials file) to the Master module in order to establish a connection between the indices of volumes and the materials of the mold. This is done in the Start Simulation and Parameters window on the General tab in the Mold Properties line(see details in the Сommon Simulation Parameters section). System templates for different technological processes are located in the %POLIGONSOFT%\Database folder.

In order to reassign indices or volume types, the model tree is used, which is displayed by clicking the Model Tree button in the Standard toolbar (button pressed by default).

The tree contains a list of all the volumes available in the simulation area (see the figure below). The volumes are distributed in three tree branches according to the type (Casting, Mold or Empty). By default, all volumes have names according to the type to which they belong and consecutive numbering (for example, "Casting 1" "Mold 2", "Mold 3", "Empty 4").


Model Tree in the Master module

Each position in the tree has a switch that controls the visibility of the volume or the entire group of volumes. When the switch is off, the selected volume or group will not be displayed in the working area of the module. When the switch is on, a volume or a group of volumes is displayed on the workspace. The display mode is set by the buttons on the Modes panel (see details in the Model Display Modes section). The display of the volumes of the Empty group is always disabled and cannot be enabled. To display volumes from the "Empty" group, you need to change their type to Casting or Mold. Each volume in the tree has a color indicator indicating the volume index in accordance with the color scale adopted by POLIGONSOFT (for more information on the indices of elements (volumes), see the Basic concepts: Volumes, Types, Indices section).

To the right of the name are condition indicators. They become active if any conditions associated with its index are assigned to the volume. By clicking on the active indicator, you can edit the specified condition, this is described in the corresponding sections below.

Some of the commands described below are available through the context menu of the tree (see the figure below), the content of which varies depending on the mode of object selection - one at a time or a group one.


Context menu of the Model Tree

To go to volume editing, double-click with the left mouse button on its name in the model tree, or on the volume itself in the working area of the module. This will open the Volume Parameters window (see the figure below). You can change the parameters of not only each volume separately, but also a group of volumes. To do this, select the required volumes in the model tree, and then open the volume parameters window through the context menu of the model tree or by pressing the Enter on the keyboard.


Volume Parameters window
(Investment casting technology template used)

At the top left, in the Name field, the name of the edited volume is indicated, which can be changed. Three radio buttons on the left side of the window allow you to change the type of volume, and the nine radio buttons on the right - its index. The plus sign (+) to the right of the index number indicates that this index has already been assigned to some volumes of the model. This does not prevent you from applying this index again, the model may contain several volumes with the same index (for example, if they are made of the same material).

The figure above shows the view of the window for editing the volume parameters when the volume type "Mold" is selected. In this case, the volume index always specifies the volume material of the mold, so the material names are written next to the indexes. The list of materials and the correspondence of these materials to indexes are set by the user on the Project Overview tab independently by loading materials from .pmat files or by loading a process template, i.e. mold materials file .pmms (for details, see the Materials Assignment section). If the volume type is Cast or Empty, the volume index list does not contain useful information. In general, casting bodies can have any indices , because usually they are not given special conditions or special technological parameters. However, in some situations (for example, when using internal chills), the volume index of the casting can be important. This is discussed in the relevant sections.

Changing the volume parameters is completed by pressing the OK button. If the volume type was changed during editing, a message will be displayed prompting you to assign default boundaries (see the figure below). This happens because after the volume type changing , the numbering of its boundaries contradicts the accepted "standard" indexing of the outer borders (by default, these are indices 17 and 18). The default numbering of boundaries is discussed in the Basic concepts: Volumes, Types, Indices section. The operation of assigning the default boundaries usually facilitates the preparation of the model, but in some cases it is undesirable (for example, when deleting "bad" elements, see the section "Mesh Quality Control" about this). If you agree to assign the default boundaries, you need to click the Yes button, otherwise, click No.


Message prompting to assign default boundaries

If you want to change only the type of volume or volumes (without changing their indices), you just select these volumes in the tree and drag them to the desired group.

By default, the "Master" establishes none-ideal contact at the casting-mold interface and ideal contact between the mold volumes. In other words, the casting and mold are FE mesh with coincident nodes at the contact boundary. Bodies of the "Mold" type represent a equivalent mesh, the elements of which can have different properties, depending on which index is assigned to them. You can split the meshes of mold bodies so that a mold-to-mold boundary forms between them, where you can set the heat transfer conditions (none-ideal contact). To do this, select the detachable volume in the model tree, double-click the window to edit the volume parameters, and tick the Detach Index checkbox located under the Volume Index group. After pressing the OK button, indices corresponding to the required heat transfer parameters should be assigned to the boundaries of the separated volume (the procedure for assigning these parameters is described in the Heat Transfer Parameters section). The separated volume has an interactive indicator (see the figure below), i.e. you can separate volumes directly in the model tree by clicking the mouse in the corresponding cell opposite the volume.

It should be noted that when detaching the mesh of a volume with a selected index, all volumes with the same index and type will be automatically detached (in the example under consideration, all volumes of the mold with index five (yellow) will be detached). In this case, the boundaries between the separated volumes and the remaining bodies of the mold will remain zero (they have index 0 by default). If such boundaries are not assigned new indices with the required heat transfer parameters, i.e. index 0 remains assigned to them, they will be processed by the solvers as adiabatic (condition of no heat transfer).

Casting volumes can also be detached to solve highly specialized problems. In this case, the same rules for assigning boundaries described above apply to them.


View of the Model Tree with separated volume (yellow)

Additional volume control commands are available through the context menu of the tree:

  • Merge Volumes. The command combines the selected volumes into one, provided that there is a contact between them. If there are the volumes among the selected one that do not touch any of the group's volumes, a corresponding message will be displayed. At the same time, if the merger process continues, the adjacent volumes will be merged. During the merge process, a window for editing volume parameters opens, where you can specify the name, type and index of the new volume, and, if necessary, separate it. By default, the new volume is being assigned a name consisting of the names of the volumes to be merged, for example, mold1-mold2-mold3. You can use the key combination Ctrl + M on the keyboard to call the command.
  • Rename. The command opens access to editing the name of the volume selected in the tree. You can use the F2 key on the keyboard to call the command.
  • Delete. The command removes the selected volume from the model permanently. You can use the Delete key on the keyboard to call the command. If you are not sure about the need to delete the volume, you can temporarily move it to the Empty group, it will remain in it as long as the file is open.


To see the results of the Type assignment, you should enable the coloring of the model by volume types by clicking the Show Volume Types button on the Modes panel.

To see the results of assigning the Element Indices, you should enable the model coloring by the volume indices by clicking the Show Volume Indices button on the Modes panel.

More details about display modes see in the Model Display Modes section.


Editing Element Indices in the Arbitrary Areas


Sometimes it is necessary to change the model, for example, delete fragments of some volumes or assign new properties to some arbitrary area of it. Typically, this requires changes in CAD and create a new FE model in the mesh generator, but in some cases, these operations can be avoided. To do this, it is possible to create a new volume from an arbitrary group of elements and assign the necessary parameters to it: type and index, in Master module. The concepts of "type" and "index" of an element are described in the Basic concepts: Volumes, Types, Indices section.

To select elements, use a group of buttons on the Edit toolbar (see the figure below). Some of the buttons in this group may be inactive if there are no selected elements in the model.


Buttons for element indices edition

To create a new volume within the existing model, you need to select the required group of elements and assign new parameters to it.


The easy way to select elements is to use Select Elements button. This opens the dialog box (see the figure below), where you can use the checkboxes and radio buttons to describe the group of elements that you want to select.

Select Elements dialog

To select elements of one or more indices, you need to check the corresponding lines in the right part of the window (1). Only those volume indices that are used in the model are available for selection. The radio buttons (2) indicate where to select the elements of the selected indices - in the casting or in the mold. Depending on the selected volume type (2), the list of volume indices available for selection (1) changes. Radio buttons (3) set how to select elements, selected (1) indices and located in specified (2) volumes, in relation to the section plane. Either all elements satisfying conditions (1) and (2), or their part located before or after a given section plane can be marked. To use condition (3), you must first set the section plane to the desired position. How to control the section is described in the Section Plane section.

Finally, the radio buttons (4) set the operation on a given group of elements - select or deselect.

It should note, that an element is considered to be located on any side of the section plane if all four of its vertices are located there.

After the group of elements is set, click the OK button in the Select Elements dialog to apply the selection operation.


Another quicker way to select elements is to select it by frame. To activate this method, you should click Select Elements by Frame button. The button will remain pressed.

It should be noted that the selection of elements in an arbitrary way, without relying on the types and volume indices assigned to them, can be a rather laborious procedure.

When the pressed button Select Elements by Frame, the combination Shift+LMB allows you to stretch a frame of the required size in the working area (see the figure below). All elements that are entirely within the frame area will be selected. The selection result is shown in Fig. below. This operation can be repeated several times to obtain the desired result.


a)
b)
Selection of a group of elements with a frame: a) selection area; b) result

The Shift+RMB combination allows you to set an area with a frame to deselect a part or all selected elements (see the figure below). All elements that fall entirely within the frame area will be deselected. This operation can be repeated several times to obtain the desired result.


a)
b)
Canceling the selection of a group of elements with a frame: a) selection area; b) result

The appearance of the marked elements (tags) is determined by the settings (see the Mark Parameters section).


The selected elements can be assigned a new index and other parameters. To do this, click the Assign Volume Index button on the editing panel. The Volume Parameters dialog will open (see the figure below).

Volume Parameters Tool

In the Volume Parameters dialog, you can assign a name, type and index to a group of selected elements. After pressing the OK button, a new volume will appear in the model tree, and the selection will be canceled.

Other operations with volumes are described in the Editing Volume Types and Indices section.


All available marks can be removed by clicking the Deselect All Elements button on the Edit toolbar or using the Ctrl+E hotkeys.

To see the result of assigning indices to selected elements, you can turn on the GM visualization by volume indices by clicking the Show Volume Indices button on the Modes panel. This is described in more detail in the Model Display Modes section.

Editing Boundary Indices


All faces of the elements of each GM volume numbers or indices are assigned . Faces forming surfaces (boundaries) with the same heat transfer conditions are combined under one index. It can be both internal body-body boundaries, and external body-environment. In total, 19 indices can be assigned to the faces of finite elements, they are numbered from 0 to 18. Index 0 is automatically assigned to the inner faces of elements of the same type (casting-casting and mold-mold), it is also assigned to the boundaries of symmetry. Indexes in the range from 1 to 9 can be assigned to the boundaries between bodies with imperfect contact (casting-mold, mold-mold) or external boundaries. Indexes in the range from 10 to 18 can only be assigned to the outer boundaries of bodies (casting-environment and mold-environment boundaries). For more details see the section «Basic concepts: Volumes, Types, Indices».

According to the index, boundary conditions are assigned at the border, which are connected to the project through heat transfer files, temperature modes and other files of simulation conditions. For example, to simulate any casting technology, you need a heat transfer parameters file that defines the heat transfer conditions at the boundaries between bodies and at the outer boundaries. When loading the model, the Master module automatically assigns indices to the boundaries according to a certain scheme, which is suitable for modeling most casting technologies: the boundary between the casting and the mold gets the index 1, the boundary between the casting and the medium gets the index 17, the boundary between the mold and the environment - the index 18. The foundry technology templates supplied with POLIGONSOFT are configured for this scheme. To understand the relationship between the boundary index and the corresponding boundary condition, before you start editing the boundaries, you need to connect the * .afo heat transfer file from the POLIGONSOFT database to the project. The %POLIGONSOFT%\Database folder contains such files with heat transfer settings for different technological processes. See Сommon Simulation Parameters section for details how to connect a heat transfer file to a project.

Quite often, automatically assigned boundary conditions (boundary indices) do not reflect the real situation. For example, for sand mold casting, the boundary between the casting and the mold will automatically get index 1. When using a heat transfer file for this type of casting, the heat transfer conditions at boundary 1 will be set correctly. But if a steel or cast-iron chill is installed in the mold, the boundary between the chill and the casting will also have an index 1 by default. Obviously, this is wrong and this situation should be corrected. In the file of heat transfer parameters * .afo for sand casting, there are preset conditions for heat transfer between the casting and the chill, but they have a different index. Therefore, you should select the surface between the casting and the chill and change its index to the correct one.

Operations related to the selection of elements faces are performed using the group of buttons on the Edit toolbar.


Buttons for working with boundaries

Automatic boundary indexing can be performed by clicking the Set Default Boundary Indices button. For more information on default (or standard) boundaries, see Basic concepts: Volumes, Types, Indices section. After pressing the button, a warning is displayed (see the figure below).


Warning when assigning default boundaries

The offer to load default boundaries automatically appears when changing the volume type.

The user can make changes to the default boundary indexing scheme by introducing boundaries with other indices, if necessary, in order to set special boundary conditions for some sections of the boundaries. There are two ways to do this in the Master module.


The easiest way to select borders is to use the Select Boundaries button. This opens a dialog box (see the figure below), in which you need to set the signs of the boundaries that need to be selected.


Select Boundaries dialog

To select the boundaries of one or more indices, check the corresponding lines in the right part of the window (1). Only those boundary indices that are used in the model are available for selection. Borders with indices 10-18 are checked using a special edit box with arrows. Radio buttons (2) indicate where to select the boundaries of the selected indices - in the casting or in the mold. Depending on the selected volume type (2), the list of boundaries indices available for highlighting (1) changes. Radio buttons (3) set how to select boundaries, selected in (1) indices and located in specified (2) volumes, in relation to the section plane. Either all boundaries satisfying conditions (1) and (2), or their part located before or after a given section plane can be selected. To use condition (3), you must first set the section plane to the desired position. How to control the section plane is described in the Section Plane section. Finally, the radio buttons (4) specify the operation itself over the given group of borders – select or cancel selection.

The additional checkbox Only visible (5) applies all filter settings only to those borders that are currently displayed in the workspace. If the setting is enabled, the borders of hidden model elements (whose display is disabled) will not be selected.

It should be note that the face of an element is considered to be located on one or another side of the section plane if all three of its vertices are located there.

After the conditions are set, click the OK button in the Select Boundaries dialog to perform a selection of the boundaries.

Another way to select boundaries is interactive, using the left mouse button and holding down the Ctrl key. In this way, you can select either single element faces or entire surfaces.


Switching between selection modes is carried out by the Select Contiguous button on the Edit toolbar.

If the button is not pressed, the outer edges of the elements are selected one at a time (see the figure below). When the Select Contiguous button is pressed, a mouse click with the Ctrl button pressed will select the entire surface (see the figure below). To deselect individual faces or surfaces use the right mouse button while holding down the Ctrl key.


Selecting individual faces
Surface selection
Boundaries selection

When the Select Contiguous button is pressed, a surface is considered to be a set of adjacent faces, where the dihedral angle between adjacent faces is not less than the specified value in the range from 0 to 180 degrees. For example, if you set this angle to 180°, the adjacent faces that lie in the same plane will be selected, and if the value is 0°, all the faces of the volume will be selected. The default surface angle is 150°.


The angle value for the surface, as well as the appearance of the marked boundaries, can be set in the Settings window on the Mark Parameters tab, which is opened by the Mark Parameters button on the Edit toolbar (for more details see Mark Parameters section).

The marked boundaries can be assigned a new index. To do this, click the Assign Boundary Index button on the Edit toolbar. The Boundary Index dialog will open (see the figure below).


Boundary Index dialog

Select the desired item in the list, then press the OK button. The marked boundaries will be assigned the selected index and deselected (see the figure below).


The result of replacing the index of the selected boundary part

All borders can be deselected either by clicking the Deselect All Boundaries button on the Edit toolbar, or by using the Ctrl+G keyboard shortcut.

To see the result of assigning boundary indices, you can turn on the GM rendering by boundary indices by clicking the Show Boundary Indices button on the Modes toolbar (see details in Model Display Modes section).

Mesh Quality Control


Mathematical models used in POLIGONSOFT impose rather stringent requirements on the quality of meshes, therefore, quality control of elements is an important task of preparing a model for calculation. Ideally, all elements of a tetrahedral finite element mesh should have the shape of a regular tetrahedron, i.e. to the tetrahedral pyramids, consisting of equal equilateral triangles. In practice, this is unattainable, but it should be ensured that the grid does not contain elements with a highly distorted shape - the grid must be fairly regular. In addition, elements whose volume is zero or close to zero - the so-called degenerate elements - are unacceptable. The presence of even one degenerate element can make the calculation impossible. Therefore, the Master module contains functions for improving the overall quality of the mesh and correcting individual "bad" and degenerate elements. In practice, correction of the shape of an element can almost always be done only at the expense of degrading the quality of the elements surrounding it. Therefore, the algorithms for correction cannot be applied to a large number of elements at the same time. And vice versa - regularization algorithms are effective only when there are a large number of vertex elements (i.e. mesh nodes) which can be freely moved.

Mesh quality control and the correctness of the elements is maintained by a group of buttons on the Edit toolbar.


Buttons for the work with the mesh quality

A "bad" element is one that has too small and/or too large dihedral angles. Both mean that the shape of the element is far from correct, it is degenerated almost to a flat figure. Such an element will not be accepted for calculation in the Fourier and Hooke solvers. In addition, solvers impose a number of size requirements on finite elements and the ratio of volume and surface values. These requirements are dictated by the numerical calculation method, they cannot be changed.


To edit the limit values for the angles of elements that determine the quality of the mesh, click the Mark Parameters button on the Edit toolbar (for more details, see the Mark Parameters section).

The simplest and at the same time favor way to improve the quality of the mesh as a whole is the performing the Regularization function on the Edit toolbar.

With regularization, an attempt is made to bring the shape of all finite elements of the model closer to a regular tetrahedron. If the task of the functions of correcting "bad" elements (discussed below) is to avoid their presence even at the expense of some deterioration in the shape of neighboring elements, then the main task of regularization is to improve the quality of as many elements as possible without deteriorating the quality of others. In addition, during regularization, the nodes lying on the boundaries of the volumes, on the position of which its shape depends, remain stationary.

Regularization can be performed multiple times and the mesh quality can be improved. The function works automatically, without user intervention. You can evaluate the results of regularization by the quality histogram, which appears automatically upon completion of the operation (see the figure below).

The mesh quality histogram shows the percentage distribution of finite element dihedral angles over angle intervals. The width of the interval can be changed using the field located above the histogram from 1 to 60 degrees.


      

Mesh quality histogram with 10 degree and 1 degree angles interval

The filled columns show the distribution of dihedral angles after regularization, and the outlines show the distribution of dihedral angles before regularization. The legend is shown in the upper right corner of the diagram. In a good quality mesh, most of the dihedral angles should be close to 70 degrees in magnitude. The diagram is refreshed each time the regularization is performed.

The color of the columns of the chart is the same as the color of the element marks. To change the color of the diagram, use the Special Colors button on the Standard toolbar (for more details see Special Colors section).


The angle distribution histogram can be viewed by clicking the Mesh Quality button on the Edit toolbar.

Another way to improve the quality of the mesh is to reshape only the bad elements. In this case, the algorithm automatically moves the mesh nodes, trying to eliminate invalid geometric parameters in individual, selected mesh elements.


To find and select bad elements with angles outside the limit values, you should execute the Check Mesh Quality command on the Edit toolbar.

When the command is executed, it starts checking the angles and volumes of the casting and mold elements. Then the Check Mesh Quality dialog opens (see the figure below) and the check results are displayed in the table.


Check Mesh Quality dialog

The first column shows information about the mesh quality of the "casting" type bodies , the middle column shows information about the mesh quality of the "mold" type bodies, the third (right) column of the table shows summary information.

For each type of volume, the table shows the following information:

  • the number of elements with dihedral angles outside the range of permissible values;
  • minimum and maximum values of dihedral angles in elements;
  • how many elements with dihedral angles that go beyond the limits of permissible values are located on the outer boundaries of the volumes;
  • the number of elements with degenerate volumes;
  • how many degenerate elements are located on the outer boundaries of the volumes.

Found bad elements will be selected (marked) in accordance with the previously set label parameters (see the figure below).


An example of highlighting bad elements after checking the mesh quality

The correction of the marked bad elements is started by clicking on the Repair button in the Check Mesh Quality dialog.

The correction process is a responsible operation, since the movement of grid nodes applies to any "bad" elements of the model, including those located on the boundaries. This means that the operation may distort the shape of the object. Therefore, in the lower part of the Check Mesh Quality dialog, there are options for the algorithm for correcting elements (more precisely, shifting mesh nodes). By enabling and disabling the corresponding parameters, you can choose where to allow the displacement of nodes (in the casting, in the mold) and enable or disable the displacement of nodes at the boundaries. The last option allows you to keep the shape of the bodies without distortion.

It should be note that the maximum effect can be obtained when all types of fixes are allowed (enabled by default). The only drawback is the possible distortion of the boundaries of the casting or mold. It is better to run the simulation with some geometry inaccuracies than rebuilding the finite element mesh.

Correction of elements ends with a message (see the figure below). The information in the table of the Check Mesh Quality dialog will be updated.


Mesh repair completion message

If, after correcting the elements, the information in the table shows that not all elements with bad angles or degenerate volumes have been corrected, it is advisable to repeat the correction procedure again. It should be repeated as long as the number of detected "bad" elements continues to decrease. If after the next cycle the number of "bad" elements has not decreased, further attempts will be useless.

Elements that could not be repaired can be removed from the model. To do this, press the Delete button in the Check Mesh Quality dialog. After confirming the operation (see the figure below), the marked elements are deleted from the geometric model without the possibility of restoration.


Deleting elements operation confirmation

Thus, there will be no more bad elements in the model, but the model itself will slightly lose its geometric accuracy.

Attention! Removing bad elements from the model should be done only after the volume types have been specified and the correct bounding indices have been assigned.

When bad features are removed, the faces of adjacent features become external. However, these faces will have an index of 0 (inner face, see the figure below), which contradicts the accepted "standard" indexing of outer boundaries (usually these are indices 17 and 18). The standard numbering of borders is discussed in the Basic concepts: Volumes, Types, Indices section.


Boundaries with index 0 (gray) in places where bad elements are removed

If, after removing bad elements from the model, you change the type of any volume, you will be prompted to assign “standard” boundary indices to all volumes. This can have irreversible consequences. you cannot assign standard boundaries after deleting "bad" elements, so you should click the No button in the message box. If you click Yes, outer boundaries appear in the body of the casting or mold where the feature was removed, resulting in incorrect calculation of temperature fields in that area. It is important to preserve the outer boundaries with the index 0 in such places, this minimizes the calculation error, the system will assume the presence of a virtual element beyond such a boundary.

If, after removing "bad" features, gray faces appear on the outer boundaries of the casting or mold, you can optionally assign them indices of the boundaries with the environment using the tools described in the Editing Boundary Indices section.


Additional options for transforming finite elements are provided by the Split Elements command on the Editing toolbar. It becomes available after selecting any group of finite elements and opens the Split Elements dialog (see the figure below).


Split Elements dialog

In the window, you can select one of two schemes for the selected mesh elements.

  • Create nodes on element edges option divides each selected element into 8 smaller elements. The new nodes are located at the midpoints of the edges of the selection. This scheme is basic and cannot be turned off.
  • Create a node in the element center option expands the capabilities of the basic scheme. It adds another node in the center of each selected element. As a result, each selected element is divided into 12 elements. The option can be turned off. This option is enabled by default.

After pressing the OK button, the mesh refinement operation will be performed. An example of mesh refinement in the casting block feeders is shown in Fig. below.


a)                                                    b)                                                     c)
An example of local mesh refinement in the feeders of a casting:
a) original mesh; b) selection of elements in feeders; c) refined mesh

It should be note that the refinement, especially repeated many times, can lead to a deterioration in the quality of the mesh due to the appearance of elements with angles that go beyond the specified maximum and minimum values. In this regard, it is highly recommended to check the quality of the mesh after grinding it and correct the “bad” elements as described above.


Shells Creation


In investment casting, the melt is poured into a mold made by alternately wetting the master model with a binder and sprinkling it with bulk material. The result is a casting mold that has an arbitrary shape on the outside and an internal cavity that exactly repeats the geometry of the master model. Building such shells in a CAD system is a laborious process, so the Master module has a special tool for their automatic creation.

The shell is generated for all volumes displayed in the work area, i.e. before creating it, you must turn off the display of bodies in the model tree for which the shell will not be built. In addition, it may be necessary to set a special border index on surfaces where the shell is not needed. For more details, see the descriptions of the respective parameters below.


The Create Shell command on the Edit toolbar opens a dialog box where you can set the thickness and other parameters of the shell. The parameters of the finite element mesh of the shell can be selected automatically by the algorithm, or set by the user. By default, automatic selection of parameters is used (see the figure below).

Create Shell dialog

In automatic mode, the following shell parameters are set in the window:

  • Inlet Area Surface Index. The parameter specifies the index of the border where no shell is required. This is usually the surface through which the melt is fed. In the picture below, this is boundary 17. The parameter is enabled if there is a check mark in front of its name. In this case, the shell construction algorithm will create a separate volume above this area, which can then be removed (see the figure below). If the checkbox is off, the parameter is ignored and a closed envelope will be built for the entire volume.

  • Shell volume index. The parameter sets the volume index that will be assigned to the constructed shell. The index for which the required thermophysical properties of the material are specified in the file of the mold materials (for example, ceramics) are usually indicated here.

  • Volume Index Above Inlet Area. The parameter sets the volume index that is created above the boundary through which the form is filled with melt. It is only available when the Inlet Area Surface Index parameter is enabled. When the shell is created, its elements above the melt supply boundary will be allocated into a separate volume with a specified index, which can then be deleted to open access to mold filling and heat exchange with the environment (see the figure below).

  • Shell Thickness, mm. When creating real mold, their thickness differs in different areas, this is due to the topology of the master model. The algorithm simulates the result of a real process, therefore the parameter sets the average thickness of the shell. For example, it can be equal to the number of layers of bulk material, or measured on a test sample.

    a)                                                      b)                                                      c)
Constructio of the shell with the thickness 100 mm for the casting bush:
а) bush model prepared for the shell creation (boundaries indices mode);
b) the result of the shell construction (volume indices mode);
c) the model of the casting block in the form ready for simulating (boundaries indices mode)<

The shell generator is started by pressing the OK button. While it is running, a progress bar is displayed in the status bar. Upon completion of generation, the model will be loaded into the working area of the module.

The whole process of building the shell is recorded in the log files, in case of an error or in the event of a process failure, they can be viewed and sent to the developer. To do this, go to the folder for storing the log files through the Help → Show Logs menu, find the folder with the name containing the date and time the module was started and examine the contents of the files "shell_builder _ *. log" and "shell_builder _ *. err".

After construction the shell, it is necessary to remove the volume above the melt feeding boundary and perform the operation of assigning the standard boundaries (for more details, see the Editing Boundary Indices section).

After building the shell, its quality must be checked with the standard procedure of the Master module and, if necessary, correcting the bad elements. For more information about checking the quality of the mesh, see the Mesh Quality Control section. If the number of bad elements exceeds reasonable limits and cannot be corrected, it is recommended to delete the shell and repeat the operation of creating it, switching to the custom mode for controlling the shell mesh parameters.


The Settings button in the lower left corner of the Create Shell dialog opens the preprocessor settings window, where advanced parameters are located that control the shell creation algorithm (see the figure below).


Advanced shell creation options

These parameters control the operation of the mesh generator and affect the quality of the resulting shell. Most of them are optimally tuned and do not need to be changed unnecessarily. The exception is the Mesh Parameters group, where the Auto and User Defined switches are located (marked in the figure above). By default, the automatic selection of parameters is set. If you go to the User Defined mode and then click the OK button, the appearance of the Create Shell dialog box will change (see the figure below).


Create Shell dialog (User Defined mode)

In this mode, when creating a shell, the user can specify two more parameters. By default, both parameters are in the "auto" mode, which means that they are set by the algorithm automatically (which corresponds to the automatic mode). By setting the parameters to specific values, you can control the mesh quality and achieve the desired results.

  • Cell limit per side. This parameter adjusts the size of the shell mesh element based on the dimensions of the model. Smaller values of the parameter (from 30 to 50) give a coarser result, suitable for simple castings. Higher values of the parameter (80-120) are required when constructing a shell for shaped castings of complex shapes and allow a more accurate and high-quality mesh to be obtained. The minimum possible value of the parameter is 20. When entering a value of 0, the parameter switches to automatic mode.

  • Layer count. The parameter specifies the number of layers of elements in shell mesh. It is recommended to set this parameter less than the shell thickness. The minimum possible value of the parameter is 1. When you enter a value of 0, the parameter switches to automatic mode.

The procedure of shell creation can be repeated as many times as necessary to obtain multilayer shells with different properties. For example, in this way, first a model of a ceramic mold can be created, and then a model of thermal insulation wrapped around the mold. For more complex insulation, some layers of thermal insulation can be cut off by editing the volume indices (see details in Editing Element Indices in the Arbitrary Areas section).


a)                                      b)                                     c)                                     d)
Creation of mold and thermal insulation shells for investment casting modeling:
a) casting; b) ceramic shell; в) thermal insulation, 1 layer; г) thermal insulation, 2 layer

New Volumes Creation (Extrude Mesh)


In the Master module, you can create new volumes or modify the geometry of existing volumes by extruding a surface or part of it. This can be an extremely useful option when there is a need for a minor and relatively simple GM change. In this case, there is no need to change the original CAD model, then re-generate the FE mesh and import it into POLIGONSOFT. You can easily change the dimensions of the risers of a simple form, create layers of thermal insulation, covers, change the geometry of the casting, etc.

The created new volume is a mesh of the specified thickness, obtained by offsetting the selected surface. Therefore, before starting extrusion, you should mark the required surface or part of the surface (for methods of marking boundaries, see the Editing Boundary Indices section). As an example, let's increase the outer diameter of the pipe by 2 mm. In fig. below the outer cylindrical surface of the pipe is highlighted (the body has the "casting" type).


Selecting the outer cylindrical surface of a pipe

The Extrude Mesh button on the Edit toolbar calls a dialog box in which the thickness of the future shell is set in mm (see the figure below).


Extrude Mesh dialog

The window specifies the distance in mm by which the surface mesh should be extruded and the number of mesh layers that will be created in this case. By default, the second parameter is set to "auto", which means that the required number of layers will be calculated based on the size of the elements of the original mesh and the distance by which it will be extruded. The "auto" value can be changed to the desired one, then the layer thickness will be recalculated and the value will be shown in the information field.

After clicking the OK button, the operation will be performed, a new volume of the "mold" type will be created with a certain volume index that is different from the parent body (see the figure below).


Extrude Mesh Result

Now you can change the type, volume index and other parameters of the new volume, achieving the desired result. For example, following the purpose of the example, you can set the new body to a "casting" type and merge it with the pipe body, which will change the geometry of the casting.


Model Information


Model Information button opens a window that contains information about the loaded geometric model. To call the command, you can use the key combination Ctrl + I (see the figure below).


Information panel

The window displays the following information:

  • dimensions of the geometric model;
  • the number of nodes and elements in the casting and mold;
  • volumes of bodies included in the model, summed up by volume indices;
  • areas of boundaries, castings and shapes, summed by boundary indices.

The units in which the information is displayed are specified in the list in the upper left corner of the window (by default - meters). If necessary, the values of the table cells can be copied to the clipboard. To do this, you need to make two clicks on the cell to make it active, after which the copy function will become available.

You can close the information window by clicking the Close button in the lower right corner of the window.

Mold Filling

POLIGONSOFT can simulate the filling of a casting mold with melt and other flow-related problems. The calculation is performed by the Euler solver, which solves the problems of melt flow and its temperature change during interaction with the mold and the environment, including solidification. The following results can be obtained as a result of the calculation:

  • dynamics of the filling process of the mold with melt in time;
  • velocity field of the moving melt;
  • temperature field of the melt and the mold;
  • hydrostatic pressure field in the casting;
  • information about possible misruns (stopping the flow).

Additionally, Euler solver can be used to evaluate the gating system operation without considering the casting body.

The problem of predicting shrinkage cavities, macro- and microporosity is not solved at this stage.

To perform the calculation of mold filling, the following data are used, prepared in the Master module:

Preparing the Modeling Area


The Preparing the Modeling Area section describes the general principles of preparing the model for calculating the temperature field and shrinkage porosity. The case of mold filling simulation requires additional preparation of the geometric model in the Master module.


Melt Inlet


Usually, the mold is filled with melt from a ladle. To simulate this process on the casting surface, it is necessary to specify the area (or areas) through which the melt will flow. This can be done using the special tool Inlets.


The Inlets button on the Edit toolbar opens the editor (see the figure below), in which the parameters and the order of filling the casting are set. This data can be applied to solve one specific problem or exported to a bdv file for use in other calculations, including as a process template.


Inlets Editor

The window contains a table of regions (inlets) through which the melt is fed into the mold cavity (more precisely, into the casting body). Under the table there is an area for editing curves and parameters that control the stop of the solver.

To create a inlet, press the "+" button above the table on the left. The first row appears in the table. By pressing the "+" button you can create no more than 9 inlets. To delete an unnecessary inlet, select it in the table and press the "-" button above the table on the left, the line will be deleted from the table.

After the inlet is created, it needs to be configured. Since the inlet is a virtual stream of the melt, you need to set its dimensions and the entry point on the casting surface. It is assumed that the jet is in the shape of a cylinder, so its size is given by the diameter of the circle. Left-clicking with the Ctrl key pressed on the casting body starts drawing a circle, which denotes the melt stream (see the figure below). The center of the circle is set by the mouse pointer.


The melt inlet specified at the upper end of the downsprue

To the right of the table is the Snap to mesh node option. When this parameter is enabled, the center of the stream is always at the mesh node. When the parameter is deactivated, the center of the stream is displayed at random. It is recommended to enable this parameter if the melt inlet is specified at the symmetry boundary, this will avoid errors when calculating the area of the stream.

In the inlet line, you can set the parameters that determine the operation of the inlet (see the figure below).


Inlet parameters

  • Index. This is the number of the inlet. By default, the indices are assigned to the inlet in ascending order, but the user can independently assign any index from 1 to 9 to the inlet. To do this, double-click the left mouse button in the cell and enter a new value. You should be aware that multiple inlet with the same index cannot be specified. The checkbutton next to the color indicator shows whether the inlet will be used in the calculation or not. Thus, several inlets can be set, but only those that are checked will be used in the calculation.


  • Type. The parameter sets the Inlet operation mode. It is set by a value from the drop-down list:


    • Inlet. This is the default lnlet type, it means that through the specified lnlet, the melt enters the casting volume (i.e., into the mold cavity).


    • Outlet. This type of Inlet sets the condition for the outflow of the melt from the casting volume (i.e., the mold cavity). Sometimes the task of checking the operation of the gating-feeding system (GFS) is considered separately from the task of filling the casting. In this case, the body of the casting is the GFS itself, in which there are places for supplying the melt and places from where the melt flows out into the intended casting. In this case, it is necessary to specify the regions of the melt inlet and the regions of its outflow from the GFS.


  • Diameter, mm. Specifies the stream diameter in millimeters. To edit, you need to double-click the left mouse button in the cell and enter a new value.


  • Velocity Mode. From the drop-down list, you can choose how to set the velocity of the melt entering through the inlet: Value or Curve. If Value is selected, the speed is given by a constant in the Velocity ​​cell. This means that the melt velocity will be constant throughout the calculation. If the Curve option is selected, a table and a curve drawing area become active at the bottom of the editor, where you can set the curve points in the velocity-time coordinates. If the calculation time exceeds the time set for the last point of the curve, the velocity remains the same as the velocity set for the last point of the curve. For more information about editing curves, see the Editing Constants and Functions section.


  • Turn On. The parameter sets the moment when the Inlet starts. Double-clicking in a cell shows a drop-down list where you can select one of two conditions (see the figure below):


    • Time. Inlet starts to work at a certain moment of time, which is set in seconds from the beginning of the process in the Time, s field. If one Inlet is used, then the beginning of its operation, as a rule, coincides with the beginning of the calculation, i.e. starts at 0 second.


    • Other Inlet. The Inlet starts working when another Inlet with the number specified in the Inlet Index field finishes. For example, if you need to simulate adding melt to the riser, then on the riser surface, inlet # 2 is set, which starts working after the main inlet # 1, set on the surface of the downsprue, finishes working. If desired, in the Pause, s field, you can set a pause in seconds between the end of the first inlet and the start of the inlet approach. Such a pause may correspond to the movement of the ladle from the downsprue to the riser.

    

Inlet turn on conditions

  • Turn Off. The parameter sets the moment when the Inlet operation stops. Double-clicking in a cell shows a drop-down list where you can select one of two conditions (see the figure below):


    • Time. The Inlet finishes working when a certain moment of time of the filling process is reached, which is set in seconds in the Time, s field.


    • Percentage of filled. The Inlet finishes working when the total volume of all volumes of the Casting type is filled to the level specified in the Filled,% field. If one Inlet is used, then the end of its work, as a rule, is set by the condition of 100% filling of the casting.

    

Inlet turn off conditions

After the Inlets are set, press the OK button. The dialog box will be closed, but the Inlets outlines will continue to be displayed on the model.

As an example of the use of several inlets, consider the technology of filling a form through a riser with topping up in two risers (see the figure below).

The main inlet # 1 supplies the melt to the downsprue from the beginning of the calculation (from the zero second of the process). When the form is half full (the condition for turning off the feed - 50% of the filled volume), inlet # 1 is turned off, followed by a pause of 2 seconds, during which the ladle moves towards the inlet # 2, set on the surface of the first riser. In 2 seconds after the end of the first inlet , inlet # 2 starts to work - the melt is added to the first riser. It ends when the form is 3/4 (75%) full, then a pause of 2 seconds follows and inlet # 3, set on the surface of the second riser, begins to work. Inlet # 3 finishes work when the form is completely filled out.


Example of using multiple inlets

For each inlet, the velocity of the incoming melt is set. The velocity can be set by a constant value or by a speed-time curve.

In most cases, the flow calculation is stopped (followed by switching to the solidification and porosity calculation) when the mold is completely filled or when all Inlets are turned off. This situation is automatically tracked by the preprocessor and the solver and no additional parameters need to be specified. In cases where the moment of stopping the flow calculation cannot be determined automatically, it is required to specify it manually in the Stop Flow Simulation group. See the section below for more information.

Prefilled Volumes


In some problems, at the beginning of the calculation, there are volumes filled with melt. For example, it can be a pre-filled sprue bowl from which the melt is poured into a mold or crucible of a melting plant during vacuum casting. In such a situation, it can be inconvenient to set the flow of the melt into the mold using the inlet. To simulate the melt flow, an appropriate geometric model must be prepared, in which the pre-filled volume or volumes will be present.

If the container, for example, the ladle at the beginning of the calculation is not completely filled, its model should be prepared, consisting of two parts - filled and empty, or only the filled part of the bucket should be used in the calculation (see the figure below). In any case, the pre-filled volume must be of the "casting" type and have a separate index (in the figure below it is the volume with index 9, red).

The selected volume index should be marked as “filled” in the model tree. To do this, left-click on the empty glass icon to the right of the bucket volume (in this example, this is index 9). In this case, the icon will change to the image of a full glass. To undo the changes, you need to click on the icon again. If the model contains multiple casting volumes with the selected index, all of these volumes will be marked as filled, so it is important to set a unique index to the filled volume. In general, you can specify several filled volumes, and they can have different indices, but you cannot specify as filled one volume from a group of volumes with the same index.


a)                                                 b)                                                    c)         
Ways to set the pre-filled volume:
a) complete geometry; b) only the filled part; c) Model Tree

If a prefilled volume is specified in the model, this is a sufficient condition to run the Euler flow solver; you do not need to specify melt inlets. However, if the problem statement requires it, leads can be specified as usual in the required quantity as described in the previous section.

In most cases, the flow calculation is stopped (followed by switching to the solidification and porosity calculation) when the mold is completely filled or when all the inlets are turned off. This situation is automatically tracked by the preprocessor and the solver and no additional parameters need to be specified. In cases when the moment of stopping the flow calculation cannot be determined automatically, it is required to specify it manually in the Stop Flow Simulation group. See the section below for more on this.


Solver Stop Parameters


At a certain moment in time, the Euler flow solver stops working and further calculations are performed by other solvers (if specified in the simulation parameters). In most cases, this happens when all inlets finish their work, for example, when the mold cavity is completely filled with melt. There are situations when the desired moment of switching between solvers does not coincide with filling the mold, here are some of them:

  • Centrifugal casting, where the melt enters a container from which it is forced out by centrifugal force into molds. In this situation, the cast volumes are never completely filled.
  • Large ingots casting. Switching to the solidification calculation immediately after filling the ingot turns off the convection of the melt, which leads to the "freezing" of the velocity and the temperature fields. As a result, the calculation of porosity may not be accurate. It takes some time after the end of filling the mold for the velocity field in the ingot to come to a more or less stationary state.

In order to flexibly control the moment of stopping the flow solver and switching to the calculation of solidification and porosity (if specified) in the Inlets dialog there is a Stop Flow Simulation group, where you can select one of two criteria that stops the work of the Euler solver (see the figure below).


Stop Flow Simulation Group in the Inlets Dialog

Consider the criteria for stopping the Euler flow solver.

  • Inlet Stopping. When creating an inlet or several inlets that stop working at a given percentage of filling the casting, the Master by default includes the criterion for stopping the Euler solver when the last lead is turned off. This will be the inlet with the maximum fill percentage is set in the Turn Off Condition. In the example with three leads described in the section Melt Inlet this is inlet number 3, it has the criterion for stopping 100% filling of the casting. This is the simplest and most common case when modeling a flow; in this case, you do not need to perform any special actions to adjust the parameters. If there are several inlets and some of them stop working by the percentage of filling, and some by time, a situation arises when the solver cannot automatically decide the moment of stopping the flow calculation. In this case, editing of the approach index in the parameter line becomes available, and disabling the specified approach turns off the flow solver. Additionally, you can set a pause during which the solver will continue to count the flow of the melt after turning off the specified inlet.
  • Time. If this parameter is selected, the calculation of the flow will run over time regardless of the mode and disconnection of the inlets. This parameter must be set if the casting is filled by draining the melt from the prefilled volume (see section above). If this parameter is selected, the calculation of the flow will run over time regardless of the mode and disconnection of the inlets. This parameter must be set if the casting is filled by draining the melt from the prefilled volume (see Section above). Also, this parameter must be set, for example, for vacuum centrifugal casting, when the melt is poured into the distribution volume, and then flow into the mold under the centrifugal force.

Solidification and Porosity

Calculation of casting solidification taking into account heat exchange with mold elements and the environment, as well as pipe shrinkage and macro- and microporosity prediction - the main task of the casting process simulation, performed in the Fourier solver. As a result of the calculation, the following results can be obtained:

  • casting and mold temperature fields;
  • hydrostatic pressure field in the casting;
  • casting phase field;
  • shrinkage pipes in the casting;
  • macroporosity and microporosity in the casting.

Note. At the beginning of the solidification calculation, the mold is filled with melt. The filling problem is solved by the Euler module (see the corresponding section).

To perform the majority of temperature fields, solidification and shrinkage porosity calculations, in addition to the prepared geometric model, the following data are required, which are need to be prepared in the Master module:

  • initial casting temperature;
  • initial mold temperature;
  • ambient temperature;
  • simulation stop criterion;
  • alloy properties (see the Materials Assignment section);
  • mold thermal properties (see the Materials Assignment section);
  • heat transfer parameters (see the Heat Transfer Parameters section);
  • other conditions depending on the process being modeled.

A detailed description of all parameters is given in chapter 6 in Thermal and Porosity Simulation section.

To simulate special casting processes (vacuum casting, continuous casting, high and low pressure die casting, etc.) additional data may required (see relevant sections).

Stresses and Deformations

POLIGONSOFT allows modeling the stress-strain state (SSS) of a casting and a mold, taking into account their contact interaction. The calculation of stresses and strains is performed by the Hooke solver based on the results of a previously made thermal calculation using the Fourier solver. The following results can be obtained as a result of the calculation:

  • casting and mold effective stress fields;
  • casting and mold effective strain fields;
  • displacement field in the casting (warping);
  • failure indicator field in the casting;
  • casting and mold effective plastic strain fields;
  • gap field between the casting and mold;
  • hot tearing indicator field;
  • principal stress fields in the casting;
  • etc.

As a result of the work of the Fourier thermal solver, files are written that contain the temperature fields of the casting and the mold. These files are part of the initial data for the work of the Hooke solver.

To calculate SSS, you need the following data, which are prepared in Master and Mirage (optional) modules:


Preparing the Modeling Area


Preparing the Modeling Area section describes the general principles of preparing a model for calculating the temperature and porosity field. The task of modeling the SSS requires additional preparation of the geometric model in the Master module.

First of all, thermomechanical properties of materials must be specified for all volumes involved in the stress calculation. They are specified in the Material Properties editor on the Stress tab. The set of thermomechanical characteristics describes the behavior of the material when its temperature field changes and is called the Thermomechanical Model. Two thermomechanical models are available for the casting material, they are selected from the list in the first row of the table (see the figure below):

  • Linear Elastic;
  • Elasto-plastic.

Material Properties editor, Stress tab.
Elasto-plastic model

For most cases, the Elasto-plastic model is suitable; it requires specifying the entire set of parameters. In some cases, it is permissible to use a Linear Elastic model of material behavior for casting (see the figure below). This is a simpler model, it is suitable for cases when the body is always in a solid state and its temperature changes insignificantly and the heating-cooling cycle does not lead to an irreversible change in shape and size.


Linear Elastic model

In the current version of POLIGONSOFT, only one thermomechanical model can be used for a mold, regardless of the number of materials used in it. The behavior of the mold is set by the calculation parameter Contact Interface with Mold in the Start Simulation and Parameters window on the Stress tab (for more details, see the Stress Simulation Parameters section). There are several options available for the interaction of the casting with the mold:

  • Not simulated;
  • Rigid Mold;
  • Linear Elastic Mold;
  • Elasto-plastic Mold.

The first two models do not require any thermomechanical properties. This means that even if such properties are specified, they will be ignored. The other two models require that thermomechanical properties be specified for the mold. When selecting the type of contact interaction with the mold in the Start Simulation and Parameters window, you will also need to set the Mold Stress Properties Volume Index parameter there. The parameter specifies the volume index of the mold to which the material properties are tied and its thermomechanical properties will be used when calculating stresses. Thus, when selecting the interaction of a casting with a linear-elastic or elastic-plastic mold, thermomechanical properties are used only from one material specified in the Mold Stress Properties Volume Index parameter. In this case, it is not necessary to switch the thermomechanical model in the material properties; it is enough that all the necessary properties are specified.

Some mold materials can be excluded from the calculation. For example, if the interaction of a casting with an absolutely rigid or linear-elastic metal or sand mold is calculated, then the mechanical impact of such materials as fibrous refractories, wools, burnt-out exothermic sleeves can be ignored. Such mold materials can (and should!) be excluded from the SSS calculation, while their role in the formation of the thermal picture is preserved. To do this, check the box Empty for Stress for the selected volume indices in the Volume Parameters window, or click the mouse in the model tree in the corresponding cell opposite the desired index (see the figure below). Volumes not involved in the SSS calculation are considered "empty" from the point of view of the thermomechanical problem and are marked with an icon in the model tree.

Note: Only volumes of the Mold type can be marked empty.


Assigning the Empty for Stress option
in the Model Tree (left) and in the Volume Parameters window (right)


Removal of Gating System Elements


The stress-strain state of the casting at the end of solidification is generally determined by the history of the temperature field of the casting and its contact with the mold. This state is not final and may differ significantly from the stress-strain state of the casting removed from the mold and freed from the elements of the gating-feed system (GFS). Thus, in order to understand what the final geometry of the casting will be and how the residual stresses will be distributed in it, it is necessary to calculate the stress-strain state after removing the casting block from the mold and then after removing the GFS elements.

To calculate the change in SSS when removing collectors, risers, etc., these elements of the model must have unique volume indices. It is best to split the model into separate bodies at the stage of its preparation in the CAD system, but in extreme cases, you can split the casting mesh into parts using the Master tools described in the Editing Element Indices in the Arbitrary Areas. section. As a result, the casting will consist of bodies with different indices, while index 1 is always assigned to the casting itself, i.e. the part that remains and cannot be deleted (see the figure below).


Casting prepared for SSS calculation after GFS elements are removed. Volume indices: 1 (white) – casting (part); 4 (cyan) – riser; 5 (yellow) – gating system

The moment of removing the casting from the mold and the order of removing the elements of the gating system are specified in the calculation parameters (for more details, see the sections Stress Simulation Parameters and Stress Simulation at subsequent technological stages).

Radiation Heat Transfer

Radiation heat transfer can be considered on the basis of a simple expression for the heat flux (Stefan-Boltzmann law) or calculated using a complex algorithm that takes into account the effects of radiation energy reflection and shading of some bodies in the computational domain by others. Such a model for calculating radiation heat transfer used in POLIGONSOFT does not take into account the influence of the environment, i.e. the heat transfer coefficient to the environment and its temperature (including the dynamic medium) are not taken into account.

A problem involving complex radiation heat transfer requires special preparation of the geometric model and, in some cases, special process conditions, for example, the translation mode that sets the movement of bodies when modeling directional crystallization (for more details, see the Translation Modes section).

Since radiation heat transfer involves the calculation of the interaction of parts of the geometric model (casting and mold) with the surrounding bodies (for example, a melting unit, a mold heating furnace,...), it is necessary to include these bodies into a geometric model. You should remember that it is recommended to create a closed computational domain. For example, if the casting is placed in a vacuum melting plant, the computational domain will include: a casting, a mold and a vacuum melting plant body surrounding them.


a)                                           b)

A simple example of a computational domain for calculating radiation heat transfer:
a) geometric pattern. Different colors are shown: white - casting, purple - mold, yellow - imitates the surrounding body that closes the computational domain; b) thermal fields showing the influence of the "adjacent" casting

In the above example, the entire geometric model is included in the computational domain. However, it can only include part of the geometric model. Consider a more general (and more complex) case: a mold with a melt is placed in a mold heating furnace. Part of the form goes beyond its limits. It can be assumed that for the part of the mold that is in the furnace, the influence of the environment (air) is insignificant. Most of the heat is transferred by radiation. On the other hand, this cannot be said for the part of the form that goes beyond its boundaries. The task can become more difficult if the furnace has a movable bottom. Then, when the casting moves in the vertical direction, the area located in the zone of radiation heat transfer changes in time. Consider two options for configuring the computational domain.

Option 1. The geometric model (file g3d), in addition to the casting and the mold, includes a model of the mold heating furnace and the enclosing body – Enclosure, which simulates the walls of the vacuum melting plant or other surrounding space that closes the computational domain (see figure below). In this case, the boundary of the radiation heat transfer computational domain by is set according to the Enclosure. If the mold heating furnace and the casting move relative to each other, a curve (or curves) of their movement is specified (for more details, see the Movement Modes section).


a)                                           b)

Computational domain (option 1): a) geometric model; b) the calculated area of radiation heat transfer (denoted by a parallelepiped). Different colors are shown: white - casting, purple - mold, yellow - mold heating furnace, green - melting unit or a body that closes the computational domain

Note. The shape and size of the Enclosure can significantly affect the accuracy of the simulation results. For example, when modeling radiation heat transfer in a vacuum melting plant, it is important to accurately convey the shape and dimensions of its internal space.

Option 2. In the geometric model (.g3d file), in addition to the casting and the mold, only the model of the mold heating furnace is included. In this case, the boundary of the radiation heat transfer region is set by the heating furnace (see the figure below). If the mold heating furnace and the casting move relative to each other, a curve (or curves) of their movement is specified (for more details, see the Movement Modes section). The difference from the first option is that for those parts of the mold and casting that go beyond the boundary of the computational domain (over the heating furnace border), the calculation will be carried out according to the usual heat exchange rules with the environment. The coefficient of heat transfer to the environment, emissivity and the environment temperature, specified in the heat transfer parameters, will be taken into account (for more details, see Heat Transfer Parameters section). It should be note that with this version of the problem setup, the part of the casting block located outside the ray tracing boundaries (outside the mold heating furnace) does not participate in radiative heat transfer with energy reflection effects, which can be a rather serious deviation from the real situation. For example, in the described example, the part of the mold that is outside the heating furnace will not “feel” its (furnace) influence.


а)                                           б)

Computational domain (option 2): а) geometric model; б) computational domain (indicated by a parallelepiped). Different colors are shown: white - casting, purple - mold, yellow - mold heating furnace

The first option corresponds to a more accurate formulation of the problem, since in reality the radiation of the heating furnace also affects the parts of the block protruding from it. The second option can be used when it is necessary to take into account the exit of the mold from the radiative heat exchange zone with a change in the heat transfer model, for example, when the mold moves in a vacuum melting unit and is immersed in a liquid metal cooler or other medium. In this case, the accuracy of the calculation may suffer due to the creation of a clear virtual boundary that the radiation cannot "surpass". (see the figure below).


a)                   b)

Simulatiuon results (temporary) when the computational area is defined: a) option 1; b) option 2

Let us give an example of an incorrectly specified computational domain. In fig. below, the dimensions of the computational domain are given according to the casting mold. As a result, an open computational domain was obtained - the upper end of the bowl does not have radiation heat exchange with any other surfaces. Due to the specified dimensions of the radiation area, the possibility of heat transfer from the top of the bowl to the environment is excluded. If you carry out a calculation with such a calculated area, then the results obtained will be absolutely incorrect (cooling will go from bottom to top, the metal in the bowl will solidified last).


An example of an incorrectly defined computational domain

To set the computational domain for radiative heat transfer, you need to specify a body that plays the role of Enclosure. To do this check the Enclosure checkbox (1) in the Volume Parameters dialog of the selected body (see the figure below). In this case, an interactive icon (2) will appear in the Model Tree opposite the marked body. You can control the enable/disable Enclosure by right-clicking on the icon in the Model Tree.


Assignment of Enclosure to the body of the melting plant (volume index 3)

It should be note that Enclosure, like other volume parameters, is assigned to the volume index, so all other bodies with the same index will also be marked as Enclosure.

The volume indices for radiation heat transfer are assigned in accordance with the idea of the calculation. For example, if you want to calculate the cooling of a casting block installed in a metal box and the temperature distribution in it is of interest, you should assign a volume index according to the required material properties (see the Materials Assignment section). If the calculation of the temperature in the Enclosure is not of interest, you can assign a constant or variable temperature to it (for more details see the Temperature Modes section).

Calculation of radiation heat transfer (ray tracing) is performed only between the boundaries with the "default" environment, i.e. between boundaries 17 and 18 in any combination. For each of these boundaries it is necessary to set the emissivity E in the heat transfer file and the parameter Surface Fraction > 0 (for more details see the Heat Transfer Parameters section). If the computational domain includes many bodies with different emissivity: casting mold, sand, thermal insulation, steel flask, etc.), in this case you can set the values ​​of the effective total emissivity, depending on the emissivity of each surface that falls within the raytraced region for boundaries 17 and 18 in the heat transfer file. All types of boundary heat transfer can be taken into account in the heat transfer settings. This will make it possible to solve "mixed" problems that have areas where heat transfer by radiation is calculated taking into account the effects of shading (there is no environment), and areas in which the calculation is carried out taking into account the environment.

If other external boundaries except 17 and 18 are assigned in the radiative heat transfer computation domain, they do not participate in the complex radiative heat transfer, i.e. for them, heat exchange with the medium can be set, including taking into account radiation according to the Stefan-Boltzmann law. This can be used, for example, if it is important to take into account the heating of the walls of the vacuum melting plant from the inside and their heat exchange with the environment from the outside. On the fig. below on the outer walls of the melting plant (given by the Enclosure for the area of radiative heat exchange), boundary 3 is set for which convective heat exchange with the medium is set.


An example of assigning boundaries for the case when the calculation of temperatures in the Enclosure is performed

To reduce the calculation time, it is reasonable to assign the outer boundaries of the Enclosure body index 0 if the calculation of temperatures in the Enclosure is not performed (see the figure below). You can also assign boundaries with index 0 to some surfaces of the computational domain that do not play a significant role in heat transfer.


An example of a GM with "zero" external boundaty indices

Before starting the calculation, you should enable Radiation parameter on the Solidification tab in the Start Simulation and Parameters window (for more details, see the Thermal Simulation Parameters section). This parameter launches a special module that is used to calculate ray tracing between objects, incl. moving. The module uses as source information a .g3d geometry file, in which at least one volume index is set as Enclosure, a .mve movement modes file (if specified) and creates a .rad file based on them, which contains not only the FE mesh, but also information about which surfaces (triangular facets) "see" each other and which "do not see" in each specified moment in time (the so-called view factors). If the volumes of the model do not move, the view factors do not change during the calculation. However, if the bodies are assigned movement modes, the view factors must be calculated for different positions of the bodies relative to each other. Further, the .rad file is used by the Fourier solver.

Exothermic Sleeves

The solution to the problem of saving metal while maintaining or increasing the efficiency of casting feeding can be achieved by using mold materials that insulate or give off heat to the riser zone. When modeling the conditions for manufacturing a casting in a mold containing special heat-insulating materials, to take into account their effect, it is enough to set the necessary thermophysical properties in the corresponding volume of the 3D mesh model. The case of modeling the casting formation in the presence of materials from an exothermic (thermite) mixture requires the use of a special model, which includes a description of the conditions for starting a reaction during thermal interaction of a material with a liquid metal, its dynamics with active heat release.

The decision on the need to use exothermic materials is made when developing or optimizing a foundry technology, based on considerations of economy and the possibility of ensuring the required quality of castings. In a number of cases, the installation, in particular, of exothermic sleeves of the risers may be the only available option to ensure continuity in the feeding zone of the casting due to the structurally determined insufficient size of the platform for placing the risers. When initially choosing the size of the riser sleeve (see figure below), the concept of its module (equivalent to ordinary riser) is used to select a specific model from the proposed size range. Subsequent adjustment of the decision made can be carried out based on the results of the analysis of its effectiveness in POLIGONSOFT.


Making identical castings using risers with comparable feed efficiency:
1 – heated reduced size; 2 – traditional type

Riser sleeves can be supplied complete with sandy breaker cores attached to the bottom of the sleeve to form a thinned section of metal - the riser neck. This reduces the labor intensity of separating the riser from the casting and reduces the amount of machining in the corresponding area. With the correct choice of the width of the riser neck, formed by the central channel of the breaker core, the liquid metal feeding does not stop during the entire time of operation of the riser. In this case, premature solidification of the riser neck does not occur due to its strong heating and low heat removal from the breaker core.

In most cases, exothermic sleeves are made from a bonded mixture of powders or shavings of aluminum and iron oxide with additives. In contact with the melt, such a material ignites when relatively low temperatures of several hundred degrees are reached. The oxidation reaction is accompanied by intense heat release (Q): Fe2O3+2Al®Al2O3+2Fe+Q. For a relatively short time of the reaction in the sleeve volume, the released heat makes it possible to heat the riser, create a strongly heated zone around it, including the area filled with reaction products, and the adjacent mold. In most cases, the sleeve retains its coherence after completion of the reaction, which allows it to work already as a loose, effectively insulating material, and also facilitates measures to prevent clogging of the recycled mixture.


Exothermic Mixture Model


To simulate the exothermic sleeves behavior during casting solidification, POLIGONSOFT implemented a common adapted model, when the sleeve material has a set of effective thermophysical properties, including the heat of transformation. When using this approach, it is possible to simulate with sufficient accuracy the conditions of thermal behavior of the material without calculating the reaction of interaction between the components of the exothermic mixture. In this case, the effective thermal properties of the material are tied to the model and size of a particular sleeve or their family and include a “correction” associated with phenomena that cannot be taken into account when solving the Fourier equation (1).


(1)

where t – temperature;
          t – time;
          C(t), r(t), l(t) – temperature-dependent values of specific heat,
                density and thermal conductivity of the exothermic mixture, respectively;
          I – heat source function, written in the following form:



(2)

where m – converted substance fraction, which varies from 0 to 1 for a fixed time for the
                 complete reaction tign;
          L – thermal effect of transformation.


Initially I=0. The moment of activation of the heat source in the sleeve mesh node is determined by the ignition temperature tign, which starts the counter of the increment of the current reaction time τr (from 0 to τign) and the fraction of the converted substance. Heat release stops at m=1, the mixture in this area is considered "burnt out" (3).


(3)


The initial data for solving equations (1) and (2) are the given properties of the material C(t), r(t), l(t), L, and m=m(tr); m set on the interval [0, 1] as a function dependent on the reaction time tr on the interval [0,tign].

Type III boundary conditions apply at the contact between the sleeve and the metal, the mold, and the environment. The heat transfer coefficients along the boundaries with the metal and the mold can also be included in a set of individually specified special data due to the special characteristics of the sleeve surface.


Preparing the Modeling Area


To calculate the foundry technology using exothermic materials, the model prepared in the Master module must contain their volumes (see the figure below), which are assigned indices in accordance with the specified materials or the loaded process template (mold materials properties file). For more information about indices and materials, see the Basic concepts: Volumes, Types, Indices и Materials Assignment sections .


Casting model with exothermic sleeves
(volume index 7, dark green)

In the example shown in Fig. above, the volumes of sleeves are assigned volume index 7, for which the properties of the exothermic material should be specified. For such a material, in addition to the usual set of properties specified on the Thermal tab in the Material Properties editor, you should enter parameters that specify the conditions for the thermite reaction to occur (see the figure below). They are set on the Exothermic tab (for more details, see the Material Properties section and the corresponding subsection).


Thermal properties of exothermic material,
given for volumes with index 7


Exothermic material properties,
given for volumes with index 7

The exothermic materials properties can be set by the user based on their own data or loaded from the database, from the %POLIGONSOFT%\Database\Exothermic folder.

By default, the heat transfer between the casting and the exothermic sleeve is set in the same way as between the casting and the sand mold, i.e. the boundary with index 1 is assigned. There are no boundaries between the sleeve and the mold by default, an equivalent mesh is used, which corresponds to the conditions of ideal contact. To set heat transfer conditions at the casting-sleeve boundary other than those accepted by default, select the corresponding interfaces and set a unique index for them, for which the required conditions must be entered in the heat transfer parameters (for more details, see the Editing Boundary Indices and Heat Transfer Parameters sections). The figure below an example of assigning special indexes on the casting-sleeve and sleeve-mold surfaces is shown.


An example of a model with special heat transfer conditions at the boundaries
casting-sleeve (index 2) and sleeve-mold (index 3)

For the boundaries selected in this way, the appropriate heat transfer conditions should be assigned (see the example in the figure below) in the Heat Transfer editor. Parameters can be set independently or imported from templates stored in the %POLIGONSOFT%\Database\Exothermic folder.


Heat transfer parameters for assigned boundaries
casting-sleeve (index 2) and sleeve-mold (index 3)

It should be note that in order to set the sleeve-mold interface, it is necessary to separate the sleeve volume index from the mold and set the contact condition between the parts of the mold (for more details, see the Editing Volume Types and Indices). In most cases, you can skip this and work with ideal contact between mold elements.

After entering other parameters of the casting process, set in a standard way, and performing calculations, it becomes possible to observe the thermal history of exothermic sleeves and analyze the effectiveness of their use to provide improved casting feeding (see figure below).


Temperature-phase fields in the casting and the mold after 30 minutes from the start of solidification:
using exothermic sleeves (left); without sleeves (right)

Virtual Mold

POLIGONSOFT allows you to simulate thermal processes in a mold without its finite element (FE) model using a virtual semi-infinite mold (or virtual mold ). This can be useful, for example, when simulating large sand casting or investment casting using a casting box with sand. The virtual mold can also be used for all types of casting, if the main purpose of modeling is to fill the mold with melt.

When using a virtual mold, it is necessary to specify the properties of its material and the heat transfer coefficient at the casting-virtual mold or mold-virtual mold interface.

Since the volume (FE mesh) of the mold is absent in the model, its thermophysical properties cannot be specified as usual - using the volume index. Instead, the boundary index specified on the outer surfaces of the casting or mold is used, behind which the virtual volume is supposed to be located. In accordance with the rules for using boundary indices (see Basic concepts: Volumes, Types, indices section for more details) indices from 1 to 9 can be assigned to the boundaries between a casting and a mold or a mold and a mold . If such an index is assigned to the outer boundary of the volume, i.e. a surface beyond which there is no other volume, the POLIGONSOFT solvers will assume that a virtual volume of the Mold type is located beyond this boundary. The material properties of the virtual mold will be used according to the boundary index. This rule is true for both casting boundaries and mold boundaries - in any case, the virtual volume will a Mold type. Below are two examples of using a virtual mold.

Example 1. Sand casting. The figure below shows an example of a sand casting in which the volume of the mold is replaced with a virtual one. The full case is shown next to it for comparison. The complete model (right side of the figure) consists of a casting, a mold and a chill plate, a model with a virtual model (the left side of the figure) consists of a casting and a chill plate. In the full case, the volume of the mold with index 3 is given, for which the properties of the sand are set, and the heat transfer parameters are set by the standard boundary with the index 1. In the problem with the virtual mold, sand properties and the heat transfer parameters are set by the boundary with the index 3. In this case, you need to make sure that the heat transfer for boundary 3 is set correctly.


a)
b)
Volume Indices: 3 – Sand Mold; 5 – Iron Chill

a)
b)
Boundary Indices: 1 - Casting-Mold Interface; 3 - Casting-Virtual Mold Interface
Case setup for sand casting: a) a case with a real shape; b) a case with a virtual mold

Example 2. Investment casting. In the manufacture of steel castings by the Investment Casting method, the shell mold is usually placed in a box with a sand. Geometry design of a shell mold in a CAD system is a rather laborious task, therefore the preprocessor has a special tool that allows you to create a realistic mesh model of a shell mold (see details in Shells Creation section). However, in this case, the problem of joining the mesh of the sand, which can be designed in the CAD with the shell mesh, which is built without reference to geometric objects. The figure below shows hot to solve this problem using a virtual mold.

In variant a, it is assumed to use completely virtual sand through the designation of a border with an index of 6 on the outer surface of the form (border 6 is selected in accordance with the attached Mold Properties File). The heat transfer coefficient for this boundary (shall-sand) is specified in the heat transfer file.

In variant b, some layer of sand was created using the Create Shell tool. In this case, the thickness of the shell with the sand properties (volume index 6 in accordance with the attached Mold Properties File) is chosen arbitrarily, for example, it can be equal to the distance from the ceramic mold to the bottom of the box, or you can focus on the heating thickness of the sand. In any case, the second shell simulates the heating of real sand, and on its outer boundary heat transfer with virtual sand is specified, as in variant a.


a)
b)
Boundary Indices: 2 – Ceramic Shell; 6 – Silica Sand

a)
b)
Boundary Indices: 6 – Mold-Virtual Mold Interface
Case setup for investment casting: a) virual sand case; b) semi-virual sand case

Centrifugal Casting


Introduction


The principle of centrifugal casting is that the mold is filled with melt and the casting is formed by rotating the mold around a horizontal, vertical or inclined axis. This achieves additional impact on the melt and the solidifying casting by the field of centrifugal forces. The process is implemented on special centrifugal machines.

The main feature of the castings formation during centrifugal casting is that the filling of the mold with melt and the solidification of the casting occur under the action of centrifugal forces, which, as a rule, are many times greater than the force of gravity.

There are many different types of centrifugal casting units. In general, centrifugal casting allows for denser castings. In gravity casting, the acceleration vectors of gravity are always directed parallel to each other in one direction. In centrifugal casting, the acceleration vectors are directed from the rotation axis in all directions. This allows for the organization of such formation of thermal hot spots that will more effectively remove shrinkage defects from the body of the part. In addition, since centrifugal forces are significantly greater than gravitational forces, centrifugal casting improves the filling of the mold, provides increased piezometric pressure in the body of the solidifying casting, and a finer structure on the outer surface due to a better contact and more intensive cooling along the outer surface.

Limitations of the Model


The centrifugal casting model used in POLIGONSOFT has the following limitations:

  • When modelling the centrifugal casting you should orient the bodies in the Master module in such a way that the rotation axis direction coincides with the Y axis of coordinates.
  • The inlet should be set so that its center is located on the axis of rotation (or close to it).
  • When modelling solidification and porosity, it is believed that centrifugal force is far superior to gravity. Such permission allows us to neglect the effect of gravity on the course of solidification and shrinkage. It means that melt mirrors are always oriented parallel to the rotation axis. At the same time, this allows us to consider the centrifugal casting regardless of the orientation of the axis of rotation in space, which makes it possible to simulate the process with a vertical, horizontal or inclined axis of rotation.

Preparing the Modeling Area


To start the calculation of centrifugal casting, the same data set is required as for gravity casting (sand casting or investment casting).

In this case, the orientation of the model must be set in such a way that the axis of rotation coincides with the direction of the Y axis of the coordinate system. To change the orientation of the model, see the Change Model Orientation, Position and Size) section.

The melt inlet should be set so that its center is located on the axis of rotation (or close to it). For more details on the inlets, see the Melt Inlet section.


It is also necessary to set the mold rotation axis and process parameters. To do this, in the Master module there is a special tool Rotation, which is launched with the button from the Process toolbar.

When the Rotation tool dialog box is open, left-clicking with the Ctrl key pressed on the model starts drawing a circle, the center of which marks the point of the axis of rotation (see the figure below). The direction of the axis of rotation corresponds to the normal to the surface on which the point is located. An arrow on the circle indicates the direction of rotation.


Rotation Tool

Rotation parameters are set in the dialog box:

  • Rotation speed, r/min. Sets the speed of rotation of the mold.

  • Coordinates, mm. This parameter displays the coordinates of a point lying on the axis of rotation, which is set interactively by clicking the left mouse button while holding down the Ctrl key. If necessary, you can edit this line manually to more accurately define the position of the point. Coordinates x y z are set sequentially in decimal numbers separated by a space.

  • Axis. This parameter is not editable, it shows the direction of the rotation axis. By default, it corresponds to the outward normal to the surface on which the point is specified. The direction of rotation of the shape will correspond to the rule of the right hand: if the thumb of the right hand corresponds to the direction of the axis of rotation, the rest of the fingers clenched into a fist indicate the direction of rotation of the mold. For ease of perception, an arrow is drawn on the circle showing this direction. To change the direction of rotation of the form to the opposite, it is necessary to press the button to the right of the output field of the direction of the rotation axis vector.

After all the rotation parameters are set, press the OK button. The dialog box closes, but the revolve path continues to display on the model.

Mold Gases

From the first moments of pouring the molten metal into a sand mold, an uneven pressure field develops within it. In the simplest case (casting into a non-gassing mold), this occurs due to the varying temperatures of the air within the pores of the mixture, located at different distances from the interface with the metal. The situation becomes more complex when the mixture has an unstable structure due to the presence of thermally unstable binders, burnable additives, and moisture, the gaseous products of combustion and evaporation of which are transported into adjacent layers by the pressure gradient.

When unfavorable conditions for gas evolution arise (the presence of cores surrounded by metal without gas venting channels, suboptimal density distribution in the mold with increased density in layers adjacent to gas emission zones, insufficient gas permeability of the mixture as a whole, etc.), a pressure gradient develops, directed toward the casting rather than into the mixture. If at this point there is no solid crust on the casting surface and the resulting pressure is capable of overcoming surface tension forces, gases are embedded in the metal. If they cannot be removed in the liquid phase of the casting, this leads to the risk of gas defects due to mold failure.

Modeling the mold's gas flow allows to identify the potential for gas defects in castings. Calculations allows to monitor the conditions for gas formation and release, allowing us to make decisions about process adjustments, core or mold ventilation, and density control during manufacturing.

In its current implementation, the model allows for the study of gas transfer conditions in a non-gas-forming mold and in a mold with a gas-forming binder (chemical source). Gas transfer processes in the mold are determined in parallel with thermal calculations based on the numerical solution of a differential equation describing the pressure field.

The following fields are available based on the calculation results:

  • time distribution of excess (over initial) gas pressure in the mold volume;
  • time distribution of the degree of transformation of the binder in the mold volume.

When viewing the calculation results, the following assumptions made in the first version of the gas exchange model implementation in the mold should be taken into account.

  • Gases are not removed through the inner surface of the mold, overcoming the backpressure of the liquid metal. The tendency for gas defects to occur should be judged by the accumulation of gases near the surface of the casting during the initial period of its formation, their increased pressure, and a weak tendency to migrate toward the outer surface of the mold.
  • An ideal contact is assigned between different mold volumes without the possibility of taking into account gaps, core prints, etc.

To calculate the gas regime of the mold, the following data is used, prepared in the Master module:

  • a model prepared in the same way as when performing solidification calculations (for more details, see the Preparing the Modeling Area section);
  • for mold materials that are capable of filtering and releasing gases, assign the corresponding characteristics to the mold material properties (see the Preparing the Modeling Area section below and Material Properties section);
  • initial conditions of the process and solver parameters (see the Mold Gas Simulation section).

Gas Transport Model in the Mold Volume


The calculation of gas transport in the mold is based on the numerical solution of equation, obtained on the basis of the equations of continuity, Darcy and the ideal gas state:


(1)

where P – pressure in the pores of the mixture;
          T – temperature;
          t – time;
          Kp – permeability coefficient;
          m – dynamic coefficient of gas viscosity;
          m – porosity of the mold material;
          M – molar mass of gas;
          R – universal gas constant;
          I – specific mass of gas released per unit volume of the mixture.

The derivative expresses the mass power of the gas release source and, if present, takes on a nonzero value. The description of the conditions for thermal degradation of the binder accompanied by gas release is based on the model:


(2)

where g – degree of conversion of the binder;
           f – specific gas-producing capacity;
           Mb – the initial mass of the binder per unit volume of the mixture (calculated based on its known
                    addition to the mixture and the density of the mixture in the mold);
           rg – gas density.

The degree of conversion is calculated according to the equation:


(3)

where n – chemical reaction order indicator;
           K – reaction rate according to the Arrhenius equation:


(4)

where Ea – activation energy of the thermal destruction process.

The analysis of the pressure field and the transformation of the substance to be gasified in different areas of the mold in the presence of chemical sources is performed on the basis of a joint numerical solution of the equations shown above and the Fourier equation.


Preparing the Modeling Area


The geometric model required for performing a numerical analysis of the mold gases is no different from the model for calculating the solidification of a casting in a mold, except in the case where the mold contains volumes that are impermeable to the gas flow. Different mold volumes must be assigned indices based on the materials used. For more information on indices and materials, see the Basic concepts: Volumes, Types, Indices and Materials Assignment sections.

For gas-producing or filtering materials, in addition to the standard set of properties specified in the Material Properties editor on the Thermal tab, you should specify parameters characterizing the filtration and gas emission conditions on the Gases tab. For more details see the Material Properties section.

If the mold contains impermeable materials (e.g., metal refrigerators) or if, for some reason, certain volumes in the mold are designated as impermeable, the Gas tab should be left blank when defining the properties of such a material. These volumes will be automatically recognized as impermeable during the calculation. When preparing the geometric model, you will also need to detuch such volumes (see the figure below) in the Volume Parameters window. Boundary 8, which is used to mark the interface between permeable and impermeable materials, must be assigned to the contact boundary of the detuched volumes with the remaining part of the mold. Boundary heat transfer conditions for the boundaty 8 should be assigned in the Heat Transfer editor (e.g., 5000-10000 W/m2 to simulate conditions close to ideal contact). For more information on separating volume indices and assigning boundary indices, see the Editing Volume Types and Indices and Editing Boundary Indices sections. The remaining boundaries are set as usual, including the default boundaries that can be used.


Detuching a volume with an impermeable material in the Volume Parameters window

It is necessary to ensure that the mold material density is specified in the Thermal tab of the property editor for the volumes involved in the outgassing calculation (see the figure below). This means that the specific heat capacity is entered in combination with the density, rather than just the volumetric heat capacity. Density is required as a parameter for calculating the initial (i.e., before the onset of thermal destruction) specific mass of the binder in the mixture. Specifying the density as a constant value will directly apply to this calculation. If a density curve is specified, the value corresponding to the minimum specified temperature will be used.


Mixture density data presented in the Thermal tab

Databases

In addition to preparing a geometric model, a lot of data is required for the calculation, such as: thermophysical and mechanical properties of materials and alloys, heat transfer parameters between bodies and with the environment, temperature, speed conditions, etc. Some of this data: material properties and heat transfer parameters, is mandatory for any calculation; others are used depending on the case. Editing, storage and other operations with this data are performed in the same type of editors, a description of which is given below in the relevant sections.


Project Overview Tab


The Project Data tree is a convenient tool for quickly viewing boundary conditions and data assigned in volumes and at the boundaries of the calculation area, as well as a starting point for assigning new conditions. The tree is located on the Project Overview tab (1), which is located on the left side of the Master module (see the figure below).


Project Overview Tab

The top of the tab contains filters that can be used to change the amount of information displayed in the tree (2). When loading a project, the In Project filter is automatically set to display only information about the project data in the tree. The volume and boundary indices assigned in the project are indicated in the tree by the "(+)" sign at the beginning of the line.

Each tree position has a button (3), at the end of the line, which launches the corresponding data editor. To delete data for a selected tree line, for example, a material, click the button  (4). If this icon is missing, the data has not yet been created.

For many calculation conditions, such as Mold Materials, Heat Transfer, Translations, the data can be loaded from previously saved files, or, conversely, exported, in order to use them in other projects. To export data for all conditions, use the button , located in the data set line (5), and to import, use the button located there. The POLIGONSOFT database contains ready-made settings for parameters and materials for different casting technologies, which can be loaded into the project in a similar way and, if necessary, edited for a specific task. Some conditions, such as "Casting Material" or "Dynamic Environment" do not have import and export buttons (6), they can only be imported or exported directly from the editor.


Data Editors. General Controls and Features


The material properties, boundary conditions and other data editors, despite the fundamental differences, have the same interface. The window of any editor contains a toolbar (1), a table of editable parameters (2) and an area for editing curves (3) (see the figure below). If there is a * symbol at the end of the editor's caption, it means that the parameters have been changed. When you close the editor, you will be prompted to save your changes.


Temperature Modes Editor

In some editors under the toolbar (1) there is a line where you can enter the name of the data set. This name will then be displayed in the corresponding line in the parameters table (Start Simulation and Parameters window). Below is a line for entering a comment on the data set. In most cases, the dataset name is edited in the tree on the Project Overview tab.

When creating or editing data for a specific volume index, boundary, etc., the editor is launched from the corresponding tree on the Project Overview tab. If the data has not yet been created, all elements of the editor interface are inactive (see the figure below). To start working in the editor, click the New button, which will activate its interface.


Empty Translation editor

Further in this section, interface elements and buttons common to all editors are considered. The following sections give more detailed descriptions of all editors and the data edited in them.


Toolbar Buttons


All editors have a toolbar at the top with buttons for general data operations (see figure below). The number of commands may vary depending on the functionality of the editor.


Editor toolbar



The New button creates a new data set with "default" parameters. The data is created for use in the current project.



The Import button opens a standard dialog for loading data into the project. It is possible to import data created in older versions of POLIGONSOFT.



The Export button export data from the project to files that can form a user database and then be used in other projects. This command is also used if you want to save edited data in the original file.



The Apply button saves the imported or/and edited data in the project. The original data file remains unchanged. To update it, use the Export button.

The Close button closes the editor. If the data in the editor has been changed, a message will be displayed asking you to save it.

On the right side of the toolbar are filters that change the set of displayed data:

  • Current. Shows only the currently selected position in the tree, for example, the mold material for the selected volume index. This filter is set automatically when you launch the editor from the tree in the Project Overview panel.
  • In Project. Shows only project data, including empty ones.
  • With Data. Shows all positions where data is defined, including those that are not assigned in the project. This filter is set automatically when you launch the editor from the menu or the Process toolbar.

Editing Constants and Functions


Many material characteristics and simulation conditions can be specified either by constants or by curves. In this case, the drop-down list is active in the parameters table, in the Mode column (see the figure below).


Editing the ambient temperature in the Heat Transfer Editor

If Value is selected, then the parameter should be set as a constant in the Value column. If Curve is selected, the curve editing area at the bottom of the editor becomes active. It contains: on the left - a table with coordinates of points and a toolbar for working with a curve; on the right is the curve drawing area (see the figure below).


Curve Editing Area

Above the table with the coordinates of the curve points there is a toolbar with buttons that control the data in the table.




The Add Row button creates a new row in the table. The row is created under the current (active) table row. New row cells default values based on the values set in the previous and next cells. At the same time, a new point is displayed in the curve drawing area, it is highlighted in yellow.



The Delete Row button deletes the selected row in table.



The Import Curve Points button opens a standard dialog for loading data into a table. It is possible to import data from .csv (MS Excel), .lin (POLIGONSOFT) files and any other text files containing curve points in the "x y" format.



The Export Curve Points button exports the data from the table to a .lin text file (POLIGONSOFT) that contains the curve points in the "x y" format. The file can be imported into other spreadsheet editors such as MS Excel.



The Copy Table button copies the table to the clipboard. The data can be pasted into another table or into another application.



The Paste Table button pastes data from the clipboard.

The Delete Table button deletes all table rows. To delete one or more selected table rows, use the button.

Editing the table results in an immediate change in the appearance of the curve in the drawing area on the right. The curve must contain at least two points. The table of editable parameters will display the minimum and maximum values of the parameter.

The drawing area is interactive, the user can change the appearance of the curve and edit it with the mouse.

  • Clicking the left mouse button on a point highlights it in yellow. At the same time, the corresponding row in the table is highlighted. You can select any number of points by clicking on them with the mouse and holding down the Ctrl button on the keyboard.
  • The selected point can be moved by holding down the left mouse button. At the same time, the values of the coordinates of the point in the corresponding line of the table change.
  • By holding down the left mouse button, you can move the entire curve over the drawing area.
  • Using the mouse wheel, you can change the scale of the curve display.



The Fit button returns the curve display to its default state.

The Line button allows you to create a linear section of the curve between two selected points. When two points are selected, the button becomes active. On fig. below is an example of how this command works.


Initial curve
Curve with linear section
Example of creating a linear section of a curve


Computed Parameters


Some editors have parameters whose values ​​can not only be set by the user, but also calculated using other set parameters. For example, some characteristics of a casting alloy, necessary for the shrinkage porosity model, can be calculated using other properties, such as density. If a parameter can be calculated, there is a button with a calculator in the right side of the cell where the value is usually entered (see the figure below).


Button in the parameter cell to calculate it

The calculated values ​​can be edited or replaced with other values ​​at the user's discretion. The requirements for the calculated parameters and the calculation methodology are described in the relevant sections.


Material Properties


The Material Properties editor (see figure below) is used for creating and editing materials used in the calculation. Materials are stored in the %POLIGONSOFT%\Database folder in pmat files (older versions have bd, fil and bdd extensions). Each material must contain a set of characteristics required for POLIGONSOFT solvers: Euler, Fourier, Hooke, etc.


The button to launch the Material Properties editor is located on the Process toolbar and on the Main tab of the ribbon menu. The editor can also be opened from the Materials tree on the Project Overview tab and from the Start Simulation and Parameters window. See Materials Assignment section for more details.


Material Properties Editor

Under the toolbar there is a line with the name of the material, which can be changed. The comment box is located below. The following is a table of material properties, with properties divided into tabs:

  • Composition;
  • Flow;
  • Thermal;
  • Porosity;
  • Stress;
  • Exothermic;
  • Grain Growth;
  • Hardenability & CCT.

Tabs can be of three colors (see the figure below): gray - the tab does not contain data; green - the tab contains the minimum set of material properties that are sufficient to run the solver; pink - the set of properties is not enough to run the solver.


Different tabs colors

On the tab marked in pink, the properties that are not set but required to run the solver are also highlighted in pink (see the figure below).


Highlighting undefined properties

The ability to import material properties in a legacy format is retained. Older versions of materials have the extensions bd (thermal properties), fil (shrinkage properties) and bdd (mechanical properties). When importing a bd (thermal properties) file, other file types with the same name are searched. If other alloy properties are found, the Master offers to load them into the appropriate tabs (see the figure below).


Message suggesting importing material properties

When importing fil (shrinkage properties) and bdd (mechanical properties) files, the properties are loaded into the respective tabs, but other file types are not searched.


Composition Tab


The tab contains a table of chemical composition (see figure below). In the first line of the table, you can select the base element of the alloy (Fe, Al, Ni, Ti, etc.) and then enter the content of the remaining elements in weight percent.


Material Properties, Composition Tab

The tab turns green if the content of at least one chemical element or a basic element is specified in the table.

The data in this tab can be used as input for the Heat Treatment solver (for more details see the Heat Treatment Simulation) section).


Flow Tab


The tab contains material characteristics related to its behavior in the liquid state (see figure below).


Material Properties, Flow Tab

Dynamic Viscosity is required for the continuous casting hot tearing model. It can also be used to calculate the Kinematic Melt Viscosity parameter, which is used to model microporosity (see the Porosity Tab for more details).


Thermal Tab


The tab contains a set of material characteristics required to calculate the temperature and phase fields (see the figure below). To calculate temperature fields in bodies where there is no phase transition (usually it is a mold), the minimum set of thermophysical properties required for calculation includes:

  • Thermal conductivity;
  • Specific heat;
  • Density.

These properties can be specified either as constants or as functions of temperature.

In the presence of a phase transition, i.e. solidification (usually these are bodies of the “casting” type), you should additionally specify the following characteristics:

  • Latent Heat;
  • Fraction of Solid.

Material Properties, Thermal Tab

The Solid Fraction is a function of temperature. Its values should start from 0 at the liquidus temperature and increase to 1 at the solidus temperature. Thus, the fraction of the solid phase is determined by a falling curve strictly between 0.0 and 1.0. Isothermal transformation (for example, eutectic) is recommended to be "stretched" by an interval, for example, of 2-3 degrees. Latent Heat is given by a constant.


It should be noted that the POLIGONSOFT algorithms allow you to set the fraction of the solid phase with flat sections. The recommendation to "stretch" them by 2-3 degrees is not mandatory.


An alternative way of defining the phase transition is possible. Instead of Specific Heat and Latent Heat, the corresponding Enthalpy curve can be given as a function of temperature H(T), which is defined as follows:


where cp(T) – specific heat; Q – latent heat; fs(T) – fraction of solid.

Since there are two ways for defining a phase transition, the editor automatically avoids conflict between the following options:

  • Specific Heatcp(T);
  • Latent Heat Q;

or

  • Enthalpy H(T).

If specific heat or latent heat is entered, the enthalpy entry is blocked, and vice versa.

Thermal conductivity and enthalpy must always be given as a temperature curve.


Porosity Tab


The tab contains a set of material properties needed to calculate pipe shrinkage, macro- and microporosity in a casting (see the figure below).


Material Properties, Porosity Tab

The set of parameters required for the "standard" macroporosity model used in POLIGONSOFT:

  • Melting Volume Change;
  • Liquid Phase Volume Change;
  • Pl1 – fraction of liquid phase at the beginning of linear shrinkage;
  • Pl2 – fraction of liquid phase at the end of the gravity flow;
  • Pl3 – fraction of liquid phase when channels are fully blocked;
  • Volume Change Spectrum dV = f(fs).

Some parameters can be calculated based on the properties specified in the Thermal tab. This mainly requires the Density specified as a function of temperature and the Fraction Solid curve. The calculated values ​​can be edited or replaced with other values ​​at the user's request.

To calculate the parameter, you should click the calculator button in the right side of the cell (see the figure below).


Button in the parameter cell to calculate it

Melting Volume Change β, %. The parameter, better known as “volume shrinkage”, is calculated using the equation:


where ρliq and ρsol – density in kg/m3 at liquidus and solidus temperatures, respectively.

In this case, the density values ​​in the liquid and solid state are automatically calculated using the specified density curve at the solidus and liquidus points. If the required temperature values ​​are not among the points of the density curve, they will be found using the linear interpolation method for the neighboring points. This approach is used when calculating all other parameters.

Liquid Phase Volume Change βL, 10-6/K. The parameter is the coefficient of thermal linear expansion of the liquid phase above the liquidus temperature Tliq. It is calculated using two points of the density curve: one point, for example, is the last point of the curve at the maximum temperature Tmax, the second point is the density at Tliq:


where ρliq – density in kg/m3 at liquidus temperature; ρTmax – density at temperature Tmax; Tmax > Tliq.

Volume Change Spectrum dV(fs). The curve characterizes the change in volume during the solidification process. The horizontal axis shows the values ​​of the Fraction Solid, specified by the fs(T) curve on the Thermal tab. The vertical axis shows the corresponding Volume Change, calculated using the equation:


where ρfs – density in kg/m3 at a temperature corresponding to the current fraction of the solid phase.

It is easy to see that when fs = 0 the density ρfs = ρliq and dV = 0, and when fs = 1 the density ρfs = ρsol and dV = 1.

To calculate microporosity according to the "standard" POLIGONSOFT model, in addition to the parameters listed above, it is required to set the parameters:

  • Kinematic Melt Viscosity;
  • Liquid Strength.

To calculate the filtration coefficient using the equation:


where
fL и fs – fractions of liquid and solid phases, respectively; λII – secondary dendrite arm spacing

need to define:

  • coefficient Ко;
  • power n.

Kinematic Melt Viscosityγ, 10-5 m2/s can be calculated if the dynamic viscosity curve is specified in the Flow tab and the density curve in the Thermal tab. The parameter characterizes the melt flow in the two-phase region near the solidus temperature and is calculated using the equation:


where Ny – melt Dynamic Viscosity, Pa∙s. Since the dynamic viscosity is specified in centipoise, a coefficient of 0.001 is used to convert it to Pa∙s; ρT – is the density in kg/m3 at a temperature corresponding to the fraction of the solid phase at which microporosity is formed. The temperature Tsol+10° is used in the calculation.

To calculate the parameter, you should click the calculator button in the right side of the cell.

At the end of the table, the set of parameters required for the new porosity model is listed, which extends the capabilities of the "standard" macro- and microporosity models:

  • Fraction of liquid phase at the beginning of linear shrinkage;
  • Fraction of liquid phase when channels are fully blocked;
  • Critical Pressure Pcrit;
  • Melt Compressibility Modulus E;
  • Surface Tension Coefficient σ;
  • Secondary Dendrite Arm Spacing λII.

Some parameters from this list duplicate the parameters of the "standard" macroporosity, but nevertheless must be set again. This is because the parameter values are different when used in different models and their separation is convenient for switching between new and standard porosity models.

For more information on the new model of porosity and its parameters, see the New Pipe Shrinkage and Macroporosity Model section.


Stress Tab


The tab contains a set of material properties required to calculate the stress-strain state of the casting and the mold (see the figure below).


Material Properties, Stress Tab

Different stress models can be used in the calculation, each requires its own set of mechanical characteristics. The stress model is selected from the list, which is located in the first line of the parameter table. The following options are available:

  • Linear Elastic;
  • Elasto-plastic.

The Linear Elastic model is characterized mainly by Young's Modulus. In addition, the Poisson's Ratio and the Thermal Expansion Coefficient must also be specified. The value of Poisson's Ratio for metals is approximately 0.3.

For an Elasto-plastic model, the properties described above for an elastic model, i.e. Young's Modulus, Poisson's Ratio and Thermal Expansion Coefficient must also be given. In addition to this, the Yield stress and Hardening Coefficient must be specified. Hardening Coefficient, as a first approximation, can be calculated using a simple formula:


H(T) = E(T) / 20,

where E(T) is Young's Modulus, defined as a function of temperature.

To calculate the parameter, you should click the calculator button in the right side of the cell.

All listed properties can be specified by constants or temperature curves. Of course, it is extremely undesirable to set the properties with a constant, since in this case the accuracy of the calculation decreases. For example, if all properties are given by constants, then the result of the calculation will in no way depend on the change in the temperature field of the casting and the mold.

The source of stress properties can be: literature, the Internet and scientific research. If these sources cannot provide the necessary data, then they can be obtained empirically, i.e. by carrying out the necessary mechanical tests.


Exothermic Tab


The tab contains a set of parameters (see the figure below) necessary for heat release simulation when using mold elements made of exothermic (thermitic) mixtures (sleeves, covers, mats, etc.). Such a material must be specified in the model as a separate volume of the “mold” type with the necessary thermophysical properties are specified on the Thermal tab (for more details see the Thermal Tab section). Then the required parameters for the exothermic reaction should be specified:

  • Exothermic Energy. The amount of heat released by 1 kg of material during the period from the beginning to the completion of the thermite reaction.
    Specified as an integer in kJ/kg


  • Burnt Fraction. Sets the total duration and rate of the reaction (see the figure below). Generaly, the reaction may proceed unevenly with time intervals of deceleration and acceleration of the process, which can be taken into account in the calculation.
    It is given by a time - reaction fraction curve, where the time is given in seconds.


  • Ignition Temperature. The parameter sets the temperature of the thermite reaction (ignition), which "switch on" the calculation of heat release in the corresponding mesh node of the exothermic material.
    Specified as an integer in Celsius.

Material Properties, Exothermic Tab

If at least one of the three parameters described above is not set, the tab will be colored pink. This means that there is not enough data to run the calculation. The corresponding parameter lines will also be pink. If all parameters are set, the tab is colored green.

Additionally (but not necessarily) the thermophysical properties of the material after completion of the reaction can be specified. They differ from the initial exothermic material properties, specified on the Thermal tab, and can be used in the calculation to take into account the heat-insulating ability of burnt material. The following properties will be used when calculating temperatures for those mesh nodes where the burnt fraction is equal to 1:

  • Conductivity;
  • Density;
  • Specific Heat.

These properties can be specified by constants, or temperature curves. All three properties must be set. If the properties of the material after the completion of the reaction are not set, then properties from the Thermal tab will be used for the heat transfer calculation in the areas, where the fraction of burnt material is equal to 1.

When using exothermic materials, it may be necessary to enter specific heat transfer parameters at the casting-material and material-mould interfaces. These parameters are set in the Heat Transfer editor (see the Heat Transfer Parameters section). For more information on modeling foundry technology using exothermic materials, see the Exothermic Sleeves section.


Grain Growth Tab


The tab contains coefficients determined by the chemical composition of the alloy, which are used in the equation for the growth rate of the primary axis of the dendrite in the grain growth model (see figure below). This data is used by the Macrostructure solver. For more information about macrostructure simulation, see the Macrostructure Simulation section.


Material Properties, Grain Growth Tab


Hardenability & CCT Tab


The tab contains the data required to calculate the microstructure and properties in the model bodies of the "casting" type, using the following models:

  • Hardenability-based heat treatment model for steel;
  • groups of physically based models (GPBM).

For more information on modeling heat treatment processes using different models, including those listed above, see the Heat Treatment Simulation section.

For calculations using the model based on hardenability, it is necessary to specify data characterizing the conditions and results of standard (EN ISO 642) tests for hardenability (see the figure below):


Material properties for calculations using the hardenability model

  • Austenitizing Temperature. Heating temperature of a standard sample for hardening in Celsius.
    Enter an integer value in Celsius degrees.


  • Austenitizing Time. The holding time of a standard sample for quenching at the austenitizing temperature. As a rule, this is 0.5 hours.
    Enter a float value in hours.
    Default value: 0.5.


  • Grain Size Number. Additionally (but not necessarily), the austenite grain size number obtained with the specified austenitization parameters can be specified. The presence of this parameter will allow for a correction to be made for the actual austenitization conditions when modeling the HT in the user's task and thereby increase the accuracy of the calculation.
    Enter an integer value.


  • Hardness. Hardenability curve – distances from the hardened end of the sample in millimetres and the corresponding hardness values ​​in Rockwell units (HRC).

For calculations using GPBM, it is necessary to enter the thermokinetic diagram (CCT) of austenite decomposition during continuous cooling and the conditions for obtaining it (see the figure below). The conditions in which the CCT was obtained are described by the following parameters:

  • Austenitizing temperature – the temperature to which the samples were heated during testing to obtain the CCT diagram.
    Enter an integer value in Celsius degrees.


  • Austenitizing Time.The holding time of a sample for quenching at the austenitizing temperature. As a rule, this is 0.5 hours.
    Enter a float value in hours.
    Default value: 0.5.


  • Grain Size Number. Additionally (but not necessarily), the austenite grain size number obtained with the specified austenitization parameters can be specified. The presence of this parameter will allow for a correction to be made for the actual austenitization conditions when modeling the HT in the user's task and thereby increase the accuracy of the calculation.
    Enter an integer value.


  • A3 Temperature. Temperature value for critical point A3. If the steel has no ferritic transformation, enter zero.
    Enter an integer value in Celsius degrees.


  • A1 Temperature. Temperature value for critical point A1. If the steel has no pearlitic transformation, enter zero.
    Enter an integer value in Celsius degrees.


  • Threshold for initiating austenite decomposition. The value corresponding to the volume fraction of decomposed austenite, taken as the beginning of the transformation.
    Enter in %.
    Default value: 1.


  • Equilibrium Ferrite. The volume fraction of ferrite that is obtained under equilibrium conditions. Can be taken approximately as corresponding to conditions of very slow cooling (e.g. at 0.01°C/s).
    Enter in %.


  • Equilibrium Pearlite. As in the case of ferrite, the corresponding lowest cooling rate (e.g. 0.01°C/s) can be taken as approximately.
    Enter in %.

The sum of the proportions of ferrite and pearlite under equilibrium conditions should be 100%.

To draw the CCT diagram, a separate data entry will be required for each of its branches in the form of pairs of values ​​“cooling rate in °C/s – transformation temperature in °C”.

  • Ferrite Formation Start Temperature. Entering this curve is mandatory. The curve must be described by at least 3 points. If the steel does not have a ferrite transformation, enter 3 any cooling rate values and set the corresponding temperatures equal to zero.


  • Ferrite Volume Fraction. Volume fractions of ferrite corresponding to the cooling rates introduced above. Entering this data is optional.


  • Pearlite Formation Start Temperature. Entering this curve is mandatory. The curve must be described by at least 3 points. If the steel does not have a pearlite transformation, enter 3 any cooling rate values and set the corresponding temperatures equal to zero; in this case, it is permissible not to enter the curve for the end of the transformation.


  • Pearlite Formation Finish Temperature. Curve of pearlite transformation end. If data is available, input is recommended. When inputting, the curve must be specified by at least 3 points.


  • Pearlite Volume Fraction. Volume fractions of pearlite corresponding to the cooling rates introduced above. Entering this data is optional.


  • Bainite Formation Start Temperature. Entering this curve is mandatory. The curve must be described by at least 3 points.


  • Bainite Formation Finish Temperature. Curve of bainite transformation end. If data is available, input is recommended. When inputting, the curve must be specified by at least 3 points.


  • Bainite Volume Fraction. Volume fractions of bainite corresponding to the cooling rates introduced above. Entering this data is optional.


  • Martensite Formation Start Temperature. Entering this curve is mandatory. The curve can be degenerate and, as an exception, defined by only one point.


  • Martensite Formation Finish Temperature. Volume fractions of martensite corresponding to the cooling rates introduced above. Entering this data is optional.


  • Martensite Volume Fraction. Volume fractions of martensite corresponding to the cooling rates introduced above. Entering this data is optional.


  • Hardness. The obtained hardness values ​​in Vickers units correspond to the cooling rates introduced above. The cooling rate column is not editable and is generated automatically based on all the values ​​entered above when specifying the curves for the beginning and end of phase formation. Entering hardness values ​​is optional.

When entering data to draw the CCT diagram, the following rules should be followed:

  1. The temperatures of the curves of the beginning and end of one transformation must be specified for the same cooling rates. For example, if the temperature value for the rate of 0.5 °C/s is specified for the Pearlite Formation Start Temperature curve, then the temperature value on the Pearlite Formation Finish Temperature curve must be specified for the same cooling rate. If the value of the end of transformation temperature on the original CCT is missing, 0 should be specified.
  2. The curves of the beginning and end of transformations should be specified by at least three points.
  3. When entering volume fractions, it is necessary to ensure that the sum of the structural components for one cooling rate is equal to 100%.

Input control can be performed using the diagram constantly displayed at the bottom (see the figure below). For convenience, branches on the diagram are also plotted in the coordinates "cooling rate - temperature" (default display type of CCT diagram in the editor window).

Usually, the curves of the beginning and end of the transformation (except for martensitic) are convex upwards. This corresponds to the fact that the lower the cooling rate, the higher the temperature of the beginning or end of the transformation. This allows the most correct transformation of the CCT diagram into an isothermal decomposition diagram (TTT) and the description of the family of C-shaped curves by Zener-Hillert dependencies.

The input of the curve of the end of the ferrite transformation is not provided, since it coincides with the curve of the beginning of the pearlite transformation.


Material properties for calculations according to GPBM

To display the CCT diagram in the form familiar to a thermal engineer in the “temperature-time” coordinates, you should check the CCT in Temperature-Time coordinates checkbox, it is located under the diagram.

The hardenability curve and thermokinetic diagram can be calculated. To access the corresponding generator (see the figure below), click the calculator button in the parameter group name line. Calculations are based on the steel's chemical composition and austenite grain size, which can also be calculated based on the austenitizing temperature and time. The data obtained in this way can be used in heat treatment process analysis using the hardenability-based model and GPBM, and are also generally useful for designing heat treatment modes. For more information on performing calculations, see the Heat Treatment Simulation section.


a)
b)
Dialogs for entering initial data for generating: а) hardenability; б) CCT

The hardenability curve is calculated using the Grossman method [1], modified to take into account the austenite grain size. The result of the calculation is a hardenability curve in the coordinates "distance from the hardened end of the sample" (mm) – "hardness" (HRC), which is obtained in practice during standard testing.

The CCT calculation is based on the initial calculation of the TTT diagram using the metod [2] followed by its transformation into the CCT using the Sheil additivity rule [3]. A series of other models were also used in the complex calculation:

  • when calculating critical points A3, A1, Bs, Ms;
  • model for calculating the equilibrium content of ferrite and pearlite;
  • Johnson-Mehl-Avrami-Kolmogorov equation;
  • Creusot-Loir formulas for calculating hardness.

The result of the CCT calculation is a series of cooling rates and the corresponding temperatures of the beginning and end of the transformations of ferrite, pearlite, the temperature of the beginning of the transformation of martensite, bainite (if the transformation occurs at a given speed), the volume fractions of ferrite, pearlite, bainite, martensite (if the transformation occurs) and hardness.

Generators can only be used to calculate the properties of hypoeutectoid steels. Chemical composition restrictions are listed in the table below.


Table. Permissible limits for chemical element content when using generators.
Element Content limits, %
Basic elements
C 0.05-0.60
Si ≤1.37
Mn ≤2.00
Cr ≤3.00
Ni ≤4.00
Mo ≤1.00
Residual levels
Cu ≤0.15
Al ≤0.06
Nb ≤0.02
V ≤0.20
Ti ≤0.20
W ≤0.20
B ≤0.0008
N ≤0.012
Co ≤0.10
Mg ≤0.005
Sn ≤0.01
Zr ≤0.005
S ≤0.035
P ≤0.035

After launching the hardenability or CCT generator, you must assign the initial data for the calculation in its window. After the generator completes, the calculation results will be saved in the corresponding section (Hardenability or CCT) of the editor. The calculations are performed independently: only hardenability, only CCT, or both can be generated.

To obtain the result, you should act in the following order.

  1. In the Composition tab, ensure that Fe is selected as the base element, and that the steel composition is entered correctly. If you are entering a new alloy, fill out the tab. The steel composition must be compatible with the elemental content restrictions specified above.
  2. Open the generator settings window by clicking the button with the calculator icon in the desired section (Hardenability or CCT).
  • By default, the Grain Size Number parameter is set to 7 (can be edited).
  • Only in the CCT Calculation window: the "Threshold for initial austenite decomposition" parameter is additionally set at the default level of 1% (can be edited).
  • The input fields for other parameters (Austenitizing Temperature and Austenitizing Time) are empty by default (can be edited).
  • Any parameter values ​​entered by the user during previous launches are not restored when the generators are reopened.
  1. Please fill in the required input fields.
  • To perform a standard calculation, the values ​​in the default window can be accepted, then during the generation process the values ​​for the parameters Austenitizing Temperature and Austenitizing Time will be calculated automatically.
  • When entering the Austenitizing Temperature and Austenitizing Time parameters, you can either specify a value for the Grain Size Number parameter or leave the parameter field blank. In the latter case, the required value will be selected automatically during generation.
  • The "Threshold for initial austenite decomposition" parameter (in the CCT Calculation window) can be either specified or not - in the latter case, it will be assigned automatically during generation.
  • If one of the parameters Austenitizing Temperature or Austenitizing Time is specified, then it is always necessary to enter a value for the second one.
  1. Input data control in the generator window is maintained automatically.
  • If the entered parameter combinations are invalid, the lines that need to be filled in are highlighted and the Continue button is blocked.
  • Parameters can only be entered within the permitted limits (see table below).

Table. Restrictions on the input of initial data for generating hardenability and CCT
Parameter Range

Austenitizing Temperature, °C

750 - 1200

Austenitizing Time, h

0.1 - 5

Grain Size Number

-7 17

Threshold for initial austenite decomposition

0.1 - 1

  1. To start the calculation, after entering the initial data, click the Continue button..
  • If no errors are found, the required properties will be generated and a confirmation prompt will appear asking to fill in the corresponding lines of the Hardenability and CCT tab, replacing the previously entered values. Upon confirmation, the calculated data will be transferred to the material editor tab. If you cancel, the calculated data will be deleted, and the Material Editor tab will remain unchanged.
  • If the alloy type does not match and the chemical composition is outside the permissible limits, the calculation is interrupted with an error message.
  • If the user has specified the Austenitizing Temperature parameter, but the value is below the A3 point, the calculation is rejected: an error message will be displayed indicating the A3 temperature value for the entered chemical composition.
  • If the alloying element combination classifies the steel as hypereutectoid, a corresponding message will be displayed and the calculation will be rejected. To obtain the result, the C, Cr, and W content will need to be reduced, or the Si and Ni content increased.
  1. ASTM A 255-2014 – Standard Test Methods for Determining Hardenability of Steel.
  2. Li M., Niebuhr D., Meekisho L., Atteridge D.G. A computational model for the prediction of steel hardenability. //Metallurgical and Materials Transactions B, 29, 1998, pp. 661-672.
  3. Chen X., Xiao N., Li D., Li G., Sun G. The finite element analysis of austenite decomposition during continuous cooling in 22MnB5 steel. //Modelling and Simulation in Materials Science and Engineering, V. 22, Nr. 6, 2014.

Gas Tab


The Gas tab contains a set of properties (see the figure below) necessary for modeling the mold's gas regime. The parameters on this tab are divided into two groups. The properties in the Filtration Properties of Compacted Sand Mixture group are needed to calculate the transport of gases driven by pressure differences due to uneven temperature distribution within the mold and gas formation during binder gasification. Gasification parameters are specified in the Gas Evolution Due to Binder Decomposition group. If no characteristics are entered in the second group, the material will be considered non-gassy. The movement of gases in a non-gassy form is determined solely by their thermal expansion within the pores during mixture heating.


Properties for mold gas calculation

In the Filtration Properties of Compacted Sand Mixture group the following properties are specified.

  • Permeability Coefficient. Specifies the permeability coefficient of the porous mold material.
    Specified in m2 by a constant value or a temperature curve.


  • Dynamic Viscosity. Viscosity of the gas filling the pores.
    Specified in Ps·s by a constant value or a temperature curve.


  • Gas Permeability. The gas permeability of the molding sand, determined by a widely used testing method using a standard laboratory apparatus. Gas permeability Krel is related to the permeability coefficient Kp by the formula Kp = 3.094·10-13Krel.This row is available for entry if no data is entered in the Permeability Coefficient and Gas Permeability Coefficient rows.
    Enter a floating (real) value in std. units.


  • Gas Permeability Coefficient. This parameter can be specified instead of the Permeability Coefficient, Dynamic Viscosity and Gas Permeability parameters.
    Specified in m4/(H·s) by a constant value or a temperature curve.


  • Porosity of Compacted Sand. Specifies the mold material porosity after compaction.
    Enter a floating (real) value in %.

Depending on the available data, the following combinations of parameters can be entered:

  • Permeability Coefficient
  • Dynamic Viscosity
  • Porosity of Compacted Sand,

or:

  • Gas Permeability Coefficient
  • Porosity of Compacted Sand,

or:

  • Dynamic Viscosity
  • Gas Permeability
  • Porosity of Compacted Sand

In the Gas Evolution Due to Binder Decomposition group the following properties are specified.

  • Molar Mass. The molar mass of the gas in the pores of the mixture is specified.
    Enter a floating (real) value in kg/mol.


  • Binder Content. The proportion of thermally unstable binder is specified.
    Enter a floating (real) value in %.


  • Specific Gas Evolution. The gas formation capacity of a unit mass of a thermodestructing component is specified.
    Enter a floating (real) value in m3/kg.


  • Initiation Temperature for Binder Thermodestruction. The temperature at which the thermodestruction reaction begins. If the temperature in some area of the mixture drops below the specified threshold, the reaction stops.
    Enter a floating (real) value in Celsius degrees


  • Degree of Binder Decomposition Chemical Reaction. Used in the equation for calculating the degree of binder conversion (for more details, see the Gas Transport Model in the Mold Volume section).
    Enter a floating (real) value.


  • Pre-Exponential Factor for Arrhenius equation (for more details, see the Gas Transport Model in the Mold Volume section).
    Enter a floating (real) value in s-1.


  • Activation Energy of Binder Decomposition Reaction. Activation energy in the Arrhenius equation (for more details, see the Gas Transport Model in the Mold Volume section).
    Enter a floating (real) value in kJ/mol.

Heat Transfer Parameters


The Heat Transfer editor is designed to create and edit heat transfer settings at the boundaries between volumes (Body-to-Body) and at the outer boundaries (Body-Environment). The settings can be applied to solve one specific problem or exported to a file with the afo extension for use in other calculations, including as a casting process template.


The button to launch the Heat Transfer editor is located on the Process toolbar and on the Main tab of the ribbon menu. The editor can also be opened from the Start Simulation and Parameters window and from the Heat Transfer tree on the Project Overview tab.

In the editor table, heat exchange conditions are specified for all surfaces of each volume of the computational domain. Heat exchange parameters for each surface (boundary) are determined by a specific task and depend on where this boundary is located. The list of boundaries does not include a boundary with index 0, since this boundary is either internal or external with adiabatic conditions and no parameters need to be specified for it. To display in the table only those boundaries that are specified in the model, you can use the In Project filter located on the toolbar (see the figure below). To display a list of boundaries where data is specified (not all boundaries may be present in the project), use the With Data filter. The Current filter shows only the boundary selected in the tree on the Project Overview tab.


Heat Transfer editor, In Project filter

The list of heat transfer parameters for each boundary begins with a line that specifies its index and color, which is used for visualization in the pre- and postprocessor in the boundary indices mode. In the same line, in the Name column, you can enter an explanatory text describing, for example, the heat transfer conditions. This text will be displayed in preprocessor dialogs when editing boundaries. The next column Mode indicates one of the two heat transfer modes. For boundaries with indices 1-9, it can be changed. Сlicking on a line brings up a drop-down list in which you need to select the desired option (see the figure below):

  • Body-to-Body mode sets the heat transfer between the casting and the mold, or between different parts of the mold. In this mode, heat is transferred across a selected boundary from a body with a higher temperature to a body with a lower temperature. It is assumed that the boundary with this index (and this color) is set on the contact surfaces of both bodies. If in the computational domain one of the outer surfaces of the body has a boundary for which this mode is specified, but the second body is absent, then the heat transfer will go between the real and the virtual body. For example, if the boundary with this mode is set on the outer surface of the casting, and the body of the mold is absent, then the heat transfer will go into a virtual mold of infinite thickness. In this case, the thermal activity of the mold will be taken into account (the root of the product of the volumetric heat capacity and thermal conductivity). The characteristics of the mold material from the specified set of mold materials with an index equal to the boundary index will be used. For more information on the virtual mold, see the Virtual Mold section.


  • Body-Environment mode sets the heat exchange between the volume (casting or mold) and the environment with a given temperature.

Changing the heat transfer mode for the boundary

As indicated in the Basic concepts: Volumes, Types, indices section, the first nine boundaries, starting from index 1, are universal and can be used as boundaries between bodies of the computational domain or as external boundaries with the environment. All other boundaries with indices 10-18 set the parameters of heat transfer between the body and the environment.

In the Body-Environment boundary mode, the Ambient Temperature line becomes available.The ambient temperature can be set by a constant value or by a time-temperature curve. To select, you need to click the left mouse button in the Mode column and select the desired option from the drop-down list (see the figure below).


Editing the ambient temperature for the boundary with index 2

If the ambient temperature is constant, select the Value mode and set the temperature in the next cell Value. If the ambient temperature changes depending on time, select the Curve mode. In this case, the table and the curve drawing area in the lower part of the editor window become active, where you can set the points of the curve in the time-temperature coordinates (see the figure below). For more information on editing curves, see the Editing Constants and Functions section.


Ambient Temperature curve

It should be note in mind that for boundaries 17 and 18, the ambient temperature is set by a constant in the Start Simulation and Parameters window.

The next four lines set the heat transfer parameters, in accordance with the POLIGONSOFT model . This model can be represented as a scheme (see figure below) and describes three mechanisms of heat transfer: heat conduction, convection and radiation.


Scheme of the heat transfer model in POLIGONSOFT

It is assumed that, in the general case, at the boundary between volume 1 and volume 2, two successive layers can exist: let us designate them as I and II. Layer I is characterized by the heat transfer coefficient of the layer α, which is set by the user. For example, this could be a coating applied to a mold. If the thermal conductivity of the coating material λI and its thickness δ are known, then the heat transfer coefficient of the layer α can be calculated by the equation:


(1)


Layer II is characterized by a calculated heat transfer coefficient β. In general, the contact between volumes 1 and 2 is not ideal. It is assumed that heat transfer in layer II can take place in three ways: heat transfer in the zone of direct contact of volumes, convection and radiation where there is no contact. Each mechanism acts on a part of the boundary, which is specified by the relative areas S1, S2 and S3 respectively. Thermal conductivity and convection are given by the heat transfer coefficients (let's call them β1 and β2). The radiation is characterized by a given emissivity E. Then the calculated heat transfer coefficient β can be calculated as:


(2)

where – radiation heat transfer coefficient.


Knowing the coefficients α and β, one can write down the equation for a heat flux density:


(3)

where .


At the outer boundaries of volumes, heat transfer in the general case can be carried out by heat transfer, convection and radiation. In addition, a constant heat flux can be specified at the outer boundary. The heat flux density in general is determined by the expression:


(4)

where q – heat flux;
          Т – temperature on the outer surface of the body;
          Тa – ambient temperature.


Using the described model, the user can specify any mechanism for transferring heat from body to body or from body to environment. The following lines are used for this:

  • Layer Heat Transfer Coef. (α, W/m2/K). This row specifies the heat transfer coefficient of the layer α, which can be used to account for the effect of the non-stick coating of the mold. It can be calculated using the formula (1). Specified by a constant value in the Value cell. If you set the parameter to 0, it will mean that layer I is missing.
    Enter a floating (real) value.
    Default value: 0.


  • Heat Transfer (β1, W/m2/K). The line specifies the heat transfer coefficient βI, which then participates in the calculation of the total coefficient according to the formula (3). This parameter can be used to define the heat transfer between two bodies (for example, a casting-mold). Set as a constant value or curve by selecting the appropriate mode. A constant coefficient is set in the Value cell, the curve can be set in the lower part of the window. In the cell Surface Fraction, [0 ... 1] it is indicated which part of the surface is involved in heat exchange: 0 - does not participate; 0.5 - 50% percent of the surface is involved; 1 - the entire surface is involved.
    Enter a floating (real) value.
    Default value: 1000 for boundary 1 and 0 for other boundaries.


  • Convecton (β2, W/m2/K). The line specifies the heat transfer coefficient βII, which then participates in the calculation of the total coefficient according to formula (3). The parameter can be used to define heat transfer at the outer boundary (for example, mold-environment).. Set as a constant value or curve by selecting the appropriate mode. A constant coefficient is set in the Value cell, the curve can be set in the lower part of the window. In the cell Surface Fraction, [0 ... 1] it is indicated which part of the surface is involved in heat exchange: 0 - does not participate; 0.5 - 50% percent of the surface is involved; 1 - the entire surface is involved.
    Enter a floating (real) value.
    Default value: 10 for boundaries 17, 18 and 0 for other boundaries.


  • Emissivity. The line specifies the emissivity E, which then participates in the calculation of the total coefficient according to formula (3). The parameter can be used to define radiation at the outer boundary (for example, mold-environment). Specified by a constant value in the Value cell. In the cell Surface Fraction, [0 ... 1] it is indicated which part of the surface is involved in heat exchange: 0 - does not participate; 0.5 - 50% percent of the surface is involved; 1 - the entire surface is involved.
    Enter a floating (real) value.
    Default value: 0.8 for boundary 17, 0.5 for boundary 18 and 0 for other boundaries.

To switch off one or two heat transfer mechanisms from the calculation, it is necessary to set the corresponding fraction of the surface equal to 0. At least one of the three mechanisms of the second layer must be enabled.

  • Heat Flux, W/m2. In some cases, it is necessary to specify the heat input or loss across the boundary as a given heat flow. The heat flow can be specified by a constant value or a time-heat flux curve. Heat input across the boundary is specified by a positive value, and heat loss by a negative heat flow value. The line is only available in the Body-Environment boundary mode.
    Enter a floating (real) value.
    Default value: 0.


  • Interface. The parameter specifies the mesh type on the contact boundary of two volimes. It is recommended to create mesh models where the meshes are joined at common boundaries with the "node-node" condition, i.e. for each mesh node 1 lying on the contact surface, there is a coincident mesh node 2 with the same coordinates. Using coincident meshes gives the most accurate result when calculating the heat flow between bodies. However, there are situations when this condition cannot be met. For example, if the bodies move relative to each other or the generation of a coincident mesh is impossible for technical reasons. In such problems, a more complex algorithm for finding contact surfaces and calculating their common area is used. The parameter takes the following values ​​(selected from the list):
    • Coincident Mesh. The most common and recommended interface type.
    • Noncoincident Mesh. Selecting this value makes available additional parameters responsible for the operation of a more complex algorithm.
    Default value: Coincident Mesh.


  • Maximum Gap. In general, it is assumed that there may be a gap between noncoinciding meshes. The parameter specifies the search area for the "adjacent" face through which the heat flow can be directed. The higher the value of the parameter, the greater the probability of finding contact between two meshes. However, you should not set the value too large, exceeding the characteristic size of the element edge at the contact boundary.
    Enter a floating (real) value in millimeters.
    Default value: 1.


  • Filter. To limit the search for contacts, you can use one of the filters from the list:
    • Casting-Mold. Search for contact only at the boundary between different volume types .
    • Casting-Casting or Mold-Mold. Search for contact only on the boundary between the same types of volumes.
    • Any-Any. Search for contact at the boundary between all types of volumes.
    Default value: Casting-Mold.


  • Heat transfer under the condition. The parameter specifies how to assign heat transfer parameters if the contact conditions have changed. In the Mode column, you can select one of two conditions:
    • Contact not found. This is the default condition. It is triggered if the bodies had contact along the specified boundary, and at some point in time this contact is lost (completely or partially) and the boundary becomes external. For example, a slab moves inside the mold (there is contact), and then part of it ends up outside the mold in the secondary cooling zone (contact with the crystallizer is partially lost). In this case, the heat transfer parameters specified for the boundary selected from the list in the Value column (by default, this is boundary 0) will be used in the contact loss area.
    • Contact found. This is an inverse condition that will start to operate when contact between bodies occurs along a given boundary (completely or partially). This condition is convenient to use, for example, when a slab passes through different cooling zones with different conditions of heat exchange with the environment, but some of its sections are in contact with rolls and the heat transfer conditions at the "slab-rolls" boundary are the same for all rollers. In this case, the heat transfer parameters specified for the boundary selected from the list in the Value column (by default, this is boundary 0) will be used in the contact occurrence area.

Data from files of outdated formats that were used in previous versions of POLIGONSOFT can be imported into the project:

  • files with afo extension with heat transfer parameters;
  • files with qgr extension with heat flux data.

Temperature Modes


Temperature Modes editor is intended for creating and editing temperature modes of the casting and mold volumes. The settings can be applied to solve one specific problem or exported to a file with the .tel extension for use in other calculations. The ".tel" temperature file can be a part of a casting process template.


The button to launch the Temperature Modes editor is located on the Process toolbar and on the Main tab of the ribbon menu. The editor can also be opened from the Start Simulation and Parameters window.

In the table that occupies the upper part of the editor (see the figure below), temperature modes are set for all volumes of the model, for which you need to set special conditions that differ from those specified in the Start Simulation and Parameters window on the General tab (see deatils in Сommon Simulation Parameters section). Regardless of the problem being solved, the table always contains data on all 9 volume indices, which in general can all be specified in the model (more details about volume indices see in Basic concepts: Volumes, Types, Indices section). To display in the table only those volume indices that are specified in the model, you can use the In Project filter located on the toolbar (see the figure below). To display only volume indices with specified temperature modes, use the With Modes filter. Use the All filter to display the entire list of volume indices.


Temperature Modes Editor

For each volume index, the Name column contains the index, its color designation and explanatory text (name of the volume). In the next column Mode there is a switch for the set temperature mode. From the drop-down list it is proposed to select:

  • [not assigned] – the temperature mode for the volume(s) is not set, the temperature for the volumes with this index will be assigned in the Start Simulation and Parameters window.


  • Initial Temperature – the initial temperature of the volume(s) is set, which is different from the one specified in the Start Simulation and Parameters window.


  • Temperature vs Time – the temperature of the volume(s) is given by a constant or a temperature-time curve. If the calculation time exceeds the time set for the last point of the curve, the temperature of the volume (s) remains equal to the temperature set for the last point of the curve.


  • Temperature vs Time with Cutoff – the temperature of the volume(s) is set by the temperature-time curve. If the calculation time exceeds the time specified for the last point of the curve, the mode is "turned off", i.e. the temperature of the volume(s) begins to be calculated as usual.

Please note that the temperature set for the volume index will be applied to all volumes with this index regardless of their type (casting or mold). Pay attention to the temperature icons in the Model Tree to avoid mistakes.

The selected temperature mode is set in the Temperature line. This line is visible only if the temperature mode is set for the volume index. In the Initial Temperature mode, the temperature can only be set by a constant in the Value column. In the Temperature vs Time with Cutoff mode, the temperature can only be set by the curve in temperature-time coordinates.

If the Temperature vs Time mode is set, in the Mode column you can select one of the options: Value or Curve. If Value is selected, the temperature of the volume(s) of the selected index is set by a constant in the Value column. This means that the temperature in the volume(s) will be constant throughout the entire calculation

If the Curve option is selected, the table and the curve drawing area become active at the bottom of the editor, where you can specify the points of the curve in temperature-time coordinates. If the calculation time exceeds the time set for the last point of the curve, the temperature of the volume(s) remains equal to the temperature set for the last point of the curve. For more information on editing graphs, see the Editing Constants and Functions section.

For the procedure for saving temperature modes and descriptions of other buttons located on the panel, see Toolbar Buttons section.


Dynamic Environment


The Dynamic Environment editor is designed to create and edit an environment consisting of several layers with changing parameters. The settings can be applied to solve one specific problem or exported to a file with the pdes extension for use in other calculations (import of older files with the air extension is also supported).


The button for launching the Dynamic Environment editor is located on the Process toolbar and on the Main tab of the ribbon menu. The editor can also be opened from the Start Simulation and Parameters window and from the tree on the Project Overview tab.

The dynamic environment is an environment that is divided vertically into 5 layers. Each layer has its own characteristics (height, temperature, and other parameters described below), some of which can change during the simulation. These settings overlap the parameters of the "normal" environment, which are set in the Start Simulation and Parameters window on the General tab (for more details, see the Сommon Simulation Parameters section) and in the Heat Transfer editor (for more details, see the Heat Transfer Parameters section). The dimensions of each layer are determined only by its height or thickness; in the other two directions, the layer is infinite.

During the calculation process, the casting and mold can be stationary in one or more layers simultaneously, can move through these layers and go beyond them. Outside the dynamic environment, normal environmental conditions apply.

In the table that occupies the upper part of the editor (see the figure below), the parameters of the layers of the dynamic environment are specified, for which it is necessary to set special conditions that differ from those specified in the "Parameters and run calculation" window on the "General" tab (for more information, see the Сommon Simulation Parameters section). Regardless of the problem being solved, the table always contains data on all 5 layers of the environment, but some layers can be switched off and thus the dynamic environment will have gaps where heat exchange will proceed according to the rules of the "normal" environment.


Dynamic Environment editor

For each layer, the Name column by default indicates its name, number and color designation (the name can be replaced with a custom one). The next column, Mode, contains the layer activity switch. From the drop-down list, you can select:

  • Enabled – the layer is active, its parameters will be used in the calculation.


  • Disabled – the layer is set, but not active, in the space limited by the layer height, the parameters of the "normal" environment are in effect.

Below are the layer parameters:

  • Layer height. Sets the height of the selected layer. The layers of the dynamic environment are counted from the bottom up from the lower point of the calculation area (mesh) in the direction of the Y axis.
    Enter a floating (real) value in millimeters.
    Default value: 100.

It should be remembered that even if the dynamic environment layer is disabled, its height matters, since it sets the gap between adjacent layers or an additional layer of "normal" environment on the edge.

  • Layer Temperature. This parameter specifies the temperature of the environment in the layer. In the Mode column, you can select one of the options: Value or Curve. If you select Value, the temperature of the environment in the selected layer is set by the constant in the Value column. This means that the temperature of the environment will be constant throughout the calculation. If you select Curve, the table in the lower part of the editor where you can set curve points in time-temperature coordinates and the curve drawing area become active. If the calculation time exceeds the time specified for the last point of the curve, the temperature of the environment in the layer remains equal to the temperature specified for the last point of the curve. For more information on editing graphs, see the Editing Constants and Functions section.
    Enter a floating (real) value in Celsius degrees.
    Default value: 20.


  • Volume Indices. The parameter defines the indices of the volumes for which the temperature specified by the parameter "Layer temperature" (see above) will be assigned at each calculation step. A text string consisting of the numbers (indices) of the volumes separated by spaces (e.g. 2 6 4) is entered in the field. If one of the specified volumes crosses the layer boundary, the temperature will be assigned only to the part of the volume that is within the range of this parameter.
    Default value: empty string.


  • Heat Transfer Coefficient. The parameter specifies the conditions of heat exchange with the environment in the layer. The heat transfer parameters specified by the boundary of the computational domain or the parameters of any other boundary can be used. This allows changing the heat exchange conditions when the domain enters the layer of the dynamic environment. For example, immersion of the mold in a liquid metal cooler can be simulated in this way. In the Mode column, you can select one of the options:
    • As on the domain boundary. In this case, the heat transfer parameters for the boundary indices specified in the calculation domain will be used in the layer area.
    • As specified for the boundary. In this case, the heat transfer parameters set for the boundary selected from the list in the Value column will be used in the layer area. These parameters will be applied to all boundaries of the calculation area (or their parts) that fall within the layer.

For information on the procedure for importing/exporting a dynamic environment and a description of other buttons located on the toolbar, see the Toolbar Buttons section.


Translation Modes


When modeling some foundry processes, it is necessary to take into account the movement of casting and/or mold volumes. For example, when simulating vacuum casting, you may need to take into account the movement of the mold within the furnace. To create and edit volume movements, use the Translation editor and the Translations tree in the Project Overview panel. Translation modes can be configured and used to solve one specific problem, or imported (and exported) into the project from files with the pmov extension (importing files of older versions with the mve extension is also supported).


The Translation editor button is located on the Process toolbar and on the Main tab of the ribbon menu. The editor can also be opened from the Translation tree on the Project Overview tab.

The editor title indicates the volume index the translation mode will be created for, for example, in Fig. Below is the editor for moving volumes with index 2. Please keep in mind that the specified mode will be applied to all volumes with index 2, regardless of their type (casting or mold). If the translation mode has not yet been specified, then all interface elements are inactive. To start working in the editor, you need to create a translation mode for the selected index - click the New button.

Unlike other editors, there are no Import and Export commands on the toolbar, since these operations are performed through the Translations tree (for more details, see the Project Overview Tab section). Otherwise, the principles of operation follow the general rules described in Data Editors. General Controls and Features and Editing Constants and Functions sections.

The translation mode for the selected volume index is set by the table at the bottom of the window (see the figure below). For each line, the first column Time indicates the moment in time, the next three columns indicate the shift in x, y, z coordinates relative to the previous position described in the previous line. The fourth column “L” is not editable, it shows the shift value, which is calculated using the equation:


,

where Xt, Yt, Zt – shift coordinates at time t;
           Xt-1, Yt-1, Zt-1 – shift coordinates at time t-1;
           t > t-1.


The L(t) graph is displayed to the right of the coordinate table.


Translation mode for volumes with index 2

Unlike most other editors, here you can set a single-point movement schedule. In this case, bodies with the specified index will be instantly moved in the specified direction. This can be useful for modeling processes that do not require the action to be performed in time, such as opening a mold in high-pressure casting.

Often, translation is used when simulating vacuum casting in combination with radiation heat transfer (for more details, see the Radiation Heat Transfer section). Often, translation is used when simulating vacuum casting in combination with radiation heat transfer (for more details, see the “Radiation Heat Transfer” section). If translations are specified, a ray tracing calculation will be performed for each point in the graph, and linear interpolation between two adjacent times for all intermediate points in time will be used. For such calculations, you should specify the number of points necessary for a qualitative description of the visibility of objects during the entire calculation. But it should be remembered that a large number of points in the translation curve significantly slows down the calculation with radiation heat transfer.

All specified translation are displayed in the Translations tree on the Project Overview panel. They can be exported to a .pmov (for more details, see the Project Overview Tab section). This file can then be imported into another project (importing files of older versions with the mve extension is also supported).


Surface Nucleation


For some castings, requirements may be imposed on the quality of the macrostructure: size, shape, quantity and orientation of grains. POLIGONSOFT allows you to model the processes associated with the nucleation and growth of the grain structure of the metal (for more details, see the Macrostructure Simulation section). To create and edit boundary conditions describing the process of grain nucleation at the casting-mold boundaries, use the Surface Nucleation editor and the tree of the same name in the Project Overview panel. Nucleation conditions can be configured and used to solve one specific problem, or imported (and exported) into the project from files with the psn extension.


The button to launch the Surface Nucleation editor is located on the Process toolbar and on the Main tab of the ribbon menu. The editor can also be opened from the Surface Inception tree on the Project Overview tab.

The editor header indicates the index of the boundary for which the generation parameters will be created, for example, in the figure below the editor is open for the boundary with index 1. If the parameters have not yet been set, the editor looks like in the figure below, all interface elements are inactive. To start working in the editor, you need to create generation parameters for the selected boundary - click the Create button.

Unlike other editors, there are no Import and Export commands on the toolbar, since these operations are performed via the Surface Nucleation tree (for more details, see the Project Overview Tab section). Otherwise, the operating principles follow the general rules described in the Data Editors. General Controls and Features and Editing Constants and Functions sections.

The parameters for the grain structure generation for the selected boundary index are specified in the table (see the figure below). In the first line of the editor, you can enter a name for the data set being created. By default, the boundary identifier is used.


Nucleation parameters for boundary with index 1

The appearance of solid phase nuclei in a supercooled melt is a random process. The number of nuclei that appear in the melt when the supercooling changes by a small dT) value is described by the Gaussian distribution (see equation (3) in the Grain Nucleation Model section). The nucleation of the solid phase for the casting-mold boundaries is determined by the following parameters:

  • Structure Type. Specifies the type of structure that will be generated on the surface:
    • Equiaxed. Suitable for most cases, since it is difficult to know in advance what structure will begin to form on the surface, this is the purpose of the calculation. When choosing this type of structure, grains of arbitrary orientation will form on the surface.
    • Texture. It is assumed that grains of a given orientation will nucleate on the surface. It is suitable for simulating the interaction of a melt with a seed having a polycrystalline structure with a limited range of grain orientation angles.
    • Single Crystal. It is assumed that one or more grains of a given orientation will be generated on the surface. Suitable for simulating the interaction of the melt with a seed having a single-crystal structure with a given orientation.
    Default value: Equiaxed.


  • Orientation. When using the 3D solver, it is specified by three Euler angles a, b, g for all types of structure except equiaxed. When using the 3D solver, it is specified by three Euler angles a, b, g for all types of structure except equiaxed. For a textured structure, the orientation of the nucleated grains is specified by the ranges of values ​​of each angle. To do this, a pair of numbers separated by a space is entered in each line. For a single-crystal structure, the orientation is specified by three Euler angles; entering angle ranges is not allowed. When using the 2D solver, only the angle a is specified.
    Specified as integers in the range from -45 to 44 degrees.


  • Maximum Density of Nuclei. The parameter (corresponds to nmax in equation (3)) sets the density of solid phase nuclei formed in the melt upon reaching undercooling many times greater than DTmax (see equation (3)). The choice of the parameter value depends on the choice of the 2d or 3d version of the macrostructure solver (see equation (4) in the Grain Nucleation Model section).
    Specified in m-2 or m-1 if 3d or 2d solver is selected respectively.
    Default value: 3000 if 3d solver is specified.


  • Average Undercooling. The parameter determines the position of the maximum of the Gaussian curve (corresponds to DTmax in equation (3)).
    Specified in degrees.
    Default value: 0.1.


  • Standard Deviation. The parameter defines the peak width of the Gaussian curve (corresponding to DTs in equation (3)).
    Specified in degrees.
    Default value: 0.1.


  • Integration Range Parameter. The parameter sets the melt undercooling range, at which the nucleation of the solid phase can occur. Outside this range, the appearance of solid phase nuclei is impossible. Thus, the subcooling range is a multiple of the standard deviation DTs.
    Default value: 5.

Structure Inheritance


Macrostructure simulation (for more details, see the Macrostructure Simulation section) can be performed either in the entire casting or in a certain specified area, which is set in the Start Simulation and Parameters window on the Structure tab (for more details, see the Macrostructure Simulation Parameters section). To obtain results corresponding to the specified technology, it may be necessary to set the conditions for inheriting the structure on the faces of the parallelepiped that limits the calculation area.

For example, the figure below shows the structure modeling area, in which the structure will be generated on the side faces in accordance with the surface nucleation conditions at the casting-mold boundary. In this case, the upper and lower faces of the area border on the melt, where the nucleation conditions cannot be specified, since there is no boundary. To take into account the influence of the structure formed outside the calculation area, it is possible to apply the mechanism of its inheritance from neighboring areas of the casting.

To do this, use the Structure Inheritance editor and the tree of the same name in the Project Overview panel. The inheritance conditions can be configured and used within the framework of solving one specific problem, or imported (and exported) into the project from files with the pinh extension.


Example of computational domain
for macrostructure simulation


The button for launching the Structure Inheritance editor is located on the Process toolbar and on the Main tab of the ribbon menu. The editor can also be opened from the Structure Inheritance tree on the Project Overview tab.

The Structure Inheritance tree contains 6 lines corresponding to the six faces of the parallelepiped that limits the calculation area. Since the faces of the parallelepiped are parallel to the coordinate planes, they are designated by icons with filled faces. The face color corresponds to the color of the axis perpendicular to which the face is located. When calling the editor, the face is identified in its header by the axis and position relative to the origin: Xmin, Xmax, Ymin, Ymax, Zmin and Zmax. If the parameters have not yet been set, all interface elements are inactive. To start working in the editor, you need to create inheritance parameters for the selected face - click the Create button.

Unlike other editors, the Import and Export commands, as well as filters, are inactive here, since these operations are performed via the Surface Nucleation tree (for more details, see the Project Overview Tab section). Otherwise, the operating principles follow the general rules described in the Data Editors. General Controls and Features and Editing Constants and Functions sections.

The parameters of grain structure inheritance for the selected face of the parallelepiped that limits the modeling area are specified in the table (see the figure below). In the first line of the editor, you can enter a name for the data set being created. By default, the plane identifier is used (for example, Xmin Plane in the figure below).


Inheritance parameters for the Xmin plane

The appearance of solid phase nuclei in a supercooled melt is a random process. The number of nuclei that appear in the melt when the supercooling changes by a small dT) value is described by the Gaussian distribution (see equation (3) in the Grain Nucleation Model section). The nucleation of the solid phase for the casting-mold boundaries is determined by the following parameters:

  • Structure Type. Specifies the type of structure that will be generated on the face:
    • Equiaxed. Suitable for casting cases where the structure formation is not controlled by the technological process. When choosing this type of structure, grains of arbitrary orientation will be formed on the face of the calculation domain.
    • Texture. It is assumed that grains of a given orientation will be generated on the boundary of the computational domain. Suitable for simulating the interaction of the melt with the region where a structure with a limited range of grain orientation angles has already formed.
    • Монокристалл. It is assumed that one or more grains of a given orientation will be generated on the boundary of the computational domain. Suitable for simulating the interaction of the melt with the region where a single-crystal structure with a given orientation has already formed.
    Default value: Equiaxed.


  • Orientation. When using the 3D solver, it is specified by three Euler angles a, b, g for all types of structure except equiaxed. When using the 3D solver, it is specified by three Euler angles a, b, g for all types of structure except equiaxed. For a textured structure, the orientation of the nucleated grains is specified by the ranges of values ​​of each angle. To do this, a pair of numbers separated by a space is entered in each line. For a single-crystal structure, the orientation is specified by three Euler angles; entering angle ranges is not allowed. When using the 2D solver, only the angle a is specified.
    Specified as integers in the range from -45 to 44 degrees.


  • Maximum Density of Nuclei. The parameter (corresponds to nmax in equation (3)) sets the density of solid phase nuclei formed in the melt upon reaching undercooling many times greater than DTmax (see equation (3)). The choice of the parameter value depends on the choice of the 2d or 3d version of the macrostructure solver (see equation (4) in the Grain Nucleation Model section).
    Specified in m-2 or m-1 if 3d or 2d solver is selected respectively.
    Default value: 3000 if 3d solver is specified.