# hrbeu221's blog

## rotation spring setting

Submitted by hrbeu221 on Thu, 2012-07-26 13:35.Dear all, these days I use combin14 to simulate the boundary conditions. I use two types of springs to simulate support. One is linear displacement spring, the other one is rotational spring. When the spring is on the radial edges, the rotation is about radius direction. When the spring is on the circumferential direction, the rotation is about circumferential direction. So i translate all the nodes to local cylndrial coordinate system. When the spring constants are large enough, the simulation results do not meet well with the clamped boundary conditions. I am not sure for the rotation spring in cylindrical coordinate system. But i have already tried on rectangular plate, it is wright. So please help me on this problem. Thank you for your concern

## Vibration analysis with uncertainty properties

Submitted by hrbeu221 on Wed, 2012-07-11 20:56.Dear all,

These days I am doing some research on vibration analysis with uncertain properties. In traditional, Monte Carlo Simulation is used to do this simulation. But the drawback is time-consuming for the determinstic analysis.

For this, how to solve this problem? Please help me on it. Please give me some suggestion on it.

Thank you for your concern.

- 1 comment
- 677 reads

## compiler for FORTRAN

Submitted by hrbeu221 on Mon, 2012-05-21 16:21.Hello, everyone,

These days, I am studying the Fortran. I want to know which complier is good for fortran. My computer system is 64 bit Win7 system. I want to use fortran to do some numerical calculation about the structure dynamics.

Thank you for your sugestion.

- 2 comments
- 1584 reads

## orthotropic plates analysis

Submitted by hrbeu221 on Sun, 2012-04-29 19:44.
Hi, please help me on the simulation on the orthotropic material in ABAQUS. I have the parameters for the material density=1584;

Ex=127.9e9

Ey=10.27e9

muxy=0.22;

muyx=0.0177;

Gxy=7.312e9;How can i use the orthotropic or engineering constant to simulate the material. I want to calculate frequency for orthotropic plates.

Graphite epoxy plates

Graphite/epoxy consisting of 12 plies oriented at 0 and fabricated from Fiberite

Hy-E1034C prepreg

Thank you for your concern.

## Explicit Integral

Submitted by hrbeu221 on Tue, 2012-03-06 15:51.When use matlab or mathematica integrate this y = int(cos(m*pi*x/a)*cos(n*pi*x/a)/(x+a0)^2,x,0,a); (m,n=0,1,2,3,4,5...)

There is some explicit integral. I try to solve this problem. But the result is complex.

Is there some effective way to solve this problem.

Thanks!

## Vibration of Circular Plates

Submitted by hrbeu221 on Thu, 2012-02-16 23:40.Dear All,

These days I am do some research on vibration of circular plates (out-of-plane and in-plane). In the reference, the displacement in radial direction can be expressed with Bessel polynomial. I want to use the Rayleigh-Ritz method to solve these problems. The strain energy and kinetic energy can obtain with integral. If there is the bessel function, how to calculate the integral. Also I want to consider some spring in the boundary to change the boundary conditions.

Or there are some other method to do? Express the displacement in the radial direction with trigonometric function?

Hope all can give me some advice on it.

Thanks.

- 1 comment
- 1773 reads

## Circle transform to square

Submitted by hrbeu221 on Fri, 2012-01-20 11:32.
**Hi, everyone, Is there some methd for the coordinate transform from circle to square.
**

Thank you for your suggestion.

- 2 comments
- 1661 reads

## Ask for help for one paper

Submitted by hrbeu221 on Fri, 2011-10-28 21:18.Ask help for one paper in the inter noise conference 2011

Title: Application of virtual mode synthesis method for ship compartment impact noise

Page: 147-156

Thank you very much.

## draw mode shape with matlab

Submitted by hrbeu221 on Tue, 2011-09-13 15:04.
if M is your mass matrix and K is the stiffness matrix, use

[phi w2]=eig(K,M);

phi contains the eigenvectors (each column is a modeshape)

w2 is a diagonal matrix, each diagonal entry contains wn(i)^2

(i.e. the i'th natural frequency squared).

Then how to draw the mode shape with eigenvectors?

Thank you.

- 3 comments
- 6970 reads

## Recent comments

19 hours 18 min ago

1 day 2 hours ago

1 day 7 hours ago

1 day 21 hours ago

1 day 23 hours ago

1 day 23 hours ago

2 days 2 hours ago

2 days 4 hours ago

2 days 5 hours ago

3 days 7 hours ago