d0ab6854-0abc-48de-9d35-4f86a037d0b6 3 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

xyz# pvtm# z.f# BCy# PilesPVTM PilesXYZ.seat PVTM.group Z.p_fix BC.y x y e.x e.y x' y' I.x I.y I.xy θ I.x' I.y' x'1 y'1 PVTM.pile PVTM.pile_corner i N e.x' e.y' P.z V.x V.y T.z M.y M.x V.x' V.y' M.x' M.y' Vx' Vy' Mx' My' Pz Tz ID.XmaxYmax ID.XmaxYmin ID.XminYmax ID.XminYmin Clear PilesXYZ.seat xyz# : PVTM.group pvtm# : P.z PVTM.group 1 el 1 * 0.01 N * + : V.x PVTM.group 2 el 1 * 0.02 N * + : V.y PVTM.group 3 el 1 * 0.03 N * + : T.z PVTM.group 4 el 1 * 0.04 N * m * + : M.y PVTM.group 5 el 1 - (* 0.05 N * m * + : M.x PVTM.group 6 el 1 * 0.06 N * m * + : Z.p_fix z.f# : BC.y BCy# : x PilesXYZ.seat 1 col : y PilesXYZ.seat 2 col : N y length : e.x x sum N / : e.y y sum N / : x' x e.x - : y' y e.y - : I.x y' i el (2^(i 1 N sum : I.y x' i el (2^(i 1 N sum : I.x 1 mm 2^* I.x + : I.y 2 mm 2^* I.y + : I.xy x' i el y' i el * (i 1 N sum : θ 12 / 2 I.xy * I.y I.x - / atan * : I.x' I.x θ cos (2^* I.y θ sin (2^* + 2 I.xy * θ sin * θ cos * - : I.y' I.x θ sin (2^* I.y θ cos (2^* + 2 I.xy * θ sin * θ cos * + : I.z x' i el (2^y' i el (2^+ (i 1 N sum : x'1 x' θ cos * y' θ sin * + : y'1 y' θ cos * x' θ sin * - : e.x' e.x θ cos * e.y θ sin * + : e.y' e.y θ cos * e.x θ sin * - : V.x' V.x θ cos * V.y θ sin * - : V.y' V.x θ sin * V.y θ cos * + : M.x' M.x θ cos * M.y θ sin * + : M.y' M.x θ sin * - M.y θ cos * + : Pz N Zeros 1 + : Vx' N Zeros 1 + : Vy' N Zeros 1 + : Tz N Zeros 1 + : My' N Zeros 1 + : Mx' N Zeros 1 + : Tz Tz 0.000001 * kN * m * : Vx' V.x' N / y' T.z I.z / * - : Vy' V.y' N / x' T.z I.z / * + : Mx' Vy' Z.p_fix 2 / * : Pz P.z N / x'1 P.z e.x' * M.y' - V.x' Z.p_fix * 2 / + I.y' / * + y'1 P.z e.y' * M.x' + V.y' Z.p_fix * 2 / + I.x' / * + : BC.y pinned_c ≡ My' M.y' N / My' * : My' Vx' Z.p_fix 2 / * : if PVTM.pile 1 N Zeros 1 + : j 1 N range PVTM.pile 1 j el 6 Zeros 1 + : PVTM.pile 1 j el Pz j el Vx' j el Vy' j el Tz j el My' j el Mx' j el 61 mat : 21 line for PVTM.pile_corner 24 Zeros 1 + : ID.XmaxYmin 0 : i 1 N range x'1 i el m / 3 round x'1 Max m / 3 round ≡ y'1 i el m / 3 round y'1 Min m / 3 round ≡ ID.XmaxYmin i : continue if continue if for ID.XmaxYmax 0 : i 1 N range x'1 i el m / 3 round x'1 Max m / 3 round ≡ y'1 i el m / 3 round y'1 Max m / 3 round ≡ ID.XmaxYmax i : continue if continue if for ID.XminYmax 0 : i 1 N range x'1 i el m / 3 round x'1 Min m / 3 round ≡ y'1 i el m / 3 round y'1 Max m / 3 round ≡ ID.XminYmax i : continue if continue if for ID.XminYmin 0 : i 1 N range x'1 i el m / 3 round x'1 Min m / 3 round ≡ y'1 i el m / 3 round y'1 Min m / 3 round ≡ ID.XminYmin i : continue if continue if for PVTM.pile_corner 11 el ID.XmaxYmin kN * : PVTM.pile_corner 21 el PVTM.pile ID.XmaxYmin el : PVTM.pile_corner 12 el ID.XmaxYmax kN * : PVTM.pile_corner 22 el PVTM.pile ID.XmaxYmax el : PVTM.pile_corner 13 el ID.XminYmax kN * : PVTM.pile_corner 23 el PVTM.pile ID.XminYmax el : PVTM.pile_corner 14 el ID.XminYmin kN * : PVTM.pile_corner 24 el PVTM.pile ID.XminYmin el : PilesPVTMresults# PVTM.pile_corner : PilesXYZ.seat PVTM.group Z.p_fix BC.y x y e.x e.y x' y' I.x I.y I.xy θ I.x' I.y' x'1 y'1 PVTM.pile PVTM.pile_corner i N e.x' e.y' P.z V.x V.y T.z M.y M.x V.x' V.y' M.x' M.y' Vx' Vy' Mx' My' Pz Tz ID.XmaxYmax ID.XmaxYmin ID.XminYmax ID.XminYmin Clear PilesPVTMresults# 1001 line :

coordinates of abutment seat piles are defined below. Some input information:

]]>

abutment seat skew angle

θ.skew 19 ° * :

L.seat 33480 mm * θ.skew cos / : 35.409 m *

Length of abutment seat

distance from seat edge to CL of first pile

L.edge 1000 mm * :

number of pile rows along the X' axis

No.p_x' 16 :

number of pile rows along the Y' axis

No.p_y' 2 :

total number of piles

No.p No.p_x' No.p_y' * : 32

Pile spacing along X' axis - piles are centered on abutment seat CL

sp.p_x' 2000 mm * :

Lets calculate the pile spacing along Y' axis:

sp.p_y' L.seat 2 L.edge * - No.p_x' 1 - / : mm 2227.276

Now we will define the pile coordinates in a system X'Y'Z:

coordinate of the 1st pile:

Pile1.X' sp.p_x' 2 / : 1 m *

Pile1.Y' L.seat 2 / L.edge - (- : 16.705 m * -

i 1 No.p range i No.p_x' ≤ PilesX'Y'Z.seat i 1 el Pile1.X' : PilesX'Y'Z.seat i 2 el Pile1.Y' i 1 - (sp.p_y' * + : 21 line PilesX'Y'Z.seat i 1 el Pile1.X' sp.p_x' - : PilesX'Y'Z.seat i 2 el PilesX'Y'Z.seat No.p_x' 2 el i No.p_x' - 1 - (sp.p_y' * - : 21 line if for

Pile coordinates; both X'Y'Z and XYZ:

rot.z θ.skew cos θ.skew sin - θ.skew sin θ.skew cos 22 mat : 0.946 0.326 - 0.326 0.946 22 mat

rotation matrix

PilesX'Y'Z.seat 1 m * 16.7046 m * - 1 m * 14.4773 m * - 1 m * 12.25 m * - 1 m * 10.0227 m * - 1 m * 7.7955 m * - 1 m * 5.5682 m * - 1 m * 3.3409 m * - 1 m * 1.1136 m * - 1 m * 1.1136 m * 1 m * 3.3409 m * 1 m * 5.5682 m * 1 m * 7.7955 m * 1 m * 10.0227 m * 1 m * 12.25 m * 1 m * 14.4773 m * 1 m * 16.7046 m * 1 m * - 16.7046 m * 1 m * - 14.4773 m * 1 m * - 12.25 m * 1 m * - 10.0227 m * 1 m * - 7.7955 m * 1 m * - 5.5682 m * 1 m * - 3.3409 m * 1 m * - 1.1136 m * 1 m * - 1.1136 m * - 1 m * - 3.3409 m * - 1 m * - 5.5682 m * - 1 m * - 7.7955 m * - 1 m * - 10.0227 m * - 1 m * - 12.25 m * - 1 m * - 14.4773 m * - 1 m * - 16.7046 m * - 322 mat

PilesXYZ.seat PilesX'Y'Z.seat rot.z * : 4.493 m * - 16.12 m * - 3.768 m * - 14.014 m * - 3.043 m * - 11.908 m * - 2.318 m * - 9.802 m * - 1.592 m * - 7.696 m * - 0.867 m * - 5.59 m * - 0.142 m * - 3.484 m * - 0.583 m * 1.379 m * - 1.308 m * 0.727 m * 2.033 m * 2.833 m * 2.758 m * 4.939 m * 3.483 m * 7.045 m * 4.209 m * 9.151 m * 4.934 m * 11.257 m * 5.659 m * 13.363 m * 6.384 m * 15.469 m * 4.493 m * 16.12 m * 3.768 m * 14.014 m * 3.043 m * 11.908 m * 2.318 m * 9.802 m * 1.592 m * 7.696 m * 0.867 m * 5.59 m * 0.142 m * 3.484 m * 0.583 m * - 1.379 m * 1.308 m * - 0.727 m * - 2.033 m * - 2.833 m * - 2.758 m * - 4.939 m * - 3.483 m * - 7.045 m * - 4.209 m * - 9.151 m * - 4.934 m * - 11.257 m * - 5.659 m * - 13.363 m * - 6.384 m * - 15.469 m * - 322 mat

PilesXYZ.seat 3 m * 3 m * - 3 m * 3 m * 3 m * - 3 m * 3 m * - 3 m * - 42 mat :

if i delete this, the function will use larger matrix and will crash!

PilesXYZ.seat 33763561232633330 m * 13815 m * 38427 m * 16 mat kN * 5 m * frame PilesPVTM kN 1234687915429550 m * 3854 m * 2386 m * 61 mat 1464312649550 m * 3161 m * 2386 m * 61 mat 1000212646770 m * 3161 m * 1693 m * 61 mat 223815426770 m * 3854 m * 1693 m * 61 mat 24 mat