49 lines
1.6 KiB
Python
49 lines
1.6 KiB
Python
|
from abaqusConstants import *
|
||
|
|
from abaqus import *
|
||
|
|
import random
|
||
|
|
incFile=open('Ori29.dat','w')
|
||
|
|
# seeds 1-10
|
||
|
|
#random.seed(79789)
|
||
|
|
#random.seed(673673)
|
||
|
|
#random.seed(868256)
|
||
|
|
#random.seed(15675676)
|
||
|
|
#random.seed(8836572)
|
||
|
|
# seeds 20-29
|
||
|
|
random.seed(6753679)
|
||
|
|
random.seed(88822657)
|
||
|
|
random.seed(267567)
|
||
|
|
random.seed(833683655)
|
||
|
|
random.seed(566445)
|
||
|
|
for i in range(1,1000):
|
||
|
|
# Generate 'c' direction
|
||
|
|
angle1=random.random()*2.*pi
|
||
|
|
sigma=0.02 # extruded rod
|
||
|
|
a2=random.normalvariate(1., sigma)
|
||
|
|
if a2>1.:
|
||
|
|
a2=1.-(a2-1.)
|
||
|
|
angle2=acos(a2)
|
||
|
|
# angle2=acos(random.random()) # homogeneous
|
||
|
|
cx=cos(angle1)*cos(angle2)
|
||
|
|
cy=sin(angle1)*cos(angle2)
|
||
|
|
cz=sin(angle2)
|
||
|
|
# Get arbitary normal to C
|
||
|
|
nx=cos(angle1)*cos(angle2+pi/2.)
|
||
|
|
ny=sin(angle1)*cos(angle2+pi/2.)
|
||
|
|
nz=sin(angle2+pi/2.)
|
||
|
|
# Rotate about random angle
|
||
|
|
a3=random.random()*2.*pi
|
||
|
|
xx=nx*(cos(a3)+cx*cx*(1.-cos(a3)))+ny*(cx*cy*(1.-cos(a3))-cz*sin(a3))+nz*(cx*cz*(1.-cos(a3))+cy*sin(a3))
|
||
|
|
xy=nx*(cx*cy*(1.-cos(a3))+cz*sin(a3))+ny*(cos(a3)+cy*cy*(1.-cos(a3)))+nz*(cy*cz*(1.-cos(a3))-cx*sin(a3))
|
||
|
|
xz=nx*(cz*cx*(1.-cos(a3))-cy*sin(a3))+ny*(cz*cy*(1.-cos(a3))+cx*sin(a3))+nz*(cos(a3)+cz*cz*(1.-cos(a3)))
|
||
|
|
# Y-Axis is normal to C and X
|
||
|
|
yx=xy*cz-xz*cy
|
||
|
|
yy=xz*cx-xx*cz
|
||
|
|
yz=xx*cy-xy*cx
|
||
|
|
incFile.write('*MATERIAL, NAME=Mat '+str(i)+ '\n')
|
||
|
|
incFile.write('*USER MATERIAL, CONSTANTS=22, UNSYMM \n')
|
||
|
|
incFile.write('45000.,0.3,'+str(xx)+','+str(xy)+','+str(xz)+','+str(yx)+','+str(yy)+','+str(yz)+' \n')
|
||
|
|
incFile.write('4.,1.,0.,10.,20.,150.,7500.,40., \n')
|
||
|
|
incFile.write('260.,7500.,5.,200.,0.11,10. \n')
|
||
|
|
incFile.write('*Depvar \n')
|
||
|
|
incFile.write(str(163)+' \n')
|
||
|
|
incFile.close()
|