*HEADING DSLOAD AND TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS [C3D6] AND MATRIX INPUT *NODE,NSET=ALL 1 , 0. , 0. , 0. 2 , 7. , 0. , 0. 3 , 7. , 0. , -7. 4 , 0. , 0. , -7. 5 , 0. , 7. , 0. 6 , 7. , 7. , 0. 7 , 7. , 7. , -7. 8 , 0. , 7. , -7. *NSET,NSET=NS1 1,4,5,8 *NSET,NSET=NS2 1,2,3,4 *NSET,NSET=NS3 1,2,5,6 *ELEMENT,TYPE=C3D6,ELSET=ELEMS 1, 4,3,8,1,2,5 ** The following element is modeled thorugh *MATRIX INPUT ** 2, 3,7,8,2,6,5 ** ** Surface elements added to apply loads on matrix nodes ** *ELEMENT, TYPE=SFM3D3, ELSET=ESURF1 3,3,8,7 *ELEMENT, TYPE=SFM3D3, ELSET=ESURF2 4,2,6,5 *ELEMENT, TYPE=SFM3D4, ELSET=ESURF3 5,2,3,7,6 *ELEMENT, TYPE=SFM3D4, ELSET=ESURF4 6,5,6,7,8 *ELSET,ELSET=ESURF ESURF1,ESURF2,ESURF3,ESURF4 *SURFACE SECTION,ELSET=ESURF,DENSITY=0.0 *surface, NAME=SIDE_1 ELEMS, S1 ESURF1, SPOS *surface, NAME=SIDE_2 ELEMS, S2 ESURF2, SPOS *surface, NAME=SIDE_3 1, S3 *surface, NAME=SIDE_4 ESURF3, SPOS *surface, NAME=SIDE_5 ESURF4, SPOS *surface, NAME=SIDE_6 1, S5 *SOLID SECTION,MATERIAL=ELMAT, ELSET=ELEMS *MATERIAL,NAME=ELMAT *ELASTIC 3.E6,0.3 *DENSITY 10.0, *EXPANSION .0001, ** Matrix data for element 2 *MATRIX INPUT,NAME=MATRIX1,TYPE=SYMMETRIC,INPUT=tempdsl_stiff.inp *MATRIX ASSEMBLE,STIFFNESS=MATRIX1 *MATRIX INPUT,NAME=MATRIX2,TYPE=SYMMETRIC,INPUT=tempdsl_mass.inp *MATRIX ASSEMBLE,MASS=MATRIX2 *INITIAL CONDITIONS,TYPE=TEMPERATURE ALL,-10. *BOUNDARY NS1,1 NS2,2 NS3,3 **************************** *STEP,PERTURBATION *STATIC *DSLOAD SIDE_1,P,3000. SIDE_3,P,3000. SIDE_5,P,3000. ***NODE PRINT ** RF, ***EL PRINT ** LOADS, *OUTPUT,FIELD *ELEMENT OUTPUT S, *NODE OUTPUT RF, U, *OUTPUT,HISTORY *ENERGY OUTPUT *END STEP **************************** *STEP,PERTURBATION *STATIC *DSLOAD SIDE_2,P,3000. SIDE_4,P,3000. SIDE_6,P,3000. *END STEP *************************** *STEP *Frequency 5, *End step *************************** *STEP,PERTURBATION *STATIC *DSLOAD SIDE_1,HP,3000.,0.,-7. SIDE_3,HP,3000.,0.,-7. SIDE_5,HP,3000.,0.,-7. *END STEP *************************** *STEP,PERTURBATION *STATIC *DSLOAD SIDE_2,HP,3000.,0.,-7. SIDE_4,HP,3000.,0.,-7. SIDE_6,HP,3000.,0.,-7. *END STEP *************************** *STEP,PERTURBATION *STATIC *TEMPERATURE ALL,40. ***EL PRINT ** S, ***NODE PRINT ** RF, *OUTPUT,FIELD *ELEMENT OUTPUT S, *NODE OUTPUT RF, U, *OUTPUT,HISTORY *ENERGY OUTPUT *END STEP