language: AWK (mawk) (mawk-1.3.3)
date: 882 days 16 hours ago
link:
visibility: public 
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
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177

/FILNAM,CARBON
/PREP7
ET,1,42             !ELEMENT TYPE 1 IS PLANE42
ET,2,45             !ELEMENT TYPE 2 IS SOLID45
ET,3,TARGE170        !ELEMENT TYPE 3 IS TARGE170
ET,4,CONTA174         !ELEMENT TYPE 3 IS CONTA174
MP,EX,1,75000       !GLASS SUBSTRATE MAT PROPS
MP,NUXY,1,0.23        
MP,EX,2,4750        !CARBON FILM MAT PROPS
MP,NUXY,2,0.22
!MP,DENS,2,12000
TB,BISO,2,1         !BILINEAR ISOTROPIC HARDENING RULE
TBTEMP,0
TBDATA,1,1.15,0      !YIELD STRENGTH
TBDATA,2,560     !Tangent mod. Ni (Mpa)
MP,EX,3,1141000     !DIAMOND INDENTER MAT PROPS
MP,NUXY,3,0.07
K,1                  !CREATE KEYPOINTS
K,2,4
K,3,10
K,4,,6
K,5,4,6
K,6,10,6
K,7,,7.2
K,8,4,7.2
K,9,10,7.2
K,10,,8
K,11,4,8
K,12,10,8
L,1,2       !CREATE LINES
L,2,3
L,4,5
L,5,6
L,7,8
L,8,9
L,10,11
L,11,12
L,1,4
L,2,5
L,3,6
L,4,7
L,5,8
L,6,9
L,7,10
L,8,11
L,9,12
AL,1,3,9,10         !CREATE AREAS
AL,2,4,10,11
AL,3,5,12,13
AL,4,6,13,14
AL,5,7,15,16
AL,6,8,16,17
LESIZE,7,,,15       ! SPECIFY LINE DIVISIONS
LESIZE,5,,,15
LESIZE,3,,,15
LESIZE,1,,,15
LESIZE,8,,,10,4
LESIZE,6,,,10,4
LESIZE,4,,,10,4
LESIZE,2,,,10,4
LESIZE,9,,,10,0.25
LESIZE,10,,,10,0.25
LESIZE,11,,,10,0.25
LESIZE,12,,,8
LESIZE,13,,,8
LESIZE,14,,,8
LESIZE,15,,,8
LESIZE,16,,,8
LESIZE,17,,,8
MSHKEY,1                         ! ENFORCE MAPPED MESHING
TYPE,1                           ! SWITCH TO ET 1
MAT,1                            ! SWITCH TO MATERIAL 1
AMESH,1,6                        ! CREATE MESH
TYPE,2                           ! SWITCH TO ET 2
ESIZE,,4                         ! SPECIFY # OF ELEMS FOR SWEEP
VROTAT,1,2,3,4,5,6,1,10,45       ! SWEEP AREAS TO CREATE VOLUME
NSEL,S,LOC,Y,6,8              !SELECT NODES AT 6um<=y<= 8um
ESLN,S,1                    !SELECT ELEMENTS WITH SELECTED NODES
EMODIF,ALL,MAT,2            !MODIFY ELEMS TO BE GLASS
TYPE,4                      !SWITCH TO ET 4
AMESH,22                    !CREATE CONTACT AREA MESH
ESEL,S,TYPE,,4              !SELECT ELEMENTS WITH ET 4
ESURF,,REVE                 !ADJUST OUTWARD NORMAL OF ELEMENTS
                            !GENERATE SOLID MODEL FOR INDENTER
K,23,,8
K,26,,10.5
K,24,3.7,9.5
K,25,3.7,10.5
L,23,24
L,24,25
L,25,26
L,26,23
AL,40,41,42,43
LESIZE,41,,,5
LESIZE,43,,,5
LESIZE,40,,,15
LESIZE,42,,,15
TYPE,1
MAT, 3
AMESH,27                        !CREATE MESH FOR INDENTER
TYPE,2                          !SWITCH TO ET 2
ESIZE,,4                        !SPECIFY # OF ELEMS FOR SWEEP
VROTAT,27,,,,,,23,26,45         !SWEEP AREAS TO CREATE VOLUME
TYPE,3                          !SWITCH TO ET 3
AMESH,28                        !CREATE TARGET AREA MESH
ALLSEL,ALL
FINISH                          !EXIT PREPROCESSOR
/SOLU                           !ENTER SOLUTION PROCESSOR
SOLCONTROL,0
ANTYPE,STATIC
NSEL,S,LOC,Z                    !SELECT NODES AT Z = 0
DSYM,SYMM,Z,0                   !APPLY SYMMETRY BG ABOUT X-Y PLANE
CLOCAL,11,,,,,,,45              !DEFINE LOCAL CS
NSEL,S,LOC,Z                   !SELECT NODES AT Z = 0(IN LOCAL CS)
DSYM,SYMM,Z,11
CSYS,0                          !SWITCH TO GLOBAL CARTESIAN CS
NSEL,S,LOC,Y,10.5
D,ALL,UY,0
NSEL,S,LOC,X
NSEL,R,LOC,Z
D,ALL,UX
D,ALL,UZ
*DIM,DIS,ARRAY,57
*DIM,DISR,TABLE,57
*DIM,SUMS,TABLE,57
A=0
*DO,I,1,36
DIS(I)=A
A=A+0.0156
*ENDDO
B=0.015
*DO,I,37,57
DIS(I)=A-B
B=B+0.015
*ENDDO
NLGEOM,ON
AUTOTS,ON
OUTRES,,1
NROPT,FULL,,
NEQIT,100
ALLSEL
*DO,I,1,57
NSEL,S,LOC,Y
D,ALL,UX
D,ALL,UZ
D,ALL,UY,DIS(I)
ALLSEL
LSWRITE,I
*ENDDO
ESEL,U,TYPE,,1
LSSOLVE,1,57
FINISH
/POST1
*DO,J,1,57
SET,J
NSEL,S,LOC,Y
*GET,NUMNOD,NODE,0,COUNT
*GET,CURNOD,NODE,0,NUM,MIN
SUM=0
*DO,I,1,NUMNOD
*GET,RFY,NODE,CURNOD,RF,FY
SUM=SUM+RFY
CURNOD=NDNEXT(CURNOD)
*ENDDO
DISJ=DIS(J)
DISR(J)=DIS(J)
SUM=SUM/100
SUMS(J)=SUM
/OUTPUT,NANO_RF_D,OUT,,APPEND
*VWRITE,DISJ,SUM
(E16.8,5X,E16.8)
/OUTPUT
*ENDDO
*VPLOT,DISR,SUMS
/AXLAB,X,Indentation Depth (micrometer) 
/AXLAB,Y,Forc(mN)