�� sr !com.femlab.server.ModelFileHeader�D���% L tagst Lcom/femlab/util/FlStringList;L typesq ~ L vrsnt Lcom/femlab/util/FlVersion;xpw sr com.femlab.util.FlVersion��%�/B = I buildI majorL datet Ljava/lang/String;L extq ~ L nameq ~ L rcsq ~ L reactionExtq ~ L reactionNameq ~ xpw t q ~ t #COMSOL Reaction Engineering Lab 3.5t at COMSOL 3.5q ~ w [t $Name: $t $Date: 2008/12/03 17:02:19 $xur [Ljava.lang.String;��V��{G xp t modelinfot xfemt guit fem0t g4t g5t fem1t fem1.0t g3q ~ q ~ t fem1416t fem1416.0q ~ q ~ uq ~ q ~ t femstructt guistructq ~ t drawq ~ t geomt meshq ~ t solutiont xmeshq ~ q ~ q ~ q ~ !xsr com.femlab.api.client.ModelInfo�^���% L descrq ~ L docURLq ~ [ imaget [Bxpw pt pxuq ~ t�`
COMSOL 3.5
a
0
603
$Name: $
$Date: 2008/12/03 17:02:19 $
1408
geom
on
on
off
on
on
on
on
0.1
TIME
SmeSolid3
cartesian
{"u","v","w","p"}
{"X","Y","Z"}
smsld
SME
{"shlag(2,'u')","shlag(2,'v')","shlag(2,'w')"}
4
2
2
off
_smsld
EX_smsld
EY_smsld
EZ_smsld
FX_smsld
FXAmp_smsld
FXPh_smsld
FXg_smsld
FY_smsld
FYAmp_smsld
FYPh_smsld
FYg_smsld
FZ_smsld
FZAmp_smsld
FZPh_smsld
FZg_smsld
GXY_smsld
GXZ_smsld
GYZ_smsld
TaX_smsld
TaY_smsld
TaZ_smsld
alphaX_smsld
alphaY_smsld
alphaZ_smsld
eX_smsld
eXi_smsld
eXY_smsld
eXYi_smsld
eXZ_smsld
eXZi_smsld
eY_smsld
eYi_smsld
eYZ_smsld
eYZi_smsld
eZ_smsld
eZi_smsld
nuXY_smsld
nuXZ_smsld
nuYZ_smsld
refpntX_smsld
refpntY_smsld
refpntZ_smsld
sX_smsld
sXi_smsld
sXY_smsld
sXYi_smsld
sXZ_smsld
sXZi_smsld
sY_smsld
sYi_smsld
sYZ_smsld
sYZi_smsld
sZ_smsld
sZi_smsld
Lag2
static
on
freq
off
ref
ref
off
{"lm1","lm2","lm3"}
ideal
0
0
0
0
0
0
global
global
standard
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
{1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}
0
0
0
0
0
0
global
global
standard
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
{1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}
Fluid load
-T_x_ns*dvol_ale/dvol
0
0
0
-T_y_ns*dvol_ale/dvol
0
0
0
-T_z_ns*dvol_ale/dvol
0
0
0
global
0
0
distr_force
global
standard
free
free
fixed
0
0
0
0
0
0
0
0
0
0
0
0
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
0
E_smsld/h*min(1e-3*5^auglagiter,1)
E_smsld/h*min(1e-3*5^auglagiter,1)
nofriction
fast
auto
1e-2
auto
1e-6
0
0
0
0
0
Inf
1e-3
0
1e-3
0
1e-3
0
1e-3
0
{2,2,2,2,2,1,1,3,1,1,1,2}
1 2 3
1 2 3
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1
{"0","0","0","0"}
{"0","0","0","0"}
0
1
Fluid domain
Solid domain
2.0e11
8e6
7850
nodamping
0
0
0
{"0","0","0","0"}
{"0","0","0","0"}
global
iso
0
neo_hookean
8e5
2e5
2e5
1e10
-3e11
-5e11
-6e11
1e11
8e10
{"2e10"}
{"3000"}
4e10
3e10
iso
mises
tangent
2.0e10
2.0e10
2.0e8
mises_smsld
misesKin_smsld
2.0e10[Pa]/(1-2.0e10[Pa]/2.0e11[Pa])*epe_smsld
0
0
0
1.2e-5
0.33
2.0e11
0.33
7.52e10
1.2e-5
0
0
0
0
2.0e11
0.33
7.52e10
1.2e-5
0
0
0
0
2.0e11
0.33
7.52e10
1.2e-5
0
0
0
0
0
{"2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)"}
{"2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)"}
{"1.2e-5","1.2e-5","1.2e-5","0","0","0"}
{"1.2e-5","1.2e-5","1.2e-5","0","0","0"}
0
0
0
0
0
0
global
0
293.15
293.15
0
500
17.44
51.6
global
standard
free
0
0
0
0
0
0
0
0
0
0
0
0
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
dx_guess_smsld
dx_smsld
dy_guess_smsld
dy_smsld
dz_guess_smsld
dz_smsld
dr_guess_smsld
dr_smsld
R0_guess_smsld
X0_guess_smsld
Y0_guess_smsld
Z0_guess_smsld
{"X","Y","Z"}
{"X","Y","Z"}
none
0
0
0
0
0
0
{"0","0","1"}
{"0","0","1"}
global
0
1
{1,2}
-1
0
0
0
MovingMesh
cartesian
{}
{"X","Y","Z"}
ale
{"shlag(2,'x')","shlag(2,'y')","shlag(2,'z')"}
4
2
2
off
_ale
Lag2
winslow
transient
off
ale
ale
ref
off
{"lm4","lm5","lm6"}
non-ideal
Structural displacement
Fixed
global
def
{"u","v","w"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"1","1","1"}
{"1","1","1"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{2,2,2,2,2,1,1,3,1,1,1,2}
1 2 3
1 2 3
1 1 1
1 1 1
1 1 1
1 1 1
{"xinit_ale","yinit_ale","zinit_ale"}
{"xinit_ale","yinit_ale","zinit_ale"}
1
Fluid domain
Solid domain
free
phys
{"0","0","0"}
{"u","v","w"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{1,2}
FlNavierStokes
cartesian
{"u2","v2","w2","p2","nxw","nyw","nzw"}
{"x","y","z"}
ns
{"shlag(2,'u2')","shlag(2,'v2')","shlag(2,'w2')","shlag(1,'p2')"}
4
2
2
1
2
off
_ns
Lagp2p1
static
Off
ale
off
{"lm7","lm8","lm9","lm10"}
non-ideal
0
0
{1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}
Structural velocity
inlet
walltype
int
outlet
walltype
uv
p
noslip
noslip
noslip
noslip
lwall
ntotstress
totstress
cont
u0
U0in
U0in
U0in
U0in
16*umax*y*(width-y)*z*(height-z)/(width^2*height^2)
0
0
0
0
0
0
0
0
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
1
0
0
0
0
0
ut
0
0
0
0
vt
0
0
0
0
wt
{1,2,2,2,2,2,2,3,2,2,2,4}
1 2 3 4
1 2 3 4
1 1 1 2
1 1 1 2
1 1 1 2
1 1 1 2
{"0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0"}
1
0
Fluid domain
Solid domain
1e3
1
1e-3
1
0
0
0
0
0.5
1
gls
0.25
1
sc
0.35
0
1
0.1
0
navierstokes
usesGMG
0
Lagp2p1
0
Lagp2p1
0
0
0
{1,2}
10
{"X","Y","Z"}
{"x","y","z"}
{"ref","ale"}
{"shlag('order',2,'basename','x','frame','ale')","shlag('order',2,'basename','y','frame','ale')","shlag('order',2,'basename','z','frame','ale')","shlag('order',2,'basename','u2','frame','ale')","shlag('order',2,'basename','v2','frame','ale')","shlag('order',2,'basename','w2','frame','ale')","shlag('order',1,'basename','p2','frame','ale')","shlag(2,'u')","shlag(2,'v')","shlag(2,'w')"}
4
2
2
1
dvol
ref
2
fixed
dvol_ale
ale
2
moving_abs
ref
dvol_ale
ale
2
moving_rel
{"u","v","w"}
ref
on
1
SI
{"0","0","0","0","0","0","0","0","0","0","0"}
1 1 1 1 1 1 1 1 1 1 2
{1,1}
1 2 3 4 5 6 7
8 9 10
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
{"0","0","0","0","xinit_ale","yinit_ale","zinit_ale","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-Xx*test(xx)-Xy*test(xy)-Xz*test(xz)","-Yx*test(yx)-Yy*test(yy)-Yz*test(yz)","-Zx*test(zx)-Zy*test(zy)-Zz*test(zz)","(-2*eta_ns*u2x+p2)*test(u2x)-eta_ns*(u2y+v2x)*test(u2y)-eta_ns*(u2z+w2x)*test(u2z)-rho_ns*(u2*u2x+v2*u2y+w2*u2z)*test(u2)+gls_ns+shock_capt_ns","-eta_ns*(v2x+u2y)*test(v2x)+(-2*eta_ns*v2y+p2)*test(v2y)-eta_ns*(v2z+w2y)*test(v2z)-rho_ns*(u2*v2x+v2*v2y+w2*v2z)*test(v2)","-eta_ns*(w2x+u2z)*test(w2x)-eta_ns*(w2y+v2z)*test(w2y)+(-2*eta_ns*w2z+p2)*test(w2z)-rho_ns*(u2*w2x+v2*w2y+w2*w2z)*test(w2)","-divU_ns*test(p2)"}
{"-eX_smsld_test*SX_smsld-eY_smsld_test*SY_smsld-eZ_smsld_test*SZ_smsld-2*eXY_smsld_test*SXY_smsld-2*eYZ_smsld_test*SYZ_smsld-2*eXZ_smsld_test*SXZ_smsld","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
1 2
1 3
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 3 3 3 3 3 3 3
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(v)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
0
0
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
(F11_smsld*(SX_smsld*F11_smsld+SXY_smsld*F12_smsld+SXZ_smsld*F13_smsld)+F12_smsld*(SXY_smsld*F11_smsld+SY_smsld*F12_smsld+SYZ_smsld*F13_smsld)+F13_smsld*(SXZ_smsld*F11_smsld+SYZ_smsld*F12_smsld+SZ_smsld*F13_smsld))/J_smsld
(F21_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F22_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F23_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld
(F31_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F32_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F33_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld
(F11_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F12_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F13_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld
(F21_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F22_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F23_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld
(F11_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F12_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F13_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld
uX+0.5*(uX^2+vX^2+wX^2)
vY+0.5*(uY^2+vY^2+wY^2)
wZ+0.5*(uZ^2+vZ^2+wZ^2)
0.5*(uY+vX+uX*uY+vX*vY+wX*wY)
0.5*(vZ+wY+uY*uZ+vY*vZ+wY*wZ)
0.5*(uZ+wX+uX*uZ+vX*vZ+wX*wZ)
E_smsld*((1-nu_smsld)*eX_smsld+nu_smsld*eY_smsld+nu_smsld*eZ_smsld)/((1+nu_smsld)*(1-2*nu_smsld))
E_smsld*(nu_smsld*eX_smsld+(1-nu_smsld)*eY_smsld+nu_smsld*eZ_smsld)/((1+nu_smsld)*(1-2*nu_smsld))
E_smsld*(nu_smsld*eX_smsld+nu_smsld*eY_smsld+(1-nu_smsld)*eZ_smsld)/((1+nu_smsld)*(1-2*nu_smsld))
E_smsld*eXY_smsld/(1+nu_smsld)
E_smsld*eYZ_smsld/(1+nu_smsld)
E_smsld*eXZ_smsld/(1+nu_smsld)
-K_smsld*evol_smsld
sqrt((K_smsld+4*G_smsld/3)/rho_smsld)
sqrt(G_smsld/rho_smsld)
sqrt(sX_smsld^2+sY_smsld^2+sZ_smsld^2-sX_smsld*sY_smsld-sY_smsld*sZ_smsld-sX_smsld*sZ_smsld+3*sXY_smsld^2+3*sYZ_smsld^2+3*sXZ_smsld^2)
0.5*(sX_smsld*eX_smsld+sY_smsld*eY_smsld+sZ_smsld*eZ_smsld+2*sXY_smsld*eXY_smsld+2*sYZ_smsld*eYZ_smsld+2*sXZ_smsld*eXZ_smsld)
-1+Jel_smsld
1+uX
uY
uZ
vX
1+vY
vZ
wX
wY
1+wZ
F11_smsld*F22_smsld*F33_smsld+F12_smsld*F23_smsld*F31_smsld+F13_smsld*F21_smsld*F32_smsld-F11_smsld*F23_smsld*F32_smsld-F12_smsld*F21_smsld*F33_smsld-F13_smsld*F22_smsld*F31_smsld
detF_smsld
J_smsld
(F22_smsld*F33_smsld-F23_smsld*F32_smsld)/detF_smsld
(F13_smsld*F32_smsld-F12_smsld*F33_smsld)/detF_smsld
(F12_smsld*F23_smsld-F13_smsld*F22_smsld)/detF_smsld
(F31_smsld*F23_smsld-F21_smsld*F33_smsld)/detF_smsld
(F11_smsld*F33_smsld-F31_smsld*F13_smsld)/detF_smsld
(F21_smsld*F13_smsld-F11_smsld*F23_smsld)/detF_smsld
(F21_smsld*F32_smsld-F31_smsld*F22_smsld)/detF_smsld
(F31_smsld*F12_smsld-F11_smsld*F32_smsld)/detF_smsld
(F11_smsld*F22_smsld-F21_smsld*F12_smsld)/detF_smsld
max(max(abs(s1_smsld-s2_smsld),abs(s2_smsld-s3_smsld)),abs(s1_smsld-s3_smsld))
X
Y
Z
x-X
y-Y
z-Z
sqrt(u2^2+v2^2+w2^2)
w2y-v2z
u2z-w2x
v2x-u2y
sqrt(Vx_ns^2+Vy_ns^2+Vz_ns^2)
u2x+v2y+w2z
rho_ns*U_ns*h_ale/eta_ns
rho_ns*(u2*u2x+v2*u2y+w2*u2z)+p2x-F_x_ns+if(gmg_level>0,0,-eta_ns*(2*u2xx+u2yy+v2xy+u2zz+w2xz))
rho_ns*(u2*v2x+v2*v2y+w2*v2z)+p2y-F_y_ns+if(gmg_level>0,0,-eta_ns*(v2xx+u2yx+2*v2yy+v2zz+w2yz))
rho_ns*(u2*w2x+v2*w2y+w2*w2z)+p2z-F_z_ns+if(gmg_level>0,0,-eta_ns*(w2xx+u2zx+w2yy+v2zy+2*w2zz))
rho_ns*u2
rho_ns*v2
rho_ns*w2
eta_ns
rho_ns
nojac(1/max(2*rho_ns*sqrt(emetric_ale(u2-xt,v2-yt,w2-zt)),48*eta_ns/h_ale^2))
0.5*nojac(if((u2-xt)^2+(v2-yt)^2+(w2-zt)^2
rho_ns*divU_ns
8000000
7850
0
0 t�`
0
0
800000
200000
200000
10000000000
-300000000000
-500000000000
-600000000000
E_smsld*nu_smsld/((1+nu_smsld)*(1-2*nu_smsld))
0.5*E_smsld/(1+nu_smsld)
E_smsld/(3*(1-2*nu_smsld))
0.5*E_smsld/(1+nu_smsld)
20000000000
20000000000
200000000
mises_smsld
misesKin_smsld
20000000000[Pa]*epe_smsld/(1-20000000000[Pa]/200000000000[Pa])
1.2e-005
0.33
200000000000
0.33
75200000000
1.2e-005
0
0
0
0
200000000000
0.33
75200000000
1.2e-005
0
0
0
0
200000000000
0.33
75200000000
1.2e-005
0
0
0
0
0
0
0
0
0
0
0
293.15
293.15
500
17.44
51.6
0
0
0
0
0
0
0
0
0
dx_guess_smsld
dx_smsld
dy_guess_smsld
dy_smsld
dz_guess_smsld
dz_smsld
dr_guess_smsld
dr_smsld
R0_guess_smsld
X0_guess_smsld
Y0_guess_smsld
Z0_guess_smsld
0
0
0
1
0
0
0
0
0
0
0
0
0
1000
0.001
0
0
0
0.5
0.25
sc
0.35
1
0.1
0
Lagp2p1
Lagp2p1
-test(nojac(rho_ns)*(nojac(u2-xt)*u2x+nojac(v2-yt)*u2y+nojac(w2-zt)*u2z)+p2x+if(gmg_level>0,0,-nojac(eta_ns)*(2*u2xx+u2yy+v2xy+u2zz+w2xz)))*taum_ns*res_u2_ns-test(nojac(rho_ns)*(nojac(u2-xt)*v2x+nojac(v2-yt)*v2y+nojac(w2-zt)*v2z)+p2y+if(gmg_level>0,0,-nojac(eta_ns)*(v2xx+u2yx+2*v2yy+v2zz+w2yz)))*taum_ns*res_v2_ns-test(nojac(rho_ns)*(nojac(u2-xt)*w2x+nojac(v2-yt)*w2y+nojac(w2-zt)*w2z)+p2z+if(gmg_level>0,0,-nojac(eta_ns)*(w2xx+u2zx+w2yy+v2zy+2*w2zz)))*taum_ns*res_w2_ns-test(divU_ns)*tauc_ns*res_p2_ns
0
(d(x,xi1)^2+d(x,xi2)^2+d(x,xi3)^2)^2+2*(d(y,xi1)*d(x,xi1)+d(y,xi2)*d(x,xi2)+d(y,xi3)*d(x,xi3))^2+(d(y,xi1)^2+d(y,xi2)^2+d(y,xi3)^2)^2+2*(d(z,xi1)*d(x,xi1)+d(z,xi2)*d(x,xi2)+d(z,xi3)*d(x,xi3))^2+2*(d(z,xi1)*d(y,xi1)+d(z,xi2)*d(y,xi2)+d(z,xi3)*d(y,xi3))^2+(d(z,xi1)^2+d(z,xi2)^2+d(z,xi3)^2)^2
sqrt((u2-xt)^2+(v2-yt)^2+(w2-zt)^2)
d(x,xi1)^2+d(x,xi2)^2+d(x,xi3)^2
d(y,xi1)*d(x,xi1)+d(y,xi2)*d(x,xi2)+d(y,xi3)*d(x,xi3)
d(y,xi1)^2+d(y,xi2)^2+d(y,xi3)^2
d(z,xi1)*d(x,xi1)+d(z,xi2)*d(x,xi2)+d(z,xi3)*d(x,xi3)
d(z,xi1)*d(y,xi1)+d(z,xi2)*d(y,xi2)+d(z,xi3)*d(y,xi3)
d(z,xi1)^2+d(z,xi2)^2+d(z,xi3)^2
-nojac(0.5*min(ck_ns*sqrt((tauc_ns*res_p2_ns/(nojac(rhoint_g1_ns/vol_g1_ns)*nojac(sqrt(u2int_g1_ns/vol_g1_ns))))^2+(taum_ns*res_u2_ns)^2+(taum_ns*res_v2_ns)^2+(taum_ns*res_w2_ns)^2),U_ale_ns)*min(gijgij_ns^(-0.25),rho_ns*U_ale_ns/(48*eta_ns)))*rho_ns*(g11_ns*(u2x*test(u2x)+v2x*test(v2x)+w2x*test(w2x))+g21_ns*(u2y*test(u2x)+v2y*test(v2x)+w2y*test(w2x))+g31_ns*(u2z*test(u2x)+v2z*test(v2x)+w2z*test(w2x))+g21_ns*(u2x*test(u2y)+v2x*test(v2y)+w2x*test(w2y))+g22_ns*(u2y*test(u2y)+v2y*test(v2y)+w2y*test(w2y))+g32_ns*(u2z*test(u2y)+v2z*test(v2y)+w2z*test(w2y))+g31_ns*(u2x*test(u2z)+v2x*test(v2z)+w2x*test(w2z))+g32_ns*(u2y*test(u2z)+v2y*test(v2z)+w2y*test(w2z))+g33_ns*(u2z*test(u2z)+v2z*test(v2z)+w2z*test(w2z)))
{0,0}
{0,0}
{1,2}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"FXg_smsld*u_test","FYg_smsld*v_test","FZg_smsld*w_test","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","-nx_ns*p2*test(u2)","-ny_ns*p2*test(v2)","-nz_ns*p2*test(w2)","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-x+X","-y+Y","-z+Z","-u2+u0_ns","-v2","-w2","0"}
{"0","0","0","0","-x+X","-y+Y","-z+Z","-u2","-v2","-w2","0"}
{"0","0","0","0","u-x+X","v-y+Y","w-z+Z","-u2","-v2","-w2","0"}
{"-u","-v","-w","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-x+X","-y+Y","-z+Z","0","0","0","-p2"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","test(-u2)","test(-v2)","test(-w2)","0"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","test(-u2)","test(-v2)","test(-w2)","0"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","test(-u2)","test(-v2)","test(-w2)","0"}
{"test(-u)","test(-v)","test(-w)","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","0","0","0","test(-p2)"}
1 2 3 4 5 6 7
1 2 3 4 5 6 7
8 9 10 1 2 3 4 5 6 7
8 9 10
1 2 3 4 5 6 7
1 2
1 2
1 2
1 3
1 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 3 3 3 3 3 3 3
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(u)
reacf(v)
reacf(v)
reacf(w)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
FX_smsld
0
FY_smsld
0
FZ_smsld
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
sqrt(real(u)^2+real(v)^2+real(w)^2)
(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld
(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld
(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld
(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld
(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld
(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld
X
Y
Z
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-T_x_ns*dvol_ale/dvol
0
0
0
-T_y_ns*dvol_ale/dvol
0
0
0
-T_z_ns*dvol_ale/dvol
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
dnX
nX
dnY
nY
dnZ
nZ
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
nX
nX
unX
nX
nY
nY
unY
nY
nZ
nZ
unZ
nZ
16*umax*y*(width-y)*z*(height-z)/(width^2*height^2)
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
1
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
nx
nx
unx
nx
ny
ny
uny
ny
nz
nz
unz
nz
0
0
0
0
0
0
0
0
0
0
0
0
{0,0,0,0,0,0,0,0,0,0,0,0}
{0,0,0,0,0,0,0,0,0,0,0,0}
{1,2,2,2,2,3,3,4,3,3,3,5}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
1 2 3 4 5 6 7
8 9 10 1 2 3 4 5 6 7
1 2
1 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(v)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
0
0
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
X
Y
Z
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
t�`
0
0
0
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
1 2 3 4 5 6 7
8 9 10 1 2 3 4 5 6 7
1 2
1 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(v)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
0
0
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
X
Y
Z
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{1,1,1,1,2,2,2,2,2,2,2,2,1,1,1,1}
-1
0
0
0
uTIME
uTIMETIME
uXTIME
uYTIME
uZTIME
vTIME
vTIMETIME
vXTIME
vYTIME
vZTIME
wTIME
wTIMETIME
wXTIME
wYTIME
wZTIME
xTIME
xTIMETIME
xxTIME
xyTIME
xzTIME
yTIME
yTIMETIME
yxTIME
yyTIME
yzTIME
zTIME
zTIMETIME
zxTIME
zyTIME
zzTIME
10
u2TIME-xTIME*u2x-yTIME*u2y-zTIME*u2z
u2TIMETIME
u2xTIME-xTIME*u2xx-yTIME*u2xy-zTIME*u2xz
u2yTIME-xTIME*u2yx-yTIME*u2yy-zTIME*u2yz
u2zTIME-xTIME*u2zx-yTIME*u2zy-zTIME*u2zz
v2TIME-xTIME*v2x-yTIME*v2y-zTIME*v2z
v2TIMETIME
v2xTIME-xTIME*v2xx-yTIME*v2xy-zTIME*v2xz
v2yTIME-xTIME*v2yx-yTIME*v2yy-zTIME*v2yz
v2zTIME-xTIME*v2zx-yTIME*v2zy-zTIME*v2zz
w2TIME-xTIME*w2x-yTIME*w2y-zTIME*w2z
w2TIMETIME
w2xTIME-xTIME*w2xx-yTIME*w2xy-zTIME*w2xz
w2yTIME-xTIME*w2yx-yTIME*w2yy-zTIME*w2yz
w2zTIME-xTIME*w2zx-yTIME*w2zy-zTIME*w2zz
p2TIME-xTIME*p2x-yTIME*p2y-zTIME*p2z
p2TIMETIME
p2xTIME-xTIME*p2xx-yTIME*p2xy-zTIME*p2xz
p2yTIME-xTIME*p2yx-yTIME*p2yy-zTIME*p2yz
p2zTIME-xTIME*p2zx-yTIME*p2zy-zTIME*p2zz
elpric
{"1"}
{"real(sX_smsld)","real(sXY_smsld)","real(sY_smsld)","real(sXZ_smsld)","real(sYZ_smsld)","real(sZ_smsld)"}
{"2"}
s
smsld
{"X","Y","Z"}
elpric
{"1"}
{"real(eX_smsld)","real(eXY_smsld)","real(eY_smsld)","real(eXZ_smsld)","real(eYZ_smsld)","real(eZ_smsld)"}
{"2"}
e
smsld
{"X","Y","Z"}
elcplscalar
{"1"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1","2","3","4","13","14","15","16"}
{"5","6","7","8","9","10","11","12"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1","2","3","4","5","6","7","8","21","22","23","24"}
{"9","10","11","12","13","14","15","16","17","18","19","20"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1","2","3","4","5","12"}
{"6","7","8","9","10","11"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1"}
{"2"}
{"FXtotpnt_g1_smsld","FYtotpnt_g1_smsld","FZtotpnt_g1_smsld","MXtotpnt_g1_smsld","MYtotpnt_g1_smsld","MZtotpnt_g1_smsld","FXtotedg_g1_smsld","FYtotedg_g1_smsld","FZtotedg_g1_smsld","MXtotedg_g1_smsld","MYtotedg_g1_smsld","MZtotedg_g1_smsld","FXtotbnd_g1_smsld","FYtotbnd_g1_smsld","FZtotbnd_g1_smsld","MXtotbnd_g1_smsld","MYtotbnd_g1_smsld","MZtotbnd_g1_smsld","FXtotsub_g1_smsld","FYtotsub_g1_smsld","FZtotsub_g1_smsld","MXtotsub_g1_smsld","MYtotsub_g1_smsld","MZtotsub_g1_smsld"}
{"1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24"}
elcplsum
{"1"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMXpnt_smsld
RMYpnt_smsld
RMZpnt_smsld
{"1","2","3","4","13","14","15","16"}
{"5","6","7","8","9","10","11","12"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMX_smsld
RMY_smsld
RMZ_smsld
{"1","2","3","4","5","6","7","8","21","22","23","24"}
{"9","10","11","12","13","14","15","16","17","18","19","20"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMX_smsld
RMY_smsld
RMZ_smsld
{"1","2","3","4","5","12"}
{"6","7","8","9","10","11"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMX_smsld
RMY_smsld
RMZ_smsld
{"1"}
{"2"}
{"RFXtotpnt_g1_smsld","RFYtotpnt_g1_smsld","RFZtotpnt_g1_smsld","RMXtotpnt_g1_smsld","RMYtotpnt_g1_smsld","RMZtotpnt_g1_smsld","RFXtotedg_g1_smsld","RFYtotedg_g1_smsld","RFZtotedg_g1_smsld","RMXtotedg_g1_smsld","RMYtotedg_g1_smsld","RMZtotedg_g1_smsld","RFXtotbnd_g1_smsld","RFYtotbnd_g1_smsld","RFZtotbnd_g1_smsld","RMXtotbnd_g1_smsld","RMYtotbnd_g1_smsld","RMZtotbnd_g1_smsld","RFXtotsub_g1_smsld","RFYtotsub_g1_smsld","RFZtotsub_g1_smsld","RMXtotsub_g1_smsld","RMYtotsub_g1_smsld","RMZtotsub_g1_smsld"}
{"1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24"}
elconst
+FXtotpnt_g1_smsld+FXtotedg_g1_smsld+FXtotbnd_g1_smsld+FXtotsub_g1_smsld
RFXtotsub_g1_smsld
+FYtotpnt_g1_smsld+FYtotedg_g1_smsld+FYtotbnd_g1_smsld+FYtotsub_g1_smsld
RFYtotsub_g1_smsld
+FZtotpnt_g1_smsld+FZtotedg_g1_smsld+FZtotbnd_g1_smsld+FZtotsub_g1_smsld
RFZtotsub_g1_smsld
+MXtotpnt_g1_smsld+MXtotedg_g1_smsld+MXtotbnd_g1_smsld+MXtotsub_g1_smsld
RMXtotsub_g1_smsld
+MYtotpnt_g1_smsld+MYtotedg_g1_smsld+MYtotbnd_g1_smsld+MYtotsub_g1_smsld
RMYtotsub_g1_smsld
+MZtotpnt_g1_smsld+MZtotedg_g1_smsld+MZtotbnd_g1_smsld+MZtotsub_g1_smsld
RMZtotsub_g1_smsld
elcplscalar
{"1"}
-qual_ale
4
{"1","2"}
{"minqual_tmp1_ale"}
{"1"}
{"1"}
elconst
-minqual_tmp1_ale
elcplscalar
{"1"}
nojac(rho_ns)
nojac(1)
nojac(max((+(u2-xt)^2+(v2-yt)^2+(w2-zt)^2),(eta_ns/rho_ns/h_ale)^2))
nojac(-1/max(emetric_ale(u2-xt,v2-yt,w2-zt),(visc_vel_fact_ns*eta_ns/rho_ns)^2*emetric2_ale))
ref ref ref ref ref ref ref ref
0 0
0 0
0 0
0 0
{"1"}
{"2"}
{"rhoint_g1_ns","vol_g1_ns","u2int_g1_ns","mdt2_cfl_g1_ns"}
{"1","2","3","4"}
{"4"}
elconst
nojac(sqrt(-mdt2_cfl_g1_ns))
{"BLK1","EXT1"}
1409
geom
on
on
off
on
off
off
on
0.1
t
0
{"x","y"}
{"xy"}
1
on
1
SI
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
{"E1"}
1.5 [m/s]
1 [m]
1 [m]
{0}
{0}
smsld
weak
SI
COMSOL 3.5
a
0
603
$Name: $
$Date: 2008/12/03 17:02:19 $
1408
geom
on
on
off
on
on
on
on
0.1
TIME
SmeSolid3
cartesian
{"u","v","w","p"}
{"X","Y","Z"}
smsld
SME
{"shlag(2,'u')","shlag(2,'v')","shlag(2,'w')"}
4
2
2
off
_smsld
EX_smsld
EY_smsld
EZ_smsld
FX_smsld
FXAmp_smsld
FXPh_smsld
FXg_smsld
FY_smsld
FYAmp_smsld
FYPh_smsld
FYg_smsld
FZ_smsld
FZAmp_smsld
FZPh_smsld
FZg_smsld
GXY_smsld
GXZ_smsld
GYZ_smsld
TaX_smsld
TaY_smsld
TaZ_smsld
alphaX_smsld
alphaY_smsld
alphaZ_smsld
eX_smsld
eXi_smsld
eXY_smsld
eXYi_smsld
eXZ_smsld
eXZi_smsld
eY_smsld
eYi_smsld
eYZ_smsld
eYZi_smsld
eZ_smsld
eZi_smsld
nuXY_smsld
nuXZ_smsld
nuYZ_smsld
refpntX_smsld
refpntY_smsld
refpntZ_smsld
sX_smsld
sXi_smsld
sXY_smsld
sXYi_smsld
sXZ_smsld
sXZi_smsld
sY_smsld
sYi_smsld
sYZ_smsld
sYZi_smsld
sZ_smsld
sZi_smsld
Lag2
static
on
freq
off
ref
ref
off
{"lm1","lm2","lm3"}
ideal
0
0
0
0
0
0
global
global
standard
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
{1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}
0
0
0
0
0
0
global
global
standard
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
{1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}
Fluid load
-T_x_ns*dvol_ale/dvol
0
0
0
-T_y_ns*dvol_ale/dvol
0
0
0
-T_z_ns*dvol_ale/dvol
0
0
0
global
0
0
distr_force
global
standard
free
free
fixed
0
0
0
0
0
0
0
0
0
0
0
0
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
0
E_smsld/h*min(1e-3*5^auglagiter,1)
E_smsld/h*min(1e-3*5^auglagiter,1)
nofriction
fast
auto
1e-2
auto
1e-6
0
0
0
0
0
Inf
1e-3
0
1e-3
0
1e-3
0
1e-3
0
{2,2,2,2,2,1,1,3,1,1,1,2}
1 2 3
1 2 3
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1
{"0","0","0","0"}
{"0","0","0","0"}
0
1
Fluid domain
Solid domain
2.0e11
8e6
7850
nodamping
0
0
0
{"0","0","0","0"}
{"0","0","0","0"}
global
iso
0
neo_hookean
8e5
2e5
2e5
1e10
-3e11
-5e11
-6e11
1e11
8e10
{"2e10"}
{"3000"}
4e10
3e10
iso
mises
tangent
2.0e10
2.0e10
2.0e8
mises_smsld
misesKin_smsld
2.0e10[Pa]/(1-2.0e10[Pa]/2.0e11[Pa])*epe_smsld
0
0
0
1.2e-5
0.33
2.0e11
0.33
7.52e10
1.2e-5
0
0
0
0
2.0e11
0.33
7.52e10
1.2e-5
0
0
0
0
2.0e11
0.33
7.52e10
1.2e-5
0
0
0
0
0
{"2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)"}
{"2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","2.0e11/((1+0.33)*(1-2*0.33))*0.33","0","0","0","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*0.33","2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)","0","0","0","0","0","0","2.0e11/((1+0.33)*2)"}
{"1.2e-5","1.2e-5","1.2e-5","0","0","0"}
{"1.2e-5","1.2e-5","1.2e-5","0","0","0"}
0
0
0
0
0
0
global
0
293.15
293.15
0
500
17.44
51.6
global
standard
free
0
0
0
0
0
0
0
0
0
0
0
0
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0"}
{"0","0","0"}
{"0","0","0"}
0
0
0
0
0
0
dx_guess_smsld
dx_smsld
dy_guess_smsld
dy_smsld
dz_guess_smsld
dz_smsld
dr_guess_smsld
dr_smsld
R0_guess_smsld
X0_guess_smsld
Y0_guess_smsld
Z0_guess_smsld
{"X","Y","Z"}
{"X","Y","Z"}
none
0
0
0
0
0
0
{"0","0","1"}
{"0","0","1"}
global
0
1
{1,2}
-1
0
0
0
MovingMesh
cartesian
{}
{"X","Y","Z"}
ale
{"shlag(2,'x')","shlag(2,'y')","shlag(2,'z')"}
4
2
2
off
_ale
Lag2
winslow
transient
off
ale
ale
ref
off
{"lm4","lm5","lm6"}
non-ideal
Structural displacement
Fixed
global
def
{"u","v","w"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"1","1","1"}
{"1","1","1"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{2,2,2,2,2,1,1,3,1,1,1,2}
1 2 3
1 2 3
1 1 1
1 1 1
1 1 1
1 1 1
{"xinit_ale","yinit_ale","zinit_ale"}
{"xinit_ale","yinit_ale","zinit_ale"}
1
Fluid domain
Solid domain
free
phys
{"0","0","0"}
{"u","v","w"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{1,2}
FlNavierStokes
cartesian
{"u2","v2","w2","p2","nxw","nyw","nzw"}
{"x","y","z"}
ns
{"shlag(2,'u2')","shlag(2,'v2')","shlag(2,'w2')","shlag(1,'p2')"}
4
2
2
1
2
off
_ns
Lagp2p1
static
Off
ale
off
{"lm7","lm8","lm9","lm10"}
non-ideal
0
0
{1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}
Structural velocity
inlet
walltype
int
outlet
walltype
uv
p
noslip
noslip
noslip
noslip
lwall
ntotstress
totstress
cont
u0
U0in
U0in
U0in
U0in
16*umax*y*(width-y)*z*(height-z)/(width^2*height^2)
0
0
0
0
0
0
0
0
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
{"0","0","0"}
1
0
0
0
0
0
ut
0
0
0
0
vt
0
0
0
0
wt
{1,2,2,2,2,2,2,3,2,2,2,4}
1 2 3 4
1 2 3 4
1 1 1 2
1 1 1 2
1 1 1 2
1 1 1 2
{"0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0"}
1
0
Fluid domain
Solid domain
1e3
1
1e-3
1
0
0
0
0
0.5
1
gls
0.25
1
sc
0.35
0
1
0.1
0
navierstokes
usesGMG
0
Lagp2p1
0
Lagp2p1
0
0
0
{1,2}
10
{"X","Y","Z"}
{"x","y","z"}
{"ref","ale"}
{"shlag('order',2,'basename','x','frame','ale')","shlag('order',2,'basename','y','frame','ale')","shlag('order',2,'basename','z','frame','ale')","shlag('order',2,'basename','u2','frame','ale')","shlag('order',2,'basename','v2','frame','ale')","shlag('order',2,'basename','w2','frame','ale')","shlag('order',1,'basename','p2','frame','ale')","shlag(2,'u')","shlag(2,'v')","shlag(2,'w')"}
4
2
2
1
dvol
ref
2
fixed
dvol_ale
ale
2
moving_abs
ref
dvol_ale
ale
2
moving_rel
{"u","v","w"}
ref
on
1
SI
{"0","0","0","0","0","0","0","0","0","0","0"}
1 1 1 1 1 1 1 1 1 1 2
{1,1}
1 2 3 4 5 6 7
8 9 10
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
{"0","0","0","0","xinit_ale","yinit_ale","zinit_ale","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-Xx*test(xx)-Xy*test(xy)-Xz*test(xz)","-Yx*test(yx)-Yy*test(yy)-Yz*test(yz)","-Zx*test(zx)-Zy*test(zy)-Zz*test(zz)","(-2*eta_ns*u2x+p2)*test(u2x)-eta_ns*(u2y+v2x)*test(u2y)-eta_ns*(u2z+w2x)*test(u2z)-rho_ns*(u2*u2x+v2*u2y+w2*u2z)*test(u2)+gls_ns+shock_capt_ns","-eta_ns*(v2x+u2y)*test(v2x)+(-2*eta_ns*v2y+p2)*test(v2y)-eta_ns*(v2z+w2y)*test(v2z)-rho_ns*(u2*v2x+v2*v2y+w2*v2z)*test(v2)","-eta_ns*(w2x+u2z)*test(w2x)-eta_ns*(w2y+v2z)*test(w2y)+(-2*eta_ns*w2z+p2)*test(w2z)-rho_ns*(u2*w2x+v2*w2y+w2*w2z)*test(w2)","-divU_ns*test(p2)"}
{"-eX_smsld_test*SX_smsld-eY_smsld_test*SY_smsld-eZ_smsld_test*SZ_smsld-2*eXY_smsld_test*SXY_smsld-2*eYZ_smsld_test*SYZ_smsld-2*eXZ_smsld_test*SXZ_smsld","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
1 2
1 3
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 3 3 3 3 3 3 3
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(v)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
0
0
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
(F11_smsld*(SX_smsld*F11_smsld+SXY_smsld*F12_smsld+SXZ_smsld*F13_smsld)+F12_smsld*(SXY_smsld*F11_smsld+SY_smsld*F12_smsld+SYZ_smsld*F13_smsld)+F13_smsld*(SXZ_smsld*F11_smsld+SYZ_smsld*F12_smsld+SZ_smsld*F13_smsld))/J_smsld
(F21_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F22_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F23_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld
(F31_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F32_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F33_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld
(F11_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F12_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F13_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld
(F21_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F22_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F23_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld
(F11_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F12_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F13_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld
uX+0.5*(uX^2+vX^2+wX^2)
vY+0.5*(uY^2+vY^2+wY^2)
wZ+0.5*(uZ^2+vZ^2+wZ^2)
0.5*(uY+vX+uX*uY+vX*vY+wX*wY)
0.5*(vZ+wY+uY*uZ+vY*vZ+wY*wZ)
0.5*(uZ+wX+uX*uZ+vX*vZ+wX*wZ)
E_smsld*((1-nu_smsld)*eX_smsld+nu_smsld*eY_smsld+nu_smsld*eZ_smsld)/((1+nu_smsld)*(1-2*nu_smsld))
E_smsld*(nu_smsld*eX_smsld+(1-nu_smsld)*eY_smsld+nu_smsld*eZ_smsld)/((1+nu_smsld)*(1-2*nu_smsld))
E_smsld*(nu_smsld*eX_smsld+nu_smsld*eY_smsld+(1-nu_smsld)*eZ_smsld)/((1+nu_smsld)*(1-2*nu_smsld))
E_smsld*eXY_smsld/(1+nu_smsld)
E_smsld*eYZ_smsld/(1+nu_smsld)
E_smsld*eXZ_smsld/(1+nu_smsld)
-K_smsld*evol_smsld
sqrt((K_smsld+4*G_smsld/3)/rho_smsld)
sqrt(G_smsld/rho_smsld)
sqrt(sX_smsld^2+sY_smsld^2+sZ_smsld^2-sX_smsld*sY_smsld-sY_smsld*sZ_smsld-sX_smsld*sZ_smsld+3*sXY_smsld^2+3*sYZ_smsld^2+3*sXZ_smsld^2)
0.5*(sX_smsld*eX_smsld+sY_smsld*eY_smsld+sZ_smsld*eZ_smsld+2*sXY_smsld*eXY_smsld+2*sYZ_smsld*eYZ_smsld+2*sXZ_smsld*eXZ_smsld)
-1+Jel_smsld
1+uX
uY
uZ
vX
1+vY
vZ
wX
wY
1+wZ
F11_smsld*F22_smsld*F33_smsld+F12_smsld*F23_smsld*F31_smsld+F13_smsld*F21_smsld*F32_smsld-F11_smsld*F23_smsld*F32_smsld-F12_smsld*F21_smsld*F33_smsld-F13_smsld*F22_smsld*F31_smsld
detF_smsld
J_smsld
(F22_smsld*F33_smsld-F23_smsld*F32_smsld)/detF_smsld
(F13_smsld*F32_smsld-F12_smsld*F33_smsld)/detF_smsld
(F12_smsld*F23_smsld-F13_smsld*F22_smsld)/detF_smsld
(F31_smsld*F23_smsld-F21_smsld*F33_smsld)/detF_smsld
(F11_smsld*F33_smsld-F31_smsld*F13_smsld)/detF_smsld
(F21_smsld*F13_smsld-F11_smsld*F23_smsld)/detF_smsld
(F21_smsld*F32_smsld-F31_smsld*F22_smsld)/detF_smsld
(F31_smsld*F12_smsld-F11_smsld*F32_smsld)/detF_smsld
(F11_smsld*F22_smsld-F21_smsld*F12_smsld)/detF_smsld
max(max(abs(s1_smsld-s2_smsld),abs(s2_smsld-s3_smsld)),abs(s1_smsld-s3_smsld))
X
Y
Z
x-X
y-Y
z-Z
sqrt(u2^2+v2^2+w2^2)
w2y-v2z
u2z-w2x
v2x-u2y
sqrt(Vx_ns^2+Vy_ns^2+Vz_ns^2)
u2x+v2y+w2z
rho_ns*U_ns*h_ale/eta_ns
rho_ns*(u2*u2x+v2*u2y+w2*u2z)+p2x-F_x_ns+if(gmg_level>0,0,-eta_ns*(2*u2xx+u2yy+v2xy+u2zz+w2xz))
rho_ns*(u2*v2x+v2*v2y+w2*v2z)+p2y-F_y_ns+if(gmg_level>0,0,-eta_ns*(v2xx+u2yx+2*v2yy+v2zz+w2yz))
rho_ns*(u2*w2x+v2*w2y+w2*w2z)+p2z-F_z_ns+if(gmg_level>0,0,-eta_ns*(w2xx+u2zx+w2yy+v2zy+2*w2zz))
rho_ns*u2
rho_ns*v2
rho_ns*w2
eta_ns
rho_ns
nojac(1/max(2*rho_ns*sqrt(emetric_ale(u2-xt,v2-yt,w2-zt)),48*eta_ns/h_ale^2))
0.5*nojac(if((u2-xt)^2+(v2-yt)^2+(w2-zt)^2
rho_ns*divU_ns
8000000
7850
0
0
0
0
800000
200000
200000
10000000000
-300000000000
-500000000000
-600000000000
E_smsld*nu_smsld/((1+nu_smsld)*(1-2*nu_smsld))
0.5*E_smsld/(1+nu_smsld)
E_smsld/(3*(1-2*nu_smsld))
0.5*E_smsld/(1+nu_smsld)
20000000000
20000000000
200000000
mises_smsld
misesKin_smsld
20000000000[Pa]*epe_smsld/(1-20000000000[Pa]/200000000000[Pa])
1.2e-005
0.33
200000000000
0.33
75200000000
1.2e-005
0
0
0
0
200000000000
0.33
75200000000
1.2e-005
0
0
0
0
200000000000
0.33
75200000000
1.2e-005
0
0
0
0
0
0
0
0
0
0
0
293.15
293.15
500
17.44
51.6
0
0
0
0
0
0
0
0
0
dx_guess_smsld
dx_smsld
dy_guess_smsld
dy_smsld
dz_guess_smsld
dz_smsld
dr_guess_smsld
dr_smsld
R0_guess_smsld
X0_guess_smsld
Y0_guess_smsld
Z0_guess_smsld
0
0
0
1
0
0
0
0
0
0
0
0
0
1000
0.001
0
0
0
0.5
0.25
sc
0.35
1
0.1
0
Lagp2p1
Lagp2p1
-test(nojac(rho_ns)*(nojac(u2-xt)*u2x+nojac(v2-yt)*u2y+nojac(w2-zt)*u2z)+p2x+if(gmg_level>0,0,-nojac(eta_ns)*(2*u2xx+u2yy+v2xy+u2zz+w2xz)))*taum_ns*res_u2_ns-test(nojac(rho_ns)*(nojac(u2-xt)*v2x+nojac(v2-yt)*v2y+nojac(w2-zt)*v2z)+p2y+if(gmg_level>0,0,-nojac(eta_ns)*(v2xx+u2yx+2*v2yy+v2zz+w2yz)))*taum_ns*res_v2_ns-test(nojac(rho_ns)*(nojac(u2-xt)*w2x+nojac(v2-yt)*w2y+nojac(w2-zt)*w2z)+p2z+if(gmg_level>0,0,-nojac(eta_ns)*(w2xx+u2zx+w2yy+v2zy+2*w2zz)))*taum_ns*res_w2_ns-test(divU_ns)*tauc_ns*res_p2_ns
0
(d(x,xi1)^2+d(x,xi2)^2+d(x,xi3)^2)^2+2*(d(y,xi1)*d(x,xi1)+d(y,xi2)*d(x,xi2)+d(y,xi3)*d(x,xi3))^2+(d(y,xi1)^2+d(y,xi2)^2+d(y,xi3)^2)^2+2*(d(z,xi1)*d(x,xi1)+d(z,xi2)*d(x,xi2)+d(z,xi3)*d(x,xi3))^2+2*(d(z,xi1)*d(y,xi1)+d(z,xi2)*d(y,xi2)+d(z,xi3)*d(y,xi3))^2+(d(z,xi1)^2+d(z,xi2)^2+d(z,xi3)^2)^2
sqrt((u2-xt)^2+(v2-yt)^2+(w2-zt)^2)
d(x,xi1)^2+d(x,xi2)^2+d(x,xi3)^2
d(y,xi1)*d(x,xi1)+d(y,xi2)*d(x,xi2)+d(y,xi3)*d(x,xi3)
d(y,xi1)^2+d(y,xi2)^2+d(y,xi3)^2
d(z,xi1)*d(x,xi1)+d(z,xi2)*d(x,xi2)+d(z,xi3)*d(x,xi3)
d(z,xi1)*d(y,xi1)+d(z,xi2)*d(y,xi2)+d(z,xi3)*d(y,xi3)
d(z,xi1)^2+d(z,xi2)^2+d(z,xi3)^2
-nojac(0.5*min(ck_ns*sqrt((tauc_ns*res_p2_ns/(nojac(rhoint_g1_ns/vol_g1_ns)*nojac(sqrt(u2int_g1_ns/vol_g1_ns))))^2+(taum_ns*res_u2_ns)^2+(taum_ns*res_v2_ns)^2+(taum_ns*res_w2_ns)^2),U_ale_ns)*min(gijgij_ns^(-0.25),rho_ns*U_ale_ns/(48*eta_ns)))*rho_ns*(g11_ns*(u2x*test(u2x)+v2x*test(v2x)+w2x*test(w2x))+g21_ns*(u2y*test(u2x)+v2y*test(v2x)+w2y*test(w2x))+g31_ns*(u2z*test(u2x)+v2z*test(v2x)+w2z*test(w2x))+g21_ns*(u2x*test(u2y)+v2x*test(v2y)+w2x*test(w2y))+g22_ns*(u2y*test(u2y)+v2y*test(v2y)+w2y*test(w2y))+g32_ns*(u2z*test(u2y)+v2z*test(v2y)+w2z*test(w2y))+g31_ns*(u2x*test(u2z)+v2x*test(v2z)+w2x*test(w2z))+g32_ns*(u2y*test(u2z)+v2y*test(v2z)+w2y*test(w2z))+g33_ns*(u2z*test(u2z)+v2z*test(v2z)+w2z*test(w2z)))
{0,0}
{0,0}
{1,2}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"FXg_smsld*u_test","FYg_smsld*v_test","FZg_smsld*w_test","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","-nx_ns*p2*test(u2)","-ny_ns*p2*test(v2)","-nz_ns*p2*test(w2)","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-x+X","-y+Y","-z+Z","-u2+u0_ns","-v2","-w2","0"}
{"0","0","0","0","-x+X","-y+Y","-z+Z","-u2","-v2","-w2","0"}
{"0","0","0","0","u-x+X","v-y+Y","w-z+Z","-u2","-v2","-w2","0"}
{"-u","-v","-w","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-x+X","-y+Y","-z+Z","0","0","0","-p2"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","test(-u2)","test(-v2)","test(-w2)","0"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","test(-u2)","test(-v2)","test(-w2)","0"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","test(-u2)","test(-v2)","test(-w2)","0"}
{"test(-u)","test(-v)","test(-w)","0","0","0","0","0","0","0","0"}
{"0","0","0","0","-test(x)","-test(y)","-test(z)","0","0","0","test(-p2)"}
1 2 3 4 5 6 7
1 2 3 4 5 6 7
8 9 10 1 2 3 4 5 6 7
8 9 10
1 2 3 4 5 6 7
1 2
1 2
1 2
1 3
1 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 3 3 3 3 3 3 3
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(u)
reacf(v)
reacf(v)
reacf(w)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
FX_smsld
0
FY_smsld
0
FZ_smsld
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
sqrt(real(u)^2+real(v)^2+real(w)^2)
(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld
(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld
(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld
(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld
(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld
(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld
X
Y
Z
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
eta_ns*(2*nx_ns*u2x+ny_ns*(u2y+v2x)+nz_ns*(u2z+w2x))
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
-nx_ns*p2+2*nx_ns*eta_ns*u2x+ny_ns*eta_ns*(u2y+v2x)+nz_ns*eta_ns*(u2z+w2x)
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
eta_ns*(nx_ns*(v2x+u2y)+2*ny_ns*v2y+nz_ns*(v2z+w2y))
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
-ny_ns*p2+nx_ns*eta_ns*(v2x+u2y)+2*ny_ns*eta_ns*v2y+nz_ns*eta_ns*(v2z+w2y)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
eta_ns*(nx_ns*(w2x+u2z)+ny_ns*(w2y+v2z)+2*nz_ns*w2z)
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-nz_ns*p2+nx_ns*eta_ns*(w2x+u2z)+ny_ns*eta_ns*(w2y+v2z)+2*nz_ns*eta_ns*w2z
-T_x_ns*dvol_ale/dvol
0
0
0
-T_y_ns*dvol_ale/dvol
0
0
0
-T_z_ns*dvol_ale/dvol
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
dnX
nX
dnY
nY
dnZ
nZ
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
nX
nX
unX
nX
nY
nY
unY
nY
nZ
nZ
unZ
nZ
16*umax*y*(width-y)*z*(height-z)/(width^2*height^2)
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
1
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
nx
nx
unx
nx
ny
ny
uny
ny
nz
nz
unz
nz
0
0
0
0
0
0
0
0
0
0
0
0
{0,0,0,0,0,0,0,0,0,0,0,0}
{0,0,0,0,0,0,0,0,0,0,0,0}
{1,2,2,2,2,3,3,4,3,3,3,5}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
1 2 3 4 5 6 7
8 9 10 1 2 3 4 5 6 7
1 2
1 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 2
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(v)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
0
0
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
X
Y
Z
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,1}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
{"0","0","0","0","0","0","0","0","0","0","0"}
1 2 3 4 5 6 7
8 9 10 1 2 3 4 5 6 7
1 2
1 2
1 1 1 1 2 2 2 2 2 2 2
1 1 1 1 2 2 2 2 2 2 2
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"","","","","","","","","","",""}
{"",t�:"","","","","","","","","",""}
{"u","v","w","p","x","y","z","u2","v2","w2","p2"}
reacf(u)
reacf(v)
reacf(w)
(Y+v-refpntY_smsld)*RFZ_smsld-(Z+w-refpntZ_smsld)*RFY_smsld
(Z+w-refpntZ_smsld)*RFX_smsld-(X+u-refpntX_smsld)*RFZ_smsld
(X+u-refpntX_smsld)*RFY_smsld-(Y+v-refpntY_smsld)*RFX_smsld
0
0
0
sqrt(real(u)^2+real(v)^2+real(w)^2)
X
Y
Z
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}
{1,1,1,1,2,2,2,2,2,2,2,2,1,1,1,1}
-1
0
0
0
uTIME
uTIMETIME
uXTIME
uYTIME
uZTIME
vTIME
vTIMETIME
vXTIME
vYTIME
vZTIME
wTIME
wTIMETIME
wXTIME
wYTIME
wZTIME
xTIME
xTIMETIME
xxTIME
xyTIME
xzTIME
yTIME
yTIMETIME
yxTIME
yyTIME
yzTIME
zTIME
zTIMETIME
zxTIME
zyTIME
zzTIME
10
u2TIME-xTIME*u2x-yTIME*u2y-zTIME*u2z
u2TIMETIME
u2xTIME-xTIME*u2xx-yTIME*u2xy-zTIME*u2xz
u2yTIME-xTIME*u2yx-yTIME*u2yy-zTIME*u2yz
u2zTIME-xTIME*u2zx-yTIME*u2zy-zTIME*u2zz
v2TIME-xTIME*v2x-yTIME*v2y-zTIME*v2z
v2TIMETIME
v2xTIME-xTIME*v2xx-yTIME*v2xy-zTIME*v2xz
v2yTIME-xTIME*v2yx-yTIME*v2yy-zTIME*v2yz
v2zTIME-xTIME*v2zx-yTIME*v2zy-zTIME*v2zz
w2TIME-xTIME*w2x-yTIME*w2y-zTIME*w2z
w2TIMETIME
w2xTIME-xTIME*w2xx-yTIME*w2xy-zTIME*w2xz
w2yTIME-xTIME*w2yx-yTIME*w2yy-zTIME*w2yz
w2zTIME-xTIME*w2zx-yTIME*w2zy-zTIME*w2zz
p2TIME-xTIME*p2x-yTIME*p2y-zTIME*p2z
p2TIMETIME
p2xTIME-xTIME*p2xx-yTIME*p2xy-zTIME*p2xz
p2yTIME-xTIME*p2yx-yTIME*p2yy-zTIME*p2yz
p2zTIME-xTIME*p2zx-yTIME*p2zy-zTIME*p2zz
elpric
{"1"}
{"real(sX_smsld)","real(sXY_smsld)","real(sY_smsld)","real(sXZ_smsld)","real(sYZ_smsld)","real(sZ_smsld)"}
{"2"}
s
smsld
{"X","Y","Z"}
elpric
{"1"}
{"real(eX_smsld)","real(eXY_smsld)","real(eY_smsld)","real(eXZ_smsld)","real(eYZ_smsld)","real(eZ_smsld)"}
{"2"}
e
smsld
{"X","Y","Z"}
elcplscalar
{"1"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1","2","3","4","13","14","15","16"}
{"5","6","7","8","9","10","11","12"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1","2","3","4","5","6","7","8","21","22","23","24"}
{"9","10","11","12","13","14","15","16","17","18","19","20"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1","2","3","4","5","12"}
{"6","7","8","9","10","11"}
FXg_smsld
FYg_smsld
FZg_smsld
(((Y+v-refpntY_smsld))*(FZg_smsld)-((Z+w-refpntZ_smsld))*(FYg_smsld))
(((Z+w-refpntZ_smsld))*(FXg_smsld)-((X+u-refpntX_smsld))*(FZg_smsld))
(((X+u-refpntX_smsld))*(FYg_smsld)-((Y+v-refpntY_smsld))*(FXg_smsld))
ref ref ref ref ref ref ref ref ref ref ref ref
4
4
4
4
4
4
{"1"}
{"2"}
{"FXtotpnt_g1_smsld","FYtotpnt_g1_smsld","FZtotpnt_g1_smsld","MXtotpnt_g1_smsld","MYtotpnt_g1_smsld","MZtotpnt_g1_smsld","FXtotedg_g1_smsld","FYtotedg_g1_smsld","FZtotedg_g1_smsld","MXtotedg_g1_smsld","MYtotedg_g1_smsld","MZtotedg_g1_smsld","FXtotbnd_g1_smsld","FYtotbnd_g1_smsld","FZtotbnd_g1_smsld","MXtotbnd_g1_smsld","MYtotbnd_g1_smsld","MZtotbnd_g1_smsld","FXtotsub_g1_smsld","FYtotsub_g1_smsld","FZtotsub_g1_smsld","MXtotsub_g1_smsld","MYtotsub_g1_smsld","MZtotsub_g1_smsld"}
{"1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24"}
elcplsum
{"1"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMXpnt_smsld
RMYpnt_smsld
RMZpnt_smsld
{"1","2","3","4","13","14","15","16"}
{"5","6","7","8","9","10","11","12"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMX_smsld
RMY_smsld
RMZ_smsld
{"1","2","3","4","5","6","7","8","21","22","23","24"}
{"9","10","11","12","13","14","15","16","17","18","19","20"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMX_smsld
RMY_smsld
RMZ_smsld
{"1","2","3","4","5","12"}
{"6","7","8","9","10","11"}
RFX_smsld
RFY_smsld
RFZ_smsld
RMX_smsld
RMY_smsld
RMZ_smsld
{"1"}
{"2"}
{"RFXtotpnt_g1_smsld","RFYtotpnt_g1_smsld","RFZtotpnt_g1_smsld","RMXtotpnt_g1_smsld","RMYtotpnt_g1_smsld","RMZtotpnt_g1_smsld","RFXtotedg_g1_smsld","RFYtotedg_g1_smsld","RFZtotedg_g1_smsld","RMXtotedg_g1_smsld","RMYtotedg_g1_smsld","RMZtotedg_g1_smsld","RFXtotbnd_g1_smsld","RFYtotbnd_g1_smsld","RFZtotbnd_g1_smsld","RMXtotbnd_g1_smsld","RMYtotbnd_g1_smsld","RMZtotbnd_g1_smsld","RFXtotsub_g1_smsld","RFYtotsub_g1_smsld","RFZtotsub_g1_smsld","RMXtotsub_g1_smsld","RMYtotsub_g1_smsld","RMZtotsub_g1_smsld"}
{"1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24"}
elconst
+FXtotpnt_g1_smsld+FXtotedg_g1_smsld+FXtotbnd_g1_smsld+FXtotsub_g1_smsld
RFXtotsub_g1_smsld
+FYtotpnt_g1_smsld+FYtotedg_g1_smsld+FYtotbnd_g1_smsld+FYtotsub_g1_smsld
RFYtotsub_g1_smsld
+FZtotpnt_g1_smsld+FZtotedg_g1_smsld+FZtotbnd_g1_smsld+FZtotsub_g1_smsld
RFZtotsub_g1_smsld
+MXtotpnt_g1_smsld+MXtotedg_g1_smsld+MXtotbnd_g1_smsld+MXtotsub_g1_smsld
RMXtotsub_g1_smsld
+MYtotpnt_g1_smsld+MYtotedg_g1_smsld+MYtotbnd_g1_smsld+MYtotsub_g1_smsld
RMYtotsub_g1_smsld
+MZtotpnt_g1_smsld+MZtotedg_g1_smsld+MZtotbnd_g1_smsld+MZtotsub_g1_smsld
RMZtotsub_g1_smsld
elcplscalar
{"1"}
-qual_ale
4
{"1","2"}
{"minqual_tmp1_ale"}
{"1"}
{"1"}
elconst
-minqual_tmp1_ale
elcplscalar
{"1"}
nojac(rho_ns)
nojac(1)
nojac(max((+(u2-xt)^2+(v2-yt)^2+(w2-zt)^2),(eta_ns/rho_ns/h_ale)^2))
nojac(-1/max(emetric_ale(u2-xt,v2-yt,w2-zt),(visc_vel_fact_ns*eta_ns/rho_ns)^2*emetric2_ale))
ref ref ref ref ref ref ref ref
0 0
0 0
0 0
0 0
{"1"}
{"2"}
{"rhoint_g1_ns","vol_g1_ns","u2int_g1_ns","mdt2_cfl_g1_ns"}
{"1","2","3","4"}
{"4"}
elconst
nojac(sqrt(-mdt2_cfl_g1_ns))
1409
geom
on
on
off
on
off
off
on
0.1
t
0
{"x","y"}
{"xy"}
1
on
1
SI
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
{}
1.5 [m/s]
1 [m]
1 [m]
{0}
{0}
smsld
weak
SI
uq ~ t�`
off
timeseg
time
auto
1.0E-6
25
tol
100
1
off
on
off
1.0
1.0E-4
10.0
off
on
0.0010
25
lumped
automatic
auto
on
1
4
tol
10
1
off
const
1.0
minimal
1.0
1.0E-2
10.0
0.0010
0.01
range(0,0.1,2)
0
1
0.1
11
step
none
tlist
free
bdf_ida
5*dt_cfl_g1_ns
off
15
off
on
on
on
on
on
on
off
on
on
0.0010
5
1
0.1
0.75
linear
off
off
maybe
bweuler
0
off
6
0
300
0.0
0
lambda
0
step
none
off
1
off
0.0
0.0
0.0
on
auto
auto
auto
eliminate
auto
1000
on
none
on
off
off
on
off
currgeom
l2
10000000
longest
weak
on
0
1
2
1
fltpft
1
1.1
0.8
0.2
on
off
off
sensitivity
adjoint
on
cholesky
analytic
500
10
1.0e-6
1.0e-6
1.0e-6
1.0e-6
0
off
w,v,u,u2,p2,v2,z,y,w2,x
w,v,u,u2,p2,v2,z,y,x,w2
on
init_expr_currsol_radio
Automatic
0
u_init_radio
Automatic
0
off
off
off
off
spooles
linsolver
0.0
0.1
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
2
1.0E-6
400.0
10000
left
right
1
50
1
1.0
3
6
v
5000
512.0
off
2.0
on
droptol
gmres
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
u2 v2 w2 p2
1e-3
1
gmres
linsolver
0.0
0.1
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
2
1.0e-3
20
10000
right
right
1
100
1
1.0
3
6
v
5000
512.0
off
2.0
on
droptol
direct
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
gmg
precond
0.01
1.0
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
2
1.0e-3
20
10000
right
right
1
100
1
1.0
3
6
v
5000
512.0
off
2.0
on
droptol
direct
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
vanka
presmoother
0.01
0.1
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
2
1.0e-3
20
10000
right
right
1
100
1
1.0
3
6
v
5000
512.0
off
2.0
on
droptol
direct
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
vanka
postsmoother
0.01
0.1
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
2
1.0e-3
20
10000
right
right
1
100
1
1.0
3
6
v
5000
512.0
off
2.0
on
droptol
direct
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
pardiso
coarsesolver
0.0
0.1
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
10
0.1
1
500
right
right
1
100
1
1.0
3
6
v
5000
512.0
off
2.0
on
droptol
direct
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
u v w p x y z
1e-3
1
gmres
linsolver
0.0
0.1
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
5
1.0e-3
20
10000
right
right
1
100
1
0.8
3
6
v
5000
512.0
off
2.0
on
droptol
gmres
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
ssor
precond
0.01
1.0
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
5
1.0e-3
20
10000
right
right
1
100
1
0.8
3
6
v
5000
512.0
off
2.0
on
droptol
gmres
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1 Geom2
on
on
1
iter
0.5
1
1
20
10
1.0E-2
1
const
off
1.0
1.0E-4
1.0E-2
10.0
minimal
onevery
2
iter
0.5
1
1
20
10
1.0E-2
1
const
off
1.0
1.0E-4
1.0E-2
10.0
minimal
onevery
off
off
off
none
gmres
1.0e-3
20
100
right
gmg
shape
vanka
vanka
direct
pardiso
ssor
5
0.8
5*dt_cfl_g1_ns
none
gmres
1.0e-3
20
100
right
gmg
shape
vanka
vanka
direct
pardiso
ssor
5
0.8
5*dt_cfl_g1_ns
none
gmres
1.0e-3
20
100
right
gmg
shape
vanka
vanka
direct
pardiso
ssor
5
0.8
5*dt_cfl_g1_ns
none
gmres
1.0e-3
20
100
right
gmg
shape
vanka
vanka
direct
pardiso
ssor
5
0.8
5*dt_cfl_g1_ns
Geom1
0
off
off
off
1027175
COMSOL 3.5.0.603
1280511348908
1277744001404
off
html
off
a4
off
full
off
Geom1
post
1
-6.1257422745431E-17
2.5
0.0
1.100000023841858
0.0
1.0
0.5
0.5
0.5
-1.5508119515295125
1.5508119515295125
-1.1606183662945708
1.1606183662945708
-6.417132113755014
-9.442001095650228
7.771519808666241
1.25
0.550000011920929
0.5
0.0
0.0
1.0
point
headlight
on
255
255
255
-6.417132
-9.442001
7.7715197
directional
light 1
on
255
255
255
1.0
-1.0
1.0
directional
light 2
on
255
255
255
-1.0
1.0
-1.0
directional
light 3
on
255
255
255
-1.0
-1.0
1.0
directional
light 4
on
255
255
255
1.0
-1.0
-1.0
on
on
on
on
off
on
on
on
on
on
on
on
off
off
off
off
1.0
1.0
1.0
size
cycle
off
off
0.5
orbit
none
none
off
off
off
off
on
off
off
on
on
on
on
on
on
1.0E-6
intersection
1.0
normal
off
2
1.0
off
on
off
on
off
off
on
on
6
on
1.5
0.7
1.6
0.04
0.4
20
1.0
1.0
1.0
on
0
longest
{"","","","","","","","","","","","","","","",""}
{"","","","","","","","","","","","","","","",""}
{"","","","","","","","","","","","","","","","","","","","","","","",""}
{"","","","","","","","","","","","","","","","","","","","","","","",""}
{"","","","","","","","","","","","","","","","","","","","","","","",""}
{"","","","","","","","","","","","","","","","","","","","","","","",""}
{"","","","","","","","","","","",""}
{"","","","","","","","","","","",""}
{"","","","","","","","","","","",""}
{"","","","","","","","","","","",""}
{"tri","tri","tri","tri","tri","tri","tri","tri","tri","tri","tri","tri"}
{"",""}
{"",""}
on
{"U_ns"}
on
0.0
0.0
on
off
m/s
0
off
0
0.55
on
0
.1
on
step
none
step
none
step
none
Rainbow
off
on
on
255,0,0
interp
fill
off
{"mises_smsld"}
on
off
Pa
5
off
step
none
{"mises_smsld"}
on
off
Pa
Rainbow
off
on
on
255,0,0
interp
fill
off
1
random
{"mises_smsld"}
on
155.92214158073534
268786.2380466773
on
off
Pa
Rainbow
off
on
on
255,0,0
interp
fill
off
{"disp_smsld"}
on
0.0
0.08263585652656665
on
off
m
Rainbow
off
on
on
255,0,0
interp
fill
off
{"disp_smsld"}
on
0.0
0.08263585652656665
on
off
m
Rainbow
off
on
on
255,0,0
interp
off
sub
{"u","v","w"}
off
m
{"u","v","w"}
off
m
{"u","v","w"}
off
m
7
off
7
off
7
off
step
none
step
none
step
none
cone
proportional
255,0,0
on
1
off
sub
{"s1_smsld","s1X_smsld","s1Y_smsld","s1Z_smsld","s2_smsld","s2X_smsld","s2Y_smsld","s2Z_smsld","s3_smsld","s3X_smsld","s3Y_smsld","s3Z_smsld"}
off
Pa
{"s1_smsld","s1X_smsld","s1Y_smsld","s1Z_smsld","s2_smsld","s2X_smsld","s2Y_smsld","s2Z_smsld","s3_smsld","s3X_smsld","s3Y_smsld","s3Z_smsld"}
off
Pa
5
off
5
off
5
off
step
none
step
none
step
none
cone
proportional
0,153,0
on
1
off
{"u","v","w"}
m
off
255,0,0
{"mises_smsld"}
Pa
Rainbow
off
on
line
{"mises_smsld"}
off
Pa
on
1
medium
20
sub
off
step
none
step
none
step
none
0.01
0.01
0.01
Inf
5000
on
off
0.15
0.05
0.15
off
1
0.5
1.3
0.5
none
off
mass
lines
khanforce_ns
radius,partr,1e-4
4*pi/3*1e-9
{"sign(u2-partu)*pi*partr^2*rho_ns*(1.84*(abs(u2-partu))^(2/3.45-0.31)*(2*partr*rho_ns/eta_ns)^(-0.31)+0.293*(abs(u2-partu))^(2/3.45+0.06)*(2*partr*rho_ns/eta_ns)^0.06)^3.45","sign(v2-partv)*pi*partr^2*rho_ns*(1.84*(abs(v2-partv))^(2/3.45-0.31)*(2*partr*rho_ns/eta_ns)^(-0.31)+0.293*(abs(v2-partv))^(2/3.45+0.06)*(2*partr*rho_ns/eta_ns)^0.06)^3.45","sign(w2-partw)*pi*partr^2*rho_ns*(1.84*(abs(w2-partw))^(2/3.45-0.31)*(2*partr*rho_ns/eta_ns)^(-0.31)+0.293*(abs(w2-partw))^(2/3.45+0.06)*(2*partr*rho_ns/eta_ns)^0.06)^3.45"}
{"u","v","w"}
m
on
{"0","0","0"}
{"0","0","0"}
step
none
step
none
step
none
off
255,0,0
{"mises_smsld"}
Pa
Rainbow
off
on
255,0,0
on
1
once
step
none
stick
line
{"mises_smsld"}
off
Pa
on
1
medium
0.001
off
{""}
{""}
off
on
on
0.001
partt
off
5
0.001
off
{"",""}
{"",""}
off
on
on
0.001
{"partu","partv","partw"}
partt
off
5
off
{"mises_smsld"}
off
on
Pa
{"disp_smsld"}
off
off
m
{"disp_smsld"}
off
off
m
on
off
on
5
{"Geom1"}
0
ale
0
stored
step
none
on
on
off
all
on
off
off
on
on
on
{"u","v","w"}
m
{"u","v","w"}
m
{"u","v","w"}
m
on
0.293668
0
stored
step
none
moviefile
AVI
640
480
10
full
11
off
off
off
{"disp_smsld"}
off
m
0
0
stored
step
none
off
auto
4
on
{"disp_smsld"}
off
m^2
0
ale
0
stored
step
none
off
auto
4
on
{"disp_smsld"}
off
m^3
0
ale
0
stored
step
none
off
auto
4
on
{"mises_smsld"}
off
J
0
ale
0
stored
step
none
{"disp_smsld"}
on
on
off
m
Rainbow
off
on
on
255,0,0
interp
fill
arc
off
{"disp_smsld"}
m
{"disp_smsld"}
m
off
{"disp_smsld"}
on
off
m
arc
off
{"disp_smsld"}
m
solid
cyclecolor
none
off
off
255,0,0
Rainbow
off
on
on
255,0,0
interp
fill
{"disp_smsld"}
off
m
off
{"disp_smsld"}
m
solid
cyclecolor
none
off
off
255,0,0
255,0,0
0
ale
0
stored
step
none
on
on
on
{"lin","lin"}
on
on
1
0
1
on
{"mises_smsld"}
on
on
off
Pa
{"0","1","0","0","0","1","0","0","0"}
on
global
5
off
off
step
none
Rainbow
off
on
on
255,0,0
interp
fill
off
{"mises_smsld"}
off
Pa
arc
off
{"mises_smsld"}
Pa
solid
cyclecolor
none
off
off
255,0,0
Rainbow
off
on
on
255,0,0
interp
fill
200
{"0","1","0","0","0","0"}
{"mises_smsld"}
off
Pa
off
{"mises_smsld"}
J
solid
cyclecolor
none
off
off
255,0,0
{"0","0","0"}
step
none
step
none
step
none
255,0,0
0
ale
0
stored
step
none
on
on
on
{"lin","lin"}
on
1
{"disp_smsld"}
off
m
element
2
{"disp_smsld"}
off
off
m
element
2
{"disp_smsld"}
off
off
m
element
2
{"mises_smsld"}
off
off
Pa
10
off
10
off
10
off
step
none
step
none
step
none
element
2
0
ale
0
stored
step
none
on
ptd
off
3
coorddata
off
off
on
0
ale
0
stored
step
none
step
none
step
none
{"mises_smsld"}
off
Pa
{"0","0","0"}
{}
{}
off
solid
cyclecolor
none
off
off
255,0,0
on
on
on
{"lin","lin"}
0
0
stored
step
none
{}
{}
{}
{}
off
off
192,192,192
on
on
0,0,0
on
off
128,128,128
on
on
0,0,0
on
0,0,0
off
all
1
random
Rainbow
off
on
Geom2
draw
-1
-0.125
2.625
-0.46621621621621623
1.4662162162162162
-1.0
1.0
0.025
0.025
0.2
-1.0
1.0
-1.0
1.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
1.0
0.0
on
on
on
on
off
on
on
on
off
on
on
on
off
off
off
off
1.0
1.0
1.0
size
cycle
off
off
0.5
orbit
none
none
off
off
off
off
on
off
off
on
on
on
on
on
on
intersection
1.0
normal
off
2
10.0
off
on
off
on
off
off
on
on
5
on
1
0.3
1.3
0.001
1
10
1.0
1.0
on
0
regular
{}
{}
{}
{}
{}
{}
{}
{}
{}
on
on
on
off
off
Rainbow
off
on
on
255,0,0
interp
fill
off
on
off
off
on
off
Rainbow
off
on
on
255,0,0
20
off
off
step
none
off
off
on
on
off
off
Rainbow
off
on
on
255,0,0
interp
off
sub
off
off
off
15
off
15
off
step
none
step
none
arrow
proportional
255,0,0
on
1
off
off
off
8
off
8
off
step
none
step
none
arrow
proportional
0,153,0
on
1
off
off
255,0,0
Rainbow
off
on
off
20
sub
off
step
none
step
none
0.001
0.01
0.01
Inf
5000
on
off
0.05
20
off
1
0.5
1.3
0.5
none
off
mass
lines
1
{"","",""}
on
{"0","0","0"}
partstart_explicit
10
off
step
none
{"0","0"}
step
none
step
none
off
255,0,0
Rainbow
off
on
255,0,0
on
1
once
step
none
stick
0.001
off
{""}
{""}
off
on
on
0.001
partt
off
5
0.001
off
{"",""}
{"",""}
off
on
on
0.001
{"partu","partv","partw"}
partt
off
5
off
off
on
off
off
on
off
on
1
{"Geom2"}
0
0
stored
step
none
on
on
off
all
on
off
off
on
on
on
1
0
stored
step
none
moviefile
AVI
640
480
10
full
11
off
off
off
off
0
0
stored
step
none
off
auto
4
on
off
0
0
stored
step
none
off
auto
4
on
off
0
0
stored
step
none
on
on
off
Rainbow
off
on
on
255,0,0
interp
fill
off
on
off
arc
off
solid
cyclecolor
none
off
off
255,0,0
Rainbow
off
on
on
255,0,0
interp
fill
off
off
solid
cyclecolor
none
off
off
255,0,0
255,0,0
0
0
stored
step
none
on
on
on
{"lin","lin"}
on
on
1
0
off
off
arc
off
solid
cyclecolor
none
off
off
255,0,0
Rainbow
off
on
on
255,0,0
interp
fill
200
{"0","1","0","0"}
5
off
off
step
none
off
off
solid
cyclecolor
none
off
off
255,0,0
{"0","0"}
step
none
step
none
255,0,0
0
0
stored
step
none
on
on
on
{"lin","lin"}
on
0
off
element
2
off
off
element
2
off
off
10
off
10
off
step
none
step
none
element
2
0
0
stored
step
none
on
ptd
off
2
coorddata
off
off
on
0
0
stored
step
none
step
none
off
{"0","0"}
{}
{}
off
solid
cyclecolor
none
off
off
255,0,0
on
on
on
{"lin","lin"}
0
0
stored
step
none
0
0,1,0,0,0,1,0,0,0
0,1,0,0,0,0,1,0,0,0,0,1
Geom2
Geom1
10
0
up
on
off
128,128,128
off
on
192,192,192
on
0,0,0
off
all
1
random
Rainbow
off
on
uq ~ q ~ sr com.femlab.geom.Blockx��a��E� D lxD lyD lzL baseq ~ L lxExprt !Lcom/femlab/geom/Geom$Expression;L lyExprq ~ 2L lzExprq ~ 2xr com.femlab.geom.Prim3�n.D��� D rot[ axist [D[ axisExprt "[Lcom/femlab/geom/Geom$Expression;[ posq ~ 4[ posExprq ~ 5L rotExprq ~ 2xr com.femlab.geom.Geom3]�(�.� xr com.femlab.geom.Geom�6��{ L geomAssoct Lcom/femlab/geom/GeomAssoc;L versionq ~ xpw q ~ w !�ur [B���T� xp !� G e o m 3 �����|�=-C��6? ���������|�= @ �<���������|�= �������? ���������|�= @�������? �<���������|�= ��� �?���������|�= @ �?���������|�= ����������? �?���������|�= @�������? �?���������|�= ���� � ���� � ���� � �? ���� � �? �?���� � �?���� � �?���� � �? ���� � �? �?���� � �? �?���� � �? ���� � �?���� � �? ���� � �? �?���� � �?���� � �?���� � �? ���� � �? �?���� � �? �?���� � �? ���� � �?���� � ���� � ���� � ���� � �?���� � �?���� � �?���� � �?���� � �?���� � �?���� � �?���� � �?���� � �?���� � �?���� � �?���� � �?���� � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � � � � � � � B e z i e r T r i �? @ �<