-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathCoF_modification_script.py
More file actions
59 lines (40 loc) · 2.52 KB
/
Copy pathCoF_modification_script.py
File metadata and controls
59 lines (40 loc) · 2.52 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
# This is just a demo to show how the code is structured. For the full, working version, please visit our website.
# ######################################################################
# ################# CAE Assistant Company ################
# ############## www CAEassistant com #############
# ########### Copy right by CAE Assistant Company ###############
## ######################################################################
# ONLY the BUYER of this package has permission to use its codes.
# Any distribution of this subroutine is illegal and will be prosecuted
# ######################################################################
# ######################################################################
# CAE Assisitant Services:
# Toturial Packages,Consultancy,Articles,Q&A,Video Gallery,Online Course
# ######################################################################
# Need help with your project?
# You can get initial free consultation from (Support CAEassistant com)
# ######################################################################
from abaqusConstants import SCALAR, CENTROID
from numpy import sin, cos, sqrt
COF = [0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65]
Sf = 835 # MPa
damage = []
angle = []
for cof in COF:
#Hidden Content Please purchase the full product to access the complete code.
mdb.jobs['Fretting'].submit()
mdb.jobs['Fretting'].waitForCompletion()
odb = session.openOdb("Fretting.odb")
results = odb.rootAssembly.instances['PLANE-1'].elementSets['RESULTS']
#Hidden Content Please purchase the full product to access the complete code.
for i, j, k in zip(Emax_sub.values, Emin_sub.values, Smax_sub.values):
Ex_max = i.data[0]
Ey_max = i.data[1]
#Hidden Content Please purchase the full product to access the complete code.
for o in range(-90, 91):
#Hidden Content Please purchase the full product to access the complete code.
Ymin = (-(Ex_min - Ey_min)*sin(2*radians(o))/2 + Exy_min*cos(2*radians(o))/2)*2
#Hidden Content Please purchase the full product to access the complete code.
if FS > FS_local:
FS_local = FS
#Hidden Content Please purchase the full product to access the complete code.