In CalculiX, how can we get the total force carried by some nodes representing for instance a section through the body? *NODE PRINT and RF give only the reaction forces on nodes, as I understand it, only at restrained nodes, in the restrained directions.
Below is an image showing my intent. I also provided an example case with beam elements. So, I want to find the force carried at a section through the body (a minimum section for instance, or a certain location of interest).
# ## Geometry # ## Points pnt p1 0 0 0 pnt p2 0.000000 160.000000 0 ## Lines line L1 p1 p2 40 # Line set definitions. seta Lset1 l L1 elty Lset1 be3r mesh all # All nodes are stored in a set. seta nodes n all # Save nodes of minimum neck radiusd lines. # Boundary condition nodes are selected from the all nodes and stored in sets. enq nodes fix rec 0 0 0 0.001 a enq nodes tip rec 0 160 _ 0.001 a enq nodes mid rec 0 80 _ 0.001 a # Node sets containing boundary condition nodes are written to files. send fix abq nam send tip abq nam send mid abq nam # Export nodal information to a file. send nodes abq # Export elemental information to a file. seta belems e Lset1 send belems abq seta fs p p1 seta fse n fs comp fse u send fse abq nam rot -z # Save image of the geometry with BCs and contact surfaces. ulin Geometry with boundary conditions and contact surfaces. plot l all plot m all view edge off view elem all plus n fix r 6 plus n tip b 6 plus n mid t 6 frame hcpy png geom
*INCLUDE, INPUT=nodes.msh *INCLUDE, INPUT=belems.msh *INCLUDE, INPUT=fix.nam *INCLUDE, INPUT=tip.nam *INCLUDE, INPUT=mid.nam *INCLUDE, INPUT=fse.nam *BOUNDARY Nfix,1,3 *MATERIAL, NAME=Mat *ELASTIC 185000.000000,0.300000 *DENSITY 0.000000008220000000000000 *BEAM SECTION,MATERIAL=Mat,ELSET=Ebelems,SECTION=RECT 1.,1. *STEP, NLGEOM, INC=200 *STATIC *CLOAD Ntip, 2, 100 **DLOAD **Ebelems, CENTRIF, 274155.677808, 0., 0., 0., 0., 0., 1. *EL FILE, ELSET=Ebelems U,S *NODE PRINT, NSET=Nfix RF *NODE PRINT, NSET=Nmid RF **EL PRINT, ELSET=Efse **S *END STEP