Giới thiệu dự án

Nghiên cứu ứng dụng "Các yếu tố ảnh hưởng đến động lực làm việc của người lao động tại Công ty Cổ phần Phương Minh" giải quyết bài toán cốt lõi trong quản trị nguồn nhân lực (Human Resource Management - HRM) tại các doanh nghiệp sản xuất và thi công xây lắp kỹ thuật. Trong bối cảnh kinh tế tri thức và công nghiệp hóa, nguồn nhân lực kỹ thuật cao quyết định trực tiếp đến 85% năng lực cạnh tranh và hiệu quả thi công công trình. Tuy nhiên, tình trạng suy giảm động lực, biến động lao động kỹ thuật và hiệu suất không đồng đều đang là thách thức lớn đối với doanh nghiệp.

+-------------------------------------------------------------------------------+
|                       MÔ HÌNH HỆ THỐNG ĐỘNG LỰC LAO ĐỘNG                      |
|                                                                               |
|  [Đặc điểm công việc (DDCV)]      ---\                                        |
|  [Điều kiện làm việc (DKLV)]       ----\                                      |
|  [Lương, thưởng & phúc lợi (LTPL)] -----> [ĐỘNG LỰC LÀM VIỆC CỦA NLĐ (DLLV)]  |
|  [Quan hệ đồng nghiệp/sếp (QH)]    ----/                                      |
|  [Đào tạo & thăng tiến (DTTT)]    ---/                                        |
+-------------------------------------------------------------------------------+

Bối cảnh và vấn đề nghiên cứu

Công ty Cổ phần Phương Minh (FUMICO) hoạt động chính trong lĩnh vực sản xuất cọc bê tông ly tâm dự ứng lực, đúc kim loại và thi công xây lắp điện đến 220KV tại Thừa Thiên Huế. Với quy mô 383 lao động (năm 2019), trong đó lao động trực tiếp chiếm 76,11% và công nhân kỹ thuật chiếm 83,60%, công ty đối mặt với các điểm nghẽn nghiêm trọng:

  • Áp lực nợ ngắn hạn chiếm tỷ trọng lớn trên tổng nguồn vốn (72,47 tỷ VNĐ trên 92,88 tỷ VNĐ tổng tài sản năm 2019, tương đương 78,02%), gây sức ép lên chính sách dòng tiền chi trả phúc lợi.
  • Tỷ lệ hao mòn tài sản cố định và điều kiện thi công nặng nhọc ngoài hiện trường tác động trực tiếp đến tâm lý và thể lực người lao động (lao động nam chiếm 84,45%).
  • Thiếu hụt công cụ định lượng chuẩn hóa để đo lường chính xác mức độ thỏa mãn và động lực nội tại của công nhân viên.

Mục tiêu nghiên cứu

  1. Hệ thống hóa cơ sở lý luận về động lực lao động dựa trên Tháp nhu cầu Maslow (1943), Thuyết kỳ vọng Vroom (1964) và Thuyết công bằng Adams (1963).
  2. Xây dựng và kiểm định mô hình định lượng 5 nhóm nhân tố với 20 biến quan sát tác động đến Động lực làm việc ($DLLV$).
  3. Thu thập dữ liệu thực nghiệm ($n = 150$), xử lý qua bộ công cụ trích xuất thống kê chuyên sâu.
  4. Đề xuất ma trận giải pháp điều hành nhân sự thực thi theo từng nhóm biến có hệ số tác động chuẩn hóa cao nhất.

Phạm vi và giới hạn

  • Không gian: Toàn bộ các phân xưởng sản xuất (Bê tông, Đúc kim loại, Cơ khí) và khối phòng ban chức năng tại Công ty Cổ phần Phương Minh.
  • Thời gian: Phân tích dữ liệu thứ cấp giai đoạn 2017–2019; khảo sát dữ liệu sơ cấp từ tháng 10/2020 đến tháng 01/2021.
  • Giới hạn: Mẫu nghiên cứu lấy theo phương pháp chọn mẫu thuận tiện; giới hạn trong mô hình hồi quy tuyến tính cổ điển OLS.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp quản trị động lực truyền thống tại doanh nghiệp chủ yếu mang tính định tính, phụ thuộc vào cảm tính của cán bộ quản lý phân xưởng, dẫn đến sự bất bình đẳng trong phân phối thu nhập và đánh giá hiệu suất.

Tiêu chí so sánh Phương pháp quản trị truyền thống Mô hình định lượng ứng dụng (Đề tài)
Cơ sở đánh giá Báo cáo cảm tính, quan sát hành vi rời rạc Thang đo Likert 5 điểm chuẩn hóa, kiểm định EFA/OLS
Độ tin cậy dữ liệu Không có chỉ số đo lường sai số Cronbach's Alpha > 0.7; KMO > 0.5; Sig. < 0.05
Tính cá nhân hóa Áp dụng cào bằng toàn bộ công ty Phân tầng rõ rệt giữa khối sản xuất và khối văn phòng
Chi phí vận hành Thấp nhưng rủi ro thất thoát nhân sự cao Tối ưu hóa chi phí chính sách nhờ trúng điểm nghẽn

Ma trận ưu tiên yêu cầu nghiên cứu (MoSCoW):

  • Must have: Kiểm định độ tin cậy thang đo Cronbach's Alpha; kiểm định độ hội tụ và phân biệt EFA; kiểm định đa cộng tuyến (VIF < 2.0).
  • Should have: Phân tích tương quan Pearson; kiểm định hồi quy đa biến OLS; phân tích phương sai ANOVA.
  • Could have: Kiểm định trung bình mẫu One-Sample T-test so sánh với ngưỡng kỳ vọng $\mu = 3.0$.
  • Won't have: Mô hình cấu trúc tuyến tính SEM bậc cao (giữ nguyên phạm vi xử lý SPSS để tối ưu tính thực thi cho doanh nghiệp).

Thiết kế hệ thống và mô hình nghiên cứu

Hệ thống biến quan sát được mã hóa thành 5 nhóm biến độc lập và 1 biến phụ thuộc:

$$\text{DLLV} = \beta_0 + \beta_1 \cdot \text{DKLV} + \beta_2 \cdot \text{DDCV} + \beta_3 \cdot \text{LTPL} + \beta_4 \cdot \text{QH} + \beta_5 \cdot \text{DTTT} + \varepsilon$$

   [BẢNG MÃ HÓA THANG ĐO CHI TIẾT (20 BIẾN ĐỘC LẬP + 3 BIẾN PHỤ THUỘC)]

+-----------------------------------------------------------------------------------+
|               KIẾN TRÚC LUỒNG XỬ LÝ DỮ LIỆU ĐỊNH LƯỢNG (DATA PIPELINE)            |
|                                                                                   |
|  [Phiếu khảo sát Likert 5] -> [Mã hóa Excel 2019] -> [SPSS v20.0 / Python Engine]  |
|                                                              |                    |
|       +------------------------------------------------------+                    |
|       |                                                                           |
|       v                                                                           |
|  [Cronbach's Alpha] -> [EFA Factor Extraction] -> [Pearson Correlation]           |
|                                                          |                        |
|                                                          v                        |
|  [Ma trận giải pháp quản trị] <--- [One-Sample T-test] <--- [OLS Multiple Regress]|
+-----------------------------------------------------------------------------------+

Phương pháp nghiên cứu (Methodology)

Quy trình nghiên cứu kết hợp phương pháp nghiên cứu định tính (phỏng vấn sâu chuyên gia, hiệu chỉnh thang đo) và nghiên cứu định lượng:

  1. Xác định cỡ mẫu:
    • Theo công thức Hair et al. (1998): $n \ge 5 \times m = 5 \times 20 = 100$ quan sát.
    • Theo Tabachnick & Fidell (2001): $n \ge 50 + 8 \times k = 50 + 8 \times 5 = 90$ quan sát.
    • Mẫu thực tế phát ra và thu về hợp lệ: $n = 150$ (đạt 150% ngưỡng tối thiểu, đảm bảo độ tin cậy thống kê $\alpha = 0.05$).
  2. Ngăn ngừa rủi ro: Loại bỏ phiếu lỗi bằng kiểm tra logic tự động trên dải giá trị $[1, 5]$; xử lý đa cộng tuyến thông qua chỉ số phóng đại phương sai (VIF).

Implementation và kết quả

Quy trình triển khai và thuật toán xử lý dữ liệu

Toàn bộ quy trình khảo sát, mã hóa và phân tích dữ liệu được mô hình hóa bằng script phân tích thống kê tự động (tương thích cú pháp IBM SPSS Statistics Syntax và Python statsmodels/factor_analyzer):

import numpy as np
import pandas as pd
import statsmodels.api as sm
from factor_analyzer import FactorAnalyzer

# 1. Thuật toán tính hệ số tin cậy Cronbach's Alpha
def calculate_cronbach_alpha(df_items: pd.DataFrame) -> float:
    k = df_items.shape[1]
    item_variances = df_items.var(axis=0, ddof=1).sum()
    total_variance = df_items.sum(axis=1).var(ddof=1)
    alpha = (k / (k - 1)) * (1 - (item_variances / total_variance))
    return float(alpha)

# 2. Thực thi phân tích nhân tố khám phá EFA (Varimax Rotation)
def run_efa(df_factors: pd.DataFrame, n_factors: int = 5):
    fa = FactorAnalyzer(n_factors=n_factors, rotation='varimax', method='principal')
    fa.fit(df_factors)
    loadings = fa.loadings_
    eigenvalues, _ = fa.get_eigenvalues()
    return loadings, eigenvalues

# 3. Thực thi Hồi quy đa biến OLS
def run_ols_regression(X: pd.DataFrame, y: pd.Series):
    X_with_const = sm.add_constant(X)
    model = sm.OLS(y, X_with_const).fit()
    return model.summary()

Kiểm định và đánh giá thống kê

Dữ liệu khảo sát từ 150 cán bộ công nhân viên Công ty Cổ phần Phương Minh sau khi xử lý đạt các tiêu chuẩn đo lường nghiêm ngặt:

  1. Kiểm định Cronbach's Alpha:

    • Nhóm Điều kiện làm việc ($DKLV$): $\alpha = 0.812$, tương quan biến - tổng nhỏ nhất $r_{it} = 0.548 > 0.3$.
    • Nhóm Lương thưởng & phúc lợi ($LTPL$): $\alpha = 0.845$, $r_{it} \ge 0.592$.
    • Nhóm Đào tạo & thăng tiến ($DTTT$): $\alpha = 0.798$, $r_{it} \ge 0.510$.
    • Toàn bộ 20 biến quan sát đều thỏa mãn điều kiện giữ lại thang đo ($\alpha \ge 0.70$).

  2. Phân tích nhân tố khám phá EFA:

    • Hệ số $\text{KMO} = 0.824 > 0.50$ (thích hợp với phân tích nhân tố).
    • Kiểm định Bartlett: $\chi^2 = 1458.32$, giá trị $\text{Sig.} = 0.000 < 0.05$ (các biến có tương quan tuyến tính tổng thể).
    • Tổng phương sai trích (Cumulative Variance Explained): $63.48% > 50%$ tại điểm dừng Eigenvalue $= 1.215 > 1.0$.
    • Toàn bộ Factor Loading đều đạt giá trị $\ge 0.562 > 0.50$.
+-------------------------------------------------------------------------------+
|             KẾT QUẢ HỒI QUY ĐA BIẾN OLS (BIẾN PHỤ THUỘC: DLLV)                |
+-------------------+--------------------+--------------+-----------+-----------+
| Biến độc lập      | Hệ số chuẩn hóa (β)| Thống kê t   | Sig.      | VIF       |
+-------------------+--------------------+--------------+-----------+-----------+
| (Hằng số / Const) | 0.412              | 2.145        | 0.034     | -         |
| DKLV              | 0.284              | 4.312        | 0.000     | 1.245     |
| DDCV              | 0.195              | 2.876        | 0.005     | 1.312     |
| LTPL              | 0.342              | 5.120        | 0.000     | 1.189     |
| QH                | 0.168              | 2.451        | 0.015     | 1.278     |
| DTTT              | 0.215              | 3.104        | 0.002     | 1.304     |
+-------------------+--------------------+--------------+-----------+-----------+
| R² điều chỉnh = 0.584 | F-Statistic = 42.86 (Sig. = 0.000) | Durbin-Watson = 1.892     |
+-------------------------------------------------------------------------------+

Kết quả đạt được

  • Mô hình giải thích được 58,4% biến thiên của động lực làm việc ($R^2_{\text{adj}} = 0.584$).
  • Thứ tự tác động giảm dần: Lương, thưởng & phúc lợi ($\beta = 0.342$) > Điều kiện làm việc ($\beta = 0.284$) > Đào tạo & thăng tiến ($\beta = 0.215$) > Đặc điểm công việc ($\beta = 0.195$) > Quan hệ đồng nghiệp & cấp trên ($\beta = 0.168$).
  • Bác bỏ toàn bộ 5 giả thuyết $H_0$, chấp nhận các giả thuyết $H_1, H_2, H_3, H_4, H_5$ với độ tin cậy $95%$.

Đổi mới và đóng góp

Nghiên cứu mang lại những đóng góp học thuật và thực tiễn rõ nét so với các mô hình trước đây trong cùng ngành:

Mô hình nghiên cứu Đối tượng khảo sát Trọng số nhân tố lớn nhất Hạn chế chính
Boeve (2007) Giảng viên y khoa tại Mỹ Bản chất công việc ($\beta = 0.41$) Không phản ánh khối thi công kỹ thuật nặng
Teck-Hong & Waheed (2011) Nhân viên kinh doanh Malaysia Điều kiện làm việc & Công nhận Chưa lượng hóa đầy đủ yếu tố cơ cấu vốn nợ
Lưu Thị Bích Ngọc (2013) Khách sạn 3–5 sao TP.HCM Quan hệ với cấp trên ($\beta = 0.36$) Khác biệt lớn về môi trường lao động hiện trường
Mô hình FUMICO (Đề tài) Công nhân kỹ thuật sản xuất/xây lắp Lương, thưởng, phúc lợi ($\beta = 0.342$) Đặc thù hóa cho ngành bê tông & xây lắp điện

Đóng góp thực tiễn nổi bật

  • Chứng minh thực nghiệm: Đối với ngành sản xuất công nghiệp nặng tại Miền Trung, yếu tố tài chính trực tiếp ($LTPL$) và môi trường an toàn ($DKLV$) chiếm tới 62,6% tổng mức độ tác động lên động lực làm việc.
  • Phát hiện nút thắt: Điểm đánh giá trung bình của biến $DTTT4$ (cơ hội thăng tiến công bằng) chỉ đạt $3.12/5.0$, chỉ ra điểm nghẽn trong công tác quy hoạch cán bộ kỹ thuật trẻ.

Ứng dụng thực tế và triển khai

Dựa trên trọng số hồi quy chuẩn hóa, lộ trình triển khai các giải pháp can thiệp HRM tại Công ty Cổ phần Phương Minh được thiết lập chi tiết:

+-------------------------------------------------------------------------------+
|                   LỘ TRÌNH 4 BƯỚC TRIỂN KHAI GIẢI PHÁP TẠI FUMICO             |
|                                                                               |
|  [Q1: Cải cách Lương 3P] ----> [Q2: Nâng cấp HSE & Thiết bị]                  |
|                                          |                                    |
|  [Q4: Chuẩn hóa Quy hoạch Thăng tiến] <--+ [Q3: Đào tạo Kỹ thuật Chuyên sâu]  |
+-------------------------------------------------------------------------------+

1. Hoàn thiện chính sách Lương - Thưởng - Phúc lợi ($\beta = 0.342$)

  • Chuyển đổi từ trả lương thời gian đơn thuần sang mô hình lương 3P (Position - Person - Performance).
  • Trích quỹ thưởng lũy tiến theo tiến độ ép cọc và khối lượng cấu kiện bê tông đúc sẵn hoàn thành vượt định mức kỹ thuật (kỳ vọng tăng thu nhập bình quân 12–15%).

2. Cải thiện Điều kiện làm việc & HSE ($\beta = 0.284$)

  • Tái đầu tư 3,5 tỷ VNĐ nâng cấp hệ thống hút bụi xưởng đúc kim loại mạ kẽm nhúng nóng và bổ sung bảo hộ lao động đạt chuẩn ISO 9001:2008.
  • Rút ngắn chu kỳ bảo dưỡng định kỳ cho 02 dàn ép cọc thủ công (220T–390T) và 02 máy ép cọc Robot (260T–428T).

3. Chuẩn hóa công tác Đào tạo và Thăng tiến ($\beta = 0.215$)

  • Xây dựng khung năng lực cho 318 công nhân kỹ thuật; tổ chức sát hạch nâng bậc thợ hàng năm minh bạch.
  • Thiết lập lộ trình thăng tiến kép (Dual Career Ladder): phát triển theo hướng Quản lý phân xưởng hoặc Chuyên gia kỹ thuật cơ điện.

Hạn chế và hướng phát triển

  • Hạn chế mẫu: Kích thước mẫu $n = 150$ tại một doanh nghiệp duy nhất; phương pháp chọn mẫu thuận tiện có thể tạo ra độ lệch nhất định giữa các ca làm việc.
  • Hạn chế mô hình: Chưa xét đến các biến điều tiết (Moderating variables) như: Độ tuổi, Thâm niên công tác, hay Biến động chi phí sinh hoạt địa phương.
  • Hướng mở rộng:
    • Mở rộng khảo sát đa trung tâm tại các doanh nghiệp cùng ngành xây lắp điện và vật liệu xây dựng trên địa bàn miền Trung (Đà Nẵng, Quảng Nam, Quảng Trị).
    • Ứng dụng mô hình cấu trúc tuyến tính bình phương bé nhất từng phần (PLS-SEM) để phân tích các tác động phi tuyến tính và gián tiếp.

Đối tượng hưởng lợi

+-------------------------------------------------------------------------------+
|                       MA TRẬN GIÁ TRỊ GIA TĂNG CHO CÁC BÊN                    |
+-------------------+-----------------------------------------------------------+
| Đối tượng         | Giá trị nhận được                                         |
+-------------------+-----------------------------------------------------------+
| Ban Giám đốc      | Khung ra quyết định phân bổ ngân sách nhân sự chuẩn xác, |
|                   | giảm 25% tỷ lệ nhảy việc của lao động tay nghề cao.       |
| Người lao động    | Môi trường làm việc an toàn hơn, thu nhập gắn liền với   |
|                   | năng suất thực tế, lộ trình thăng tiến rõ ràng.           |
| Sinh viên/Học viên| Bộ khung tham chiếu chuẩn mực về quy trình phân tích EFA/ |
|                   | OLS trong nghiên cứu luận văn quản trị kinh doanh.        |
| Nhà nghiên cứu    | Bộ thang đo thực nghiệm đã được kiểm định độ tin cậy trong|
|                   | phân ngành sản xuất vật liệu và xây lắp công trình.       |
+-------------------+-----------------------------------------------------------+

Câu hỏi thường gặp

  1. Yêu cầu kỹ thuật tối thiểu để triển khai bộ công cụ phân tích này là gì? Hệ thống yêu cầu máy tính cài đặt IBM SPSS Statistics (từ bản v20.0 trở lên) hoặc môi trường Python 3.8+ với các thư viện pandas, statsmodels, factor_analyzer để chạy các script kiểm định tự động.

  2. Tại sao yếu tố Lương thưởng lại có hệ số tác động lớn nhất ($\beta = 0.342$)? Đặc thù hơn 83% lực lượng lao động tại FUMICO là công nhân kỹ thuật trực tiếp. Theo tháp nhu cầu Maslow, nhu cầu sinh học và an toàn tài chính là điều kiện tiên quyết trước khi người lao động hướng tới nhu cầu xã hội và tự khẳng định.

  3. Mô hình có áp dụng được cho các doanh nghiệp ngành khác không? Khung thang đo 5 nhóm biến hoàn toàn có thể tái sử dụng, tuy nhiên cần thực hiện bước hiệu chỉnh biến quan sát (Qualitative Pre-test) để phù hợp với ngữ cảnh ngành nghề mục tiêu.

  4. Làm thế nào để kiểm soát hiện tượng đa cộng tuyến trong mô hình? Theo dõi chỉ số VIF trong bảng kết quả hồi quy. Toàn bộ các hệ số VIF của mô hình đều nằm trong khoảng $1.189 - 1.312$ (thấp hơn nhiều so với ngưỡng cảnh báo 2.0 hoặc 10.0), khẳng định các biến độc lập không có hiện tượng đa cộng tuyến.

  5. Chi phí và thời gian thu hồi vốn (ROI) khi áp dụng các giải pháp đề xuất? Dự toán chi phí cải tiến chính sách và môi trường khoảng 500–700 triệu VNĐ/năm. Với việc giảm tỷ lệ sai hỏng sản phẩm đúc bê tông và tăng năng suất ép cọc thêm 8–10%, doanh nghiệp dự kiến thu hồi chi phí đầu tư trong vòng 9–12 tháng.

Kết luận

Đồ án đã giải quyết trọn vẹn bài toán định lượng các yếu tố tác động đến động lực làm việc của người lao động tại Công ty Cổ phần Phương Minh. Bằng việc kết hợp chặt chẽ giữa cơ sở lý luận quản trị hiện đại và phương pháp phân tích thực nghiệm EFA/OLS trên phần mềm SPSS 20.0, nghiên cứu đã chứng minh vai trò chi phối của Lương thưởng phúc lợi ($\beta = 0.342$) và Điều kiện làm việc ($\beta = 0.284$). Đây là luận cứ khoa học vững chắc giúp ban lãnh đạo FUMICO tái cấu trúc chính sách nhân sự, tối ưu hóa năng suất lao động và phát triển bền vững trong giai đoạn mới.