文件列表:

S_e.m (407, 2018-05-10)
ancf_beam.m (1212, 2018-05-10)
ancf_params.m (635, 2018-05-10)
animate.m (764, 2018-05-10)
calc_points.m (879, 2018-05-10)
constraints.m (1139, 2018-05-10)
elastic_forces.m (2230, 2018-05-10)
gauss_points.m (1073, 2018-05-10)
gravitational_forces.m (799, 2018-05-10)
init_cond.m (331, 2018-05-10)
integrators (0, 2018-05-10)
integrators\DAE_newmark.m (3195, 2018-05-10)
integrators\ODE_euler.m (1324, 2018-05-10)
integrators\ODE_rk2.m (1816, 2018-05-10)
mass_matrix.m (1185, 2018-05-10)
plot_length.m (925, 2018-05-10)
plot_points.m (1888, 2018-05-10)
shape_fun.m (883, 2018-05-10)
snapshots.m (1361, 2018-05-10)
sym (0, 2018-05-10)
sym\Generate_eps_eps_e.m (729, 2018-05-10)
sym\Generate_grav.m (411, 2018-05-10)
sym\Generate_kappa_kappa_e.m (1337, 2018-05-10)
sym\Generate_mass.m (463, 2018-05-10)
sym\eps_eps_e.m (936, 2018-05-10)
sym\grav.m (378, 2018-05-10)
sym\kappa_kappa_e.m (1529, 2018-05-10)
sym\mass.m (1039, 2018-05-10)
sym\tilde.m (149, 2018-05-10)

# ANCF_beam Matlab simulation of a flexible beam modeled with gradient-deficient ANCF beam elements. ## 1. Set problem parameters by modifying the file `ancf_params.m` These include: - total length (L) - radius (r) - number of ANCF cable elements (ne) as well as material properties: - elasticity modulus (p.E) - material density (p.rho) and the end constraints (0: free, 1: ball joint, 2: weld joint) - left side (p.leftCnstr) - rightSide (p.rightCnstr) The default parameter values correspond to a cable of length 3 and radius 0.02, clamped at both ends and modeled using 5 ANCF cable elements. The elasticity modulus is 2e7 and density 7200. ## 2. Call the `ancf_beam` function Specify - duration (TEND) - step-size (H) - integration method (METHOD) ```matlab ancf_beam Simulate a flexible noodle modeled with ANCF beam elements [data,params] = ancf_beam(TEND, H, METHOD) simulates a noodle over the time interval [0,TEND] using a step-size H and the specified method: 'newmark' - use a DAE formulation and Newmark integration 'rk2' - use an ODE formulation and an RK2 method (change value 'a' in ODE_rk2 to select a different RK2 method; default is midpoint) 'euler' - use an ODE formulation and forward Euler integration ``` ------- Sample results ``` >> [data, params] = ancf_beam(1, 1e-3, 'newmark') data = struct with fields: t: [1×1001 double] e: [36×1001 double] ed: [36×1001 double] edd: [36×1001 double] lam: [12×1001 double] params = struct with fields: ne: 5 n: 36 L: 0.6000 E: 20000000 A: 0.0013 I: 1.2566e-07 rho: 7200 g: [3×1 double] leftCnstr: 2 rightCnstr: 2 ```

近期下载者：

相关文件：

评论：[我要评论] [举报此文件]

收藏者：