Can you elaborate a little bit on that, more precisely, your modulus of elasticity will be variable with respect to what ?
To load your model with a static sinusoidal shape, the simplest way should be to store the nodes to be loaded in an array, then use their coordinates (using NX, NY and NZ) to build the load shape.
! Here select your nodes, and
! the number of the lowest
last_node=lowest_node_number
*DO,I,1,NbNodes,
curr_node=NDNEXT(last_node)
curr_x=nx(curr_node)
curr_y=ny(curr_node)
curr_z=nz(curr_node)
f,curr_node,fx,COS(2*PI*curr_x/lambda_x)
f,curr_node,fy,COS(2*PI*curr_y/lambda_y)
f,curr_node,fy,COS(2*PI*curr_z/lambda_z)
last_node=curr_node
*ENDDO
See at the end of the *GET function, if you don't know these commands.
Good Luck
Nicolas JOBERT