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:

Testing status by file:

• cmiss_test.log
  

  i686-linux
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  mips-irix
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  Success	cm-debug-clear-malloc:	cmiss_test.log.retain.
  Success	cm-debug-clear-malloc7:	cmiss_test.log.retain.
  Success	cm64:	cmiss_test.log.retain.
  Success	cm64-debug:	cmiss_test.log.retain.
  Success	cm64-debug-clear-malloc:	cmiss_test.log.retain.
  rs6000-aix
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  Success	cm64:	cmiss_test.log.retain.
  Success	cm64-debug:	cmiss_test.log.retain.
  x86_64-linux
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.

  
  
• cmiss_test.out
  

  i686-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  mips-irix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
  rs6000-aix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  x86_64-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.

  
  

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

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