Free edges between areas

I’m modelling two walls which share an edge. That edge should not transfer shear between the walls.
I suspect I need the “SARB FIX” syntax, but I can’t seem to get it to work?

Please see the following code sample:

+prog Sofimshc urs:21.1
SYST 3D GDIV 100000 GDIR NEGZ
CTRL MESH 1

SPT 1 X 0       Y 0    Z 0
SPT 2 X 1000    Y 0    Z 0
SPT 3 X 1000    Y 0    Z 1000
SPT 4 X 0       Y 0    Z 1000

SPT 5 X 2000    Y 0    Z 0
SPT 6 X 2000    Y 0    Z 1000

SAR 1 MNO 11  GRP 11  T 100
    SARB TYPE OUT NA 1 NE 2
    SARB TYPE OUT NA 2 NE 3 ! shared edge
    SARB TYPE OUT NA 3 NE 4
    SARB TYPE OUT NA 4 NE 1

SAR 2 MNO 11  GRP 11  T 100
    SARB TYPE OUT NA 2 NE 5
    SARB TYPE OUT NA 5 NE 6
    SARB TYPE OUT NA 6 NE 3
    SARB TYPE OUT NA 3 NE 2 ! shared edge

SLN 1 NPA 1 NPE 5 FIX F
END

+prog sofiload urs:50.2
ACT 'G' GAMU 1

LC 1 TYPE 'G' TITL 'test'
    LINE AUTO PROJ N TYPE PZZ P1 -10  X1 1000   Y1 0 Z1 1000 $$
                              P2 -10  X2 2000   Y2 0 Z2 1000
end

+PROG ASE urs:51
LC ALL
END 

Thank you

Ah I figures it out:

SARB TYPE OUT NA 3 NE 2 FIX PX

Is there any way to make PX unidirectional though?

I’m guessing you mean e.g. force in only positive x-direction not negative.

I think you will have to define a non-linear spring between the elements.
I suggest you model the boundary as two seperate sln’s (use distinct spt’s for each line), e.g.:
Spt 1 x 0 0 0
Spt 2 x 1 0 0
Spt 3 x 0 0 0
Spt 4 x 1 0 0
Sln 1 Np 1 2
Sln 2 Np 3 4
And then use Sarb NL to use the lines
Connect the two with slns.
Maybe check out coupling_beam_01.dat (open teddy-> file -> examples -> sofimshc)