Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant 1 x z y Q Descrizione del problema di de Saint Venant G V S -solido cilindrico con una dimensione predominante; -materiale omogeneo, elastico, lineare, isotropo; -assenza di forze di volume; -distribuzioni di sforzi applicati solo sulle basi.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
1
x
zy
Q
Descrizione del problema di de Saint Venant
GV
S
-solido cilindrico con una dimensione predominante;-materiale omogeneo, elastico, lineare, isotropo;-assenza di forze di volume;-distribuzioni di sforzi applicati solo sulle basi.
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
2
x z
y
Q
carichi applicati solo sulle basi del cilindro
B
N
Mx
Ty
My
Mz
Tx
B
N
Mz
Ty My
Mx
Tx
risultanti delle distribuzioni di tensioni superficiali
N
Mz
TyMy
Mx
Tx
0
0
0
0
0 0
BB
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
3
Equazioni governanti [1/2]
Equilibrio
( )0
1 00
, , ,
, , ,
, , ,
x x yx y zx z
xy x y y zy z
xz x yz y z z
aσ τ ττ σ ττ τ σ
⎧ + + =⎪ + + =⎨⎪ + + =⎩
in V
( )0 0
1 0 00
0
zx x
n zy y
xz yz z
xz x yz y
nb n
n n
ττ
τ τ σ
τ τ
⎡ ⎤ ⎡ ⎤⎢ ⎥ ⎢ ⎥= = ⇒⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥⎣ ⎦⎣ ⎦⇒ + =
t 0
su S
( ) ( )
11
1 122 1
1 1
xy xyx x y z
y x y z yz yz
z x y z zx zx
GE E EEG
E E E G
E E E G
ν ν γ τε σ σ σ
ν νε σ σ σ γ τν
ν νε σ σ σ γ τ
⎧⎧ == + − − ⎪⎪⎪⎪⎪⎪ = − + − = =⎨⎨
+⎪⎪⎪⎪ = − − + =⎪⎪ ⎩⎩
Legame costitutivo:
in V
x z
y
Q
V
S
n
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
4
Congruenza
( )( )Trot rot =ε 0
N.B.: sono condizioni sufficienti solo per domini semplicemente connessi
( ) ( )( )( )
3 2
2
2
, , ,
, , ,
, , ,
, , , , ,
, , , , ,
, , , , ,
x yy y xx xy xy
y zz z yy yz yz
z xx x zz zx zx
x yz xy z zx y yz x x
y zx yz x xy z zx y y
z xy zx y yz x xy z z
ε ε γ
ε ε γ
ε ε γ
ε γ γ γ
ε γ γ γ
ε γ γ γ
⇓
+ =⎧⎪
+ =⎪⎪ + =⎪⎪⎨ = + −⎪⎪ = + −⎪⎪ = + −⎪⎩
in V
x z
y
Q
V
S
n
Equazioni governanti [2/2]
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
5
Metodo semi-inverso: si assume un andamento noto per alcuni campi, salvo poi verificare che questo sia coerente con la soluzione
Per le equazioni di equilibrio si ricava:
( )( )
, , ,
,,
zx zx
zy zy
z z xz x yz y
x yx y
τ ττ τσ τ τ
⎧ =⎪⇒ =⎨⎪ = − −⎩
0x y xyσ σ τ≡ ≡ ≡Si assume: in V
( )0
1 00
,
,
, , ,
zx z
zy z
xz x yz y z z
aτττ τ σ
⎧ =⎪ =⎨⎪ + + =⎩
( )( )( ) ( ), ,
,,
,
zx zx
zy zy
z xz x yz y
x yx y
z h x y
τ ττ τ
σ τ τ
⎧ =⎪⎪⇒ =⎨⎪ = − + +⎪⎩
mentre le equazioni costitutive divengono: ( ) ( )
( )
012
1 1
,
,
x zxy
y z zx zx
z z zy zy
E
x yE G
x yE G
νε σ γνε σ γ τ
ε σ γ τ
⎧ = −⎪ =⎧⎪ ⎪⎪ ⎪= − =⎨ ⎨⎪ ⎪⎪ ⎪= =⎪ ⎩⎩
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
6
( )( )( )( ) ( )
1
, ,
,,
,
zx zx
zy zy
z zx x zy y
x ya x y
z h x y
τ ττ τ
σ τ τ
⎧ =⎪⎪ =⎨⎪ = − + +⎪⎩
( ) ( )
( )
012
1 1
,
,
x zxy
y z zx zx
z z zy zy
E
x yE G
x yE G
νε σ γνε σ γ τ
ε σ γ τ
⎧ = −⎪ =⎧⎪ ⎪⎪ ⎪= − =⎨ ⎨⎪ ⎪⎪ ⎪= =⎪ ⎩⎩
( ) ( )( )( )
3 2
2
2
, , ,
, , ,
, , ,
, , , , ,
, , , , ,
, , , , ,
x yy y xx xy xy
y zz z yy yz yz
z xx x zz zx zx
x yz xy z zx y yz x x
y zx yz x xy z zx y y
z xy zx y yz x xy z z
ε ε γ
ε ε γ
ε ε γ
ε γ γ γ
ε γ γ γ
ε γ γ γ
+ =⎧⎪
+ =⎪⎪ + =⎪⎪⎨ = + −⎪⎪ = + −⎪⎪ = + −⎪⎩
sostituendo le precedentinelle equazioni di congruenza
si ottengonole seguenti relazioni: ( ) ( )
( )
000
32
20
,
,
,
, , , ,
, , , ,
,
z xx
z yy
z zz
z yz zx y zy x x
z zx zx y zy x y
z xy
E
E
σσσ
σ γ γν
σ γ γν
σ
=⎧⎪ =⎪⎪ =⎪⎪
= − −⎨⎪⎪
= −⎪⎪
=⎪⎩
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
7
( )( )
( )
00003
2
2
,
,
,
,
, , , ,
, , , ,
z xx
z yy
z zz
z xy
z yz zx y zy x x
z zx zx y zy x y
E
E
σσσσ
σ γ γν
σ γ γν
=⎧⎪ =⎪⎪ =⎪⎪ =⎨⎪ = − −⎪⎪⎪ = −⎪⎩
Lo sforzo σz può assumere solo l’espressione:
( )z a by cx z d ey fxσ = + + + + +
( )
( )
( )
( )3 2
2
, , , ,
, , , ,
z
z yz zx y zy x x
z zx zx y zy x y
a by cx z d ey fxEe
Ef
σ
σ γ γν
σ γ γν
⎧ = + + + + +⎪⎪ = = − −⎨⎪⎪ = = −⎩
Il sistema è equivalente al seguente:
( )2 , ,zx y zy xEe x f y m γ γν
− + + = −
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
8
( )( )
( )3 2
, ,
z
zx y zy x
a by cx z d ey fx
e x f y mE
σνγ γ
⎧ = + + + + +⎪⎨
− = − + +⎪⎩
L’equilibrio globale richiede:
0 0
0 0 0
0 0 0
cos . cos .cos .cos .
B Bz z
x x x x y
y y y y x
N N t M M tT T t M M T zT T t M M T z
⎧ ⎧= = = =⎪ ⎪= = = +⎨ ⎨⎪ ⎪= = = −⎩ ⎩
l coefficienti possono essere legati alle azioni interne.
N
Mz
Ty My
Mx
TxN
Mz
TyMy
Mx
Tx
0
0
0
0
0 0
z
B B
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
9
( ) 0
0
0
z x y x yA
x y
x y
N dA aA bS cS z d A eS f S N
aA bS cS Nd A eS f S
σ= = + + + + + ≡ ⇒
+ + =⎧⇒ ⎨ + + =⎩
∫
( ) 0 0
0
0
x z x x xy x x xy x yA
x x xy x
x x xy y
M y dA aS bJ cJ z d S e J f J M zT
aS bJ cJ Md S e J f J T
σ= = + + + + + ≡ + ⇒
⎧ + + =⎪⇒ ⎨ + + =⎪⎩
∫
( ) 0 0
0
0
y z y xy y y xy y y xA
y xy y y
y xy y x
M x dA aS bJ cJ z d S eJ f J M zT
aS bJ cJ Md S eJ f J T
σ− = = + + + + + ≡ − + ⇒
⎧ + + = −⎪⇒ ⎨ + + =⎪⎩
∫
00 0
0 0
0,
x y x y
x x xy x x x xy y
y xy y y y xy y x
A S S a N A S S dS J J b M S J J e TS J J c M S J J f T
Con assi baricentrici (“centrali”, Sx=Sy=0)e principali d’inerzia (Jxy=0):
( ) ( )1 1 1 1 2 2, ,
, ,
, , , ,
yxz
x y
yxx y x y zx y zy x
y x
yxxz x yz y
y x
MMN y xA J J
TTdiag A J J diag A J J y x mE J J E
TT x yJ J
σ
ν νγ γ
τ τ
− − − −
⎧= + −⎪
⎪⎪ ⎛ ⎞⎪⎡ ⎤⎡ ⎤= ⇒ = ⇒ − = − +⎜ ⎟⎨⎣ ⎦ ⎣ ⎦ ⎜ ⎟
⎝ ⎠⎪⎪ ⎛ ⎞⎪ + = − +⎜ ⎟⎪ ⎝ ⎠⎩
J J
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
13
Postulato di de Saint Venant
La soluzione così determinata non soddisfa, in generale, le equazioni locali di equilibrioin corrispondenza delle basi del cilindro.
Peraltro, le distribuzioni di sforzi su tali basi sono quasi sempre incognite.
Si osserva che le distribuzioni di sforzi trovate sono in equilibriocon i risultanti di quelle applicate sulle basi.
Di conseguenza, le distribuzioni di sforzi ottenute come differenza tra quelle trovate, calcolate sulle basi del cilindro, e quelle realmente applicate su tali basi sono auto-equilibrate.
N
Mz
Ty My
Mx
Tx=
N=0
Mz=0
Ty=0
My=0
Mx=0
Tx=0N
Mz
Ty My
Mx
Tx+
Si assume che l’effetto di distribuzioni auto-equilibrate di sforzi sia significativo
solo a una distanza dalle basi inferiore alla dimensione massima della sezione.
Tale ipotesi è ben verificata sperimentalmente nella maggior parte delle applicazionidella meccanica dei solidi.
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
19
Interpretazione cinematica delle costanti di integrazione
( )
( )
0
0
,
,
u y z ny m zy r z qE
v x z nx m zx s z pE
ν
ν
⎧ = + + + +⎪⎪⎨⎪ = − − + +⎪⎩
( )
( )
0
0
,
,
u y z n m z y r z qE
v x z n m z x s z pE
ν
ν
⎧ ⎛ ⎞= + + + =⎜ ⎟⎪ ⎝ ⎠⎪⎪⎪⎨⎪ ⎛ ⎞⎪ = − + + + =⎜ ⎟
⎝ ⎠⎪⎪⎩
Q
y, v
x, uz, w
ϕz
ϕz
CIR=C(z)
u
v
se cost.z z= =
( ) ( )
( ) ( )
z Q
z Q
z y u z
z x v z
ϕ
ϕ
= − +
= + +
( )( )
0 0
cost.
z z
zz
nd z m
dz E
ϕ ϕϕ νϑ
⎧ = = −⎪⇒ ⎨
= = − =⎪⎩
( ) ( )( ) ( )
00
00
,
,z z
z z
u y z z y r z q
v x z z x s z p
ϕ ϑ
ϕ ϑ
⎧ = − + + +⎪⇒ ⎨= + + + +⎪⎩
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
20
( ) ( )( ) ( )
( ) ( )( ) ( )
0 0 0
0 0 0
,
,
z Q
Qz
z z Q y
z u z
z z Q x
v zz
u y z z y u z
v x z z x v z
ϕ
ϕ
ϕ ϑ ϕ
ϕ ϑ ϕ
⎧ = − + + +⎪⎪⎪⇒ ⎨⎪ = + + + −⎪⎪⎩
Qy, v
x, u
z, w
ϕy
( ) ( )( ) ( )
00
0
,
,z z
y
u y z r y z y q
y z u y
ϑ ϕ
ϕ
= − − + =
= +
u
Q
y, v
x, uz, w
( ) ( )( ) ( )
00
0
,
,z z
x
v x z s x z x p
x z v x
ϑ ϕ
ϕ
= + + + =
= − +
v
ϕx
per si ha:cost.x x= = per si ha:cost.y y= =
( ) ( )
( ) ( )
0
0
,
,
x z x z
y z y z
dv x zx s x xdz
du y zy r y y
dz
ϕ ϑ ϕ ϑ
ϕ ϑ ϕ ϑ
⎧= − = − − = −⎪⎪
⎨⎪ = + = + − = −⎪⎩
NB: descrivono una roto-traslazione rigida del solido cilindrico0 0 0 0 0, , , ,x y z Q Qu vϕ ϕ ϕ
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
21
Osservazione sui Centri di Istantanea Rotazione (CIR)
Q
y, v
x, uz, w
ϕz
ϕz
CIR=C(z)
u
v
( ) ( )( ) ( )
( ) ( )( ) ( )
0 0 0
0 0 0
,
,
z Q
Qz
z z Q y
z u z
z z Q x
v zz
u y z z y u z
v x z z x v z
ϕ
ϕ
ϕ ϑ ϕ
ϕ ϑ ϕ
⎧ = − + + +⎪⎪⎪⎨⎪ = + + + −⎪⎪⎩
la posizione di C si ottiene dal sistema:
( ) ( )
( ) ( )
0 0 0
0 0 0
0
0
,
,
C z z C Q y
C z z C Q x
u y z z y u z
v x z z x v z
ϕ ϑ ϕ
ϕ ϑ ϕ
⎧ = = − + + +⎪⎨⎪ = = + + + −⎩
( )0 00 0
0 0, :Q yQ xC C
z z z z
u zv zx y C z
z zϕϕ
ϕ ϑ ϕ ϑ⎧ ⎫+−
⇒ = − = ⇒⎨ ⎬+ +⎩ ⎭
la posizione del CIR varia con z
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
22
Esempio: lamina rettangolare [1/2]
Q≡G
y, v
x, uz, w
Mt
0 0 0
0 0
00
x y z
G Gu vϕ ϕ ϕ⎧ = = =⎨
= =⎩Se si assumono:
(roto-traslazione
nulla della sezione in z=0):
Mt Mt
y, v
z, w
b/2 b/2
( )( )
,,
z
z
u y z z yv x z z x
ϑϑ
⎧ = −⎨
= +⎩
gli spostamentisono i seguenti: G ≡ C
v
y, v
x, u
b2z zϕ =ϑ
Mt
Mt
x
y z
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
23
0 0 0
0 0
00 0;
x y z
G Gu vϕ ϕ ϕ⎧ = = =⎨
= ≠⎩
se, invece, si aggiunge una traslazione verticale durante il processo deformativo:
( ) ( )( ) ( ) 0
,,
z
z G
u y z z yv x z z x v
ϑϑ
⎧ = −⎨ = + +⎩
gli spostamentisono i seguenti:
Esempio: lamina rettangolare [2/2]
0 1 0,GC C
z
vx yzϑ
⎧ ⎫= − =⎨ ⎬
⎩ ⎭e si ha:
G
y, v
x, uv
C(z)
b2z zϕ =ϑ
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
24
Le posizioni dei CIR delle varie sezioni dipendono dalle roto-traslazioni rigide del solido;è però sempre possibile definire una opportuna roto-traslazione del solido,in modo da ottenere un unico CIR, arbitrariamente scelto, comune a tutte le sezioni:
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
0 0 0 0 0 0
0 0 0 0 0 0
0
0
,
,
C z z C Q y Q z C y z C
C z z C Q x Q z C x z C
u y z z y u z u y y z
z
v x z z x v z v x x z
ϕ ϑ ϕ ϕ ϕ ϑ
ϕ ϑ ϕ ϕ ϕ ϑ
⎧ = − + + + = − + − =⎪⎪
∀⎨⎪
= + + + − = + − − =⎪⎩
( ) ( )( )
( )
( ) ( ) ( )
0
0
,
,z
z z C
z
z z C
u y z z y y
v x z z x xϕ
ϕ ϑ
ϕ ϑ
⎧ = − + −⎪
⇒ ⎨⎪ = + + −⎩
0 0
0
0 0
0
Q z C
y z C
Q z C
x z C
u y
y
v x
x
ϕ
ϕ ϑ
ϕ
ϕ ϑ
⎧ =⎪
=⎪⎨
= −⎪⎪ =⎩
ne consegue:
Q
y, v
x, u z, w
ϕz
ϕz
C
u
vNB: la posizione di C è definibile in modo arbitrario
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
25
( ) ( ) ( )( ) ( ) ( )
( )( ) ( ) ( )( ) ( ) ( )
0
0
00
, ,
,
,,
, , , ,
, , , ,
z z C
z z C
zx zx x z
zy zy y z
xz x yz y
x y xy
u y z z y y
v x z z x xw w x y
x y G x y G w u
x y G x y G w v
ϕ ϑ
ϕ ϑ
τ γ
τ γ
τ τσ σ τ
⎧ = − + −⎪
= + + −⎪⎪ =⎪⎪ = = +⎨⎪ = = +⎪⎪ + =⎪
= = =⎪⎩
La soluzione si completa con le restanti equazioni:
( ) ( ) ( )( ) ( ) ( )
, , , ,
, , , ,zx x z x z C
zy y z y z C
x y G w u G w y y
x y G w v G w x x
τ ϑ
τ ϑ
⎧ ⎡ ⎤= + = − −⎣ ⎦⎪⇒ ⎨⎡ ⎤= + = + −⎪ ⎣ ⎦⎩
( ) ( ),,z
w x yx y
ϑΨsi definisce, per ϑz ≠ 0, la funzione di ingobbamento Ψ:
( ) ( ), ,zw x y x yϑ⇒ = Ψ( ) ( )( ) ( )
, ,
, ,zx z x C
zy z y C
x y G y y
x y G x x
τ ϑ
τ ϑ
⎧ ⎡ ⎤= Ψ − −⎣ ⎦⎪⇒ ⎨⎡ ⎤= Ψ + −⎪ ⎣ ⎦⎩
Inoltre, l’equilibrio indefinito fornisce:
0, , , ,xz x yz y z xx yy zG Gτ τ ϑ ϑ⎡ ⎤+ = Ψ +Ψ = ∆Ψ =⎣ ⎦
( ) ( ), , , , , ,xz x yz y z x C x z y C yG y y G x xτ τ ϑ ϑ ⎡ ⎤⎡ ⎤+ = Ψ − − + Ψ + −⎣ ⎦ ⎣ ⎦
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
26
( ) ( )( ) ( )
, ,
, ,zx z x C
zy z y C
x y G y y
x y G x x
τ ϑ
τ ϑ
⎧ ⎡ ⎤= Ψ − −⎣ ⎦⎪⎨
⎡ ⎤= Ψ + −⎪ ⎣ ⎦⎩
In ultimo, l’equilibrio (1b) sul contorno Γ fornisce:
0 in A∆Ψ =
( ) ( ) 0, ,n xz x yz y
z x C x z y C y
n n
G y y n G x x n
τ τ τ
ϑ ϑ
= + =
⎡ ⎤⎡ ⎤= Ψ − − + Ψ + − =⎣ ⎦ ⎣ ⎦
, ,x x y yx y n n
n x n y n∂Ψ ∂Ψ ∂ ∂Ψ ∂
= + = Ψ + Ψ∂ ∂ ∂ ∂ ∂
Sfruttando la derivazione di funzioni composte:
x z
y
Q
AΓ
n τP
C
si ricava la condizione al contorno: ( ) ( ) suC x C yy y n x x nn
∂Ψ= − − − Γ
∂
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
27
Dunque, la funzione di ingobbamento è soluzione del problema di Neumann−Dini:
( )( )
( ) ( ) ( )
4 04
4
in
suC x C y
a A
b y y n x x nn
⎧ ∆Ψ =⎪⎨ ∂Ψ
= − − − Γ⎪ ∂⎩
( ),x zx z x CA AT dA G y y dAτ ϑ= = Ψ − −∫ ∫
Si noti che, come atteso, la soluzione del problema (4) fornisce (s ascissa curvilinea di Γ):
( ) ( ), , ,xx yy y yA Ax dA x dA x n ds
ΓΨ = −Ψ = −Ψ∫ ∫ ∫
( ), ,x z x x y y C xT G x n n y y n dsϑΓ
⎡ ⎤= Ψ + Ψ − −⎣ ⎦∫
( ), ,z x C x z xxAG x y y n ds G x dAϑ ϑ
Γ⎡ ⎤= Ψ − − − Ψ⎣ ⎦∫ ∫
Utilizzando (4a):
( )
( ) ( ) ( ) ( )2 0
z C x
z C y z C x z C
G x y y n dsn
G x x x n ds G x x x t ds G x xx dx
ϑ
ϑ ϑ ϑ
Γ
Γ Γ Γ
∂Ψ⎡ ⎤= − − =⎢ ⎥∂⎣ ⎦
= − − = − = − ≡
∫
∫ ∫ ∫
(4b)
( ) ( )2 0,y zy z y C z CA AT dA G x x dA G y yy dyτ ϑ ϑ
Γ= = Ψ + − = − ≡∫ ∫ ∫
Analogo risultatosi ha per Ty:
(4a)
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
28
OSSERVAZIONI
zx zyAy x dAτ τ= − +∫
Da ultimo, per il momento torcente si ricava:
( ) ( )inerzia geometrica torsionale
, ,
t
t z x C y CA
J
M G y y y x x x dAϑ ⎡ ⎤⎡ ⎤= − Ψ − − + Ψ + −⎣ ⎦ ⎣ ⎦∫ ( )4 tz
t
McGJ
ϑ⇒ =
( ) ( )t B zx B zyAM y y x x dAτ τ⎡ ⎤= − − + −⎣ ⎦∫ zx zy B x B yA
y x dA y T x Tτ τ= − + + −∫0 0
1) La soluzione del problema (4) non dipende dal sistema di riferimento scelto.
Infatti, l’operatore di Laplace e la derivata direzionale
non dipendono dal tipo di riferimento utilizzato.
2 2
2 2x y∂ ∂
∆ = +∂ ∂ n
∂∂
Inoltre, il termine noto può scriversi come (k versore dell’asse z):
( ) ( ) ( )C x C yy y n x x n P C
n∂Ψ
= − − − = − − ∧∂
n ki
in cui il vettore e i prodotti vettore (∧) e scalare (•) non dipendono dal riferimento scelto.
xz
y
Q
Γ
nP
C
( )P C−
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
29
2) Esistenza della soluzione del problema (4)
Dall’equazione (4a) si ricava l’identità seguente (teorema della divergenza):
0 0A
dA ds dsnΓ Γ
∂Ψ∆Ψ = ⇒ ∆Ψ = ∇Ψ = =
∂∫ ∫ ∫ni
( ) ( ) ( ) 0C x C yds y y n x x n ds P C dsnΓ Γ Γ
∂Ψ= − − − = − − ∧ =
∂∫ ∫ ∫ n ki
L’esistenza della soluzione è allora subordinata al rispetto della condizione:
( ) ( ) ( ) ( ) ( )d P CP C P C P C P Cds−
−− − ∧ = − − ∧ = − = −i i i i
tn k n k t
Utilizzando le proprietà del prodotto misto
e la definizione di t (versore tangente):
( ) ( ) ( ) ( )0
1 02
is L
s
d P Cds P C ds P C P Cn ds
=
Γ Γ=
∂Ψ −= − = − − ≡
∂∫ ∫ i i
si ottiene (Li lunghezza del contorno i-esimo):
xz
y
Q
nP
C
t
s
s
Γ
s
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
30
3) Molteplicità della soluzione del problema (4)
La funzione Ψ(x,y) che risolve il problema (4) è definita a meno di una costante arbitraria:
( )
( )
( ) ( ) ( )
( ) ( ) ( )
4 0
44
4
ˆ inˆ
su
ˆˆ ˆ ˆ; , ,ˆ
C x C y
tz z t x C y C tA
t
a A
b y y n x x nn n
Mc J y y y x x x dA JGJ
ϑ ϑ
⎧ ∆Ψ = ∆Ψ =⎪
∂Ψ ∂Ψ⎪⎪ = = − − − Γ⇒ ⎨ ∂ ∂⎪
⎡ ⎤ ⎡ ⎤⎪ = = = − Ψ − − + Ψ + − =⎣ ⎦ ⎣ ⎦⎪⎩∫
se si pone: ( ) ( )ˆ , ,x y x y kΨ = Ψ + ⇒
La costante k concorre a definire una traslazione arbitraria del solido in direzione dell’asse z:
( ) ( ) ( ) ( ) 0ˆˆ , , , ,z z z Qw x y x y x y k w x y wϑ ϑ ϑ= Ψ = Ψ + = +
Essendo il solido non vincolato si può, ad esempio, scegliere k in modo da annullarelo spostamento medio:
( ) ( )10ˆ , ,A Aw x y dA k x y dA
A= ⇒ = − Ψ∫ ∫
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
31
4) Dipendenza della soluzione dal centro di istantanea rotazione C [1/3]
La funzione Ψ(x,y) dipende esplicitamentedalla scelta (arbitraria) del punto C.
( )
( ) ( ) ( )
4 0
4
in
suC x C y
a A
b y y n x x nn
⎧ ∆Ψ =⎪⎨ ∂Ψ
= − − − Γ⎪ ∂⎩
Sia Ψ'(x,y) la soluzione associata al punto C' ≠ C.( )
( ) ( ) ( )
4 0
4
in
suC x C y
a A
b y y n x x nn ′ ′
′⎧ ∆Ψ =⎪
′⎨ ∂Ψ= − − − Γ⎪ ∂⎩
Si consideri la funzione F = Ψ'(x,y) − Ψ(x,y) + (yC' − yC) x − (xC' − xC) y
( ) ( )
0 0
0
in
suC C C C
F AF x yy y x xn n n n n′ ′
′∆ = ∆Ψ − ∆Ψ + =⎧⎪
′∂ ∂Ψ ∂Ψ ∂ ∂⎨= − + − − − = Γ⎪ ∂ ∂ ∂ ∂ ∂⎩
Si ha:
La soluzione è: F=k (=cost.); ne consegue: Ψ'(x,y) = Ψ(x,y) − (yC' − yC) x + (xC' − xC) y + k
inoltre:
I due campi di spostamento w' e w differiscono per una roto-traslazione rigida,con rotazioni attorno agli assi x e y.
Problema di Neumann−Diniomogeneo
( ) ( ) ( ) ( ) ( ), , ,y x
z z C C z C C zw x y x y w x y y y x x x y kα α
ϑ ϑ ϑ ϑ′ ′′ ′= Ψ = − − + − +
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
32
4) Dipendenza della soluzione dal centro di istantanea rotazione C [2/3]
Mt
Mt
x
y z
C
v
y, v
x, u
b2z zϕ =ϑ
−u
Esempio: sezione a T
( )( )( ) ( )
,,
0, y , 0 0
z
z
u y z z yv x z z xw w x
ϑϑ
⎧ = −⎪ = +⎨⎪ = =⎩
Mt
Mt
x
y z C'
v
y, v
x, u
b2z zϕ =ϑ
u ( ) ( )( )( )( )
,,
0, 0,0
y
z C
z
z C
u y z z y yv x z z xw yw x y x
α
ϑϑ
ϑ
′
′
′⎧ = − −⎪ ′ = +⎪⎪
′⎨ =⎪ ′ = −⎪⎪⎩
x
zαy
αy
x
z
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
33
⇒
4) Dipendenza della soluzione dal centro di istantanea rotazione C [3/3]
t x C y C x yA A AJ y y y x x x dA x y dA y x dA⎡ ⎤⎡ ⎤= − Ψ − − + Ψ + − = + + − Ψ + Ψ⎣ ⎦ ⎣ ⎦∫ ∫ ∫
Si definisca, per semplicità operativa, la funzione: ( ) ( ), , C Cx y x y xy yxΨ Ψ + −
Si può allora scrivere:
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
36
( ) ( )2 220, , , ,x y x yA A A
y x dA dA dA− Ψ + Ψ ≡ − Ψ + Ψ = − ∇Ψ ≤∫ ∫ ∫
Se ne deduce la notevole proprietà:
Dunque, per la funzione Ψ(x,y) valgono l’identità e la delimitazione seguenti:
2
t Q QAJ J dA J= − ∇Ψ ≤∫
6) Ulteriori proprietà di Jt [2/3]
Si noti che, considerando il teorema di trasporto per il momento polare d’inerzia
( ), la migliore delimitazione per Jt si ottiene nel caso Q=G:2
Q GJ J A GQ= +
( )2mint G G QA
J J dA J J= − ∇Ψ ≤ =∫
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
37
6) Ulteriori proprietà di Jt [3/3]
Infine, Jt è sempre positivo dato che, sfruttando l’identità ricavata in precedenza:
2 2 , ,t x yAJ x y y x dA= + − Ψ + Ψ =∫
( ) ( )2 2, , 0y xA
x y dA= + Ψ + −Ψ >∫Quest’ultima espressione è necessariamente positiva, poiché il suo annullamento
richiederebbe l’esistenza di una soluzione per il problema:, , 1, , 1x xy
y yx
yx
⎧ ⎧Ψ = Ψ = +⎪ ⎪⇒⎨ ⎨Ψ = − Ψ = −⎪ ⎪⎩ ⎩
( ) ( )22, , , ,x y x yA A
y x dA dA− Ψ + Ψ ≡ − Ψ + Ψ∫ ∫
si può scrivere:
( ) ( ) ( )222 2
0
, , , , , ,x y x y x yAx y y x y x dA⎡ ⎤= + − Ψ + Ψ + − Ψ + Ψ + Ψ + Ψ =⎣ ⎦∫
( ) ( )222 2 2 , 2 , , ,x y x yAx y y x dA= + − Ψ + Ψ + Ψ + Ψ =∫
e ciò, data la regolarità di , è in palese contraddizione col teorema di Schwarz.( ),x yΨ
Si può allora definire il fattore di torsione come: 1G
t
JqJ
= ≥
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
38
Centro di torsione [1/5]
Mt Mt
Q
y
x
Mt
z Q
y
x z
w(x,y)
βx
La rotazione della sezione attorno all’asse z è ben definita,ma l’ingobbamento non permette di definire le rotazionidella sezione attorno agli assi x e y.
Allora, si utilizzano le rotazioni medie (βx e βy).
Tra le varie definizioni possibili, risulta conveniente adottare quella coerentecon un’interpretazione energetica.
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
39
My
Mx
My
MxQ
y
x z
Centro di torsione [2/5]
σz
Visto che gli enti statici che lavorano per le rotazioni sonole coppie, si considera la distribuzione di sforzi generata nelcaso di flessione deviata.
[ ] 110
z x
y
y x MM
σ −
⎡ ⎤⎢ ⎥= ⎢ ⎥⎢ ⎥−⎣ ⎦
J
Le rotazioni medie (βx e βy) della sezione sono quelle cherendono identicamente soddisfatta la seguente relazione:
( ) ( ), , ,x x y y z x yAM M x y w x y dA M Mβ β σ+ = ∀∫
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
40
Centro di torsione [3/5]
Sostituendo le precedenti relazioni in quest’ultima si ricava:
( )[ ] 110
, ,xx x y y x y z x x yA
yy
M M M M x y y x dA M M MM
ββ β ϑ
β−
⎡ ⎤⎡ ⎤ ⎢ ⎥+ = = Ψ ∀⎡ ⎤ ⎢ ⎥⎣ ⎦ ⎢ ⎥⎣ ⎦ ⎢ ⎥−⎣ ⎦
∫ J
Sia ora una funzione d’ingobbamento soluzione di (4) con C≡Q.( ),x yΨ
( ) ( ) [ ]1
, , , ,C Cx y x y k x y yx
⎡ ⎤⎢ ⎥Ψ = Ψ − − ⎢ ⎥⎢ ⎥⎣ ⎦
Dalla precedente osservazione n.4 si deduce: ( ) ( ) ( ) ( ), , Q C Q Cx y x y y y x x x y kΨ =Ψ − − + − +
per semplicità operativa, e senza perdita di generalità, si sceglie “e” in modo che sia α0=0.
(i = 0 indica il contornoesterno della sezione)
potenzialedegli sforzi
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
47
( ) ( )Bt z zy B zx BA
M M x x y y dAτ τ= = − − −∫
Inoltre, per le azioni interne si ha:
x zxAT dAτ= ∫
x z
y
Q
A Γ
n tP
s
s
s
nt
n
t
Sezione con N fori
1,
i
N
y y i yAi
dx
dA n ds n dsφ φ αΓ Γ
=−
= = = ∑∫ ∫ ∫ = 00
y zyAT dAτ= ∫
1,
i
N
x x i xAi dy
dA n ds n dsφ φ αΓ Γ
== − = − = −∑∫ ∫ ∫ = 0
0
zy zx B zy B zxA A Ax y dA x dA y dAτ τ τ τ
=
= − − +∫ ∫ ∫0 0
( ) ( ), ,x yAx y dAφ φ
=
= − −∫ x y An x n y ds dAφ φ φ φ
Γ= − + + + =∫ ∫
( )( ) ( )
2x y A
P Q h s
n x n y ds dAφ φΓ
− • =
= − + +∫ ∫n
0
0
01
22
2 2 2i
i
N
iAi
AA
dA d dφ α ω α ωΓ Γ
=−+
= − −∑∫ ∫ ∫
12 2
N
t i iAi
M dA Aφ α=
= + ∑∫
dω = ds h(s)/2y
x z
n
s
t
Q
Pds
h(s)
dωMt = 2 (volume di φ estesa ai fori)
φ α1
α2
0
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
48
( ) 2, ,yy xxG m
Eνφ φ⇒ − − =
00
,in
,
susu 1...
zx y
zy x
i i
A
i N
τ φτ φ
φφ α
⎧ +⎡ ⎤ ⎡ ⎤= =⎪ ⎢ ⎥ ⎢ ⎥−⎪ ⎣ ⎦⎣ ⎦⎨= Γ⎪
⎪ = Γ =⎩
τ
Dunque, l’equilibrio è soddisfatto in A e su Γ mediante le assunzioni seguenti:
Inoltre, il momento torcente è espresso tramite la relazione:1
2 2N
t i iAi
M dA Aφ α=
= + ∑∫
φ α1
α2
Si noti che le costanti αi possono essere scelte in modo del tutto arbitrario nel rispetto
dell’equilibrio: esse assumono la funzione di “variabili iperstatiche”.
A questo punto, attraverso le equazioni costitutive, si sostituiscono gli sforzi tangenzialinell’equazione di congruenza:
2, ,, ,
zx y zy xzx y zy x m
G G Eτ τ νγ γ− = − =
1, ,xx yy mνφ φ φ
ν⇒ + = ∆ =
+
L’insieme delle equazioni finora ottenute forma un problema ben posto
solo nel caso in cui la sezione sia semplicemente connessa, cioè
senza fori (N=0; la costante m si determina dall’espressione di Mt).0
10
2su
t A
m
M dA
νφν
φφ
⎧∆ =⎪ +⎪⎨ = Γ⎪
=⎪⎩ ∫
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
49
Infatti, le equazioni di congruenza utilizzate (cioè quelle nella forma “interna”: rot[rot(ε)T]=0)
sono sufficienti solo nel caso N=0.
Per sezioni molteplicemente connesse (N>0) è necessario aggiungere altre equazioni dicongruenza: per ricavarle si utilizza il PLV (come condizione sufficiente di congruenza).
Q
y, v
x, uz, w
ϕz
ϕz
C
u
v
Γi
Γ0
Sistema reale (di cui si vuoleimporre la congruenza)
Sistema equilibrato
Q
y
xz
Γi
Γ0
Mt*
0
1
0
2 2
,in
,
susu 1...
zx y
zy x
i iN
t i iAi
A
i N
M dA A
τ φτ φ
φφ α
φ α
∗ ∗∗
∗ ∗
∗
∗ ∗
∗ ∗ ∗
=
⎧ ⎡ ⎤ ⎡ ⎤+= =⎪ ⎢ ⎥ ⎢ ⎥−⎪ ⎣ ⎦⎣ ⎦
⎪⎪ = Γ⎨= Γ =⎪
⎪⎪ = +⎪⎩
∑∫
τ
ϕz(b) = ϕz + ϑz b0
γzx , γzy = ??
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
50
0
, , , ,
e t z t z t z
i zx zx zy zy y zx x zy y zx x zyV V A
L M M M b
L dV dV b dA
ϕ ϕ ϑ
τ γ τ γ φ γ φ γ φ γ φ γ
∗ ∗ ∗
∗ ∗ ∗ ∗ ∗ ∗
⎧ = − =⎪
⇒⎨⎪ = + = − = −⎩ ∫ ∫ ∫
0
1
0
2 2
susu 1..., , , , con
ei
i it z y zx x zy t iA N
L b L b t i iAi
i NM dA MM dA A
φφ αϑ φ γ φ γ φ α
φ α
∗
∗ ∗∗ ∗ ∗ ∗ ∗ ∗
∗ ∗ ∗
=
⎧ = Γ⎪⎪ = Γ =⇒ = − ∀ ⎨⎪ = +⎪⎩
∫∑∫
, ,, ,t z y zx x zy zx y zy x zx y zy xA AM dA n n ds dAϑ φ γ φ γ φ γ φ γ φ γ φ γ∗ ∗ ∗ ∗ ∗ ∗ ∗
Γ= − = − − −∫ ∫ ∫
Applicando il lemma di Green si ricava:
e sostituendo le espressioni che definiscono φ* e raccogliendo i termini omogenei si ottiene:
( ) ( ) 1
2 2 0, , ,i
N
zx y zy x z i z i zx y zy x iAi
dA A n n dsγ γ ϑ φ α ϑ γ γ φ α∗ ∗ ∗ ∗
Γ=
− + + − − = ∀∑∫ ∫
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
51
Questa relazione deve valere per ogni arbitraria scelta della funzione φ* e dei parametri αi*.
( ) 2 1su ...k
zx y zy x z k kn n ds A k Nγ γ ϑΓ
− = Γ =∫
( ) ( ) 1
2 2 0, , ,i
N
zx y zy x z i z i zx y zy x iAi
dA A n n dsγ γ ϑ φ α ϑ γ γ φ α∗ ∗ ∗ ∗
Γ=
− + + − − = ∀∑∫ ∫
Se si assume αi* = 0 per i=1...N, si ha: ( )2 0, ,zx y zy x zAdAγ γ ϑ φ φ∗ ∗− + = ∀∫
che, per l’arbitrarietà di φ*, implica: 2, , inzx y zy x z Aγ γ ϑ− = −
Ne consegue che deve risultare: ( ) 1
2 0i
N
i z i zx y zy x ii
A n n dsα ϑ γ γ α∗ ∗
Γ=
− − = ∀∑ ∫
Assumendo αi* = 0 per i≠k, e αk* = 1 si ha:
Ricordando che dall’approccio negli spostamenti si era ricavato ϑz = −νm/E, è immediato
identificare la prima equazione ottenuta come quella indefinita di congruenza.
Invece, le ulteriori N equazioni di circuitazione rappresentano relazioni che completano le
condizioni di congruenza nella forma “interna”, per il problema della torsione, nel caso di domini
pluri-connessi.
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
52
Interpretazione delle equazioni di circuitazione [1/2]
Le equazioni di circuitazione sono suscettibili di un’interpretazione cinematica.
Se si esprimono gli scorrimenti angolari in funzione degli spostamenti e si utilizzano le
equazioni che esprimono gli spostamenti u e v nel piano x-y, per il k-esimo contorno si ricava:
( ) ( ) ( ), , , ,k k
zx y zy x x z y y z xn n ds w u n w v n dsγ γΓ Γ
− = + − +∫ ∫
( ) 2 1su ...k
zs
zx y zy x z k kn n ds A k Nγ
γ γ ϑΓ
− = Γ =∫
( ) ( )( )
, ,k k
C
x x y y z
h s
w t w t ds P C dsϑΓ Γ
=
= − − − −∫ ∫ in
( ) ( ), ,k
x y y x z C y z C xw n w n y y n x x n dsϑ ϑΓ
=
= − − − − −∫
y
( ) ( ) 20
2k k
kk
z C
Aw L w
w ds ds
ϑ ωΓ Γ
−−
∂= − −
∂∫ ∫
( ) ( )0 2 2k z k z kw w L A Aϑ ϑ
=
= − + =
( ) ( )0kw L w⇒ =
dωC = ds hC(s)/2 < 0
x z
ns
tC
Pds
hC(s) < 0
dωC
Q
Γk
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
53
Interpretazione delle equazioni di circuitazione [2/2]
Si riconduca la sezione pluri-connessa ad una
semplicemente connessa introducendo N tagli:
Per questo solido le equazioni di congruenza interna (rot[rot(ε)T]=0) sono sufficienti,
Γk
w(Lk)−w(0)
Dunque, ogni equazione di circuitazione, unitamente alla regolarità delle deformazioni
stesse sul dominio, garantisce la continuità degli spostamenti attraverso i tagli e
riconduce la sezione semplicemente connessa (con N tagli) alla sezione iniziale
molteplicemente connessa (con N fori, senza tagli).
ma tra i lembi opposti, separati dai tagli, si possono avere discontinuità negli
spostamenti in direzione assiale w(x,y):
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
Riepilogo delle equazioni per il problema della torsione secondo l’approccio negli sforzi.
Si noti che, secondo un bilancio grossolano, le incognite sono la funzione φ(x,y), le N costantiαk e la curvatura torsionale ϑz; corrispondentemente, nel problema l’equazione di Poisson èarricchita dalle N equazioni di circuitazione e dalla relazione che esprime Mt.
Le equazioni di circuitazione possono essere espresse, via legame costitutivo, in funzione di φ.
( )k
zx y zy xn n dsγ γΓ
−∫
2 1su ...k
z k kds G A k Nnφ ϑ
Γ
∂⇒ = Γ =
∂∫
1k
zx y zy xn n dsG
τ τΓ
= −∫ ( )1 , ,k
y y x xn n dsG
φ φΓ
= − −∫1
k
n
dsG
φ
φΓ
∂∂
= ∇∫ in
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
55
Se, per comodità di rappresentazione, si introducono la funzione: ( ) ( ),,z
x yx y
Gφ
φϑ
=
il sistema di equazionidiviene il seguente: ( )
0
1
20
52 1
2 2
insusu 1...
su ...
,
k
k k
k k
Nt
z t k kAkt
A
k N
ds A k NnM J dA AGJ
φφφ α
φ
ϑ φ α
Γ
=
⎧∆ = −⎪ = Γ⎪⎪ = Γ =⎪⎨ ∂
= Γ =⎪ ∂⎪⎪ = +⎪⎩
∫
∑∫
e le costanti: kk
zGααϑ
=
che, nel caso di sezione semplicementeconnessa (N=0), si riduce ad un problemadi Dirichlet per l’equazione di Poisson: 0
20
insu
Aφφ⎧∆ = −⎨
= Γ⎩
2,tz t A
t
M J dAGJ
ϑ φ= = ∫a cui segue il calcolo della curvatura torsionale:
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
56
OSSERVAZIONI
1) La soluzione del problema (5) non dipende dal sistema di riferimento scelto.
Infatti, l’operatore di Laplace e la derivata direzionale
non dipendono dal tipo di riferimento utilizzato.
2 2
2 2x y∂ ∂
∆ = +∂ ∂ n
∂∂
2) Esistenza e unicità della soluzione del problema (5).
Il problema (5) può essere inteso come una generalizzazione di un problema di Poissoncon condizioni al contorno di Dirichlet: esso è sempre “ben posto”, e la soluzione esistesempre ed è unica.
Le N costanti αk e la curvatura torsionale ϑz sono determinate mediante le N equazionidi circuitazione e la relazione che esprime Mt.
Inoltre, la funzione , le costanti e la costante dipendono solo dallageometria del dominio.
( ),x yφ kα tJ
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
57
3) Legame tra funzione potenziale degli sforzi e funzione di ingobbamento.
Il confronto tra le espressioni degli sforzi nei due approcci permette di desumerela relazione esistente tra le funzioni φ e Ψ:
La simmetria del dominio e delle condizioni al contornosuggeriscono l’utilizzo delle coordinate polarie l’indipendenza dall’anomalia θ : ( ) ( ),r rφ φ θ φ= =
λ = 1/100λ = 1/40 λ = 1/20 λ = 1/10 λ = 1/5 Se si utilizzano le quantitàcon indice “0” si compiono,al variare della spessore,i modesti errori relativiriportati in tabella.Si definisce (ad es. per Jt):
( )0
0 t tt
t
J Jerr JJ−
0+
0−
∼0
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
77
Torsione di profili aperti di piccolo spessore
Q
y
x
Mt
zL’analisi del tubo aperto di piccolo spessore ha mostrato che le quantità contrassegnate con l’indice “0” forniscono una stimasufficientemente accurata della soluzione.
Questa osservazione suggerisce l’utilizzo della funzione comeapprossimazione del potenziale degli sforzi, definito come:
0φ
Al fine di garantire un’adeguata accuratezza dei risultati, è necessario che lo spessore t(s) siaregolare e debolmente variabile ( |t'(s)|<<1 e |t''(s)|<<1/|t(s)|): in prossimità delle estremità dellasezione, di zone caratterizzate da brusche variazioni di spessore o di forti curvature della lineamedia, la soluzione approssimata basata su richiede, in generale, significative correzioni.
( )( )
0
22,
4t sn s nφ = −
t(s)
n
0φ
0φ
ns
t(s)
n = coordinata locale in direzione normalealla linea media del profilo
s = ascissa curvilinea lungo la linea media del profilo
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
78
( )( )
0
220,0 0
4tn nφ = − ≠
L’espressione utilizzata soddisfa la condizione di annullamentosul contorno, salvo che nelle estremità della sezione.
n
0φ
s
Si è però visto, negli esempi, che questa violazione dellacondizione al contorno porta a errori che rimangono limitatialle estremità della sezione (ovvero, la soluzione approssimatanon coglie gli “effetti di estremità”), ma che comunquenon influiscono sull’accuratezza della soluzione nel resto della sezione.
Inoltre, se lo spessore rispetta le condizioni esposte in precedenza e la linea media non ècaratterizzata da forti curvature si ha:
( ) ( ) ( )
( ) ( ) ( ) ( )
0
220
2
, 1 02
, 1 1 02 2
n s t s t ssn s t s t s t s
s
φ
φ
⎧∂ ′= ≅⎪⎪ ∂⎨∂⎪ ′ ′′⎡ ⎤= + ≅⎣ ⎦⎪ ∂⎩
( )2 2
0 00 2 2, 0 2 2n s
s nφ φ
φ∂ ∂
⇒ ∆ ≅ + ≅ − = −∂ ∂
Dunque, nell’ambito delle ipotesi fatte, soddisfa il problema di Poisson (5).( )( )
0
22,
4t sn s nφ = −
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
79
Si può allora procedere al calcolo della rigidezza geometrica torsionale (Jt):
( )02
22 ,
m
t
tn s dn dsφ
+
Γ −≅ ∫ ∫0
02t t AJ J dAφ= = ∫
( )( )
( ) 22 22
24m
t s
t s
t s n dn ds+
Γ −
⎡ ⎤= −⎢ ⎥⎣ ⎦∫ ∫ ( )31
3 mt s ds
Γ= ∫
( )313 m
tJ t s dsΓ
⇒ = ∫( Γm rappresenta la linea
media del profilo)
e alla determinazione della distribuzione di sforzi:
( )( )
0
22,
4t sn s nφ = − ⇒
t(s)n
s
τzsτzn
τzn τzs
( ) ( )0
0
1,2
, 2
zn s tz
zs tn
t s t sMGJ n
τ φϑ
τ φ
⎡ ⎤⎡ ⎤ ′+⎡ ⎤ ⎢ ⎥= = =⎢ ⎥⎢ ⎥ ⎢ ⎥−⎣ ⎦ ⎣ ⎦ ⎣ ⎦
τ
( ) ( )( ) 2
,2 1
t
t
t s t sMs n t sJ
⎛ ⎞ ′⎡ ⎤⇒ =± =⎜ ⎟ ⎢ ⎥±⎣ ⎦⎝ ⎠
τ
τzs
τznτ
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
80
Determinazione della funzione di ingobbamento .( ),n sΨ
Si osserva che il campo degli spostamenti assiali non varia in modo significativo nello
spessore del profilo (“piccolo” se confrontato con le dimensioni caratteristiche della sezione)
rispetto a come può variare lungo il profilo.
Dunque, si determina l’ingobbamento solo lungo la linea media del profilo: ( ) ( )0 00,s sΨ Ψ = Ψ
0 00
02 0
m
nn
nn nφ φ
=Γ =
∂ ∂= = − ≡
∂ ∂La simmetria del potenziale nello spessore fornisce:( )0 ,n sφ
00 0, , , ,x x y y x x y yn n n n
nφ
φ φ φ φ∂
= + ≅ +∂
Utilizzando le relazioni che legano e si ricava:φ Ψ
( ) ( ), ,y x x yx n y n= −Ψ − + Ψ − =
( ) ( ), ,x y y x x yn n xn yn= Ψ −Ψ − + ( ) ( ), ,x x y y x yt t xn yn= − Ψ + Ψ − + ( )P Qs
∂Ψ= − − −
∂in
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
81
ds( )
( )
0
0 0 00
mn n n
h s
P Qn n sφ φ
Γ= = =
∂ ∂ ∂Ψ≅ = − − − ≡
∂ ∂ ∂ni
( )( )0d s h s
dsΨ
⇒ = − ( ) ( )( )
2
0
2
0
ds
s
s h s ds kω
ω
⇒ Ψ = − +∫
= area settoriale
Q
y
xz
s
P
( )sω
( ) ( ) ( )0 2 0s sω ω⎡ ⎤Ψ = −⎣ ⎦
n
Q
y
xz
sPt
h(s)
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
83
0 20 40 60 80 100
0
20
40
60
80
100 Esempio: profilo a L [2/3]
0 20 40 60 80 1000
20
40
60
80
100
0
10
20
30
-2 0 2 4 6 8 10 12 14 16
0
2
4
6
8
10
12 concentrazione di sforzi(|τ|→∞ in angolo per soluzione analitica) Effetti di
estremità
( ),x yφ
τ
( ),x yφ
[mm]
[mm]
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
84
0 20 40 60 80 1000
10
20
30
40
50
60
70
80
90
100
Esempio: profilo a L [3/3]
0
50
100
0
50
100
-40-20
0
2040
( ),x y∗Ψ
( ),x y∗Ψ
N.B.: xC* = yC* ≅ 5.3 mm
( )00 0s∗Ψ ≡
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
85
Torsione di solidi con sezioni molteplicemente connessee pareti di piccolo spessore (“profili pluri-cellulari”)
Per piccoli spessori, gli N fori divengno N celle.
Così come fatto per i profili aperti, il potenziale degli sforziviene approssimato da un’espressione quadratica in n.Se si considera la parete comune alle celle k-esima e j-esimasi ha:
( ) ( ) ( ) ( )0, , , ,j kn s n s n s n sφ φ φ φ≅ + +
φ αk
αj
k j
k j
φk
φj
∼φ0( ) ( ) ( )
( ) ( )( )
( ) ( )( )
20 2
0
,,
4
, 1,2
, 1,2
z
jj j
z
kk k
z
n s t sn s n
G
n s nn sG t s
n s nn sG t s
φφ
ϑ
φφ α
ϑ
φφ α
ϑ
⎧⎪ = = −⎪⎪
⎡ ⎤⎪ = = −⎢ ⎥⎨⎢ ⎥⎪ ⎣ ⎦
⎪ ⎡ ⎤⎪ = = +⎢ ⎥⎪ ⎢ ⎥⎣ ⎦⎩
( ), 2n sφ⇒ ∆ ≅ −
nt
αk
αj
n
φ
s
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant
86
φ
k j ≅ ++
φ0
k j
φk
k j
φj
k j
Equazioni di circuitazione
nt
kΓ
k
kA
L’interpretazione fornita per le equazioni di congruenza di circuitazione(oppure, in alternativa, il PLV scritto con riferimento a un sistemaequilibrato scelto in modo opportuno) garantisce la validità delle stesseanche con riferimento a percorsi che, pur non coincidendonecessariamente coi bordi dei fori, contengono gli stessi.
2k
kds Anφ
Γ
∂=
∂∫
( ) ( ) ( ) ( )0, , , ,j kn s n s n s n sφ φ φ φ≅ + +
Corso di Scienza delle Costruzioni II - A.A. 2007-2008 - Appunti sul problema di De Saint-Venant