|
|
[ÇóÖú]
"С³æÇóÖú¡° 䶨Ò庯Êý»ò±äÁ¿ 'Beam_InputData547'¡£ ÒÑÓÐ1È˲ÎÓë
С³æ¸ÕѧϰMATLAB£¬ÏÖÕýÔÚѧϰ¡¶MatlabÓÐÏÞÔª½á¹¹¶¯Á¦Ñ§·ÖÎöÓ빤³ÌÓ¦Óá·£¬ÔÚ5.4.3˲̬ÎÊÌâ·ÖÎö¡£Àý5.7Öа´ÕÕÊé±¾¸½´øµÄÔ´³ÌÐòÔËÐÐÌáʾ
䶨Ò庯Êý»ò±äÁ¿ 'Beam_InputData547'¡£
³ö´í EX547 (line 28)
Beam_InputData547; % import the input data for the information of
´úÂ뿪ͷ²¿·ÖÈçÏ£º
%-------------------------------------------------------------------------%
% Example 5.4.7
% To solve transient response of a 2-d frame structure.
% The solution methods are: 1) central difference scheme. 3) Houbolt integration scheme.
% 4) Wilson ?? integration scheme. 5) Newmark integration scheme
% nodal dof: {u1 v1 w1 x1 y1 z1 u2 v2 w2 x2 y2 z2}
% Problem description
% Find the response of a frame structure which is made of three beams of lengths of 4 m,
% 3 m and 4 m, respectively. All beams have cross- section of 0.10 m height by 0.05 m
% width. The elastic modulus is 2.10x10^11 Pa. The frame is subjected to an impulse load
% of amplitude 500 N in the middle of the top beam. One end of the each vertical beam is
% fixed. (see Fig. 5-9 for the element discretization)
% Variable descriptions
% k, m - element stiffness matrix and mass matrix
% kk, mm - system stiffness matrix and mass matrix
% ff - system force vector
% index - a vector containing system dofs associated with each element
% bcdof - a vector containing dofs associated with boundary conditions
% bcval - a vector containing boundary condition values associated with the dofs in 'bcdof'
% dsp - displacement matrix
% vel - velocity matrix
% acc - acceleartion matrix
%--------------------------------------------------------------------------
% (0) input control data
%--------------------------------------------------------------------------
clear; clc;
Beam_InputData547; % import the input data for the information of
% nodes, elements, loads, constraints and materials
Opt_beam=1; % option for type of the beam
% =1 Euler Bernoulli beam
% =2 Timoshenko beam
Opt_mass=2; % option for mass matrix
% =1 consistent mass matrix
% =2 lumped mass matrix
Opt_section=1; % option for type of cross-section
% = 1 rectangular cross-section
% = 2 circular cross-section
TypeMethod=1; % option for selecting the solution method
% = 1 central difference scheme
% = 3 Houbolt integration scheme
% = 4 Wilson integration scheme
% = 5 Newmark integration scheme
Typeload=1; % option for selecting the load type
% = 1 impulse load
% = 2 step load
% = 3 Harmonic load
dt=0.00001; % time step size
ti=0; % initial time
tf=0.200; % final time
nt=fix((tf-ti)/dt); % number of time steps
tt=ti:dt:ti+nt*dt; % generate time samples vector
ac=0.00002; bc=0.00008; % Parameters for proportional damping
al=0; % angle between the reference coordinate system and
% the local coordinate system for the space element
%--------------------------------------------------------------------------
% (1) initialization of matrices and vectors to zero
%--------------------------------------------------------------------------
k=zeros(No_nel*No_dof,No_nel*No_dof); % element stiffness matrix
m=zeros(No_nel*No_dof,No_nel*No_dof); % element mass matrix
kk=zeros(Sys_dof,Sys_dof); % initialization of system stiffness matrix
mm=zeros(Sys_dof,Sys_dof); % initialization of system mass matrix
ff=zeros(Sys_dof,1); % initialization of system force vector
bcdof=zeros(Sys_dof,1); % initializing the vector bcdof
bcval=zeros(Sys_dof,1); % initializing the vector bcval
index=zeros(No_nel*No_dof,1); % initialization of index vector
------------------------------------------------------------------------------------------------------
䶨Òå±äÁ¿ÊÇ·ñ¿ÉÒÔÀí½âΪBeam_InputData547=£¿£¿£¿
»òBeam=input£¨data547£©£¬µ«ÔÚÔ´ÎļþÖв¢Î´ÕÒµ½Ïà¹ØÐÅÏ¢¡£
ÁíÍâ±¾ÕÂÔ´³ÌÐò¿ªÊ¼µÚÒ»¾ä¾ùΪ¡±Beam_InputData ¡±
¾ù³öÏÖÈçÉÏÎÊÌ⣬¿ÒÇë´óÉñ´Í½Ì£¡£¡
ÏÖ½«µÚÎåÕÂÏà¹ØÄÚÈݼ°³ÌÐòÔÚ¸½ÉÏ¡£¡£¡£!["С³æÇóÖú¡° 䶨Ò庯Êý»ò±äÁ¿ 'Beam_InputData547'¡£]()
´íÎóÌáʾ.png |
» ±¾Ìû¸½¼þ×ÊÔ´Áбí
-
»¶Ó¼à¶½ºÍ·´À¡£ºÐ¡Ä¾³æ½öÌṩ½»Á÷ƽ̨£¬²»¶Ô¸ÃÄÚÈݸºÔð¡£
±¾ÄÚÈÝÓÉÓû§×ÔÖ÷·¢²¼£¬Èç¹ûÆäÄÚÈÝÉæ¼°µ½ÖªÊ¶²úȨÎÊÌ⣬ÆäÔðÈÎÔÚÓÚÓû§±¾ÈË£¬Èç¶Ô°æÈ¨ÓÐÒìÒ飬ÇëÁªÏµÓÊÏ䣺xiaomuchong@tal.com
- ¸½¼þ 1 : µÚÎåÕÂ_Ô´³ÌÐò.rar
2015-08-31 15:07:57, 40.55 K
» ²ÂÄãϲ»¶
|