Example 235: Gauss point fitting to obtain a stress field


The comfile run by this example is as follows:

# Example_235 Fit a nodal strain energy field from Gauss point values 
#             and plot
                                   
fem                            #Set the FEM environment   
fem define node;r;gaussfit;example #Define building nodes
fem define base;r;;example         #Read in a bilinear basis function
fem define element;r;;example      #Define building elements
fem define fibre;d                 #Default fibre angle
fem define element;r;;example fibre
fem define window;c            #Calculate the window size
fem draw line                  #Draw the mesh
fem define equation;r;;example lock #Define a linear elasticity problem with
#                              !  Plane stress elements
fem define material;r;;example     #Define steel material properties
fem define initial;r;;example      #Define the loading on the building
fem draw node                  #Draw the nodes and show which nodes
#                              !  are fixed.
fem draw reaction sc 0.01      #Draw the initial reactions on the 
#                              !  building
fem define solve;r;;example        #Define the solve parameters
fem solve                      #Solve the linear elastic problem
fem draw line deformed dotted scale 0.01 #Draw the deformed lines
fem hide node                  #Hide the nodes
fem hide reaction              #Hide the reactions
fem draw gauss                 #Draw the gauss points
fem update gauss stress principal #Puts 1st princ stress 
#                              !  into YG(ng,1,ne)
#                              !  puts 2nd princ stress into YG(ng,2,ne)
#                              !  puts principal  angle into YG(ng,3,ne)
#                              !  puts strain energy    into YG(ng,4,ne)
fem define field;r;;example    #Define a field to hold the fitted 
#                              !  solution
fem define element;r;;example field #Define field elements
fem define fit;r;;example gauss lock #Define a gauss point fit
fem list fit                   #List the fit
fem list node                  #List the initial field
fem fit                        #Fit the first principal stress into
#                              !  field 1 and the second principal
#                              !  stress into field 2.
fem list node                  #List the fitted fields
fem draw field field 1         #Draw the first principal stress field
fem draw field field 2         #Draw the second principal stress field

Additional testing commands:

fem list node

Files used by this example are:

Name             Modified     Size

example_235.com 10-Apr-2000 2.4k gaussfit.ipbase 10-Apr-2000 1.1k gaussfit.ipelem 10-Apr-2000 6.7k gaussfit.ipelfb 10-Apr-2000 3.7k gaussfit.ipelfd 10-Apr-2000 5.1k gaussfit.ipequa 26-May-2003 1.7k gaussfit.ipfiel 10-Apr-2000 7.9k gaussfit.ipfit 13-Apr-2007 1.7k gaussfit.ipinit 10-Apr-2000 757 gaussfit.ipmate 10-Apr-2000 1.4k gaussfit.ipnode 10-Apr-2000 6.8k gaussfit.ipsolv 13-Apr-2007 1.1k test_output.com 10-Apr-2000 14

Download the entire example:

Name                      Modified     Size

examples_2_23_235.tar.gz 14-Apr-2007 8.0k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:40 20161
cm-debugSuccessSat Mar 5 00:01:57 20161
mips-irix
cmSuccessSun Aug 19 01:33:49 20077
cm-debugSuccessWed Aug 15 01:30:10 200711
cm-debug-clear-mallocSuccessSat Aug 18 01:33:52 200717
cm-debug-clear-malloc7SuccessMon Aug 20 01:32:09 200716
cm64SuccessSun Aug 19 01:34:20 20077
cm64-debugSuccessTue Aug 21 01:26:56 200712
cm64-debug-clear-mallocSuccessThu Apr 1 10:12:04 20045
rs6000-aix
cmSuccessWed Mar 4 01:08:00 20091
cm-debugSuccessMon Mar 2 01:07:13 20092
cm64SuccessWed Mar 4 01:08:00 20091
cm64-debugSuccessTue Mar 3 01:12:05 20092
x86_64-linux
cmSuccessSun Mar 6 00:01:01 20161
cm-debugSuccessSat Mar 5 00:01:12 20161

Testing status by file:


Html last generated: Sun Mar 6 05:50:12 2016

Input last modified: Fri Apr 13 10:31:22 2007


CMISS Help / Examples / 2 / 23 / 235