sbosse/simndt2b
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fri 29 Nov 2024 04:30:50 PM CET

Browse Source
This commit is contained in:
sbosse 2024-11-29 16:33:22 +01:00
parent 3748d0b7aa
commit 54752865ff
1 changed files with 60 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

60
src/SimNDT/engine/EFIT2Dnumpy.py Normal file
View File

@ -0,0 +1,60 @@
def sourceVel(Vx, Vy, XL, YL, source, Stype, win, dtdxx):
if Stype == 0:
Vx[XL[:], YL[:]] -= source * dtdxx * win[:]
elif Stype == 1:
Vy[XL[:], YL[:]] -= source * dtdxx * win[:]
elif Stype == 2:
Vx[XL[:], YL[:]] -= source * dtdxx * win[:]
Vy[XL[:], YL[:]] -= source * dtdxx * win[:]
return (Vx, Vy)
def sourceStress(Txx, Tyy, XL, YL, source, Stype, win, dtdxx):
if Stype == 0:
Txx[XL[:], YL[:]] -= source * dtdxx * win[:]
elif Stype == 1:
Tyy[XL[:], YL[:]] -= source * dtdxx * win[:]
elif Stype == 2:
Txx[XL[:], YL[:]] -= source * dtdxx * win[:]
Tyy[XL[:], YL[:]] -= source * dtdxx * win[:]
return (Txx, Tyy)
def velocityVoigt(Txx, Txy, Tyy, vx, vy, d_vx, d_vy, BX, BY, ABS, ddx, dt):
d_vx[:-1, 1:] = (ddx * BX[:-1, 1:]) * (Txx[1:, 1:] - Txx[:-1, 1:] + Txy[:-1, 1:] - Txy[:-1, 0:-1]);
vx[:-1, 1:] += dt * d_vx[:-1, 1:]
d_vy[1:, :-1] = (ddx * BY[1:, :-1]) * (Txy[1:, :-1] - Txy[0:-1, :-1] + Tyy[1:, 1:] - Tyy[1:, :-1]);
vy[1:, :-1] += dt * d_vy[1:, :-1]
vx[:, :] *= ABS[:, :]
vy[:, :] *= ABS[:, :]
return (vx, vy, d_vx, d_vy)
def stressVoigt(Txx, Txy, Tyy, vx, vy, d_vx, d_vy, C11, C12, C44, ETA_VS, ETA_S, ETA_SS, ABS, dtx):
Txx[1:, 1:] += ((dtx * C11[1:, 1:]) * (vx[1:, 1:] - vx[0:-1, 1:]) +
(dtx * C12[1:, 1:]) * (vy[1:, 1:] - vy[1:, 0:-1]) +
(dtx * ETA_VS[1:, 1:]) * (d_vx[1:, 1:] - d_vx[0:-1, 1:]) +
(dtx * ETA_S[1:, 1:]) * (d_vy[1:, 1:] - d_vy[1:, 0:-1]))
Tyy[1:, 1:] += ((dtx * C12[1:, 1:]) * (vx[1:, 1:] - vx[0:-1, 1:]) +
(dtx * C11[1:, 1:]) * (vy[1:, 1:] - vy[1:, 0:-1]) +
(dtx * ETA_S[1:, 1:]) * (d_vx[1:, 1:] - d_vx[0:-1, 1:]) +
(dtx * ETA_VS[1:, 1:]) * (d_vy[1:, 1:] - d_vy[1:, 0:-1]))
Txy[:-1, :-1] += ((dtx * C44[:-1, :-1]) * (vx[:-1, 1:] - vx[:-1, :-1] + vy[1:, :-1] - vy[:-1, :-1]) +
(dtx * ETA_SS[:-1, :-1]) * (d_vx[:-1, 1:] - d_vx[:-1, :-1] + d_vy[1:, :-1] - d_vy[:-1, :-1]))
Txx[:, :] *= ABS[:, :]
Tyy[:, :] *= ABS[:, :]
Txy[:, :] *= ABS[:, :]
return (Txx, Tyy, Txy)