Impact simulation
Hi there,
Please I am trying to impact a steel ball on the surface area of a cylindrical rod. I want to get the force-time history during the impact and also need an output file to view in paraview. Please I need help with the following:
1) Specification of material micro properties such as young modulus and poisson's ratio for the steel-ball and steel rod . The impactor is a steel and my rod is also a steel
2)Inserting the cylindrical (0.02m diamter and 1.5m length) rod into the domain
3)Basis of getting the normalK based on the material properties.
i was able to modify part of the tutorial code to suit my work. Here is the code and response that I am getting after attempting to run it.
#Import the appropriate Esys-Particle module
from esys.lsm import*
from esys.lsm.util import Vec3, BoundingBox
from esys.lsm.geometry import*
#instantiate a simulation object
#number of processor for simulation
sim = LsmMpi (numWorkerProce
#initialise the neighbour search algorithm:
sim.initNeighbo
particleType=
gridSpacing=
verletDist=0.005 #i.e, 0.2 Xply by min. particle radius which is 0.025m
)
#setting number of timestep and timestep size
sim.setNumTimeS
sim.setTimeStep
#Domain of simulation: here is the cube with 40m each size
domain = BoundingBox(
sim.setSpatialD
#create a steel ball with NRotSphere type, id=0 position at that vector posn, with the specified radius and mass
steel_ball=
#initialise the velocity of the steel_ball Vx, Vy and Vz (m/s)
steel_ball.
#add the newly created steel_ball
sim.createParti
#initialise gravity in the domain:
sim.createInter
GravityPrms(
)
#add a horizontal wall to act like the steel rod surface:
sim.createWall(
name="steel_rod",
posn=Vec3(0,0,0),
normal=Vec3(0,1,0)
)
#specify the type of interaction between wall and steel ball
sim.createInter
NRotElasticWal
name=
wallName=
normalK=10000.0
)
)
#add local viscosity to simulate air resistance
sim.createInter
LinDampingPrms(
name=
viscosity=0.1,
maxIterations=100
)
)
#create a FieldSaver to store forces on the steelrod:
force_saver=
wallName=
fieldName="Force",
fileName=
YbeginTimeStep=0,
endTimeStep=100,
timeStepIncr=0.001
)
sim.createField
#add check pointer to store simulation data for paraview
sim.createCheck
CheckPointPrms (
fileNamePrefix ="snapshot",
beginTimeStep =0,
endTime=1000,
timeStepIncr=4
)
)
#run the simulation
sim.run()
adminuser@
adminuser@
adminuser@
adminuser@
adminuser@
adminuser@
CSubLatticeCont
slave started at local/global rank 0 / 1
Traceback (most recent call last):
File "calibration.py", line 71, in <module>
timeStepIncr=4
Boost.Python.
CheckPointP
did not match C++ signature:
__init_
Kind Regards
Tope
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Temitope
- Solved:
- Last query:
- Last reply: