In this section, a cantilever beam loaded by point forces at its free end is analyzed and a thermal calculation of a furnace is performed.
The geometry, loading and boundary conditions of the cantilever beam are shown in Figure
1. The size of the beam is 1x1x8 , the loading consists of a
point force of
N and the beam is completely fixed (in all directions) on the left end. Let us take 1 m and 1 MN as units of length and force, respectively. Assume that the beam geometry was generated and meshed with CalculiX GraphiX (cgx) resulting in the mesh in Figure 2. For reasons of clarity, only element labels are displayed.
A CalculiX input deck basically consists of a model definition section describing the geometry and boundary conditions of the problem and one or more steps (Figure 3) defining the loads.
The model definition section starts at the beginning of the file and ends at the occurrence of the first *STEP card. All input is preceded by keyword cards, which all start with an asterisk (*), indicating the kind of data which follows. *STEP is such a keyword card. Most keyword cards are either model definition cards (i.e. they can only occur before the first *STEP card) or step cards (i.e. they can only occur between *STEP and *END STEP cards). A few can be both.
In our example (Figure 4), the first keyword card is *HEADING, followed by a short description of the problem. This has no effect on the output and only serves for identification. Then, the coordinates are given as triplets preceded by the *NODE keyword. Notice that data on the same line are separated by comma's and must not exceed a record length of 132 columns. A keyword card can be repeated as often as needed. For instance, each node could have been preceded by its own *NODE keyword card.
Next, the topology is defined by use of the keyword card *ELEMENT. Defining the topology means listing for each element its type, which nodes belong to the element and in what order. The element type is a parameter on the keyword card. In the beam case 20-node brick elements with reduced integration have been used, abbreviated as C3D20R. In addition, by adding ELSET=Eall, all elements following the *ELEMENT card are stored in set Eall. This set will be later referred to in the material definition. Now, each element is listed followed by the 20 node numbers defining it. With *NODE and *ELEMENT, the core of the geometry description is finished. Remaining model definition items are geometric boundary conditions and the material description.
The only geometric boundary condition in the beam problem is the fixation at z=0. To this end, the nodes at z=0 are collected and stored in node set FIX defined by the keyword card *NSET. The nodes belonging to the set follow on the lines underneath the keyword card. By means of the card *BOUNDARY, the nodes belonging to set FIX are subsequently fixed in 1, 2 and 3-direction, corresponding to x,y and z. The three *BOUNDARY statements in Figure 4 can actually be grouped yielding:
*BOUNDARY FIX,1 FIX,2 FIX,3
or even shorter:
*BOUNDARY FIX,1,3
meaning that degrees of freedom 1 through 3 are to be fixed (i.e. set to zero).
The next section in the input deck is the material description. This section is special since the cards describing one and the same material must be grouped together, although the section itself can occur anywhere before the first *STEP card. A material section is always started by a *MATERIAL card defining the name of the material by means of the parameter NAME. Depending on the kind of material several keyword cards can follow. Here, the material is linear elastic, characterized by a Young's modulus of 210,000.0 and a Poisson coefficient of 0.3 (steel). These properties are stored beneath the *ELASTIC keyword card, which here concludes the material definition. Next, the material is assigned to the element set Eall by means of the keyword card *SOLID SECTION.
Finally, the last card in the model definition section defines a node set LOAD which will be needed to define the load. The card starting with two asterisks in between the model definition section and the first step section is a comment line. A comment line can be introduced at any place. It is completely ignored by CalculiX and serves for input deck clarity only.
In the present problem, only one step is needed. A step always starts with a *STEP card and concludes with a *END STEP card. The keyword card *STATIC defines the procedure. The *STATIC card indicates that the load is applied in a quasi-static way, i.e. so slow that mass inertia does not play a role. Other procedures are *FREQUENCY, *BUCKLE, *MODAL DYNAMIC and *DYNAMIC. Next, the concentrated load is applied (keyword *CLOAD) to node set LOAD. The forces act in y-direction and their magnitude is 1, yielding a total load of 9.
Finally, the printing and file storage cards allow for user-directed output generation. The print cards (*NODE PRINT and *EL PRINT) lead to an ASCII file with extension .dat. If they are not selected, no .dat file is generated. The *NODE PRINT and *EL PRINT cards must be followed by the node and element sets for which output is required, respectively. Element information is stored at the integration points.
The *NODE FILE and *EL FILE cards, on the other hand, govern the output written to an ASCII file with extension .frd. The results in this file can be viewed with CalculiX GraphiX (cgx). Quantities selected by the *NODE FILE and *EL FILE cards are always stored for the complete model. Element quantities are extrapolated to the nodes, and all contributions in the same node are averaged. Selection of fields for the *NODE PRINT, *EL PRINT, *NODE FILE and *EL FILE cards is made by character codes: for instance, U are the displacements and S are the (Cauchy) stresses.
The input deck is concluded with an *END STEP card.
The output files for the beam problem consist of file beam.dat and beam.frd. The beam.dat file contains the displacements for set Nall and the stresses in the integration points for set Eall. The file beam.frd contains the displacements and extrapolated stresses in all nodes. It is the input for the visualisation program CalculiX GraphiX (cgx). An impression of the capabilities of cgx can be obtained by looking at Figures 5, 6 and 7.
Figure 5 shows the deformation of the beam under the prevailing loads. As expected, the beam bends due to the lateral force at its end. Figure 6 shows the normal stress in axial direction. Due to the bending moment one obtains a nearly linear distribution across the height of the beam. Finally, Figure 7 shows the Von Mises stress in the beam.
The second problem involves a thermal calculation of the furnace depicted in
Figure 8. The furnace consists of a bottom plate at a temperature
, which is prescribed. It changes linearly in an extremely short time from
300 K to 1000 K after which it remains constant. The side walls of the furnace
are isolated from the outer world, but exchange heat through radiation with
the other walls of the furnace. The emissivity of the side walls and bottom
is
. The top of the furnace exchanges heat through radiation with
the other walls and with the environmental temperature which is fixed at 300
K. The emissivity of the top is
. Furthermore, the top exhanges
heat through convection with a fluid (air) moving at the constant rate of
0.001 kg/s. The temperature of the fluid at the right upper corner is 300
K. The walls of the oven are made of 10 cm steel. The material
constants for steels are: heat conductivity
,
specific
heat
and density
. The material
constants for
air are : specific heat
and density
. The
convection coefficient is
. The dimensions of the furnace
are
(cube). At
all parts are at
. We would like to know the temperature at locations A,B,C,D and E as a
function of time.
The input deck is listed in Figure 9. It starts with the node definitions. The highest node number in the structure is 602. The nodes 603 up to 608 are gas nodes, i.e. in the gas extra nodes were defined (z=0.3 corresponds with the top of the furnace, z=0 with the bottom). Gas node 603 corresponds to the location where the gas temperature is 300 K (``inlet''), node 608 corresponds to the ``outlet'', the other nodes are located in between. The coordinates of the gas nodes actually do not enter the calculations. Only the convective definitions with the keyword *FILM govern the exchange between furnace and fluid. With the *ELEMENT card the 6-node shell elements making up the furnace walls are defined. Furthermore, the gas nodes are also assigned each to one element (element type D), a so-called gas element. These elements are needed for the assignment of material properties to the gas. Indeed, traditionally material properties are assigned to elements and not to nodes. The gas nodes 603 up to 608 are assigned to the gas elements 301 up to 306.
Next, two node sets are defined: GAS contains all gas nodes, Ndown contains all nodes on the bottom of the furnace.
The *PHYSICAL CONSTANTS card is needed in those analyses in which radiation plays a role. It defines absolute zero, here 0 since we work in Kelvin, and the Stefan Boltzmann constant. In the present input deck SI units are used throughout.
Next, the material constants for STEEL are defined. For thermal analyses the conductivity, specific heat and density must be defined. The *SHELL SECTION card assigns the STEEL material to the element set FURNACE, defined by the *ELEMENT statement before. It contains all elements belonging to the furnace. Furthermore, a thickness of 0.01 m is assigned.
The material constants for material GAS consist of the density and the specific heat. These are the constants for the fluid. Conduction in the fluid is not considered. The material GAS is assigned to element set EGAS containing all gas elements.
The *INITIAL CONDITIONS card defines an initial temperature of 300 K for all nodes, i.e. furnace nodes AND gas nodes. The *AMPLITUDE card defines a ramp function starting at 0.3 at 0.0 and increasing linearly to 1.0 at 1.0. It will be used to define the temperature boundary conditions at the bottom of the furnace. This ends the model definition.
The first step describes the linear increase of the temperature boundary
condition between and
. The INC=100 parameter on the *STEP card
allows for 100 increments in this step. The procedure is *HEAT TRANSFER,
i.e. we would like to perform a purely thermal analysis: the only unknowns are
the temperature and there are no mechanical unknowns (e.g. displacements). The
step time is 1., the initial increment size is 0.1. Both appear on the line
underneath the *HEAT TRANSFER card. The absense of the parameter STEADY STATE
on the *HEAT TRANSFER card indicates that this is a transient analysis.
Next come the temperature boundary conditions: the bottom plate of the furnace is kept at 1000 K, but is modulated by amplitude A1. The result is that the temperature boundary condition starts at 0.3 x 1000 = 300K and increases linearly to reach 1000 K at t=1 s. The second boundary conditions specifies that the temperature of (gas) node 603 is kept at 300 K. This is the inlet temperature. Notice that ``11'' is the temperature degree of freedom.
The mass flow rate in the fluid is defined with the *MASS FLOW RATE card. The first line tells us that the mass flow rate between (gas)node 603 and (gas)node 604 is 0.001. Since this rate is positive the gas flows from node 603 towards node 604. The user must assure conservation of mass (this is actually also checked by the program).
The first set of radiation boundary conditions specifies that the top face of the bottom of the furnace radiates through cavity radiation with an emissivity of 1 and an environment temperature of 1000 K. For cavity radiation the environment temperature is used in case the view factor at some location does not amount to 1. What is short of 1 radiates towards the environment. The first number in each line is the element, the number in the label (the second entry in each line) is the face of the element exposed to radiation. In general, these lines are generated automatically in cgx (CalculiX GraphiX).
The second and third block define the internal cavity radiation in the furnace for the top and the sides. The fourth block defines the radiation of the top face of the top plate of the furnace towards the environment, which is kept at 300 K. The emissivity of the top plate is 0.8.
Next come the film conditions. Forced convection is defined for the top face
of the top plate of the furnace with a convection coefficient
. The
first line underneath the *FILM keyword indicates that the second face of
element 51 interacts through forced convection with (gas)node 604. The last
entry in this line is the convection coefficient. So for each face interacting
with the fluid an appropriate gas node must be specified with which the
interaction takes place.
Finally, the *NODE FILE card makes sure that the temperature is stored in the .frd file and the *NODE PRINT card takes care that the gas temperature is stored in the .dat file.
The complete input deck is part of the test examples of CalculiX (furnace.inp). For the present analysis a second step was appended keeping the bottom temperature constant for an additional 3000 seconds.
What happens during the calculation? The walls and top of the furnace heat up
due to conduction in the walls and radiation from the bottom. However, the top
of the furnace also loses heat through radiation with the environment and
convection with the gas. Due to the interaction with the gas the temperature
is asymmetric: at the inlet the gas is cool and the furnace will lose more
heat than at the outlet, where the temperature of the gas is higher and the
temperature difference with the furnace is smaller. So due to convection we
expect a temperature increase from inlet to outlet. Due to conduction we
expect a temperature minimum in the middle of the top. Both effects are
superimposed. The temperature distribution at
is shown in
Figure 10. There is a temperature gradient from the bottom of the
furnace towards the top. At the top the temperature is indeed not
symmetric. This is also shown in Figure 11, where the temperature of
locations A, B, C, D and E is plotted as a function of time.
Notice that steady state conditions have not been reached yet. Also note that 2-D elements (such as shell elements) are automatically expanded into 3-D elements with the right thickness. Therefore, the pictures, which were plotted from within CalculiX GraphiX, show 3-D elements.