Znn3bq.jpeg
²é¿´: 4356  |  »Ø¸´: 3
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

tyutzxy

½ð³æ (СÓÐÃûÆø)

[ÇóÖú] "С³æÇóÖú¡° 䶨Ò庯Êý»ò±äÁ¿ '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

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ashton

гæ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by lovellbaby at 2018-08-02 16:45:09
Ô´³ÌÐòÀïûÓÐÌṩÏà¹ØµÄԭʼÊý¾Ý£¿

ûÓÐBeam_InputData-----һϵÁÐÊý¾Ý
4Â¥2021-06-11 12:19:40
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 4 ¸ö»Ø´ð

lovellbaby

¾èÖú¹ó±ö (ÕýʽдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

Ô´³ÌÐòÀïûÓÐÌṩÏà¹ØµÄԭʼÊý¾Ý£¿
dajiahaohaohao
2Â¥2018-08-02 16:45:09
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 282£¬Çóµ÷¼Á +8 jggshjkkm 2026-04-09 9/450 2026-04-10 02:10 by belle15
[¿¼ÑÐ] ¿¼Ñе÷¼Á +18 ˶ÐǸ° 2026-04-09 19/950 2026-04-10 02:08 by belle15
[¿¼ÑÐ] 314Çóµ÷¼Á +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[¿¼ÑÐ] 266Çóµ÷¼Á +25 ÑôÑôÍÛÈû 2026-04-07 25/1250 2026-04-09 20:32 by onlyÖÜ
[¿¼ÑÐ] µ÷¼Á +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[¿¼ÑÐ] һ־Ը³¶«´óѧ071000ÉúÎïѧѧ˶³õÊÔ·ÖÊý276Çóµ÷¼Á +19 Ľ¾øcc 2026-04-09 23/1150 2026-04-09 17:14 by ¶ÎΰÑÞ
[¿¼ÑÐ] Ò»Ö¾Ô¸ÎäÀí³µÁ¾ 281 Çóµ÷¼Á +5 Éϰ¶Ñо¿Éú. 2026-04-07 5/250 2026-04-09 15:56 by onlyÖÜ
[¿¼ÑÐ] 280Çóµ÷¼Á +7 wzzzÍõ 2026-04-09 7/350 2026-04-09 15:32 by ÊÍ·ÅÌìÐÔ
[¿¼ÑÐ] 296Çóµ÷¼Á +3 Íô£¡£¿£¡ 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[¿¼ÑÐ] 318Çóµ÷¼Á +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[¿¼ÑÐ] ±¾¿Æ211£¬293·ÖÇëÇóµ÷¼Á +12 Á«²Ë¾ÍÊÇź°É 2026-04-03 13/650 2026-04-08 20:30 by ±³¶Ô´óº£³ö·¢
[¿¼ÑÐ] Çóµ÷¼Á +11 wwwwabcde 2026-04-07 11/550 2026-04-07 23:16 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϺ½£¬ÊýһӢһѧ˶317Çóµ÷¼Á£¡£¡ +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by ¿¼ÑÐѧУÕеãÈË
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +6 Ò»ÑùYWY 2026-04-05 6/300 2026-04-05 20:30 by ÄϺ½~ÍòÀÏʦ
[¿¼ÑÐ] 272Çóµ÷¼Á +4 µçÆøÀî 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
[¿¼ÑÐ] ¿¼Ñе÷¼Á +11 СsunÒªºÃÔË 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
[¿¼ÑÐ] Ò»Ö¾Ô¸¶«±±´óѧ085901ÍÁľר˶345Çóµ÷¼Á +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 350Ò»Ö¾Ô¸±±¾©º½¿Õº½Ìì´óѧ08500²ÄÁÏ¿ÆÑ§Ó빤³ÌÇóµ÷¼Á +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 320Çóµ÷¼Á +3 ũҵ¹¤³ÌÓëÐÅÏ¢¼ 2026-04-03 3/150 2026-04-03 11:40 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +4 Ò»ÑùYWY 2026-04-03 4/200 2026-04-03 09:48 by À¶ÔÆË¼Óê
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û