function in_data = beam_3d_a(in_data)
%unités mm N MPa
%définition des noeuds

in_data.ND = [   
1	  0	    1	1 ;
2	  10	1	1 ;
3	  20	1	1 ;
4	  30	1	1 ;
5	  40	1	1 ;
6	  50	1	1 ;
7	  60	1	1 ;
8	  70	1	1 ;
9	  80	1	1 ;
10	  90	1	1 ;
11	  100	1	1 ;
12	  110	1	1 ;
13	  120	1	1 ;
14	  130	1	1 ;
15	  140	1	1 ;
16	  150	1	1 ;
17	  160	1	1 ;
18	  170	1	1 ;
19	  180	1	1 ;
20	  190	1	1 ;
21	  200	1	1 ;
22	  210	1	1 ;
23	  220	1	1 ;
24	  230	1	1 ;
25	  240	1	1 ;
26	  250	1	1 ;
27	  260	1	1 ;
28	  270	1	1 ;
29	  280	1	1 ;
30	  290	1	1 ;
31	  300	1	1 ;
32	  310	1	1 ;
33	  320	1	1 ;
34	  330	1	1 ;
35	  340	1	1 ;
36	  350	1	1 ;
37	  360	1	1 ;
38	  370	1	1 ;
39	  380	1	1 ;
40	  390	1	1 ;
41	  400	1	1 ];

%paramètres matériaux et géométrie poutre
% AISI_Steel_1005

E = 200000;
nu=0.25;
G = E/(2*(1+nu));
b = 30;
ho = 30;
rho = 25000;

%E   = 8e9;
%G   = 3e9;
%b   = 0.3;
%ho  = 0.4;
%rho = 1000;
%nu  = 0.33;


%initialise les carac du matériau à 0 : inutile ??

in_data.mater.rho=0; in_data.mater.b=0; in_data.mater.ho=0;
in_data.mater.E=0; in_data.mater.A=0; in_data.mater.I=0; in_data.mater.G=0;

in_data.nu = nu;
%définitions des barres

in_data.EL = [  
1	3	1	2	E	G	b	ho	rho	nu;
2	3	2	3	E	G	b	ho	rho	nu;
3	3	3	4	E	G	b	ho	rho	nu;
4	3	4	5	E	G	b	ho	rho	nu;
5	3	5	6	E	G	b	ho	rho	nu;
6	3	6	7	E	G	b	ho	rho	nu;
7	3	7	8	E	G	b	ho	rho	nu;
8	3	8	9	E	G	b	ho	rho	nu;
9	3	9	10	E	G	b	ho	rho	nu;
10	3	10	11	E	G	b	ho	rho	nu;
11	3	11	12	E	G	b	ho	rho	nu;
12	3	12	13	E	G	b	ho	rho	nu;
13	3	13	14	E	G	b	ho	rho	nu;
14	3	14	15	E	G	b	ho	rho	nu;
15	3	15	16	E	G	b	ho	rho	nu;
16	3	16	17	E	G	b	ho	rho	nu;
17	3	17	18	E	G	b	ho	rho	nu;
18	3	18	19	E	G	b	ho	rho	nu;
19	3	19	20	E	G	b	ho	rho	nu;
20	3	20	21	E	G	b	ho	rho	nu;
21	3	21	22	E	G	b	ho	rho	nu;
22	3	22	23	E	G	b	ho	rho	nu;
23	3	23	24	E	G	b	ho	rho	nu;
24	3	24	25	E	G	b	ho	rho	nu;
25	3	25	26	E	G	b	ho	rho	nu;
26	3	26	27	E	G	b	ho	rho	nu;
27	3	27	28	E	G	b	ho	rho	nu;
28	3	28	29	E	G	b	ho	rho	nu;
29	3	29	30	E	G	b	ho	rho	nu;
30	3	30	31	E	G	b	ho	rho	nu;
31	3	31	32	E	G	b	ho	rho	nu;
32	3	32	33	E	G	b	ho	rho	nu;
33	3	33	34	E	G	b	ho	rho	nu;
34	3	34	35	E	G	b	ho	rho	nu;
35	3	35	36	E	G	b	ho	rho	nu;
36	3	36	37	E	G	b	ho	rho	nu;
37	3	37	38	E	G	b	ho	rho	nu;
38	3	38	39	E	G	b	ho	rho	nu;
39	3	39	40	E	G	b	ho	rho	nu;
40	3	40	41	E	G	b	ho	rho	nu];

%définition contraintes

in_data.CON = [1  0  0  0  0  0  0];

%définition chargement

in_data.LOAD_ = [  
41   0   1000  0   0   0   0];

%paramètres masse et dynamiques

in_data.MASS = [ ];

in_data.dynam.TIMEH     = [ 'bedr.txt' ];
in_data.dynam.delta_tm  = [0.0079];
in_data.dynam.TIMEHDIR  = [1 0 0 0 0 0];
in_data.dynam.TIMEHM    = [1:1:length(in_data.ND)*6];
in_data.dynam.TIMEHPL   = [19*6-2];
in_data.dynam.DAMP_C    = [0.06 0.06];
in_data.dynam.DAMP_F    = [3];        
in_data.dynam.ab        = [0.04 0.001]; 
in_data.dynam.MODA      = [1];          
