gen zone brick size 10 5 5
mod el
prop shear 1e8 bulk 2e8
fix x y z ran x 0 y 0 5 z 0 5
fix x y z ran x 0 10 y 0 5 z 0
fix x y z ran x 10 y 0 5 z 0 5
table 1
def find_add
head = null
pnt = gp_head
loop while pnt # null
x_pos = gp_xpos(pnt)
if x_pos = width then
new = get_mem(2)
mem(new) = head
mem(new+1) = pnt
head = new
endif
pnt = gp_next(pnt)
endloop
end
set width = 10.0
find_add
def app_vel
while_stepping
ad = head
loop while ad # null
pnt = mem(ad+1)
gp_xvel(pnt) = vel_max * gp_zpos(pnt) / height
gp_zvel(pnt) = -vel_max * (gp_xpos(pnt) - width) / height
ad = mem(ad)
endloop
end
set large vel_max = 1e-2 height = 5.0
step 100
plot grid black vel red
上面的head=new 是什么意思那？起到什么作用那？请哪位大侠指导一下？
