Skip to content

Commit 1818ce5

Browse files
committed
Update ASDConcrete1D_MakeLaws.py
1 parent 3fbeb25 commit 1818ce5

File tree

1 file changed

+4
-2
lines changed

1 file changed

+4
-2
lines changed

source/user/manual/material/uniaxial/examples/ASDConcrete1D_MakeLaws.py

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -70,6 +70,7 @@ def _make_tension(E, ft, Gt):
7070
Td[i] = 1.0-si/qi # compute damage
7171
return (Te, Ts, Td)
7272

73+
7374
def _make_compression(E, fc, ec, Gc):
7475
'''
7576
a quadratic hardening followed by linear softening for compressive response
@@ -82,8 +83,8 @@ def _make_compression(E, fc, ec, Gc):
8283
Gc1 = fc*(ec-ec_pl)/2.0
8384
Gc2 = max(Gc1*1.0e-2, Gc-Gc1)
8485
ecr = ec + 2.0*Gc2/(fc+fcr)
85-
Ce = [0.0, ec0] # total strain points
86-
Cs = [0.0, fc0] # nominal stress points
86+
Ce = [0.0, ec0] # total strain points
87+
Cs = [0.0, fc0] # nominal stress points
8788
Cpl = [0.0, 0.0] # desired values of equivalent plastic strains
8889
nc = 10
8990
dec = (ec-ec0)/(nc-1)
@@ -92,6 +93,7 @@ def _make_compression(E, fc, ec, Gc):
9293
Ce.append(iec)
9394
Cs.append(_bezier3(iec, ec0, ec1, ec, fc0, fc, fc))
9495
Cpl.append(Cpl[-1]+(iec-Cpl[-1])*0.7)
96+
9597
# end of linear softening - begin residual plateau
9698
Ce.append(ecr)
9799
Cs.append(fcr)

0 commit comments

Comments
 (0)