https://nemesis-project.org/index.php?action=history&feed=atom&title=Example%3ATransient_03 2026-04-21T17:04:12Z Revision history for this page on the wiki MediaWiki 1.19.24 https://nemesis-project.org/index.php?title=Example:Transient_03&diff=1874&oldid=prev 2007-01-31T18:43:27Z

<p></p> <p><b>New page</b></p><div>&lt;pre&gt;<br /> from numpy import array,arange<br /> from matplotlib import use<br /> from matplotlib.pylab import plot,xlabel,ylabel,title,grid,savefig<br /> <br /> steps=1000<br /> timestep=0.01<br /> for k in range(0,8):<br /> <br /> domain.dim(1)<br /> material.uniElastic(1,1000.,0.2,10.)<br /> section.user(1,1.)<br /> <br /> node.add(1,0.)<br /> node.add(2,1.)<br /> <br /> element.bar(1,1,2,1,1)<br /> node.fix(1,1)<br /> <br /> lc.define(1)<br /> load.sin(2,1,10.,150.,0.)<br /> <br /> tracker.node(1,2)<br /> domain.RayleighDamping(0.0001*k,0.25*k)<br /> analysis.transient()<br /> imposer.elimination()<br /> algorithm.linear()<br /> control.Newmark(0.25,0.5,timestep)<br /> soe.full()<br /> analysis.run(1,steps)<br /> <br /> t=arange(0,(steps+1)*timestep,timestep)<br /> s=[]<br /> for i in range(tracker.steps(1)):<br /> s.append(tracker.data(1,i)[ 5])<br /> s=array(s) <br /> domain.clear()<br /> plot(t, s, linewidth=1.0)<br /> xlabel('time')<br /> ylabel('u')<br /> title('Rayleight damping')<br /> grid(True)<br /> savefig('RayleighDamping.png')<br /> &lt;/pre&gt;<br /> <br /> [[Category:Examples]]</div>

Fotis