Dynamic and static cases added with correct mesh ordering (f047426e) · Commits · Asutosh Padhy / ABAQUS SMA

Lagoudas/Own/Old/1_Elasticity.f→Elastic/Python_old/1_Elasticity.f

+0 −0

File moved.

Lagoudas/Own/Old/3D_UM_old.yaml→Elastic/Python_old/3D_UM_old.yaml

+3 −3

Original line number	Diff line number	Diff line
		geometry:
		film_length : 0.5 #x length of film
		film_height : 0.5 # y height of film
		film_width : 1.0 #z width of film
		film_length : 0.75 #x length of film
		film_height : 0.75 # y height of film
		film_width : 0.75 #z width of film
		mesh: 0.25 #Mesh size in mm

		UMAT: "1_Elasticity.f"

Lagoudas/Own/Old/3D_film_old.py→Elastic/Python_old/3D_film_old.py

+16 −15

Original line number	Diff line number	Diff line
		@@ -117,7 +117,7 @@ def seeding_old(part,geometry):
		#################################################################################
		def mesh_old(part,geometry):
		seeding_old(part,geometry)
		elemtype1 = msh.ElemType(elemCode=C3D8, elemLibrary=STANDARD)
		elemtype1 = msh.ElemType(elemCode=C3D20, elemLibrary=STANDARD)
		reg = part.cells
		picked = (reg,)
		part.setElementType(regions=picked, elemTypes=(elemtype1,))
		@@ -167,7 +167,6 @@ def BC(loading, geometry, model):

		#Left nodes fixed in X
		#####################################

		n0= a.instances["film-1"].nodes.getByBoundingBox(-mesh,-mesh,-mesh,mesh,h+mesh,w+mesh)
		region = regionToolset.Region(nodes=n0)
		a.Set(nodes=n0, name="BC-fix-X")
		@@ -178,18 +177,25 @@ def BC(loading, geometry, model):

		####################################################
		# Bottom and top edge nodes fixed in Y and Z
		n0= a.instances["film-1"].nodes.getByBoundingCylinder((0,0,0),(l-mesh,0,0),mesh)
		n1= a.instances["film-1"].nodes.getByBoundingCylinder((0,h,w),(l-mesh,h,w),mesh)
		region = regionToolset.Region(nodes=n1+n0)
		a.Set(nodes=n1+n0, name="BC-fix-YZ")

		# n0= a.instances["film-1"].nodes.getByBoundingCylinder((0,0,0),(l-mesh,0,0),mesh)
		# n1= a.instances["film-1"].nodes.getByBoundingCylinder((0,h,w),(l-mesh,h,w),mesh)
		n0= a.instances["film-1"].nodes.getByBoundingSphere((0,0,0),mesh)
		region = regionToolset.Region(nodes=n0)
		a.Set(nodes=n0, name="BC-fix-YZ")
		model.DisplacementBC(name="BC-YZ",
		region=region, u1=UNSET, u2=0.0,createStepName='Initial',
		u3=0.0, ur1=0.0, ur2=UNSET, ur3=UNSET, amplitude='', fixed=OFF,
		distributionType=UNIFORM, fieldName="", localCsys=None) # Displacements set
		u3=0.0, ur1=UNSET, ur2=UNSET, ur3=UNSET, amplitude='', fixed=OFF,
		distributionType=UNIFORM, fieldName="", localCsys=None)

		# Displacements set
		n1= a.instances["film-1"].nodes.getByBoundingSphere((0,h,0),mesh)
		region = regionToolset.Region(nodes=n1)
		a.Set(nodes=n1, name="BC-fix-Z")
		model.DisplacementBC(name="BC-Z",
		region=region, u1=UNSET, u2=UNSET,createStepName='Initial',
		u3=0.0, ur1=UNSET, ur2=UNSET, ur3=UNSET, amplitude='', fixed=OFF,
		distributionType=UNIFORM, fieldName="", localCsys=None)

		# Displacements set
		###############################################
		# Initial temperature assign to all nodes

		@@ -201,13 +207,10 @@ def BC(loading, geometry, model):
		crossSectionDistribution=CONSTANT_THROUGH_THICKNESS, magnitudes=(325.0,
		))


		# right faces prescribed displacement in X
		#side1Faces = p.faces.getByBoundingBox(-0.01+l,-0.01,-0.01,l+0.01,h+0.01,w+0.01)
		s = a.instances["film-1"].nodes.getByBoundingBox(l-mesh,-mesh,-mesh,l+mesh,h+mesh,w+mesh)



		if loading["type"]==1:
		region = regionToolset.Region(nodes=s)
		a.Set(nodes=s, name="BC-F")
		@@ -218,7 +221,6 @@ def BC(loading, geometry, model):
		# mdb.models['Model-1'].boundaryConditions['Bearing1'].setValuesInStep(
		# stepName='Unload', amplitude='Amp-2')


		#model.boundaryConditions['Disp-1'].setValuesInStep(stepName='Unload', amplitude='')

		elif loading["type"]==2:
		@@ -234,7 +236,6 @@ def BC(loading, geometry, model):
		distributionType=UNIFORM, fieldName="", localCsys=None)

		#model.boundaryConditions['Disp-1'].setValuesInStep(stepName='Unload', amplitude='')

		# ############### set Output requests ###############
		def outputreq(model):
		"""Defines output variables to be shown in Visualization after solving

Lagoudas/Own/Old/JOB.inp

deleted100755 → 0

+0 −200

Original line number	Diff line number	Diff line
		*Heading
		** Job name: JOB Model name: Model-1
		** Generated by: Abaqus/CAE 2018
		*Preprint, echo=NO, model=NO, history=NO, contact=NO
		**
		** PARTS
		**
		*Part, name=film
		*Node
		1, 0.5, 0.5, 1.
		2, 0.5, 0.25, 1.
		3, 0.5, 0., 1.
		4, 0.5, 0.5, 0.75
		5, 0.5, 0.25, 0.75
		6, 0.5, 0., 0.75
		7, 0.5, 0.5, 0.5
		8, 0.5, 0.25, 0.5
		9, 0.5, 0., 0.5
		10, 0.5, 0.5, 0.25
		11, 0.5, 0.25, 0.25
		12, 0.5, 0., 0.25
		13, 0.5, 0.5, 0.
		14, 0.5, 0.25, 0.
		15, 0.5, 0., 0.
		16, 0.25, 0.5, 1.
		17, 0.25, 0.25, 1.
		18, 0.25, 0., 1.
		19, 0.25, 0.5, 0.75
		20, 0.25, 0.25, 0.75
		21, 0.25, 0., 0.75
		22, 0.25, 0.5, 0.5
		23, 0.25, 0.25, 0.5
		24, 0.25, 0., 0.5
		25, 0.25, 0.5, 0.25
		26, 0.25, 0.25, 0.25
		27, 0.25, 0., 0.25
		28, 0.25, 0.5, 0.
		29, 0.25, 0.25, 0.
		30, 0.25, 0., 0.
		31, 0., 0.5, 1.
		32, 0., 0.25, 1.
		33, 0., 0., 1.
		34, 0., 0.5, 0.75
		35, 0., 0.25, 0.75
		36, 0., 0., 0.75
		37, 0., 0.5, 0.5
		38, 0., 0.25, 0.5
		39, 0., 0., 0.5
		40, 0., 0.5, 0.25
		41, 0., 0.25, 0.25
		42, 0., 0., 0.25
		43, 0., 0.5, 0.
		44, 0., 0.25, 0.
		45, 0., 0., 0.
		*Element, type=C3D8
		1, 16, 17, 20, 19, 1, 2, 5, 4
		2, 17, 18, 21, 20, 2, 3, 6, 5
		3, 19, 20, 23, 22, 4, 5, 8, 7
		4, 20, 21, 24, 23, 5, 6, 9, 8
		5, 22, 23, 26, 25, 7, 8, 11, 10
		6, 23, 24, 27, 26, 8, 9, 12, 11
		7, 25, 26, 29, 28, 10, 11, 14, 13
		8, 26, 27, 30, 29, 11, 12, 15, 14
		9, 31, 32, 35, 34, 16, 17, 20, 19
		10, 32, 33, 36, 35, 17, 18, 21, 20
		11, 34, 35, 38, 37, 19, 20, 23, 22
		12, 35, 36, 39, 38, 20, 21, 24, 23
		13, 37, 38, 41, 40, 22, 23, 26, 25
		14, 38, 39, 42, 41, 23, 24, 27, 26
		15, 40, 41, 44, 43, 25, 26, 29, 28
		16, 41, 42, 45, 44, 26, 27, 30, 29
		*Nset, nset=_PickedSet2, internal, generate
		1, 45, 1
		*Elset, elset=_PickedSet2, internal, generate
		1, 16, 1
		** Section: SMA_Section
		*Solid Section, elset=_PickedSet2, material=SMA_UM, CONTROLS=HG-1
		*End Part
		*SECTION CONTROLS, NAME=HG-1, HOURGLASS=ENHANCED
		**
		**
		** ASSEMBLY
		**
		*Assembly, name=Assembly
		**
		*Instance, name=film-1, part=film
		*End Instance
		**
		*Nset, nset=BC-F, instance=film-1, generate
		1, 15, 1
		*Nset, nset=BC-Temp, instance=film-1, generate
		1, 45, 1
		*Nset, nset=BC-fix-X, instance=film-1, generate
		31, 45, 1
		*Nset, nset=BC-fix-YZ, instance=film-1
		16, 30, 31, 45
		*Nset, nset=_PickedSet5, internal, instance=film-1, generate
		31, 45, 1
		*Nset, nset=_PickedSet7, internal, instance=film-1
		16, 30, 31, 45
		*Nset, nset=_PickedSet9, internal, instance=film-1, generate
		1, 45, 1
		*Nset, nset=_PickedSet11, internal, instance=film-1, generate
		1, 15, 1
		*End Assembly
		**
		** MATERIALS
		**
		*Material, name=SMA_UM
		*Conductivity
		22.,
		*Density
		6450.,
		*Depvar
		100,
		*Specific Heat
		329.,
		*User Material, constants=2
		7e+10, 0.33
		**
		** BOUNDARY CONDITIONS
		**
		** Name: BC-YZ Type: Displacement/Rotation
		*Boundary
		_PickedSet7, 2, 2
		_PickedSet7, 3, 3
		_PickedSet7, 4, 4
		** Name: BC-fix-X Type: Displacement/Rotation
		*Boundary
		_PickedSet5, 1, 1
		**
		** PREDEFINED FIELDS
		**
		** Name: Predefined Field-1 Type: Temperature
		*Initial Conditions, type=TEMPERATURE
		_PickedSet9, 325.
		** ----------------------------------------------------------------
		**
		** STEP: Loading
		**
		*Step, name=Loading, nlgeom=NO, inc=1000
		*Static
		1., 100., 0.001, 1.
		**
		** LOADS
		**
		** Name: Load-1 Type: Concentrated force
		*Cload
		_PickedSet11, 1, 150000.
		**
		** OUTPUT REQUESTS
		**
		*Restart, write, frequency=0
		**
		** FIELD OUTPUT: F-Output-1
		**
		*Output, field
		*Node Output
		RF, U
		*Element Output, directions=YES
		S, SDV
		**
		** HISTORY OUTPUT: H-Output-1
		**
		*Output, history
		*Element Output
		IRA1, IRA2, IRA3, IRAR1, IRAR2, IRAR3
		*End Step
		** ----------------------------------------------------------------
		**
		** STEP: Unload
		**
		*Step, name=Unload, nlgeom=NO, inc=1000
		*Static
		1., 100., 0.001, 1.
		**
		** LOADS
		**
		** Name: Load-1 Type: Concentrated force
		*Cload
		_PickedSet11, 1, 0.
		**
		** OUTPUT REQUESTS
		**
		*Restart, write, frequency=0
		**
		** FIELD OUTPUT: F-Output-1, F-Output-2
		**
		*Output, field
		*Node Output
		RF, U
		*Element Output, directions=YES
		S, SDV
		**
		** HISTORY OUTPUT: H-Output-1, H-Output-2
		**
		*Output, history
		*Element Output
		IRA1, IRA2, IRA3, IRAR1, IRAR2, IRAR3
		*End Step

Lagoudas/Own/my_sme/3D_SMA_UM.f→Lagoudas/Own/my_sme_dyn/3D_SMA_UM.f

+0 −0

File moved.

View file