Hajakala
Mechanical
- Mar 23, 2017
- 16
I am trying to initially constrain my model using MPC and then release node by node when the field variable switches from 0 to 1.
In the odb file I can see the FV changing from 0 to 1 but I am not sure how to use the field in the mpc subroutine.
In the subroutine mpc then field looks different to that in usdfld and ufield:
USDFLD Field(1)
UFIELD Field(1,1)
MPC Field(1,1,1)
like above the field has different dimensions and it doesn't seem to match in each subroutine.
In my MPC subroutine when I try to do as below, it gives me an error saying that "The number of subscripts is incorrect.":
IF (Field(1).LT.1) THEN
UE(1) = U(2,1)
A(1,1,1) = 1.
A(1,1,2) = 0.0
ELSE
LMPC=0.0
UE(1) = 0.0
A(1,1,1) = 0.01
A(1,1,2) = 0.0
END IF
JDOF(1,1) = 2
JDOF(1,2) = 2
RETURN
END
Please help.
Thanks.
Joe
In the odb file I can see the FV changing from 0 to 1 but I am not sure how to use the field in the mpc subroutine.
In the subroutine mpc then field looks different to that in usdfld and ufield:
USDFLD Field(1)
UFIELD Field(1,1)
MPC Field(1,1,1)
like above the field has different dimensions and it doesn't seem to match in each subroutine.
In my MPC subroutine when I try to do as below, it gives me an error saying that "The number of subscripts is incorrect.":
IF (Field(1).LT.1) THEN
UE(1) = U(2,1)
A(1,1,1) = 1.
A(1,1,2) = 0.0
ELSE
LMPC=0.0
UE(1) = 0.0
A(1,1,1) = 0.01
A(1,1,2) = 0.0
END IF
JDOF(1,1) = 2
JDOF(1,2) = 2
RETURN
END
Please help.
Thanks.
Joe