Re: 3D modelisation - seismic aggression

Hoi Magic

did you extract the displacement from your calculations, something along the following line (it is not in your comm file):
(replace MODELE etc by your own defintion):

MACRO_MATR_ASSE(MODELE=modmod,
CHAM_MATER=material,
CARA_ELEM=shellch,
CHARGE=Loaddyn,
NUME_DDL=CO('NUMEDDL'),
MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
OPTION='RIGI_MECA',),
_F(MATRICE=CO('MASSE'),
OPTION='MASS_MECA',),),);

modes=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
MATR_B=MASSE,
CALC_FREQ=_F(OPTION='BANDE',
FREQ=(1.0,50.0,),),);



result=CALC_ELEM(reuse =modes,
MODELE=Cmod,
CHAM_MATER=Amat,
RESULTAT=result,
OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
EXCIT=(_F(CHARGE=bcforce,),
_F(CHARGE=cylco2,),),);

result=CALC_NO(reuse =result,
RESULTAT=result,
OPTION=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',),);


IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=Cmesh,
RESULTAT=result,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','DEPL',),),);

(this is for use version CA10 and below, CA11 changed some of the NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','DEPL',),),) parameters.


Edit sunday 2012-02-05:
I didnot download your last file, where you do obtain your results.
So sorry for the mistake.

Hoi Magic

magic3306 wrote: Volume basis : 8 x 8 x 9 - 6.4 x 6.4 x 9 = 207.36 m3
Door : 1.2 x 2.7 x .8 = 2.592
Window x 2 : 1.7 x 2.2 x .8 x 2 = 5.984
Ventilation : .65 x 1.2 x 0.8 = 0.624
Roof : 2 x (7 x .25 x .3) + 2 x (6.4 x .25 x .3) = 2.01

VOLUME = Volume basis - Door - Window x 2 - Ventilation - Roof = 196.15 m3
Mass : 2000 x VOLUME = 392300 kg

Oh... I found only 104 kg in more than the automated calculation!!!
So there is an error of 0.026 % !! it's nothing!!!

I guess 0.026 % is too much.
I recalculated your volume and got 196.098 m3, and in the order of e-10 mismatch for the mass in [kg].

Try the following matlab/octave program and you can see that the volume of the ventilation is not the same:

dx1 = abs(5.05-6.75)
dy1 = abs(3.3-1.1)
A1 = dx1*dy1
V1 = A1*0.8

dx2 = abs(5.6-6.8)
dy2 = abs(2.7-0.0)
A2 = dx2*dy2
V2 = A2*0.8

dx3 = abs(1.9-3.2) %% different
dy3 = abs(1.1-0.45)
A3 = dx3*dy3
V3= A3*0.8

dx4 = abs(6.75-5.05)
dy4 = abs(1.1-3.3)
A4 = dx4*dy4
V4 = A4*0.8

V51 = 2*7.0*0.3*0.25
V52 = 2*6.4*0.3*0.25
V5 = V51+V52

Aopen = dx1*dy1+dx2*dy2+dx3*dy3+dx4*dy4
Vopen = Aopen*0.8

Vwand2 = 8*8*9
Vwand1 = (8-2*0.8)*(8-2*0.8)*9
Vwand = Vwand2-Vwand1

V1
V2
V3
V4
V5

Vtot = Vwand- V1-V2-V3-V4-V5

mtot = Vtot*2000
dm = mtot-3.92196E+05 %% result: dm = -1.1642e-10