from abaqus import *
from abaqusConstants import *
import part
import material
import section
import assembly
import step
import load
import regionToolset
import mesh
import job
import visualization

#...................model.......................
mymodel=[]
for i in range(4):
    m=mdb.Model(name="model{}".format(i+1))
    mymodel.append(m)
    mysketch=mymodel[i].ConstrainedSketch(name='sketch', sheetSize=200)
    mymodel[i].rectangle(point1=(-20,-5), point2=(20,5))
    mypart=mymodel[i].Part(name='Part', dimensionality=THREE_D, type=DEFORMABLE_BODY)
    mypart.BaseSolidExtrude(sketch=mysketch, depth=100)
    mymodel.append(m)

E=[210,150,70,40]
Nu=[0.3,0.2,0.15,0.1]
Dens=[7800,6500,4600,2200]
name=[Steel,Aluminum,Nickle,Alloy]
for i in range(4):
    mat_name = mymodel[i].Material(name='{}'.format(name[i]))
    E = mat_name.Elastic(table=((E[i-1],Nu[i-1]),))
    D = mat_name.Density(table=((Dens[i-1],),))